JP2013085965A - 円軌道コーンビームコンピュータ断層撮影（ｃｔ）においてアーチファクトを大幅に軽減するための方法およびシステム - Google Patents

Abstract

【課題】円軌道コーンビームコンピュータ断層撮影（ＣＴ）においてアーチファクトを大幅に軽減すること。
【解決手段】円軌道コンピュータ断層撮影においてコーンビームアーチファクトを補正する方法は、ａ）所定の円形線源軌道における実測の投影データから参照画像を再構成するステップと、ｂ）前記円形線源軌道に対応する理論的な完全円形軌道に従って前記参照画像を順投影することにより投影データを生成するステップと、ｃ）前記生成された投影データから補正画像を再構成するステップと、ｄ）前記参照画像と前記補正画像とを用いて最終画像を生成することによって前記コーンビームアーチファクトを大幅に軽減するステップとを備える。
【選択図】 図３

Description

本発明は、一般に画像処理およびシステムに関し、より詳細には、円軌道コーンビームコンピュータ断層撮影（ＣＴ）におけるアーチファクトの大幅な軽減に関する。

コーンビームアーチファクトは、円軌道コンピュータ断層撮影ではよく知られている問題である。東芝のＡｑｕｉｌｉｏｎＯＮＥなどのほとんどの高度なＣＴシステムにおけるＸ線源のコーン角は非常に大きく、ラドン領域のデータが欠損しているために、ＣＦＫ画像でコーンビームアーチファクトが生じる。

従来技術による１つの手法では、円軌道コーンビームＣＴにおけるアーチファクトは、円形軌道と線軌道の組み合わせなどの理論的に完全な軌道に基づく正確な再構成を適用することによって大幅に除去される。追加の線スキャンによって、正確な再構成のための円軌道を有する理論的に完全な軌道が得られるが、追加のスキャンは、収集が使用不可能かまたは非実用的のいずれかであることが多い。そのうえ、円形データと線データは同時に得られないので、２回のスキャンの間でモーションまたは造影剤の増強効果（agent enhancement）が変化すると、２つのデータセット間でデータの不一致が発生し、画像の精度が影響を受ける。最後に、追加スキャンにより、患者が曝露される放射線量が増加する。

上記の理由から、円形データのみからの画像ボリュームの正確な再構成が特に注目されている。従来技術による別の手法では、線データを推定するためにスキャノグラムが用いられる。この手法では患者の放射線量は増加しないが、推定された線データがアーチファクトの大半を軽減するのに役立つ場合でも、コーンビームアーチファクトが依然として認められる。同時にまた、モーションまたは造影剤の増強効果が変化すると、結果として得られる画像がいくらか不正確になる。

目的は、円軌道コーンビームコンピュータ断層撮影（ＣＴ）においてアーチファクトを大幅に軽減することにある。

本実施形態による円軌道コンピュータ断層撮影においてコーンビームアーチファクトを補正する方法は、
ａ）所定の円形線源軌道における実測の投影データから参照画像を再構成するステップと、
ｂ）前記円形線源軌道に対応する理論的な完全円形軌道に従って前記参照画像を順投影することにより投影データを生成するステップと、
ｃ）前記生成された投影データから補正画像を再構成するステップと、
ｄ）前記参照画像と前記補正画像とを用いて最終画像を生成することによって前記コーンビームアーチファクトを大幅に軽減するステップとを備える。

本実施形態により円軌道上で収集されたデータから再構成された画像においてコーンビームアーチファクトを大幅に軽減するためのマルチスライスＸ線ＣＴ装置またはスキャナの一実施形態を示す図。 本実施形態によってコーンビーム（ＣＢ）アーチファクトが大幅に軽減される理由の一側面を示す図。 本実施形態によりコーンビーム（ＣＢ）アーチファクトを大幅に軽減する例示的な処理に関与する全般的なステップを示す流れ図。 本実施形態によってコーンビームアーチファクトを大幅に軽減するための特定の概念的解決策を示す図の集まり。 本実施形態におけるコーンビームアーチファクト軽減処理に関与するステップを示すフロー図。 本実施形態における大幅なコーンビームアーチファクト軽減の効果を評価するために用いられるディスクファントムの側面図を示す図。 本実施形態における大幅なコーンビームアーチファクト軽減の効果を評価するために用いられるディスクファントムの上面図を示す図。 コーンビームアーチファクト補正を行わずにディスクファントムから再構成された従来のＦＤＫ画像。 コーンビームアーチファクト補正を行わずに５００回反復した後のディスクファントムから再構成されたＳＡＲＴ画像。 本実施形態によりコーンビームアーチファクトを大幅に軽減するための１つの例示的処理によりディスクファントムから再構成された補正された画像。 リビニング（rebinning）を行わずにＣＦＫを用いた再構成画像。 リビニングは行わないがＯＳＲフィルタ処理を行ったＣＦＫを用いた再構成画像。 比較のためにＨＦＫを用いた再構成画像。 リビニングを行ったＣＦＫを用いた再構成画像。 本実施形態のコーンビームアーチファクトを大幅に軽減する処理によりリビニングと１回順投影された線データ補正とを行ったＣＦＫを用いた再構成画像。 本実施形態のコーンビームアーチファクトを大幅に軽減する処理によりリビニングと３回反復順投影された線データ補正とを行ったＣＦＫを用いた再構成画像。

ここで図面を参照すると、同じ参照番号はすべての図面を通して対応する構造を示すが、特に図１を参照すると、図は、ガントリ１００と他のデバイスまたはユニットとを含む、本実施形態によるマルチスライスＸ線ＣＴ装置またはスキャナの一実施形態を示す。ガントリ１００は、側面図にて示されており、Ｘ線管１０１と、環状フレーム１０２と、多列または２次元配列型のＸ線検出器１０３とをさらに含む。Ｘ線管１０１およびＸ線検出器１０３は、環状フレーム１０２上の被検体（subject）Ｓを横切る対角線上に取り付けられ、環状フレーム１０２は、回転軸ＲＡのまわりに回転可能に支持される。回転ユニット１０７は、フレーム１０２を０．４秒／回転などの高速で回転させ、被検体Ｓは、軸ＲＡに沿って、図示の紙面に対する奥行き方向または手前方向に移動される。

マルチスライスＸ線ＣＴ装置は、Ｘ線管１０１がＸ線を生成するようにスリップリング１０８を介してＸ線管１０１に印加されるべき管電圧を生成する高電圧発生装置１０９をさらに含む。Ｘ線は被検体Ｓに向かって照射され、被検体Ｓの断面領域は円によって表されている。Ｘ線検出器１０３は、被検体Ｓを透過した照射Ｘ線を検出するために、被検体Ｓを挟んでＸ線管１０１の反対側にある。

さらに図１を参照すると、Ｘ線ＣＴ装置またはスキャナは、Ｘ線検出器１０３からの検出された信号を処理するための他のデバイスをさらに含む。データ収集回路またはデータ収集システム（ＤＡＳ）１０４は、各チャンネルに対してＸ線検出器１０３から出力された信号を電圧信号に変換し、それを増幅して、さらにデジタル信号に変換する。Ｘ線検出器１０３およびＤＡＳ１０４は、１回転あたりの所定の全投影数（ＴＰＰＲ：total number of projections per rotation）を、最大で９００ＴＰＰＲ、９００ＴＰＰＲから１８００ＴＰＰＲの間、および９００ＴＰＰＲから３６００ＴＰＰＲの間にできるよう処理するように構成される。

上述のデータは、非接触データ伝送装置１０５を通して、ガントリ１００の外部のコンソールに収容された前処理デバイス１０６に送られる。前処理デバイス１０６は、生データに対して感度補正などの特定の補正を実行する。次に、その結果として生成されるデータ（投影データとも呼ばれるデータ）を、記憶デバイス１１２が再構成処理の直前の段階で記憶する。記憶デバイス１１２は、データ／制御バスを介して、再構成デバイス１１４、入力デバイス１１５、表示デバイス１１６、マルチスケール処理デバイス１１７、およびスキャンプラン支援装置２００と共に、システムコントローラ１１０に接続される。スキャンプラン支援装置２００は、撮像技師がスキャンプランを策定するのを支援するための機能を含む。

再構成デバイス１１４の一実施形態は、種々のソフトウェア構成要素とハードウェア構成要素とをさらに含み、投影データに対する所定の解析再構成処理を実行する。本実施形態によれば、ＣＴ装置の再構成デバイス１１４は、有利には、所定のフィルタ補正逆投影（ＦＢＰ）法を用いて画像ボリュームを再構成する。

他の実施形態によれば、ＣＴ装置の再構成デバイス１１４は、有利には、逐次再構成法を用いて全変動（ＴＶ）を最小化する。一般に、本実施形態の再構成デバイス１１４は、全変動逐次再構成（ＴＶＩＲ）アルゴリズムを実施する。ＴＶＩＲは、投影データに対して、オーダードサブセット同時代数的再構成法（ＯＳ−ＳＡＲＴ）ステップなどの同時的代数的再構成法と、ＴＶ最小化ステップなどの調整とを実行する。一実施形態では、この２つのステップは、反復回数があらかじめ設定されているメインループにおいて連続して実施される。

ＴＶ最小化ステップの前に、投影データは、オーダードサブセット同時代数的再構成法（ＯＳ−ＳＡＲＴ）に供される。投影データは、それぞれ特定数のビューを有する所定数のサブセットＮにグループ化される。一実施形態では、オーダードサブセット同時代数的再構成法（ＯＳ−ＳＡＲＴ）の実施中に、各サブセットが連続して処理されてもよい。別の実施形態では、複数のサブセットは、複数の中央処理装置（ＣＰＵ）または１つのグラフィック処理装置（ＧＰＵ）などの特定のマイクロプロセッサを利用して並列に処理されてもよい。全変動（ＴＶ）最小化ステップにおいて、再構成デバイス１１４の一実施形態では、現在の画像ボリュームの目的関数が前の画像ボリュームの目的関数より小さくなるように正のステップ幅（step size）を探索するために直線探索方式（line search strategy）を用いる。

再構成デバイス１１４はまた、オーダードサブセット同時代数的再構成法（ＯＳ−ＳＡＲＴ）の実施中に、２つの主要な演算を実行する。すなわち、再構成デバイス１１４は、各サブセットＮに対して、計算上の投影データを生成するために画像ボリュームを再投影し、更新画像ボリュームを再構成するために、計測された投影データと計算上の投影データとの正規化された差分を逆投影する。さらに詳細には、再構成デバイス１１４の一実施形態では、システムマトリックスの係数がキャッシュされないレイトレーシング法を用いることによって、画像ボリュームを再投影する。さらに、再構成デバイス１１４の一実施形態では、サブセット内のすべてのレイを同時に再投影する。この処理は任意選択で並列に実施される。逆投影において、再構成デバイス１１４の一実施形態では、所望の更新画像ボリュームを生成する目的でサブセット内のすべての正規化された差分投影データを逆投影するためにピクセルドリブン（pixel-driven）法を用いる。再構成デバイス１１４が、画像ボリュームを生成するためにサブセット内のすべてのレイサム（すなわち差分投影データ）を逆投影するので、この演算も任意選択で並列に実施される。これらの演算は、単一のＯＳ−ＳＡＲＴステップを完了するために、すべてのサブセットＮに適用される。さらに、ＡＷＡＤは任意選択で組み合わされる。

本実施形態は、上述の構成要素に加えて、コーンビームアーチファクト軽減を実行するための種々の他のソフトウェアモジュールとハードウェア構成要素とをさらに含む。本実施形態によれば、ＣＴ装置のコーンビーム（ＣＢ）アーチファクト軽減デバイス１１７は、有利には、特定の状況下でコーンビームアーチファクトを大幅に軽減するためのコーンビームアーチファクト軽減機能を実行する。一般に、ＣＢアーチファクトは、画質を劣化させる、シェーディング（shading）コーンビームと高コントラストコーンビームとを含む２つの成分を有する。１つの例示的な処理では、画質を向上させるため、シェーディングはリビニングをフィルタ処理することによって補正され、高コントラストは線画像によって補正される。上述のアーチファクトについては、本願の他の図を参照してさらに説明する。

他の実施形態では、ＣＴ装置のコーンビーム（ＣＢ）アーチファクト軽減デバイス１１７は、有利には、コーンビームアーチファクトを大幅に軽減するために正確な再構成の要素と、ＳＡＲＴなどの逐次反復再構成の要素とを結合する。より詳細に説明するように、再構成デバイス１１４は、円形線源軌道上のコーンビーム線源を用いて収集された計測投影データから円形画像を再構成する。その後、ＣＢアーチファクト軽減デバイス１１７が、この円形画像からの線データを順投影し、次に、再構成デバイス１１４が、この順投影された線データに基づいて線画像を再構成する。上記のステップに基づいて、ＣＢアーチファクト軽減デバイス１１７が円形画像と線画像とを合成する。この合成画像は、ＣＢアーチファクトは、従来の方法による再構成画像に存在するはずのＣＢアーチファクトに対し、大幅に軽減されたＣＢアーチファクトを含む。最後に、ＣＢアーチファクト軽減デバイス１１７が、補正された画像を出力する。

本実施形態では、コーンビームアーチファクト軽減デバイス１１７は、他のソフトウェアモジュールおよび／または記憶デバイス１１２、再構成デバイス１１４、表示デバイス１１６、および入力デバイス１１５などのシステム構成要素にデータ／制御バスを介して動作的に接続される。この点に関して、他の実施形態では、コーンビームアーチファクト軽減デバイス１１７単独では、必ずしもコーンビームアーチファクト機能および／またはそれに関連するタスクを実行しない。そのうえ、他の実施形態では、コーンビームアーチファクト軽減デバイス１１７は、任意選択で再構成デバイス１１４などの他のデバイスの一部である。

図２は、本実施形態によってコーンビーム（ＣＢ）アーチファクトを大幅に軽減する理由の１つの側面を示す図である。この図は、一連の細い対象０が、所定の円軌道ＣＴ上にある線源に対して相対的な位置にある例示的な状況を示す。線源位置Ｓ１にあるコーンビームは、所定のコーンビーム角で対象０に向かって照射する。この例示的な状況は、画像は診断のために良好な分解能を有するズームされたまたは所望の撮像視野（ＦＯＶ）内で再構成されることが多いことも示す。順投影には、Ｘ線ビームを減衰させる対象全体に関する情報が必要なので、図２に示すように２つのボリューム画像ＣＦＫ＿ＡおよびＣＦＫ＿Ｂが生成される。網掛け領域は、不十分に収集された実測された投影データに対応する。画像ＣＦＫ＿Ａは、最終結果に用いられるべき所望のＦＯＶを有し、画像ＣＦＫ＿Ｂは、欠損している線データを生成するために全（full）ＦＯＶを有する。必要な範囲２Ｈは、この線データを生成するための線データ方向ＬＴＤに沿って示されている。

次に図３を参照すると、流れ図は、本実施形態によりコーンビーム（ＣＢ）アーチファクトを大幅に軽減する例示的な処理に関与する全般的なステップを示す。実際、この流れ図は、本実施形態によりコーンビーム（ＣＢ）アーチファクトを大幅に軽減する例示的な処理の概念的なスキームであり、本発明は、必ずしも以下の説明に限定されない。ステップＳ１００では、実測された投影データが収集される。一般に、収集された実測された投影データは、コーンビーム円形スキャンジオメトリが原因で不十分であるので、従来の方法で再構成された場合、コーンビームアーチファクトを発生しやすい。ステップＳ１１０では、円形画像が、円形線源軌道上のコーンビーム線源を用いて収集された実測された投影データから再構成される。再構成された円形画像は、任意選択で、後で取り出すために所定の記憶域に記憶される。ステップＳ１２０では、ステップＳ１１０で再構成された円形画像から、線データが順投影または再投影される。ステップＳ１２０での線データ生成の詳細については、後で他の流れ図を参照して説明する。ステップＳ１３０では、線画像が、ステップＳ１２０で生成された順投影された線データに基づいて再構成される。上記のステップＳ１１０、Ｓ１２０、およびＳ１３０に基づいて、ステップＳ１４０で、円形画像と線画像がここで合成される。この合成画像は、従来の方法による再構成画像に存在するはずのＣＢアーチファクトに対し、大幅に軽減されたＣＢアーチファクトを含む。最後に、ステップＳ１５０で、補正された画像が表示または解析のために出力される。ＣＢアーチファクトを大幅に軽減する別の実施形態では、上述のステップのうちいくつかは、線データの精度を向上させ、さらに出力画像のアーチファクト軽減を向上させるために、ＳＡＲＴなどの既知の逐次的技法を用いて反復的に繰り返される。

図４は、本実施形態によってコーンビームアーチファクトを大幅に軽減するための特定の概念的解決策を示す図の集まりである。図２に関して既に説明したように、例示的な状況では、診断のために良好な分解能を有するズームされたまたは所望の撮像視野（ＦＯＶ）が必要である。順投影には、Ｘ線ビームを減衰させる対象全体に関する情報が必要なので、図４Ａおよび図４Ｂに示すような２つの画像ＣＦＫ＿ＡおよびＣＦＫ＿Ｂは、図４Ｃに示すように円軌道ＣＴ上のコーンビーム線源ＣＢＳを用いて収集された対応する円形データから円形フェルドカンプ（ＣＦＫ：circular Feld-Kamp）法を用いて再構成される。図４Ｃは、画像ボリュームＣＦＫ＿Ａのための円形データが画像ＣＦＫ＿Ｂのための円形データのズームされた一部分であることも示す。つまり、画像ＣＦＫ＿Ａは、画像ＣＦＫ＿Ｂの全ＦＯＶまたは大きなＦＯＶの中に所望のＦＯＶを有し、この全ＦＯＶはＺ軸方向に沿って延びる。

本実施形態によれば、図４Ｄに示すような線データは、フィルタ補正逆投影の行われたボリューム画像ＣＦＫ＿Ｂの順投影から生成される。この点に関して、本願では、再投影は、上述の線データ生成において順投影と同義的に用いられる。上述のように、画像ＣＦＫ＿Ｂは、全ＦＯＶを有する。大きなコーン角を有するＸ線は、再構成画像を越えてＺ方向の空間を透過できるので、画像ＣＦＫ＿Ｂは、任意選択で、スキャンステムのコーン角によって決まるあるＺ範囲を越えて延びる。他の実施形態では、画像ＣＦＫ＿Ｂにおけるコーンビームアーチファクトを軽減するため、線軌道ＬＴに沿った順投影の前に、所定の適応ローパス３Ｄフィルタが画像ＣＦＫ＿Ｂに適用される。あるいは、別の実施形態では、線軌道ＬＴに沿った順投影の前に、所定の因数分解手法が画像ＣＦＫ＿Ｂに適用される。

本実施形態によるコーンビームアーチファクト軽減処理の一実施形態では、線データが、画像ＣＦＫ＿Ｂの再投影によって線軌道から得られた後、図４Ｅに示すような所望のＦＯＶを有する線画像が、図４Ｄに示すような線データから再構成される。最後に、図４Ｅに示すような所望のＦＯＶを有する線画像が、＋記号によって表されるように図４Ａに示すような画像ボリュームＣＦＫ＿Ａと合成され、図４Ｆに示すような補正された画像を生成する。この合成画像は、大幅に軽減されたコーンビームアーチファクトを有する。一実施形態では、画像ボリュームＣＦＫ＿Ｂは、任意選択で、全ＦＯＶを有する線データから再構成された画像を追加することによって更新される。そのうえ、この線データは、アーチファクトを有するボリューム画像から順投影されるので、近似されたデータである。このため、別の実施形態では、線データを改良し、かつコーンビームアーチファクトの軽減を向上させるために、反復手法が用いられる。

図５は、本実施形態におけるコーンビームアーチファクト軽減処理に関与するステップを示すフロー図である。コーンビームアーチファクトを大幅に軽減する例示的な処理では、実測された投影データは、特定のコーンビーム角を用いて所定の円形データ収集技法によって収集されたと仮定される。一般に、実施形態におけるコーンビームアーチファクト軽減処理では、所望の小さな撮像視野と、対応する円軌道コーンビームデータからの全撮像視野とをそれぞれ有する２つの参照画像ＣＦＫ＿ＡとＣＦＫ＿Ｂとを生成する。実施形態におけるコーンビームアーチファクト軽減処理では、所定の線源軌道に沿って全撮像視野を有する参照画像ＣＦＫ＿Ｂを順投影することによって投影データを生成し、次に、この投影データから補正画像を再構成する。この補正は、線軌道などの所定の線源軌道が参照画像ＣＦＫ＿Ｂの円形線源軌道を補うので達成される。そのうえ、実施形態におけるコーンビームアーチファクト軽減処理では、ステップＳ３６０の後の補正された画像が本実施形態により大幅に軽減されたコーンビームアーチファクトを有するように、線画像に関してステップＳ２３０からＳ２９０などの特定のステップを反復的に繰り返す。代替実施形態では、コーンビームアーチファクト軽減処理において、反復の代わりに所定の方法で線画像の操作を実行する。

さらに図５を参照して、コーンビームアーチファクトを大幅に軽減する処理について詳細にさらに説明する。一般に、参照画像は、Ｃｓｅｇ＋ｚ／ｃｏｓ（ｇａｍｍａ）によって与えられる所定のフィルタ処理方向に沿ってフィルタ処理することによりフィルタ補正逆投影アルゴリズムに基づいて得られ、ここで、ｇａｍｍａはコーン角であり、ｚはＣｓｅｇからの垂直方向の距離であり、Ｃｓｅｇは、検出器列数−１を２で除算することによって規定される。ステップＳ２１０では、実測された投影データは、所定のＨｃｏｎｖステップに供され、このステップでは、２つの参照画像を再構成するためにコンボリューションデータを出力する。一実施形態では、Ｈｃｏｎｖステップでは、Ｒａｍｐカーネル＋Ｈｉｌｂｅｒｔカーネルのハイブリッドカーネルを利用する。別の実施形態では、２つの参照画像を生成するために、円形フェルドカンプ（ＣＦＫ）法が用いられる。さらに、頭部イメージングのための脳のシェーディングなどにおいて画質をある程度向上させるために、任意選択でリビニングステップがＨｃｏｎｖステップに追加される。代替実施形態では、コンボリューションデータのリビニングステップおよび逆リビニング（inverse rebinning）はそれぞれ、画質をある程度向上させるためにＨｃｏｎｖステップの前と後に任意選択で実行される。ステップＳ３００では、第１の参照画像ＣＦＫ＿Ａは、円形フェルドカンプ（ＣＦＫ）法を用いて円軌道コーンビームデータから再構成され、第１の参照ボリューム画像ＣＦＫ＿Ａは第１の撮像視野（ＦＦＯＶ）を有する。ＦＦＯＶは、一般に、所望の逆投影視野（back projection field of view）である。同様に、ステップＳ２２０では、第２の参照画像ＣＦＫ＿Ｂが円軌道コーンビームデータから再構成され、第２の参照ボリューム画像ＣＦＫ＿Ｂは第２の撮像視野（ＳＦＯＶ）を有し、ＳＦＯＶはＦＦＯＶより大きく、任意選択でＣＴ画像システムのガントリを対象とする。次に、両方の参照画像ＣＦＫ＿ＡおよびＣＦＫ＿Ｂは、任意選択で、後で取り出すためにそれぞれのステップＳ３００およびＳ２２０において記憶される。そのうえ、ステップＳ２２０ではまた、反復に関与するステップのその後のインスタンスのために反復カウンタＩｔｅｒを初期化する。

さらに図５を参照して、コーンビームアーチファクトを大幅に軽減する処理の反復ステップについて詳細にさらに説明する。ステップＳ２３０では、所定数の反復を追跡するために、反復カウンタＩｔｅｒが１だけ増分される。ステップＳ２４０では、ＯＳＲフィルタなどの所定のフィルタが、第２の参照ボリューム画像ＣＦＫ＿Ｂに適用され、画像ＣＦＫ＿Ｃを生成するが、この画像ＣＦＫ＿Ｃは、ＬＬＦＤＫ補正されたＣＦＫ画像である。その後、次にステップＳ２５０において、画像ＣＦＫ＿Ｃが、順投影のために拡張スライスＣＦＫ＿Ｄにボリューム拡張される（volume expanded）。ステップＳ２６０では、スライスＣＦＫ＿Ｄが線データに順投影され、ステップＳ３１０およびステップＳ２７０において、同じ線データがそれぞれ線画像Ａおよび線画像Ｂに再構成される。その後、線画像Ｂが、ステップＳ２７０において反復ごとに線画像スケーリングファクタによりスケーリングされ、ステップＳ２８０において第２の参照画像ＣＦＫ＿Ｂと合成されて、ステップＳ２９０でＬＣ画像Ｂを生成してから、ステップＳ２３０においてＬＣ画像Ｂを用いて次の反復を開始し、このステップで反復カウンタＩｔｅｒが増分される。

同様に、線画像Ａも、Ｓ３１０において線画像スケーリングファクタによりスケーリングされ、ステップＳ３２０において第１の参照画像ＣＦＫ＿Ａと合成されて、ステップＳ３３０においてＬＣ画像Ａを生成する。ステップＳ３４０において、反復カウンタＩｔｅｒの値が所定の反復回数Ｎｉｔｅｒより大きいかどうかが判断される。ステップＳ３４０でＩｔｅｒカウンタの値が最大反復値Ｎｉｔｅｒより大きくないと判断された場合、所定の反復回数はまだ完了しておらず、コーンビームアーチファクトを大幅に軽減する処理は、さらなる反復のためにステップＳ２６０に進む。一方、ステップＳ３４０でＩｔｅｒカウンタ値が最大反復値Ｎｉｔｅｒより大きいと判断された場合、所定の反復回数は完了しており、コーンビームアーチファクトを大幅に軽減する処理は、任意選択で、ステップＳ３５０でＬＬＦＤＫ補正を適用してから、ステップＳ３６０で補正された最終画像を生成する。本実施形態によれば、コーンビームアーチファクトを大幅に軽減する処理は、上記の説明に限定されず、他の実施ステップを含む。

図６Ａおよび図６Ｂはそれぞれ、本実施形態における大幅なコーンビームアーチファクト軽減の効果を評価するために用いられるディスクファントムの側面図と上面図を示す。このファントムは、Ｚ方向に沿ってディスク状の割れ目を有する。

図７Ａ、図７Ｂ、および図７Ｃは、本実施形態よる一比較例における大幅なコーンビームアーチファクト軽減の効果を示す。図７Ａは、コーンビームアーチファクト補正を行わずにディスクファントムから再構成された従来のＦＤＫ画像を示す。著しいコーンビームアーチファクトが認められる。図７Ｂは、コーンビームアーチファクト補正を行わない５００回反復後のディスクファントムから再構成されたＳＡＲＴ画像を示す。図７Ｃは、本実施形態によりコーンビームアーチファクトを大幅に軽減するための１つの例示的処理を行ったディスクファントムから再構成された補正された画像を示す。著しいコーンビームアーチファクトは大幅に軽減されている。

図８Ａから図８Ｆは、本実施形態による別の比較例における大幅なコーンビームアーチファクト軽減の効果を示す。図８Ａから図８Ｆにおける例示的な画像は、臨床で得られた頭部データの厚いＭＰＲスライスからＣＦＫ再構成法またはＨＦＫ再構成法を用いて生成されたものである。そのうえ、図８Ａから図８Ｆにおける例示的な画像は、線補正、データのリビニング、およびＯＳＲ補正の組み合わせを用いて生成されている。さらに、図８Ａは、リビニングを行わずにＣＦＫを用いた再構成画像を示す。図８Ｂは、リビニングは行わないがＯＳＲフィルタ処理を行ったＣＦＫを用いた再構成画像を示す。図８Ｃは、比較のためにＨＦＫを用いた再構成画像を示す。図８Ｄは、リビニングを行ったＣＦＫを用いた再構成画像を示す。図８Ｅは、本実施形態のコーンビームアーチファクトを大幅に軽減する処理によりリビニングと１回順投影された線データ補正とを行ったＣＦＫを用いた再構成画像を示す。図８Ｆは、本実施形態のコーンビームアーチファクトを大幅に軽減する処理によりリビニングと３回反復順投影された線データ補正とを行ったＣＦＫを用いた再構成画像を示す。

しかし、本実施形態の多数の特徴および利点が構造および機能の詳細と共に前述の説明に記載されているが、本開示は例示に過ぎないこと、さらに、詳細、特に部品の形状、大きさ、および構成、ならびにソフトウェア、ハードウェア、またはその両方の組み合わせにおける実装形態に関して、変更を加えることができるが、この変更は、添付の特許請求の範囲が表現される用語の広い一般的な意味によって最大限示される本発明の原理内に含まれることを理解されたい。

Claims (26)

1. 円軌道コンピュータ断層撮影においてコーンビームアーチファクトを補正する方法であって、
   ａ）所定の円形線源軌道における実測の投影データから参照画像を再構成するステップと、
   ｂ）前記円形線源軌道に対応する理論的な完全円形軌道に従って前記参照画像を順投影することにより投影データを生成するステップと、
   ｃ）前記生成された投影データから補正画像を再構成するステップと、
   ｄ）前記参照画像と前記補正画像とを用いて最終画像を生成することによって前記コーンビームアーチファクトを大幅に軽減するステップと
   を備える方法。
2. 前記参照画像が、フィルタ補正逆投影アルゴリズムに基づいて所定のフィルタ処理方向に沿ってフィルタ処理することによって得られる、請求項１に記載の円軌道コンピュータ断層撮影においてコーンビームアーチファクトを補正する方法。
3. 前記所定のフィルタ処理方向がＣｓｅｇ＋ｚ／ｃｏｓ（ｇａｍｍａ）によって与えられ、ここでｇａｍｍａはコーン角であり、ｚはＣｓｅｇからの垂直方向の距離であり、Ｃｓｅｇは、検出器列の数−１を２で除算することによって規定される、請求項２に記載の円軌道コンピュータ断層撮影においてコーンビームアーチファクトを補正する方法。
4. 前記参照画像が、円形フェルドカンプ（ＣＦＫ）画像とハイブリッドコンボリューション再構成（Ｈｃｏｎｖ）画像のうちの１つである、請求項２に記載の円軌道コンピュータ断層撮影においてコーンビームアーチファクトを補正する方法。
5. 前記補正画像が線スキャン再構成画像である、請求項１に記載の円軌道コンピュータ断層撮影においてコーンビームアーチファクトを補正する方法。
6. 前記補正された画像を前記ステップｂ）における前記参照画像として用いて前記ステップｂ）からｄ）を所定の回数反復する追加のステップｆ）をさらに備える、請求項１に記載の円軌道コンピュータ断層撮影においてコーンビームアーチファクトを補正する方法。
7. 前記ステップｂ）においてローパスフィルタが適用される、請求項１に記載の円軌道コンピュータ断層撮影においてコーンビームアーチファクトを補正する方法。
8. 前記ステップａ）において、前記計測円形データから第１の円形画像と第２の円形画像とを生成し、前記第１の円形画像が全撮像視野を有し、前記第２の円形画像が所望の撮像視野である、請求項１に記載の円軌道コンピュータ断層撮影においてコーンビームアーチファクトを補正する方法。
9. 前記ステップｂ）において、前記第１の円形画像に所定のＯＳＲフィルタを適用する、請求項８に記載の円軌道コンピュータ断層撮影においてコーンビームアーチファクトを補正する方法。
10. 前記ステップｂ）において、順投影された線データを生成するために前記第１の円形画像を順投影する、請求項９に記載の円軌道コンピュータ断層撮影においてコーンビームアーチファクトを補正する方法。
11. 前記ステップｃ）において、前記順投影された線データから線画像を再構成する、請求項１０に記載の円軌道コンピュータ断層撮影においてコーンビームアーチファクトを補正する方法。
12. 前記ステップｄ）において、前記第２の円形画像と前記線画像とを用いて前記補正された画像を生成することによって前記コーンビームアーチファクトを大幅に軽減する、請求項１１に記載の円軌道コンピュータ断層撮影においてコーンビームアーチファクトを補正する方法。
13. 前記補正された画像を前記ステップｂ）における前記参照画像として用いて前記ステップｂ）からｄ）を所定の回数反復する追加のステップｆ）をさらに備える、請求項１２に記載の円軌道コンピュータ断層撮影においてコーンビームアーチファクトを補正する方法。
14. 円軌道コンピュータ断層撮影においてコーンビームアーチファクトを補正するシステムであって、
    所定の円形線源軌道における実測された投影データから参照画像を再構成するための再構成デバイスと、
    前記円形線源軌道に対応する理論的な完全円形軌道に従って前記参照画像を順投影することにより投影データを生成するための、前記再構成デバイスに接続されたコーンビームアーチファクト軽減デバイスとを備え、
    前記再構成デバイスは、前記生成された投影データから補正画像を再構成し、
    前記コーンビームアーチファクト軽減デバイスは、前記参照画像と前記補正画像とを用いて最終画像を生成することによって前記コーンビームアーチファクトを大幅に軽減する、システム。
15. 前記参照画像が、フィルタ補正逆投影アルゴリズムに基づいて所定のフィルタ処理方向に沿ってフィルタ処理することによって得られる、請求項１４に記載の円軌道コンピュータ断層撮影においてコーンビームアーチファクトを補正するシステム。
16. 前記所定のフィルタ処理方向がＣｓｅｇ＋ｚ／ｃｏｓ（ｇａｍｍａ）によって与えられ、ここでｇａｍｍａはコーン角であり、ｚはＣｓｅｇからの垂直方向の距離であり、Ｃｓｅｇは、検出器列の数−１を２で除算することによって規定される、請求項１５に記載の円軌道コンピュータ断層撮影においてコーンビームアーチファクトを補正するシステム。
17. 前記参照画像が、円形フェルドカンプ（ＣＦＫ）画像とハイブリッドコンボリューション再構成（Ｈｃｏｎｖ）画像のうちの１つである、請求項１５に記載の円軌道コンピュータ断層撮影においてコーンビームアーチファクトを補正するシステム。
18. 前記補正画像が線スキャン再構成画像である、請求項１４に記載の円軌道コンピュータ断層撮影においてコーンビームアーチファクトを補正するシステム。
19. 前記再構成デバイスおよび前記コーンビームアーチファクト軽減デバイスが、前記補正された画像を前記参照画像として用いて所定の反復回数反復する、請求項１４に記載の円軌道コンピュータ断層撮影においてコーンビームアーチファクトを補正するシステム。
20. 前記再構成デバイスが、ローパスフィルタを利用する、請求項１４に記載の円軌道コンピュータ断層撮影においてコーンビームアーチファクトを補正するシステム。
21. 前記再構成デバイスが、前記計測円形データから第１の円形画像と第２の円形画像とを生成し、前記第１の円形画像が全撮像視野を有し、前記第２の円形画像が所望の撮像視野である、請求項１４に記載の円軌道コンピュータ断層撮影においてコーンビームアーチファクトを補正するシステム。
22. 前記コーンビームアーチファクト軽減デバイスが、前記第１の円形画像に所定のＯＳＲフィルタを適用する、請求項２１に記載の円軌道コンピュータ断層撮影においてコーンビームアーチファクトを補正するシステム。
23. 前記コーンビームアーチファクト軽減デバイスが、順投影された線データを生成するために前記第１の円形画像を順投影する、請求項２２に記載の円軌道コンピュータ断層撮影においてコーンビームアーチファクトを補正するシステム。
24. 前記再構成デバイスが、前記順投影された線データから線画像を再構成する、請求項２３に記載の円軌道コンピュータ断層撮影においてコーンビームアーチファクトを補正するシステム。
25. 前記コーンビームアーチファクト軽減デバイスが、前記第２の円形画像と前記線画像とを用いて前記補正された画像を生成することによって前記コーンビームアーチファクトを大幅に軽減する、請求項２４に記載の円軌道コンピュータ断層撮影においてコーンビームアーチファクトを補正するシステム。
26. 前記再構成デバイスおよび前記コーンビームアーチファクト軽減デバイスが、前記補正された画像を前記参照画像として用いて所定の反復回数反復する、請求項２５に記載の円軌道コンピュータ断層撮影においてコーンビームアーチファクトを補正するシステム。

Images