GITC-7067: shift negative bounds to handle antimeridian crossing (#48)

Author: flamingbear

CHANGELOG.md

@@ -4,6 +4,13 @@ HyBIG follows semantic versioning. All notable changes to this project will be
 documented in this file. The format is based on [Keep a
 Changelog](http://keepachangelog.com/en/1.0.0/).
 
+
+## [v2.3.0] - 2025-02-26
+
+### Changed
+
+* Fix images that cross the antimeridian. Target extents are corrected when the bounding box crosses the dateline. [[#48](https://github.com/nasa/harmony-browse-image-generator/pull/48)]
+
 ## [v2.2.0] - 2024-12-19
 
 ### Changed
@@ -105,6 +112,7 @@ For more information on internal releases prior to NASA open-source approval,
 see legacy-CHANGELOG.md.
 
 [unreleased]: https://github.com/nasa/harmony-browse-image-generator/
+[v2.3.0]: https://github.com/nasa/harmony-browse-image-generator/releases/tag/2.3.0
 [v2.2.0]: https://github.com/nasa/harmony-browse-image-generator/releases/tag/2.2.0
 [v2.1.0]: https://github.com/nasa/harmony-browse-image-generator/releases/tag/2.1.0
 [v2.0.2]: https://github.com/nasa/harmony-browse-image-generator/releases/tag/2.0.2

docker/service_version.txt

@@ -1 +1 @@
-2.2.0
+2.3.0

hybig/sizes.py

@@ -153,6 +153,10 @@ def choose_scale_extent(
         )
     else:
         left, bottom, right, top = data_array.rio.transform_bounds(target_crs)
+
+        # Correct for antimeridian crossing.
+        if left > right:
+            right = right + 360
         scale_extent = ScaleExtent(
             {'xmin': left, 'ymin': bottom, 'xmax': right, 'ymax': top}
         )

tests/fixtures/split-dateline-sample.tif

@@ -458,6 +458,29 @@ def test_scale_extent_from_input_image_with_crs_transformation(self):
                 actual_scale_extent, rel=1e-12
             )
 
+    def test_scale_extent_from_input_image_that_crosses_antimeridian(self):
+        """Verify that bounding box is correctly shifted at dateline.
+
+        Notice that the xmax value is > 180.
+        """
+        target_crs = CRS.from_string(PREFERRED_CRS['global'])
+        with open_rasterio(
+            self.fixtures / 'split-dateline-sample.tif', mode='r', mask_and_scale=True
+        ) as in_array:
+            expected_scale_extent = ScaleExtent(
+                {
+                    'xmin': 179.25694918947116,
+                    'ymin': -16.351723987726583,
+                    'xmax': 180.29656915724874,
+                    'ymax': -15.345453705696787,
+                }
+            )
+
+            actual_scale_extent = choose_scale_extent({}, target_crs, in_array)
+            assert expected_scale_extent == pytest.approx(
+                actual_scale_extent, rel=1e-12
+            )
+
 
 class TestChooseTargetDimensions(TestCase):
     def test_message_has_dimensions(self):