Jak znaleźć maksymalną wartość komórki ze stosu rastrowego. Rmax <- maxValue(RAD1998.all[[1]]) działa dobrze, ale Rmax <- maxValue(RAD1998.all[[2]]) daje NA. Oczywiście nie w ustawionym w stos rastrze. Oto mój kod: RAD1998 <- raster(paste(getwd(), "/1998bil/1998ASC5min_ppt_spas1214_0001_19980202_0810_UTC.asc.bil", sep = "")) list.ras <- mixedsort(list.files(paste(getwd(), "/1998bil/", sep = ""), full.names = T, pattern = ".asc.bil")) RAD1998.all …
Obliczyłem powierzchnię rozmieszczenia gatunków (łączenie wielokątów z plików kształtów), ale ponieważ obszar ten może składać się z dość odległych wielokątów, chciałbym obliczyć pewną miarę dyspersji. Do tej pory udało mi się pobrać centroidy każdego wielokąta, obliczyć odległość między nimi i użyć ich do obliczenia współczynnika zmienności, jak w poniższym przykładzie …
wszystkim przepraszam, że przeszkadzam, ale jestem całkiem nowy, ponieważ stanąłem przed kluczową trudnością: chcę stworzyć animowaną mapę Russina ze zmianami bezrobocia w różnych latach. Na zdjęciu widać dane z jednego roku require(sp) require(maptools) require(RColorBrewer) require(rgdal) rus<-url("http://www.filefactory.com/file/4h1hb5c1cw7r/n/RUS_adm1_RData") print(load(rus)) unempl <- read.delim2(file="C:\\unempl1.txt", header = TRUE, sep = ";",quote = "", dec=",", stringsAsFactors=F) …
Interesuje mnie stworzenie powierzchni rastrowej, która przedstawia średnią intensywność promieniowania słonecznego na jednostkę powierzchni dla małej działki ogrodowej. Wyobrażam sobie, że dane wejściowe do modelu byłyby następujące: Powierzchnia działki - dla uproszczenia załóżmy, że działka ogrodowa jest okrągła. Średnia wysokość drzewa otaczającego skraju lasu Intensywność promieniowania słonecznego pod różnymi kątami …
Załączony obraz przedstawia szczelinę leśną z czerwoną sosną reprezentowaną przez koła i białą sosną reprezentowaną przez krzyże. Interesuje mnie ustalenie, czy istnieje pozytywny lub negatywny związek między dwoma gatunkami drzew sosny (tj. Czy rosną one na tych samych obszarach). Znam Kcrossa i Kmulti w pakiecie spatstat R. Ponieważ jednak mam …
W dokumentacji spwidziałem, że coordsargument SpatialPolygonsfunkcji akceptuje tylko 2 kolumny. Czy nie można tworzyć trójwymiarowych wielokątów? Robię przestrzenno-czasowe kriging za pomocą pakietu gstat, dodając czas jako trzeci wymiar przestrzenny SpatialPointsDataFrame. Chciałbym to rozszerzyć, aby zablokować kriging z wielokątami jako blokami. Czy istnieje jakiś sposób, aby to zrobić?
Potrzebuję siatki przestrzennej jako siatki głównej dla różnych map tematycznych. Jak utworzyć siatkę przestrzenną z rastra odrzucającego wszystkie piksele NA?
W moim obszarze roboczym jest 100 SpatialPointsDataFrames. Załadowałem je w ten sposób: filenames <- list.files(path="", pattern="XYhectareTravelTimes_ez+.*shp") for(i in filenames){ filepath <- file.path("/",i) assign(i, readShapePoints(filepath)) Nazywa się je tak: XYhectareTravelTimes_ez10.*shp XYhectareTravelTimes_ez11.*shp XYhectareTravelTimes_ez12.*shp itp. Jak przekonwertować je na rastry zapętlające się w przestrzeni roboczej? Jestem bardzo nowym użytkownikiem R i mam nadzieję …
Nawiązując do mojego poprzedniego pytania , w jaki sposób „geokodować” / przekonwertować - używając R - wygenerowanego obrazu rastrowego na format, który można wyświetlić jako warstwę na mapie? Krótko mówiąc, zamierzam odtworzyć wyniki komercyjnego oprogramowania SpatialKey ( zrzut ekranu ). Korzystam z instancji GeoServer do hostowania danych map, ale jak …
Prawdopodobnie bardzo proste pytanie. Mam listę około tysiąca potencjalnych lokalizacji geograficznych (lat lon) i spośród nich muszę wybrać 200 lokalizacji, które są „najbardziej rozproszone”. Myślę, że byłoby to 200 punktów z najwyższym całkowitym średnim dystansem. Pomyśl o sklepach w mieście. Czy istnieje określona metoda, aby to zrobić? Może w pakiecie …
To jest crossspost z stackoverflow . Bardzo mało wiem o oprogramowaniu GIS i robię wszystkie moje mapowania z góry R. Przeprosiny, jeśli jest to zbyt podstawowe pytanie. Załóżmy, że mam dwa pliki kształtów z różnych źródeł, ale o różnych atrybutach. Powiedzmy, że jedna dotyczy granic administracyjnych w Teksasie ( boundaries.shp), …
Mam zestaw danych w formie csv, na przykład: latitude, longitude, value -45, 45, 10 .... Na podstawie tych danych chciałbym obliczyć średnie na poziomie hrabstwa, korzystając z ArcGIS lub R + grass. Mam warstwę ArcGIS USA Counties, ale nie udało mi się ukończyć połączenia przestrzennego. W jaki sposób mogę to …
Używamy plików cookie i innych technologii śledzenia w celu poprawy komfortu przeglądania naszej witryny, aby wyświetlać spersonalizowane treści i ukierunkowane reklamy, analizować ruch w naszej witrynie, i zrozumieć, skąd pochodzą nasi goście.
Kontynuując, wyrażasz zgodę na korzystanie z plików cookie i innych technologii śledzenia oraz potwierdzasz, że masz co najmniej 16 lat lub zgodę rodzica lub opiekuna.