VI‑PRISM: 저선량 퍼퓨전 CT를 위한 단일 물질 재구성 혁신

초록

본 논문은 기존 단일 물질 폴리크로매틱 포톤 카운팅 CT 재구성 알고리즘을 퍼퓨전 CT에 적용한 VI‑PRISM을 제안한다. 정적 배경 조직을 사전 스캔으로 알려진 상태에서 요오드 농도 맵만을 추정하며, 다양한 광자 예산과 투시 각도 수에 대해 필터드 백프라젝션(FBP)과 비교하였다. 실험 결과, 광자 수가 10²10⁵ 수준, 투시 수가 8984 범위에서도 VI‑PRISM은 RMSE와 SNR 면에서 FBP보다 우수하고, 요오드 농도 오차를 0.4 mg/ml 이하로 유지한다.

상세 분석

본 연구는 폴리크로매틱 포톤 카운팅 CT(PCCT)에서 물질 특성을 정확히 분리하기 위해 제안된 단일 물질 변분 불평등(VI) 기반 재구성 알고리즘을 퍼퓨전 CT에 맞게 확장한 VI‑PRISM을 소개한다. 핵심 아이디어는 두 개의 물질(배경 조직과 요오드) 중 배경을 사전 스캔으로 고정하고, 동적 요오드 농도만을 변수로 두어 비선형 포톤 카운팅 모델을 직접 해결한다는 점이다.

1. 수학적 모델링

   • 포톤 카운팅 데이터는 포아송 분포를 따르지만, 알고리즘은 기대값 관계식(1)을 기반으로 한다. 기대값은 에너지별 감쇠 계수와 투사 행렬을 이용한 지수식 형태이며, 배경 물질의 기여는 고정된 파라미터로 처리한다.
   • 요오드 농도 맵 (x_{\text{iodine}})에 대한 연산자는 식(4)의 형태로 정의된 단조 연산자 (F(x))이며, 이는 손실 함수의 그래디언트가 아니라 데이터와 모델 간 불일치를 측정한다. 단조성 보장은 적절한 스텝 사이즈 (\alpha_t) 하에서 수렴성을 이론적으로 확보한다.

2. 알고리즘 구조

   • 업데이트 식 (3)은 단순한 전진-후진 형태로, 현재 추정치에서 연산자 (F)를 빼고, 비음수와 TV 제약을 만족하도록 투사한다. TV 제약은 실제 요오드 분포가 거의 조각상(piecewise‑constant)임을 이용한 정규화이며, TV 반경은 실제 지상 진실 TV와 동일하게 설정해 과도한 평활화를 방지한다.
   • 기존 EXACT 알고리즘과 비교했을 때, VI‑PRISM은 업데이트가 간단해 경험적으로 빠른 수렴을 보이며, 다중 물질(배경) 존재를 허용한다는 점에서 퍼퓨전 CT에 적합하다.

3. 시뮬레이션 설계

   • 디지털 팬빔 CT 시뮬레이션은 TIGRE 프레임워크를 이용해 ACR 팬형 팬빔 기준으로 구현했으며, 100 keV까지의 스펙트럼을 2 keV 간격으로 이산화했다. 포톤 카운트는 총 예산 (B)를 투시 수 (n_s)에 균등 분배해 평균 (\bar I = B/n_s) 로 설정, 10²~10⁵ 범위의 저선량 상황을 재현한다.
   • 투시 수는 984부터 8까지 10단계로 감소시켜, 각 예산에서 뷰와 광자 수 간 트레이드오프를 정량화했다. 각 조건을 9번의 무작위 시드로 반복해 통계적 신뢰성을 확보했다.

4. 비교 대상 및 평가 지표

   • 기준은 전통적인 FBP이며, FBP는 포톤 카운트를 로그 변환 후 Hann 필터링을 적용한다. FBP는 에너지별 재구성 후 스펙트럼 가중 평균을 수행한다.
   • 평가 지표는 전역 RMSE(HU), ROI별 RMSE(mg/ml), SNR, 그리고 노이즈 레벨이며, ROI는 8개의 요오드 삽입물(0.05~2.47 mg/ml) 내부의 중심 원을 사용했다.

5. 핵심 결과

   • 모든 광자 예산에서 VI‑PRISM은 FBP 대비 RMSE가 평균 3050 % 감소하고, SNR은 24 dB 향상했다. 특히 10⁴~10⁵ 광자 수준에서 984뷰 전체를 사용할 경우, 요오드 농도 오차는 0.12 mg/ml 이하로 매우 정확했다.
   • 극단적인 저선량(10² 광자)·희소 뷰(8뷰) 상황에서도 VI‑PRISM은 요오드 오차를 0.38 mg/ml 이하로 유지했으며, FBP는 오차가 1 mg/ml을 초과해 실용적 사용이 어려웠다.
   • 통계적 유의성 검증(Bonferroni 보정) 결과, 동일 예산에서 뷰 수를 감소시켜도 VI‑PRISM의 성능 저하가 미미했으며, 오히려 TV 제약이 희소 데이터에서 과적합을 방지하는 효과가 확인되었다.

6. 의의와 한계

   • VI‑PRISM은 물리적 제약(비음수, TV)과 단조 연산자를 결합해 비선형 포톤 카운팅 모델을 직접 해결함으로써, 저선량·희소 뷰 환경에서도 정량적 요오드 맵을 정확히 복원한다. 이는 퍼퓨전 CT에서 방사선 피폭을 크게 감소시킬 수 있는 실질적 근거를 제공한다.
   • 현재는 시뮬레이션 기반 평가에 머물며, 실제 PCD 하드웨어의 에너지 비선형성, 교차 톤, 그리고 움직임 아티팩트 등에 대한 검증이 필요하다. 또한, 다중 시간 프레임을 동시에 재구성하는 동적 모델링(예: Kalman 필터와 결합)으로 확장하면 실제 임상 퍼퓨전 워크플로우에 더 가까워질 것이다.