如何使用真实坐标绘制FITS文件中的图像？

Question

astropy documentation有一个关于如何从FITS文件中绘制图像的示例。结果图上的轴对应于像素数或数据点。我希望轴代表银河坐标。

请注意，我的FITS文件的标题表示有关天空中图像位置的所有信息。如何让图表显示真实坐标而不是像素数？

Answer 1

WCS framework provided by astropy可用于绘制世界坐标。

from astropy.wcs import WCS
from astropy.io import fits
from astropy.utils.data import get_pkg_data_filename
import matplotlib.pyplot as plt

image_file = get_pkg_data_filename('tutorials/FITS-images/HorseHead.fits')
hdu = fits.open(image_file)[0]
wcs = WCS(hdu.header)

plt.subplot(projection=wcs) 
plt.imshow(hdu.data, origin='lower') 
plt.grid(color='white', ls='solid')
plt.show()