我使用下面的python脚本读取形成配合文件的世界坐标系(WCS)并将其绘制。
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中。我想知道是否有一种方法可以使它以方位角和高程显示图像。
谢谢。