我在像素坐标中有FITS文件。我想基本上阅读其纬度和经度。
我知道我们可以使用以下方法:
from astropy import WCS
w = WCS('sample.fits')
lat, long = w.all_pix2world(90, 38, 1)
但是DEC(= 38)值不是一个RA(= 90),而是
w.all_pix2world
我想获取 sample.fits 整个字段的纬度和经度(在银河系坐标系中)。
header
NAXIS = 2
NAXIS1 = 150
NAXIS2 = 150
请让我知道是否需要其他信息。任何帮助表示赞赏。
谢谢。