IRAF ABC

IRAF Manual I
IRAF ABC
Last updated: 2006 10 13 (2004 4 1)
,
hshwang@astro.snu.ac.kr
כ
1.
IRAF (Image Reduction and Analysis Facilities) CCD
. IRAF (http://noao.iraf.edu)
כ .
2 IRAF , 3 IRAF imexam
כ. 4 (preprocessing)
כ, 5 IRAF כ .
, DAOPHOT PSF
http://astro.snu.ac.kr/ ∼ hshwang DAOPHOT , PHOT
http://astro.snu.ac.kr/ ∼ hshwang כ.
PHOTCAL http://astro.snu.ac.kr/ ∼ obscos
PHOTCAL TUTORIAL (, , 2001) ELLIPSE
. (spectroscopy) http://astro.snu.ac.kr/ ∼ hshwang
Introduction to analysis of single dispersion spectra with IRAF .
e-mail כ.
2.
2.1. IRAF
IRAF כ
(Fedora Core 5) כ.
http//iraf.noao.edu PC-IRAF pciraf.ps.gz כ כ.
1. ’iraf’ . fedora ’’ ⇒
’’ ⇒ ’ ’ /bin/csh
/iraf/iraf/local .
2. (root)
.
[hshwang@localhost ~/utils]$ su
Password:
[hshwang@localhost /]# mkdir /iraf
[hshwang@localhost /]# mkdir /iraf/iraf
[hshwang@localhost /]# mkdir /iraf/iraf/local
[hshwang@localhost /]# mkdir /iraf/irafbin
1

2 HWANG & LEE
[hshwang@localhost /]# mkdir /iraf/irafbin/bin.redhat
[hshwang@localhost /]# mkdir /iraf/irafbin/noao.bin.redhat
[hshwang@localhost /]# mkdir /iraf/x11iraf
[hshwang@localhost /]# mkdir /iraf/extern
[hshwang@localhost /]# mkdir /iraf/extern/stsdas
[hshwang@localhost /]# mkdir /iraf/extern/tables
iraf .
[hshwang@localhost /]# chown -R iraf.iraf /iraf
3. http//iraf.noao.edu PC-IRAF FTP Dir
as.pcix.gen.gz, ib.rhux.x86.gz, nb.rhux.x86.gz (
/home/hshwang/utils/ ) .
.
4. ’iraf’
.
[hshwang@localhost ~]$ su iraf
[iraf@localhost ~]$ whoami
iraf
[iraf@localhost ~]$ cd /iraf/iraf
[iraf@localhost iraf]$ cat /home/hshwang/utils/as.pcix.gen.gz | zcat | tar -xpf -
[iraf@localhost iraf]$ cd /iraf/iraf/bin.redhat/
[iraf@localhost bin.redhat]$ cat /home/hshwang/utils/ib.rhux.x86.gz | zcat | tar -xpf -
[iraf@localhost bin.redhat]$ cd /iraf/iraf/noao/bin.redhat/
[iraf@localhost bin.redhat]$ cat /home/hshwang/utils/nb.rhux.x86.gz | zcat | tar -xpf -
5. script . root .
[iraf@localhost bin.redhat]$ su
Password:
[iraf@localhost bin.redhat]# cd /iraf/iraf/unix/hlib/
[iraf@localhost hlib]# setenv iraf /iraf/iraf
[iraf@localhost hlib]# source irafuser.csh
[iraf@localhost hlib]# ./install
install enter , install כ.
2.2. xgterm
xgterm http//iraf.noao.edu X11IRAF FTP Dir
x11iraf-v1.3.1-bin.redhat.tar.gz ( /home/hshwang/utils/
) . root .

IRAF ABC 3
[hshwang@localhost x11iraf]$ whoami
hshwang
[hshwang@localhost ~/utils]$ pwd
/home/hshwang/utils
[hshwang@localhost ~/utils]$ mkdir x11iraf
[hshwang@localhost ~/utils]$ cd x11iraf
[hshwang@localhost x11iraf]$ tar xzvf ../x11iraf-v1.3.1-bin.redhat.tar.gz
[hshwang@localhost x11iraf]$ su
Password:
[hshwang@localhost x11iraf]# ./install
enter OK. xgterm http//iraf.noao.edu Fedora/Mandrake
Patched XGTerm Available Download from ftp://iraf.noao.edu/pub/xgterm.fedora
xgterm.fedora ( /home/hshwang/utils/ )
. root .
[hshwang@localhost utils]# mv xgterm.fedora /usr/local/bin/xgterm
mv: overwrite ‘/usr/local/bin/xgterm’ yes
[hshwang@localhost bin]# cd /usr/local/bin
[hshwang@localhost bin]# chown root.root xgterm
[hshwang@localhost bin]# chmod 755 xgterm
[hshwang@localhost ~]$ /usr/local/bin/xgterm -sb&
xgterm כ .
2.3. ds9
ds9 display http://hea-www.harvard.edu/RD/ds9/ Linux (
) ds9.linux.4.0b10.tar.gz . root
.
[hshwang@localhost utils]# whoami
root
[hshwang@localhost utils]# cd /usr/local/bin
[hshwang@localhost bin]# tar xzvf /home/hshwang/utils/ds9.linux.4.0b9.tar.gz
[hshwang@localhost bin]# ./ds9 &
2.4. external packages
IRAF כ external package
כ STSDAS TABLES כ .
http://www.stsci.edu/resources/software hardware/stsdas/download STSDAS SOURCE
(stsdas35.tar.gz) Binaries Redhat (stsdas35.bin.rh.tar.gz) (
/home/hshwang/utils/ ) . TABLES SOURCE
(tables35.tar.gz) Binaries Redhat (tables35.bin.rh.tar.gz) . iraf
.

4 HWANG & LEE
1. Binary
[hshwang@localhost ~]$ su iraf
Password:
iraf@localhost% mkdir /iraf/extern/stsdas/bin.redhat
iraf@localhost% cd /iraf/extern/stsdas/bin.redhat
iraf@localhost% tar xzvf /home/hshwang/utils/stsdas35.bin.rh.tar.gz
iraf@localhost% mkdir /iraf/extern/tables/bin.redhat
iraf@localhost% cd /iraf/extern/tables/bin.redhat/
iraf@localhost% tar xzvf /home/hshwang/utils/tables35.bin.rh.tar.gz
2. extern.pkg
iraf@localhost% vi /iraf/iraf/unix/hlib/extern.pkg
# External (non core-system) packages. To install a new package, add the
# two statements to define the package root directory and package task,
# then add the package helpdb to the ‘helpdb’ list.
reset noao = iraf$noao/
task noao.pkg = noao$noao.cl
reset tables = /iraf/extern/tables/ ##
task tables.pkg = tables$tables.cl ##
reset stsdas = /iraf/extern/stsdas/ ##
task stsdas.pkg = stsdas$stsdas.cl ##
reset helpdb = "lib$helpdb.mip\
,noao$lib/helpdb.mip\
,tables$lib/helpdb.mip\ ##
,stsdas$lib/helpdb.mip\ ##
"
keep
3. Source
iraf@localhost% cd /iraf/extern/stsdas/
iraf@localhost% tar xzvf /home/hshwang/utils/stsdas35.tar.gz
iraf@localhost% cd /iraf/extern/tables/
iraf@localhost% tar xzvf /home/hshwang/utils/tables35.tar.gz
. 2 IRAF .
3.
3.1.
IRAF imexam
IRAF
. ’iraf’ ’root’ .

IRAF ABC 5
[hshwang@localhost ~]$ whoami
hshwang
[hshwang@localhost ~]$ pwd
/home/hshwang
[hshwang@localhost ~]$ mkdir iraf
[hshwang@localhost ~]$ cd iraf
[hshwang@localhost ~/iraf]$ mkiraf
-- creating a new uparm directory
Terminal types: xgterm,xterm,gterm,vt640,vt100,etc.
Enter terminal type: xgterm
A new LOGIN.CL file has been created in the current directory.
You may wish to review and edit this file to change the defaults.
, iraf xgterm כ,
login.cl .
[hshwang@localhost ~/iraf]$ vi login.cl
set imtype = "fits" -> set imtype = "fits" (# )
3.2.
login.cl xgterm
[hshwang@localhost ~/iraf]$ xgterm -sb&
xgterm
[hshwang@localhost ~/iraf]$ ds9 &
[hshwang@localhost ~/iraf]$ cl ; IRAF
cf) cl> logout ; IRAF
1: IRAF
1 stsdas tables ( ) כ.
display 2 IRAF כ.
cl> display dev$pix 1 fi+
[] display (image name) (frame number: 1-4) fi+
2: DS9 M51

6 HWANG & LEE
3.3. imexam
ds9 ().
imexam tool .
cl> imexam
ds9 M51(dev$pix) כ.
.
r :
s : surface plot
e : contour plot
a : ,
h :
m :
z :
q :
4.
Preprocessing
Bias, Dark, Flat, Object . Object imexam
.
/home/hshwang/IRAFman/work כ .
cl> noao.imred
im> ccdred
bias כ . imstat []
image section כ .
cc> pwd
/home/hshwang/iraf
cc> cd /home/hshwang/IRAFman/work
cc> imstat BiAS*[101:900,101:900]
-> .
כ כ image .
cc> imstat BiAS*[101:900,101:900]
# IMAGE NPIX MEAN STDDEV MIN MAX
BiAS01.fits[101:900,101:900] 640000 519.6 1027.33 0. 65535.
BiAS010.fits[101:900,101:900] 640000 519.9 27.18 377. 1489.
BiAS02.fits[101:900,101:900] 640000 519.5 27.26 372. 986.
BiAS03.fits[101:900,101:900] 640000 519.5 27.31 365. 942.
BiAS04.fits[101:900,101:900] 640000 519.5 27.29 386. 2275.

IRAF ABC 7
*****************NOTE**
CCD readout noise gain כ image header
.
cc> imhead BiAS23 lo+ | page
***********************
bias bias list .
[hshwang@localhost work]$ ls BiAS*fits > bias.list
bias image master bias image .
cc> zerocomb @bias.list out=Zero.fits combine=average reject=minmax \
ccdtype="" scale=none rdnoise=7.6 gain=5.04
(or ccdtype="" scale=none rdnoise=RDNOISE gain=GAIN)
IRAF task parameter כ
epar (taskname) parameter , IRAF
epar parameter כ כ .
, epar
( כ).
PACKAGE = ccdred
TASK = zerocombine
input = @bias.list List of zero level images to combine
(output =
Zero) Output zero level name
(combine=
average) Type of combine operation
(reject =
minmax) Type of rejection
(ccdtype=
zero) CCD image type to combine
(process=
no) Process images before combining
(delete =
no) Delete input images after combining
(clobber=
no) Clobber existing output image
(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=
7.6) ccdclip: CCD readout noise (electrons)
(gain = 5.04) ccdclip: CCD gain (electrons/DN)
(snoise =
0.) ccdclip: Sensitivity noise (fraction)

8 HWANG & LEE
(pclip = -0.5) pclip: Percentile clipping parameter
(blank = 0.) Value if there are no pixels
(mode = ql)
:go כ .
(cf. Ctrl-C: ‏,‏כ Ctrl-D: )
*****************NOTE** : IRAF task help
cc> help zerocombine
[] help (task name)
help !!!
**********************
4.1. Dark
cc> imstat DARK*[101:900,101:900]
[hshwang@localhost work]$ ls DARK* > dark.list
bias imstat , list . dark image combine bias
combine Zero.fits ccdproc task .
cc> epar ccdproc
images = @dark.list List of CCD images to correct
(fixpix =
no) Fix bad CCD lines and columns
(oversca=
no) Apply overscan strip correction
(trim = no) Trim the image
(zerocor=
yes) Apply zero level correction
(darkcor=
no) Apply dark count correction
(flatcor=
no) Apply flat field correction
(illumco=
no) Apply illumination correction
(zero = Zero) Zero level calibration image
:Ctrl-D () .
darkcomb task .
cc> darkcomb @dark.list out=Dark.fits combine=average reject=minmax \
ccdtype="" process=yes scale=exposure rdnoise=7.6 gain=5.04
4.2. Flat
cc> imstat f*[101:900,101:900]
[hshwang@localhost work]$ ls f*fits > flat.list

IRAF ABC 9
flat image filter combine
( U,B,V,I) .
cc> ccdhedit fU* subset U
cc> ccdhedit fB* subset B
cc> ccdhedit fV* subset V
cc> ccdhedit fI* subset I
*****************NOTE** : image header
cc> ccdhedit M67i_1 imagetype object
[] ccdhedit (file name)(keyword) (name)
***********************
darkcomb task ccdproc task flatcomb
task .
cc> epar ccdproc
images = @flat.list List of CCD images to correct
(fixpix =
no) Fix bad CCD lines and columns
(oversca=
no) Apply overscan strip correction
(trim = no) Trim the image
(zerocor=
yes) Apply zero level correction
(darkcor=
yes) Apply dark count correction
(flatcor=
no) Apply flat field correction
(illumco=
no) Apply illumination correction
(zero = Zero) Zero level calibration image
(dark = Dark) Dark count calibration image
:Ctrl-D () .
cc> flatcomb @flat.list out=Flat combine=median \
reject=avsigclip ccdtype="" process+ subsets+ scale=mean \
statsec=[101:900,101:900] rdnoise=7.6 gain=5.04
4.3. Object
[hshwang@localhost work]$ ls *67*fits > obj.list
cc> ccdhedit U67* subset U
cc> ccdhedit B67* subset B
cc> ccdhedit V67* subset V
cc> ccdhedit i67* subset I
object image Zero, Dark, Flat* image object
image correction . zerocor+ zerocor=yes ccdproc
task .

10 HWANG & LEE
cc> ccdproc @obj.list out=p//@obj.list ccdtype="" \
fixpix- oversc- trim- zerocor+ darkcor+ flatcor+ \
zero=Zero dark=Dark flat=Flat*.fits
p file כ.
cc> ccdlist p//@obj.list
pB67_011.fits[1024,1024][real][object][B][ZDF]:m67b
pB67_012.fits[1024,1024][real][object][B][ZDF]:m67b
pB67_013.fits[1024,1024][real][object][B][ZDF]:m67b
...
Z(Zero correction), D(Dark correction), F(Flat correction) .
4.4. Image Combine
combine image כ .
,
כ כ task imalign.
object image .
pV67 011.fits, pV67 012.fits, pV67 014.fits .
[hshwang@localhost work]$ ls pV67_011.fits pV67_012.fits pV67_014.fits > V67.list
image ( V67 011.fits) . imexam
image imexam 5−10 a
.
cc> epar imexam
(logfile=
V67.cod) logfile
(keeplog=
yes) log output results
cc> disp pV67_011 1 fi+
cc> imexam
3: imalign display
2,3 frame display imexam ,
imexam keeplog parameter no , ds9 frame tile
3 .
cc> epar imexam
(logfile=
(keeplog=
cc> disp pV67_012 2 fi+
cc> disp pV67_014 3 fi+
cc> imexam
) logfile
) log output results

IRAF ABC 11
image a , frame
a , q . xgterm
x, y shift כ . x(ref)-x(img), y(ref)-y(img)
. כ shift כ כ . כ
.
[hshwang@localhost work]$ vi V67_ini.shf
0. 0.
2.69 -1.32
8.32 -5.06
, image image shiftכ
0. image image shift כ כ.
. align .
cc> epar imalign
input = @V67.list Input images
referenc= pV67_011 Reference image
coords = V67.cod Reference coordinates file
output = al//@V67.list Output images
(shifts =
V67_ini.shf) Initial shifts file
:go shift כ
alpV67 *fits image כ. imcombine preprocessing
, . align filter כ
cc> imcombine al//@V67.list alpV67.fits combine=average \
reject=avsigclip scale=none rdnoise=RDNOISE gain=GAIN
cc> disp alpV67 3 fi+
4.5.
color image photoshop, IDL DS9
. red, green, blue fits
( align ), . V
r,g,b .
[hshwang@localhost work]$ ds9 -rgb -red -green -blue

12 HWANG & LEE
Current Green Green fits . Red
.
כ, File ⇒
Save Image As כ
. (r,g,b ) Photoshop כ.
4.6. Manuals
IRAF . IRAF
help , (
) .
1. :
− http://atro.snu.ac.kr/ ∼ obscos/
IRAF/PHOTCAL Tutorial - כ
ALLFRAME Tutorial - כ
IRAF/ELLIPSE Tutorial - כ
− http://atro.snu.ac.kr/ ∼ hshwang/
Preprocess - image כ
DAOPHOT - psf כ
PHOT - כ
spectroscopy - כ כ
2. :
− http://iraf.noao.edu/iraf/web/docs/docmain.html
Everything about IRAF