Merge pull request #651 from mraspaud/fix_freeze_epsg

Fix dynamic area freezing not preserving epsg code

Author: djhoese

pyresample/geometry.py

@@ -22,6 +22,7 @@
 import math
 import warnings
 from collections import OrderedDict
+from contextlib import suppress
 from functools import partial, wraps
 from logging import getLogger
 from pathlib import Path
@@ -66,6 +67,7 @@
 
 from pyproj import CRS
 from pyproj.enums import TransformDirection
+from pyproj.exceptions import CRSError
 
 logger = getLogger(__name__)
 
@@ -1284,6 +1286,8 @@ def freeze(self, lonslats=None, resolution=None, shape=None, proj_info=None,
         area_extent = self.area_extent
         if not area_extent or not width or not height:
             projection, corners = self._compute_bound_centers(proj_dict, lonslats, antimeridian_mode=antimeridian_mode)
+            with suppress(CRSError):
+                projection = CRS(CRS(projection).to_epsg())
             area_extent, width, height = self.compute_domain(corners, resolution, shape, projection)
         return AreaDefinition(self.area_id, self.description, '',
                               projection, width, height,

pyresample/test/test_geometry/test_area.py

@@ -1763,3 +1763,12 @@ def test_dynamic_area_can_use_bounding_box_attribute():
     swath_def_to_freeze_on = SwathDefinition(None, None, attrs=dict(bounding_box=[[0, 20, 20, 0], [55, 55, 45, 45]]))
     res_area = area_def.freeze(swath_def_to_freeze_on)
     assert res_area.area_extent == (3533500, 2484500, 5108500, 3588500)
+    assert res_area.crs == CRS("epsg:3035")
+
+
+def test_dynamic_area_preserves_epsg_code():
+    """Test that area freezing preserves epsg code."""
+    area_def = DynamicAreaDefinition("test_area", "", "epsg:3035", resolution=500)
+    swath_def_to_freeze_on = SwathDefinition(None, None, attrs=dict(bounding_box=[[0, 20, 20, 0], [55, 55, 45, 45]]))
+    res_area = area_def.freeze(swath_def_to_freeze_on)
+    assert res_area.crs == CRS("epsg:3035")