我需要对最初使用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数据？
非常感谢
奈杰尔
答案 0 :(得分：0)
除SIP系数外，astropy.wcs方法还将与TPV distortion coefficients一起使用。这意味着您可以将SCAMP astrometric solver的输出直接用于astropy.wcs。如果您希望将TPV系数转换为SIP格式，则可以使用我作为主要撰稿人的sip_tpv package。我不知道有一个包装SCAMP的Python软件包-我已经将其包装为Zwicky Transient Facility pipeline，但该代码不是公开的。
答案 1 :(得分：0)
您可以这样做：
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'))