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

Presentation on theme: "2010. 11.09 함지범 hjb0220@snu.ac.kr hjb0220@astro.snu.ac.kr 관측2 분광파트 전처리 실습(2) 2010. 11.09 함지범 hjb0220@snu.ac.kr hjb0220@astro.snu.ac.kr."— Presentation transcript:

1 2010. 11.09 함지범 hjb0220@snu.ac.kr hjb0220@astro.snu.ac.kr
관측2 분광파트 전처리 실습(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 Aperture Extraction – (1)
apall task (noao.twodspec.apextract) Determining dispersion axis => two > epar apextract dispersion axis =2 ap > epar apall input (format = onedspec) nfind = 1

4 Aperture Extraction – (2)
- 이렇게 apall을 수행하면 iraf 창이 하나가 뜬다. Determine aperture boundary – l, u Enter the background mode – b Determine background region – s Increase fitting order – : order 4 Fit the function – f 위의 정보를 참고해 구경추출을 한다.

5

6 Wavelength calibration – (1)
1) 비교광원의 선 테이블을 먼저 준비(TU manual) 2) 비교광원을 구경추출(apall task) 3) identify task를 이용해서 파장 정보를 입력 4) Science image의 헤더를 약간 수정한 후(hedit task), 5) dispcor task 를 이용해서 파장을 일치시키면 끝!

7 Wavelength calibration – (2)
ap > epar apall input => 비교광원 이미지 (output => 아웃풋 이미지 (referen (profile => Science 이미지 interact, find, recente, resize, edit, trace fittrac = no extract = yes extras, review = no, background = none

8 Wavelength calibration – (3)
ap > epar identify images = ~~~~~~~.0001.fits (order = 4) => 맘대루.. 이 task를 실행하면 또 그림이 나온다. m => 라인의 파장 값을 입력 l => 자동 파장 맞춤 f => Fitting d => 맘에 안드는 거 지우기

9

10 Wavelength calibration – (4)
위의 과정을 반복하거나, reidentify 명령어를 이용해 다른 calibration lamp에도 적용한다. 이제 Science image에 이 보정 파장 값을 적용할 차례. (헤더에 파장 정보 입력) hedit %#$%.0001 refspec1 add+ ver- Science image calibraton frame linearlize = no

11 Flux calibration – (1) 어려움. 얻은 데이터들을 모두 w*.fits로 저장하고 a.lst로 리스트를 만든다. 이 작업은 모두 eso.set이라고 하는 스크립트를 이용하기 위함이다. 스크립트 내용은 각자 해석하도록 합니다. cl > cl < eso.set 스크립트가 실행된다.

12 Flux calibration – (2) standard task (noao.onedspec)
표준성 정보를 실제 데이터베이스에서 불러옴. on> epar standd input = w (wavelength calibration까지 다 된 표준성) output = std (extinct = onedstds$ctioextinct.dat) (caldir = onedstds$/spec50cal/ starname = feige67 Page onedstds$README 쳐 보세요~

13 Flux calibration – (3) sensfunc task =>파장별 sensitivity보정을 해주기 위해 모델 만듬. on > sensfunc standard = std (extinct =onedstds$ctioextinct.dat) (newexti=extinct.dat) (graphs = srei)

14 Flux calibration – (4) calibrate task
최종적으로 Science 이미지에 보정된 정보를 입력하 는 과정 on > epar calibrate input output (extinct = onedstds$/ctioextinct.dat) Fnally, x-axis : wavelength [A] y-axis : flux [erg/cm^2/s/A]

15 Plotting spectroscopic data
noao.onedspec splot aaa.fits

16 Other correction Heliocentric correction, Dereddening, Removing sky line TU manual 참조. Heliocentric correction Rvcorrect ima=f0020 imupdate+ observa=eso Dopcor f0020 h0020 –vhelio isvel+

17 모든 자료는 TU 매뉴얼을 기초로 하였고, TU매뉴얼을 기초로 한 한글 매뉴얼을 참고하였음. 한글 매뉴얼대로 하다보면 잘 안되기도 함. 어쨌든 나누어준 두 개의 매뉴얼(iraf abc, an introduction to analysis of single dispersion spectra with IRAF)과 TU매뉴얼을 이용, 잘 연구해보세요. SNUO관측 데이터 처리는 TU데이터 처리와 다소 다 를 것임. 나중에 데이터 처리할 때 같이 연구해보아 요… 모르면 항상 Help command 쳐 보는 습관을 기르시길.