[論文レビュー] 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

ひとこと要約

本稿では、FFTに基づく離散サイン変換と線形補間を用いて、ラドンデータからの画像再構成のための高速なOPED（Orthogonal Polynomial Expansion on the Disk）アルゴリズムの実装を提示する。この手法は、元の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アルゴリズムの実装を開発すること。
• FFTおよび補間技術を活用することで、OPEDをFused Back Projection（FBP）と同等の速度にまで向上させること。
• 画像関数のやや滑らかな条件下で、高速OPEDアルゴリズムの収束性を証明すること。
• 標準的なアーチファクト（phantoms）を用いた数値的検証を通じて、高速OPEDの精度と性能を検証すること。

提案手法

• 高速OPEDアルゴリズムは、直交多項式展開の直接計算に代わり、FFTを用いて離散サイン変換を計算する。
• 非一様なラドンプロジェクションのサンプリングに対応するため、線形補間ステップを組み込むことで、計算を効率化する。
• この手法は、チェビシェフノードにおけるガウス求積に基づいており、再構成式は角度方向のビューと径方向の求積点の和として表現される。
• チェビシェフ多項式の第二種の性質を活用することで、FFTによる高速化を実現する。
• 高速実装は、元のOPED式における連続積分を離散求積と補間によって近似することで導出される。
• 理論的収束性は、再構成作用素の作用素ノルムを評価することで確立され、補間のもとでもO(m log m)のまま保たれることを示している。

実験結果

リサーチクエスチョン

• RQ1 OPEDアルゴリズムを、再構成精度を維持したままFBPと同等の計算速度にまで高速化できるか？
• RQ2 OPEDにFFTと線形補間を適用することで、元のアルゴリズムの収束性が保持されるか？
• RQ3 滑らかでない画像に対して、補間によって生じる誤差は、元のOPEDと比べてどの程度か？
• RQ4 高速OPEDは、元のOPEDと比較して、定量的な速度向上と誤差のトレードオフはどのようになるか？
• RQ5 連続的な2階または4階微分可能な関数に対して、高速OPEDアルゴリズムは一様収束するか？

主な発見

• 高速OPEDアルゴリズムは、元のOPEDに比べ26倍の高速化を達成しており、m=512の場合、再構成時間を344秒から13秒に短縮した。
• 元のシェップ・ローガン・アーチファクトと高速OPED再構成の間の相対二乗誤差（RSE）は0.00249574であり、元のOPEDの0.00239702と比べてわずかに高いが、差は無視できる。
• 高速OPEDの平均誤差（ME）は0.00981329であり、元のOPEDの0.0129175よりも低く、平均誤差が改善されていることを示している。
• 両再構成の差分画像は、視覚的および定量的にほとんど差がなく、差分RSEはたった0.000515499にとどまる。
• 理論的解析により、補間を伴う高速OPEDがC⁴(B)に属する関数に対してO(m log m)の収束性を維持することが確認され、元のアルゴリズムと同等の収束速度が保たれている。
• 数値的結果は、高速OPEDがシャープな不連続性を持つ画像（例えばシェップ・ローガン・アーチファクト）に対しても、元のOPEDとほぼ同等の性能を示すことを示している。