熵WCS转换矩阵

如何解决熵WCS转换矩阵

我正在尝试创建一个自定义WCS，以将图像的像素坐标转换为世界坐标。

给出一个带有星星的图像，我已经确定了2个星星，因此我可以将图像中的两个点的像素（x，y）匹配到（RA，DEC）。

我现在想要的是使用适当的变换矩阵创建一个自定义WCS，因此当我提供任何像素坐标时，它将返回相应的RA和DEC。

我知道astrometry.net做到了，并用适当的转换矩阵写了一个fits标头。

我的问题是，如何获得此转换矩阵并创建自定义WCS对象？

谢谢。

编辑： 这是我正在尝试的代码：

from astropy.coordinates import SkyCoord
import astropy.units as u
from astropy.wcs import WCS
from astropy.wcs.utils import fit_wcs_from_points
import numpy as np
# I have the following stars identified in my image (andromeda):
#  (X,Y)       --> (RA in degrees,DEC in degrees) --> HIP_ID
#  (640,555)   --> (17.43421495,35.61993419)      --> 5447
#  (1076,32)  --> (2.09777329,29.08952671)        --> 607
#  (161,903)  --> (30.9751282,42.32944223)        --> 9640
#  (932,327)  --> (9.83272908,30.86056254)        --> 3092
stars = SkyCoord(ra=[17.43421495,2.09777329,30.9751282,9.83272908],dec=[35.61993419,29.08952671,42.32944223,30.86056254],unit=u.deg)
stars
pixels_x = np.array([640,1076,161,932])
pixels_y = np.array([555,32,903,327])
wcs = fit_wcs_from_points((pixels_x,pixels_y),stars); wcs
wcs.wcs_pix2world(np.array([[640]]),np.array([555]),0)

为什么我什至没有正确设置参考点？

解决方法

代码正确，问题在于Astropy库的版本，请参见错误：https://github.com/astropy/astropy/pull/10155。

如果有人遇到此问题，请确保您未使用Astropy 4.0.1。