Problem with FXCOR
Pasit wrote on Apr 08, 2009
Hello,
I downloaded spectral data from ELODIE archive and tried to find radial velocity by FXCOR. But, I received an error message.
ERROR: Cannot convert to specified units
Then I tried to solve this problem by editing the header.
This is the header of my data.
noao> imhead /home/thagoon/Desktop/51peg-t1.fit l+
/home/thagoon/Desktop/51peg-t1.fit[56000][real]: HD217014
No bad pixels, min=0., max=0. (old)
Line storage mode, physdim [56000], length of user area 5630 s.u.
Created Mon 15:01:59 06-Apr-2009, Last modified Mon 15:01:59 06-Apr-2009
Pixel file "/home/thagoon/Desktop/51peg-t1.fit" [ok]
EXTEND = T / File may contain extensions
ORIGIN = 'NOAO-IRAF FITS Image Kernel July 2003' / FITS file originator
DATE = '2009-04-06T08:01:59' / Date FITS file was generated
IRAF-TLM= '15:01:59 (06/04/2009)' / Time of last modification
OBJECT = 'HD217014' / Name of the object observed
COMMENT FITS (Flexible Image Transport System) format is defined in 'Astronomy
COMMENT and Astrophysics', volume 376, page 359; bibcode: 2001A&A...376..359H
COMMENT FITS (Flexible Image Transport System) format is defined in 'Astronomy
COMMENT and Astrophysics', volume 376, page 359; bibcode: 2001A&A...376..359H
COMMENT FITS (Flexible Image Transport System) format is defined in 'Astronomy
COMMENT and Astrophysics', volume 376, page 359; bibcode: 2001A&A...376..359H
CRPIX1 = 1. / Reference pixel
CRVAL1 = 4000. / Coordinate at reference pixel
CDELT1 = 0.0500000007 / Coordinate increment (prefer CD)
IMATYP = 'OBTH/S2D/CAL' / TACOS_IMATYP
VERSION = '1.3 ' / TACOS_VERSION
LASTCOMP= 'Fri Jun 24 13:05:00 MET DST 1994' / TACOS_LASTCOMP
NOCCD = 1.00000000E+00 / ELODIE_NOCCD
NOCONFIG= 1.00000000E+00 / ELODIE_NOCONFIG
CCDGAIN = 2.00000000E+00 / ELODIE_CCDGAIN
IMABIAS = 'n19940914/bias0001' / ELODIE_IMABIAS
TEMP1 = -9.00400009E+01 / ELODIE_TEMP
TEMP2 = 2.08999996E+01 / ELODIE_TEMP
TEMP3 = 2.12000008E+01 / ELODIE_TEMP
TEMP4 = 2.12000008E+01 / ELODIE_TEMP
TEMP5 = 1.16000004E+01 / ELODIE_TEMP
TEMP6 = 1.16000004E+01 / ELODIE_TEMP
TEMP7 = 1.16000004E+01 / ELODIE_TEMP
TEMP8 = 1.16000004E+01 / ELODIE_TEMP
NOCORREF= 2.00000000E+00 / ELODIE_NOCORREF
ANGCREF = 3.30000000E+02 / ELODIE_ANGCREF
NONUIT = 2.20000000E+01 / NONUIT
NOETER = 3.80500000E+03 / NOETER
NOOBJET1= 1 / CAT_NOOBJET
NOOBJET2= 0 / CAT_NOOBJET
NOOBJET3= 217014 / CAT_NOOBJET
NOOBJET4= 0 / CAT_NOOBJET
INSTRUME= 'ELODIE ' / INSTRUME
CUMUL = 1.00000000E+00 / O_CUMUL
TIME1 = 1.99400000E+03 / O_TIME
TIME2 = 9.00000000E+00 / O_TIME
TIME3 = 1.50000000E+01 / O_TIME
TIME4 = 5.18041968E-01 / O_TIME
TIME5 = 3.01760010E+02 / O_TIME
POS1 = 2.29527779E+01 / TEL_POS
POS2 = 2.07380562E+01 / TEL_POS
OBSERVER= 'quidam ' / OBSERVER
REFNOCOD= '* ' / CAT_REFNOCODE
EQUINOX = 1.95000000E+03 / CAT_EQUINOX
ALPHA = '22:55:00' / CAT_ALPHA
DELTA = '+20:30:00' / CAT_DELTA
PROPERM1= 2.02999994E-01 / CAT_PROPERMVT
PROPERM2= 4.69999984E-02 / CAT_PROPERMVT
COMMENT = '* ' / COMMENT
NOPROG = 104 / CAT_NOPROG
IMALOC = 'n19940914/s2d/loco0017' / TACOS_IMALOC
NORDER = 6.70000000E+01 / TACOS_NORDER
DECLOC = 0.00000000E+00 / TACOS_DECLOC
DEGLOC1 = 4.00000000E+00 / TACOS_DEGLOC
DEGLOC2 = 2.50000000E+01 / TACOS_DEGLOC
DEGSIG1 = 3.00000000E+00 / TACOS_DEGSIG
DEGSIG2 = 4.00000000E+00 / TACOS_DEGSIG
OLIMLS = 1.00000000E+00 / TACOS_OLIMLS
EXTSTR = 3.00000000E+00 / TACOS_EXTSTR
JDB1 = 2.44961000E+06 / TACOS_JDB
JDB2 = 5.26821136E-01 / TACOS_JDB
BERV = 4.90955710E-01 / TACOS_BERV
BERVMX = 2.89625320E+01 / TACOS_BERVMX
DELCAL = 1.81319099E-03 / TACOS_DELCAL
DEGXLL = 4.00000000E+00 / TACOS_DEGXLL
DEGOLL = 5.00000000E+00 / TACOS_DEGOLL
COELL25 = 4.1849897120096E-01 / TACOS_COELL
COELL26 = 5.0761034225724E-02 / TACOS_COELL
COELL27 = 2.6488538167900E-02 / TACOS_COELL
COELL28 = 8.0478483402008E-02 / TACOS_COELL
COELL29 = 4.2389516158353E-02 / TACOS_COELL
COELL30 = 6.0074712830782E-02 / TACOS_COELL
OLIMLL = 1.00000000E+00 / TACOS_OLIMLL
LLOFFSET= 1.57015152E+02 / TACOS_LLOFFSET
LLPARS = -1.00000000E+00 / TACOS_LLPARS
LLCMOY = 0.00000000E+00 / TACOS_LLCMOY
LLCSIG = 0.00000000E+00 / TACOS_LLCSIG
IMACAL = 'n19940914/llo0021' / TACOS_IMACAL
LLCNRA = 1.05000000E+03 / TACOS_LLCNRA
FILERC = '/home/queloz/src/tacos/elodie/elodiehp' / TACOS_FILERC
PROXIM = 2.50000000E-01 / TACOS_PROXIM
FILENAME= 'elodie:19940914/0022' / Original filename
CD1_1 = 0.0500000007 / Coordinate increment per pixel
CREATOR = 'Pleinpot 1' / Written by Pleinpot (Hypercat)
DATE-OBS= '1994-09-15T00:31:04' / [CCYY-MM-DDThh:mm:ss] Date observation
MJD-OBS = 49610.021585 / [d] Modified JD when observ. start
AIRMASS = 1.16051221 / Airmass
CTYPE1 = 'LINEAR ' / Axis type
EXPTIME = 301.76001 / [s] Exposure time
H_CONTXT= 'SPE ' / Spectrum
H_STATUS= 0 / No error on pipeline
COMMENT The ELODIE archive (OHP)
COMMENT Note: The telescope position (POS1,POS2) and the catalogue information
COMMENT (NOOBETn)
COMMENT may not be accurate. Please use instead the designation from keyword
COMMENT (OBJECT)
H_AFR000= 'elodie:19940914/0022' / Archive identifier of current frame
H_TYPE = 1 / Science object
SN = 65.1 / S/N-per-pixel of original spectrum at 550 nm
DIVUL = 0 / data fully public
CUNIT1 = '0.1 nm ' / Unit of coordinate
BTYPE = 'FLUX-INS' / Type of data
BUNIT = 'instrumental' / Unit of data
H_DBIAS = 0. / [ADU] Bias on the frame
H_WRESOL= 0.1321429 / [0.1nm] FWHM wavelength resolution
H_ASP002= 'fa:LL_ELODIE/99002' / File used for "strengthening"
H_HISTO = 'FFC000/WCA000/'
H_FCANOR= '[FCANOR]' / Pointer to the Flux calibration relation
APNUM1 = '1 1 '
DC-FLAG = 0
WCSDIM = 1
LTM1_1 = 1.
WAT0_001= 'system=equispec'
WAT1_001= 'wtype=linear'
These are my questions:
1. Could you tell me whether my solution is right or wrong?
2. How can I solve it?
Thank a lot,
Pasit
I downloaded spectral data from ELODIE archive and tried to find radial velocity by FXCOR. But, I received an error message.
ERROR: Cannot convert to specified units
Then I tried to solve this problem by editing the header.
This is the header of my data.
noao> imhead /home/thagoon/Desktop/51peg-t1.fit l+
/home/thagoon/Desktop/51peg-t1.fit[56000][real]: HD217014
No bad pixels, min=0., max=0. (old)
Line storage mode, physdim [56000], length of user area 5630 s.u.
Created Mon 15:01:59 06-Apr-2009, Last modified Mon 15:01:59 06-Apr-2009
Pixel file "/home/thagoon/Desktop/51peg-t1.fit" [ok]
EXTEND = T / File may contain extensions
ORIGIN = 'NOAO-IRAF FITS Image Kernel July 2003' / FITS file originator
DATE = '2009-04-06T08:01:59' / Date FITS file was generated
IRAF-TLM= '15:01:59 (06/04/2009)' / Time of last modification
OBJECT = 'HD217014' / Name of the object observed
COMMENT FITS (Flexible Image Transport System) format is defined in 'Astronomy
COMMENT and Astrophysics', volume 376, page 359; bibcode: 2001A&A...376..359H
COMMENT FITS (Flexible Image Transport System) format is defined in 'Astronomy
COMMENT and Astrophysics', volume 376, page 359; bibcode: 2001A&A...376..359H
COMMENT FITS (Flexible Image Transport System) format is defined in 'Astronomy
COMMENT and Astrophysics', volume 376, page 359; bibcode: 2001A&A...376..359H
CRPIX1 = 1. / Reference pixel
CRVAL1 = 4000. / Coordinate at reference pixel
CDELT1 = 0.0500000007 / Coordinate increment (prefer CD)
IMATYP = 'OBTH/S2D/CAL' / TACOS_IMATYP
VERSION = '1.3 ' / TACOS_VERSION
LASTCOMP= 'Fri Jun 24 13:05:00 MET DST 1994' / TACOS_LASTCOMP
NOCCD = 1.00000000E+00 / ELODIE_NOCCD
NOCONFIG= 1.00000000E+00 / ELODIE_NOCONFIG
CCDGAIN = 2.00000000E+00 / ELODIE_CCDGAIN
IMABIAS = 'n19940914/bias0001' / ELODIE_IMABIAS
TEMP1 = -9.00400009E+01 / ELODIE_TEMP
TEMP2 = 2.08999996E+01 / ELODIE_TEMP
TEMP3 = 2.12000008E+01 / ELODIE_TEMP
TEMP4 = 2.12000008E+01 / ELODIE_TEMP
TEMP5 = 1.16000004E+01 / ELODIE_TEMP
TEMP6 = 1.16000004E+01 / ELODIE_TEMP
TEMP7 = 1.16000004E+01 / ELODIE_TEMP
TEMP8 = 1.16000004E+01 / ELODIE_TEMP
NOCORREF= 2.00000000E+00 / ELODIE_NOCORREF
ANGCREF = 3.30000000E+02 / ELODIE_ANGCREF
NONUIT = 2.20000000E+01 / NONUIT
NOETER = 3.80500000E+03 / NOETER
NOOBJET1= 1 / CAT_NOOBJET
NOOBJET2= 0 / CAT_NOOBJET
NOOBJET3= 217014 / CAT_NOOBJET
NOOBJET4= 0 / CAT_NOOBJET
INSTRUME= 'ELODIE ' / INSTRUME
CUMUL = 1.00000000E+00 / O_CUMUL
TIME1 = 1.99400000E+03 / O_TIME
TIME2 = 9.00000000E+00 / O_TIME
TIME3 = 1.50000000E+01 / O_TIME
TIME4 = 5.18041968E-01 / O_TIME
TIME5 = 3.01760010E+02 / O_TIME
POS1 = 2.29527779E+01 / TEL_POS
POS2 = 2.07380562E+01 / TEL_POS
OBSERVER= 'quidam ' / OBSERVER
REFNOCOD= '* ' / CAT_REFNOCODE
EQUINOX = 1.95000000E+03 / CAT_EQUINOX
ALPHA = '22:55:00' / CAT_ALPHA
DELTA = '+20:30:00' / CAT_DELTA
PROPERM1= 2.02999994E-01 / CAT_PROPERMVT
PROPERM2= 4.69999984E-02 / CAT_PROPERMVT
COMMENT = '* ' / COMMENT
NOPROG = 104 / CAT_NOPROG
IMALOC = 'n19940914/s2d/loco0017' / TACOS_IMALOC
NORDER = 6.70000000E+01 / TACOS_NORDER
DECLOC = 0.00000000E+00 / TACOS_DECLOC
DEGLOC1 = 4.00000000E+00 / TACOS_DEGLOC
DEGLOC2 = 2.50000000E+01 / TACOS_DEGLOC
DEGSIG1 = 3.00000000E+00 / TACOS_DEGSIG
DEGSIG2 = 4.00000000E+00 / TACOS_DEGSIG
OLIMLS = 1.00000000E+00 / TACOS_OLIMLS
EXTSTR = 3.00000000E+00 / TACOS_EXTSTR
JDB1 = 2.44961000E+06 / TACOS_JDB
JDB2 = 5.26821136E-01 / TACOS_JDB
BERV = 4.90955710E-01 / TACOS_BERV
BERVMX = 2.89625320E+01 / TACOS_BERVMX
DELCAL = 1.81319099E-03 / TACOS_DELCAL
DEGXLL = 4.00000000E+00 / TACOS_DEGXLL
DEGOLL = 5.00000000E+00 / TACOS_DEGOLL
COELL25 = 4.1849897120096E-01 / TACOS_COELL
COELL26 = 5.0761034225724E-02 / TACOS_COELL
COELL27 = 2.6488538167900E-02 / TACOS_COELL
COELL28 = 8.0478483402008E-02 / TACOS_COELL
COELL29 = 4.2389516158353E-02 / TACOS_COELL
COELL30 = 6.0074712830782E-02 / TACOS_COELL
OLIMLL = 1.00000000E+00 / TACOS_OLIMLL
LLOFFSET= 1.57015152E+02 / TACOS_LLOFFSET
LLPARS = -1.00000000E+00 / TACOS_LLPARS
LLCMOY = 0.00000000E+00 / TACOS_LLCMOY
LLCSIG = 0.00000000E+00 / TACOS_LLCSIG
IMACAL = 'n19940914/llo0021' / TACOS_IMACAL
LLCNRA = 1.05000000E+03 / TACOS_LLCNRA
FILERC = '/home/queloz/src/tacos/elodie/elodiehp' / TACOS_FILERC
PROXIM = 2.50000000E-01 / TACOS_PROXIM
FILENAME= 'elodie:19940914/0022' / Original filename
CD1_1 = 0.0500000007 / Coordinate increment per pixel
CREATOR = 'Pleinpot 1' / Written by Pleinpot (Hypercat)
DATE-OBS= '1994-09-15T00:31:04' / [CCYY-MM-DDThh:mm:ss] Date observation
MJD-OBS = 49610.021585 / [d] Modified JD when observ. start
AIRMASS = 1.16051221 / Airmass
CTYPE1 = 'LINEAR ' / Axis type
EXPTIME = 301.76001 / [s] Exposure time
H_CONTXT= 'SPE ' / Spectrum
H_STATUS= 0 / No error on pipeline
COMMENT The ELODIE archive (OHP)
COMMENT Note: The telescope position (POS1,POS2) and the catalogue information
COMMENT (NOOBETn)
COMMENT may not be accurate. Please use instead the designation from keyword
COMMENT (OBJECT)
H_AFR000= 'elodie:19940914/0022' / Archive identifier of current frame
H_TYPE = 1 / Science object
SN = 65.1 / S/N-per-pixel of original spectrum at 550 nm
DIVUL = 0 / data fully public
CUNIT1 = '0.1 nm ' / Unit of coordinate
BTYPE = 'FLUX-INS' / Type of data
BUNIT = 'instrumental' / Unit of data
H_DBIAS = 0. / [ADU] Bias on the frame
H_WRESOL= 0.1321429 / [0.1nm] FWHM wavelength resolution
H_ASP002= 'fa:LL_ELODIE/99002' / File used for "strengthening"
H_HISTO = 'FFC000/WCA000/'
H_FCANOR= '[FCANOR]' / Pointer to the Flux calibration relation
APNUM1 = '1 1 '
DC-FLAG = 0
WCSDIM = 1
LTM1_1 = 1.
WAT0_001= 'system=equispec'
WAT1_001= 'wtype=linear'
These are my questions:
1. Could you tell me whether my solution is right or wrong?
2. How can I solve it?
Thank a lot,
Pasit
Mike Fitzpatrick wrote on Apr 08, 2009
As far as the spectral WCS goes, I think what you need to do is change the CUNIT1 keywords to be 'angstroms' and add a DC-FLAG keyword to the image you get from ELODIE. For example,
The headers are also missing the RA/DEC/EPOCH keywords, but you may be able to use the KEYWPARS psets to use alternate keywords in the header (e.g. ALPHA/DELTA for RA/DEC and add an EPOCH keyword based on the DATE-OBS). The APNUM and WAT keywords shouldn't be needed.
-Mike
cl> hedit spec cunit1 angstroms
cl> hedit spec dc-flag 0 add+
The headers are also missing the RA/DEC/EPOCH keywords, but you may be able to use the KEYWPARS psets to use alternate keywords in the header (e.g. ALPHA/DELTA for RA/DEC and add an EPOCH keyword based on the DATE-OBS). The APNUM and WAT keywords shouldn't be needed.
-Mike
Pasit wrote on Apr 08, 2009
First,We thank you for your answer,
but!!! I have new problem in FXCOR again.
When I use the reaI template from SDSS, can't open "IRAFterm" because of the error.
error: Aperture list cannot be applied to template image.
Could you help me please?
but!!! I have new problem in FXCOR again.
When I use the reaI template from SDSS, can't open "IRAFterm" because of the error.
error: Aperture list cannot be applied to template image.
rv> fxcor
List of object spectra (/home/thagoon/Desktop/51peg-t2.fit):
List of template spectra (/home/thagoon/Desktop/spDR2-008.fit):
Aperture list cannot be applied to template image.
kpno: Kitt Peak National Observatory
ctio: Cerro Tololo Interamerican Observatory
lasilla: European Southern Observatory: La Silla
paranal: European Southern Observatory: Paranal
lick: Lick Observatory
mmto: MMT Observatory
cfht: Canada-France-Hawaii Telescope
lapalma: Roque de los Muchachos, La Palma
mso: Mt. Stromlo Observatory
sso: Siding Spring Observatory
aao: Anglo-Australian Observatory
mcdonald: McDonald Observatory
lco: Las Campanas Observatory
mtbigelow: Catalina Observatory: 61 inch telescope
dao: Dominion Astrophysical Observatory
spm: Observatorio Astronomico Nacional, San Pedro Martir
tona: Observatorio Astronomico Nacional, Tonantzintla
Palomar: The Hale Telescope
mdm: Michigan-Dartmouth-MIT Observatory
NOV: National Observatory of Venezuela
bmo: Black Moshannon Observatory
BAO: Beijing XingLong Observatory
keck: W. M. Keck Observatory
ekar: Mt. Ekar 182 cm. Telescope
apo: Apache Point Observatory
lowell: Lowell Observatory
vbo: Vainu Bappu Observatory
iao: Indian Astronomical Observatory, Hanle
flwo: Whipple Observatory
flwo1: Whipple Observatory
oro: Oak Ridge Observatory
lna: Laboratorio Nacional de Astrofisica - Brazil
saao: South African Astronomical Observatory
casleo: Complejo Astronomico El Leoncito, San Juan
bosque: Estacion Astrofisica Bosque Alegre, Cordoba
rozhen: National Astronomical Observatory Rozhen - Bulgaria
irtf: NASA Infrared Telescope Facility
bgsuo: Bowling Green State Univ Observatory
DSAZ: Deutsch-Spanisches Observatorium Calar Alto - Spain
ca: Calar Alto Observatory
holi: Observatorium Hoher List (Universitaet Bonn) - Germany
lmo: Leander McCormick Observatory
fmo: Fan Mountain Observatory
whitin: Whitin Observatory, Wellesley College
osn: Observatorio de Sierra Nevada
gemini-north: Gemini North Observatory
gemini-south: Gemini South Observatory
lasilla: European Southern Observatory: La Silla
paranal: European Southern Observatory: Paranal
esontt: European Southern Observatory, NTT, La Silla
eso36m: European Southern Observatory, 3.6m Telescope, La Silla
esovlt: European Southern Observatory, VLT, Paranal
sln: SLN - Catania Astrophysical Observatory
euo: Ege University Observatory
tug: TUBITAK National Observatory, Turkey
mgo: Mount Graham Observatory
aries: Aryabhatta Research Institute of Observational Sciences
oalp: Observatorio Astronómico de La Plata
olin: Connecticut
boyden: Boyden
obspars: Use parameters from OBSERVATORY task
Observatory identification (dao):
Could you help me please?
Mike Fitzpatrick wrote on Apr 08, 2009
The task assumes the object and template data were taken with the same instrument and so only provides a single aperture parameter that is applied to both spectra. It would make sense that having a 1-D template should be allowed against a 2-D object, but that would have to be a future addition. For now, use the SCOPY task to copy out the apertures to a new image and leave the apertures parameter at '*' to use everything.
However, the ELODIE data I got were just 1-D. Did you do something else to the header or else how did you set the aperture parameter?
The other output comes because it can't make sense of the observatory name, either the one you set in the 'observatory' parameter or one named in the OBSERVAT keyword. For Sloan spectra, and appropriate OBSERVAT keyword will be 'apo', for the ELODIE data I don't think we have the location in the observatory database. See the OBSERVATORY help page for how to add a new location to the database.
-Mike
However, the ELODIE data I got were just 1-D. Did you do something else to the header or else how did you set the aperture parameter?
The other output comes because it can't make sense of the observatory name, either the one you set in the 'observatory' parameter or one named in the OBSERVAT keyword. For Sloan spectra, and appropriate OBSERVAT keyword will be 'apo', for the ELODIE data I don't think we have the location in the observatory database. See the OBSERVATORY help page for how to add a new location to the database.
-Mike
Last post on Apr 08, 2009