Take the 2-minute tour ×
Geographic Information Systems Stack Exchange is a question and answer site for cartographers, geographers and GIS professionals. It's 100% free, no registration required.

I am trying to access a vector dataset in a ESRI File Geodatabase using Python + GDAL. I have successfully compiled GDAL with the file geodatabase API. The FileGDB driver is working correctly since entering

ogrinfo --formats

shows the FileGDB driver and entering

ogrinfo myfilegdb.gdb 

gives me the correct information on the content of the database.

However, I can't find out how to access the content itself in Python. For accessing a shapefile, I would write:

driver = ogr.GetDriverByName('ESRI Shapefile')
ds = driver.Open('shapefile.shp', 0)

When accessing a FileGDB feature class I would assume using the commands:

driver = ogr.GetDriverByName('FileGDB')
ds = driver.Open('myfilegdb.gdb/feature_class', 0)

but this doesn't seem to work since it can not identify/locate the data set. Does anyone know how to call individual feature classes from a ESRI FileGDB.

I'm using Python 2.7, GDAL 1.9.1, filegdb api 1.2 on Ubuntu 12.04 x64. Thanks for any suggestions!

share|improve this question
    
Can you give me some idea to install the OGR for FileGDB driver? –  giser Dec 11 '13 at 4:06

1 Answer 1

up vote 10 down vote accepted

You're almost there. This is on Windows 7, Python 2.6.5 32bit, and GDAL 1.9.0:

>>> from osgeo import ogr
>>> driver = ogr.GetDriverByName("FileGDB")
>>> ds = driver.Open(r"C:\temp\buildings.gdb", 0)
>>> ds
<osgeo.ogr.DataSource; proxy of <Swig Object of type 'OGRDataSourceShadow *' at 0x02BB7038> >
>>> ds.GetLayer("buildings")
<osgeo.ogr.Layer; proxy of <Swig Object of type 'OGRLayerShadow *' at 0x02BB7050> >
>>> b = ds.GetLayer("buildings")
>>> sr = b.GetSpatialRef()
>>> sr
<osgeo.osr.SpatialReference; proxy of <Swig Object of type 'OSRSpatialReferenceShadow *' at 0x02BB7080> >
>>> sr.ExportToProj4()
'+proj=utm +zone=15 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs '
>>>

Once you open the FGDB, then use GetLayer to get at your featureclass.

share|improve this answer
    
Once you know it it seems so logical :-) thanks a lot, your solution does the trick. –  Niels Sep 4 '12 at 13:46

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.