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.

Presentation on theme: "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."— Presentation transcript:

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

2 분광자료처리 전처리 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

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

4

5 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“ 로 수정!

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

7 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 .....

8 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> \ oversca+ biassec=[*,2049:2060] ccdred> ccdlist ot*fits .... otbias1.fits[317,2060][real][none][][O]:BIAS ....

9 Bias subtraction – (1) # IMAGE NPIX MEAN STDDEV MIN MAX
ccdred> imstat otbias*.fits # IMAGE NPIX MEAN STDDEV MIN MAX otbias1.fits otbias2.fits otbias3.fits : 값이 이상한 것이 있으면 그 영상은 제거 Combine Bias images ccdred> epar zerocombine (combine= average) Type of combine operation (ccdtype= ) CCD image type to combine : wq ccdred> zerocombine output=mbias

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

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

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

13 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 output=mflat

14 Flat combine – (3) ecl>noao >twodspec>longslit response task longslit> imarith tmflat 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)‏

15 Bias Flat Flat after response Object

16 Trim the images ecl Contents C0001[*,10:2040] E0001 C0002[*,10:2040] E0002 C0003[*,10:2040] E0003 C0004[*,10:2040] E0004

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