@@ -1379,7 +1379,7 @@ QgsRasterHistogram QgsGdalProvider::histogram( int theBandNo,
13791379
13801380 // Min/max, if not specified, are set by histogramDefaults, it does not
13811381 // set however min/max shifted to avoid rounding errors
1382-
1382+
13831383 double myMinVal = myHistogram.minimum ;
13841384 double myMaxVal = myHistogram.maximum ;
13851385
@@ -1388,8 +1388,8 @@ QgsRasterHistogram QgsGdalProvider::histogram( int theBandNo,
13881388 double myOffset = bandOffset ( theBandNo );
13891389 if ( myScale != 1.0 || myOffset != 0 . )
13901390 {
1391- myMinVal = (myHistogram.minimum - myOffset) / myScale;
1392- myMaxVal = (myHistogram.maximum - myOffset) / myScale;
1391+ myMinVal = ( myHistogram.minimum - myOffset ) / myScale;
1392+ myMaxVal = ( myHistogram.maximum - myOffset ) / myScale;
13931393 }
13941394
13951395 double dfHalfBucket = ( myMaxVal - myMinVal ) / ( 2 * myHistogram.binCount );
@@ -2385,34 +2385,34 @@ QgsRasterBandStats QgsGdalProvider::bandStatistics( int theBandNo, int theStats,
23852385 myRasterBandStats.statsGathered = QgsRasterBandStats::Min | QgsRasterBandStats::Max
23862386 | QgsRasterBandStats::Range | QgsRasterBandStats::Mean
23872387 | QgsRasterBandStats::StdDev;
2388-
2388+
23892389 // define if the band has scale and offset to apply
23902390 double myScale = bandScale ( theBandNo );
23912391 double myOffset = bandOffset ( theBandNo );
23922392 if ( myScale != 1.0 || myOffset != 0.0 )
23932393 {
2394- if ( myScale < 0.0 )
2395- {
2396- // update Min and Max value
2397- myRasterBandStats.minimumValue = pdfMax * myScale + myOffset;
2398- myRasterBandStats.maximumValue = pdfMin * myScale + myOffset;
2399- // update the range
2400- myRasterBandStats.range = ( pdfMin - pdfMax) * myScale;
2401- // update standard deviation
2402- myRasterBandStats.stdDev = -1.0 * pdfStdDev * myScale;
2403- }
2404- else
2405- {
2406- // update Min and Max value
2407- myRasterBandStats.minimumValue = pdfMin * myScale + myOffset;
2408- myRasterBandStats.maximumValue = pdfMax * myScale + myOffset;
2409- // update the range
2410- myRasterBandStats.range = ( pdfMax - pdfMin) * myScale;
2411- // update standard deviation
2412- myRasterBandStats.stdDev = pdfStdDev * myScale;
2413- }
2414- // update the mean
2415- myRasterBandStats.mean = pdfMean * myScale + myOffset;
2394+ if ( myScale < 0.0 )
2395+ {
2396+ // update Min and Max value
2397+ myRasterBandStats.minimumValue = pdfMax * myScale + myOffset;
2398+ myRasterBandStats.maximumValue = pdfMin * myScale + myOffset;
2399+ // update the range
2400+ myRasterBandStats.range = ( pdfMin - pdfMax ) * myScale;
2401+ // update standard deviation
2402+ myRasterBandStats.stdDev = -1.0 * pdfStdDev * myScale;
2403+ }
2404+ else
2405+ {
2406+ // update Min and Max value
2407+ myRasterBandStats.minimumValue = pdfMin * myScale + myOffset;
2408+ myRasterBandStats.maximumValue = pdfMax * myScale + myOffset;
2409+ // update the range
2410+ myRasterBandStats.range = ( pdfMax - pdfMin ) * myScale;
2411+ // update standard deviation
2412+ myRasterBandStats.stdDev = pdfStdDev * myScale;
2413+ }
2414+ // update the mean
2415+ myRasterBandStats.mean = pdfMean * myScale + myOffset;
24162416 }
24172417
24182418#ifdef QGISDEBUG
0 commit comments