Advertisement
Guest User

Untitled

a guest
Jun 28th, 2016
57
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.59 KB | None | 0 0
  1. def get_subsel(layer,filter):
  2. layers = iface.mapCanvas().layers()
  3. for lyr in layers:
  4. if lyr.name() == layer:
  5. request = QgsFeatureRequest().setFilterExpression(filter)
  6.  
  7. iter = lyr.getFeatures(request)
  8. ids = [i.id() for i in iter]
  9. lyr.setSelectedFeatures(ids)
  10.  
  11. selected_feats = lyr.selectedFeatures()
  12.  
  13. new_layer = QgsVectorLayer("Point?crs=epsg:277000", "acc_subselection", "memory")
  14. data_prov = new_layer.dataProvider()
  15. data_prov.addFeatures(selected_feats)
  16. QgsMapLayerRegistry().instance().addMapLayer(new_layer)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement