슬릿리스 스펙트럼 오염 제거를 위한 새로운 방법

초록

본 논문은 Euclid 등 슬릿리스 분광을 사용하는 우주 임무에서 발생하는 스펙트럼 오염을 해소하기 위해, 두 가지 물리 모델(즉시 선형 모델과 컨볼루션 모델)을 기반으로 한 네 가지 새로운 로컬 탈오염 기법을 제안한다. 다중 분산 방향(0°, 180°, 184°, −4°)을 동시에 활용하고, 직접 이미지(포톤미터) 정보를 이용해 혼합 행렬을 추정한다. 시뮬레이션 실험을 통해 제안 방법이 기존 NMF‑기반, GRIZLI, LINEAR 방식보다 빠르고, 초기값에 민감하지 않으며, 미검출 오염원에 대해서도 강인함을 보였다.

상세 분석

이 연구는 슬릿리스 분광에서 발생하는 ‘오염(Contamination)’ 문제를 신호 분리(Source Separation) 관점에서 재정의한다. 먼저, 관측된 스펙트럼을 두 단계로 모델링한다. (1) 즉시 선형 모델에서는 물체의 공간 프로파일 I_j(x,y)를 x와 y 방향으로 분리 가능한 형태 I_j1(x)·I_j2(y)라고 가정하고, 각 픽셀의 관측값을 a_j(m)·e_j(n) 형태의 곱으로 표현한다. 여기서 a_j는 세로 방향(교차 분산) 프로파일, e_j는 파장에 따라 변하는 스펙트럼 컨볼루션을 포함한다. 이 곱셈 구조는 관측 행렬 X(d_i)=A(i)·E(i) 로 정리될 수 있으며, A(i)와 E(i)는 각각 혼합 행렬과 스펙트럼 행렬이다.

(2) 컨볼루션 모델은 파장 의존 필터(분산 함수 D(λ)=αλ+β)와 PSF h(x,y)를 명시적으로 포함한다. 식 (4)에서 보듯이, 관측 스펙트럼은 λ에 대한 적분 형태의 컨볼루션으로 기술되며, 이는 푸리에 영역에서 곱셈으로 변환될 수 있다. 따라서 Fourier 도메인에서 혼합 모델을 구성하면, 각 파장대의 필터링 효과를 정확히 보정하면서 동시에 스펙트럼을 복원할 수 있다.

두 모델에 기반한 네 가지 탈오염 기법은 다음과 같다.

1. 즉시 모델 – 비음수 행렬 분해(NMF) 기반: 직접 이미지에서 얻은 a_j를 초기값으로 사용하고, 비음수 제약 하에 교대 최소화로 A와 E를 추정한다.
2. 즉시 모델 – 교대 최소화(ALS) with 정규화: 초기값에 민감한 NMF를 보완하기 위해, a_j를 고정하고 e_j를 선형 최소제곱으로 풀며, 정규화 파라미터로 잡음에 대한 강인성을 강화한다.
3. 컨볼루션 모델 – 푸리에 기반 직접 해법: 관측 스펙트럼을 푸리에 변환한 뒤, 각 방향별 필터 H_i(λ)와 혼합 행렬 A(i)를 이용해 역변환 없이 스펙트럼을 복원한다.
4. 컨볼루션 모델 – 정규화된 최소제곱(Constrained LS): 푸리에 영역에서 L2 정규화를 추가해 잡음 억제와 미검출 오염원에 대한 완화 효과를 얻는다.

핵심적인 기술적 통찰은 다중 분산 방향을 동시에 활용한다는 점이다. 기존 방법들은 보통 하나의 방향만 처리하거나, 방향별 결과를 후처리한다. 여기서는 네 방향의 관측 행렬을 하나의 큰 블록 행렬로 결합해, 각 방향의 a_j와 e_j가 공유되는 구조를 만들었다. 이는 동일 물체에 대한 서로 다른 오염 패턴을 상호 보완적으로 이용함으로써, 특히 겹쳐진 스펙트럼이 심한 밀집 필드에서 복원 정확도를 크게 향상시킨다.

또한, 직접 이미지(포톤미터) 활용이 중요한 역할을 한다. 직접 이미지에서 추출한 공간 프로파일 a_j는 물체의 위치와 형태 정보를 제공하며, 이는 혼합 행렬 A(i)의 초기값이자 제약조건으로 작용한다. 따라서 모델은 ‘데이터‑드리븐’ 방식이면서도 물리적 제약을 유지한다.

실험에서는 Euclid‑like 시뮬레이션 데이터를 사용해, 신호‑대‑잡음비(SNR) 520 dB 범위에서 네 방법을 평가하였다. 결과는 평균 스펙트럼 오차가 기존 NMF‑기반 방법 대비 30 % 이상 감소하고, GRIZLI와 LINEAR에 비해 연산 시간은 510배 빠르며, 미검출 오염원에 대한 잔차는 2배 이하로 감소함을 보여준다.

전반적으로 이 논문은 로컬 접근법을 통해 병렬 처리 가능성을 확보하고, 컨볼루션 모델을 푸리에 영역에서 효율적으로 풀어내어, 기존 전역적 방법들의 계산 복잡도와 초기값 의존성을 극복한다는 점에서 큰 의의를 가진다.