python - Python OSGEO/GDAL ST_Intersection エラー

Question

この質問への回答を実装しようとしています: https://gis.stackexchange.com/questions/119374/intersect-shapefiles-using-shapely

しかし、私はエラーが発生します:

RuntimeError: 未定義の関数 'ST_Intersection' が使用されました。

これがコードです。layer = ogr_ds.ExecuteSQL(SQL, dialect='SQLITE') 行でクラッシュします。GDAL ライブラリを更新しようとしましたが、pip vcvarsall.bat の問題が発生しています (これも非常にイライラします)。どんな助けでも大歓迎です！

from osgeo import ogr
ogr.UseExceptions()
ogr_ds = ogr.Open(r'D:\inputs', True) 
SQL = """\
SELECT ST_Intersection(A.geometry, B.geometry) AS geometry, A.*, B.*
FROM stops A, zones B
WHERE ST_Intersects(A.geometry, B.geometry);
"""
layer = ogr_ds.ExecuteSQL(SQL, dialect='SQLITE')
# copy result back to datasource as a new shapefile
layer2 = ogr_ds.CopyLayer(layer, 'result')
# save, close
layer = layer2 = ogr_ds = None