Python gdal在使用ogr in、Contains或oth时停止工作

Question

更新：通过进一步调查，我发现这一定是由某种损坏的几何结构引起的。但即使我跑了：

if drillhole[1].IsValid():

它导致了车祸。所以不知道怎么检查。在

我试着测试线条是否在多边形内。使用shapely它对我很好，除了速度-我有一万条线，多边形也是一样。我只是想测试一下ogr是否可以更快地完成这个任务，但是没有运气让它工作。在

^{pr2}$

知道线路有什么问题吗？公司名称：

if cell[1].Contains(drillhole[1]):

在我的两台Windows7电脑上python.exe停止工作。。。 问题签名：

  Problem Event Name:   APPCRASH
  Application Name: python.exe
  Application Version:  0.0.0.0
  Application Timestamp:    5193f3af
  Fault Module Name:    gdal111.dll
  Fault Module Version: 1.11.2.0
  Fault Module Timestamp:   54e65215
  Exception Code:   c0000005
  Exception Offset: 00000000005e5fb3
  OS Version:   6.1.7601.2.1.0.768.3

或者有没有其他更快的方法或者方法来提高shapely的使用速度？在

Answer 1

这次崩溃很可能是一个有充分记录的gotcha。在

类似的工作也可以在Shapely身上完成，而Shapely没有这个问题。另外，for循环中的查询可以通过使用Rtree索引更快地完成，这里有一些与shapely+Rtree相关的问答。在