With count data, like population, it only makes sense to aggregate using a sum resampling function (like you did with block_reduce). Using rasterio you can perform a decimated read (setting out_shape)

Note that setting out_shape to (1200, 1200) will likely read a (1200, 1199) shaped chunk of data and resample it to (1200, 1200); which is not exactly the same as reading a 1200*1200 window. Kind rega

Many thanks Sean and Even, So with rasterio compiled against gdal 3.3 I'll be able to force ignore overviews? If that deserves a new rasterio unit test I'm happy to send a pull request in the coming d

Hi Denis, Thanks for pointing to that issue. I may have misunderstood it but I don't think it helps answering my question. Cheers, Loïc ________________________________________ Sent: 10 November 2020

Hi everyone, I'm trying to perform a decimated read of a dataset that already contains overviews. Regardless of which value I pass to the resampling= argument of the read method, it seems that the alr

Hi Sean, Thank you for your response. Is it correct to interpret from your explanation that: 1- currently the reproject function will never use overviews of the source file, even when it involves down

Hi everyone, I read that the gdalwarp command line defaults to using the overview level nearest to target resolution. What about `rasterio.warp.reproject()`? Does it use overview at all? Thank you and