Disable bounding box in calc_pixmap

spacetelescope · Dec 16, 2024 · c5fff59 · c5fff59
1 parent 8978081
commit c5fff59
Show file tree

Hide file tree

Showing 5 changed files with 25 additions and 5 deletions.
diff --git a/.github/workflows/ci.yml b/.github/workflows/ci.yml
@@ -22,15 +22,13 @@ jobs:
       envs: |
         - linux: check-style
         - linux: check-security
-
         - linux: py310-xdist
         - linux: py311-xdist
         - linux: py313-xdist
         - macos: py312-xdist
         - windows: py312-xdist
         - linux: py312-xdist-cov
           coverage: codecov
-
         - linux: py312-xdist-devdeps
     secrets:
       CODECOV_TOKEN: ${{ secrets.CODECOV_TOKEN }}
diff --git a/CHANGES.rst b/CHANGES.rst
@@ -4,6 +4,14 @@
 Release Notes
 =============
 
+2.0.1 (2024-12-16)
+==================
+
+- Update ``utils.calc_pixmap`` code to be ready for upcoming changes in GWCS
+  due to which inverse WCS transformations will respect bounding box by
+  disabling the bounding box when computing pixel map. [#164]
+
+
 2.0.0 (2024-10-23)
 ==================
 

diff --git a/drizzle/resample.py b/drizzle/resample.py
@@ -6,7 +6,7 @@
 
 from drizzle import cdrizzle
 
-__all__ = ["blot_image", "Drizzle"]
+__all__ = ["Drizzle", "blot_image"]
 
 SUPPORTED_DRIZZLE_KERNELS = [
     "square",

diff --git a/drizzle/utils.py b/drizzle/utils.py
@@ -74,7 +74,21 @@ def calc_pixmap(wcs_from, wcs_to, shape=None):
         )
 
     y, x = np.indices(shape, dtype=np.float64)
-    x, y = wcs_to.world_to_pixel_values(*wcs_from.pixel_to_world_values(x, y))
+
+    # temporarily disable the bounding box:
+    if hasattr(wcs_to, "bounding_box"):
+        orig_bbox = getattr(wcs_to, "bounding_box", None)
+        wcs_to.bounding_box = None
+    else:
+        orig_bbox = None
+    try:
+        x, y = wcs_to.world_to_pixel_values(
+            *wcs_from.pixel_to_world_values(x, y)
+        )
+    finally:
+        if orig_bbox is not None:
+            wcs_to.bounding_box = orig_bbox
+
     pixmap = np.dstack([x, y])
     return pixmap
 

diff --git a/pyproject.toml b/pyproject.toml
@@ -6,8 +6,8 @@ authors = [
     { name = "STScI", email = "[email protected]" },
 ]
 dependencies = [
-    "numpy",
     "astropy",
+    "numpy",
 ]
 dynamic = [
     "version",