Re: [gdal-dev] Prevent extrapolating with ReadAsArray when resampling

2018-01-05 Thread Even Rouault
Rutger, ReadAsArray() / IRasterIO() with non-nearest resampling use the interpolation code of overview building (except for Lanczos where it uses the warping code), which was of course targetted at downsampling rather than upsampling. Nodata is taken into account when mixing together source

[gdal-dev] Prevent extrapolating with ReadAsArray when resampling

2018-01-04 Thread Rutger
Dear list, I regularly use (and enjoy!) the resampling option when using the ReadAsArray method of a Dataset. However, i noticed that in the case of nodata values in the raster, the method extrapolates "into" those nodata values. This is not always desirable in my case, especially when using any r