我正在尝试从ESRI地理数据库(.gdb)中读取一个没有使用R的几何体的表.readOGR抛出错误,因为没有定义几何体,哪种有意义:
library(rgdal)
readOGR("gSSURGO_CO.gdb","mutext")
#Error in readOGR("gSSURGO_CO.gdb", "mutext") : no features found
#In addition: Warning message:
#In ogrInfo(dsn = dsn, layer = layer, encoding = encoding, use_iconv = use_iconv, :
#ogrInfo: all features NULL
我目前的解决方案是直接访问GDAL并将我读回的表的cvs转储回R:
# Here example with the same layer
system("ogr2ogr -f CSV mutext.csv gSSURGO_CO.gdb mutext")
mutext <- read.csv("mutext.csv")
它有效,但还有更好的方法吗?
谢谢!