QGIS 3.8 Python Code for calculate area

Asked Oct 06 '19 at 14:39

Active Oct 06 '19 at 21:52

Viewed 368 times

Is this the the best solution to calculate the area of the features?

parent = iface.mainWindow()
mc = iface.mapCanvas()
fldname = "m2" #field which has the values for squaremeters
myLyr = mc.currentLayer() 
d = QgsDistanceArea()
d.setSourceCrs(myLyr.crs(), QgsProject.instance().transformContext())
d.setEllipsoid(QgsProject.instance().ellipsoid())
#https://gis.stackexchange.com/questions/309721/creating-and-updating-a-field-from-calculation?rq=1
#https://gis.stackexchange.com/questions/289601/how-to-fill-fields-with-layer-name-in-pyqgis
with edit(myLyr):
   fldindex = myLyr.fields().indexOf(fldname)
   if fldindex == -1:
       print("Das gesuchte Feld {} ist nicht im Layer {}!".format(fldname, myLyr.name()))
   else:
       fts = myLyr.getFeatures()
       for f in fts:
           geom = f.geometry()
           v=(d.measureArea(geom))
           #print("Fläche in qm ist {}!".format(v))
           f[fldindex]=v
           myLyr.updateFeature(f)

edited Oct 06 '19 at 21:52

PolyGeo

65,136
29
109
338

asked Oct 06 '19 at 14:39

klaus k.

2

For review of working code there is the [codereview.se] Stack Exchange. – PolyGeo Oct 06 '19 at 21:53
1

Looks fine, are you having issues? – Joseph Oct 07 '19 at 10:59
No there aren't issue. My hope was to make the code shorter – klaus k. Oct 07 '19 at 19:25
Years later, this still helped me getting the right area (like the expression $area would). Thank you! – Homerun_ Sep 21 '23 at 09:32

QGIS 3.8 Python Code for calculate area

0 Answers0