KMTNet
Korea Microlensing Telescope Network   


모자이크 CCD 영상 WCS 업데이트 방법

2015-09-04 / 오전 11:09:52
[첨부] KMTNet 데이터처리 스크립트
1) 키워드 업데이트
– kmtc*.fits화일을 모아놓고 perl upddetsec2를 실행한다.

2) login.cl 화일을 열고 아래를 추가한다.
– set stdimage=imt55
– set disable_wcs_maps=””
– flpr

3) crval1, crval2 업데이트
– filename은 MEF 영상 1개를 말함.
– mscred> mscsetwcs filename database=18kb.db

4) 카탈로그 만들기
– mscred> mscgetcatalog filename output=filename.cat magmin=0 magmax=20
– mscred> !awk ‘{if($3>14 && $3<16) print $0}’ filename.cat > filename.bri

5) zero offset 값 계산
– mscred> msczero filename ext=im8 coords=filename.bri update+ label+
– 모자이크영상의 offset을 계산합니다. MEF extension im8의 맨 아래쪽이 광축과 가장 가까운 부분입니다. (m) 카탈로그를 읽는다. (s) 일치 시키려는 별에 가장 가까운 좌표를 확인한다. (z) 별의 중심좌표를 카탈로그의 좌표로 재설정한다. 별을 2-3개 선택하여 offset을 구하고, (q)로 종료한다.

6) WCS 업데이트
– mscred> msccmatch filename coords=filename.cat usebpm-
– q로 나온후, 엔터를 쳐서 solution을 업데이트한다.

18k-ccd

7) WCS 확인
– ds9 –mosaicimage iraf filename 으로 영상을 로드한다
– ds9 | Analysis | Image Servers | SAO-DSS 에서 DSS 영상을 얻은 후 좌표를 비교한다.

현재 KMTNet 관측 영상은 한개의 칩에서 8개의 앰프를 이용하여 영상을 읽어내기 때문에 관측된 영상은 총 32개의 영상조각으로 구성됩니다. 종전의 디스토션 계수 데이타베이스 화일은 32개 각각에 대한 정보를 담고 있습니다. 그러다 보니 한개의 칩에서 (8개의 영상조각을 합쳤을 때) 물리적으로는 각 조각영상 간에 간격이 없음에도 불구하고, WCS 업데이트를 마친 영상을 ds9 등에서 WCS좌표계를 도입하여 출력했을때 각 조각영상에 적용한 WCS 계수의 미세한 차이로 인하여 경계부분에서 영상이 불연속적으로 약간씩 벌어지는 현상이 발견되었습니다.
따라서, 이런 부분을 개선하고자 전처리가 끝난 영상을 하나의 물리적인 칩에 대해서는 하나의 영상으로 만든 후 WCS 정보를 업데이트할 수 있도록 WCS 데이타베이스 화일을 다시 구하였습니다. 이를 위해 8개의 앰프로 읽어낸 영상을 하나의 영상으로 합치기 위한 과정을 간단한 스크립트로 작성하였습니다.(상단 첨부)
화일에 표기된 M, K, N, T 는 각각 4개의 9k 칩을 나타내며, DS9 에서 영상을 출력하였을때 보이는 위치별로 칩의 이름을 아래와 같이 표기하였습니다.
+—–+—–+
| M | K |
+—–+—–+
| N | T |
+—–+—–+
스크립 실행 후 msczero와 msccmatch를 실행시키면 9k x 9k 한장에 대해서 WCS정보를 업데이트할 수 있습니다.
  




  


[참고]   정밀 WCS 결정방법 (양유진 박사 제공)      


원자료 위치 : 양유진 박사 홈페이지






* 잘 안되는 부분이나 오류가 있으면 이충욱(leecu@kasi.re.kr)에게 알려 주세요.