[논문 리뷰] Non-Invertibility of spectral x-ray photon counting data with pileup

Robert Alvarez|arXiv (Cornell University)|2017. 02. 16.

Advanced X-ray and CT Imaging참고 문헌 11인용 수 4

한 줄 요약

이 논문은 고역반응 및 세 에너지 범위를 가진 광자 수세기 스펙트럼 CT 데이터가 악조건이 되어, 기저 계수의 선적분 추정치에서 노이즈 분산이 10⁵ 배 증가할 수 있음을 보여준다. 크래머-라오 하한을 사용하여 이 불안정성은 특정 감쇠 수준에서 행렬의 악조건성 때문임을 밝혀내었으며, 이는 네 범위 데이터나 저역반응 조건과는 다름을 보여주며, 스펙트럼 CT 시스템의 핵심 설계 취약점을 드러낸다.

ABSTRACT

-In the Alvarez-Macovski method [R.E. Alvarez and A. Macovski, Phys. Med. Biol., 1976, 733-744], the attenuation coefficient is approximated as a linear combination of functions of energy multiplied by coefficients that depend on the material composition at points within the object. The method then computes the line integrals of the basis set coefficient from measurements with different x-ray spectra. This paper shows that the transformation from photon counting detector data with pileup to the line integrals can become ill-conditioned under some circumstances leading to highly increased noise. Methods: An idealized model that includes pileup and quantum noise is used. The noise variance of the line integral estimates is computed using the Cram{è}r-Rao lower bound (CRLB). The CRLB is computed as a function of object thickness for photon counting detector data with three and four bin pulse height analysis (PHA) and low and high pileup. Results: With four bin PHA data the transformation is well conditioned with either high or low pileup. With three bin PHA and high pileup, the transformation becomes ill-conditioned for specific values of object attenuation. At these values the CRLB variance increases by approximately 10 5 compared with the four bin PHA or low pileup results. The condition number of the forward transformation matrix also shows a spike at those attenuation values. Conclusion: Designers of systems using counting detectors should study the stability of the line integral estimator output with their data.

연구 동기 및 목표

• 광자 수세기 검출기를 사용한 스펙트럼 X선 영상에서 고역반응이 있는 경우 선적분 추정기의 통계적 안정성을 조사하기 위해.
• 측정된 광자 수에서 기저 집합 계수로의 변환 과정이 악조건이 되는 조건을 규명하기 위해.
• 크래머-라오 하한(CRLB)을 사용하여 역반응 및 에너지 범위 수가 추정기 노이즈 분산에 미치는 영향을 정량화하기 위해.
• 역반응 및 에너지 범위 분할과 같은 실제 검출기 제약 조건 하에서 추정기 안정성을 평가하는 데 시스템 설계자들을 안내하기 위해.

제안 방법

• 실제 X선 스펙트럼과 감쇠 계수를 사용하여 X선 투과 모델링을 수행하고, 비가속성 정지시간 모델을 사용해 역반응을 시뮬레이션함.
• 이deal 직사각형 에너지 범위를 가진 펄스 높이 분석(PHA)을 사용해 각 에너지 범위의 기록된 광자 수 계산.
• 양자 노이즈와 역반응 효과를 포함한 다변량 정규분포로 PHA 범위 수의 확률 분포 모델링.
• 기저 집합 계수(A-벡터)에 대한 비편향 추정기의 이론적 최소 분산을 계산하기 위해 크래머-라오 하한(CRLB) 적용.
• 역행렬 과정의 수치적 안정성을 평가하기 위해 전방 변환 행렬 M의 조건수 분석.
• 다양한 물체 두께에서 저역반응 및 고역반응 조건 하에서 세 범위 및 네 범위 PHA 시스템 간 결과 비교.

실험 결과

연구 질문

• RQ1 세 범위 광자 수세기 스펙트럼 CT에서 고역반응이 기저 계수 선적분 추정치에 심각한 불안정성을 초래하는가?
• RQ2 고역반응 조건에서 물체 두께에 따라 A-벡터 추정기의 노이즈 분산은 저역반응 또는 네 범위 시스템과 비교해 어떻게 변화하는가?
• RQ3 어떤 특정 감쇠 수준에서 전방 변환 행렬이 악조건이 되며, 그 원인은 무엇인가?
• RQ4 왜 세 범위, 고역반응 데이터에서는 크래머-라오 하한 분산이 급격히 증가하지만, 네 범위 또는 저역반응 경우는 그렇지 않은가?
• RQ5 스펙트럼 왜곡이 존재할 때 역반응 매개변수 η가 조건수와 추정기 노이즈에 얼마나 영향을 미치는가?

주요 결과

• 세 범위 PHA 및 고역반응 조건에서, 특정 물체 두께에서 크래머-라오 하한(CRLB) 분산이 네 범위 PHA 또는 저역반응 조건보다 약 10⁵ 배 증가함.
• 이 불안정성은 동일한 감쇠 수준에서 전방 변환 행렬 M의 조건수에 급격한 피크가 나타나면서 악조건이 되었음을 나타내며 관련이 있음.
• 최대 분산과 조건수 피크는 일관되게 η ≈ 0.73 ± 0.04에서 발생하여, 데이터가 노이즈에 매우 민감해지는 특정 물리적 영역임을 시사함.
• 세 범위 데이터의 저역반응 조건에서는 불안정성이 나타나지 않으며, 네 범위 데이터의 고역반응 조건은 여전히 잘 조건화되어 있어, 범위 수와 역반응 수준이 안정성에 공동으로 영향을 미침.
• PHA 범위 수 곡선은 최대 분산 지점에서 기울기가 0임을 보여주며, 이는 해당 감쇠 수준에서 정보 손실과 데이터의 역행렬 불가능성이 발생했음을 의미함.
• 물체 두께가 증가함에 따라 초기에는 역반응 감소(낮은 η)로 인해 노이즈 분산이 감소하지만, 나중에는 조건수 향상 속도보다 광자 수 감소 속도가 더 빨라져 다시 증가함.