Guest User

Untitled

a guest
Nov 23rd, 2017
60
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.07 KB | None | 0 0
  1. def splitting():
  2. selectedFeatures = []
  3. for l_feat in l_layer.getFeatures():
  4. l_geom = l_feat.geometry().asPolyline()
  5. l_featid = l_feat.id()
  6. print('l_geom is',l_geom)
  7. g_layer.select(l_featid)
  8. sel_feat = g_layer.selectedFeatures()[0]
  9. sel_feat_geom = sel_feat.geometry()
  10. for p_feat in p_layer.getFeatures():
  11. if p_feat.geometry().intersects(sel_feat_geom):
  12. intersection_geom = p_feat.geometry().intersection(sel_feat_geom)
  13. if intersection_geom.asPoint() != l_geom[0] and intersection_geom.asPoint() !=l_geom[-1]:
  14. print '------', p_feat.id()
  15. l_layer.startEditing()
  16. # l_layer.select(l_feat.id())
  17. t = l_layer.splitFeatures([intersection_geom.asPoint()],True)
  18. print('t is',t)
  19. l_layer.commitChanges()
  20. # l_layer.removeSelection()
  21. else:
  22. pass
  23. g_layer.removeSelection()
  24. qgis.utils.iface.mapCanvas().refresh()
Add Comment
Please, Sign In to add comment