以下の python スクリプトを使用して、世界座標系 (WCS) を fit ファイルから読み取り、プロットします。
from astropy.io import fits
import matplotlib.pyplot as plt
from astropy.wcs import WCS
hdu = fits.open('file.fits')
header = hdu[0].header
data = hdu[0].data
wcs = WCS(header, naxis=2)
plt.figure().add_subplot(1,1,1, projection=wcs)
plt.imshow(data)
plt.grid()
plt.show()
ただし、上記のスクリプトは画像を RA と DEC に投影します。方位角と仰角で画像を表示する方法があるかどうか疑問に思っていました。
ありがとう。