«rgdal» etiketlenmiş sorular

Dışa aktarılan bir şekil dosyası olarak 2 GB'den büyük bir coğrafi veritabanında bulunan bir özelliğe sahibim. Çokgenleri bir raster dosyasındaki verilerle ilişkilendirmek için R'deki bir çıkarma işlevini çalıştırmam gerekiyor. Özelliğin bir tablo olarak dışa aktarılması bir çözüm değildir. Bir Esri dosyası coğrafi veritabanında bulunan özellik sınıflarını nasıl okuyabilirim?

27 r  features  file-geodatabase  feature-class  rgdal 

Ben indirilebilir R. Bir örnek veri dosyası içine bir ESRI dosya coğrafi veritabanından doğrudan bir tablo okumaya çalışıyorum burada . Veritabanı bir nokta özellik sınıfı (Zone9_2014_01_Broadcast) ve iki bağlantılı tablo (Zone9_2014_01_Vessel ve Zone9_2014_01_Voyage) içerir. R'deki şekil readOGRdosyasını rgeospaketten kullanarak okuyabilirsiniz : 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 = ".") …

21 r  file-geodatabase  table  rgdal 

Zillow, büyük ABD şehirlerinin farklı mahalleleri için bir dizi şekil dosyasına sahiptir . R kullanarak belirli mahallelerde belirli binaların olup olmadığını kontrol etmek istedim: 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 

İki çokgenim var. Biri tarlaları (X, Y, Z), diğeri toprak türlerini (A, B, C, D) içerir. Her alanın hangi alanında hangi toprak türünü içerdiğini bilmek istiyorum. Aşağıdakileri denedim: 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 TRUE TRUE …

19 shapefile  polygon  r  intersection  rgdal 

İki çokgen oluşturmak istiyorum. Bunlardan biri rasterimin dikdörtgen uzantılarından. Kullanarak raster boyutunu elde edebilirsiniz biliyorum r <- raster("band5.tif") e <- extent(r) plot(e)bana bunu verir: Aşağıda gösterildiği gibi raster sınırının BİR poligonunu nasıl oluşturabilirim? Rasterim şöyle görünüyor:

13 raster  r  rgdal 

gSimplifyBir şekil dosyasının geometrilerini basitleştirmek için (rgeos paketi) kullanıyorum . İşlev iyi çalışıyor, ancak şimdi çıktıyı yeni bir şekil dosyasına yazamıyorum. Bazı yollar denedim: writeOGR(simplyshape, file, driver="ESRI Shapefile", layer='test') alırım obj bir SpatialPointsDataFrame, SpatialLinesDataFrame veya SpatialPolygonsDataFrame olmalıdır Ve birlikte: writePolyShape(simplyshape, file) Alırım: Hata: (x, "SpatialPolygonsDataFrame") DOĞRU değil

12 shapefile  r  simplify  rgdal 

Aynı GeoPackage R birden çok katman yazmaya çalışıyorum, ama bir hata alıyorum Creation of output file failed. RGDAL ile .gpkg dosyalarına okuma ve yazma belgelerini aramaya çalıştım, örneğin, writeOGR()aslında çok sayıda katmanı destekleyip desteklemediğini anlamak için . Bu mümkün mü, eğer öyleyse, nasıl yapılır? Minimum çalışma örneği: library(sp) library(maptools) library(rgdal) …

12 gdal  r  rgdal  geopackage  sf 

Bende NE lat lngve var SW lat lng. Amacım basit, RStudio kullanıyorum ve yukarıdaki iki ağıt ve uzuntan bir sınırlayıcı kutu oluşturmak ve daha sonra Çokgen oluşturmak ve kare çokgene rastgele noktalar eklemek istiyorum. Bir data.frame içine aktarabileceğiniz bir csv dosyasında NELat / Lng ve SW-Lat / Lng var. Başlamak …

11 r  rgdal 

Windows 7'de R'ye okumak istediğim çok büyük bir ESRI şekil dosyası var (> 8000 çokgen, .shp dosyası> 32MB). Rgdal ve readOGR () yöntemine aşinayım. Bununla birlikte, bu şekil dosyası boyutuyla, tüm şekil dosyasını belleğe okumak çok yavaş bir işlemdir. Ayrıca, şekil dosyasında bulunan 8.000'in yalnızca birkaç kümelenmiş çokgeniyle (<100) ilgileniyorum. …

10 shapefile  r  rgdal 

R'de raster vektörü dönüştürüyorum. Ancak süreç çok uzundu. Daha hızlı yapmak için betiği çok iş parçacığına veya GPU işlemeye koyma olasılığı var mı? Rasterleştirilmiş vektöre benim senaryom. r.raster = raster() extent(r.raster) = extent(setor) #definindo o extent do raster res(r.raster) = 10 #definindo o tamanho do pixel setor.r = rasterize(setor, r.raster, …

9 raster  r  rasterization  rgdal  parallel