sharpetronics/QGIS
1
0
Fork 0
mirror of https://github.com/qgis/QGIS.git synced 2025-10-04 00:04:03 -04:00
Code Issues Packages Projects Releases Wiki Activity

[gdal] Histogram: adjust bin count according to scale

Browse Source 
Fix #59461
This commit is contained in:
Alessandro Pasotti 2025-09-22 12:10:06 +02:00
parent 0b1ba631b3
commit 37a09e66ff
2 changed files with 33 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
src/core/providers/gdal/qgsgdalprovider.cpp
View File

@ -2040,6 +2040,11 @@ QgsRasterHistogram QgsGdalProvider::histogram( int bandNo,
}
}
const double myScale { bandScale( bandNo ) };
// Adjust bin count according to scale
// Fix issue GH #59461
myHistogram.binCount = static_cast<int>( myHistogram.binCount * myScale );
if ( ( sourceHasNoDataValue( bandNo ) && !useSourceNoDataValue( bandNo ) ) ||
!userNoDataValues( bandNo ).isEmpty() )
{
@ -2092,9 +2097,8 @@ QgsRasterHistogram QgsGdalProvider::histogram( int bandNo,
double myMinVal = myHistogram.minimum;
double myMaxVal = myHistogram.maximum;
// unapply scale anf offset for min and max
double myScale = bandScale( bandNo );
double myOffset = bandOffset( bandNo );
// unapply scale and offset for min and max
const double myOffset = bandOffset( bandNo );
if ( myScale != 1.0 || myOffset != 0. )
{
myMinVal = ( myHistogram.minimum - myOffset ) / myScale;

26
tests/src/python/test_provider_gdal.py
View File

@ -550,6 +550,32 @@ class PyQgsGdalProvider(QgisTestCase, RasterProviderTestCase):
# This was triggering another crash:
self.assertEqual(rl.dataProvider().colorInterpretationName(1), "Undefined")
def testHistogramBinCountWithScale(self):
"""Test issue GH #59461"""
# Create raster with 2x2 pixels, int type, values 0-10000, scale 0.001
tmp_dir = QTemporaryDir()
tmpfile = os.path.join(tmp_dir.path(), "testBinCountWithScale.tif")
ds = gdal.GetDriverByName("GTiff").Create(tmpfile, 2, 2, 1, gdal.GDT_Int32)
ds.SetGeoTransform([0, 1, 0, 0, 0, -1])
ds.SetProjection("EPSG:4326")
ds.WriteRaster(0, 0, 2, 2, struct.pack("i" * 4, 0, 5000, 10000, -10000))
# Set band scale
band = ds.GetRasterBand(1)
band.SetScale(0.001)
ds = None
rl = QgsRasterLayer(tmpfile, "test")
self.assertTrue(rl.isValid())
provider = rl.dataProvider()
# Data type is changed by QGIS becase of the scale
self.assertEqual(provider.dataType(1), Qgis.DataType.Float32)
self.assertEqual(provider.bandScale(1), 0.001)
# Create histogram
hist = provider.histogram(1, 10000)
self.assertIsNotNone(hist)
self.assertEqual(hist.binCount, 10)
if __name__ == "__main__":
unittest.main()