Use astropy pixel_to_pixel in calc_pixmap

Author: mcara

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)
 ==================
 

drizzle/utils.py

@@ -2,6 +2,8 @@
 
 import numpy as np
 
+from astropy.wcs.utils import pixel_to_pixel
+
 __all__ = ["calc_pixmap", "decode_context", "estimate_pixel_scale_ratio"]
 
 _DEG2RAD = math.pi / 180.0
@@ -74,7 +76,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
 

pyproject.toml

@@ -7,7 +7,7 @@ authors = [
 ]
 dependencies = [
     "numpy",
-    "astropy",
+    "astropy>=5.0",
 ]
 dynamic = [
     "version",