Hello all,I just recently found that I can download data from the ESO WFI telescope, but on some of the data, after processing, there are still "streaks" (in the final stacked image).Please forgive me if I use the wrong terms, I don't have any formal training (everything I know about IRAF I have taught myself).I would like to know if anyone has any ideas on what I am doing wrong (I am just looking to make "pretty pictures" from the data, not scientific stuff).Here is the final image (with settings showing the "lines"): http://wierdling.longboys.net/iraf/irafImage.jpgI am running IRAF on a CENTOS os. Here are the steps that I do (it is for some data for M 83 by the way):
PHP Formatted Code
esowfi
> esohdr W
*.fits
esowfi
> esohdr W
*.fits
esowfi
> epar ccdproc
-- see
#1
esowfi
> epar zerocombine
-- see
#2
esowfi
> zerocombine W
*.fits
esowfi
> epar flatcombine
-- see
#3
esowfi
> flatcombine W
*.fits
-- take output and make a
list called flatfiles
esowfi
> epar ccdproc
-- see
#4
esowfi
> ccdproc W
*.fits
esowfi
> msczero
-- update zero reference
for each image in the
list
esowfi
> msczero
-- create
list of coords
esowfi
> epar msccmatch
-- see
#5
esowfi
> msccmatch
esowfi
> epar mscimage
-- see
#6
esowfi
> mscimage
esowfi
> epar mscimatch
-- see
#7
esowfi
> mscimatch
-- see
#8 for output
esowfi
> mscstack
#1 -- epar for ccdproc
I R A F
Image Reduction and Analysis Facility
PACKAGE
= mscred
TASK
= ccdproc
images
= W
*fits
List of Mosaic CCD images to process
(output
= ) List of output processed images
(bpmasks
= ) List of output bad pixel masks
(ccdtype
= object
) CCD image type to process
(noproc
= no
) List processing steps only?
(xtalkco
= no
) Apply crosstalk correction?
(fixpix
= yes
) Apply bad pixel mask correction?
(oversca
= yes
) Apply overscan strip correction?
(trim = yes
) Trim the image?
(zerocor
= no
) Apply zero level correction?
(darkcor
= no
) Apply dark
count correction?
(flatcor
= no
) Apply flat field correction?
(sflatco
= no
) Apply sky flat field correction?
(split = no
) Use split images during processing?
(merge
= no
) Merge amplifiers from same CCD?
(xtalkfi
= ) Crosstalk
file
(fixfile
= @badPics
) List of bad pixel masks
(saturat
= INDEF
) Saturated pixel threshold
(sgrow
= 0) Saturated pixel grow radius
(bleed
= INDEF
) Bleed pixel threshold
(btrail
= 20) Bleed trail minimum length
(bgrow
= 0) Bleed pixel grow radius
(biassec
= !biassec
) Overscan strip image section
(trimsec
= !trimsec
) Trim data section
(zero
= ) List of zero level calibration images
(dark
= ) List of dark
count calibration images
(flat
= ) List of flat field images
(sflat
= ) List of secondary flat field images
(minrepl
= 1.) Minimum flat field value
(interac
= no
) Fit overscan interactively?
(functio
= legendre
) Fitting
function
(order
= 4) 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
= ql
)
#2 -- epar for zerocombine
I R A F
Image Reduction and Analysis Facility
PACKAGE
= mscred
TASK
= zerocombine
input
= W
*.fits
List of zero level images to combine
(output
= Zero
) Output zero level name
(combine
= average
) Type of combine operation
(reject
= avsigclip
) Type of rejection
(ccdtype
= zero
) CCD image type to combine
(process
= yes
) Process images before combining?
(delete
= no
) Delete input images after combining?
(scale
= none
) Image scaling
(statsec
= ) Image section
for computing statistics
(nlow
= 0) 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
= 0.) ccdclip
: CCD readout noise
(electrons
)
(gain
= 1.) ccdclip
: CCD gain
(electrons
/DN
)
(snoise
= 0.) ccdclip
: Sensitivity noise
(fraction
)
(pclip
= -0.5) pclip
: Percentile clipping parameter
(blank
= 0.) Value
if there are no pixels
(mode
= ql
)
#3 -- epar for flatcombine
I R A F
Image Reduction and Analysis Facility
PACKAGE
= mscred
TASK
= flatcombine
input
= W
*.fits
List of flat field images to combine
(output
= F
) Output flat field root name
(combine
= average
) Type of combine operation
(reject
= avsigclip
) 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?
(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
= 0.) ccdclip
: CCD readout noise
(electrons
)
(gain
= 1.) ccdclip
: CCD gain
(electrons
/DN
)
(snoise
= 0.) ccdclip
: Sensitivity noise
(fraction
)
(pclip
= -0.5) pclip
: Percentile clipping parameter
(blank
= 1.) Value
if there are no pixels
(mode
= ql
)
#4 -- epar for ccdproc
I R A F
Image Reduction and Analysis Facility
PACKAGE
= mscred
TASK
= ccdproc
images
= W
*fits
List of Mosaic CCD images to process
(output
= ) List of output processed images
(bpmasks
= ) List of output bad pixel masks
(ccdtype
= object
) CCD image type to process
(noproc
= no
) List processing steps only?
(xtalkco
= no
) Apply crosstalk correction?
(fixpix
= yes
) Apply bad pixel mask correction?
(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?
(sflatco
= no
) Apply sky flat field correction?
(split = no
) Use split images during processing?
(merge
= no
) Merge amplifiers from same CCD?
(xtalkfi
= ) Crosstalk
file
(fixfile
= @badPics
) List of bad pixel masks
(saturat
= INDEF
) Saturated pixel threshold
(sgrow
= 0) Saturated pixel grow radius
(bleed
= INDEF
) Bleed pixel threshold
(btrail
= 20) Bleed trail minimum length
(bgrow
= 0) Bleed pixel grow radius
(biassec
= !biassec
) Overscan strip image section
(trimsec
= !trimsec
) Trim data section
(zero
= Zero
) List of zero level calibration images
(dark
= ) List of dark
count calibration images
(flat
= @flatfiles
) List of flat field images
(sflat
= ) List of secondary flat field images
(minrepl
= 1.) Minimum flat field value
(interac
= no
) Fit overscan interactively?
(functio
= legendre
) Fitting
function
(order
= 4) 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
#5 -- epar for msccmatch
I R A F
Image Reduction and Analysis Facility
PACKAGE
= mscred
TASK
= msccmatch
input
= @rcIn
List of input mosaic exposures
coords
= rccoords Coordinate
file (ra
/dec
)
(outcoor
= ) List of updated coordinate files
(usebpm
= yes
) Use bad pixel masks?
(verbose
= no
) Verbose?
# Coarse Search
(nsearch
= 10) Maximum number of positions to
use in search
(search
= 3.) Translation search radius
(arcsec
)
(rsearch
= 5.) Rotation search radius
(deg
)
# Fine Centroiding
(cbox
= 11) Centering box
(pixels
)
(maxshif
= 5.) Maximum centering shift to accept
(arcsec
)
(csig
= 0.1) Maximum centering uncertainty to accept
(arcsec
)
(cfrac
= 0.5) Minimum fraction of accepted centers
(listcoo
= yes
) List centered coordinates in verbose mode?
# WCS Fitting
(nfit
= 4) Min for fit
(>0) or
max not found
(<=0)
(rms
= 2.) Maximum fit RMS to accept
(arcsec
)
(fitgeom
= general
) Fitting geometry
(reject
= 3.) Fitting rejection limit
(sigma
)
(update
= yes
) Update coordinate systems?
(interac
= no
) Interactive?
(fit
= yes
) Interactive fitting?
(graphic
= stdgraph
) Graphics device
(cursor
= ) Graphics cursor
accept
= yes Accept solution?
(mode
= ql
)
#6 -- epar for mscimage
I R A F
Image Reduction and Analysis Facility
PACKAGE
= mscred
TASK
= mscimage
input
= @rcIn
List of input mosaic exposures
output
= @rcOut
List of output images
(format
= image
) Output format
(image
|mef
)
(pixmask
= yes
) Create pixel mask?
(verbose
= )_
.verbose
) Verbose output?
# Output WCS parameters
(wcssour
= image
) Output WCS source
(image
|parameters
)
(referen
= ) Reference image
(ra
= INDEF
) RA of tangent point
(hours
)
(dec
= INDEF
) DEC of tangent point
(degrees
)
(scale
= INDEF
) Scale
(arcsec
/pixel
)
(rotatio
= INDEF
) Rotation of DEC from N to E
(degrees
)
# Resampling parmeters
(blank
= 0.) Blank value
(interpo
= linear
) Interpolant
for data
(minterp
= linear
) Interpolant
for mask
(boundar
= reflect
) Boundary extension
(constan
= 0.) Constant boundary extension value
(fluxcon
= no
) Preserve flux per unit area?
(ntrim
= 8) Edge
trim in
each extension
(nxblock
= 2048) X dimension of working block size in pixels
(nyblock
= 1024) Y dimension of working block size in pixels
# Geometric mapping parameters
(interac
= no
) Fit mapping interactively?
(nx
= 10) Number of x grid points
(ny
= 20) Number of y grid points
(fitgeom
= general
) Fitting geometry
(xxorder
= 4) Order of x fit in x
(xyorder
= 4) Order of x fit in y
(xxterms
= half
) X fit cross terms type
(yxorder
= 4) Order of y fit in x
(yyorder
= 4) Order of y fit in y
(yxterms
= half
) Y fit cross terms type
(fd_in
= )
(fd_ext
= )
(fd_coor
= )
(mode
= ql
)
#7 -- epar for mscimatch
I R A F
Image Reduction and Analysis Facility
PACKAGE
= mscred
TASK
= mscimatch
input
= @rcOut
List of images
coords
= rccoords
File of coordinates
(bpm
= @rcbpm
) List of bad pixel masks
(measure
= ) Measurment
file
(scale
= yes
) Determine scale?
(zero
= yes
) Determine zero offset?
(box1
= 21) Inner box size
for statistics
(box2
= 51) Outer box size
for statistics
(lower
= 1.) Lower limit
for good data
(upper
= INDEF
) Upper limit
for good data
(niterat
= 3) Number of sigma clipping iterations
(sigma
= 3.) Sigma clipping factor
(interac
= no
) Interactive?
(verbose
= yes
) Verbose?
accept
= yes Accept scaling and update images?
(mode
= ql
)
#8 -- output for mscimatch
MSCIMATCH
:
Reading region coordinates from rccoords
48 coordinates read
Measuring regions in rc1
.fits
...
Using bad pixel mask rc1_bpm
.pl
...
48 good regions measured
Measuring regions in rc2
.fits
...
Using bad pixel mask rc2_bpm
.pl
...
40 good regions measured
Measuring regions in rc3
.fits
...
Using bad pixel mask rc3_bpm
.pl
...
47 good regions measured
Measuring regions in rc4
.fits
...
Using bad pixel mask rc4_bpm
.pl
...
45 good regions measured
Measuring regions in rc5
.fits
...
Using bad pixel mask rc5_bpm
.pl
...
40 good regions measured
Determining scale factors
...
rc1
.fits
: 1.0000 (0.0000) 0.00 (0.00)
rc2
.fits
: 1.1855 (0.0802) 336.62 (858.67)
rc3
.fits
: 1.0588 (0.0380) 108.96 (246.84)
rc4
.fits
: 0.9974 (0.0454) 171.42 (232.07)
rc5
.fits
: 1.0047 (0.0445) -192.14 (776.87)
Accept scaling and update images?
(yes
):
#9 -- epar mscstack
I R A F
Image Reduction and Analysis Facility
PACKAGE
= mscred
TASK
= mscstack
input
= @rcOut
List of images to combine
output
= m83_rc5
.fits Output image
(headers
= ) List of
header files
(optional
)
(bpmasks
= ) List of bad pixel masks
(optional
)
(rejmask
= ) List of rejection masks
(optional
)
(nrejmas
= ) List of number rejected masks
(optional
)
(expmask
= ) List of exposure masks
(optional
)
(sigmas
= ) List of sigma images
(optional
)
(combine
= median
) Type of combine operation
(median
|average
)
(reject
= avsigclip
) Type of rejection
(masktyp
= goodvalue
) Mask type
(maskval
= 0.) Mask value
(blank
= 0.) Value
if there are no pixels
(scale
= !mscscale
) Image scaling
(zero
= !msczero
) Image zero point offset
(weight
= none
) Image weights
(statsec
= ) Image section
for computing statistics
(lthresh
= 1.) Lower threshold
(hthresh
= INDEF
) Upper threshold
(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
= 0.) ccdclip
: CCD readout noise
(electrons
)
(gain
= 1.) ccdclip
: CCD gain
(electrons
/DN
)
(snoise
= 0.) ccdclip
: Sensitivity noise
(fraction
)
(sigscal
= 0.1) Tolerance
for sigma clipping scaling corrections
(pclip
= -0.5) pclip
: Percentile clipping parameter
(grow
= 0.) Radius
(pixels
) for neighbor rejection
(mode
= ql
)