CT 手法

CT の逐次近似再構成とは?

逐次近似再構成(反復再構成)は CT ボリュームを段階的に計算します。現在の推定から投影をシミュレーションし、測定した投影と比較して、両者が一致するまでボリュームを修正します。Voxray QRT のようなモデルベース手法はスキャンの物理もシミュレーションするため、ビームハードニングや金属アーティファクトがそもそも生じません。

逐次近似再構成の仕組み

多くの CT 装置はフィルタ補正逆投影(FDK)で再構成します。フィルタ処理した投影から 1 回の処理でボリュームを計算するため高速ですが、理想的なデータを前提としています:単一の X 線エネルギー、正確なジオメトリ、円軌道、そしてノイズのない十分な数の投影です。

逐次近似再構成は、再構成を最適化問題として扱います。ボリュームの推定から始め、そのボリュームが生じる投影をシミュレーションし、測定した投影と比較して、差が小さくなるようにボリュームを更新します。これを何度か繰り返すと、ボリュームは測定データを可能な限りよく説明するようになります。

シミュレーションの段階で測定を必要なだけ詳細にモデル化できるため、逐次近似法は FDK が無視せざるを得ないものを考慮できます:各検出器画素のノイズ、X 線スペクトル、線源と検出器の実際の位置、そして部品に関する事前知識です。

逐次近似再構成が役立つ場面

  • 投影数が少ない、または角度が限られる場合。スキャン時間を短くしたいときや、部品をすべての方向から撮影できないとき。

  • 短い露光、高密度の材料、小さなボクセルによるノイズの多いデータ。

  • ビームハードニングや金属アーティファクトを伴う複数材料の部品(手法が X 線スペクトルをモデル化する場合)。

  • ヘリカル、ラミノグラフィ、ロボット CT などの円軌道以外の軌道、およびジオメトリ誤差のあるスキャン。

CT 再構成手法の画質と計算時間フィルタ補正逆投影(FDK)は最も高速ですが、画質は最も低くなります。SIRT や ART などの従来の逐次近似法は時間がかかります。Voxray QRT も FDK より時間がかかりますが、最も高い画質が得られます。画質計算時間FDKSIRT / ARTVoxray QRT
模式図です。逐次近似法は FDK より時間がかかります。どれだけ長くかかるかは、スキャンとそのパラメータによって異なります。

従来の逐次近似再構成と QRT の比較

オイルセパレーターのスキャンからの 2 つのビュー。スライダーをドラッグして比較してください。

従来の逐次近似再構成によるオイルセパレーター電子部の断面
従来の逐次近似法
Voxray QRT で再構成した同じ断面
Voxray QRT
オイルセパレーターの電子部の断面。従来の逐次近似再構成では、高密度部分から基板と背景にストリークが広がります。QRT では背景がきれいなまま、小さな部品もシャープに見えます。
従来の逐次近似再構成によるオイルセパレーターの 3D 表示
従来の逐次近似法
Voxray QRT による同じ 3D 表示
Voxray QRT
同じ部品の 3D 表示。従来の逐次近似再構成ではコイルとリングの表面が粗く一部が欠けていますが、QRT では閉じた滑らかな表面になります。

逐次近似再構成の手法の種類

手法の違いは、測定のシミュレーションに何を含めるかにあります。モデル化する範囲が広いほど、避けられるアーティファクトも増えます。

代数的手法(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 ピン電源コネクタの 3D 表示

ケーススタディ

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 で再構成するか、部品とスキャナーについてご相談ください。

お問い合わせ