1
我想繪製擬合與aplpy文件適合圖像與aplpy但我得到以下錯誤:繪製使用正確的WCS改造
INFO: The WCS transformation has more axes (3) than the image it is associated with (2) [astropy.wcs.wcs]
ERROR: IndexError: list index out of range [aplpy.wcs_util]
如果我刪除標題關鍵詞CTYPE3,CD3_3,WAT3_001然後我再也看不到這個錯誤了,但是當繪製合適的文件時,我發現WCS轉換並沒有在RA和DEC座標上正確完成。
這是標題:
HCG89_halpha_plot.fits[2929,2881][real]: HCG89-SII
No bad pixels, min=0., max=0. (old)
Line storage mode, physdim [2929,2881], length of user area 5144 s.u.
Created Tue 15:40:54 28-Jan-2014, Last modified Thu 22:20:03 30-Jan-2014
Pixel file "HCG89_halpha_plot.fits" [ok]
ORIGIN = 'NOAO-IRAF FITS Image Kernel July 2003'/FITS file originator
IRAF-TLM= '2014-01-30T19:20:03'/Time of last modification
OBJECT = 'HCG89-SII' /Name of the object observed
DATE-OBS= '2013-10-11T00:14:11.700'/Date of observation
DATE = '2014-01-28T12:40:54'/Date Format is YYYY-MM-DD
TIME = '00:14:11.07 to 00:19:11.40'/~ Start & Stop of Exposure
N_PARAM = 80/Number of Parameters
TELESCOP= 'SOAR 4.1m'
INSTRUME= 'Goodman Spectrograph'
NOTES = ' '
RA = '21:20:09.042' /right ascension [hh:mm:ss.sss]
DEC = '-3:55:23.160' /declination [dd:mm:ss.sss]
AIRMASS = 1.13/airmass at approx. start of exposure
UT = '00:14:10.9' /time at approx. start of exposure [UTC]
FOCUS = -1021.0/SOAR telescope focus
MOUNT_AZ= 17.1482/SOAR mount azimuth
MOUNT_EL= 62.2622/SOAR mount elevation
ROTATOR = 228.003/Nasymth cage rotator angle [deg]
POSANGLE= 360.0/position angle [deg. E of N]
SEEING = -1.0/DIMM seeing [arcsec]
LST = '20:50:17.4' /Local Sidereal Time [hh:mm:ss.sss]
OBSRA = '21:20:09.946' /target right ascension [hh:mm:ss.sss]
OBSDEC = '-3:55:32.480' /target declination [hh:mm:ss.sss]
CAM_ANG = 0.014461/camera angle [deg]
GRT_ANG = 0.661/grating angle [deg]
CAM_TARG= 0.0/camera target [deg]
GRT_TARG= 0.661/grating target [deg]
CAM_FOC = 660/camera focus
COLL_FOC= 1002/collimator focus
FILTER = 'SII ' /primary filter wheel
FILTER2 = '<NO FILTER>' /secondary filter wheel
GRATING = '<NO GRATING>' /VPH grating [1/mm]
SLIT = '<NO MASK>' /slit [arcsec]
COL_TEMP= 10.421875/coll ext temp (deg C)
CAM_TEMP= 14.234375/cam ext temp (deg C)
EXPTIME = 300.0/integration time
RDNOISE = 4.74/CCD readnoise [e-]
GAIN = 1.4/CCD gain [e-/ADU]
OBSTYPE = 'OBJECT ' /observation type
OBSERVER= ' '
PROPOSAL= ' '
EQUINOX = 2000.0/equinox of coordinates
CRPIX1 = 1477.39599990845/Reference pixel on axis 1
CRPIX2 = 1488.23082170899/Reference pixel on axis 2
CRVAL1 = 320.041401353/Value at ref. pixel on axis 1
CRVAL2 = -3.91652009277/Value at ref. pixel on axis 2
CTYPE1 = 'RA---TAN' /Type of co-ordinate on axis 1
CTYPE2 = 'DEC--TAN' /Type of co-ordinate on axis 2
PC1_1 = 0.999884048946387/Transformation matrix element
PC1_2 = 0.0152278909432213/Transformation matrix element
PC2_1 = -0.0140395410920601/Transformation matrix element
PC2_2 = 1.0/Transformation matrix element
MJD-OBS = 56576.0098576389/Modified Julian Date of observation
RADESYS = 'FK5 ' /Reference frame for RA/DEC values
DISPAXIS= 1
DETSIZE = '[1:4096,1:4096]'
TRIMSEC = '[1:3096,1:3096]'
CCDSIZE = '[1:4096,1:4096]'
CCDSUM = '1 1 '
OPENTIME= '00:14:11.611' /GPS-Synched Time for Shutter Open (UT)
INSTRUME= 'Goodman Spectrograph'
COMMENT Triggered Acquisition, Exp Time= 05:00
WCSDIM = 3
LTM1_1 = 1.0
LTM2_2 = 1.0
LTM3_3 = 1.0
WAXMAP01= '1 0 2 0 0 0 '
WAT0_001= 'system=physical'
WAT1_001= 'wtype=linear axtype=ra'
WAT2_001= 'wtype=linear axtype=dec'
WAT3_001= 'wtype=linear'
LTV1 = -159.604000091553
LTV2 = -134.766998291016
IMCMB001= '0073.HCG89-SII_scaled.fits'
IMCMB002= '0074.HCG89-SII_scaled.fits'
IMCMB003= '0079.HCG89-SII_scaled.fits'
IMCMB004= '0080.HCG89-SII_scaled.fits'
IMCMB005= '0083.HCG89-SII_scaled.fits'
IMCMB006= '0084.HCG89-SII_scaled.fits'
NCOMBINE= 6
CTYPE3 = 'LINEAR '
CD1_1 = 4.03250815042544E-5
CD2_2 = 4.04084336509694E-5
CD3_3 = 1.