JP2024523733A - 動きアーティファクト補償による改善された金属アーティファクト低減 - Google Patents

Abstract

本発明は、歯科分野におけるデジタルボリュームトモグラフィ（ＤＶＴ）イメージングを再構成するための方法であって、方法が、（Ｓ１）動きアーティファクト補償（ＭＡＣ）を使用して、患者イメージングのサイノグラムから投影ジオメトリを推定するステップと、（Ｓ２）推定された投影ジオメトリを使用して、第１のボリュームを再構成するステップと、（Ｓ３）第１のボリュームと推定された投影ジオメトリとを使用して、サイノグラム中の金属領域と第１のボリューム（１）中の金属領域とを検出するステップと、（Ｓ４）サイノグラム中の金属領域を補正し、補正されたサイノグラムを生成するステップと、（Ｓ５）推定された投影ジオメトリと補正されたサイノグラムとを用いて、第２のボリュームを再構成するステップと、（Ｓ６）第２のボリューム中の金属領域を補正するステップとを備えることを特徴とする、方法に関する。

Description

本発明は、歯科分野におけるデジタルボリュームトモグラフィ（ＤＶＴ：digital volume tomography）イメージングを再構成するための方法に関する。詳細には、本発明は、ＤＶＴイメージングにおける金属アーティファクトの低減に関する。

ＤＶＴイメージングでは、イメージング構成要素のＸ線源とＸ線検出器とが、互いに対向し、患者の周りで回転する。サイノグラムを形成する、一連のＸ線投影画像が生成される。投影ジオメトリを知ることにより、サイノグラムからボリュームが再構成される。投影ジオメトリは、ＤＶＴデバイスの幾何学的特性と曝露中のＤＶＴデバイスの軌道とを記述する。投影ジオメトリは、投影行列によって表され得る。

ＤＶＴイメージングでは、金属などのＸ線不透過構造が、再構成されたボリューム中の画像アーティファクトにつながる。画像アーティファクトは、Ｘ線検出器の感度が、Ｘ線減衰を物理的に十分正確にイメージングするのに十分でないとき、発生する。これは特に、ノイズがそこで優勢である、強く吸収するいわゆるＸ線不透過構造の背後にある問題につながる。これは、ボリューム中のストライプアーティファクト（stripe artifact）としてしばしば現れる、再構成プロセスにおける不整合な値につながる。

簡単のためにここでは金属アーティファクト低減（ＭＡＲ：metal artifact reduction）プロシージャと呼ばれる、これらの画像アーティファクトを補正するためのよく知られているソフトウェアプロシージャがある。それらのプロシージャは、再構成されたボリューム全体にわたって金属アーティファクトが低減されるように、サイノグラム中の金属領域を、妥当と思われる値と置き換える。再構成の後に、ボリューム中の金属領域が、妥当と思われる値と置き換えられる。

知られている金属アーティファクト低減技法は、患者がＤＶＴ収集中に移動するか、またはデバイスが不十分に較正された場合、低減された品質を与える。これの理由は、ＭＡＲ方法の多くが、ボリューム中で検出された金属エリアをサイノグラム中に投影することを行うか、またはその逆を行うからである。患者の移動または不十分なデバイス較正の場合、アーティファクトの影響を受けたボリューム中の金属検出と金属領域の投影の両方が、投影ジオメトリの正しくない考慮により不正確である。再構成における不整合は、ボリューム中の顕著な残存金属アーティファクトにつながる。

不正確な投影ジオメトリにより、多すぎるまたは少なすぎる補正が、サイノグラム中でならびにボリューム中で行われるか、または補正が、間違った場所で行われる。不正確な投影ジオメトリに関する前述の問題を低減するために、検出された金属領域の拡大、または検出された金属領域の周りの不確実性範囲の拡大が、通常、実施される。これは、サイノグラムの物理的に厳密に測定されたエリアも置き換えられるので、ＭＡＲのさらなる不正確さにつながる。拡大強度は、投影ジオメトリの最大許容不正確さから得られる。

本発明の目的は、デジタルボリュームトモグラフィ（ＤＶＴ）イメージングをＸ線コーンビームコンピュータトモグラフィにおいて再構成するための方法を提供することであり、それにより従来技術の上述の欠点を克服する。

この目的は、請求項１に記載の方法によって達成される。従属請求項の主題は、さらなる発展および好ましい実施形態に関する。

本発明による方法は、ＤＶＴイメージングを再構成するためのものである。本方法は、動きアーティファクト補償（ＭＡＣ：motion artifact compensation）を使用して、患者イメージングのサイノグラムから投影ジオメトリを推定するステップと、推定された投影ジオメトリを用いて、第１のボリュームを再構成するステップと、第１のボリュームと推定された投影ジオメトリとを使用して、サイノグラム中の金属領域と第１のボリューム中の金属領域とを検出するステップと、サイノグラム中の金属領域を補正し、補正されたサイノグラムを生成するステップと、推定された投影ジオメトリと補正されたサイノグラムとを用いて、第２のボリュームを再構成するステップと、第２のボリューム中の金属領域を補正するステップとを備える。

本発明の顕著な有益な効果は、補正のより高い正確さによるものである金属アーティファクトの改善された補正と、再構成されたボリュームの結果として生じる改善された画像品質とである。

以下の説明では、本発明は、例示的な実施形態を参照しながらおよび図面を参照しながら、より詳細に説明される。

本発明の一実施形態による、フロー図。 本発明による方法が行われ得るコンピュータ援用ＤＶＴシステムを示す図。 従来技術からの、不十分に補正された金属アーティファクトをもつ、ＭＡＣを用いないＤＶＴボリュームの軸方向スライスを示す図。 本発明の一実施形態による、明確に輪郭が描かれた「金属」充填を有する、ＭＡＣを用いたＤＶＴボリュームの軸方向スライスを示す図。

図面に示されている参照番号は、例示的な実施形態の以下の説明において言及される、以下に記載される要素を指定する。

１． ＤＶＴシステム
２． Ｘ線デバイス
３． Ｘ線源
４． Ｘ線検出器
５． 制御ユニット
６． 頭部固定
７． バイト（bite）
８． コンピュータ
９． ディスプレイ
ａ． 顎アーチ
ｂ． 歯
ｃ． 金属充填
ｄ． 金属アーティファクト
本発明による方法（図１参照）は、歯科領域におけるデジタルボリュームトモグラフィ（ＤＶＴ）イメージングを再構成するために使用される。本方法は、（Ｓ１）動きアーティファクト補償（ＭＡＣ）を使用して、患者イメージングのサイノグラムから投影ジオメトリを推定するステップと、（Ｓ２）推定された投影ジオメトリを使用して、第１のボリュームを再構成するステップと、（Ｓ３）第１のボリュームと推定された投影ジオメトリとを使用して、サイノグラム中の金属領域と第１のボリューム中の金属領域とを検出するステップと、（Ｓ４）サイノグラム中の金属領域を補正し、補正されたサイノグラムを生成するステップと、（Ｓ５）推定された投影ジオメトリと補正されたサイノグラムとを使用して、第２のボリュームを再構成するステップと、（Ｓ６）第２のボリューム中の金属領域を補正するステップとを備える。

図３ａは、ＭＡＣを用いない従来技術に従って再構成された、ＤＶＴボリュームの軸方向スライスを示す。図３ｂは、ＭＡＣを用いる本発明による方法によって生成された、ＤＶＴボリュームの同じ軸方向スライスを示す。図３ｂでは、金属充填（ｃ）が、明確な輪郭により明白である。顎アーチ（ａ）と個々の歯（ｂ）の両方にわたる図３ａからの金属アーティファクト（ｄ）は、図３ｂ中で明らかに低減されている。

好ましい実施形態では、第１のボリュームは、分解能および／または使用された投影の数および／または使用された画像処理に関して、第２のボリュームとは異なる。ここで、第１のボリュームはもっぱら、ＭＡＲのための金属領域を検出するために使用され、第２のボリュームは、可視化およびさらなる診断のため使用される。したがって、第１のボリュームは、算出時間とメモリ消費とを節約するために、得られた画像品質に関してより低い要求で生成され得る。詳細には、第１のボリュームを算出するとき、ボクセル分解能はより低くなり得、投影のサブサンプリングが実施され得、投影および／またはボリュームに関する画像処理ステップが省略され得る。

別の好ましい実施形態では、ステップ（Ｓ３）におけるサイノグラム中の金属領域の検出は、２つのサブステップによって実施される。これらは、（Ｓ３．１）好ましくはしきい値および／または金属構造の形状を使用することによって、第１のボリューム中の金属領域を検出することと、（Ｓ３．２）推定された投影ジオメトリを用いて、検出された金属領域をサイノグラム中に投影することとである。金属領域の検出は、サイノグラム中の構造の重畳により、ボリューム中で、サイノグラム中よりも容易である。推定された投影ジオメトリを用いて、第１のボリューム中で検出された金属領域をサイノグラム中に投影することによって、サイノグラム中の金属領域の検出が達成される。

好ましい代替実施形態では、ステップ（Ｓ３）におけるサイノグラム中の金属領域の検出（４）は、３つのさらなるサブステップによって実施される。これらは、（Ｓ３．１）推定された投影ジオメトリを使用して、第１のボリュームのシミュレートされたサイノグラムを生成することと、（Ｓ３．２）シミュレートされたサイノグラム中の金属領域を検出し、検出された金属領域をサイノグラムに転送することと、（Ｓ３．３）推定された投影ジオメトリとサイノグラム中の検出された金属領域とを使用して、第１のボリューム中の金属領域を決定することとである。

推定された投影ジオメトリが、第１のボリュームの再構成においておよびシミュレートされたサイノグラムの生成において使用された場合、シミュレートされたサイノグラムからサイノグラムへの検出された金属領域の転送は、画像領域が同一であるので、単純である。シミュレートされたサイノグラムを生成するときに重畳構造の簡略表現が達成されるので、金属領域の検出は、シミュレートされたサイノグラム中で、サイノグラム中よりも容易である。サイノグラムについて推定された投影ジオメトリが、シミュレートされたサイノグラムの生成において使用されたので、シミュレートされたサイノグラム中で検出された金属領域は、サイノグラムに転送され得る。推定された投影ジオメトリを用いて、シミュレートされたサイノグラム中で検出された金属領域を第１のボリュームに投影することによって、第１のボリューム中の金属領域の検出が達成される。

さらなる好ましい実施形態では、ステップ（Ｓ４）におけるサイノグラム中の金属領域の補正は、（Ｓ４．１）サイノグラム中の金属領域を、新しいピクセル値で充填するサブステップと、新しいピクセル値が、隣接ピクセルから計算されるか、または人工ピクセル値に対応するかのいずれかである、（Ｓ４．２）サイノグラムのピクセル値と、前のステップ（Ｓ４．１）からの新しいピクセル値を加重ブレンドするサブステップとのうちの１つまたは複数によって実施される。サイノグラム中の金属領域は、強い吸収により、正しく物理的に測定され得ない。サイノグラム中の金属領域を、妥当と思われるピクセル値と置き換えることによって、再構成中の金属アーティファクトの発生は回避され得る。新しい妥当と思われるピクセル値は、人工値であるか、またはサイノグラム中の隣接ピクセルから計算され得る。元のピクセル値と、新しいピクセル値を加重ブレンドすることが、不確実な物理情報の一部を維持する。

さらなる好ましい実施形態では、ステップ（Ｓ６）における第２のボリューム中の金属領域の補正は、（Ｓ６．１）第２のボリューム中の金属領域を、再構成された第１のボリュームからの値または人工値で充填するサブステップと、（Ｓ６．２）加重様式で、第２のボリュームの値と、前のステップ（Ｓ６．１）からの新しい値をブレンドするサブステップとのうちの１つまたは複数によって実施される。

サイノグラム中の金属領域の置き換えにより、再構成された第２のボリューム中の金属領域中の値は、強く変造され、それらは、明らかに低すぎる吸収を示す。第１のボリュームからの値を再挿入すること、またはそれらを、妥当と思われる人工値と置き換えることによって、第２のボリューム中の金属領域の存在が可視化される。元のボクセル値と、新しいボクセル値を加重ブレンドすることによって、満足な画像印象が達成される。

別の好ましい実施形態では、ＭＡＣは、ステップ（Ｓ１）において、（Ｓ１．１）推定された投影ジオメトリを用いて、第３のボリュームを再構成するサブステップと、（Ｓ１．２）サイノグラムの投影画像を第３のボリュームと位置合わせすることによって、投影ジオメトリを推定するサブステップとを備え、それらは、好ましくは、収束基準が達せられるまで反復的に繰り返される。

患者の移動またはずれるデバイスの移動が、再構成されたボリュームとサイノグラムの投影画像との間のずれにつながる。位置合わせは、これらのずれを検出し、投影ジオメトリを調整することによってこれらのずれを補償する。投影ジオメトリの調整は、幾何学的パラメータを変動させることによって行われる。幾何学的パラメータは、好ましくは、内因性パラメータと外因性パラメータとを備え、ここで、内因性パラメータは、互いに対するＸ線検出器とＸ線源との投影特性を記述し、外因性パラメータは、ボリュームの投影画像および選択されたサブ領域ごとの回転と並進とからなる変換を記述する。幾何学的パラメータを推定するためのプロシージャとして、いわゆる「最適化方法」が使用され得る。調整された投影ジオメトリを用いた第３のボリュームの繰り返される再構成は、位置合わせのための基準ボリュームを更新する。可能な収束基準は、ａ）投影ジオメトリの変化がしきい値よりも小さいかどうか、ｂ）位置合わせの品質がしきい値よりも大きいかどうか、ｃ）反復ステップの数がしきい値よりも大きいかどうか、ｄ）算出時間がしきい値よりも大きいかどうかである。

好ましくは、ステップ（Ｓ１）におけるＭＡＣは、投影行列を使用して推定されるべき投影ジオメトリを記述し得る。好ましくは、ステップ（Ｓ１）において、ＭＡＣは、既存の幾何学的デバイス較正に基づいて投影ジオメトリを推定し、それを、サイノグラムの投影画像ごとの３つの回転パラメータと３つの並進パラメータとによって記述されるさらなる剛体運動（rigid motion）で補足し得る。随意に、内因性パラメータがさらに推定され得る。

別の好ましい実施形態では、ステップ（Ｓ１．２）における投影画像の位置合わせが、推定された投影ジオメトリを使用して、シミュレートされたサイノグラムを生成することと、類似性測度を使用して、シミュレートされたサイノグラムをサイノグラムと比較することとによって実施される。シミュレートされたサイノグラムを生成する際に、推定された投影ジオメトリを使用して、第３のボリュームの順投影が実施される。正しい投影ジオメトリが適用されたとき、シミュレートされたサイノグラムとサイノグラムとの間の類似性は最大である。いくつかの類似性測度が、類似性を評価するための文献において知られている。類似性測度の代わりに、否定差測度（negated difference measure）も使用され得る。以下の類似性測度または誤差測度、すなわち、平均２乗誤差、平均絶対差、正規化相互相関、勾配相関、勾配情報、線形スケーリングを伴う勾配情報、勾配配向、相互情報量（Mutual Information）が、位置合わせプロシージャにおいて適用され得る。

別の好ましい代替実施形態では、ステップ（Ｓ１）におけるＭＡＣは、サイノグラムにおける１つまたは複数の整合性制約に関するメトリックを評価することによって、投影ジオメトリを推定することを備える。サイノグラムにおける整合性制約に関するメトリックを評価することによって投影ジオメトリを推定することは、ボリュームの再構成が必要とされないという利点を有する。

以下の整合性制約、すなわち、データ整合性、エピポーラ整合性、フーリエ整合性、Ｇｒａｎｇｅａｔの定理、相互相関が、投影ジオメトリ推定プロシージャにおいて適用され得る。

別の好ましい代替実施形態では、ステップ（Ｓ１）におけるＭＡＣは、（Ｓ１．１）推定された投影ジオメトリを用いて、第３のボリュームを再構成するサブステップと、（Ｓ１．２）シャープネス測度、または滑らかさ条件を評価することなど、第３のボリュームに関する画像品質メトリックを算出するサブステップと、（Ｓ１．３）推定された投影ジオメトリを、画像品質メトリックを改善するように調整するサブステップとを備え、それらは、好ましくは反復的に繰り返される。画像品質メトリックは、しばしば、動きアーティファクトの存在と相関する。投影ジオメトリを調整することによって画像品質メトリックを改善することは、動きアーティファクトを低減することができる。投影ジオメトリの適応は、好ましくは反復的プロシージャにおいて、幾何学的パラメータを変動させることによって行われる。

画像品質メトリックは、勾配分散または勾配ノルムなど、ボリュームに関するシャープネス測度であり、あるいは、グレー値のエントロピー、グレー値の分散、または総変動など、ボリュームに関する滑らかさ制約によって定義され得る。

本発明による方法は、コンピュータ実装可能方法であり、コンピュータ化ＤＶＴシステム（１）上で実行され得る。図２は、ＤＶＴシステム（１）のための一実施形態の一例を示す。この点について、本発明はまた、コンピュータ可読コードを有するコンピュータプログラムを含む。コンピュータプログラムは、データストレージデバイス上に設けられ得る。コンピュータ支援ＤＶＴシステム（１）は、患者イメージングを実施するためのＸ線デバイス（２）を備え、それにより、サイノグラムが生成される。Ｘ線デバイス（２）は、曝露中に患者の頭部の周りで回転される、Ｘ線源（３）とＸ線検出器（４）とを有する。曝露中のＸ線源（３）とＸ線検出器（４）との軌道は、円形経路を描き得る。代替的に、軌道は、円形経路からずれる形状を仮定し得る。いくつかのアクチュエータが同時に制御される場合、純粋な円形経路からずれる患者の頭部の周りのデバイス軌道が、達成され得る。患者の頭部は、バイト（７）および頭部固定（６）によりＸ線デバイス中に配置される。コンピュータ化ＤＶＴシステム（１）は、制御ユニット（５）と、好ましくは、Ｘ線デバイス（２）に接続可能なコンピュータ（８）またはコンピューティングユニットと、好ましくは、とりわけデータセットを可視化するためのディスプレイ（９）とを備える。コンピュータ（８）は、ローカルエリアネットワーク（図示せず）を介して、または代替的に、インターネットを介して、Ｘ線デバイス（２）に接続され得る。コンピュータ（８）は、クラウドの一部であり得る。代替的に、コンピュータ（８）は、Ｘ線デバイス（２）に組み込まれ得る。算出は、代替的に、クラウドにおいて行われ得る。コンピュータ（８）は、コンピュータプログラムを実行し、ディスプレイ（９）上での可視化のためのものを含む、データセットを提供する。ディスプレイ（９）は、Ｘ線デバイス（２）から空間的に分離され得る。好ましくは、コンピュータ（８）はまた、Ｘ線デバイス（２）を制御し得る。代替的に、別個のコンピュータが、制御および再構成のために使用され得る。

本発明によれば、上記の実施形態によって生成されたデータセットは、好ましくはディスプレイ（９）またはプリントアウトによって、可視化のために、特に診断目的のために、医師に提示され得る。

本発明によれば、上記の実施形態によって生成されたデータセットは、好ましくはディスプレイ（９）またはプリントアウトによって、可視化のために、特に診断目的のために、医師に提示され得る。
以下に本願の出願当初の特許請求の範囲に記載された発明を付記する。
［Ｃ１］
歯科領域におけるデジタルボリュームトモグラフィ（ＤＶＴ）イメージングを再構成するための方法であって、前記方法が、
（Ｓ１）動きアーティファクト補償（ＭＡＣ）を使用して、患者イメージングのサイノグラムから投影ジオメトリを推定するステップと、
（Ｓ２）前記推定された投影ジオメトリを用いて、第１のボリュームを再構成するステップと、
（Ｓ３）前記第１のボリュームと前記推定された投影ジオメトリとを使用して、前記サイノグラム中の金属領域と前記第１のボリューム（１）中の金属領域とを検出するステップと、
（Ｓ４）前記サイノグラム中の前記金属領域を補正し、補正されたサイノグラムを生成するステップと、
（Ｓ５）前記推定された投影ジオメトリと前記補正されたサイノグラムとを用いて、第２のボリュームを再構成するステップと、
（Ｓ６）前記第２のボリューム中の前記金属領域を補正するステップと
を備えることを特徴とする、方法。
［Ｃ２］
前記第１のボリュームが、分解能および／または使用された投影の数および／または使用された画像処理に関して、前記第２のボリュームとは異なることを特徴とする、Ｃ１に記載の方法。
［Ｃ３］
ステップ（Ｓ３）における前記サイノグラム中の前記金属領域の前記検出が、
（Ｓ３．１）好ましくはしきい値および／または金属構造の形状を使用することによって、前記第１のボリューム中の前記金属領域を検出することと、
（Ｓ３．２）前記推定された投影ジオメトリを用いて、前記検出された金属領域を前記サイノグラム中に投影することと
からなることを特徴とする、Ｃ１または２に記載の方法。
［Ｃ４］
ステップ（Ｓ３）における前記サイノグラム中の前記金属領域の前記検出が、
（Ｓ３．１）前記推定された投影ジオメトリを使用して、前記第１のボリュームのシミュレートされたサイノグラムを生成することと、
（Ｓ３．２）前記シミュレートされたサイノグラム中の前記金属領域を検出し、前記検出された金属領域を前記サイノグラムに転送することと、
（Ｓ３．３）前記推定された投影ジオメトリと前記サイノグラム中の前記検出された金属領域とを使用して、前記第１のボリューム中の前記金属領域を決定することと
からなることを特徴とする、Ｃ１または２に記載の方法。
［Ｃ５］
ステップ（Ｓ４）における前記サイノグラム中の前記金属領域の前記補正は、
（Ｓ４．１）前記サイノグラム中の前記金属領域を、新しいピクセル値で充填するステップと、前記新しいピクセル値が、隣接ピクセルから計算されるか、または人工ピクセル値に対応するかのいずれかである、
（Ｓ４．２）前記サイノグラムのピクセル値と、前のステップ（Ｓ４．１）からの前記新しいピクセル値を加重ブレンドするステップと
のうちの１つまたは複数からなることを特徴とする、Ｃ１から４のいずれか一項に記載の方法。
［Ｃ６］
ステップ（Ｓ６）における前記第２のボリューム中の前記金属領域の前記補正が、
（Ｓ６．１）前記第２のボリューム中の前記金属領域を、前記再構成された第１のボリュームからの値または人工値で充填するステップと、
（Ｓ６．２）前記第２のボリュームの値と、前のステップ（Ｓ６．１）からの前記新しい値を加重ブレンドするステップと
のうちの１つまたは複数からなることを特徴とする、Ｃ１から５のいずれか一項に記載の方法。
［Ｃ７］
ステップ（Ｓ１）における前記ＭＡＣは、
（Ｓ１．１）推定された投影ジオメトリを用いて、第３のボリュームを再構成するステップと、
（Ｓ１．２）前記サイノグラムの投影画像を前記第３のボリュームと位置合わせすることによって、前記投影ジオメトリを推定するステップと
を備え、好ましくは、収束基準が達せられるまで反復的に繰り返すことを特徴とする、Ｃ１から６のいずれか一項に記載の方法。
［Ｃ８］
ステップ（Ｓ１．２）における前記投影画像の前記位置合わせが、前記推定された投影ジオメトリを用いて、シミュレートされたサイノグラムを生成することと、類似性測度によって、前記シミュレートされたサイノグラムを前記サイノグラムと比較することとによって実施されることを特徴とする、Ｃ７に記載の方法。
［Ｃ９］
ステップ（Ｓ１）における前記ＭＡＣが、前記サイノグラムにおける１つまたは複数の整合性制約に関するメトリックを評価することによって、前記投影ジオメトリの推定を実施することを特徴とする、Ｃ１から６のいずれか一項に記載の方法。
［Ｃ１０］
ステップ（Ｓ１）における前記ＭＡＣが、
（Ｓ１．１）前記推定された投影ジオメトリを用いて、第３のボリュームを再構成するステップと、
（Ｓ１．２）シャープネス測度、または滑らかさ条件の評価など、前記第３のボリュームに関する画像品質メトリックを算出するステップと、
（Ｓ１．３）前記推定された投影ジオメトリを、前記画像品質メトリックを改善するように調整するステップと
を備え、好ましくは、反復的に繰り返すことを特徴とする、Ｃ１から６のいずれか一項に記載の方法。
［Ｃ１１］
コンピュータ支援ＤＶＴシステム（１）によって実行されたとき、前記コンピュータ支援ＤＶＴシステム（１）に、Ｃ１から１０のいずれか一項に記載の方法ステップを実施させるコンピュータ可読コードを備えるコンピュータプログラム。
［Ｃ１２］
Ｘ線デバイス（２）と、Ｃ１１に記載のコンピュータプログラムを実行するように構成されたコンピューティングユニット（８）とを備えるコンピュータ化ＤＶＴシステム（１）。
［Ｃ１３］
Ｃ１から１０のいずれか一項によって提供される、可視化のためのデータセットの使用。

Claims (13)

1. 歯科領域におけるデジタルボリュームトモグラフィ（ＤＶＴ）イメージングを再構成するための方法であって、前記方法が、
   （Ｓ１）動きアーティファクト補償（ＭＡＣ）を使用して、患者イメージングのサイノグラムから投影ジオメトリを推定するステップと、
   （Ｓ２）前記推定された投影ジオメトリを用いて、第１のボリュームを再構成するステップと、
   （Ｓ３）前記第１のボリュームと前記推定された投影ジオメトリとを使用して、前記サイノグラム中の金属領域と前記第１のボリューム（１）中の金属領域とを検出するステップと、
   （Ｓ４）前記サイノグラム中の前記金属領域を補正し、補正されたサイノグラムを生成するステップと、
   （Ｓ５）前記推定された投影ジオメトリと前記補正されたサイノグラムとを用いて、第２のボリュームを再構成するステップと、
   （Ｓ６）前記第２のボリューム中の前記金属領域を補正するステップと
   を備えることを特徴とする、方法。
2. 前記第１のボリュームが、分解能および／または使用された投影の数および／または使用された画像処理に関して、前記第２のボリュームとは異なることを特徴とする、請求項１に記載の方法。
3. ステップ（Ｓ３）における前記サイノグラム中の前記金属領域の前記検出が、
   （Ｓ３．１）好ましくはしきい値および／または金属構造の形状を使用することによって、前記第１のボリューム中の前記金属領域を検出することと、
   （Ｓ３．２）前記推定された投影ジオメトリを用いて、前記検出された金属領域を前記サイノグラム中に投影することと
   からなることを特徴とする、請求項１または２に記載の方法。
4. ステップ（Ｓ３）における前記サイノグラム中の前記金属領域の前記検出が、
   （Ｓ３．１）前記推定された投影ジオメトリを使用して、前記第１のボリュームのシミュレートされたサイノグラムを生成することと、
   （Ｓ３．２）前記シミュレートされたサイノグラム中の前記金属領域を検出し、前記検出された金属領域を前記サイノグラムに転送することと、
   （Ｓ３．３）前記推定された投影ジオメトリと前記サイノグラム中の前記検出された金属領域とを使用して、前記第１のボリューム中の前記金属領域を決定することと
   からなることを特徴とする、請求項１または２に記載の方法。
5. ステップ（Ｓ４）における前記サイノグラム中の前記金属領域の前記補正は、
   （Ｓ４．１）前記サイノグラム中の前記金属領域を、新しいピクセル値で充填するステップと、前記新しいピクセル値が、隣接ピクセルから計算されるか、または人工ピクセル値に対応するかのいずれかである、
   （Ｓ４．２）前記サイノグラムのピクセル値と、前のステップ（Ｓ４．１）からの前記新しいピクセル値を加重ブレンドするステップと
   のうちの１つまたは複数からなることを特徴とする、請求項１から４のいずれか一項に記載の方法。
6. ステップ（Ｓ６）における前記第２のボリューム中の前記金属領域の前記補正が、
   （Ｓ６．１）前記第２のボリューム中の前記金属領域を、前記再構成された第１のボリュームからの値または人工値で充填するステップと、
   （Ｓ６．２）前記第２のボリュームの値と、前のステップ（Ｓ６．１）からの前記新しい値を加重ブレンドするステップと
   のうちの１つまたは複数からなることを特徴とする、請求項１から５のいずれか一項に記載の方法。
7. ステップ（Ｓ１）における前記ＭＡＣは、
   （Ｓ１．１）推定された投影ジオメトリを用いて、第３のボリュームを再構成するステップと、
   （Ｓ１．２）前記サイノグラムの投影画像を前記第３のボリュームと位置合わせすることによって、前記投影ジオメトリを推定するステップと
   を備え、好ましくは、収束基準が達せられるまで反復的に繰り返すことを特徴とする、請求項１から６のいずれか一項に記載の方法。
8. ステップ（Ｓ１．２）における前記投影画像の前記位置合わせが、前記推定された投影ジオメトリを用いて、シミュレートされたサイノグラムを生成することと、類似性測度によって、前記シミュレートされたサイノグラムを前記サイノグラムと比較することとによって実施されることを特徴とする、請求項７に記載の方法。
9. ステップ（Ｓ１）における前記ＭＡＣが、前記サイノグラムにおける１つまたは複数の整合性制約に関するメトリックを評価することによって、前記投影ジオメトリの推定を実施することを特徴とする、請求項１から６のいずれか一項に記載の方法。
10. ステップ（Ｓ１）における前記ＭＡＣが、
    （Ｓ１．１）前記推定された投影ジオメトリを用いて、第３のボリュームを再構成するステップと、
    （Ｓ１．２）シャープネス測度、または滑らかさ条件の評価など、前記第３のボリュームに関する画像品質メトリックを算出するステップと、
    （Ｓ１．３）前記推定された投影ジオメトリを、前記画像品質メトリックを改善するように調整するステップと
    を備え、好ましくは、反復的に繰り返すことを特徴とする、請求項１から６のいずれか一項に記載の方法。
11. コンピュータ支援ＤＶＴシステム（１）によって実行されたとき、前記コンピュータ支援ＤＶＴシステム（１）に、請求項１から１０のいずれか一項に記載の方法ステップを実施させるコンピュータ可読コードを備えるコンピュータプログラム。
12. Ｘ線デバイス（２）と、請求項１１に記載のコンピュータプログラムを実行するように構成されたコンピューティングユニット（８）とを備えるコンピュータ化ＤＶＴシステム（１）。
13. 請求項１から１０のいずれか一項によって提供される、可視化のためのデータセットの使用。

Images