Python OSGEO/GDAL ST_Intersection error -

i trying implement answer question: https://gis.stackexchange.com/questions/119374/intersect-shapefiles-using-shapely

but error:

runtimeerror: undefined function 'st_intersection' used.

here code. crashes on layer = ogr_ds.executesql(sql, dialect='sqlite') line. tried updating gdal library running pip vcvarsall.bat issue (also frustrating). appreciated!

from osgeo import ogr ogr.useexceptions() ogr_ds = ogr.open(r'd:\inputs', true)  sql = """\ select st_intersection(a.geometry, b.geometry) geometry, a.*, b.* stops a, zones b st_intersects(a.geometry, b.geometry); """ layer = ogr_ds.executesql(sql, dialect='sqlite') # copy result datasource new shapefile layer2 = ogr_ds.copylayer(layer, 'result') # save, close layer = layer2 = ogr_ds = none 

i think error arises because need sqlite driver compiled sqlite sql dialect. see: http://www.gdal.org/ogr_sql_sqlite.html

as updating , compiling gdal, can annoying. operating system running?