Re: Rasterio 1.2.4

Even Rouault

Hi Sean

regarding the change in web mercator, I guess this is due to a subtle change where we don't check anymore if latitude is close enough of PI/2 to decide the point doesn't project. Now we only test lat < PI/2. And as 90° converted in radian is slightly < PI/2, it now projects to a finite (although very big) northing value. So this changes the aspect of how a raster with world wide bounds in EPSG:4326 reprojects in EPSG:3857. Neither the previous (which truncated a few pixels) nor the new behavior (which results in an image with a big height/width ratio) are ideal, but I don't think there isn't an ideal solution given the characteristics of the mercator projection. The commit that modified this is

Regarding discovery of proj.db, I can only think to which might have effect on this. It should mostly be a no-op for default builds, although I see that line might have subtle effects if the proj_lib directory hardcoded at build time doesn't exist at runtime. Previously it will have been used even if not existing, whereas now it is not considered (and thus a proj.db in the current working directory might be used). The new behavior makes more sense to me.


Le 01/06/2021 à 18:37, Sean Gillies via a écrit :
Hi all,

Here is the list of bug fixes in rasterio 1.2.4: 

The wheels on PyPI for this version now include GDAL 3.3.0 and GEOS 3.9.1. The PROJ version remains the same, 7.2.1 until we get more experience with changes to web mercator projection (I believe I see expanded limits) and discovery of the proj.db file.

Thanks for the bug reports and discussions, everyone! 

Sean Gillies
My software is free, but my time generally not.

Join to automatically receive all group messages.