세분화된 HPGe 검출기에서 감마 상호작용 위치 추정을 위한 행렬 및 특이값 분해 기법

초록

본 논문은 현대 코액셜·플래너 HPGe 검출기의 신호를 행렬 형태로 모델링하고, 특이값 분해(SVD)를 이용해 겹치는 전류 파형으로부터 감마 상호작용 위치를 빠르고 정확하게 복원하는 방법을 제시한다. 일반화된 행렬 방법과 최적화된 트렁케이션 전략을 통해 잡음에 강인하면서도 실시간 처리에 적합한 알고리즘을 설계하였다.

상세 분석

이 연구는 고해상도 게르마늄 검출기(HPGe)의 전기적 응답을 선형 시스템으로 가정하고, 각 가능한 상호작용 지점에 대한 기본 파형(basis pulse)을 사전 계산하여 응답 행렬 R을 구성한다. 실제 측정 신호 s는 R과 상호작용 강도 벡터 a의 선형 결합 s = R·a 로 표현되며, 여기서 a는 희소(sparse)하고 비음수인 물리량이다. 문제는 R이 매우 큰 차원을 갖고, 특히 세분화된 전극 구조에서는 열악한 조건수(ill‑conditioned)인 경우가 많아 직접 역연산이 불안정해지는 점이다.

특이값 분해(R = U·Σ·Vᵀ)를 적용하면 행렬의 고유 모드와 대응 특이값을 명시적으로 얻을 수 있다. 작은 특이값에 해당하는 모드는 잡음 증폭을 야기하므로, 저특이값을 차단(truncation)하고 유의한 모드만 보존하는 ‘트렁케이션 SVD’가 핵심이다. 저차원 서브스페이스에 투영한 후 최소제곱(L₂) 혹은 L₁ 정규화를 결합하면, 과도한 진동 없이 희소한 상호작용 위치를 복원할 수 있다.

또한 논문은 ‘일반화 행렬 방법(generic matrix method)’을 제시한다. 이는 검출기마다 별도의 파라미터 튜닝 없이, 사전 계산된 응답 행렬과 SVD 트렁케이션 파라미터만으로 실시간 위치 추정이 가능하도록 설계되었다. 알고리즘 흐름은 (1) 신호 전처리·노이즈 필터링, (2) 응답 행렬의 SVD 수행(오프라인), (3) 실시간 신호를 트렁케이션된 Vᵀ에 투영, (4) 비음수 제약을 갖는 최적화 문제 해결, (5) 복원된 강도 벡터에서 임계값 이상인 원소를 상호작용 지점으로 판정하는 단계로 구성된다.

시뮬레이션 결과는 전형적인 코액셜 및 플래너 HPGe 검출기에서 1 mm 이하의 위치 정확도와 1 keV 수준의 에너지 분해능을 동시에 달성함을 보여준다. 특히, 두 개 이상의 상호작용이 시간적으로 겹쳐 있을 때도, 전통적인 피크 피팅 방식보다 3배 이상 빠른 처리 속도와 20 % 향상된 위치 재현성을 기록하였다. 잡음 민감도 분석에서는 SNR ≥ 10 dB 이상이면 트렁케이션 차수 선택에 크게 좌우되지 않으며, 실험적 검증에서도 동일한 성능이 확인되었다.

핵심 인사이트는 (①) 응답 행렬을 사전 계산하고 SVD로 고유 모드를 분리함으로써 잡음 증폭을 억제하고, (②) 트렁케이션 차수를 물리적 해상도와 잡음 수준에 맞게 조정함으로써 실시간 처리와 정확도 사이의 최적 균형을 찾을 수 있다는 점이다. 또한, 비음수 및 희소성 제약을 결합한 최적화는 물리적으로 의미 있는 해를 보장하면서도 계산량을 크게 줄인다. 이러한 방법론은 다중 세그먼트 HPGe 뿐 아니라, 실시간 감마‑카메라, 핵의학 영상, 방사선 방어 시스템 등 다양한 분야에 적용 가능성을 시사한다.