Merge pull request #2516 from radosuav/zonal_stats

[processing] use raster band scailing factor and offset in Zonal Statistics algorithm
qgis · Dec 1, 2015 · 17bd46e · 17bd46e
2 parents bb9eb79 + 34df72c
commit 17bd46e
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/python/plugins/processing/algs/qgis/ZonalStatistics.py b/python/plugins/processing/algs/qgis/ZonalStatistics.py
@@ -118,6 +118,7 @@ def processAlgorithm(self, progress):
 
             srcOffset = (startColumn, startRow, width, height)
             srcArray = rasterBand.ReadAsArray(*srcOffset)
+            srcArray = srcArray*rasterBand.GetScale()+rasterBand.GetOffset()
 
             newGeoTransform = (
                 geoTransform[0] + srcOffset[0] * geoTransform[1],
@@ -192,6 +193,7 @@ def processAlgorithm(self, progress):
 
                 srcOffset = (startColumn, startRow, width, height)
                 srcArray = rasterBand.ReadAsArray(*srcOffset)
+                srcArray = srcArray*rasterBand.GetScale()+rasterBand.GetOffset()
 
                 newGeoTransform = (
                     geoTransform[0] + srcOffset[0] * geoTransform[1],