【发布时间】:2020-12-04 05:42:29
【问题描述】:
我正在使用 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 访问它,但都返回错误。
【问题讨论】: