diff --git a/src/core/providers/gdal/qgsgdalprovider.cpp b/src/core/providers/gdal/qgsgdalprovider.cpp
index d4796522c58..3c142e12f07 100644
--- a/src/core/providers/gdal/qgsgdalprovider.cpp
+++ b/src/core/providers/gdal/qgsgdalprovider.cpp
@@ -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;
diff --git a/tests/src/python/test_provider_gdal.py b/tests/src/python/test_provider_gdal.py
index a9f75844205..124f9938b79 100644
--- a/tests/src/python/test_provider_gdal.py
+++ b/tests/src/python/test_provider_gdal.py
@@ -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()