I want to find the bounding box of several features selected by their attributes. Here's what I've been trying:
tNumList = [8,9,10,11]
selectGeom = QgsGeometry()
for tNum in tNumList:
request = QgsFeatureRequest().setFilterExpression( u'"Territory" = {}'.format(tNum))
iter = territoryBorders.getFeatures(request)
feat = QgsFeature()
iter.nextFeature(feat)
selectGeom.combine(feat.geometry())
print selectGeom.boundingBox().width()
Is there a better way to get a single feature from a layer then using the getFeatures(request)?
This just gives me empty rectangles. What am I doing wrong?
Best Answer
I've optimized and corrected your code.
Optimized code: