ENGLISH
모자이크 CCD 영상 WCS 업데이트 방법
본문
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을 업데이트한다.
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정보를 업데이트할 수 있습니다.
