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"

もえりかんざとばる
7 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

Microsoft PowerPoint - テキスト-Web掲載版.pptx

Microsoft PowerPoint - テキスト-Web掲載版.pptx 1 IRAF/PyRAF インストール講習会磯貝瑞希国立天文台天文データセンター 2017.12.19 ウェブ掲載版 2 はじめに本日の講習会の目的 : IRAF/PyRAFをシステム (Linux, CentOS 7) にインストールできるようになること講習内容 : 1. IRAF (v2.16.1) のインストール IRAF, STSDAS/TABLES, x11iraf, SAOimageds9