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

Answer 1

我认为出现错误是因为您需要使用 SQLite SQL 方言编译 SQLite 驱动程序。见：http://www.gdal.org/ogr_sql_sqlite.html

至于更新和编译 gdal，那可能很烦人。你运行的是什么操作系统？