[論文レビュー] Tomographic X-ray data of a walnut

K. Hämäläinen, Lauri Harhanen+4 more|arXiv (Cornell University)|Feb 11, 2015

Nuclear Physics and Applications被引用数 19

ひとこと要約

本論文は、高分解能のシンチグラム、測定行列、および真値FBP再構成を含む、入手可能なトモグラフィックX線データセットを提示している。このデータセットにより、スパースデータおよび正則化に基づくCT再構成アルゴリズムのベンチマーク評価が可能となる。本データは、ウェーブレットに基づくスパarsity手法、全変動正則化、および低投影条件下におけるS曲線パラメータ選択のテストを支援する。

ABSTRACT

This is the documentation of the tomographic X-ray data of a walnut made available at http://www.fips.fi/dataset.php . The data can be freely used for scientific purposes with appropriate references to the data and to this document in arXiv. The data set consists of (1) the X-ray sinogram of a single 2D slice of the walnut with three different resolutions and (2) the corresponding measurement matrices modeling the linear operation of the X-ray transform. Each of these sinograms was obtained from a measured 120-projection fan-beam sinogram by down-sampling and taking logarithms. The original (measured) sinogram is also provided in its original form and resolution. In addition, a larger set of 1200 projections of the same walnut was measured and a high-resolution filtered back-projection reconstruction was computed from this data; both the sinogram and the FBP reconstruction are included in the data set, the latter serving as a ground truth reconstruction.

研究の動機と目的

• 医療および産業イメージング分野におけるスパースデータ再構成アルゴリズムのテストに向けた現実世界のCTデータを提供すること。
• トモグラフィーにおける投影データの制限という課題に対処するため、高コントラストで構造的に非凸性を示す複雑な生物学的サンプルを提供すること。
• 低信号領域におけるTikhonovおよび全変動正則化などの正則化手法の開発と検証を支援すること。
• 1200投影からの真値再構成を提供することで、再構成精度の定量的評価を可能にすること。
• 既知の測定幾何学的構成と解像度レベルを備えた標準的かつ再現可能なデータセットを提供することで、逆問題分野における手法的研究を促進すること。

提案手法

• 80 kV、200 µA、50 µm画素分解能を備えた自作µCTシステムを用いて、くるみの120投影のファンビームCTデータを取得した。
• 2次元投影画像の中央行を抽出し、2296×120のファンビームシンチグラムを形成した。これは中心2次元断面を表す。
• 元のシンチグラムをボーリングによりダウンサンプリングし、対数変換を適用することで、3つの解像度レベル（82×82、164×164、328×328）を生成した。
• 線形X線変換をモデル化するため、サイズN²×Mのスパース測定行列Aを構築した。ここでMは投影数、Nは再構成画素数を表す。
• より高い精度を得るため、0.3°の角度ステップで1200投影の第二の高分解能データセットを取得し、フィルタドバックプロジェクション（FBP）を用いて高分解能真値再構成を計算した。
• MATLAB互換のデータファイル（シンチグラム、測定行列、FBP再構成）を提供し、アルゴリズムのテストおよび検証に直接利用可能にした。

実験結果

リサーチクエスチョン

• RQ1 スパース性を促進する再構成手法は、限られた投影数のもとで、現実的で複雑な生物学的データに対してどの程度の性能を示すのか？
• RQ2 低データ領域における正則化パラメータの選択が、再構成品質に与える影響は何か？
• RQ3 測定行列の解像度の選択が、反復的再構成アルゴリズムの精度および安定性に与える影響は何か？
• RQ4 1200投影からの真値再構成が、スパースデータ手法の評価に信頼できるベンチマークとしてどの程度有効なのか？
• RQ5 層状の殻や非凸性を示す内部構造を有するくるみの構造的複雑性は、臨床的または産業的CTの課題を効果的に模倣できるか？

主な発見

• 3つのダウンサンプリングされたシンチグラム（82×82、164×164、328×328）と、それに対応する測定行列およびシンチグラムデータを含む。
• 元の120投影シンチグラムのサイズは2296×120であり、3つの解像度レベルにおける画素サイズはそれぞれ0.513 mm、0.257 mm、0.128 mmである。
• 1200投影（0.3°角度ステップ）から得られた高分解能FBP再構成により、2296×2296の真値画像が得られた。
• 測定行列Aはスパースであり、A*x = m(:) の線形系を表しており、Tikhonovおよび全変動正則化フレームワークでの利用が可能である。
• 本データセットは、ウェーブレットに基づくスパース性手法および正則化パラメータ選択のS曲線法のテストに成功して使用された。
• 適切なクレジット表記を伴えば、科学的利用に無料で提供されており、逆問題研究における再現性とベンチマーク評価を支援する。