Advertisement
Guest User

Untitled

a guest
Sep 23rd, 2017
61
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.83 KB | None | 0 0
  1. layer = iface.activeLayer()
  2. areatot = 0
  3. for elem in layer.getFeatures()
  4. geom = elem.geometry()
  5. areatot += geom.area()
  6.  
  7. aeratot = [sum(elem.geometry().area() for elem in layer.getFeatures())]
  8.  
  9. # creation of a empty geometry for unioning
  10. geomtot = QgsGeometry.fromWkt('GEOMETRYCOLLECTION EMPTY')
  11. # union
  12. for elem in layer.getFeatures():
  13. geomtot = geomtot.combine(elem.geometry())
  14. area = geomtot.area()
  15.  
  16. layer = qgis.utils.iface.activeLayer()
  17. areatot = 0
  18. layer.select()
  19. for elem in layer:
  20. geom = elem.geometry()
  21. areatot += geom.area()
  22.  
  23. geomtot = QgsGeometry.fromWkt('GEOMETRYCOLLECTION EMPTY')
  24. layer.select()
  25. for elem in layer:
  26. geomtot = geomtot.combine(elem.geometry())
  27. area = geomtot.area()
  28.  
  29. vlayer = QgsVectorLayer("/path/to/shapefile/file.shp", "layer_name_you_like", "ogr")
  30.  
  31. vlayer.extent().toString()
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement