SIP兼容Python板解包答案

【问题标题】：SIP-compatible Python plate solving packageSIP兼容Python板解包
【发布时间】：2019-01-09 16:06:42
【问题描述】：

我需要对最初使用 Maxim DL 创建的 FITS 文件执行精确的像素到世界坐标转换。 Maxim 使用 Pinpoint 进行板求解，生成 TRi_j 失真系数。这些与我建议使用的 astropy.wcs 坐标变换函数不兼容，因为它们假设 SIP 失真系数。

因此，我正在寻找对 FITS 文件进行重新板解以生成 SIP 系数的选项。

到目前为止，我发现的只是 astrometry.net，但这是一项在线服务。我真的在寻找可以作为我的应用工作流程的一部分同步执行的离线溶版（最好针对 GSC 的本地副本）。

是否有任何 Astropy 附属（或其他）Python 包可以针对 GSC 执行 SIP 兼容的溶板？

或者，是否有任何与 wcs.all_pix2world 等效的可以使用 TRi_j 失真系数的方法，以便我可以使用 Maxim DL 数据？

非常感谢

奈杰尔

【问题讨论】：

标签： astropy

【解决方案1】：

除了 SIP 系数之外，astropy.wcs 方法还适用于 TPV distortion coefficients。这意味着您可以将SCAMP astrometric solver 的输出直接与astropy.wcs 一起使用。如果您希望将 TPV 系数转换为 SIP 形式，您可以使用我是主要贡献者的sip_tpv package。我不知道包装 SCAMP 的 Python 包——我已经为 Zwicky Transient Facility pipeline 包装了它，但该代码不公开。

【讨论】：

【解决方案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'))

(Based on this Q.)

【讨论】：