Disable bounding box in calc_pixmap

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

Hide file tree

Showing 4 changed files with 21 additions and 3 deletions.
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,17 @@ 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:
+    orig_bbox = wcs_to.bounding_box
+    wcs_to.bounding_box = None
+    try:
+        x, y = wcs_to.world_to_pixel_values(
+            *wcs_from.pixel_to_world_values(x, y)
+        )
+    finally:
+        wcs_to.bounding_box = orig_bbox
+
     pixmap = np.dstack([x, y])
     return pixmap
 

diff --git a/pyproject.toml b/pyproject.toml
@@ -7,7 +7,6 @@ authors = [
 ]
 dependencies = [
     "numpy",
-    "astropy",
 ]
 dynamic = [
     "version",
@@ -29,6 +28,7 @@ Documentation = "http://spacetelescope.github.io/drizzle/"
 
 [project.optional-dependencies]
 test = [
+    "astropy",
     "pytest",
     "pytest-cov",
     "pytest-doctestplus",