Preconditioned alternating projection algorithm for solving the penalized-likelihood SPECT reconstruction problem.
Preconditioned alternating projection algorithm for solving the penalized-likelihood SPECT reconstruction problem.
复制标题
求解惩罚似然 SPECT 重建问题的预处理交替投影算法
DOI：
10.1016/j.ejmp.2017.05.001
复制
发表时间：
2017-06
期刊：
影响因子：
--
通讯作者：
Xu Y
中科院分区：
文献类型：
--
作者：
Li S;Zhang J;Krol A;Schmidtlein CR;Feiglin D;Xu Y
作者最近开发了一种预处理交替投影算法（PAPA）解决惩罚似然SPECT重建问题。该算法可以解决各种各样的不可微优化模型。本文研究了基于全变差正则化的PAPA算法（TV-PAPA）和基于嵌套期望最大化（EM）-TV迭代的前向-后向算法（FB-EM-TV）的性能。采用蒙特卡罗技术模拟了扇形束准直SPECT数据的多种噪声实现，分别为具有暖背景的分段恒定体模和具有均匀活动的热球和冷球，在两个噪声水平下。使用上述算法进行重建，并进行衰减、散射、距离相关准直器模糊和灵敏度校正。评价了图像的噪声抑制性能、病灶检测能力、病灶对比度、对比度恢复系数、收敛速度和最佳参数的选择。以传统的带TV后置滤波器的EM算法（TVPF-EM）和带高斯后置滤波器的EM算法（GPF-EM）为基准。TV-PAPA和FB-EM-TV在所有调查类别中表现出相似的性能。两种算法在图像噪声抑制、病变可检测性、病变对比度和收敛速度方面均优于TVPF-EM。我们建立了最佳参数与信息密度近似遵循幂律，这为重建方法的参数选择提供了指导。对于模拟SPECT数据，TV-PAPA和FB-EM-TV产生定性和定量相似的图像。它们的表现优于基准TVPF-EM和GPF-EM，病变对比度损失有限。
The authors recently developed a preconditioned alternating projection algorithm (PAPA) for solving the penalized-likelihood SPECT reconstruction problem. The proposed algorithm can solve a wide variety of non-differentiable optimization models. This work is dedicated to comparing the performance of PAPA with total variation (TV) regularization (TV-PAPA) and a novel forward-backward algorithm with nested expectation maximization (EM)-TV iteration scheme (FB-EM-TV). Monte Carlo technique was used to simulate multiple noise realizations of the fan-beam collimated SPECT data for a piecewise constant phantom with warm background, and hot and cold spheres with uniform activities at two noise levels. They were reconstructed using the aforementioned algorithms with attenuation, scatter, distance-dependent collimator blurring and sensitivity corrections. Noise suppressing performance, lesion detectability, lesion contrast, contrast recovery coefficient, convergence speed and selection of optimal parameters were evaluated. The conventional EM algorithms with TV post-filter (TVPF-EM) and Gaussian post-filter (GPF-EM) were used as benchmarks. The TV-PAPA and FB-EM-TV demonstrated similar performance in all investigated categories. Both algorithms outperformed TVPF-EM in terms of image noise suppression, lesion detectability, lesion contrast and convergence speed. We established that the optimal parameters versus information density approximately followed power laws, which offers a guidance in parameter selection for reconstruction methods. For the simulated SPECT data, TV-PAPA and FB-EM-TV produced qualitatively and quantitatively similar images. They performed better than the benchmark TVPF-EM and GPF-EM, with only limited loss of lesion contrast.
登录
查看更多内容
用于最大后验 ECT 重建的预处理交替投影算法
影响因子：
2.1
作者：
Krol A;Li S;Shen L;Xu Y
通讯作者：
Xu Y
DOI：
10.1073/pnas.90.21.9758
发表时间：
1993-11-01
影响因子：
11.1
作者：
BARRETT, HH;YAO, J;MYERS, KJ
通讯作者：
MYERS, KJ
影响因子：
10.6
作者：
Frey, EC;Gilland, KL;Tsui, BMW
通讯作者：
Tsui, BMW
影响因子：
10.6
作者：
Shepp, L A;Vardi, Y
通讯作者：
Vardi, Y
影响因子：
4
作者：
RUDIN, LI;OSHER, S;FATEMI, E
通讯作者：
FATEMI, E