SQL queries in shapefiles using python/gdal return none result

Question

I want try use Python and GDAL to create some SQL queries in my shapefiles. I try something and it looks good, but doesn't work. I get a none result.

Here is the code :

from osgeo import ogr
ogr_ds = ogr.Open('line.shp')
TEST=3
sql = "SELECT id FROM {} where id='{}'".format(ogr_ds,TEST)
layer = ogr_ds.ExecuteSQL(sql)
print layer

• user30184: Does your SQL find something when you run it with ogrinfo?
• jessie jes: no I need like this in editor with python but I sure to have value
• user30184: I meant that do you know that the SQL part is OK? Test with ogrinfo -sql "SELECT id FROM line where id=3" line.shp. If you get "3" it is OK and the issue can only be in your Python code.
• jessie jes: yes I get it I take correct result now how to do in python ?
• user30184: I guess that you have only opened the GDAL datasource but not the layer and you should get first a layer to query with something like layer = dataSource.GetLayer(). Taken from an example pcjericks.github.io/py-gdalogr-cookbook/….

Answer 1

You are probably using the input layer in a wrong way; furthermore, also the .format{} operation doesn't seem correct.

You may try the following code:

from osgeo import ogr
filepath = 'C:/Users/path_to_the_shapefile/line.shp' # set the filepath
layer_name = filepath[:-4].split('/')[-1] # get the layer name
driver = ogr.GetDriverByName ("ESRI Shapefile")
ogr_ds = driver.Open(filepath)
TEST=3
sql = "SELECT id FROM %s WHERE id=%s" %(layer_name,TEST)
layer = ogr_ds.ExecuteSQL(sql)
print layer

I tested it with another shapefile (and a similar query) and it returns what expected:

print layer
<osgeo.ogr.Layer; proxy of <Swig Object of type 'OGRLayerShadow *' at 0x18A4D8F0> >

EDIT

If you want to print the id value, you may add these lines:

feat = layer.GetNextFeature()
val = feat.GetField(0)
print val

but this will only returns one feature and not all the queried features (this kind of operation would be useful if you are only interested in knowing if a specific value for the id field is stored in the input layer). Furthermore, you will get an error if there isn't any id value equals to the TEST variable.

• jessie jes: show me driver not define
• mgri: @jessiejes Sorry, there was a missing line. I edited the answer.
• jessie jes: now I take this <osgeo.ogr.Layer; proxy of <Swig Object of type 'OGRLayerShadow *' at 0x03166D88> > but I need the id for example if select id is 6 then I need to take only this value 6 in the print layer