6

我正在尝试将地理数据库文件 (.gdb) 中包含的 shapefile 读取到 R 中。.gdb 包含两个要素数据集，每个要素类中包含多个要素类。

问题是只有两个特征数据集之一被读取。具体来说，我正在使用所有州的 NHD 数据集ftp://nhdftp.usgs.gov/DataSets/Staged/States/FileGDB/HighResolution/ 在每个州.gdb 中有两个特征数据集，“WBD”和“水文”。仅读取“WBD”。使用 ogrListLayers 仅从“WBD”返回层。如何指定“水文”数据集并访问其中的要素类？

任何建议将不胜感激。R 版本 3.2.0。OSX v.10.10.3

编辑 06/16/15：当我运行 orgListLayers 时，它返回：

ogrListLayers('NHDH_VI.gdb')
[1] "ExternalCrosswalk"       "NHDFCode"                "NHDFeatureToMetadata"   
[4] "NHDFlow"                 "NHDFlowlineVAA"          "NHDMetadata"            
[7] "NHDProcessingParameters" "NHDReachCodeMaintenance" "NHDReachCrossReference" 
[10] "NHDSourceCitation"       "NHDStatus"               "NHDVerticalRelationship"
[13] "WBDHU14"                 "WBDHU8"                  "WBDHU2"                 
[16] "WBDHU4"                  "WBDHU6"                  "WBDHU10"                
[19] "WBDHU12"                 "WBDHU16"                 "HYDRO_NET_Junctions"

这 21 层是不同的，然后由 Mike T 和 hrbrmstr 表示。具体来说，我正在寻找：图层名称：NHDWaterbody。

当我ogrinfo -ro NHDH_VI.gdb从终端运行时。

ERROR 1: Error: Failed to open Geodatabase (This release of the GeoDatabase is either invalid or out of date.)
FAILURE: Unable to open datasource `NHDH_VI.gdb' with the following drivers.
      -> FileGDB
      -> OpenFileGDB
      ...
4

2 回答 2

6

您可能正在使用ESRI 文件地理数据库 (OpenFileGDB)驱动程序读取文件。OGR 不保留或使用要素数据集，因此所有要素类都混合在同一个平面命名空间中。在 ArcCatalog 中查看（例如）NHDH_VI.gdb：

弧形目录

然后从命令行提示符读取相同的文件ogrinfo -ro NHDH_VI.gdb

INFO: Open of `NHDH_VI.gdb'
      using driver `OpenFileGDB' successful.
1: NHDPoint (Point)
2: NHDFlowline (Multi Line String)
3: NHDLine (Multi Line String)
4: NHDArea (Multi Polygon)
5: NHDWaterbody (Multi Polygon)
6: NHDAreaEventFC (Multi Polygon)
7: NHDLineEventFC (Multi Line String)
8: NHDPointEventFC (Point)
9: WBDLine (Multi Line String)
10: NonContributingDrainageArea (Multi Polygon)
11: NWISBoundary (Multi Line String)
12: NWISDrainageArea (Multi Polygon)
13: WBDHU14 (Multi Polygon)
14: WBDHU8 (Multi Polygon)
15: WBDHU2 (Multi Polygon)
16: WBDHU4 (Multi Polygon)
17: WBDHU6 (Multi Polygon)
18: WBDHU10 (Multi Polygon)
19: WBDHU12 (Multi Polygon)
20: WBDHU16 (Multi Polygon)
21: HYDRO_NET_Junctions (Point)

同样可以从 R 中获得：

> library(rgdal)
> ogrListLayers("NHDH_VI.gdb")
 [1] "NHDPoint"                    "NHDFlowline"                
 [3] "NHDLine"                     "NHDArea"                    
 [5] "NHDWaterbody"                "NHDAreaEventFC"             
 [7] "NHDLineEventFC"              "NHDPointEventFC"            
 [9] "WBDLine"                     "NonContributingDrainageArea"
[11] "NWISBoundary"                "NWISDrainageArea"           
[13] "WBDHU14"                     "WBDHU8"                     
[15] "WBDHU2"                      "WBDHU4"                     
[17] "WBDHU6"                      "WBDHU10"                    
[19] "WBDHU12"                     "WBDHU16"                    
[21] "HYDRO_NET_Junctions"        
attr(,"driver")
[1] "OpenFileGDB"
attr(,"nlayers")
[1] 21

因此，您需要手动将数据集从 ArcCatalog 层次结构过滤到您可以从 OGR 中找到的内容。并非所有类都可用于 OGR（例如非空间表、关系类）。

于 2015-06-16T03:57:25.440 回答
0

这不是答案，但需要更多空间来澄清，因为 OP 没有提供可重复的示例。一旦 OP q 更完整，我将删除它。

当我ogrinfo在NHDH_CT.gdb（较小的档案之一）上运行时，我得到：

Layer name: NHDPoint
Layer name: NHDFlowline
Layer name: NHDLine
Layer name: NHDArea
Layer name: NHDWaterbody
Layer name: NHDAreaEventFC
Layer name: NHDLineEventFC
Layer name: NHDPointEventFC
Layer name: WBDLine
Layer name: NonContributingDrainageArea
Layer name: NWISBoundary
Layer name: NWISDrainageArea
Layer name: WBDHU14
Layer name: WBDHU8
Layer name: WBDHU2
Layer name: WBDHU4
Layer name: WBDHU6
Layer name: WBDHU10
Layer name: WBDHU12
Layer name: WBDHU16
Layer name: HYDRO_NET_Junctions

那匹配：

> ogrListLayers("NHDH_CT.gdb")
 [1] "NHDPoint"                    "NHDFlowline"                 "NHDLine"                    
 [4] "NHDArea"                     "NHDWaterbody"                "NHDAreaEventFC"             
 [7] "NHDLineEventFC"              "NHDPointEventFC"             "WBDLine"                    
[10] "NonContributingDrainageArea" "NWISBoundary"                "NWISDrainageArea"           
[13] "WBDHU14"                     "WBDHU8"                      "WBDHU2"                     
[16] "WBDHU4"                      "WBDHU6"                      "WBDHU10"                    
[19] "WBDHU12"                     "WBDHU16"                     "HYDRO_NET_Junctions"

你期待看到什么？你的ogrinfo节目是什么？（注意：夏威夷给出了类似的输出）。

于 2015-06-15T19:04:02.150 回答