problem with using msczero on a non mef image
paul wrote on Oct 25, 2006
I am having problems running msczero on a single image, ie non mosaic
image. MSCZERO gives the message "Cursor position not in an image"
for the "z" key. I beleive I have an approximately correct WCS in the
header, offset by a few arc minutes from reality.
The images are from the Burrell Schmidt, (single 2x4k CCD) I am trying
to use msczero to correct the WCS created by FINDER.TFINDER on one of
a set of dithered images and copied the headers of the other images.
I used IMCOORDS.CCSETWCS to copy the header.
When I run msczero the "u" key correctly displays the catalog location
of the stars, exactly in place for the image tfinder was run on. At
various offsets on the other images.
the "s" key reads the coordinates of the catalog star correctly
below is the output form msczero and the image header
-------------------------------------------------------
mscfinder> msczero
List of mosaic exposures (obj529):
Amp: Individual Display (zcombine=none)
im1: 507.0 620.8 507.0 620.8
>>u entered on the image displayed on DS9
List of coordinates ():
Magnitude limit (10.):
Radii of concentric circles (10):
Gray level of marks to be drawn (205):
Label the marked coordinates (no):
>>s entered at the tredicted location of bright
"Selected: 3:54:53.312 0:00:38.75 8.400"
>>z entered at matching star location
"Cursor position not in an image"
-------------------------------------------------------
mscfinder> msczero
Here is the image header,
obj352.fits[2046,4093][short]: None
No bad pixels, min=0., max=0. (old)
Line storage mode, physdim [2046,4093], length of user area 2673 s.u.
Created Sun 23:01:34 22-Oct-2006, Last modified Sun 23:01:38 22-Oct-2006
Pixel file "obj352.fits" [ok]
EXTEND = F / File may contain extensions
ORIGIN = 'NOAO-IRAF FITS Image Kernel July 2003' / FITS file originator
OBJECT = 'None ' / Name of the object observed
DATE = '2006-03-04T01:32:19' / Date FITS file was generated
IRAF-TLM= '23:01:38 (22/10/2006)' / Time of last modification
DATE-OBS= '2006-01-09T06:13:30.000' / Date and time of observation
TIMESYS = 'UTC ' / Default time system
OBSERVAT= 'Observatory' / Observatory
SIDETIME= 'hh:mm:ss.ss'
OBJEPOCH= 2000 / Epoch of object coordinates
AIRMASS = 99.99 / Airmass
EXPTIME = 30 / Exposure time
HA = 'hh:mm:ss.ss' / Hour angle
TELESCOP= 'Telescope' / Telescope
INSTRUME= 'Instrument' / Instrument
FILTER1 = 'None ' / Filter 1
FILTER2 = 'None ' / Filter 2
TELRA = 'hh:mm:ss.ss' / Telescope right ascension
TELDEC = 'hh:mm:ss.ss' / Telescope declination
OBSERVER= 'Observer' / Observer
COMMENT None
COMMENT None
COMMENT None
COMMENT None
FILTER = 'v '
WCSDIM = 2
LTM1_1 = 1.
LTM2_2 = 1.
WAT0_001= 'system=image'
WAT1_001= 'wtype=tan axtype=ra'
WAT2_001= 'wtype=tan axtype=dec'
TRIM = 'Jan 29 20:41 Trim data section is [29:2074,3:4095]'
OVERSCAN= 'Jan 29 20:41 Overscan section is [2100:2200,1:4200] with mean=2051.8'
ZEROCOR = 'Jan 29 20:41 Zero level correction image is ../Zero.fits'
FLATCOR = 'Jan 29 20:41 Flat field image is ../Flatv.fits with scale=17676.1'
CCDSEC = '[29:2074,3:4095]'
LTV1 = -28.
LTV2 = -2.
CCDMEANT= 823034487
CCDPROC = 'Jan 29 20:41 CCD processing done'
RADECSYS= 'FK5 '
EQUINOX = 2000.
CTYPE1 = 'RA---TAN'
CTYPE2 = 'DEC--TAN'
CRVAL1 = 58.7223124999995
CRVAL2 = 0.010894444444444
CRPIX1 = 1550.07738435592
CRPIX2 = 2243.89879673122
CD1_1 = 4.85762776557403E-5
CD1_2 = 4.03355977031450E-4
CD2_1 = -4.0340106805173E-4
CD2_2 = 4.86414776843567E-5
I also tried adding
DETSEC = '[1:2046,1:4093]'
DATASEC = '[1:2046,1:4093]'
IMAGEID = 1
to the header with hedit, but still get the same error message.
image. MSCZERO gives the message "Cursor position not in an image"
for the "z" key. I beleive I have an approximately correct WCS in the
header, offset by a few arc minutes from reality.
The images are from the Burrell Schmidt, (single 2x4k CCD) I am trying
to use msczero to correct the WCS created by FINDER.TFINDER on one of
a set of dithered images and copied the headers of the other images.
I used IMCOORDS.CCSETWCS to copy the header.
When I run msczero the "u" key correctly displays the catalog location
of the stars, exactly in place for the image tfinder was run on. At
various offsets on the other images.
the "s" key reads the coordinates of the catalog star correctly
below is the output form msczero and the image header
-------------------------------------------------------
mscfinder> msczero
List of mosaic exposures (obj529):
Amp: Individual Display (zcombine=none)
im1: 507.0 620.8 507.0 620.8
>>u entered on the image displayed on DS9
List of coordinates ():
Magnitude limit (10.):
Radii of concentric circles (10):
Gray level of marks to be drawn (205):
Label the marked coordinates (no):
>>s entered at the tredicted location of bright
"Selected: 3:54:53.312 0:00:38.75 8.400"
>>z entered at matching star location
"Cursor position not in an image"
-------------------------------------------------------
mscfinder> msczero
Here is the image header,
obj352.fits[2046,4093][short]: None
No bad pixels, min=0., max=0. (old)
Line storage mode, physdim [2046,4093], length of user area 2673 s.u.
Created Sun 23:01:34 22-Oct-2006, Last modified Sun 23:01:38 22-Oct-2006
Pixel file "obj352.fits" [ok]
EXTEND = F / File may contain extensions
ORIGIN = 'NOAO-IRAF FITS Image Kernel July 2003' / FITS file originator
OBJECT = 'None ' / Name of the object observed
DATE = '2006-03-04T01:32:19' / Date FITS file was generated
IRAF-TLM= '23:01:38 (22/10/2006)' / Time of last modification
DATE-OBS= '2006-01-09T06:13:30.000' / Date and time of observation
TIMESYS = 'UTC ' / Default time system
OBSERVAT= 'Observatory' / Observatory
SIDETIME= 'hh:mm:ss.ss'
OBJEPOCH= 2000 / Epoch of object coordinates
AIRMASS = 99.99 / Airmass
EXPTIME = 30 / Exposure time
HA = 'hh:mm:ss.ss' / Hour angle
TELESCOP= 'Telescope' / Telescope
INSTRUME= 'Instrument' / Instrument
FILTER1 = 'None ' / Filter 1
FILTER2 = 'None ' / Filter 2
TELRA = 'hh:mm:ss.ss' / Telescope right ascension
TELDEC = 'hh:mm:ss.ss' / Telescope declination
OBSERVER= 'Observer' / Observer
COMMENT None
COMMENT None
COMMENT None
COMMENT None
FILTER = 'v '
WCSDIM = 2
LTM1_1 = 1.
LTM2_2 = 1.
WAT0_001= 'system=image'
WAT1_001= 'wtype=tan axtype=ra'
WAT2_001= 'wtype=tan axtype=dec'
TRIM = 'Jan 29 20:41 Trim data section is [29:2074,3:4095]'
OVERSCAN= 'Jan 29 20:41 Overscan section is [2100:2200,1:4200] with mean=2051.8'
ZEROCOR = 'Jan 29 20:41 Zero level correction image is ../Zero.fits'
FLATCOR = 'Jan 29 20:41 Flat field image is ../Flatv.fits with scale=17676.1'
CCDSEC = '[29:2074,3:4095]'
LTV1 = -28.
LTV2 = -2.
CCDMEANT= 823034487
CCDPROC = 'Jan 29 20:41 CCD processing done'
RADECSYS= 'FK5 '
EQUINOX = 2000.
CTYPE1 = 'RA---TAN'
CTYPE2 = 'DEC--TAN'
CRVAL1 = 58.7223124999995
CRVAL2 = 0.010894444444444
CRPIX1 = 1550.07738435592
CRPIX2 = 2243.89879673122
CD1_1 = 4.85762776557403E-5
CD1_2 = 4.03355977031450E-4
CD2_1 = -4.0340106805173E-4
CD2_2 = 4.86414776843567E-5
I also tried adding
DETSEC = '[1:2046,1:4093]'
DATASEC = '[1:2046,1:4093]'
IMAGEID = 1
to the header with hedit, but still get the same error message.
Francisco Valdes wrote on Oct 25, 2006
Hello,
As a guess try
cl> set disable_wcs_maps=""
cl> flpr
This forces use of a simpler mapping between images and coordinates. Does this help?
Yours,
Frank Valdes
P.S. Make sure the image is displayed either by msczero itself or by mscdisplay. If you use display and then msczero then a file written to uparm needed by msczero will be absent or wrong and things will not work.
As a guess try
cl> set disable_wcs_maps=""
cl> flpr
This forces use of a simpler mapping between images and coordinates. Does this help?
Yours,
Frank Valdes
P.S. Make sure the image is displayed either by msczero itself or by mscdisplay. If you use display and then msczero then a file written to uparm needed by msczero will be absent or wrong and things will not work.
paul wrote on Oct 25, 2006
Thanks Frank,
You suggestion solved the problem.
MSCZERO works perfectly now.
Regards
Paul Harding
You suggestion solved the problem.
MSCZERO works perfectly now.
Regards
Paul Harding
Last post on Oct 25, 2006