Проблема с возвратом WCS из образа FITS

Я использую astropy для загрузки изображения FITS и извлечения из него WCS.

from astropy.io import fits
from astropy.wcs import WCS

with fits.open('hst_A2744_f606w_drz.fits') as hdul:
    wcs = WCS(hdul[1])
    print(wcs)

Что возвращает

WCS Keywords

Number of WCS axes: 2
CTYPE : 'RA---TAN'  'DEC--TAN'
CRVAL : 3.587755699764648  -30.39711750881429
CRPIX : 3000.4999999998081  2989.499999999809
CD1_1 CD1_2  : -1.3888888888888e-05  0.0
CD2_1 CD2_2  : 0.0  1.3888888888889599e-05
NAXIS : 6000  5978

Моя цель — вернуть CRVAL в виде строки, кортежа или массива, например: (3.587755699764648, -30.39711750881429)

Я пытался получить к нему доступ как wcs['CTYPE'] и wcs[CTYPE], и wcs.CTYPE, но все возвращают ошибки.


python astropy wcs
person theastronomist    schedule 04.12.2020    source источник

Ответы (1)


arrow_upward
1
arrow_downward

Вы хотите WCS.wcs

и, в частности, WCS.wcs.crval (или WCS.wcs.ctype, поскольку в своем вопросе вы написали, что вам нужны значения CRVAL, но вы показали попытки найти значения CTYPE).

person Iguananaut    schedule 07.12.2020