[论文解读] Non-Invertibility of spectral x-ray photon counting data with pileup

Robert Alvarez|arXiv (Cornell University)|Feb 16, 2017

Advanced X-ray and CT Imaging参考文献 11被引用 4

一句话总结

本文表明，具有高堆积效应和三个能量通道的光子计数能谱CT数据可能变得病态，导致基系数线积分估计的噪声方差增加10⁵倍。通过Cramér-Rao下界分析，发现这种不稳定性源于特定衰减水平下的矩阵病态性，而非常见的四通道数据或低堆积条件，揭示了能谱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.

研究动机与目标

• 研究使用具有堆积效应的光子计数探测器进行光子计数能谱CT时，线积分估计器的统计稳定性。
• 识别从测量光子计数到基集系数转换过程变得病态的条件。
• 利用Cramér-Rao下界（CRLB）量化堆积效应和能量通道数量对估计器噪声方差的影响。
• 为系统设计者提供指导，以评估在实际探测器限制（如堆积效应和能量分组）下的估计器稳定性。

提出的方法

• 使用真实X射线能谱和衰减系数建模X射线传输，通过非可停死时间模型模拟堆积效应。
• 利用脉冲幅度分析（PHA）计算每个能量通道的记录光子计数，采用理想矩形能量通道。
• 将PHA通道计数的概率分布建模为多元正态分布，同时包含量子噪声和堆积效应。
• 应用Cramér-Rao下界（CRLB）计算A向量（基集系数）无偏估计器的理论最小方差。
• 分析前向变换矩阵M的条件数，以评估逆运算过程的数值稳定性。
• 比较三通道与四通道PHA系统在不同物体厚度下，低堆积与高堆积条件下的结果。

实验结果

研究问题

• RQ1 在三通道光子计数能谱CT中，高堆积效应是否会导致基系数线积分估计出现显著不稳定性？
• RQ2 在高堆积条件下，与低堆积或四通道系统相比，A向量估计器的噪声方差如何随物体厚度变化？
• RQ3 在哪些特定衰减水平下，变换矩阵会变得病态？其不稳定性由何原因引起？
• RQ4 为何Cramér-Rao下界方差在三通道、高堆积数据中出现尖峰，而在四通道或低堆积情况下不会？
• RQ5 堆积参数η在存在能谱失真时，对条件数和估计器噪声的影响程度如何？

主要发现

• 在三通道PHA和高堆积条件下，Cramér-Rao下界（CRLB）方差在特定物体厚度处相比四通道PHA或低堆积条件增加了约10⁵倍。
• 该不稳定性与前向变换矩阵M在相同衰减水平下的条件数急剧上升有关，表明存在病态性。
• 峰值方差与条件数出现在一致的堆积参数η ≈ 0.73 ± 0.04处，表明在此特定物理条件下，数据对噪声极度敏感。
• 三通道低堆积数据不表现出该不稳定性，而四通道高堆积数据仍保持良好条件性，表明通道数与堆积水平共同决定系统稳定性。
• PHA通道计数曲线在峰值方差点处导数为零，表明在此衰减水平下信息丢失且数据不可逆。
• 随着物体厚度增加，噪声方差最初因堆积减少（η降低）而下降，但随后因光子计数下降速度超过条件数改善速度而再次上升。