Tách một tính năng khi giao với một tính năng của lớp khác bằng PyQGIS / Python?

Tôi có một lớp đệm (đa giác màu xanh lá cây) mà tôi muốn chia thành hai đa giác bất cứ khi nào nó vượt qua một rào cản (đường màu xanh). Tôi đã cố gắng sử dụng phương pháp "splitGeometry", nhưng tôi không thể làm cho nó hoạt động được. Mã của tôi cho đến nay là:

while ldbuffprovider.nextFeature(feat):
  while barprovider.nextFeature(feat2):
    if feat.geometry().intersects(feat2.geometry()):
        intersection = feat.geometry().intersection(feat2.geometry())
        result, newGeometries, topoTestPoints=feat.geometry().splitGeometry(intersection.asPolyline(),True) 

Trả về 1 cho kết quả (lỗi) và danh sách trống cho newGeometries. Bất kỳ trợ giúp nào cũng được đánh giá rất cao.

Bạn có thể sử dụng reshapeGeometrychức năng của QgsGeometryđối tượng cho việc này, nó cắt một đa giác dọc theo giao điểm của nó với một đường.

Sau đây sẽ giao cắt các đa giác bộ đệm với các dòng và thêm các tính năng đa giác phân chia vào một lớp bộ nhớ (cú pháp QGIS 2.0):

# Get the dataProvider objects for the layers called 'line' and 'buffer'
linepr = QgsMapLayerRegistry.instance().mapLayersByName('line')[0].dataProvider()
bufferpr = QgsMapLayerRegistry.instance().mapLayersByName('buffer')[0].dataProvider()

# Create a memory layer to store the result
resultl = QgsVectorLayer("Polygon", "result", "memory")
resultpr = resultl.dataProvider()
QgsMapLayerRegistry.instance().addMapLayer(resultl)


for feature in bufferpr.getFeatures():
  # Save the original geometry
  geometry = QgsGeometry.fromPolygon(feature.geometry().asPolygon())
  for line in linepr.getFeatures():
    # Intersect the polygon with the line. If they intersect, the feature will contain one half of the split
    t = feature.geometry().reshapeGeometry(line.geometry().asPolyline())
    if (t==0):
      # Create a new feature to hold the other half of the split
      diff = QgsFeature()
      # Calculate the difference between the original geometry and the first half of the split
      diff.setGeometry( geometry.difference(feature.geometry()))
      # Add the two halves of the split to the memory layer
      resultpr.addFeatures([feature])
      resultpr.addFeatures([diff])

Một xấp xỉ tốt với GDAL> = 1.10.0 được biên dịch bằng SQLite và SpatiaLite bao gồm việc bọc các lớp của bạn (ví dụ: poligon.shp và line.shp ) trong tệp OGR VRT (ví dụ: layer.vrt ):

<OGRVRTDataSource>
    <OGRVRTlayer name="buffer_line">
        <SrcDataSource>line.shp</SrcDataSource>
        <SrcSQL dialect="sqlite">SELECT ST_Buffer(geometry,0.000001) from line</SrcSQL>
    </OGRVRTlayer>
    <OGRVRTlayer name="polygon">
        <SrcDataSource>polygon.shp</SrcDataSource>
    </OGRVRTlayer>
</OGRVRTDataSource>

để có một bộ đệm rất nhỏ (ví dụ 1 micron) xung quanh line.shp có được lớp * buffer_line *. Sau đó, chúng ta có thể áp dụng sự khác biệt và khác biệt đối xứng trên các hình học này bằng SpatiaLite:

ogr2ogr splitted_polygons.shp layers.vrt -dialect sqlite -sql "SELECT ST_Difference(ST_SymDifference(g1.geometry,g2.geometry),g2.geometry) FROM polygon AS g1, buffer_line AS g2" -explodecollections

Rõ ràng, tất cả những thứ này hoàn toàn có thể thực hiện được từ tập lệnh Python:

os.system("some_command with args")

Hi vọng điêu nay co ich!

Câu trả lời của Jake có vấn đề tô pô với đa giác bị tách. Tôi đã thực hiện một giải pháp, hy vọng điều này có thể giúp: https://gist.github.com/RamonLopezEscudero/844c1401f5339143da1b2b5cf7ff27bd