CT の逐次近似再構成とは?
逐次近似再構成(反復再構成)は CT ボリュームを段階的に計算します。現在の推定から投影をシミュレーションし、測定した投影と比較して、両者が一致するまでボリュームを修正します。Voxray QRT のようなモデルベース手法はスキャンの物理もシミュレーションするため、ビームハードニングや金属アーティファクトがそもそも生じません。
逐次近似再構成の仕組み
多くの CT 装置はフィルタ補正逆投影(FDK)で再構成します。フィルタ処理した投影から 1 回の処理でボリュームを計算するため高速ですが、理想的なデータを前提としています:単一の X 線エネルギー、正確なジオメトリ、円軌道、そしてノイズのない十分な数の投影です。
逐次近似再構成は、再構成を最適化問題として扱います。ボリュームの推定から始め、そのボリュームが生じる投影をシミュレーションし、測定した投影と比較して、差が小さくなるようにボリュームを更新します。これを何度か繰り返すと、ボリュームは測定データを可能な限りよく説明するようになります。
シミュレーションの段階で測定を必要なだけ詳細にモデル化できるため、逐次近似法は FDK が無視せざるを得ないものを考慮できます:各検出器画素のノイズ、X 線スペクトル、線源と検出器の実際の位置、そして部品に関する事前知識です。
逐次近似再構成が役立つ場面
投影数が少ない、または角度が限られる場合。スキャン時間を短くしたいときや、部品をすべての方向から撮影できないとき。
短い露光、高密度の材料、小さなボクセルによるノイズの多いデータ。
ビームハードニングや金属アーティファクトを伴う複数材料の部品(手法が X 線スペクトルをモデル化する場合)。
ヘリカル、ラミノグラフィ、ロボット CT などの円軌道以外の軌道、およびジオメトリ誤差のあるスキャン。
従来の逐次近似再構成と QRT の比較
オイルセパレーターのスキャンからの 2 つのビュー。スライダーをドラッグして比較してください。
逐次近似再構成の手法の種類
手法の違いは、測定のシミュレーションに何を含めるかにあります。モデル化する範囲が広いほど、避けられるアーティファクトも増えます。
代数的手法(ART、SART、SIRT)
仕組み: ボリュームの投影が測定した投影に一致するよう、光線ごと、またはすべての光線について一度に、ボリュームを繰り返し修正します。
限界: FDK と同じ簡略化した物理を使うため、ビームハードニングは残ります。反復を重ねるとノイズが増えるため、結果はどこで止めるかに左右されます。
統計的手法
仕組み: 検出器の光子統計をモデル化し、高密度の材料を通過する光線など、ノイズの多い光線の重みを小さくします。
限界: ノイズは減りますが、多くの場合は依然として単一の X 線エネルギーを前提とします。ノイズモデルの強さは調整が必要です。
正則化手法(TV、圧縮センシング)
仕組み: 部品が一定密度の領域から成るといった事前知識を加えます。これにより少ない投影からの再構成が可能になります。
限界: 正則化が強すぎると細部が失われ、なだらかな変化が平らでブロック状の領域になります。
モデルベース手法
仕組み: 多色スペクトル、検出器、正確なジオメトリなど、スキャンの物理をシミュレーションします。これらの効果によるアーティファクトはそもそも生じません。
限界: 反復ごとの計算量がはるかに多く、装置の正確なモデルが必要です。GPU がなければ日常的に使うには遅すぎます。
Voxray QRT における逐次近似再構成
QRT は産業用 CT のための物理ベースのモデルベース逐次近似再構成です。多色 X 線スペクトルと測定の物理をモデル化し、記録された投影を最もよく説明するボリュームを計算します。
QRT は再構成中に各投影のジオメトリ誤差も推定するため、ロボット CT、ラミノグラフィ、その他の軌道を同じアルゴリズムで扱えます。NVIDIA GPU 上で動作し、データセットごとにパラメータを調整しなくても良好な結果が得られます。
物理モデル
ビームハードニングや金属アーティファクトを後から補正するのではなく、発生自体を防ぎます。
ジオメトリ補正
機械的な公差、軸のずれ、焦点のドリフトを投影ごとに補正します。
あらゆる軌道
円軌道、ヘリカル、ラミノグラフィ、ロボット CT を、お使いのスキャナーの生データから再構成します。
GPU による高速処理
2k 検出器の 2,500 投影は、RTX 4090 1 枚で 32 分、2 枚で 18 分です。
ケーススタディ
7 ピン電源コネクタの欠陥検出
Robert Bosch GmbH Corporate Research と共同で、銅合金ピンを樹脂ハウジングに収めたコネクタを対象に、フィルタ補正逆投影、デュアルエナジー融合、QRT を比較し、研磨断面で検証しました。ピン間のボイドは、フィルタ補正逆投影では見えず、デュアルエナジーでは一部のみ見え、QRT ではっきりと分解されました。
事例を読む
逐次近似再構成に関するよくある質問
CT の逐次近似再構成とは何ですか?
現在の推定から投影を繰り返しシミュレーションし、測定した投影と比較してボリュームを修正することで、シミュレーションと測定データが一致するまで CT ボリュームを計算する手法です。
フィルタ補正逆投影と逐次近似再構成の違いは何ですか?
フィルタ補正逆投影(FDK)は 1 回の処理でボリュームを計算し、理想的なデータを前提とします。逐次近似再構成は何度も処理を繰り返しますが、ノイズ、X 線スペクトル、実際のスキャンジオメトリをモデル化できるため、アーティファクトが減り、少ない投影やあらゆる軌道に対応できます。
逐次近似再構成は FDK より遅いですか?
はい。反復ごとにすべての投影をシミュレーションして比較します。GPU を使えば日常的に使える速さで、当社のベンチマークでは QRT は 2k 検出器の 2,500 投影を RTX 4090 1 枚で 32 分で再構成します。
逐次近似再構成自体がアーティファクトを生むことはありますか?
はい。単純な代数的手法では反復しすぎるとノイズが増幅され、正則化が強すぎると表面がブロック状に見えて細部が失われます。モデルベース手法は、データを平滑化するのではなく物理で説明するため、こうした調整の必要が少なくなります。
産業用 CT の逐次近似再構成は医療用 CT と同じですか?
原理は同じです。医療用 CT では主に被ばく線量の低減に使われます。産業用 CT では、複数材料部品のアーティファクト低減、信頼できる測定、スキャン時間の短縮、円軌道以外の軌道が目的になります。
自社のスキャンで QRT を試せますか?
はい。無料トライアルを申し込んで、ご自身のワークステーションでデータセットを再構成するか、再構成サービスを選んでデータセットをお送りいただければ、当社で再構成します。
ご自身のスキャンで逐次近似再構成を試す
無料トライアルでご自身のデータを QRT で再構成するか、部品とスキャナーについてご相談ください。