SIP兼容的Python板解决包

Question

我需要执行精确的像素为世界协调FITS原先采用Maxim DL创建的文件转换。 Maxim使用小点为板求解生成TRi_j失真系数。这些是不相容与astropy.wcs坐标我建议使用这些假设SIP畸变系数变换的功能。

因此，我正在寻找的选项重新platesolve的FITS文件生成SIP系数。

到目前为止，所有我发现是astrometry.net但是这是一个在线服务。我真的很为脱机platesolving（优选地针对GSC的本地副本），我可以为我的应用程序的工作流程的一部分同步执行。

是否有任何Astropy下属（或其他）Python包执行对GSC SIP兼容platesolving？

可替代地，是否有任何等同物wcs.all_pix2world能够使用TRi_j畸变系数，所以我可以用Maxim的DL数据工作？

非常感谢

奈杰尔

Answer 1

除了SIP系数，该astropy.wcs方法将与TPV distortion coefficients工作。这意味着你可以直接与astropy.wcs使用SCAMP astrometric solver的输出。如果你想给冠捷系数转换成SIP形式，可以使用对此我带头贡献者sip_tpv package。我不知道一个Python包包装SCAMP的 - 我已经把它包在Zwicky Transient Facility pipeline但代码不公开。

Answer 2

你可以这样做：

from astropy.io import fits
from astropy.wcs import WCS

hdul = fits.open(fitsfilename)[0]

wcs = WCS(hdul.header)

ax = fig.gca()
ax.scatter([34], [3.2], transform=ax.get_transform('world'))

（Qazxswpoi。）