Tôi có một tính năng chứa trong cơ sở dữ liệu địa lý lớn hơn 2GB dưới dạng shapefile xuất khẩu. Tôi cần chạy một hàm trích xuất trong R để gán các đa giác với dữ liệu từ tệp raster. Xuất khẩu tính năng dưới dạng bảng không phải là một giải pháp. Làm cách nào tôi có thể đọc các lớp đối tượng có trong cơ sở dữ liệu địa lý tệp Esri?
Câu trả lời:
Bạn có thể sử dụng rgdal để truy cập các lớp tính năng trong cơ sở dữ liệu địa lý tệp Esri.
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)layerđối số có thể được bỏ qua khi chỉ có một lớp tính năng trong GDB.
gdbcác tập tin, fgdbtrong câu trả lời này, đây là một thư mục và ogrListLayers()hoạt động trên thư mục này ...
Như đã được đăng trong câu trả lời này, giờ đây nó cũng hoạt động rất độc đáo với sf:
require(sf)
fc <- sf::st_read("C:/path/to/your/filegeodatabase.gdb", layer = "some_featureclass")Nhưng viết vào fgdb ist chưa được triển khai ( chưa? ), Bạn phải có Giấy phép ArcGIS / ArcMap cũng như thư viện R arcgisbinding(xem https://github.com/R-ArcGIS/r-bridge )
st_drivers()$write[st_drivers()$long_name == "ESRI FileGDB"]
#> [1] FALSE