《MEDICAL PHYSICS》：Quantifying photon counting detector (PCD) performance using PCD-CT images. Part II. detector spectral response function of clinical scanners
编辑推荐：
专业术语翻译摘要
背景：先前发表于《医学物理》（Med. Phys.）52:2796（2025）的研究报道了一种从重建的PCD-CT图像中估计光子计数探测器（photon-counting detector, PCD）平均输出计数的方法。该能力为临床PCD-C
专业术语翻译摘要
背景：先前发表于《医学物理》（Med. Phys.）52:2796（2025）的研究报道了一种从重建的PCD-CT图像中估计光子计数探测器（photon-counting detector, PCD）平均输出计数的方法。该能力为临床PCD-CT系统的终端用户在不获取原始数据的情况下表征PCD性能提供了基础。在先前工作的基础上，本手稿提出了一种图像域方法来表征PCD的光谱响应函数（spectral response function, SRF）。SRF强烈影响PCD-CT的光谱成像性能，且即使在采用相同半导体传感器材料的情况下，SRF也可能因PCD-CT平台的不同而有所差异。
目的：利用物理模体的DICOM图像，对临床全身PCD-CT系统（西门子NAEOTOM Alpha）的理论PCD SRF模型进行校准与验证。
方法：将厚度范围为1至13英寸的PMMA板置于扫描仪中作为X射线能谱修饰器，同时管电压在70至140 kV之间变化。Alpha PCD-CT在静态机架模式下运行，使得每次扫描期间入射到探测器的X射线能谱保持静态。采集的投影数据由扫描仪重建为DICOM图像，如同其是在完整的360°旋转过程中采集的一样。对于每个PCD能量仓，根据重建图像的噪声功率谱（noise power spectrum, NPS）估计探测器的平均输出计数。随后利用各能量仓之间估计平均计数的比值，来校准基于Alpha PCD半导体材料和像素几何结构所建立的理论SRF模型的自由参数。该方法类似于利用实验测量的半值层（half-value layer, HVL）或四分之一值层来将参数化X射线能谱模型校准至放射或CT成像系统的实际能谱。随后，通过使用铜（Cu）滤过器（0.5–4.5 mm）而非PMMA板所生成的不同X射线能谱条件下获取的PCD-CT能量仓图像，比较理论预测与实验估计的能量仓间计数比值，对校准后的SRF模型进行验证。
结果：对于基于PMMA的校准数据，模型预测的总能量仓与高能量仓计数比值的平均百分比误差（mean percent error, MPE）为-0.48%（95%置信区间：[-8.00%, 7.03%]）。对于独立的基于Cu的验证数据，模型预测的计数比值与实验测量结果保持高度一致，MPE为-2.27%（95%置信区间：[-10.76%, 6.22%]）。校准后的SRF模型显示，Alpha PCD的低能量仓与高能量仓之间的平均能量间隔在70、90、120和140 kV射束下分别为10.4、16.5、24.4和28.7 keV，表明在较高管电压下具有更大的光谱分离度。
结论：PCD-CT系统的终端用户在不获取原始探测器计数的情况下，从重建图像中估计SRF是可行的。估计的Alpha PCD SRF曲线已包含在MATLAB软件包PCD-SPECTRA中，该软件包已提供给研究界，以促进PCD-CT系统光谱性能的表征，并支持进一步的研究与应用开发。
论文解读文章
一、研究背景与问题
自全身光子计数探测器CT（photon-counting detector CT, PCD-CT）进入临床应用以来，大量研究致力于评估其光谱或多能量成像能力，尤其是与采用能量积分探测器（energy-integrating detector, EID）的双能CT系统进行比较。然而，CT是一个复杂的成像系统，除探测器本身外，众多软硬件组件共同影响图像质量。新型PCD-CT平台不仅引入了全新的探测器架构，还更新了数据预处理、重建及图像后处理算法。因此，要分离PCD对成像性能的贡献，需要超越简单图像质量比较的额外测量。
表征探测器光谱响应函数（spectral response function, SRF），即给定入射光子能量下的脉冲高度谱，是理解PCD在光谱成像中具体作用的最直接且权威的方法之一。SRF类似于传统辐射探测器中量子探测效率（detective quantum efficiency, DQE）所起的作用，是光谱探测器的基础性能指标。然而，目前商用PCD-CT系统的SRF尚未由制造商公开报告。
对于终端用户而言，原则上可通过扫描能量阈值来实验测量SRF，但这需要直接访问探测器计数和自由调节能量阈值的能力，而当前临床PCD-CT系统均不具备这些能力。因此，终端用户无法在这些商用平台上直接测量和表征SRF。作为替代方案，已针对采用碲化镉（cadmium telluride, CdTe）、碲锌镉（cadmium zinc telluride, CZT）或硅（silicon, Si）传感器的PCD建立了理论和数值SRF模型。然而，无论采用何种理论方法，建立适用于特定探测器系统的精确模型都需要对PCD系统（包括传感器、电接触和读出电子学）有全面了解。不同PCD系统之间不存在普遍适用的SRF模型，且读出电子学的精确配置和参数通常是专有的，探测器开发方之外无法获取。因此，即使认为已纳入所有相关物理过程，未经实验验证的SRF模型也可能无法准确反映特定PCD系统的真实光谱响应。
二、研究内容与结论
针对上述问题，研究人员开展了一项研究，旨在为临床PCD-CT系统提供一种基于DICOM图像校准和验证SRF模型的实用途径。该研究并非旨在提出新的SRF模型或展示其优于现有模型的优势，其新颖性在于提供了一条利用DICOM图像适配和验证特定SRF模型于临床PCD-CT系统的可行路径。
研究人员采用了一个先前为CdTe基PCD开发的参数化SRF模型，该模型基于并联和串联级联系统分析建立，涵盖了K荧光逃逸与再吸收、瑞利和康普顿散射、电荷收集、电荷扩散、电子噪声注入及能量阈值化等主要物理过程。模型中的部分参数由CdTe材料固有物理性质决定，部分可从厂商规格获得，而剩余参数集则取决于传感器制造质量、半导体-金属接触及专用集成电路（application-specific integrated circuit, ASIC），需针对每台PCD进行校准。
研究人员利用西门子NAEOTOM Alpha PCD-CT在静态机架模式下采集数据，通过改变PMMA板厚度（0至13英寸）和管电压（70、90、120和140 kV）生成56种不同的入射X射线能谱条件。对于每个PCD能量仓，根据重建图像的噪声功率谱（noise power spectrum, NPS）估计平均探测器输出计数，并计算总能量仓（total-energy, TE）与高能量仓（high-energy, HE）之间的计数比值。通过优化算法，将模型预测的TE/HE计数比值与实验测量值进行匹配，从而校准SRF模型的三个自由参数（电荷扩散半径、电荷收集因子和电子噪声）以及四个管电压设置下未知的HE能量阈值。随后，使用铜（Cu）滤过器（0.5–4.5 mm）替代PMMA板生成显著不同的X射线能谱，对校准后的SRF模型进行独立验证。
研究结果表明，对于PMMA校准数据，模型预测的TE/HE计数比值与测量值高度吻合，相对均方根误差（relative root-mean-square error, rRMSE）为0.51%，平均百分比误差（mean percent error, MPE）为-0.48%。对于独立的Cu验证数据，模型预测与实验测量结果同样保持良好一致，rRMSE为1.13%，MPE为-2.27%。校准后的SRF模型显示，Alpha PCD低能量仓与高能量仓之间的平均能量间隔在70、90、120和140 kV射束下分别为10.4、16.5、24.4和28.7 keV，表明较高管电压下具有更大的光谱分离度。
三、关键技术方法
研究人员采用了以下主要关键技术方法：1）基于级联系统分析的理论SRF模型，该模型解析计算光子-半导体相互作用概率、沉积能量的均值与方差以及入射或次级光子的平均路径长度，而非像蒙特卡洛模拟那样显式采样单个相互作用事件；2）基于图像域NPS的探测器计数估计方法，通过静态机架模式下采集的成对图像相减去除背景趋势，从噪声图像中计算二维NPS，进而估计每个PCD能量仓的平均探测器计数；3）参数化X射线能谱建模，采用TASMIP方法生成基础能谱，并通过迭代调整等效固有滤过使模型预测的HVL与实验测量值匹配；4）嵌套优化算法，外层采用Nelder-Mead单纯形算法优化SRF模型参数，内层通过穷举搜索在40至80 keV范围内以1 keV间隔确定各管电压下的最优HE能量阈值。
四、研究结果
模型校准结果：利用实验校准数据，确定的最优SRF自由参数为：电荷扩散半径16.4 μm，电荷收集因子0.98，电子噪声1.9 keV。70、90、120和140 kV设置下的HE能量阈值分别估计为53、64、68和68 keV。模型预测的TE/HE计数比值与校准数据集的测量值紧密匹配，rRMSE为0.51%，MPE为-0.48%。此外，将CdTe法诺因子（Fano factor）从0.1分别调整为0和0.2进行敏感性分析，结果显示优化参数和最终SRF均仅有微小变化，表明该校准方法在考察范围内对法诺因子的假设不敏感。
模型验证结果：将校准后的SRF模型应用于采用Cu滤过器获得的独立验证数据集，模型预测的TE/HE计数比值与测量值紧密匹配，rRMSE为1.13%，MPE为-2.27%。对于给定管电压，TE/HE比值随滤过增加而单调下降，这与射束硬化效应一致。校准后的SRF以二维热图形式呈现，可清晰显示光电峰和由Cd及Te K荧光逃逸产生的次级能带。四个代表性单能入射X射线（40、60、80和100 keV）的光电峰半高全宽（full width at half maximum, FWHM）分别为4.8、4.8、4.8和4.7 keV，几乎一致，因其主要由与入射能量无关的电子噪声参数决定。
输入能谱数量对SRF估计的影响：将校准数据集从56个条件逐步缩减至8个条件，估计的SRF模型参数有所变化，但验证精度保持稳定，rRMSE仅增加不到0.01%。仅在缩减至8个条件或更少时观察到rRMSE的轻微增加。这表明所提出的校准框架对校准条件数量具有鲁棒性，远少于56个X射线能谱配置可能即可达到相当的估计精度。
统计不确定性对SRF估计的影响：通过自举分析将图像切片分为五组独立进行SRF校准，各条件下TE/HE计数比值的平均相对标准差为2.8%。虽然各组间估计的模型参数存在明显变化，但所得SRF的验证性能保持相对稳定，rRMSE范围为1.13%至1.30%。这表明图像导出的计数比值不确定性可导致个别拟合参数的变化，而所得SRF的预测精度对此相对不敏感，同时支持使用更多图像切片进行NPS估计以降低统计不确定性。
PCD SRF模型的应用：利用校准后的SRF，研究人员计算了四个管电压设置下、以20 cm软组织等效材料作为射束衰减体时的低能量仓（low-energy, LE）和高能量仓（high-energy, HE）有效X射线能谱。LE与HE仓之间的平均能量间隔在70、90、120和140 kV下分别为10.4、16.5、24.3和28.3 keV，表明较高管电压下具有更大的光谱分离度，可直接可视化和量化两仓间的光谱重叠程度。
五、讨论与结论
该研究提供了一种从重建DICOM图像估计PCD SRF的实用方法，使PCD-CT系统的终端用户无需访问原始探测器数据即可表征探测器光谱性能。校准后的SRF模型在独立验证条件下与实验测量保持良好一致，证明了该方法的可行性和准确性。该方法类似于利用实验测量的半值层校准参数化X射线能谱模型的经典做法，并通过使用多种入射能谱和阈值条件显著降低了参数解非唯一性的风险。研究提供的MATLAB软件包PCD-SPECTRA可供研究界使用，以促进PCD-CT系统光谱性能的表征并支持进一步的研究与应用开发。
