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?
Yanıtlar:
Esri dosya coğrafi veritabanlarındaki özellik sınıflarına erişmek için rgdal kullanabilirsiniz .
require(rgdal)
# The input file geodatabase
fgdb <- "C:/path/to/your/filegeodatabase.gdb"
# List all feature classes in a file geodatabase
subset(ogrDrivers(), grepl("GDB", name))
fc_list <- ogrListLayers(fgdb)
print(fc_list)
# Read the feature class
fc <- readOGR(dsn=fgdb,layer="some_featureclass")
# Determine the FC extent, projection, and attribute information
summary(fc)
# View the feature class
plot(fc)layerGDB'de sadece bir özellik sınıfı olduğunda , argümanın atlanabileceğini bilmek yardımcı olur .
gdbDosyalar hakkında hiçbir şey bilmeyenler için , fgdbbu cevapta bir dizin var ve ogrListLayers()bu dizinde çalışıyor ...
Bu cevapta zaten yayınlandığı gibi , bu da şu şekilde çok iyi çalışıyor sf:
require(sf)
fc <- sf::st_read("C:/path/to/your/filegeodatabase.gdb", layer = "some_featureclass")Ancak bir fgdb'ye yazılmadıysa ( henüz? ), Bir ArcGIS / ArcMap Lisansının yanı sıra R kütüphanesine sahip olmanız gerekir arcgisbinding(bkz. Https://github.com/R-ArcGIS/r-bridge ).
st_drivers()$write[st_drivers()$long_name == "ESRI FileGDB"]
#> [1] FALSE