에너지 임계값 편향 계산기로 구현한 광자계산 CT 스펙트럼 보정

초록

본 논문은 광자계산 CT(PCCT)에서 발생하는 스펙트럼 불일치를 두 가지 물리적 요인, 즉 고정된 스펙트럼 스키우 항과 가변적인 에너지 임계값 편향(ETB) 항으로 모델링한다. 캘리브레이션 단계에서 이 두 항을 추정한 뒤, 실시간으로 픽셀별 보정 벡터와 물질 분해 벡터를 생성하는 ETB‑Cal 알고리즘을 제안한다. 실험 결과, 기존 방법 대비 비균일성이 크게 감소하였다.

상세 분석

본 연구는 PCCT 시스템에서 에너지 임계값 편향(ETB)이 스펙트럼 불일치의 주요 원인임을 명확히 규정하고, 이를 물리 기반 두 항 모델(스키우 gᵢ(E)와 편향 ΔEₖᵢ)으로 분리한다. 스키우 항은 픽셀마다 존재하는 필터 두께 편차(ΔL_b,i)를 기반으로 한 지수 함수 형태로 정의되며, 이는 CdTe와 Al 같은 기본 재료의 선형 감쇠 계수를 이용해 1차 테일러 전개로 근사한다. 편향 항은 실제 임계값이 명목값과 차이 나는 정도를 나타내며, 작은 편차 가정 하에 적분 구간을 선형화하여 ΔEₖᵢ·S₀·gᵢ·H·R 형태로 표현한다. 이러한 가정은 실험적으로 검증된 바와 같이 편향이 수 keV 이하일 때 오차가 미미함을 보여준다.

모델을 정량화하기 위해 λₖᵢ(픽셀 i, 임계값 k) 를 3개의 기본 스펙트럼 함수 F₁, F₂, F₃ 로 분해하고, 각 함수에 대한 적분값 Aₖˡ, Bₖˡ을 사전 계산한다. 이후 αᵢ, ΔDᵢ, ΔTᵢ, ΔEₖᵢ 네 개의 미지 파라미터를 다중 필터(µ₁, µ₂)와 다양한 kVp 조건에서 얻은 측정값으로 선형 연립 방정식 형태로 풀어낸다. 이 과정은 기존의 다항식 기반 보정이 요구하는 대규모 캘리브레이션 데이터와 달리, 몇 가지 알려진 필터 조합만으로 충분히 파라미터를 추정할 수 있게 한다.

알고리즘 흐름은 (1) 원시 카운트 데이터의 합-정규화, (2) 스키우 항 gᵢ(E) 계산, (3) 편향 항 ΔEₖᵢ 계산 순으로 진행되며, 각 단계는 행렬 연산 수준의 연산량으로 구현 가능해 실시간 적용이 가능하다. 실험에서는 테이블탑 PCCT 시스템에 Gammex 다중 에너지 팬텀과 교토 두부 팬텀을 사용했으며, 보정 전후의 HU 비균일성을 비교하였다. 보정 전 29.3 HU였던 Gammex 팬텀의 비균일성이 5.8 HU로, 교토 팬텀은 27.9 HU에서 3.2 HU로 각각 감소하였다. 이는 동일 조건에서 다항식 기반 모델(8.3 HU, 6.5 HU)보다 우수한 결과이며, X‑ray fluorescence 기반 캘리브레이션 없이도 높은 정확도를 달성한다는 점에서 실용성이 크게 향상된다.

또한 물질 분해 실험에서는 교토 두부 팬텀의 K‑엣지 물질(요오드) 정량이 보정 전보다 평균 12 % 감소한 반면, ETB‑Cal 적용 후 2 % 이하의 오차로 회복되었다. 이는 스키우와 편향을 독립적으로 보정함으로써 에너지 분해 정확도가 크게 개선됨을 의미한다. 마지막으로, 알고리즘의 연산 복잡도는 픽셀당 4개의 실수 파라미터만을 필요로 하며, GPU 가속 시 실시간 영상 재구성 파이프라인에 무리 없이 통합될 수 있음을 시뮬레이션으로 확인하였다.

요약하면, 본 논문은 물리적 근거에 기반한 두 항 모델을 통해 ETB를 정량적으로 계산하고, 최소한의 캘리브레이션으로 고정밀 스펙트럼 보정을 구현한다는 혁신적인 접근을 제시한다.