Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- def splitting():
- selectedFeatures = []
- for l_feat in l_layer.getFeatures():
- l_geom = l_feat.geometry().asPolyline()
- l_featid = l_feat.id()
- print('l_geom is',l_geom)
- g_layer.select(l_featid)
- sel_feat = g_layer.selectedFeatures()[0]
- sel_feat_geom = sel_feat.geometry()
- for p_feat in p_layer.getFeatures():
- if p_feat.geometry().intersects(sel_feat_geom):
- intersection_geom = p_feat.geometry().intersection(sel_feat_geom)
- if intersection_geom.asPoint() != l_geom[0] and intersection_geom.asPoint() !=l_geom[-1]:
- print '------', p_feat.id()
- l_layer.startEditing()
- # l_layer.select(l_feat.id())
- t = l_layer.splitFeatures([intersection_geom.asPoint()],True)
- print('t is',t)
- l_layer.commitChanges()
- # l_layer.removeSelection()
- else:
- pass
- g_layer.removeSelection()
- qgis.utils.iface.mapCanvas().refresh()
Add Comment
Please, Sign In to add comment