sharpetronics/QGIS
1
0
Fork 0
mirror of https://github.com/qgis/QGIS.git synced 2025-02-27 00:33:48 -05:00
Code Issues Packages Projects Releases Wiki Activity

[processing] be more flexible with multipoint layers, allow sampling for single part point

Browse Source
This commit is contained in:
nirvn 2018-07-28 10:49:06 +07:00 committed by Mathieu Pellerin
parent 3e31a6570c
commit e19e1dea9d
1 changed files with 11 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

16
python/plugins/processing/algs/qgis/RasterSampling.py
View File

@ -36,6 +36,7 @@ from qgis.core import (NULL,
QgsField,
QgsFeatureSink,
QgsRaster,
QgsPointXY,
QgsProcessing,
QgsProcessingParameterRasterLayer,
QgsProcessingParameterString,
@ -160,13 +161,18 @@ class RasterSampling(QgisAlgorithm):
attrs = i.attributes()
if i.geometry().isMultipart():
raise QgsProcessingException(self.tr('''Impossible to sample data
of a Multipart layer. Please use the Multipart to single part
algorithm to transform the layer.'''))
if i.geometry().isMultipart() and i.geometry().constGet().partCount() > 1:
sink.addFeature(i, QgsFeatureSink.FastInsert)
feedback.setProgress(int(n * total))
feedback.reportError(self.tr('Impossible to sample data of multipart feature {}.').format(i.id()))
continue
# get the feature geometry as point
point = i.geometry().asPoint()
point = QgsPointXY()
if i.geometry().isMultipart():
point = i.geometry().asMultiPoint()[0]
else:
point = i.geometry().asPoint()
# reproject to raster crs
try: