mirror of
https://github.com/qgis/QGIS.git
synced 2025-03-06 00:05:02 -05:00
06c4b07222
Since the layer has more complete knowledge of the crs (ie, when provider could not determine crs and user has selected it from the list), and also better knowledge of layer fields (virtual fields, joined fields) we should use these rather than the provider methods.
153 lines
6.3 KiB
Python
153 lines
6.3 KiB
Python
# -*- coding: utf-8 -*-
|
|
|
|
"""
|
|
***************************************************************************
|
|
SymmetricalDifference.py
|
|
---------------------
|
|
Date : September 2014
|
|
Copyright : (C) 2014 by Alexander Bruy
|
|
Email : alexander dot bruy at gmail 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__ = 'Alexander Bruy'
|
|
__date__ = 'September 2014'
|
|
__copyright__ = '(C) 2014, Alexander Bruy'
|
|
|
|
# This will get replaced with a git SHA1 when you do a git archive
|
|
|
|
__revision__ = '$Format:%H$'
|
|
|
|
import os
|
|
|
|
from qgis.PyQt.QtGui import QIcon
|
|
|
|
from qgis.core import Qgis, QgsFeature, QgsGeometry, QgsFeatureRequest, NULL, QgsWkbTypes
|
|
from processing.core.ProcessingLog import ProcessingLog
|
|
from processing.core.GeoAlgorithm import GeoAlgorithm
|
|
from processing.core.GeoAlgorithmExecutionException import GeoAlgorithmExecutionException
|
|
from processing.core.parameters import ParameterVector
|
|
from processing.core.outputs import OutputVector
|
|
from processing.tools import dataobjects, vector
|
|
|
|
pluginPath = os.path.split(os.path.split(os.path.dirname(__file__))[0])[0]
|
|
|
|
|
|
class SymmetricalDifference(GeoAlgorithm):
|
|
|
|
INPUT = 'INPUT'
|
|
OVERLAY = 'OVERLAY'
|
|
OUTPUT = 'OUTPUT'
|
|
|
|
def getIcon(self):
|
|
return QIcon(os.path.join(pluginPath, 'images', 'ftools', 'sym_difference.png'))
|
|
|
|
def defineCharacteristics(self):
|
|
self.name, self.i18n_name = self.trAlgorithm('Symmetrical difference')
|
|
self.group, self.i18n_group = self.trAlgorithm('Vector overlay tools')
|
|
self.addParameter(ParameterVector(self.INPUT,
|
|
self.tr('Input layer'), [ParameterVector.VECTOR_TYPE_ANY]))
|
|
self.addParameter(ParameterVector(self.OVERLAY,
|
|
self.tr('Difference layer'), [ParameterVector.VECTOR_TYPE_ANY]))
|
|
self.addOutput(OutputVector(self.OUTPUT,
|
|
self.tr('Symmetrical difference')))
|
|
|
|
def processAlgorithm(self, progress):
|
|
layerA = dataobjects.getObjectFromUri(
|
|
self.getParameterValue(self.INPUT))
|
|
layerB = dataobjects.getObjectFromUri(
|
|
self.getParameterValue(self.OVERLAY))
|
|
|
|
geomType = layerA.wkbType()
|
|
fields = vector.combineVectorFields(layerA, layerB)
|
|
writer = self.getOutputFromName(self.OUTPUT).getVectorWriter(
|
|
fields, geomType, layerA.crs())
|
|
|
|
featB = QgsFeature()
|
|
outFeat = QgsFeature()
|
|
|
|
indexA = vector.spatialindex(layerB)
|
|
indexB = vector.spatialindex(layerA)
|
|
|
|
featuresA = vector.features(layerA)
|
|
featuresB = vector.features(layerB)
|
|
|
|
total = 100.0 / (len(featuresA) * len(featuresB))
|
|
count = 0
|
|
|
|
for featA in featuresA:
|
|
add = True
|
|
geom = featA.geometry()
|
|
diffGeom = QgsGeometry(geom)
|
|
attrs = featA.attributes()
|
|
intersects = indexA.intersects(geom.boundingBox())
|
|
for i in intersects:
|
|
layerB.getFeatures(QgsFeatureRequest().setFilterFid(i)).nextFeature(featB)
|
|
tmpGeom = featB.geometry()
|
|
if diffGeom.intersects(tmpGeom):
|
|
diffGeom = QgsGeometry(diffGeom.difference(tmpGeom))
|
|
if not diffGeom.isGeosValid():
|
|
ProcessingLog.addToLog(ProcessingLog.LOG_ERROR,
|
|
self.tr('GEOS geoprocessing error: One or '
|
|
'more input features have invalid '
|
|
'geometry.'))
|
|
add = False
|
|
break
|
|
|
|
if add:
|
|
try:
|
|
outFeat.setGeometry(diffGeom)
|
|
outFeat.setAttributes(attrs)
|
|
writer.addFeature(outFeat)
|
|
except:
|
|
ProcessingLog.addToLog(ProcessingLog.LOG_WARNING,
|
|
self.tr('Feature geometry error: One or more output features ignored due to invalid geometry.'))
|
|
continue
|
|
|
|
count += 1
|
|
progress.setPercentage(int(count * total))
|
|
|
|
length = len(layerA.fields())
|
|
|
|
for featA in featuresB:
|
|
add = True
|
|
geom = featA.geometry()
|
|
diffGeom = QgsGeometry(geom)
|
|
attrs = featA.attributes()
|
|
attrs = [NULL] * length + attrs
|
|
intersects = indexB.intersects(geom.boundingBox())
|
|
for i in intersects:
|
|
layerA.getFeatures(QgsFeatureRequest().setFilterFid(i)).nextFeature(featB)
|
|
tmpGeom = featB.geometry()
|
|
if diffGeom.intersects(tmpGeom):
|
|
diffGeom = QgsGeometry(diffGeom.difference(tmpGeom))
|
|
if not diffGeom.isGeosValid():
|
|
ProcessingLog.addToLog(ProcessingLog.LOG_ERROR,
|
|
self.tr('GEOS geoprocessing error: One or '
|
|
'more input features have invalid '
|
|
'geometry.'))
|
|
add = False
|
|
break
|
|
|
|
if add:
|
|
try:
|
|
outFeat.setGeometry(diffGeom)
|
|
outFeat.setAttributes(attrs)
|
|
writer.addFeature(outFeat)
|
|
except:
|
|
ProcessingLog.addToLog(ProcessingLog.LOG_WARNING,
|
|
self.tr('Feature geometry error: One or more output features ignored due to invalid geometry.'))
|
|
continue
|
|
|
|
count += 1
|
|
progress.setPercentage(int(count * total))
|
|
|
|
del writer
|