Fix issues in outlier det. due to GWCS.inverse enforcing bbox

spacetelescope · Dec 16, 2024 · bae0418 · bae0418
1 parent 4953380
commit bae0418
Showing 1 changed file with 12 additions and 1 deletion.
diff --git a/src/stcal/outlier_detection/utils.py b/src/stcal/outlier_detection/utils.py
@@ -284,7 +284,18 @@ def calc_gwcs_pixmap(in_wcs, out_wcs, in_shape):
     log.debug("Bounding box from data shape: {}".format(bb))
 
     grid = gwcs.wcstools.grid_from_bounding_box(bb)
-    return np.dstack(reproject(in_wcs, out_wcs)(grid[0], grid[1]))
+
+    # temporarily disable the bounding box:
+    orig_bbox = out_wcs.bounding_box
+    out_wcs.bounding_box = None
+    try:
+        pixmap = np.dstack(reproject(in_wcs, out_wcs)(grid[0], grid[1]))
+    finally:
+        # restore bounding box:
+        if orig_bbox is not None:
+            out_wcs.bounding_box = orig_bbox
+
+    return pixmap
 
 
 def reproject(wcs1, wcs2):