Pytania otagowane jako rgdal

Wiązania R dla biblioteki abstrakcji danych geoprzestrzennych

5
Odczytaj tabelę z geobazy danych pliku ESRI (.gdb) za pomocą R.
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 = …
3
Sprawdzanie, czy punkty mieszczą się w wielokącie Shapefile
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 …
19 r  rgdal  zillow 
2
Wyodrębnianie obszarów przecięcia w R.
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 …
1
Utwórz wielokąty zasięgu i granicy danego rastra w R.
Chcę utworzyć dwa wielokąty. Jeden jest z prostokątnych rozmiarów mojego rastra. Wiem, że mogę uzyskać zasięg użycia rastra, r <- raster("band5.tif") e <- extent(r) plot(e)daje mi to: Jak mogę utworzyć JEDEN wielokąt granicy rastra, jak pokazano poniżej? Tak wygląda mój raster:
13 raster  r  rgdal 
2
writeOGR z wielokątem przestrzennym uproszczonym przez gSimplify
Używam gSimplify(pakiet rgeos), aby uprościć geometrię pliku kształtu. Funkcja działa dobrze, ale teraz nie mogę zapisać wyniku w nowym pliku kształtu. Próbowałem na kilka sposobów: writeOGR(simplyshape, file, driver="ESRI Shapefile", layer='test') dostaję obj musi być SpatialPointsDataFrame, SpatialLinesDataFrame lub SpatialPolygonsDataFrame i z: writePolyShape(simplyshape, file) Dostaję: Błąd: is (x, „SpatialPolygonsDataFrame”) nie jest PRAWDĄ
12 shapefile  r  simplify  rgdal 
2
R - Utwórz obwiednię, przekonwertuj na klasę Polygon i wydrukuj
Mam NE lat lngi SW lat lng. Mój cel jest prosty, używam RStudio i chcę utworzyć obwiednię z powyższych dwóch łat i długich, a następnie wygenerować wielokąt i chcę dodać losowe punkty do wielokąta kwadratowego. Mam NELat / Lng i SW-Lat / Lng w pliku csv, który mógłbym zaimportować do …
11 r  rgdal 
1
Zapisywanie wielu warstw w GeoPackage za pomocą writeOGR () w R?
Próbuję napisać wiele warstw do tej samej GeoPackage w badania, ale pojawia się błąd Creation of output file failed. Próbowałem wyszukać dokumentację odczytu i zapisu do plików .gpkg za pomocą RGDAL, np. Aby dowiedzieć się, czy writeOGR()rzeczywiście obsługuje wiele warstw, z niewielkim powodzeniem. Czy to w ogóle możliwe, jeśli tak, …
11 gdal  r  rgdal  geopackage  sf 
3
Czytaj tylko część pliku kształtu ESRI do R
Mam bardzo duży plik kształtu ESRI (> 8000 wielokątów, plik .shp to> 32 MB), który chcę wczytać do R na Windows 7. Znam rgdal i readOGR (). Jednak przy takim rozmiarze pliku shapefile wczytywanie całego pliku shapefile do pamięci jest procesem bardzo powolnym. Ponadto interesuje mnie tylko kilka zgrupowanych wielokątów …
10 shapefile  r  rgdal 
1
Szybsze przetwarzanie wektora na raster za pomocą R.
Przekształcam wektor na raster w R. Jednak proces ten był zbyt długi. Czy istnieje możliwość włączenia skryptu w proces wielowątkowy lub przetwarzanie GPU, aby zrobić to szybciej? Mój skrypt do zrasteryzowanego wektora. r.raster = raster() extent(r.raster) = extent(setor) #definindo o extent do raster res(r.raster) = 10 #definindo o tamanho do …
Korzystając z naszej strony potwierdzasz, że przeczytałeś(-aś) i rozumiesz nasze zasady używania plików cookie i zasady ochrony prywatności.
Licensed under cc by-sa 3.0 with attribution required.