如何把fits文件的pix坐标和world坐标相互转换呢,用astropy可以有两种方法来进行转换
from astropy.wcs import WCS
from astropy.io import fits as pf
data = pf.getdata('M31_rotated_cal.fits')
hdr = pf.getheader('M31_rotated_cal.fits')
c = SkyCoord('0h42m48s', '+41d16m48s')
w_n = WCS(hdr)
print(w_n.world_to_pixel(c))
print(w_n.wcs_pix2world([[251, 200]], 0))
然后你把world_to_pixel和wcs_pix2world中的pix和word颠倒一下也是可以用的,两种方式输出和输入的的格式不一样,可以print出来看看格式