Как построить изображение из файла FITS с реальными координатами?

Question

В документации по астропии есть пример того, как построить изображение из файла FITS.Оси на полученном графике соответствуют количеству пикселей или точек данных.Я хотел бы, чтобы ось вместо этого указывала галактические координаты.

Обратите внимание, что заголовок моего файла FITS показывает всю информацию, касающуюся положения изображения на небе.Как получить на графике реальные координаты, а не количество пикселей?

Answer 1

Каркас WCS, предоставляемый 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()