2010. 11.02 함지범 hjb0220@snu.ac.kr hjb0220@astro.snu.ac.kr 관측2 분광파트 전처리 실습 2010. 11.02 함지범 hjb0220@snu.ac.kr hjb0220@astro.snu.ac.kr

분광자료처리 전처리 Preprocessing Flux Cal. 분광관측 스펙트럼 추출 (Remove cosmic ray)‏ → Bias Dark Flat Object STD star Calibration lamp 전처리 Preprocessing Overscan correction combine bias Bias subtraction combine dark Dark subtraction combine flat (Remove cosmic ray)‏ combine object Flat fielding 스펙트럼 추출 Extract spectra :2D → 1D Wavelength calibration Extinction correction → Radial velocity Flux Cal. → Age Metallicity

샘플 데이터 준비 http://ulisse.pd.astro.it/Astro/AsMon/Vol_1/ 압축 풀기 tar xzvf spectra.tar.gz 이름 바꾸기(확장자가 fits여야 ds9이 읽음) csh rename.csh (a00* -> a00*.fits) (c-shell script. 파일은 홈페이지에 있음....) 각각의 자료에 관한 정보는 ps파일에 있음.

IRAF 준비 $ xgterm –sb & $ mkiraf Enter terminal type: xgterm -> 폴더에 login.cl 이라는 파일이 생김. $ vi login.cl 혹은 $ emacs login.cl 명령으로, #set stdimage = imt800 → set stdimage = imt4096 #set imtype = "imh" → set imtype = "fits“ 로 수정!

Bits correction 샘플 데이터의 값은 -32768~+32768의 값을 가짐. 0~+65536의 pixel value range로 설정해야 함. 모든 픽셀에 +32768을 더해줌. ecl > imarith @list_a + 32768 혹은 ecl > epar imarith에서 변수 변경 후 :go 혹은 :wq.

Overscan correction – (1) ecl> disp tflat1 3 fi+ ecl> imexam (a, e, s, r ..... q )‏ implot tflat1.fits ( :c ??? :a ?? :l ?? , :x ? ? , :y ? ? ..... )‏ → overscan range : [*,2049:2060] cf) trim section : [28:315,10:2040] IRAF/ccdred package 로 가기 ecl> noao : ...... noao> imred : ..... ccdred. ..... ccdred> ? ..... ccdproc darkcombine zerocombine flatcombine .....

Overscan correction – (2) ccdred> epar ccdproc (lpar task, help task)‏ (ccdtype= ) CCD image type to correct (fixpix = no) Fix bad CCD lines and columns? (oversca= no) Apply overscan strip correction? (trim = no) Trim the image? (zerocor= no) Apply zero level correction? (darkcor= no) Apply dark count correction? (flatcor= no) Apply flat field correction? :wq $ cat ????.list > bdfos.list ccdred> ccdproc @bdfos.list output=o//@bdfos.list \ oversca+ biassec=[*,2049:2060] ccdred> ccdlist ot*fits .... otbias1.fits[317,2060][real][none][][O]:BIAS ....

Bias subtraction – (1) # IMAGE NPIX MEAN STDDEV MIN MAX ccdred> imstat otbias*.fits # IMAGE NPIX MEAN STDDEV MIN MAX otbias1.fits 653020 -0.01422 5.153 -27.28 497.7 otbias2.fits 653020 -0.01343 5.116 -29.29 463.7 otbias3.fits 653020 -0.01282 5.347 -27.3 1287. : 값이 이상한 것이 있으면 그 영상은 제거 Combine Bias images ccdred> epar zerocombine (combine= average) Type of combine operation (ccdtype= ) CCD image type to combine : wq ccdred> zerocombine o//@bias.list output=mbias

Bias subtraction – (2) Dark, Flat, Object, STD 영상에서 Bias 제거하기 $ cp bdfos.list dfos.list vi dfos.list ( tbias*.fits 삭제 )‏ ccdred> ccdproc o//@dfos.list output=z//@dfos.list \ zerocor+ zero=mbias ccdred> ccdlist zt*.fits .... ztdark1.fits[317,2060][real][none][][OZ]:DARK .....

Dark subtraction Bias subtraction과 비슷. Dark image들을 합성할 때는 darkcombine사용 !!!!중요!!!! ccdproc task에서 zerocor을 켜고 bias 이미지를 입력하여, Dark이미지가 자동으로 bias보정되게 해야함… 혹은 ecl>imcombine 과 ecl > imarith task를 사용하여, 단순하게 dark와 bias를 만들 수 있다!

Flat combine – (1) 측광 자료 처리와는 다소 다름. long slit의 경우 공간 좌표가 1차원임. (파장정보가 1차원이 들어가기 때문에…) Flat은 원래 각각의 위치에 따른 sensitivity를 보정을 해 주는 것인데, 분광 자료의 경우 위치가 오직 1차 원적임. 따라서 파장에 따른 의존도가 생기고, 이 파장 의존 도를 없애주기 위한 작업이 필요함..

Flat combine – (2) ccdred> imstat dtflat*.fits Combine Flat images ccdred> epar flatcombine (combine= median) Type of combine operation (reject = avsigclip) Type of rejection (ccdtype= ) CCD image type to combine (process= no) Process images before combining? (subsets= no) combine images by subset parameter? (scale = none) Image scaling : wq ccdred> flatcombine d//@flat.list output=mflat

Flat combine – (3) ecl>noao >twodspec>longslit response task longslit> imarith tmflat + 20400 tmflat0 longslit> response tmflat0 tmflat0[50:250,*] rmflat Fit the normalization spectrum for tmflat0 interactively (yes): Dispersion axis(1=along lines, 2=columns, 3=z) (1:3) (2): 2 : order 20, f(fitting), q(quit)‏

Bias Flat Flat after response Object

Trim the images ecl > imcopy @list_d @list_e Contents of lists @list_d @list_e C0001[*,10:2040] E0001 C0002[*,10:2040] E0002 C0003[*,10:2040] E0003 C0004[*,10:2040] E0004

NEXT WEEK Aperture extraction(apall task) Wavelength calibration(identify task) Flux calibration(several….. Tasks…) Another corrections…… - Heliocentric correction - Dereddening - Removing sky lines…