Data reduction method for OPTICAM multiband time series of transiting exoplanets

이 논문은 OPTICAM 기기의 따뜻한 픽셀 (warm pixels) 문제를 해결하기 위해 3x3 중앙값 필터를 적용한 다대역 시간 계열 데이터 축소 방법론을 제안하고, 이를 통해 TOI-7149 b 의 통과 관측에서 가장 높은 베이지안 증거를 얻어 최적의 전처리 기법으로 검증했습니다.

핵심

1. 배경: 왜 이 연구가 필요할까?

천문학자들은 지구 밖 행성이 별 앞을 지나갈 때, 별빛이 살짝 어두워지는 현상을 관측합니다. 이때 여러 가지 색깔 (파장) 의 빛을 동시에 찍으면, 그 행성의 대기가 어떤 성분으로 이루어졌는지, 혹은 그 행성이 진짜인지 가짜인지 (예: 쌍성계일 수 있음) 를 더 정확하게 알 수 있습니다.

하지만 여기에는 큰 걸림돌이 있었습니다.

• OPTICAM 이라는 장비: 멕시코의 산 페드로 마르트르 천문대에 있는 2.1 미터 망원경에 달린 고성능 카메라입니다.
• 문제점 (따뜻한 픽셀, Warm Pixels): 이 카메라는 매우 민감해서 10 초 이상 사진을 찍으면, 화면에 **예상치 못한 '따뜻한 점들' (Warm Pixels)**이 생깁니다.
  • 비유: 마치 카메라 렌즈에 기름때나 먼지가 낀 것처럼, 사진을 찍을 때마다 그 자리에 빛나는 점들이 생기는 겁니다.
  • 특이점: 이 점들은 일반적인 '어두운 사진 (다크 프레임)'을 빼주는 방법으로도 지워지지 않습니다. 왜냐하면 매번 다른 모양과 위치로 튀어 오르기 때문입니다. 마치 귀신처럼 제자리에서 사라졌다가 다른 곳에 나타나는 것이죠.

이 '따뜻한 점들'이 빛의 흐름 (광도 곡선) 에 섞여 들어가면, 행성의 크기를 재거나 대기를 분석할 때 큰 오차가 생깁니다.

2. 해결책: 6 가지 청소 방법의 대결

연구팀은 이 '따뜻한 점'들을 지우기 위해 6 가지 다른 디지털 청소 방법을 시험해 보았습니다. TOI-7149 b 라는 행성이 별 앞을 지나가는 데이터를 가지고 실험을 했습니다.

1. 기존 방법 (Standard): 그냥 어두운 사진을 빼고 평평하게 만드는 일반적인 방법. (결과: 점들이 그대로 남음)
2. 가우시안 블러 (Gaussian Blur): 이미지를 흐릿하게 만들어 점들을 부드럽게 하는 방법.
   • 비유: 흐린 안경을 끼고 보는 것처럼, 점들을 지우려다 보니 별 자체의 모양도 뭉개져서 별빛이 주변으로 퍼져버렸습니다. 마치 별빛이 번져서 가짜 별을 만들어낸 꼴이 되었습니다.
3. 중앙값 필터 (Median Filter): 주변 픽셀들의 '중간 값'을 가져와서 이상한 점들을 교체하는 방법.
   • 비유: 반찬 그릇에 이상한 돌멩이가 섞여 있으면, 주변의 정상적인 콩들 사이에서 '가장 보통인 콩'을 골라 돌멩이를 대체하는 방식입니다.
   • 연구팀은 이 필터의 크기를 3x3, 5x5, 7x7 픽셀로 바꿔가며 실험했습니다.

3. 승자는? "3x3 중앙값 필터"

실험 결과, **3x3 픽셀 크기의 중앙값 필터 (w3)**가 가장 훌륭한 성과를 냈습니다.

• 왜 3x3 이 좋을까?
  • 너무 작으면 (1x1): 점들을 다 지우지 못합니다.
  • 너무 크면 (5x5, 7x7): 점들은 지워지지만, 별의 빛까지 너무 많이 잘라내거나 뭉개버립니다. 마치 큰 스펀지로 닦으려다 그림 전체를 지워버리는 격입니다.
  • 적당함 (3x3): 점들은 깔끔하게 지우면서, 별의 빛은 거의 손상하지 않습니다.
• 성능: 이 방법을 쓰면 데이터의 '붉은 소음 (Red Noise, 시간에 따라 변하는 잡음)'이 가장 적게 남았고, 행성 모델을 맞출 때 통계적으로 가장 확실한 결과를 얻었습니다.

4. 결론 및 제안: 새로운 청소 도구 (PROFE)

연구팀은 이 발견을 바탕으로 **새로운 자동화 프로그램 (PROFE)**을 만들었습니다.

• 어떻게 작동하나요?
  1. 원본 사진을 받습니다.
  2. 3x3 중앙값 필터를 적용해서 '따뜻한 점'들을 지웁니다.
  3. 그다음 일반적인 보정 (어두운 사진 빼기, 평평하게 만들기) 을 합니다.
  4. 최종적으로 행성이 지나가는 빛의 변화를 정밀하게 계산합니다.

이 방법은 OPTICAM 카메라를 사용하는 모든 천문학자들에게 추천됩니다. 특히 별빛이 약할 때 (어두운 별을 볼 때) 이 '따뜻한 점'들의 영향이 커지는데, 이 방법을 쓰면 그 영향을 크게 줄일 수 있습니다.

요약: 한 줄로 정리하면?

"고성능 카메라에 생기는 **예측 불가능한 잡음 (따뜻한 점)**을 지우기 위해, **3x3 크기의 '중간값 필터'**라는 정교한 디지털 청소 도구를 개발했고, 이를 통해 행성 관측의 정확도를 크게 높였다는 이야기입니다."

이 연구는 마치 사진 편집 프로그램에서 '스마트 필터'를 개발하여, 사진의 노이즈는 지우되 피사체 (별) 는 선명하게 남기는 기술을 천문학에 적용한 사례라고 볼 수 있습니다.

기술 요약

제공된 논문 "Data reduction method for OPTICAM multiband time series of transiting exoplanets"의 기술적 요약입니다.

1. 연구 배경 및 문제 제기 (Problem)

• 장비 및 환경: 멕시코 산 페드로 마르트르 천문대 (OAN-SPM) 의 2.1m 망원경에 장착된 OPTICAM 장비를 사용하여 외계행성 (TOI-7149 b) 의 통과 (transit) 관측을 수행했습니다. OPTICAM 은 3 개의 독립적인 채널에서 동시에 다대역 (multiband) 관측이 가능한 sCMOS (Andor Zyla 4.3-Plus) 검출기를 사용합니다.
• 핵심 문제: sCMOS 검출기는 장노출 (10 초 이상) 시 **워밍 픽셀 (warm pixels)**이 급격히 증가하는 문제가 발생합니다.
  • 이러한 워밍 픽셀은 표준적인 암전 (dark) 차감이나 평평화 (flat-field) 보정으로 제거되지 않습니다. 이는 워밍 픽셀의 밝기가 프레임마다 예측 불가능하게 변하기 때문입니다.
  • 워밍 픽셀은 광도 곡선 (light curve) 에 **적색 잡음 (red noise)**과 백색 잡음을 유발하여 행성 통과 모델의 정밀도와 신뢰성을 떨어뜨립니다.
  • 특히 TOI-7149 b 와 같은 어두운 항성 (M 왜성) 의 경우, 신호 대 잡음비 (S/N) 가 낮아 워밍 픽셀의 영향이 더욱 치명적입니다.

2. 방법론 (Methodology)

저자들은 워밍 픽셀의 영향을 완화하기 위해 6 가지 전처리 (pre-processing) 방법을 비교 평가했습니다. 모든 방법은 과학 이미지와 보정 이미지 (암전, 평평화) 에 적용되었습니다.

1. 표준 감축 (Standard Reduction, st): 암전 차감과 평평화 나눗셈만 수행 (기저선).
2. 가우시안 컨볼루션 (Gaussian Convolution):
   • $\sigma=1$ 픽셀 (g1) 및 $\sigma=3$ 픽셀 (g3) 커널을 사용하여 이미지를 평활화 (smoothing) 하는 방법.
3. 중앙값 필터 (Median Filter):
   • $3\times3$픽셀 (w3),$5\times5$픽셀 (w5),$7\times7$ 픽셀 (w7) 윈도우 크기를 사용하여 각 픽셀을 주변 이웃 픽셀의 중앙값으로 대체하는 방법.

평가 프로세스:

• TOI-7149 b 의 $g', r', i'$ 대역 통과 데이터를 사용하여 각 감축 방법을 적용한 후, **AstroImageJ (AIJ)**를 통해 차등 광도 측정을 수행했습니다.
• juliet 패키지를 사용하여 다대역 통과 모델을 동시에 피팅했으며, 각 밴드별로 가우시안 프로세스 (GP) 를 도입하여 상관 잡음을 모델링했습니다.
• 성능 지표: 베이지안 증거 (Bayesian Evidence, $\ln Z$), 광도 정밀도 (RMS), 10 분 구간 상관 잡음 (RMS10 min), $\chi^2_{\nu}$ 등을 비교했습니다.

3. 주요 결과 (Results)

• 워밍 픽셀의 특성: 30 초 노출 시 채널 1 은 약 17-20%, 채널 2 는 30-35%, 채널 3 은 23-60% 의 픽셀이 워밍 픽셀로 분류되었습니다. 이 픽셀들은 시간에 따라 예측 불가능하게 변화하여 고정된 픽셀 마스크로는 제거가 불가능했습니다.
• 가우시안 필터의 한계: 가우시안 컨볼루션 (특히 g3) 은 워밍 픽셀을 제거하는 것처럼 보이지만, 별의 신호를 확산시켜 인접 픽셀로 퍼뜨립니다. 이로 인해 **가짜 별 (fake stars)**이 생성되거나 광도 측정 구간에 불필요한 플럭스가 추가되어 시스템 잡음을 유발했습니다.
• 중앙값 필터의 우수성:
  • **$3\times3$ 중앙값 필터 (w3)**가 모든 대역에서 가장 우수한 성능을 보였습니다.
  • 베이지안 증거 ($\ln Z$): w3 방법이 다른 모든 방법보다 통계적으로 유의미하게 높은 증거 값을 기록했습니다 ($\Delta \ln Z \gtrsim 20$).
  • 잡음 감소: w3 은 10 분 구간 RMS (상관 잡음 지표) 를 가장 낮게 줄였으며, 시간 평균 곡선 (time-averaging curve) 에서 백색 잡음 기대치에 가장 근접했습니다.
  • 신호 보존: $3\times3$필터는 워밍 픽셀을 효과적으로 제거하면서도 별의 피크 신호를 약 9$5\times5, 7\times7$) 나 가우시안 필터보다 별의 형태 왜곡을 최소화했습니다.
• 행성 파라미터: w3 방법을 통해 얻은 행성 반지름 ($R_p/R_*$) 및 궤도 경사각 등의 파라미터는 기존 연구 (Kanodia et al. 2025) 와 일치하면서도 더 정밀한 오차 범위를 보였습니다.

4. 기여 및 제안된 파이프라인 (Key Contributions & Pipeline)

• 최적화 방법 제안: OPTICAM sCMOS 데이터 처리를 위해 $3\times3$ 중앙값 필터를 표준 감축 전의 필수 전처리 단계로 제안했습니다. 이는 워밍 픽셀 비율이 최대 60% 에 달하는 환경에서 가장 균형 잡힌 해결책입니다.
• 자동화 파이프라인 개발:
  • PROFE: Python 기반의 오픈소스 모듈로, 원시 데이터의 전처리 (중앙값 필터 적용), 메타데이터 업데이트, 진단 생성 (히스토그램, 적색 잡음 곡선 등) 을 담당합니다.
  • AstroImageJ (AIJ): 표준 보정 (마스터 암전/평평화 생성, 차등 광도 측정) 을 수행합니다.
  • 이 두 도구를 결합하여 OAN-SPM 2.1m+OPTICAM 관측 데이터를 고품질 통과 광도 곡선으로 변환하는 통합 파이프라인을 구축했습니다.

5. 의의 (Significance)

• sCMOS 기술의 천문학 적용: 고화질, 고속 판독이 가능한 sCMOS 검출기가 천문학에서 널리 사용되고 있지만, 워밍 픽셀 문제는 큰 걸림돌이었습니다. 이 연구는 이러한 검출기의 단점을 보완하여 고정밀 외계행성 통과 관측을 가능하게 하는 구체적인 데이터 감축 전략을 제시했습니다.
• 다대역 관측의 정밀도 향상: 다대역 동시 관측은 행성 대기 분석 및 위양성 (false positive) 식별에 필수적입니다. 제안된 방법은 특히 신호가 약한 어두운 항성 관측 시 잡음을 효과적으로 제어하여, 지상 기반 망원경으로도 우주 기반 망원경에 버금가는 정밀도의 대기 특성 분석을 위한 데이터 품질을 확보할 수 있음을 입증했습니다.
• 오픈소스 도구 제공: PROFE 코드를 공개하여 다른 연구자들이 OPTICAM 데이터를 효율적으로 처리하고 재현 가능한 과학 결과를 도출할 수 있도록 기여했습니다.

결론적으로, 이 논문은 OPTICAM 의 sCMOS 검출기에서 발생하는 예측 불가능한 워밍 픽셀 문제를 $3\times3$ 중앙값 필터링을 통해 효과적으로 해결함으로써, 외계행성 통과 관측의 정밀도와 신뢰성을 획기적으로 높인 방법론과 도구를 제시했습니다.