Szukam innego, bardziej eleganckiego rozwiązania problemu statystyki przestrzennej. Surowe dane składają się ze współrzędnej xy dla każdego pojedynczego drzewa (tj. Przekonwertowanego na punktowy plik .shp). Chociaż nie zostało użyte w tym przykładzie, każde drzewo ma również odpowiedni wielokąt (tj. Jako .shp), który reprezentuje średnicę korony. Dwa obrazy po lewej pokazują …
chcę podzielić plik kształtu z wieloma płytkami na wiele mniejszych płytek. Poszukuję inteligentnego rozwiązania, ponieważ ręczne edytowanie nie wchodzi w rachubę Czy ktoś może pomóc lub znać fajną funkcję?
Próbuję przekonwertować plik kształtu do rastra w R. Moje podejście polega na tym, aby czytać w rastrze w następujący sposób: library(rgdal) # Loads SP package by default demo <- readOGR('F:/data/', 'shapefile') # Creates a SpatialPolygonsDataFrame class (sp) To działa dobrze i mogę to wykreślić. Jest to jednak duży plik kształtów …
Podane punkty danych o długości i szerokości geograficznej oraz trzeciej wartości właściwości tego punktu. Jak mogę grupować punkty w grupy (podregiony geograficzne) na podstawie wartości nieruchomości? Przeszukałem google i odkryłem, że problem ten nazywa się „klastrowaniem przestrzennym” lub „regionalizacją”. Nie znam się jednak na obsłudze danych geograficznych i nie mam …
Jak uprościć sfwielokąt bez wprowadzania przerw i fragmentów? Na przykład z plikiem shapefile użyłbym rmapshaper::ms_simplify(): library("pryr") library("rgdal") library("rmapshaper") download.file("https://borders.ukdataservice.ac.uk/ukborders/easy_download/prebuilt/shape/England_gor_2011.zip", destfile = "regions.zip") unzip("regions.zip") regions <- readOGR(".", "england_gor_2011") object_size(regions) # ~13MB regions <- ms_simplify(regions) object_size(regions) # < 1MB Próbowałem, sf::st_cast()który na stronach podręcznika stwierdza: Rzuć geometrię na inny typ: uprość lub …
W pakiecie HistData dla R ( https://r-forge.r-project.org/R/?group_id=574 ) mam zestawy danych związane z mapą wybuchu cholery w Londynie w 1854 roku przez Johna Snowa. Wierzę, że są autorytatywne, ponieważ zostały starannie zdigitalizowane pod nadzorem Waltera Toblera. Niektóre szczegóły dotyczące tych zestawów danych zostały opisane przez Johna Mackenzie na stronie http://www1.udel.edu/johnmack/frec480/cholera/cholera2.html …
Jestem absolutnym początkującym użytkownikiem danych geograficznych, więc proszę wybacz mi, jeśli pytanie nie jest właściwe. Pobrałem dane z NCDC NARR i udało mi się załadować do R przy użyciu rasterpakietu. Chciałbym uzyskać listę z szerokością, długością i wartością. Rozumiem, że rasterToPoints()powinienem robić dokładnie to, co chcę, jednak moje wartości szerokości …
Odwrotny klip zapisuje tylko część obiektu przestrzennego, która znajduje się poza granicami innego obiektu, w przeciwieństwie do zwykłego klipu, który zapisuje części znajdujące się wewnątrz innego obiektu. Wykonujesz odwrotny klip w ArcMap? pokazuje, jak to zrobić w ArcMap. Jak to zrobić w R? Powtarzalny przykład (na komputerach z systemem Linux): …
Mam plik kształtu linii reprezentujący sieć dróg. Chcę zrasteryzować te dane, a wartości wynikowe w rastrze pokazują całkowitą długość linii, które mieszczą się w komórce rastra. Dane są w projekcji British National Grid, więc jednostkami będą metry. Idealnie chciałbym wykonać tę operację przy użyciu Ri domyślam się, że rasterizefunkcja z …
Uważam, że duża część dokumentacji QGIS jest ukierunkowana na użycie GUI. Chcę mieć możliwość uruchamiania skryptów R w QGIS, ale mam problem ze znalezieniem zasobu, który ugruntuje mnie w relacji między strukturami danych QGIS i R oraz w jaki sposób przekazywać struktury między nimi. Powiązane pytanie ( czy możliwe jest …
Mam plik kształtu wielokątów i inny plik CSV, który zawiera listę punktów w postaci par (Lat, Lng) .. Chcę sprawdzić dla każdej pary (łac., Lng) z pliku CSV, w którym wielokącie wpada. Plik shapefile jest wyświetlany, a plik proj wygląda następująco: PROJCS["Transverse_Mercator",GEOGCS["GCS_OSGB 1936", DATUM["D_OSGB_1936",SPHEROID["Airy_1830",6377563.396,299.3249646]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",49],PARAMETER["central_meridian",-2],PARAMETER["scale_factor",0.9996012717],PARAMETER["false_easting",400000],PARAMETER["false_northing",-100000],UNIT["Meter",1]] Mój plan jest następujący: Odczytaj plik …
Załączony zestaw danych pokazuje około 6000 sadzonek w około 50 lukach leśnych o zmiennej wielkości. Interesuje mnie nauka, w jaki sposób te sadzonki rosną w obrębie swoich odpowiednich luk (tj. Skupione, losowe, rozproszone). Jak wiecie, tradycyjnym podejściem byłoby uruchomienie Globalnego Morana I. Jednak agregacje drzew w obrębie agregatów luk wydają …
Jestem bardzo nowy w korzystaniu z danych GIS i mam jedynie niewielkie doświadczenie z R. Czytałem o tym, jak analizować dane przestrzenne za pomocą książki PDF spatial-analyst.net, więc nie jestem całkowicie zagubiony, ale pomyślałem, że mógłbym opisać mój problem i ludzie mogą sugerować pomysły. Mam zestaw danych z około 2000 …
Czytałem ostatnio o kartogramach w tle i natknąłem się na interesującą alternatywę - mapy wartości według alfa . Koncepcja jest opisana bardziej szczegółowo tutaj, a pełny papier jest dostępny [nie otwarcie] [tutaj] 2 . Istnieje kilka szczegółów dotyczących implementacji na blogu indiemaps : ... wszystkie powyższe grafiki zostały utworzone za …
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.