Mismatch between actual minimum pixel value and statistics minimum

[phargogh]

Greetings! I have a TiTiler instance (which of course is using rio-tiler under the hood), and I have a COG (https://storage.googleapis.com/natcap-data-cache/global/landscan-pop/landscan_2023.tif) that I am attempting to render tiles for. The problem is that rio-tiler appears to be incorrectly reporting the minimum pixel value.

When I inspect the raster with gdalinfo, I see the correct minimum value (0), and the raster's nodata value is correctly stated as well (2147483647.0). gdalinfo https://storage.googleapis.com/natcap-data-cache/global/landscan-pop/landscan_2023.tif will demonstrate this.

And when I download the file and read in the raster with GDAL, I can see that the minimum pixel value matches the nodata value:

>>> from osgeo import gdal
>>> gdal.Open('landscan_2023.tif').ReadAsArray().min()
-2147483647

However, when I use rio-tiler to get statistics, the minimum value is incorrectly stated:

>>> from rio_tiler.io import Reader
>>> Reader('https://storage.googleapis.com/natcap-data-cache/global/landscan-pop/landscan_2023.tif').statistics()['b1'].min
-2147483648.0

Because of this mismatch between the nodata value and the rio-tiler-identified pixel minimum, the rest of the percentile breaks are being incorrectly identified.

Any suggestions? Thanks!