Advertisement
Guest User

Untitled

a guest
Jan 22nd, 2018
62
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.02 KB | None | 0 0
  1. ##skropty=group
  2. ##Warstwa1=vector polygon
  3. ##Wyjscie=output vector
  4.  
  5. from qgis.core import *
  6. from PyQt4.QtCore import QVariant
  7.  
  8. warstwa1 = processing.getObject(Warstwa1)
  9.  
  10. pola = QgsFields(warstwa1.fields())
  11. powr = QgsField('Powierzchnia', QVariant.Double)
  12. pola.append(powr)
  13.  
  14.  
  15. writer = processing.VectorWriter(
  16. Wyjscie,
  17. warstwa1.dataProvider().encoding(),
  18. pola,
  19. QGis.WKBPolygon,
  20. warstwa1.crs())
  21.  
  22. feats = processing.features(warstwa1)
  23. n = len(feats)
  24. for i, feat in enumerate(feats):
  25. progress.setPercentage(int(100 * i / n))
  26. obiekt = QgsFeature(pola)
  27. geometria = feat.geometry()
  28. obiekt.setGeometry(geometria)
  29. gobiektu = obiekt.geometry()
  30. obiekt['ID'] = feat['ID']
  31. obiekt['KOD'] = feat['KOD']
  32. obiekt['MEZOREGION'] = feat['MEZOREGION']
  33. obiekt['MAKROREGIO'] = feat['MAKROREGIO']
  34. obiekt['PODPROWINC'] = feat['PODPROWINC']
  35. obiekt['PROWINCJA'] = feat['PROWINCJA']
  36. obiekt['Powierzchnia'] = gobiektu.area()
  37. writer.addFeature(obiekt)
  38. del writer
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement