HyBIG 2.8.0: Improved reprojection support and bug fixes (#67)

* GITC-9048: Fix latitude max < latitude min error when transforming bounds, fix polar projection, fix adapter

* GITC-9047/9048/9148: Updated unit tests, removed bad read window code, updated changelog

* Update dependency versions

Author: jackiryan

CHANGELOG.md

@@ -4,6 +4,14 @@ 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/).
 
+## [2.8.0] - 2026-07-17
+
+* [GITC-9047](https://bugs.earthdata.nasa.gov/browse/GITC-9047): Updated adapter code to now log a message when a request produces no data, but not throw a service error. Previously, jobs with no output would produce this error: `WorkItem failed: nasa/harmony-browse-image-generator:2.7.0: list index out of range`
+* [GITC-9048](https://bugs.earthdata.nasa.gov/browse/GITC-9048): Fixed a bug where some granules fail to process due to an inverted latitude coordinate system.
+* [GITC-9148](https://bugs.earthdata.nasa.gov/browse/GITC-9148): Removed an "optimization" when reading from the source data that caused an unintended bug when reprojecting some types of data from EPSG:4326 to EPSG:3413. The consequence of removing this code is that requests for downsampled browse images will use more memory, but this is not a common mode for HyBIG.
+* Update harmony-service-lib to 3.0 and pillow to 12.3.0
+
+
 ## [2.7.0] - 2026-05-11
 
 ### Changed

docker/service_version.txt

@@ -1 +1 @@
-2.7.0
+2.8.0

harmony_service/adapter.py

@@ -121,7 +121,13 @@ def process_item(self, item: Item, source: HarmonySource) -> Item:
                 # image_file_list is a list of tuples (image, world, auxiliary)
                 # we need to stage them each individually, and then add their final
                 # locations to a list before creating the stac item.
-                item_assets = []
+                item_assets: list[tuple[str, str, str]] = []
+
+                if not image_file_list:
+                    # If output is empty, log an error message but don't throw an
+                    # Exception, just return an empty stac item
+                    self.logger.error('No output assets produced from HyBIG call.')
+                    return self.create_output_stac_item(item, item_assets)
 
                 for (
                     browse_image_name,

hybig/browse.py

@@ -7,7 +7,7 @@
 
 import numpy as np
 import rasterio
-from affine import Affine, dumpsw
+from affine import dumpsw
 from harmony_service_lib.message import Message as HarmonyMessage
 from harmony_service_lib.message import Source as HarmonySource
 from matplotlib.colors import BoundaryNorm, Normalize
@@ -228,22 +228,10 @@ def process_tile(
     if src_window is None:
         return False
 
-    # Compute downsampled read dimensions: read at output resolution rather than
-    # full source resolution, avoiding loading the entire source into memory when
-    # the output is much smaller than the source (e.g. height=1280, width=2560).
-    src_pixel_x = abs(src_ds.transform.a)
-    src_pixel_y = abs(src_ds.transform.e)
-    dst_pixel_x = abs(grid_params['transform'].a)
-    dst_pixel_y = abs(grid_params['transform'].e)
+    read_width = int(src_window.width)
+    read_height = int(src_window.height)
 
-    win_height = int(src_window.height)
-    win_width = int(src_window.width)
-
-    # Cap at full window size so we never upsample during the read.
-    read_width = min(win_width, max(1, round(win_width * src_pixel_x / dst_pixel_x)))
-    read_height = min(win_height, max(1, round(win_height * src_pixel_y / dst_pixel_y)))
-
-    # Explicitly load a subset of ds at the target resolution
+    # Explicitly load a subset of the dataset needed for the browse image tile
     tile_source = read_window_with_mask_and_scale(
         src_ds, src_window, out_shape=(band_count, read_height, read_width)
     )
@@ -255,11 +243,7 @@ def process_tile(
         return False
 
     src_crs = src_ds.crs
-    # Adjust the window transform to reflect the downsampled pixel size.
-    window_transform = src_ds.window_transform(src_window)
-    src_transform = window_transform * Affine.scale(
-        win_width / read_width, win_height / read_height
-    )
+    src_transform = src_ds.window_transform(src_window)
 
     dst_nodata = TRANSPARENT
 

hybig/sizes.py

@@ -151,7 +151,18 @@ def choose_scale_extent(
             }
         )
     else:
-        left, bottom, right, top = transform_bounds(src_ds.crs, dst_crs, *src_ds.bounds)
+        # Normalize the source bounds before transforming. Some inputs have a
+        # flipped (south-up) geotransform with a positive Y resolution, which
+        # makes src_ds.bounds report top < bottom. transform_bounds requires a
+        # properly ordered box (min/max), otherwise PROJ raises
+        # "latitude max < latitude min."
+        src_left = min(src_ds.bounds.left, src_ds.bounds.right)
+        src_right = max(src_ds.bounds.left, src_ds.bounds.right)
+        src_bottom = min(src_ds.bounds.bottom, src_ds.bounds.top)
+        src_top = max(src_ds.bounds.bottom, src_ds.bounds.top)
+        left, bottom, right, top = transform_bounds(
+            src_ds.crs, dst_crs, src_left, src_bottom, src_right, src_top
+        )
 
         # Correct for antimeridian crossing.
         if left > right:

pip_requirements.txt

@@ -1,7 +1,7 @@
-harmony-service-lib~=2.11.0
+harmony-service-lib~=3.0.0
 matplotlib==3.9.0
 numpy==1.26.4
-pillow==12.2.0
+pillow~=12.3.0
 pyproj==3.6.1
 pystac~=1.0.1
 rasterio==1.3.10

tests/unit/test_browse.py

@@ -13,6 +13,7 @@
 from numpy.testing import assert_array_equal, assert_equal
 from osgeo_utils.auxiliary.color_palette import ColorPalette
 from rasterio import Affine
+from rasterio.coords import BoundingBox
 from rasterio.crs import CRS
 from rasterio.io import DatasetReader, DatasetWriter
 from rasterio.warp import Resampling
@@ -161,7 +162,7 @@ def test_create_browse_imagery_with_mocks(self, rasterio_open_mock, reproject_mo
         ds.crs = CRS.from_string('EPSG:4326')
         ds.count = 1
         ds.colormap = Mock(side_effect=ValueError)
-        ds.bounds = (-180.0, -90.0, 180.0, 90.0)
+        ds.bounds = BoundingBox(-180.0, -90.0, 180.0, 90.0)
         ds.window_transform = Mock(return_value=file_transform)
         ds.nodatavals = (255,)
         ds.scales = (1,)
@@ -313,7 +314,7 @@ def test_create_browse_imagery_excludes_all_nan_tiles(self, rasterio_open_mock):
         ds.crs = CRS.from_string('EPSG:4326')
         ds.count = 1
         ds.colormap = Mock(side_effect=ValueError)
-        ds.bounds = (-180.0, -90.0, 180.0, 90.0)
+        ds.bounds = BoundingBox(-180.0, -90.0, 180.0, 90.0)
         ds.window_transform = Mock(return_value=file_transform)
         ds.nodatavals = (None,)
         ds.scales = (1,)
@@ -368,74 +369,6 @@ def test_read_window_with_mask_and_scale_without_out_shape(self):
         ds.read.assert_called_once_with([1], window=window)
         self.assertEqual(result.shape, (1, 4, 4))
 
-    @patch('hybig.browse.reproject')
-    @patch('rasterio.open')
-    def test_process_tile_downsamples_read_for_coarser_output(
-        self, rasterio_open_mock, reproject_mock
-    ):
-        """Test process_tile reads at output resolution when output is coarser
-        than source.
-
-        This is the memory fix: a 36000x18000 source with a 1280x2560 output
-        should not load the full-resolution source into memory.
-        """
-        # Source: 1° pixels, 10x10
-        src_affine = Affine(1.0, 0.0, -5.0, 0.0, -1.0, 5.0)
-        ds = Mock(spec=DatasetReader)
-        ds.crs = CRS.from_string('EPSG:4326')
-        ds.transform = src_affine
-        ds.shape = (10, 10)
-        ds.count = 1
-        ds.nodatavals = (None,)
-        ds.scales = (1,)
-        ds.offsets = (0,)
-        # Return data at the downsampled (2x2) size
-        ds.read.return_value = np.array([[[0, 100], [200, 255]]], dtype='float64')
-        ds.window_transform.return_value = src_affine
-
-        dest_write_mock = Mock(spec=DatasetWriter)
-        rasterio_open_mock.return_value.__enter__.return_value = dest_write_mock
-
-        # Output: 5° pixels, 2x2 — coarser than source
-        out_affine = Affine(5.0, 0.0, -5.0, 0.0, -5.0, 5.0)
-        grid_params = GridParams(
-            {
-                'height': 2,
-                'width': 2,
-                'crs': CRS.from_string('EPSG:4326'),
-                'transform': out_affine,
-            }
-        )
-
-        result = process_tile(
-            ds,
-            grid_params,
-            None,
-            'PNG',
-            self.tmp_dir / 'output.png',
-            self.tmp_dir / 'output.pgw',
-            self.logger,
-        )
-
-        self.assertTrue(result)
-
-        # Source window covers the full 10x10 source (with buffer, clamped).
-        # read_width  = min(10, round(10 * 1.0 / 5.0)) = 2
-        # read_height = min(10, round(10 * 1.0 / 5.0)) = 2
-        ds.read.assert_called_once()
-        read_kwargs = ds.read.call_args.kwargs
-        self.assertIn('out_shape', read_kwargs)
-        _bands, read_height, read_width = read_kwargs['out_shape']
-        self.assertEqual(read_width, 2)
-        self.assertEqual(read_height, 2)
-
-        # src_transform passed to reproject must have 5° pixel size (scaled from 1°)
-        self.assertEqual(reproject_mock.call_count, 1)
-        actual_src_transform = reproject_mock.call_args.kwargs['src_transform']
-        self.assertAlmostEqual(actual_src_transform.a, 5.0)
-        self.assertAlmostEqual(abs(actual_src_transform.e), 5.0)
-        self.assertAlmostEqual(actual_src_transform.c, -5.0)  # origin unchanged
-
     @patch('hybig.browse.reproject')
     @patch('rasterio.open')
     def test_process_tile_full_res_when_resolutions_match(