version 1.0 November 2010 PyRAF Y. Nakajima Computer and Data Management Division Subaru Telescope NAOJ

Size: px

Start display at page:

Download "version 1.0 November 2010 PyRAF Y. Nakajima Computer and Data Management Division Subaru Telescope NAOJ"

もえりかんざとばる
5 years ago
Views:

1 version 1.0 November 2010 PyRAF Y. Nakajima Computer and Data Management Division Subaru Telescope NAOJ

2 Chapter 1 PyRAF 1.1 PyRAF PyRAF IRAF Python STScI 1998 (1) IRAF-CL (2) CL-? (3) IRAF Python wrapper / IRAF CL- PyRAF PyRAF Python 1

3 2 CHAPTER 1. PYRAF Python IRAF Python CL Python Python 1.2 PyRAF PyRAF IRAF CL pyraf IRAF PyRAF ~/iraf login.cl PyRAF iraf mkiraf IRAF-CL --> IRAF-CL --> imhead dev$pix long+ dev$pix[512,512][short]: m51 B 600s No bad pixels, min=-1., max= Line storage mode, physdim [512,512], length of user area 1621 s.u. Created Mon 23:54:13 31-Mar-1997, Last modified Sun 16:37:53 12-Mar Pixel file "HDR$pix.pix" [ok] KPNO-IRAF / / IRAF-MAX= E4 / DATA MAX IRAF-MIN= E0 / DATA MIN imstat --> imstat dev$pix # IMAGE NPIX MEAN STDDEV MIN MAX dev$pix ds9 disp -->! ds9 & --> disp dev$pix frame to be written into (1:16) (1): z1=35. z2= IRAF ecl tab

4 1.2. PYRAF 3 PyRAF -->.exit PyRAF logfile -->.logfile[filename[append overwrite]] append logfile.logfile IRAF-CL PyRAF GUI epar PyRAF epar GUI (Unlearn) Help

4 CHAPTER 1. PYRAF GUI --> imstat dev$pix fields=midpt,stddev --> iraf.imstat( dev$pix, fields= midpt,stddev ) --> iraf.imstat.fields= image,midpt,mode --> iraf.imstat( dev$pix ) 1.2.

5 4 CHAPTER 1. PYRAF GUI --> imstat dev$pix fields=midpt,stddev --> iraf.imstat( dev$pix, fields= midpt,stddev ) --> iraf.imstat.fields= image,midpt,mode --> iraf.imstat( dev$pix ) PyRAF xgterm CL xterm imexam radial profile xterm CL xgterm PyRAF OpenGL Tkinter --> prow dev$pix xgterm print save [Window] [New] > i = iraf --> i.imstat.fields= image,midpt,mode --> i.imstat( dev$pix )

6 1.2. PYRAF 5 iraf Stdout=1 --> s = iraf.imhead("dev$pix", long=yes, Stdout=1) --> print s[0] dev$pix[512,512][short]: m51 B 600s --> iraf.head(nl=3, Stdin=s) dev$pix[512,512][short]: m51 B 600s No bad pixels, min=-1., max= Line storage mode, physdim [512,512], length of user area 1621 s.u. Stdin Stderr pyraf -s s silent pyraf -h

8 Chapter 2 PyRAF 2.1 (dev$pix) imstat imstat1.py iraf.images() iraf.imstat("dev$pix") chmod +x imstat1.py PyRAF 1 python python imstat1.py python 3 PyRAF 5 images imstat 6 [ ] python 6 1 7

9 8 CHAPTER 2. PYRAF PyRAF imstat2.py def run_imstat(): iraf.images() iraf.imstat("dev$pix") PyRAF 1 import imstat2 run imstat() imstat2.run imstat() : PyRAF import def imstat-2.py imstat 2.py imstat2.py PyRAF imstat3.py def run_imstat(): iraf.images() iraf.imstat("dev$pix") if name == " main ": run_imstat() chmod +x imstat3.py PyRAF import imstat3 imstat3.run imstat() PyRAF chmod +x mkiraf login.cl login.cl images iraf.images() iraf.imstat() (mkimage.py) 1 login.cl python

10 def run_mkimage(): iraf.imred(_doprint=0) iraf.ccdred(_doprint=0) iraf.ccdtest() iraf.mkimage("test.fits","make","100","2"," ") if name == " main ": run_mkimage() doprint=0 doprint=1 2.2 imstat file.py import os,sys def imstat_file(infile): os.access sys.argv sys.argv[1] infile if os.access(infile, os.r_ok): f=open(infile) for line in f: fname, xx, yy = line[:-1].split( ) if not line.startswith("#"): fitsname=fname+.fits if os.access(fitsname, os.r_ok): iraf.imstat(fname) else: print "Error: can t read", fname f.close() else: print "Error: can t open", infile if name == " main ": imstat_file(sys.argv[1]) imstat file

11 10 CHAPTER 2. PYRAF %more hoge.list ff ff ff ff %imstat_file.py hoge.list (.fits ) imstat # imstat (imstat file2.py) import os,sys def imstat_file(infile): if os.access(infile, os.r_ok): fout=open( imstat_out.txt, w ) f=open(infile) for line in f: fname, xx, yy = line[:-1].split( ) if not line.startswith("#"): fitsname=fname+.fits if os.access(fitsname, os.r_ok): iraf.imstat.fields= image,mean,midpt,stddev out_imstat=iraf.imstat(fname,stdout=1) print >> fout, out_imstat[1] else: print "Error: can t read", fname f.close() fout.close() else: print "Error: can t open", infile if name == " main ": imstat_file(sys.argv[1]) imstat out imstat # IMAGE MEAN MIDPT STDDEV out imstat[0]

12 iraf.imstat(fname,stdout="out.txt",stderr="err.txt") Stdout="STDOUT" Stderr="STDERR" import os,sys os sys python python math, shutil, glob import os,sys,math,shutil,glob def run_myscript(): for file in glob.glob( *.fits ): os.remove(file) x=100 y=math.log10(x) prefix=[ jband, hband, kband ] for band in prefix: image=band fits shutil.copy(image, temp.fits ) os.system( myprog temp.fits ) os.rename( temp.fits, band+ temp.fits ) copy remove iraf imcopy delete os.sysyem UNIX SExtractor

13 12 CHAPTER 2. PYRAF [ ] (1) imstat file.py (2) ( ) (skylevel) (stddev) digiphot.apphot.daofind daofind datapars.fwhmpsf = 2.2 datapars.sigma = stddev datapars.datamin = skylevel-6*stddev datapars.datamax = datapars.readnoise = 30 datapars.epadu = 5 findpars.threshold = 10 (3) daofind apphot.phot daofind datapars phot centerpars.cbox = 5 centerpars.maxshift = 2 fitskypars.annulus = 10 fitskypars.dannulu = 5 photpars.aperture = 9 (4) phot.mag.1 xc, yc, mag, merr, cier, sier, pier txdump

14 python try:except: PyRAF IRAF try:except: (trytest.py) def test_try(): try: iraf.columns("garbage",37) except iraf.iraferror, e: print "error was caught" print e if name == " main ": test_try() try:except: IRAF 2.4 python PyRAF 2 epar mytask mytask.par login.cl uparm inword,s,a,"",,,"string to print" xx,r,a,5.,,,"real value" iy,i,a,1,,,"test integer" mode,s,h,"al" inword xx iy xx 5. iy 1 mytask.py 2 IRAF-CL login.cl package user task myscript ="/home/nakajima/script/myscript.cl"

15 14 CHAPTER 2. PYRAF def mytask(inword,xx,iy): print inword yy=2.0*xx ii=10*iy print "yy=%.1f ii=%d" % (yy,ii) printf parfile = iraf.osfn("home$script/mytask.par") t = iraf.iraftaskfactory(taskname="mytask", value=parfile, function=mytask) iraf.osfn() par mytask PyRAF --> pyexecute("home$script/mytask.py") PyRAF login.cl keep pyexecute("home$script/mytask.py") epar mytask GUI

version 1.0 July 2012 python Yasushi Nakajima Computer and Data Management Division Subaru Telescope NAOJ

version 1.0 July 2012 python Yasushi Nakajima Computer and Data Management Division Subaru Telescope NAOJ Contents 1 3 1.1 python................................... 3 1.2..................... 5 1.3..........................