Use astropy pixel_to_pixel in calc_pixmap

spacetelescope · Oct 25, 2024 · 99d15a7 · 99d15a7
1 parent 6b0a118
commit 99d15a7
Showing 3 changed files with 10 additions and 2 deletions.
diff --git a/CHANGES.rst b/CHANGES.rst
@@ -4,6 +4,13 @@
 Release Notes
 =============
 
+
+2.0.1 (unreleased)
+==================
+
+- Use ``astropy``'s ``pixel_to_pixel()`` in ``utils.calc_pixmap()``. [#161]
+
+
 2.0.0 (2024-10-23)
 ==================
 

diff --git a/drizzle/utils.py b/drizzle/utils.py
@@ -1,5 +1,6 @@
 import math
 
+from astropy.wcs.utils import pixel_to_pixel
 import numpy as np
 
 __all__ = ["calc_pixmap", "decode_context", "estimate_pixel_scale_ratio"]
@@ -74,7 +75,7 @@ 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))
+    x, y = pixel_to_pixel(x, y)
     pixmap = np.dstack([x, y])
     return pixmap
 

diff --git a/pyproject.toml b/pyproject.toml
@@ -7,7 +7,7 @@ authors = [
 ]
 dependencies = [
     "numpy",
-    "astropy",
+    "astropy>=5.0",
 ]
 dynamic = [
     "version",