[논문 리뷰] MIARMA: An information preserving method for filling gaps in time series. Application to CoRoT light curves

J. Pascual-Granado, R. Garrido외 1명|arXiv (Cornell University)|2014. 10. 03.

Stellar, planetary, and galactic studies참고 문헌 21인용 수 8

한 줄 요약

MIARMA는 별의 광선 곡선에서 원래 주파수 성분을 유지하기 위해 전진-후행 자동회귀이동평균(ARMA) 모델링을 사용하는 시간 시리즈의 새로운 간격 보정 방법이다. 특히 CoRoT 데이터에서 효과적으로 작동하며, 선형 보간보다 더 나은 성능을 보이며 가짜 앨리어스를 제거하고 이전에 발견되지 않은 미세 구조 신호를 드러낸다. 이는 주로 저주파수 진동을 보이는 별들, 예를 들어 Be 항성과 δ Scuti 항성에서 두드러진다.

ABSTRACT

The method here presented intends to minimize the effect of the gaps in the power spectra by gap-filling preserving the original information, that is, in the case of asteroseismology, the stellar oscillation frequency content. We make use of a forward-backward predictor based on autoregressive moving average modelling (ARMA) in the time domain. The method MIARMA is particularly suitable for replacing invalid data such as those present in the light curves of the CoRoT satellite due to the pass through the South Atlantic Anomaly, and eventually for the data gathered by the NASA planet hunter Kepler. We select a sample of stars from the ultra-precise photometry collected by the asteroseismic camera on board the CoRoT satellite: the δ Scuti star HD 174966, showing periodic variations of the same order as the CoRoT observational window, the Be star HD 51193, showing longer time variations, and the solar-like HD 49933, with rapid time variations. We showed that in some cases linear interpolations are less reliable to what was believed. In particular: the power spectrum of HD 174966 is clearly aliased when this interpolation is used for filling the gaps; the light curve of HD 51193 presents a much more aliased spectrum than expected for a low frequency harmonic signal; and finally, although the linear interpolation does not affect noticeably the power spectrum of the CoRoT light curve of the solar-like star HD 49933, the ARMA interpolation showed rapid variations previously unidentified that ARMA interprets as a signal. In any case, the ARMA interpolation method provides a cleaner power spectrum, that is, less contaminated by spurious frequencies. In conclusion, MIARMA appears to be a suitable method for filling gaps in the light curves of pulsating stars observed by CoRoT since the method preserves their frequency content, which is a necessary condition for asteroseismic studies.

연구 동기 및 목표

• 특히 CoRoT 광선 곡선에서 간격으로 인해 발생하는 위조 주파수 앨리어스 문제를 해결한다.
• 저주파수 진동 항성에서 원래 주파수 성분을 유지하지 못하는 선형 보간의 한계를 극복한다.
• 신뢰할 수 있는 항성진동학 분석을 위해 진짜 신호 성분을 유지하는 간격 보정 방법을 개발한다.
• 다양한 항성 유형, 특히 δ Scuti, Be 및 태양 유사 진동 항성에 적용 가능한 방법을 확보한다.
• 기존의 인painting 기법들이 간혹 잡음을 신호로 재구성할 수 있는 잠재적 위험을 피할 수 있는 강력한 대안을 제공한다.

제안 방법

• 시간 시리즈의 누락된 값을 추정하기 위해 자동회귀이동평균(ARMA) 모델 기반의 전진-후행 예측기법을 사용한다.
• 유효한 데이터 세그먼트를 ARMA 과정으로 모델링하여 분석 대상의 시간적 동역학을 포착하며, 해석적 형태를 가정하지 않는다.
• 원래 신호의 스펙트럼 특성을 유지하면서 간격을 반복적으로 보정하는 방식으로 적용한다.
• 모델 적합도 평가 및 잡음과 진짜 신호의 구분을 위해 아카이케 정보 기준(AIC)을 활용한다.
• 가짜 백색 잡음 데이터를 사용하여 방법의 타당성을 검증하여 인위적 에너지가 생성되지 않음을 확인한다.
• 신호 유지 능력과 스펙트럼 정밀도를 평가하기 위해 k-inpainting 및 선형 보간과의 결과를 비교한다.

실험 결과

연구 질문

• RQ1 저주파수 진동 항성의 간격이 있는 광선 곡선에서 선형 보간이 얼마나 심각하게 진짜 주파수 성분을 손상시키는가?
• RQ2 ARMA 기반 간격 보정이 CoRoT 관측 창문 함수로 인한 위조 앨리어스를 제거할 수 있는가?
• RQ3 MIARMA는 기존의 보간 방식에 의해 가려졌던 광선 곡선의 이전에 발견되지 않은 미세 구조를 드러낼 수 있는가?
• RQ4 합성 백색 잡음 데이터에서 MIARMA는 k-inpainting에 비해 잡음의 통계적 성질을 얼마나 잘 유지하는가?
• RQ5 ARMA 보간은 간격이 있는 시간 시리즈에서 내재된 신호와 잡음을 신뢰성 있게 구분할 수 있는가?

주요 결과

• 선형 보간은 δ Scuti 항성과 Be 항성인 HD 174966 및 HD 51193의 파wer 스펙트럼에 심각한 스펙트럼 앨리어스를 유도하여 진짜 진동 주파수를 왜곡한다.
• 태양 유사 항성인 HD 49933의 경우 선형 보간은 파워 스펙트럼에 뚜렷한 영향을 주지 않지만, MIARMA는 이전에 알려지지 않은 빠른 변동을 드러내며 이를 진짜 신호로 해석할 수 있다.
• MIARMA는 항성 유형에 관계없이 파워 스펙트럼의 모든 앨리어스를 성공적으로 제거한다. 특히 선형 보간이 실패하는 장주기 변수에 대해서도 마찬가지다.
• 순수한 백색 잡음 시간 시리즈에 적용했을 때, MIARMA는 표준편차 σ = 0.0753을 가진 보간 값을 생성하여 인위적 신호 재구성이 일어나지 않음을 시사한다.
• 반면 k-inpainting 알고리즘은 동일한 백색 잡음 데이터에서 표준편차 σ = 0.2730을 생성하여 잡음을 신호로 재구성함으로써 신호 유지 원칙을 위반함을 보여준다.
• 아카이케 정보 기준은 MIARMA로 보간한 HD 49933 광선 곡선의 미세 구조가 백색 잡음이 아니라는 것을 확인하여 진짜 물리적 신호로 간주할 수 있음을 뒷받침한다.