Advertisement
Guest User

Untitled

a guest
Aug 19th, 2019
98
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.11 KB | None | 0 0
  1. roadlist=[]
  2. layer = iface.activeLayer()
  3. chainage_layer=QgsVectorLayer(r'Z:My DriveMangoesmappingSpatial Projects2019DSC126_Road_UpdatesWorkingOutputDSC_Roads_10m_Chainage.tab', "Chainage", "ogr")
  4. for feature in layer.getFeatures():
  5. if feature['RoadMntnc'] not in roadlist:
  6. roadlist.append(feature['RoadMntnc'])
  7.  
  8. for road in roadlist:
  9. print (road)
  10. layer.selectByExpression('"RoadMntnc"= '{0}''.format(road), QgsVectorLayer.SetSelection)
  11. selection=layer.selectedFeatures()
  12. for feature in selection:
  13. chainage_layer.selectByExpression('"RoadMntnc"= '{0}''.format(road), QgsVectorLayer.SetSelection)
  14. chainage_selection=layer.selectedFeatures()
  15. if start_point(buffer(feature.geometry,5,10)).intersects(chainage_selection.geometry()):
  16. print (feature2['distance'])
  17. #aggregate(layer:='DSC_Roads_10m_Chainage.tab', aggregate:='max', expression:="distance" * (attribute(@parent, 'RoadMntnc') = attribute($currentfeature, 'RoadMntnc')),filter:= within(start_point(geometry(@parent)),buffer($geometry,1,10))) # process in model builder -field calculator.
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement