[论文解读] Fast OPED algorithm for reconstruction of images from Radon data

Yuan Xu, Oleg Tischenko|ArXiv.org|Mar 21, 2007

Medical Imaging Techniques and Applications参考文献 8被引用 7

一句话总结

本文提出了一种OPED（Disk上的正交多项式展开）算法的快速实现方法，通过基于FFT的离散正弦变换和线性插值，从Radon数据中重建图像。该方法在保持重建精度的同时，相比原始OPED实现了26倍的加速，数值测试显示误差差异可忽略不计，且理论证明了对适度光滑函数的收敛性。

ABSTRACT

A fast implementation of the OPED algorithm, a reconstruction algorithm for Radon data introduced recently, is proposed and tested. The new implementation uses FFT for discrete sine transform and an interpolation step. The convergence of the fast implementation is proved under the condition that the function is mildly smooth. The numerical test shows that the accuracy of the OPED algorithm changes little when the fast implementation is used.

研究动机与目标

• 开发一种计算高效的OPED算法实现，用于从Radon数据中进行图像重建。
• 通过利用FFT和插值技术，使OPED在速度上能够与FBP相媲美。
• 在图像函数具有适度光滑性条件下，证明快速OPED算法的收敛性。
• 通过在标准体模上的数值测试，验证快速OPED的精度和性能。

提出的方法

• 快速OPED算法使用FFT计算离散正弦变换，替代原始正交多项式展开的直接计算。
• 引入线性插值步骤以处理Radon投影的非均匀采样，从而实现高效计算。
• 该方法基于切比雪夫节点上的高斯求积，重建公式表示为角度视图和径向求积点的求和。
• 利用第二类切比雪夫多项式的性质，实现FFT加速。
• 通过将原始OPED公式中的连续积分近似为离散求积和插值，推导出快速实现。
• 通过有界重构算子的算子范数，建立理论收敛性，表明在插值下其仍保持O(m log m)的收敛阶。

实验结果

研究问题

• RQ1 能否在保持重建精度的同时，将OPED算法加速至与FBP相当的计算速度？
• RQ2 在OPED中使用FFT和线性插值是否保持了原始算法的收敛性质？
• RQ3 插值引入的误差与原始OPED相比，在光滑和非光滑图像上的表现如何？
• RQ4 与原始OPED相比，快速OPED在速度和误差之间的定量权衡如何？
• RQ5 对于具有连续二阶或四阶导数的函数，快速OPED算法能否实现一致收敛？

主要发现

• 快速OPED算法相比原始OPED实现了26倍的加速，当m=512时，重建时间从344秒减少至13秒。
• 原始Shepp-Logan体模与快速OPED重建结果之间的相对平方误差（RSE）为0.00249574，仅略高于原始OPED的0.00239702。
• 快速OPED的平均误差（ME）为0.00981329，低于原始OPED的0.0129175，表明平均误差有所改善。
• 两重建结果之间的差分图像在视觉和定量上均显示可忽略的差异，差分RSE仅为0.000515499。
• 理论分析证实，带插值的快速OPED在C⁴(B)类函数上仍保持O(m log m)的收敛阶，与原始算法的收敛速率一致。
• 数值结果表明，即使对于具有锐利不连续性的图像（如Shepp-Logan体模），快速OPED的表现也几乎与原始OPED相当。