Fix reproject_shapes() to avoid flipped coordinats (#462)

* Fix reproject_shapes() to avoid flipped coordinats

Fix reproject_shapes() by enforcing always_xy=True in Transformer

Previously, reprojections involving EPSG:4326 could produce flipped
coordinates because pyproj defaults to EPSG axis order (lat, lon).
This caused inconsistent results when transforming from projected
CRS (e.g. UTM) to geographic CRS.

The fix explicitly sets always_xy=True in Transformer.from_crs,
ensuring that all transformations consistently interpret coordinates
as (x=longitude, y=latitude).

* Add release note

Author: juliusmeier

RELEASE_NOTES.rst

@@ -29,6 +29,7 @@ Upcoming Release
 * Raise a `FileNotFoundError` if the `temp_dir` explicitly specified for cutout preparation does not exist instead of failing with an obfuscated error message (https://github.com/PyPSA/atlite/pull/445).
 * Addressed `rasterio` DeprecationWarning on `crs.is_valid`.
 * Fix calls to `cdsapi` for ERA5 to be compliant with current API syntax (https://github.com/PyPSA/atlite/pull/414 and https://github.com/PyPSA/atlite/pull/454).
+* Fix `reproject_shapes()` to preserve coordinate order by enforcing (x=lon, y=lat) with `always_xy=True` when transforming from one CRS to another (https://github.com/PyPSA/atlite/pull/462). 
 
 `v0.4.1 <https://github.com/PyPSA/atlite/releases/tag/v0.4.1>`__ (12th May 2025)
 =======================================================================================

atlite/gis.py

@@ -88,7 +88,7 @@ def reproject_shapes(shapes, crs1, crs2):
     """
     Project a collection of shapes from one crs to another.
     """
-    transformer = Transformer.from_crs(crs1, crs2)
+    transformer = Transformer.from_crs(crs1, crs2, always_xy=True)
 
     def _reproject_shape(shape):
         return transform(transformer.transform, shape)