[论文解读] Sinogram constrained TV-minimization for metal artifact reduction in CT

Clemens Schiffer, Kristian Bredies|arXiv (Cornell University)|Apr 26, 2014

Advanced X-ray and CT Imaging参考文献 2被引用 3

一句话总结

该论文提出了一种新颖的、基于 sinogram 约束的总变差（TV）最小化方法，通过在受损 sinogram 数据上施加不等式约束，同时利用 TV 正则化进行图像重建，从而减少 CT 成像中的金属伪影。采用 Chambolle-Pock 算法，该方法在合成 Shepp-Logan 模型中成功保留了金属植入物附近的解剖细节，硬约束下 PSNR 达到 47.6 dB，噪声数据下达到 40.1 dB。

ABSTRACT

A new method for reducing metal artifacts in X-ray computed tomography (CT) images is presented. It bases on the solution of a convex optimization problem with inequality constraints on the sinogram, and total variation regularization for the reconstructed image. The Chambolle-Pock algorithm is used to numerically solve the discretized version of the optimization problem. As proof of concept we present and discuss numerical results for synthetic data.

研究动机与目标

• 解决由高密度植入物引起的束硬化和截断效应导致的 CT 成像中金属伪影问题。
• 在金属物体导致的 sinogram 数据缺失或受损情况下，提升图像重建质量。
• 开发一种变分方法，结合总变差正则化与 sinogram 层级约束，以保留金属附近结构的细节。
• 通过合成数据和迭代优化方法，验证该方法的可行性。

提出的方法

• 该方法将凸优化问题公式化，包含未受损 sinogram 区域的数据保真项和图像上的总变差（TV）正则化项。
• 在 sinogram 上引入逐点不等式约束，确保金属影响区域的重建投影值不低于阈值 C，代表最小可检测强度。
• 使用 Chambolle-Pock 原始-对偶算法求解优化问题，该算法能高效处理非光滑 TV 和指示函数。
• 算法采用离散化的 Radon 变换及其伴随变换，并引入预条件因子以确保收敛性。
• 通过投影步骤强制执行约束集，将金属影响区域的 sinogram 值设为至少 C。
• 采用对偶公式与近端算子实现高效计算，包括 TV 的软阈值化和 sinogram 约束的截断操作。

实验结果

研究问题

• RQ1 与标准 TV 最小化相比，对 sinogram 施加不等式约束是否能提升金属影响区域 CT 扫描的图像重建质量？
• RQ2 在金属植入物附近，sinogram 层级约束的引入如何影响细微解剖结构的保留？
• RQ3 不同数据保真度和正则化参数对伪影抑制和图像质量的影响如何？
• RQ4 该方法在处理含噪声或截断的 sinogram 数据时，是否能有效保持诊断级图像质量？

主要发现

• 在含 10×10 金属插入物的合成 Shepp-Logan 幻象中，该方法在硬约束下实现了 47.6 dB 的 PSNR，显著优于从截断 sinogram 重建的标准方法。
• 在 5% 高斯噪声存在下，经过 80,000 次迭代后，该方法实现了 40.1 dB 的 PSNR，表现出对噪声的鲁棒性。
• 该方法成功保留了金属植入物附近被标准重建完全遮蔽的小型组织结构。
• Chambolle-Pock 算法实现了对非光滑、约束优化问题的稳定且收敛的求解。
• 尽管图像质量优异，该方法在 GPU 上运行 80,000 次迭代仍需约一小时计算时间，表明仍需进一步加速。
• 结果表明，该方法可推广至其他形式的投影数据损坏或缺失问题，而不仅限于金属伪影。