Użyłem integracji R z QGIS i uznałem, że jest to bardzo przydatne narzędzie. Po zapoznaniu się z rozszerzeniem R zestawu narzędzi dla 9 i rozszerzeniem dla 10 wydaje się, że oba oferują tylko klastrowanie punktowe, aby zademonstrować dowód słuszności koncepcji. Czy inne metody analizy dostępne w R można zaimplementować w …
W QGIS utworzyłem 2 wielokąty. Używając ich w R, wielokąty stają się automatycznie SpatialPolygonsDataFrame (SPDF). Chciałbym połączyć je w pojedynczy SPDF (co jest bardzo łatwe w ArcGis przy użyciu narzędzia Scalanie narzędzi ). Jestem pewien, że powinien istnieć prosty sposób na ukończenie tego w R, ale nie mogę tego znaleźć. …
Próbuję wykonać połączenie przestrzenne między danymi punktów i danymi wielokąta. Mam dane wskazujące współrzędne przestrzenne zdarzenia w moim pliku csv A i mam inny plik, plik kształtu B, który zawiera granice obszaru jako wielokąty. head(A) month longitude latitude lsoa_code crime_type 1 2014-09 -1.550626 53.59740 E01007359 Anti-social behaviour 2 2014-09 -1.550626 …
Próbuję odczytać tabelę bezpośrednio z geobazy bazy danych pliku ESRI do R. Przykładowy plik danych można pobrać tutaj . Baza danych zawiera punktową klasę obiektów (Zone9_2014_01_Broadcast) i dwie połączone tabele (Zone9_2014_01_Vessel i Zone9_2014_01_Voyage). Możesz odczytać plik kształtu w R, używając readOGRz rgeospakietu: library(rgeos) library(downloader) download("https://coast.noaa.gov/htdata/CMSP/AISDataHandler/2014/01/Zone9_2014_01.zip", dest="Zone9_2014_01.zip", mode="wb") unzip("Zone9_2014_01.zip", exdir = …
Czy istnieje funkcja przycinania obiektu mapy sf, podobna do tej maptools::pruneMap(lines, xlim= c(4, 10), ylim= c(10, 15))używanej w przypadku SpatialPolygon lub SpatialLine? Rozważam, st_intersection()ale może być właściwy sposób.
Biorąc pod uwagę dużą (~ 1 milion) próbkę nierównomiernie rozmieszczonych punktów - czy można wygenerować nieregularną siatkę (pod względem wielkości, ale może również mieć nieregularny kształt, jeśli to możliwe?), Która będzie zawierać określoną minimalną ilość n punktów? Dla mnie nie ma większego znaczenia, jeśli generetowane „komórki” takiej siatki zawierają dokładnie …
Jak przekonwertować punkty lon-lat na proste funkcje (sfg), a następnie umieścić je w prostej kolekcji funkcji (sfc)? Oto MWE, które nie działa, ale jest najbliższe, na jakie się natknąłem. library(data.table) library(sf) # The DT data.table is the data I have (but 10,000s of rows, each row is a point) DT …
Planuję niektóre geolokalizacje (lon, lat)tweetów zebranych za pomocą Stream API z 7mi. promień. Poniższy obrazek pokazuje dwie mapy tych samych danych, ale przy użyciu różnych wartości alfa / przezroczystości dla każdej obserwacji. Patrząc na mapę po lewej stronie zauważyłem, że geolokalizacje podążają głównie za wyobrażoną „siatką” z określonymi węzłami. Czy …
coast<-readShapeSpatial("coastline.shp") landc<-readShapeSpatial("landcover.shp") ro<-readShapeSpatial("roads.shp") bc<-gBuffer(ro,width=100) landc$ratings=1 landc$ratings[landc$LANDUSE_ID==4]=0 Powyżej biorę dowolną kategorię, która ma 4, i w nowej kolumnie wpisuję ją jako 0. W tym momencie chcę również nazwę kolumny ratingsdla bc, gdzie zajmie 0, jeśli jest w buforze i 1, jeśli jest na zewnątrz. Problem polega na tym, że bcjest SpatialPolygonsi …
Mam zestaw obrazów rastrowych, które reprezentują konkretny miesiąc na przestrzeni lat i chcę utworzyć oś czasu średnich wartości obszaru za pomocą pliku kształtu. Jak wyodrębnić wartości z rastrów i zaimportować je w R w sposób, w jaki mogę je wykorzystać?
I Morana , miara autokorelacji przestrzennej, nie jest szczególnie solidną statystyką (może być wrażliwa na wypaczone rozkłady atrybutów danych przestrzennych). Jakie są bardziej niezawodne techniki pomiaru autokorelacji przestrzennej? Szczególnie interesują mnie rozwiązania, które są łatwo dostępne / możliwe do wdrożenia w języku skryptowym, takim jak R. Jeśli rozwiązania dotyczą wyjątkowych …
Zillow ma zestaw plików kształtowych dla różnych dzielnic głównych amerykańskich miast. Chciałem sprawdzić, czy niektóre budynki były obecne w niektórych dzielnicach za pomocą R: library(rgeos) library(sp) library(rgdal) df <- data.frame(Latitude =c(47.591351, 47.62212,47.595152), Longitude = c(-122.332271,-122.353985,-122.331639), names = c("Safeco Field", "Key Arena", "Century Link")) coordinates(df) <- ~ Latitude + Longitude wa.map …
Mam dwa wielokąty. Jedno zawiera pola (X, Y, Z), a drugie zawiera typy gleby (A, B, C, D). Chcę wiedzieć, który obszar na każdym polu zawiera rodzaj gleby. Próbowałem następujące: library(rgdal) library(rgeos) Field<-readOGR("./","Field") Soil<-readOGR("./","Soil") Results<-gIntersects(Soil,Field,byid=TRUE) rownames(Results)<-Field@data$FieldName colnames(Results)<-Soil@data$SoilType > Results A B C D Z TRUE FALSE FALSE FALSE Y FALSE …
Próbuję wykonać połączenie na wspólnym polu po scaleniu dwóch sąsiednich plików kształtów. Pliki kształtów kończą się co najmniej jednym cienkim kawałkiem przestrzeni między nimi. Podczas próby połączenia pojawia się następujący błąd osieroconej dziury: Błąd w createPolygonsComment (p): rgeos_PolyCreateComment: osierocona dziura, nie można znaleźć zawierającego wielokąt dla dziury o indeksie 17 …
Moim celem jest zmodyfikowanie istniejącego pliku kształtu przez scalenie niektórych wielokątów. Po zaimportowaniu pliku kształtu i użyciu polecenia UnionSpatialPolygons otrzymuję kontur wielokąta, który chcę. Jednak jest to teraz obiekt SpatialPolygons, a nie SpatialPolygonsDataFrame, więc nie jestem w stanie wyeksportować go do pliku kształtu za pomocą writeOGR. Jak mogę obejść ten …
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.