我正在尝试实现这个问题的答案： https://gis.stackexchange.com/questions/119374/intersect-shapefiles-using-shapely
但我得到错误：
RuntimeError：使用了未定义的函数'ST_Intersection'。
这是代码。它在layer = ogr_ds.ExecuteSQL（SQL，dialect ='SQLITE'）行崩溃。我尝试更新GDAL库，但我遇到了一个点vcvarsall.bat问题（也非常令人沮丧）。非常感谢任何帮助！
string答案 0 :(得分：0)
我认为错误的产生是因为您需要使用SQLite SQL方法编译SQLite驱动程序。请参阅：http://www.gdal.org/ogr_sql_sqlite.html
至于更新和编译gdal，这可能很烦人。你在运行什么操作系统？