[论文解读] A statistical reconstruction algorithm for positronium lifetime imaging using time-of-flight positron emission tomography

Hsin‐Hsiung Huang, Zheyuan Zhu+4 more|arXiv (Cornell University)|Jun 13, 2022

Medical Imaging Techniques and Applications被引用 5

一句话总结

本文提出了一种用于时间飞行（TOF）正电子发射断层扫描（PET）中正电子素寿命成像（PLI）的统计重建算法，通过联合建模摄取与寿命分布的似然函数，实现高保真、定量的寿命图像恢复。在采用L-BFGS-B优化算法并结合真实TOF时间分辨率模型的最大似然框架下，该方法即使在有限时间分辨率条件下，也能从模拟的TOF-PET数据中成功重建出高保真度的寿命图像。

ABSTRACT

Positron emission tomography (PET) is an important modality for diagnosing diseases such as cancer and Alzheimer's disease, capable of revealing the uptake of radiolabeled molecules that target specific pathological markers of the diseases. Recently, positronium lifetime imaging (PLI) that adds to traditional PET the ability to explore properties of the tissue microenvironment beyond tracer uptake has been demonstrated with time-of-flight (TOF) PET and the use of non-pure positron emitters. However, achieving accurate reconstruction of lifetime images from data acquired by systems having a finite TOF resolution still presents a challenge. This paper focuses on the two-dimensional PLI, introducing a maximum likelihood estimation (MLE) method that employs an exponentially modified Gaussian (EMG) probability distribution that describes the positronium lifetime data produced by TOF PET. We evaluate the performance of our EMG-based MLE method against \st{traditional} approaches using exponential likelihood functions and penalized surrogate methods. Results from computer-simulated data reveal that the proposed EMG-MLE method can yield quantitatively accurate lifetime images. We also demonstrate that the proposed MLE formulation can be extended to handle PLI data containing multiple positron populations.

研究动机与目标

• 为实际TOF-PET系统中存在有限时间分辨率的正电子素寿命成像（PLI）缺乏完整统计模型的问题提供解决方案。
• 开发一种联合最大似然重建框架，能够从列表模式TOF-PET数据中同时估计摄取图像与寿命图像。
• 将有限时间分辨率和探测器响应的影响纳入统计模型，以避免重建过程中的信息损失。
• 通过准确恢复独立于示踪剂摄取机制的寿命分布，实现定量PLI。
• 利用真实系统参数下的计算机模拟TOF-PET数据验证算法性能。

提出的方法

• 构建一个考虑有限时间分辨率及正交正电子素（o-Ps）寿命指数衰减的二维列表模式TOF-PET数据统计模型。
• 推导出同时包含摄取与寿命参数的联合对数似然函数，整合探测器响应函数与时间飞行不确定性。
• 采用L-BFGS-B优化算法，在摄取与寿命图像均满足非负性约束的条件下最大化似然函数。
• 应用惩罚型代理（PS）方法以处理逆拉普拉斯变换问题，实现寿命组分的稳定分解。
• 推导出似然函数关于寿命参数的解析梯度，以促进高效优化。
• 通过在真实TOF时间分辨率（200–600 ps FWHM）下模拟数据，并分别使用触发计数（PC）与预设时间（PT）采集模式进行图像重建，验证模型。

实验结果

研究问题

• RQ1 能否构建一个准确考虑TOF-PET系统中有限时间分辨率的统计模型，以实现可靠的正电子素寿命重建？
• RQ2 与分别重建方法相比，摄取与寿命图像的联合估计在提升PLI准确度与定量性方面有何改进？
• RQ3 所提出的最大似然算法在具有真实探测器分辨率的模拟TOF-PET数据中，能在多大程度上恢复真实的寿命分布？
• RQ4 在考虑时间分辨率影响的前提下，该算法在不同采集模式（PC与PT）下是否保持一致的性能表现？
• RQ5 该方法能否在不依赖理想化假设（如完美时间分辨率）的前提下，准确分离寿命组分？

主要发现

• 所提出的算法即使在有限时间分辨率（200–600 ps FWHM）条件下，也能从模拟的TOF-PET数据中高精度地重建出寿命图像。
• 联合最大似然估计框架能够同时恢复摄取与寿命分布，显著提升重建保真度。
• 该算法在不同采集模式下均表现出鲁棒性，PC与PT模式下的最大似然估计结果仅相差一个正标量因子。
• 在统计模型中引入有限时间分辨率可防止信息损失，并实现准确的寿命组分分离。
• 采用L-BFGS-B优化并施加非负性约束，确保了稳定收敛与物理上合理的图像重建。
• 针对寿命参数推导出的梯度可实现高效且精确的优化，证明了该方法在实际成像场景中的可行性。