가우시안 대체 모델, 포아송 역문제에서 뛰어난 MSE 성능

초록

포아송 잡음이 지배하는 저선량 영상 복원에서, 정규화된 가우시안 2차 손실이 포아송 MAP과 동등하거나 더 낮은 평균제곱오차(MSE)를 달성한다는 이론적·실험적 결과를 제시한다.

상세 분석

본 논문은 포아송 잡음이 지배하는 영상 역문제에서 평균제곱오차(MSE)를 최소화하는 것이 포아송 가능도 최대화와 반드시 일치하지 않음을 강조한다. 저선량 영역(예상 카운트 µ≪1)에서는 포아송 최대우도추정량(MLE)이 작은 계수 a_j 때문에 1/(s a_j) 형태의 큰 스파이크를 만들어 분산이 급증하고, 이는 전체 MSE를 크게 악화시킨다. 저자는 이를 완화하기 위해 두 가지 정규화 전략을 제시한다. 첫 번째는 포아송 MAP, 즉 등방성 가우시안 사전분포에 기반한 티크노프 정규화이며, 해석적으로 각 모드별 MSE 비율이 γ_j=τ/(s a_j)^2에 의해 조절되는 식(14)를 도출한다. γ_j가 클수록(즉 a_j가 작을수록) 정규화 효과가 강해져 스파이크가 억제된다. 두 번째는 가우시안 대체 모델이다. 이중 이질분산(heteroscedastic) 가우시안 근사는 포아송을 정규분포로 근사하지만 로그항 때문에 여전히 비선형성을 유지한다. 보다 단순한 동질분산(homoscedastic) 가우시안 모델을 도입하여 y_j∼N(s(Ax)_j,1) 로 가정하고, 티크노프 정규화를 포함한 MAP 추정식을 (17)에서와 같이 닫힌 형태로 얻는다. 이 경우 각 모드의 MSE 비율은 식(18)으로, 저선량 한계에서 (1+γ_j)^{-2} 만큼 감소함을 보인다. 흥미롭게도 이 감소율은 포아송 MAP의 (2/(1+√(1+4γ_j)))^2 보다 더 작아, 동질 가우시안 MAP이 더 강한 수축 효과를 제공한다는 점을 강조한다. 이론적 결과는 A와 x*가 다항식 감쇠를 보이는 경우에도 적용 가능하며, d(s)≈s^{1/(α+β)} 로 정의되는 dose‑dependent 해상도 개념을 도입해 저선량 regime이 언제 발생하는지 정량화한다. 실험적으로는 2D 평행빔 CT 시뮬레이션에서 포아송 MAP, 정규화된 가우시안 MAP, 가중 최소제곱(PWLS) 및 단순 LS를 비교한다. 저선량(평균 카운트 10)에서는 동질 LS와 정규화된 가우시안 MAP이 포아송 MAP과 동등하거나 약간 우수한 MSE를 보이며, 고선량(평균 카운트 1000)에서는 모든 방법이 수렴한다. 이러한 결과는 저선량 상황에서 정확한 포아송 모델보다 정규화와 수축 효과가 MSE 개선에 더 핵심적임을 실증한다.