I am trying to read shapefiles contained within a geodatabase file (.gdb) into R. The .gdb contains two Feature Datasets with multiple Feature Classes within each.
The problem is only one of the two Feature Datasets is being read. Specifically, I am using the NHD dataset for all states ftp://nhdftp.usgs.gov/DataSets/Staged/States/FileGDB/HighResolution/ In each state .gdb are two feature Datasets, 'WBD' and 'Hydrography'. Only 'WBD' is being read. Using ogrListLayers only returns layers from 'WBD'. Can I specify the 'Hydrography' dataset and access the Feature Classes within it?
When I run orgListLayers, it returns:
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" Yet others get:
ogrListLayers("NHDH_VI.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" Why the difference? Specifically I am looking for: Layer name: NHDWaterbody.
Any suggestions would be much appreciated. R version 3.2.0. OSX v.10.10.3
أكثر...
The problem is only one of the two Feature Datasets is being read. Specifically, I am using the NHD dataset for all states ftp://nhdftp.usgs.gov/DataSets/Staged/States/FileGDB/HighResolution/ In each state .gdb are two feature Datasets, 'WBD' and 'Hydrography'. Only 'WBD' is being read. Using ogrListLayers only returns layers from 'WBD'. Can I specify the 'Hydrography' dataset and access the Feature Classes within it?
When I run orgListLayers, it returns:
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" Yet others get:
ogrListLayers("NHDH_VI.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" Why the difference? Specifically I am looking for: Layer name: NHDWaterbody.
Any suggestions would be much appreciated. R version 3.2.0. OSX v.10.10.3
أكثر...