Skip to content
Permalink
Browse files

refactor Geometry convert alg

  • Loading branch information
alexbruy committed Oct 3, 2014
1 parent 8e696b7 commit be0798961992492d3f285861d7d15b873fea1b09
@@ -83,6 +83,7 @@
from mmqgisx.HubDistance import HubDistance
from mmqgisx.HubLines import HubLines
from mmqgisx.Merge import Merge
from mmqgisx.GeometryConvert import GeometryConvert

from ConcaveHull import ConcaveHull
from Polygonize import Polygonize
@@ -150,7 +151,7 @@ def __init__(self):
DeleteColumn(), DeleteDuplicateGeometries(),
TextToFloat(), ExtractByAttribute(),
SelectByAttribute(), Grid(), Gridify(), HubDistance(),
HubLines(), Merge(),
HubLines(), Merge(), GeometryConvert(),
#mmqgisx_delete_duplicate_geometries_algorithm(),
#mmqgisx_geometry_convert_algorithm(),
#mmqgisx_grid_algorithm(),
@@ -0,0 +1,221 @@
# -*- coding: utf-8 -*-

"""
***************************************************************************
Gridify.py
---------------------
Date : May 2010
Copyright : (C) 2010 by Michael Minn
Email : pyqgis at michaelminn dot com
***************************************************************************
* *
* This program is free software; you can redistribute it and/or modify *
* it under the terms of the GNU General Public License as published by *
* the Free Software Foundation; either version 2 of the License, or *
* (at your option) any later version. *
* *
***************************************************************************
"""

__author__ = 'Michael Minn'
__date__ = 'May 2010'
__copyright__ = '(C) 2010, Michael Minn'

# This will get replaced with a git SHA1 when you do a git archive

__revision__ = '$Format:%H$'

from PyQt4.QtCore import *
from qgis.core import *
from processing.core.GeoAlgorithm import GeoAlgorithm
from processing.core.GeoAlgorithmExecutionException import \
GeoAlgorithmExecutionException
from processing.core.parameters import ParameterVector
from processing.core.parameters import ParameterSelection
from processing.core.outputs import OutputVector

from processing.tools import dataobjects, vector

class GeometryConvert(GeoAlgorithm):
INPUT = 'INPUT'
TYPE = 'TYPE'
OUTPUT = 'OUTPUT'

TYPES = ['Centroids',
'Nodes',
'Linestrings',
'Multilinestrings',
'Polygons'
]

def defineCharacteristics(self):
self.name = 'Convert geometry type'
self.group = 'Vector geometry tools'

self.addParameter(ParameterVector(self.INPUT,
'Input layer', [ParameterVector.VECTOR_TYPE_ANY]))
self.addParameter(ParameterSelection(self.TYPE,
'New geometry type', self.TYPES))

self.addOutput(OutputVector(self.OUTPUT, 'Output'))

def processAlgorithm(self, progress):
layer = dataobjects.getObjectFromUri(
self.getParameterValue(self.INPUT))
index = self.getParameterValue(self.TYPE)

splitNodes = False
if index == 0:
newType = QGis.WKBPoint
elif index == 1:
newType = QGis.WKBPoint
splitNodes = True
elif index == 2:
newType = QGis.WKBLineString
elif index == 3:
newType = QGis.WKBMultiLineString
elif index == 4:
newType = QGis.WKBPolygon
else:
newType = QGis.WKBPoint

writer = self.getOutputFromName(self.OUTPUT).getVectorWriter(
layer.pendingFields(), newType, layer.crs())

features = vector.features(layer)
count = len(features)
total = 100.0 / float(count)

for count, f in enumerate(features):
geom = f.geometry()
geomType = geom.wkbType()

if geomType in [QGis.WKBPoint, QGis.WKBPoint25D]:
if newType == QGis.WKBPoint:
writer.addFeature(f)
else:
raise GeoAlgorithmExecutionException(
'Cannot convert from %s to %s', geomType, newType)
elif geomType in [QGis.WKBMultiPoint, QGis.WKBMultiPoint25D]:
if newType == QGis.WKBPoint and splitNodes:
points = geom.asMultiPoint()
for p in points:
feat = QgsFeature()
feat.setAttributes(f.attributes())
feat.setGeometry(QgsGeometry.fromPoint(p))
writer.addFeature(feat)
elif newType == QGis.WKBPoint:
feat = QgsFeature()
feat.setAttributes(f.attributes())
feat.setGeometry(geom.centroid())
writer.addFeature(feat)
else:
raise GeoAlgorithmExecutionException(
'Cannot convert from %s to %s', geomType, newType)
elif geomType in [QGis.WKBLineString, QGis.WKBLineString25D]:
if newType == QGis.WKBPoint and splitNodes:
points = geom.asPolyline()
for p in points:
feat = QgsFeature()
feat.setAttributes(f.attributes())
feat.setGeometry(QgsGeometry.fromPoint(p))
writer.addFeature(feat)
elif newType == QGis.WKBPoint:
feat = QgsFeature()
feat.setAttributes(f.attributes())
feat.setGeometry(geom.centroid())
writer.addFeature(feat)
elif newType == QGis.WKBLineString:
writer.addFeature(f)
else:
raise GeoAlgorithmExecutionException(
'Cannot convert from %s to %s', geomType, newType)
elif geomType in [QGis.WKBMultiLineString, QGis.WKBMultiLineString25D]:
if newType == QGis.WKBPoint and splitNodes:
lines = geom.asMultiPolyline()
for line in lines:
for p in line:
feat = QgsFeature()
feat.setAttributes(f.attributes())
feat.setGeometry(QgsGeometry.fromPoint(p))
writer.addFeature(feat)
elif newType == QGis.WKBPoint:
feat = QgsFeature()
feat.setAttributes(f.attributes())
feat.setGeometry(geom.centroid())
writer.addFeature(feat)
elif newType == QGis.WKBLineString:
lines = geom.asMultiPolyline()
for line in lines:
feat = QgsFeature()
feat.setAttributes(f.attributes())
feat.setGeometry(QgsGeometry.fromPolyline(line))
writer.addFeature(feat)
elif newType == QGis.WKBMultiLineString:
writer.addFeature(f)
else:
raise GeoAlgorithmExecutionException(
'Cannot convert from %s to %s', geomType, newType)
elif geomType in [QGis.WKBPolygon, QGis.WKBPolygon25D]:
if newType == QGis.WKBPoint and splitNodes:
rings = geom.asPolygon()
for ring in rings:
for p in ring:
feat = QgsFeature()
feat.setAttributes(f.attributes())
feat.setGeometry(QgsGeometry.fromPoint(p))
writer.addFeature(feat)
elif newType == QGis.WKBPoint:
feat = QgsFeature()
feat.setAttributes(f.attributes())
feat.setGeometry(geom.centroid())
writer.addFeature(feat)
elif newType == QGis.WKBMultiLineString:
rings = geom.asPolygon()
feat = QgsFeature()
feat.setAttributes(f.attributes())
feat.setGeometry(QgsGeometry.fromMultiPolyline(rings))
writer.addFeature(feat)
elif newtype == QGis.WKBPolygon:
writer.addFeature(f)
else:
raise GeoAlgorithmExecutionException(
'Cannot convert from %s to %s', geomType, newType)
elif geomType in [QGis.WKBMultiPolygon, QGis.WKBMultiPolygon25D]:
if newType == QGis.WKBPoint and splitNodes:
polygons = geom.asMultiPolygon()
for polygon in polygons:
for line in polygon:
for p in line:
feat = QgsFeature()
feat.setAttributes(f.attributes())
feat.setGeometry(QgsGeometry.fromPoint(p))
writer.addFeature(feat)
elif newType == QGis.WKBPoint:
feat = QgsFeature()
feat.setAttributes(f.attributes())
feat.setGeometry(geom.centroid())
writer.addFeature(feat)
elif newType == QGis.WKBLineString:
polygons = geom.asMultiPolygon()
for polygons in polygons:
feat = QgsFeature()
feat.setAttributes(f.attributes())
feat.setGeometry(QgsGeometry.fromPolyline(polygon))
writer.addFeature(feat)
elif newType == QGis.WKBPolygon:
polugons = geom.asMultiPolygon()
for polygons in polygons:
feat = QgsFeature()
feat.setAttributes(f.attributes())
feat.setGeometry(QgsGeometry.fromPolygon(polygon))
writer.addFeature(feat)
elif newType in [QGis.WKBMultiLineString, QGis.WKBMultiPolygon]:
writer.addFeature(f)
else:
raise GeoAlgorithmExecutionException(
'Cannot convert from %s to %s', geomType, newType)

progress.setPercentage(int(count * total))

del writer

0 comments on commit be07989

Please sign in to comment.
You can’t perform that action at this time.