Flatcombine
camjh wrote on Mar 28, 2008
Hi,
I have a problem with the Iraf FLATCOMBINE task, when trying to create a combined flat file of UVES flat files (ff lamp 1) I get errors. I have tried running it with different lists and by typing in the single files. Everytime I get an error:
List of flat field images to combine (flat615.fits, flat226.fits):
ERROR on line 47: No images to combine
flatcombine ()
But I am able to combine the flats using COMBINE, and ZEROCOMBINE also works..
To solve the problem I have so far tried the following:
-removing iraf directory and creating a new using mkiraf (commented a few
suggested lines in the login.cl file in)
-flatcombine *.fits ccdtype=""
-gave all rights to the fits files
-created a new directory, where I moved the flat fits files to, renamed them and gave all rights to them
-changed setinstrument (echelle,direct)
-changed ccdproc's biassec and trimsec ([values],image)
-unlearn flatcombine and ccdproc.
My task input look like this:
PACKAGE = ccdred
TASK = setinstrument
instrume= echelle Instrument ID (type ? for a list)
(site = kpno) Site ID
(directo= ccddb$) Instrument directory
(review = yes) Review instrument parameters?
query = echelle Instrument ID (type q to quit)
(mode = al)
PACKAGE = ccdred
TASK = ccdproc
images = flat615.fits, flat226.fits List of CCD images to correct
(output = flattest2.fits) List of output CCD images
(ccdtype= ) CCD image type to correct
(max_cac= 0) Maximum image caching memory (in Mbytes)
(noproc = no) List processing steps only?
(fixpix = no) Fix bad CCD lines and columns?
(oversca= yes) Apply overscan strip correction?
(trim = yes) Trim the image?
(zerocor= yes) Apply zero level correction?
(darkcor= no) Apply dark count correction?
(flatcor= yes) Apply flat field correction?
(illumco= no) Apply illumination correction?
(fringec= no) Apply fringe correction?
(readcor= no) Convert zero level image to readout correction?
(scancor= no) Convert flat field image to scan correction?
(readaxi= line) Read out axis (column|line)
(fixfile= ) File describing the bad lines and columns
(biassec= image) Overscan strip image section
(trimsec= image) Trim data section
(zero = ) Zero level calibration image
(dark = ) Dark count calibration image
(flat = ) Flat field images
(illum = ) Illumination correction images
(fringe = ) Fringe correction images
(minrepl= 1.) Minimum flat field value
(scantyp= shortscan) Scan type (shortscan|longscan)
(nscan = 1) Number of short scan lines
(interac= yes) Fit overscan interactively?
(functio= chebyshev) Fitting function
(order = 1) Number of polynomial terms or spline pieces
(sample = *) Sample points to fit
(naverag= 1) Number of sample points to combine
(niterat= 1) Number of rejection iterations
(low_rej= 3.) Low sigma rejection factor
(high_re= 3.) High sigma rejection factor
(grow = 0.) Rejection growing radius
(mode = al)
PACKAGE = ccdred
TASK = flatcombine
input = flat615.fits, flat226.fits List of flat field images to combine
(output = flattest) Output flat field root name
(combine= average) Type of combine operation
(reject = crreject) Type of rejection
(ccdtype= flat) CCD image type to combine
(process= yes) Process images before combining?
(subsets= yes) Combine images by subset parameter?
(delete = no) Delete input images after combining?
(clobber= no) Clobber existing output image?
(scale = mode) Image scaling
(statsec= ) Image section for computing statistics
(nlow = 1) minmax: Number of low pixels to reject
(nhigh = 1) minmax: Number of high pixels to reject
(nkeep = 1) Minimum to keep (pos) or maximum to reject (neg)
(mclip = yes) Use median in sigma clipping algorithms?
(lsigma = 3.) Lower sigma clipping factor
(hsigma = 3.) Upper sigma clipping factor
(rdnoise= rdnoise) ccdclip: CCD readout noise (electrons)
(gain = gain) ccdclip: CCD gain (electrons/DN)
(snoise = 0.) ccdclip: Sensitivity noise (fraction)
(pclip = 0.) pclip: Percentile clipping parameter
(blank = 1.) Value if there are no pixels
(mode = al)
Thank you in advance
Camilla
I have a problem with the Iraf FLATCOMBINE task, when trying to create a combined flat file of UVES flat files (ff lamp 1) I get errors. I have tried running it with different lists and by typing in the single files. Everytime I get an error:
List of flat field images to combine (flat615.fits, flat226.fits):
ERROR on line 47: No images to combine
flatcombine ()
But I am able to combine the flats using COMBINE, and ZEROCOMBINE also works..
To solve the problem I have so far tried the following:
-removing iraf directory and creating a new using mkiraf (commented a few
suggested lines in the login.cl file in)
-flatcombine *.fits ccdtype=""
-gave all rights to the fits files
-created a new directory, where I moved the flat fits files to, renamed them and gave all rights to them
-changed setinstrument (echelle,direct)
-changed ccdproc's biassec and trimsec ([values],image)
-unlearn flatcombine and ccdproc.
My task input look like this:
PACKAGE = ccdred
TASK = setinstrument
instrume= echelle Instrument ID (type ? for a list)
(site = kpno) Site ID
(directo= ccddb$) Instrument directory
(review = yes) Review instrument parameters?
query = echelle Instrument ID (type q to quit)
(mode = al)
PACKAGE = ccdred
TASK = ccdproc
images = flat615.fits, flat226.fits List of CCD images to correct
(output = flattest2.fits) List of output CCD images
(ccdtype= ) CCD image type to correct
(max_cac= 0) Maximum image caching memory (in Mbytes)
(noproc = no) List processing steps only?
(fixpix = no) Fix bad CCD lines and columns?
(oversca= yes) Apply overscan strip correction?
(trim = yes) Trim the image?
(zerocor= yes) Apply zero level correction?
(darkcor= no) Apply dark count correction?
(flatcor= yes) Apply flat field correction?
(illumco= no) Apply illumination correction?
(fringec= no) Apply fringe correction?
(readcor= no) Convert zero level image to readout correction?
(scancor= no) Convert flat field image to scan correction?
(readaxi= line) Read out axis (column|line)
(fixfile= ) File describing the bad lines and columns
(biassec= image) Overscan strip image section
(trimsec= image) Trim data section
(zero = ) Zero level calibration image
(dark = ) Dark count calibration image
(flat = ) Flat field images
(illum = ) Illumination correction images
(fringe = ) Fringe correction images
(minrepl= 1.) Minimum flat field value
(scantyp= shortscan) Scan type (shortscan|longscan)
(nscan = 1) Number of short scan lines
(interac= yes) Fit overscan interactively?
(functio= chebyshev) Fitting function
(order = 1) Number of polynomial terms or spline pieces
(sample = *) Sample points to fit
(naverag= 1) Number of sample points to combine
(niterat= 1) Number of rejection iterations
(low_rej= 3.) Low sigma rejection factor
(high_re= 3.) High sigma rejection factor
(grow = 0.) Rejection growing radius
(mode = al)
PACKAGE = ccdred
TASK = flatcombine
input = flat615.fits, flat226.fits List of flat field images to combine
(output = flattest) Output flat field root name
(combine= average) Type of combine operation
(reject = crreject) Type of rejection
(ccdtype= flat) CCD image type to combine
(process= yes) Process images before combining?
(subsets= yes) Combine images by subset parameter?
(delete = no) Delete input images after combining?
(clobber= no) Clobber existing output image?
(scale = mode) Image scaling
(statsec= ) Image section for computing statistics
(nlow = 1) minmax: Number of low pixels to reject
(nhigh = 1) minmax: Number of high pixels to reject
(nkeep = 1) Minimum to keep (pos) or maximum to reject (neg)
(mclip = yes) Use median in sigma clipping algorithms?
(lsigma = 3.) Lower sigma clipping factor
(hsigma = 3.) Upper sigma clipping factor
(rdnoise= rdnoise) ccdclip: CCD readout noise (electrons)
(gain = gain) ccdclip: CCD gain (electrons/DN)
(snoise = 0.) ccdclip: Sensitivity noise (fraction)
(pclip = 0.) pclip: Percentile clipping parameter
(blank = 1.) Value if there are no pixels
(mode = al)
Thank you in advance
Camilla
hanin wrote on Mar 28, 2008
AFAIK, you must have the keyword 'IMAGETYP= flat' in the flat images' header. This is because you have 'ccdtype= flat' in the flatcombine parameter.
cmiiw..
cmiiw..
camjh wrote on Mar 28, 2008
Thank you!
Changing the IMAGETYP in the header to flat did the trick.
Cheers
Camilla
Changing the IMAGETYP in the header to flat did the trick.
Cheers
Camilla
Mike Fitzpatrick wrote on Mar 28, 2008
To solve the problem I have so far tried the following:
......
-flatcombine *.fits ccdtype=""
For the record, Hanin has the correct reply, but I was confused by this (failed) solution since it is normally what we also recommend. Setting ccdtype="" tells the task to ignore the IMAGETYP and use the images regardless, so this should have worked. Are you sure this didn't work?
-Mike
Last post on Mar 28, 2008