Hello, I am trying to reproject a population distribution raster file of China into 3857 space. I tried but the images are not aligning. Looking at the CRS it has a lot of different standards, I am

Hi, We fixed a window shift bug in 1.1.5. See https://github.com/mapbox/rasterio/pull/1932 https://github.com/mapbox/rasterio/pull/1937 https://github.com/mapbox/rasterio/pull/1938 I wonder if this is

I am using rasterio 1.1.5.

Thanks for pointing that out. But it's really odd. I've been trying to update rasterio using conda, but the version remains at v1.1.0. I just tried doing a new install using conda install -c

The latest version is 1.1.5 (Jun 3, 2020). 1.1.0 has been released in 2019.

Yep - I'm using 1.1.0

Do you use the latest version of rasterio?

I'm using merge with method = 'max' to merge two rasters with identical metadata. When I open the output file and the input files in QGIS, the top part of the output raster is perfectly merged, but

Hello, I've been having a problem using windows, probably a misunderstanding on my part. My goal is to dump a specific data from a band into a csv file, first I tried to use a window but the results

I originally had a geotiff (size= 4490 x 7341) which I tiled into tiles of size 256 x 256 to perform an operation. The result is a list of 522 numpy array. From the original geotiff I have the geo

Thanks for letting me know! Thumbs up for your work!

Hi Thomas, I think Denis found the root of the problem and we have a fix at https://github.com/mapbox/rasterio/pull/1938 that will be a part of the 1.1.5 release. Thank you for the test cases and for

The issue has been fixed via https://github.com/mapbox/rasterio/commit/15d8dd759256f12d412433f5b9824c8d16b33f8d.

FYI...I was able to achieve the desired behavior by "lying" about the nodata value. The actual nodata was nan, but when I specified something else (e.g., nodata=-9999) in the profile, the resampled

Here is an example using rioxarray: https://corteva.github.io/rioxarray/stable/examples/interpolate_na.html It uses scipy griddata.

I see what you mean. Rasterio doesn't have a built in option for this, unless there is a GDAL resampling trick that I am not aware of. -- Sean Gillies

Hi Sean. Thanks for the message. I understand your example, but I may have been confusing in my email. I think I'm asking something different. Is there a way to force the resampling process to not

Hi Nick, Ignoring nodata values is the default behavior for rasterio, as it is for the GDAL library that rasterio uses. Here's an example using one of rasterio's test files. $ rio insp

rasterio.transform.array_bounds should do the trick - https://rasterio.readthedocs.io/en/latest/api/rasterio.transform.html#rasterio.transform.array_bounds

I am not 100 sure, but I now simply rounded all windows (origin + offsets) and this seems to have fixed it.