X‑레이 CT 촬영의 확률적 모델과 기능적 중심극한정리

초록

본 논문은 단위 원판을 통과하는 모든 직선을 인덱스로 하는 이산 확률 과정을 정의하고, 이를 X‑레이 CT에서 관측되는 광자 수와 연결한다. 먼저 L² 수렴을 보이는 대수의 법칙을 증명하고, 이후 평균 광자 수와 스캔 라인 수의 성장 비율에 따라 다양한 정규화 항을 추가한 기능적 중심극한정리를 제시한다. 특히, 보정항을 단계적으로 삽입함으로써 요구되는 평균 광자 수의 차수를 임의로 낮출 수 있음을 보이며, 각 정리마다 베리‑에센 정리를 통해 수렴 속도도 정량화한다. 또한, 광자 수가 0이 되는 경우를 처리하기 위한 세 가지 정규화 방법(N1‑N3)을 제안하고, 결과에 미치는 영향을 분석한다.

상세 분석

이 연구는 X‑레이 CT의 물리적 현상을 확률론적으로 모델링함으로써 기존의 Beer‑Lambert 법칙을 평균값 수준에서만 사용하는 한계를 극복한다. 저자들은 단위 원판 D 내부의 비음수 연속함수 f를 attenuation 함수로 두고, 각 직선 y∈Z에 대해 X‑ray 변환 Xf(y)=∫_{Ly∩D}f·dx를 정의한다. 실제 촬영에서는 각 직선마다 평균 N개의 광자를 방출하지만, 실제 방출 수 V_y는 포아송(N) 분포를 따르고, 각 광자는 독립적으로 흡수 여부가 Bernoulli(p_y) (p_y=e^{-Xf(y)}) 로 결정된다. 따라서 관측된 흡수되지 않은 광자 수 S_y는 포아송(Np_y) 를 따른다.

핵심 난관은 S_y가 0이 될 확률이 양수이므로 로그를 직접 취할 수 없다는 점이다. 이를 해결하기 위해 (N1) 모든 경우에 1을 더하는 가장 단순한 정규화 방식을 채택하고, (N2) 0일 때만 1을 더하는 방식, (N3) 0이면 재시도하는 방식도 제시한다. 본문에서는 (N1)을 기본으로 삼아 Y_{n,m,N}^f(y)=−log