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

Use sample instead of identify in raster sampling alg

Browse Source 
About 25x faster on large layers
This commit is contained in:
Nyall Dawson 2018-07-16 21:12:02 +10:00
parent 0f41ca9e13
commit 8b2d5d0374
1 changed files with 8 additions and 15 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -31,7 +31,8 @@ import os
from qgis.PyQt.QtGui import QIcon
from qgis.PyQt.QtCore import QVariant
from qgis.core import (QgsApplication,
from qgis.core import (NULL,
QgsApplication,
QgsField,
QgsFeatureSink,
QgsRaster,
@ -179,20 +180,12 @@ class RasterSampling(QgisAlgorithm):
feedback.reportError(self.tr('Could not reproject feature {} to raster CRS').format(i.id()))
continue
if sampled_raster.bandCount() > 1:
for b in range(sampled_raster.bandCount()):
attrs.append(
sampled_raster.dataProvider().identify(
point,
QgsRaster.IdentifyFormatValue).results()[b + 1]
)
attrs.append(
sampled_raster.dataProvider().identify(
point,
QgsRaster.IdentifyFormatValue).results()[1]
)
for b in range(sampled_raster.bandCount()):
value, ok = sampled_raster.dataProvider().sample(point, b + 1)
if ok:
attrs.append(value)
else:
attrs.append(NULL)
i.setAttributes(attrs)