Найти физические координаты пикселя в файле соответствия с питоном
Я пытаюсь получить физические координаты неба данного пикселя из скрипта Python. Я хотел бы использовать WCS от Astropy, но я сделаю что-нибудь из Python.
Я попробовал эти два фрагмента кода.
from astropy.io import fits
from astropy.wcs import WCS
def astropymethod1(img):
# from http://astropy.readthedocs.org/en/latest/wcs/
w = WCS(img)
lon, lat = w.all_pix2world( 100., 100., 1)
print lon, lat
def astropymethod2(img):
# from http://astropy.readthedocs.org/en/latest/wcs/
hdu = fits.open(img)
w = WCS(hdu[0].header)
lon, lat = w.wcs_pix2world(100., 100., 1)
print lon, lat
Проблемы заключаются в том, что я получаю сообщение об ошибке при первой попытке использовать WCS, и в результате получаем только те значения пикселей, которые я ввел.
WARNING: FITSFixedWarning: The WCS transformation has more axes (2) than the image it is associated with (0) [astropy.wcs.wcs]