ACS/WFC 이미지 위성 궤적 자동 탐지를 위한 개선된 Median Radon 변환

초록

본 논문은 HST ACS/WFC 데이터에서 위성 궤적 및 기타 선형 인공물을 자동으로 식별하기 위해 Median Radon Transform(MRT)을 도입한다. 기존 Radon Transform의 합산 방식 대신 각 경로의 중간값을 사용함으로써 별·은하와 같은 밝은 천체의 영향을 크게 억제한다. 85 %의 회수율과 2.5 %의 오탐률(회절 스파이크 제외)을 달성했으며, 대규모 데이터셋에서 전체 오탐률은 약 10 %로 추정한다. 2002‑2022년 사이 ACS/WFC 이미지에 적용한 결과, 위성 궤적 오염 빈도는 약 2배 증가했지만 평균 밝기에는 뚜렷한 변화가 없었다. 알고리즘은 acstools 패키지에 findsatmrt 함수로 제공된다.

상세 분석

이 연구는 선형 인공물 탐지를 위한 기존 방법들의 한계를 정확히 짚어낸다. 전통적인 Radon Transform(RT)은 모든 가능한 직선을 따라 픽셀 값을 합산해 최대값을 찾지만, 이미지에 별·은하 등 고광도 점원이 존재하면 이들 역시 큰 합산값을 만들어 위성 궤적을 가릴 위험이 있다. 저자들은 이를 해결하기 위해 “Median Radon Transform”(MRT)이라는 변형을 제안한다. MRT는 각 직선 경로에 대해 합계가 아니라 중간값(median)을 계산한다. 중간값은 소수의 극단값에 민감하지 않으므로, 점광원에 의해 발생하는 가짜 신호를 자연스럽게 억제한다.

알고리즘 흐름은 다음과 같다. (1) ACS/WFC의 FL​C 파일을 로드하고, DQ 플래그를 이용해 결함 픽셀을 NaN 처리한다. (2) 배경을 중간값으로 추정해 차감하고, 2×2 재빈도를 적용해 연산량을 10배 가량 감소시킨다. (3) skimage.transform.radon을 기반으로 구현된 MRT를 계산한다. 여기서는 각 경로에 사용된 유효 픽셀 수 L을 동시에 기록해, L < 25인 경우는 신뢰도가 낮아 제외한다. (4) MRT 결과에 대해 S/N > 5인 영역을 photutils.StarFinder와 3가지 폭(15, 7, 3 픽셀) 커널로 검출한다. (5) 검출된 후보를 회전·절단해 1‑D 프로파일을 만든 뒤, 3차 다항식으로 배경을 보정하고 Gaussian 피팅으로 폭·피크 밝기를 재측정한다. 재측정된 S/N이 임계값을 넘고, 폭이 75 픽셀 이하이면 최종 후보로 채택한다. (6) 후보를 일정 길이(≥100 픽셀) 구간으로 나누어 각 구간별 S/N ≥ 3을 확인, 지속성(persistence) 점수를 부여한다. 지속성 점수가 0 이하이면 폐기한다.

이 절차는 두 가지 핵심 장점을 제공한다. 첫째, 중간값 기반이라 별·은하가 산재한 필드에서도 위성 궤적을 효과적으로 구분한다. 둘째, 길이와 S/N을 고려한 지속성 검증으로 부분적인 선형 잡음(예: charge bleeding)이나 짧은 스파이크를 걸러낸다.

성능 평가는 세 단계로 이루어진다. (i) 이론적 감도 분석에서는 경로 길이 L에 따라 최소 검출 밝기가 달라짐을 보여준다. 가장 긴 대각선(L≈4580 픽셀)에서는 배경 잡음의 0.1배 수준까지도 S/N = 5로 검출 가능하고, 코너 근처(L≈500 픽셀)에서는 0.3배 수준이다. (ii) 시뮬레이션에서는 실제 CCD 특성(읽기 잡음, CTE 손실, 플랫필드 비균일성 등)을 포함한 모의 이미지에 인공 위성 궤적을 삽입해, 0.05–0.5 배 밝기 범위에서 회수율과 오탐률을 측정했다. (iii) 실제 ACS/WFC 데이터(2002‑2022)와 Frontier Fields 프로젝트의 눈으로 식별된 궤적을 비교했을 때, 85 % 회수율과 2.5 % 오탐률을 기록했다. 다중 이미지에서 동일 궤적이 반복 나타나는 경우를 이용해 전체 오탐률을 약 10 %로 추정했다.

한계점도 명확히 제시한다. 매우 혼잡한 은하단 필드에서는 별·은하가 밀집해 중간값 자체가 상승, MRT가 여전히 가짜 신호에 민감해진다. 또한, 재빈도 과정에서 궤적 폭이 2 픽셀 이하인 경우 검출 효율이 떨어진다. 저자들은 이러한 상황에서 파라미터(θ, ρ 샘플링, 최소 L, 폭 제한 등)를 조정하거나, 별 마스크를 사전 적용하는 방안을 제안한다.

전반적으로 MRT 기반 자동 탐지는 기존 합산형 RT보다 높은 신뢰도와 낮은 오탐률을 제공하며, acstools에 통합된 findsatmrt 함수는 HST 데이터 파이프라인에 손쉽게 삽입될 수 있다.