0

更新：通过进一步的调查，我发现它一定是由某种损坏的几何图形引起的。但即使我运行：

if drillhole[1].IsValid():

它正在导致崩溃。所以不知道如何检查。

我正在尝试测试线条是否在多边形内。除了速度之外，使用 shapely 对我来说效果很好——我有一万条线，多边形也是如此。我只是想测试一下 ogr 是否可以更快地做到这一点，但没有运气让它发挥作用。

import ogr
#Load cells as polygons with slight buffer
data_source = ogr.Open(file_path_cells)
source_layer = data_source.GetLayer()
source_layer.ResetReading()
cells = []
for source_feature in source_layer:
  feature_id = source_feature.GetFieldAsInteger(0)
  feature_geometry = source_feature.geometry()
  feature_geometry = feature_geometry.Buffer(0.1, quadsecs = 3)
  cells.append((feature_id,feature_geometry))

#Load drillholes as lines wihtin cells
data_source = ogr.Open(file_path_drillholes)
source_layer = data_source.GetLayer()
source_layer.ResetReading()
drillholes = []
for source_feature in source_layer:
  feature_elev = source_feature.GetFieldAsInteger("elev")
  feature_geometry = source_feature.geometry()
  drillholes.append((feature_elev,feature_geometry))  
for cell in cells:
  for drillhole in drillholes:
    if cell[1].Contains(drillhole[1]):
      print("yes")
    else:
      print("no")

知道这条线有什么问题吗？：

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

在我的两台 Windows 7 机器上，我总是让 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

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

4

1 回答 1

2

崩溃很可能是一个有据可查的陷阱。

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

于 2015-07-05T21:21:45.253 回答