IRAF/DAOPHOT Photometry Simple Steps - BATC home page

IRAF/DAOPHOT
♦ IRAF starten (z.B. "xiraf", "xiraf -ds9" oder "cl")
♦ Benoetigte packages laden:
♦ Ubersicht ueber die main tasks und die dazugehoerigen
Parameterfiles:

main tasks:
PACKAGE = daophot
TASK = daofind
I R A F
Image Reduction and Analysis Facility
image = bild Input image(s)
output = default Output coordinate file(s) (default: image.coo.)
(starmap=
) Output density enhancement image(s)
(skymap =
) Output sky image(s)
(datapar=
) Data dependent parameters
(findpar=
) Object detection parameters
(boundar= nearest) Boundary extension (constant|nearest|reflect|wra
(constan= 0.) Constant for boundary extension
(interac= no) Interactive mode
(verify = no) Verify critical daofind parameters
(update = no) Update critical daofind parameters
(verbose= yes) Print daofind messages
(graphic= no) Graphics device
(display= yes) Display device
(icomman=
) Image cursor: [x y wcs] key [cmd]
(gcomman=
) Graphics cursor: [x y wcs] key [cmd]
(mode = ql)
PACKAGE = daophot
TASK = phot
I R A F
Image Reduction and Analysis Facility
image = bild Input image(s)
coords = default Input coordinate list(s) (default: image.coo.)
output = default Output photometry file(s) (default: image.mag.)
skyfile =
Input sky value file(s)
(plotfil=
) Output plot metacode file
(datapar=
) Data dependent parameters
(centerp=
) Centering parameters
(fitskyp=
) Sky fitting parameters
(photpar=
) Photometry parameters
(interac= no) Interactive mode
(radplot= no) Plot the radial profiles
(verify = no) Verify critical phot parameters
(update = no) Update critical phot parameters
(verbose= yes) Print phot messages
(graphic= stdgraph) Graphics device
(display= stdimage) Display device
(icomman=
) Image cursor: [x y wcs] key [cmd]
(gcomman=
) Graphics cursor: [x y wcs] key [cmd]
(mode = ql)
PACKAGE = daophot
TASK = pstselect
I R A F
Image Reduction and Analysis Facility
image = bild Image for which to build psf star list
photfile= default Photometry file (default: image.mag.)
pstfile = default Output psf star list file (default: image.pst.)
maxnpsf = 1000 Maximum number of psf stars
(plotfil=
) Output plot metacode file
(datapar= kdat.par) Data dependent parameters
(daopars= dao.par) Psf fitting parameters
(interac= yes) Select psf stars interactively
(plottyp= mesh) Default plot type (mesh|contour|radial)
(verbose= )_.verbose) Print pstselect messages
(verify = )_.verify) Verify critical pstselect parameters
(update = )_.update) Update critical pstselect parameters
(graphic= )_.graphics) Graphics device
(display= )_.display) Image display device
(icomman=
) Image cursor: [x y wcs] key [cmd]
(gcomman=
) Graphics cursor: [x y wcs] key [cmd]
(mode = ql)

PACKAGE = daophot
TASK = psf
I R A F
Image Reduction and Analysis Facility
image = bild Input image(s) for which to build PSF
photfile= default Input photometry file(s) (default: image.mag.)
pstfile = default Input psf star list(s) (default: image.pst)
psfimage= default Output PSF image(s) (default: image.psf.)
opstfile= default Output PSF star list(s) (default:image.pst.)
groupfil= default Output PSF star group file(s) (default: image.ps
(plotfil=
) Output plot metacode file
(datapar=
) Data dependent parameters
(daopars=
) Psf fitting parameters
(matchby= yes) Match psf star list to photometry file(s) by id
(interac= no) Compute the psf interactively
(showplo= no) Show plots of PSF stars
(plottyp= mesh) Default plot type (mesh|contour|radial)
(verbose= )_.verbose) Print psf messages
(verify = )_.verify) Verify critical psf parameters
(update = )_.update) Update critical psf parameters
(graphic= )_.graphics) Graphics device
(display= )_.display) Display device
(icomman=
) Image cursor: [x y wcs] key [cmd]
(gcomman=
) Graphics cursor: [x y wcs] key [cmd]
(mode = ql)
PACKAGE = daophot
TASK = allstar
image =
photfile=
psfimage=
allstarf=
rejfile =
subimage=
(cache =
(datapar=
(daopars=
(verbose=
(verify =
(update =
(mode =
I R A F
Image Reduction and Analysis Facility
bild Image corresponding to photometry
default Input photometry file (default: image.mag.)
default PSF image (default: image.psf.)
default Output photometry file (default: image.als.)
default Output rejections file (default: image.arj.)
default Subtracted image (default: image.sub.)
yes) Cache the data in memory
) Data dependent parameters
) Psf fitting parameters
yes) Print allstar messages
no) Verify critical allstar parameters
)_.update) Update critical allstar parameters
ql)

parameter files:
PACKAGE = daophot
TASK = datapars
I R A F
Image Reduction and Analysis Facility
(scale = 1.) Image scale in units per pixel
(fwhmpsf= 3.) FWHM of the PSF in scale units
(emissio= yes) Features are positive
(sigma = 0.5) Standard deviation of background in counts
(datamin= 0.) Minimum good data value
(datamax= 65000.) Maximum good data value
(noise = poisson) Noise model
(ccdread=
) CCD readout noise image header keyword
(gain =
) CCD gain image header keyword
(readnoi= 4.4) CCD readout noise in electrons
(epadu = 1.31) Gain in electrons per count
(exposur= INDEF) Exposure time image header keyword
(airmass= INDEF) Airmass image header keyword
(filter = INDEF) Filter image header keyword
(obstime= INDEF) Time of observation image header keyword
(itime = INDEF) Exposure time
(xairmas= INDEF) Airmass
(ifilter= INDEF) Filter
(otime = INDEF) Time of observation
(mode = ql)
PACKAGE = daophot
TASK = findpars
I R A F
Image Reduction and Analysis Facility
(thresho= 3.) Threshold in sigma for feature detection
(nsigma = 1.5) Width of convolution kernel in sigma
(ratio = 1.) Ratio of minor to major axis of Gaussian kernel
(theta = 0.) Position angle of major axis of Gaussian kernel
(sharplo= 0.2) Lower bound on sharpness for feature detection
(sharphi= 1.) Upper bound on sharpness for feature detection
(roundlo= -1.) Lower bound on roundness for feature detection
(roundhi= 1.) Upper bound on roundness for feature detection
(mkdetec= no) Mark detections on the image display
(mode = ql)
I R A F
Image Reduction and Analysis Facility
PACKAGE = daophot
TASK = fitskypars
(salgori=
(annulus=
(dannulu=
(skyvalu=
(smaxite=
(sloclip=
(shiclip=
(snrejec=
(sloreje=
(shireje=
(khist =
(binsize=
(smooth =
(rgrow =
(mksky =
(mode =
median) Sky fitting algorithm
10.) Inner radius of sky annulus in scale units
10.) Width of sky annulus in scale units
0.) User sky value
10) Maximum number of sky fitting iterations
0.) Lower clipping factor in percent
0.) Upper clipping factor in percent
50) Maximum number of sky fitting rejection iteratio
3.) Lower K-sigma rejection limit in sky sigma
3.) Upper K-sigma rejection limit in sky sigma
3.) Half width of histogram in sky sigma
0.1) Binsize of histogram in sky sigma
no) Boxcar smooth the histogram
0.) Region growing radius in scale units
no) Mark sky annuli on the display
ql)

I R A F
Image Reduction and Analysis Facility
PACKAGE = daophot
TASK = centerpars
(calgori=
(cbox =
(cthresh=
(minsnra=
(cmaxite=
(maxshif=
(clean =
(rclean =
(rclip =
(kclean =
(mkcente=
(mode =
centroid) Centering algorithm
3.) Centering box width in scale units
0.) Centering threshold in sigma above background
1.) Minimum signal-to-noise ratio for centering algo
10) Maximum iterations for centering algorithm
1.) Maximum center shift in scale units
no) Symmetry clean before centering
1.) Cleaning radius in scale units
2.) Clipping radius in scale units
3.) K-sigma rejection criterion in skysigma
no) Mark the computed center
ql)
PACKAGE = daophot
TASK = photpars
(weighti=
(apertur=
(zmag =
(mkapert=
(mode =
I R A F
Image Reduction and Analysis Facility
constant) Photometric weighting scheme
3) List of aperture radii in scale units
25.) Zero point of magnitude scale
no) Draw apertures on the display
ql)
PACKAGE = daophot
TASK = daopars
I R A F
Image Reduction and Analysis Facility
(functio= gauss) Form of analytic component of psf model
(varorde= 0) Order of empirical component of psf model
(nclean = 0) Number of cleaning iterations for computing psf
(saturat= no) Use wings of saturated stars in psf model comput
(matchra= 3.) Object matching radius in scale units
(psfrad = 11.) Radius of psf model in scale units
(fitrad = 3.) Fitting radius in scale units
(recente= yes) Recenter stars during fit
(fitsky = no) Recompute group sky value during fit
(groupsk= no) Use group rather than individual sky values
(sannulu= 20.) Inner radius of sky fitting annulus in scale uni
(wsannul= 10.) Width of sky fitting annulus in scale units
(flaterr= 0.75) Flat field error in percent
(proferr= 5.) Profile error in percent
(maxiter= 50) Maximum number of fitting iterations
(clipexp= 6) Bad data clipping exponent
(clipran= 2.5) Bad data clipping range in sigma
(mergera= INDEF) Critical object merging radius in scale units
(critsnr= 1.) Critical S/N ratio for group membership
(maxnsta= 10000) Maximum number of stars to fit
(maxgrou= 60) Maximum number of stars to fit per group
(mode = ql)

tvmark:
PACKAGE = tv
TASK = tvmark
frame =
coords =
(logfile=
(autolog=
(outimag=
(deletio=
(command=
(mark =
(radii =
(lengths=
(font =
(color =
(label =
(number =
(nxoffse=
(nyoffse=
(pointsi=
(txsize =
(toleran=
(interac=
(mode =
I R A F
Image Reduction and Analysis Facility
1 Default frame number for display
bild.coo.1 Input coordinate list
) Output log file
no) Automatically log each marking command
) Output snapped image
) Output coordinate deletions list
) Image cursor: [x y wcs] key [cmd]
circle) The mark type
10) Radii in image pixels of concentric circles
10) Lengths and width in image pixels of concentric
raster) Default font
202) Gray level of marks to be drawn
no) Label the marked coordinates
no) Number the marked coordinates
10) X offset in display pixels of number
10) Y offset in display pixels of number
1) Size of mark type point in display pixels
1) Size of text and numbers in font units
1.5) Tolerance for deleting coordinates in image pixel
no) Mode of use
ql)
tvmark eignet sich zum Ueberpruefen
und Ergaenzen der Ergebnisse von
daofind.
Dabei sollte man alleridngs die
"cuts" variieren. Zum einen kann
man helle Sterne am besten bei
moderaten "cuts" verifizieren
(links oben) und zum anderen kann
man schwache Sterne nur bei sehr
hohen Kontrasten (rechts unten)
erkennen.
daofind arbeitet nicht perfekt. Man
wird nie genau alle Sterne finden.
Um Fehldetektionen aus der
Koordinaten-liste zu loeschen und
nicht detektierte Objekte
hinzuzufuegen benutzt man tvmark
wie folgt:
tvmark 1 bild.coo.1 interactive=yes
oder
tvmark 1 bild.coo.1 int+
Dieser Befehl bringt einen in
einen interaktiven Mode (aehnlich
wie bei imexam). Es gibt folgende
nuetzliche Keystroke-commands:
l - plotten der Koordinaten-liste
d - loeschen (delete)
a - hinzufuegen (add)
Fuer "a" und "d" muss sich der
cursor auf der entsprechenden
Koordinate befinden.