JP3976684B2 - 画像における動きの影響を低減する方法および装置 - Google Patents

Description

本発明は、画像処理の分野に関し、特定的には、特に映されている対象の動きにより、獲得中画像に取り込まれるアーチファクトが、出力画像の生成において補正される、画像処理技術に関する。

磁気共鳴画像法（ＭＲＩ）は、医療診断画像法に幅広く使用される技術である。従来のＭＲＩスキャナでは、患者は、場の方向に対して平行か、反平行かのいずれかである、非ゼロスピン量子数を伴う核の磁気モーメントの整列を生じる、強度静磁場中に配置される。２の向きの間でのモーメントのボルツマン分布は、場の方向に沿った正味の磁化を生じる。この磁化は、検査中の核種および加えられる場の強度によって決定される周波数の無線周波数（ＲＦ）磁場をかけることによって操作されることができる。医療用画像法のほぼすべての場合において、検査される種は水素原子の核であり、主に水分子で体内に存在し、およびＲＦパルスは、これら水プロトンの共鳴周波数でかけられる。しかしながら、他の核種および他の種に結合される水素を映すことも可能である。

ＲＦ場からの核によって吸収されるエネルギはその後再放出され、および適切に調整されたアンテナで、発振電圧として、または自由誘導減衰信号として検出されてもよい。より一般的には、さらなるＲＦパルスまたは磁場勾配が、信号獲得を遅らせ、およびスピン‐エコーあるいは勾配エコー‐エコー信号を生成するために使用される。
空間情報は、エコー獲得の間に、またはエコー獲得の前にかけられる、磁場勾配として知られる追加的な線形に変化する磁場によって、エコー信号に符号化される。空間符号化の原理は、場勾配がある場合は、一定の核モーメント、そしてその共鳴周波数によって経験される正味の場は、スキャナとの位置の関数である。勾配がエコー獲得中にかけられるとき、受信された信号は、勾配方向に沿った様々な場所での核を表すある範囲の周波数成分を含む。この信号のフーリエ変換は、患者を通した１次元投影を作り出す。この技術は、周波数符号化として公知である。２次元符号化は、位相符号化勾配として知られる、周波数符号化軸に垂直にかけられる、追加的勾配の使用を必要とする。この勾配は、データ獲得前の短い時間、かけられる。獲得プロセスは、異なる強度の位相符号化勾配を用いて、何回も、おそらくは２５６または５１２回繰り返される。周波数および位相同時符号化は、２次元フーリエ変換を受けると所望の画像を提供する、２次元データ・セットを産出する。このデータのアレイは、ｋ−空間として知られるものに存在し、および画像空間のフーリエ変換である。位相符号化勾配の効果は、データ獲得の開始を（勾配強度に依存した）ｋ−空間の１つの軸に沿った特定の場所に移動することであり、周波数符号化は、他の軸に平行にｋ−空間を通過するスイープを表す。これらのスイープの各々は、“ショット”や“ビュー”として知られる。

第３次元での空間的位置限定は、追加の位相符号化勾配を用いて、又は、より一般的には、勾配および狭帯域ＲＦパルスを用いる核モーメントの初期摂動を単一の断層Ｘ線写真のスライスに制約することによって達成されることができる。この原理は、マルチスライスＭＲＩに容易に拡張することができる。
従来のＭＲＩにおいては、単一の位相符号化ビューは、各ＲＦ励起の後に獲得される。しかしながら、各励起の後に、一連の、異なって符号化されたエコーを獲得するために、さらなるＲＦパルスおよび位相符号化勾配が使用される、より早い画像化シーケンスが現在存在する。これらのエコーは、ｋ−空間のいくつかの線を横切り、および一連のエコーの長さに等しい係数だけ走査時間を低減する。極端な場合には、単一のショット・エコー・平面画像（ＥＰＩ）技術が、１００ｍｓより短く続く単一の獲得において、２次元ｋ−空間の全体を網羅するが、空間解像度および画像品質は、大幅に譲歩されるかもしれない。

ＭＲＩ画像の獲得中の患者の動きは、臨床的に関連のある情報をあいまいにしうる画像の劣化を生じる。各読み出し期間は、数ミリ秒かかる一方で、読み取りの間の時間間隔は１００から４０００ｍｓの間でもよい。患者の動きによって生じるぼけおよびゴースティング・アーチファクトの大半は、単一の読み出しの間の動きよりもむしろ、ｋ−空間における線の間の動きによる。
動きは、その結果の画像に、位相符号化方向に沿ったぼけやゴースティングとして現れる、ｋ−空間の線の間の位相エラーにつながる。これらの位相エラーは、並進動きおよび回転動きから生じうる。読み出し方向での患者の並進は、ｋ−空間の各線における周波数依存位相シフトを生じる。空間領域における回転は、ｋ−空間での回転でもあり、およびｋ−空間における位置のより複雑な関数であるｋ−空間変化を生じる。
動きを通して画像に取り込まれる画像アーチファクトを補正しようと試みるため、様々な技術が採用されてきた。しかしながら、患者の動きを補正するとして知られている技術のほとんどは、追加の走査または追加の設備さえも伴う場合もある、修正された信号獲得技術を伴う。

国際特許出願ＷＯ９８／０１８２８は、純粋にポストデータ収集信号処理効果を用いて、画像において動きが誘発したアーチファクトの影響を低減するための技術を開示する。そこに説明されている技術において、可能性のある動きが誘発するアーチファクトに対抗するためにデータを操作し、および操作されたデータは、画像品質が向上しているかどうかを見るために、焦点条件を用いて比較される。この技術は、しかしながら、大量の処理を伴いうる。さらに、該方法は、より正確に動きのパラメータを決定するために、ｋ−空間線のグループ化を伴ってもよいが、このグループ化は、見つかった動きの一時的な解像度を下げうる。
動きが誘発した画像アーチファクトを補正するための別の方法は、コンプレックス・セットに対する投影法（ＰＯＣＳ）、Ｈｅｄｌｅｙ Ｍ、Ｈｏｎｇ Ｙ ａｎｄ Ｒｏｓｅｎｆｅｌｄ Ｄ．“Ｍｏｔｉｏｎ Ａｒｔｉｆａｃｔ Ｃｏｒｒｅｃｔｉｏｎ ｉｎ ＭＲＩ ｕｓｉｎｇ ｇｅｎｅｒａｌｉｚｅｄ ｐｒｏｊｅｃｔｉｏｎｓ”ＩＥＥＥ Ｔｒａｎｓ．Ｍｅｄ．Ｉｍａｇ．，１０：４０−４６，１９９１である。これは、それによって、バイナリ・マスクを形成するために、上質の画像が使用される方法である。該マスクは、組織‐空気境界を定義し、すなわち、マスクの外側には信号があってはならない。獲得された画像における動きが誘発したアーチファクトは、空気中で見かけ上の信号を生じる。ＰＯＣＳ法は獲得された画像においてマスクの外側すべてを黒に設定する。画像データはそれからｋ−空間にフーリエ変換される。新しい複合ｋ−空間は、測定されたデータの絶対値から、および前のステップからの推定値の位相から形成される。この新しいｋ−空間は、画像領域にフーリエ変換され、該処理は繰り返される。この方法は、しかしながら、処理が繰り返されると大量のフーリエ変換を伴うので、大量の計算量および従って時間を伴う。さらに、該方法は、常に達成されるとは限らない処理の前に、獲得された画像と、バイナリ・マスクとの空間整列を求める。

拡散重みつけされたＭＲＩにおいて動きによって誘発されたアーチファクトを補正するための異なる方法が、‘Ｃｏｒｒｅｃｔｉｏｎ ｏｆ Ｍｏｔｉｏｎａｌ Ａｒｔｉｆａｃｔｓ ｉｎ Ｄｉｆｆｕｓｉｏｎ−Ｗｅｉｇｈｔｅｄ Ｉｍａｇｅｓ Ｕｓｉｎｇ ａ Ｒｅｆｅｒｅｎｃｅ Ｐｈａｓｅ Ｍａｐ’、Ｕｌｕｇ ｅｔ ａｌ． Ｍａｇｎｅｔｉｃ Ｒｅｓｏｎａｎｃｅ ｉｎ Ｍｅｄｉｃｉｎｅ、Ｖｏｌ．３４，ｎｏ．２，１９９５，ｐｐ４７６−４８０で説明されている。そこで説明された方法は、拡散重みつけのない、または大変低い重みつけで測定し、および位相マップを構築するためにこの情報を使用する。この位相マップからの位相は、それから動きによって誘発される位相成分を排除し、かつ調節されたデータ・セット形成するために、より高い重みをかけられたデータ・セットの位相を置き換えるために使用される。

Ｕｌｕｇらによって説明される方法は、動きに対して大変敏感であり、すべてのＭＲＩ画像には適用できない、拡散重みつけＭＲＩ走査に関する。さらに、説明された該方法は、映し出されている間の対象の回転を補正することはできないだろう。
本発明は、画像における、動きによって誘発されたアーチファクトを補正するための代替的な方法を提供しようと試みる。

本発明に従って、次のものが提供される
走査中に、対象の望ましくない動きによって取り込まれるアーチファクトが補正された、走査された対象の画像を作る方法であって；
該対象の走査から導出された走査ｋ−空間データ・セットを取り、
該対象の参照ｋ−空間データ・セットを取り、
該走査k-空間データ・セットにおける走査中の何らかの動きの程度を決定するために、該走査ｋ−空間データ・セットと参照ｋ−空間データ・セットとを比較し、
動きを補正された画像セットを作るために、該走査ｋ−空間データ・セットおよび該走査ｋ−空間データ・セットの動きの程度に関する情報を使用し、および
画像を作るために、該画像を補正された走査ｋ−空間データ・セットを使用するステップ
を含む方法。

参照ｋ−空間データ・セットは、対象の画像に対応する、動きの無い（または該環境においてできるだけ動きのないものに近い）ｋ−空間データである。動きの無いデータ・セットを取ることによって、崩れた画像データ・セットにつながる動きがモデル化されうる。走査ｋ−空間データ・セットのみからのデータは、モデル化された動きに基づいて補正された画像を作るために使用される。２のデータ・セットがｋ−空間または空間周波数領域で比較されるので、データ・セット全体の反復したフーリエ変換は必要ない。該方法が、決定された動きに関して走査ｋ−空間データ・セットのデータを補正し、それを他のデータとおきかえるのではないので、実際の画像からの情報が保存される。
該方法は好ましくは、走査ｋ−空間ならびに参照ｋ−空間データ・セットをｋ−空間におけるセクションに分割するステップを含み、および動きを決定するために、走査ｋ−空間データ・セットと参照データ・セットを比較するステップは、セクションごとにデータ・セットを比較する。２Ｄ画像に関して、各セクションは都合よく、スキャナによって取られた位相符号データの線でありうる。３Ｄ画像に関しては、各セクションはｋ−空間データの平面である。

２のデータ・セットを比較するステップは、何らかの回転の程度を決定するためにデータ・セットを比較するステップと、および何らかの並進の程度を決定するためにデータ・セットを比較するステップを含むと都合がよい。画像において拡張や収縮を生じさせる動きも、そのような方法によって補正されてもよい。
好ましくは、データ・セットを比較するステップは、何らかの回転動きの程度を決定するために、データ・セットにおける選択された点を比較するステップを含む。選択された点のｋ−空間絶対値が、比較のために使用されると都合よい。ｋ−空間データの絶対値は位相から独立しており、よって並進動きによって影響されない。映し出されている対象の回転はｋ−空間において回転を生じさせる一方で、並進は、ｋ−空間データへ位相傾斜を与える。点の選択は、ｋ−空間データの、ブロックへの分離を伴ってもよい。ブロック内の一定の点が、例えば各ブロックにおける最大の絶対値値を持つ点を選択することによって、選択されうる。ブロックへの分離は、ｋ−空間の合理的な範囲が達成され、十分な角解像度を与えることを確かにする。

回転動きの程度は、参照ｋ−空間データにおける選択された点の可能な角度動きを繰り返し検討することにより、および回転動きの程度を決定するために、適切な焦点基準を用いて、走査ｋ−空間データ・セットと変換された参照データ・セットとを比較することによって決定されうる。参照ｋ−空間データ・セットを変換するステップは、好ましくは、線形補間を用いて実行されるが、３次またはスプラインなど他の補間が使用されうる。
並進動きの程度を決定するためにデータ・セットを比較するステップは、Ｒｅｄｄｙ ＢＳ ａｎｄ Ｃｈａｔｔｅｒｊｉ ＢＮ、“Ａｎ ＦＦＴ−Ｂａｓｅｄ Ｔｅｃｈｎｉｑｕｅ ｆｏｒ Ｔｒａｎｓｌａｔｉｏｎ， Ｒｏｔａｔｉｏｎ ａｎｄ Ｓｃａｌｅ−Ｉｎｖａｒｉａｎｔ Ｉｍａｇｅ Ｒｅｇｉｓｔｒａｔｉｏｎ”、ＩＥＥＥ Ｔｒａｎｓ．Ｉｍａｇ．Ｐｒｏｃ．，５：１２６６−１２７１，１９９６に説明されたような、クロス・スペクトル方法に基づいた方法を用いて決定されるのが都合がよい。クロス・スペクトル解析を実行する際に、ｋ−空間は、解像度を改善するために、ゼロパッドされてもよく、または１に等しい絶対値および任意の位相を有する複素数値でパッドされてもよい。クロス・スペクトル方法は、２のデータ・セットの間の定数および線形に変化する位相差を検出する。２のデータ・セットは、後述されるように、関数Ｋ_cを形成するために結合される。Ｋ_cのフーリエ変換は、ピークを与える。中心からのピークのシフトは位相スロープと関連しており、ＦＥ方向への動きの量を示す。ピークの位相は、２つのデータ・セットの間の位相差を明らかにし、およびＰＥ方向の動きの程度を明らかにする。

並進動きの程度を決定した後、走査ｋ−空間データ・セットは、決定された動きに関して補正されてもよい。この補正されたデータ・セットの位相は、周波数符号化方向へのサブピクセルの動きを決定するために、参照ｋ−空間データの位相と比較されうる。
何らかの点が、参照ｋ−空間データ・セットまたは走査ｋ−空間データ・セットのいずれかで信頼性がないことが分かっている場合、Ｋ_cにおける対応する点は信頼性がない。そのような点は、動きの決定に影響を与える信頼性のない点を避けるため、１に等しい絶対値および任意の位相を有する点と置き換えられうる。
クロス・スペクトルに基づいた方法は、確実性の測定を含んでもよい。理想的な場合では、Ｋ_cのフーリエ変換の結果は、単一のピークになるだろう。現実には、ノイズ、異なる獲得基準等による参照および走査ｋ−空間データ・セットの間の差が、影響を及ぼす。これは、Ｋ_cのＦＴが、非ゼロである他の値を伴うピークとなることを生じる。他の値の平均値に対するピークの高さの比は、信頼性の測定値として使用されることができ、この比が高いほど、動きの決定はより信頼性がある。

参照データ・セットおよび画像データ・セットを取るステップは、データ・グリッドの中心に一致して最大（ＤＣ）点がくるように、走査ｋ−空間データ・セットと参照ｋ−空間参照データ・セットの両方を整列するステップを含んでもよい。さらに、該方法は、最大（ＤＣ）点において、２のデータ・セットが同じ複素数値を有するように、２のデータ・セットを全体的にスケール調整するステップを含みうる。
最大絶対値を伴うｋ−空間点は、０空間周波数に対応する点であると仮定されることが注記されるべきである。０空間周波数はときにＤＣと呼ばれる。
さらに、該方法は、一定の量のぼけを取り込むために、獲得された画像と参照画像の両方をフィルタにかけるステップを含んでもよい。画像にフィルタをかけることで対応するｋ−空間をぼかし、および線形補間を用いて回転動きの決定を手助けしうる。ｋ−空間をぼかすことは、それがあまり迅速に変化せず、よって点の間の補間がより簡単であることを意味する。
走査ｋ−空間データ・セットは、磁気共鳴画像（ＭＲＩ）走査手段から得られたデータ・セットでもよい。

参照ｋ−空間データ・セットは、多くの方法で得られうる。多くの走査プロトコルは、信号対ノイズ比をより良くするために、２またはそれ以上のデータ・セットおよび平均を取得する。そのような場合、データ・セットの１は、その他における何らかの動きを補正するために、参照画像として取られうる。最善のデータ・セットは、画像のエントロピの測定値、すなわち画像における乱れの程度を用いることによって決定されうる。
あるいは、対象が、造影剤の導入で連続的に映し出されている場合、１つの画像は参照ｋ−空間データ・セットとして使用されうる。
あるいは、参照データ・セットは以前に獲得された走査から取られうる。例えば進行性の疾患を監視する際、異なる時に同じ被写体への繰り返した走査が必要となりうる。本発明は、参照データ・セットからのｋ−空間情報を用いるのではなく、走査ｋ−空間データ・セットに関する何らかの動きの程度を決定するだけなので、前の画像が使用されることができ、および該技術は、画像の間での差、例えば様々なコントラスト、ノイズ、または増進する腫瘍などの特徴に対しても合理的に強くあるべきである。使用できる参照データが無い場合、参照画像は、使用可能なデータから作成されうる。

画像全体を補正するよりは、画像の領域は、局限された領域のみに補正をしてマスクされうる。この場合、該方法は、画像の所望の部分を選択し、および走査ｋ−空間データ・セットを作るために、画像をフーリエ変換する前に、画像の所望でない範囲を削除するステップを含む。
本発明の別の特徴においては、対象を走査し、ならびに該対象から受信された信号から走査ｋ−空間データ・セットを導出する走査手段と、および該走査中に、前記対象の望ましくない動きによって取り込まれるアーチファクトが補正された、走査された対象の画像を作るために、走査ｋ−空間データ・セットを操作する処理手段を含む画像システムであって、該処理手段は、対象走査から導出された走査ｋ−空間データ・セットを取り、対象の参照ｋ−空間データ・セットを取り、画像データ・セットにおける何らかの動きの程度を決定するために走査ｋ−空間データ・セットと参照ｋ−空間データ・セットとを比較し、動きを補正された走査ｋ−空間データ・セットを作るために走査ｋ−空間データ・セットの動きの程度に関する情報を使用し、および画像を作るために、動き補正された走査データ・セットを使用するように適応しているシステムが提供される。
本発明は、以下の図面を参照して、例示的にのみ説明される；

図１を参照すると、磁気共鳴画像システム１０の概略図が示されている。システム１０は、従来の種類の磁気共鳴画像スキャナ１２を含んでいる。スキャナ１２は、患者の中の原子核の場方向に沿った正味の整列を生じさせるだけ十分に強い磁場を生成する超伝導または抵抗性主磁石２０を有する。スキャナ１２は、主磁石２０の磁場における望ましくない不均一性を補正するために、シム（ｓｈｉｍ）コイル２２も含む。シムコイル２２によって作られる磁場は、シムコイル電力供給ユニット２４によって制御される。
特定の原子核の共鳴周波数は、核およびかけられる磁場の強度に特有である。空間情報を提供するために、磁場勾配は、コイル２６などの勾配コイルによって生成される。勾配コイルは、３直交方向に勾配場を生成するようにしばしば配置される。勾配コイルによって生成される磁場は、勾配コイル電力供給ユニット２８によって制御される。患者の原子核から信号を生成するために、無線周波数磁気パルスが、送信コイル３０によって生成される。このパルスは、一定の患者のスライスの体積の中で、核スピンの角度を‘フリップ’する。これらの励起されたスピンまたは磁化はそれから、送信コイル３０と同じコイルでもよい受信コイルに電流を誘導する。コイル３０は、送信ユニット３２および受信ユニット３４に接続され、その各々は、周波数源３６からも信号を受信する。

システム１０は、システム１０の構成要素の操作を制御する制御コンピュータ３８を含む。コンピュータ３８は、勾配タイミング、磁場強度、および向き制御の形式で、勾配コイル電力供給ユニット２８を制御する。加えて、該コンピュータは、送信機タイミングと一緒に受信ユニット３４から信号を受信する。
患者の器官の画像を形成するために、患者がシステム１０へと挿入され、および一連の測定値が、静的および／または可変的勾配場の様々な組み合わせを用いて取られる。患者の組織からの信号は組織の属性、磁場勾配強度、勾配の向き、およびかけられる無線周波数パルスに関連するタイミングに依存する。可変的勾配は、受信された信号の位相、周波数、および強度を符号化する。時間の関数として受信された信号は、その後の処理のためにコンピュータ３８のメモリに記憶されるきちんとした順序のセットを形成する。
その後の信号処理段階において、受信された信号のきちんとした順序のセットに関してフーリエ変換が実行されてもよく、該変換の絶対値は、画像を形成するため、グレイ・スケールに信号を割り当てるために使用される。受信された信号のセットは、ｋ−空間に存在すると言われている。

従来のＭＲＩにおいては、患者がデータの獲得中に動くと、受信される信号が影響を受け、ｋ−空間信号の部分が崩れる。画像が再構築される方法によって、この動きは画像全体に影響を与え、最終画像においてぼけおよび／またはゴースティング・アーチファクトを生じさせる。
ここで図２を参照すると、システム１０の操作の機能的ブロック図が示されている。コンピュータ３８は、スキャナ１２からの情報を制御し、かつ受信し、およびディスプレイ５０で画像を生成するためのこの情報を使用する。この画像は、初期再構築画像である。システム１０のオペレータが、初期画像が崩れていると判断する場合、追加の信号処理ルーチンが選択される。あるいは、さらなる信号処理が自動的に発生しうる。いずれの場合も、記憶された画像データは、患者の動きの影響を低減するために処理される。

画像を処理する該方法は、本発明に従った画像処理ルーチンにおける段階の流れ図１００を示す図３を参照して、次の通りである。
獲得時間はノードに分割される。２Ｄ画像に関して、ノードは位相符号化線に対応する。３Ｄ画像に関しては、ノードはｋ−空間データの平面である。
該システムの軸は、次のように定義される；Ｘは読み出し、または周波数符号化方向である。２Ｄ画像に関しては、Ｙは位相符号化勾配方向であり、およびＺはスルー・スライス方向である。３Ｄ画像においては、Ｙは遅い位相勾配符号化方向であり、およびＺは速い位相符号化勾配方向である。よって、ｋ_xが最も迅速に変化し、それからｋ_z、そしてｋ_yにおいて最も遅く変化するデータ点が獲得される。３Ｄ画像において、ノードは定数ｋ_yの平面である。
画像処理は、画像データ１０２に関して実行され、それは例えば、順次ｋ−空間線として得られてもよいが、他の獲得技術が可能である。参照画像データ１０４も使用される。
最初に、ある程度のデータ・セットの前処理が、比較を助けるために実行されてもよい。この場合、第一の段階は、各セットに関する最大（ＤＣ）値ｋ−空間データ点が、両方のデータ・セットの中心ピクセルになるように、両方のｋ−空間が整列される整列段階１０６でもよい。
２のデータ・セットは、ＤＣにおける複素数値が各データ・セットにおいて同じであることを確実にするため、正規化段階１０８においてスケール調整されてもよい。これは好ましくは参照データ・セット１０４に関して実行される。

フィルタリング段階１１０は、実際の画像および参照画像をフィルタにかけるために採用される。このフィルタは、後に比較段階を助けるある程度のｋ−空間へのぼかしを取り込む。適切なフィルタ関数は、画像の端に向かって離れていく、画像中央における、値１のガウシアン（ガウス波形）をピクセルに乗じることである。使用された１つのフィルタは、ｎがｋ−空間次元に沿った点の数である、（０．３１２５ｎ）²に等しい分散を有するのガウシアンである。
そのような方法でフィルタにかけることによって、ｋ−空間はぼかされる（および画像は中央で明るく、端に向かうにつれて暗くなる）。これは後に線形補間ステップを助ける。
所望であれば、マスキング段階１１２は、関心のないピクセルが０に設定されうる場合に適用されうる。
フィルタにかけられ、かつ適切にマスクされた画像はそれから、ｋ−空間画像データ・セットおよびｋ−空間参照データ・セットを作るために、フーリエ変換段階１１４を受ける。フィルタにかけることおよび隠すことは、実際の画像で実行されるので、ｋ−空間へのフーリエ変換が必要となる。

適切な角解像度が達成されることを確実にするために、ｋ−空間はそれから分割段階１１６において様々なブロックへと分割される。ｋ−空間をブロックに分割することで、アルゴリズムの速度の増加を可能にしうるが、正確さの損失を生じるかもしれず、この分割はすべてのｋ−空間点を別のブロックにすることによって省略することができる。２Ｄ画像に関して、線に沿って均等に間隔をあけた多くのブロックがあってもよい。線に２５６の点がある場合、３２ブロックを選択することによって１ブロック当たり８点を生じるが、他の組み合わせが使用されうる。３Ｄでは、平面は、平面の各次元において一定の数のブロックに分割される。
次に、ノード選択段階１１８において、各ノードが同様に選択される。各ノードに関して、各ブロックにおける最大絶対値を持つ点の座標ｒ₁は、段階１２０で決定されるが、他の選択基準が利用されうる。

選択段階１２２では対象の回転に対応する角変換Ｔ_aが選択され、次いで異なる変換が選択される。角変換Ｔ_aによって変換された、選択された点の座標に対応する参照ｋ−空間Ｋ_rにおける位置の絶対値が決定される。この段階１２４において、当業者には良く理解されるように、フィルタリング段階１１０が適用された時は、線形補間技術で十分である。
最も正確に回転動きを表す角変換Ｔ_aは、どの角変換が適切な焦点基準を最小化するかを見つけることによって決定されうる。
適切な焦点基準は；

であり、Ｋ_m（ｒ_i）は、選択された点に関する修正されていない画像データである。Ｔ_a（ｒ_i）が測定された参照ｋ−空間の中に入らないとき、それは焦点基準合計から排除されるべきである。

この焦点基準は、中心（ＤＣ）ノードに近いノードに適している。中心（ＤＣ）ノードから離れた、すなわち６０より多いノードだけ離れたノードに関して、焦点基準合計は、分子および分母が、比を取る前に別個に合計されるように変更されうる。焦点基準の正規化は、様々な絶対値を伴うｋ−空間点からの相対的寄与を変えるために変更されうる。絶対値が端より大きい中心近くで、各点は、最善の角解像度を与える同じ重みを運ぶ。しかしながら、より低い絶対値を有する輪郭に関しては、ノイズは相当な問題であり、および正規化を変更することは、ノイズと角解像度の間のトレードオフを許す。
その代わりに、１つのデータ・セットの中の補間点の代わりに、以前に回転した、および記憶されたデータ・セットからの相対線は、焦点基準の評価において使用されてもよい。
該ノードに関してすべての角変換が調べられ、および最小の焦点基準が決定されると、次のノードが試みられうる。
すべてのノードが処理されると、例えば、よく知られているように、見つかった回転をｋ−空間データ・セットに適用するために、剪断アルゴリズムが適用される回転補正段階が行われる。しかしながら、補間など、補正のための他の方法も適用されうる。

次に、再度ノード選択ステップ１３０で開始する並進が補正される。しかしながら代わりの実施態様において、特定のノードに関する並進補正は、ノードに関する回転補正が実行された直後に決定されうる。すべてのノードが回転補正を受けるのを待つ必要はない。
並進は、クロス・スペクトルに基づいた方法を用いて決定される。Ｋ_rが、一定のノードにおける参照データの線（または３Ｄでの平面）におけるすべてのｋ−空間を表しており、およびＫ_mが画像データにおいて同様である場合、Ｋ_cのフーリエ変換におけるピークの位置は、定義されたｘ方向（および３Ｄではｚ）に並進を与える。Ｋ_cは次の方程式によって与えられる；

アレイの中心からの、Ｋ_cのフーリエ変換のピークのシフトは、ｘおよびｚ方向（２Ｄではｘのみ）を与え、およびｙ方向は、ピークの位相ｐから決定されうる。ｙ置換Ｙは次の式によって与えられる；

ｌは位相符号化線の数で、中心（ＤＣ）位相符号化線においては０である。
それによってｘおよびｚ置換が見られる解像度の改善のために、Ｋ_cはゼロ・パッドされてもよい。あるいは、Ｋ_cは、絶対値１および任意の位相を有する複素数値でパッドされてもよい。後者は、Ｋ_cのフーリエ変換におけるピークを、より明確に区別させ、および正確に位置を特定させるかもしれない。

さらに、Ｋ_rまたはＫ_mのいずれかにおける何らかの点は信頼性がないことが分かっている場合、Ｋ_cにおける対応する点も信頼性はない。そのような場合、信頼性がない点は、絶対値１および任意の位相を有する点と置き換えられうる。Ｋ_cのすべての点が、Ｋ_cに存在する固有の正規化により絶対値１を有するので、絶対値１が選択される。或る点が、複数の受信コイルからのｋ−空間のデータ平均または計算など、他の処理ステップによって信頼性がないと決定される場合もある。
並進情報は、次のノードに移動する前に、該ノードに関してデータを補正するための補正段階１３４において使用されうる。

上述の、クロス・スペクトルに基づいた方法が、データ・パッドなく使用される場合、該方法は、一番近いピクセルへの、周波数符号化方向補正での並進動きを決定する。代わりの実施態様（図示されていない）において、該方法は、参照データ・セットおよび画像データ・セットを用いて、上述のように実行されうる。画像データ・セットは、決定された動きに関して補正されることができ、それからこの中間補正されたデータが参照データ・セットと比較される。残りの位相傾斜は、ｆｅ方向にサブピクセルシフトを与えるように決定されうる。位相傾斜は、Ｋ_rと補正されたＫ_mとの比の位相に関して、強い線適合アルゴリズムを用いて決定されてもよい。
このように、獲得された画像は参照画像と比較されてもよく、および獲得された画像の何らかの回転および／または並進動きは補正されうる。ｋ−空間絶対値データに関する回転比較を実行することは、何らかの並進動きに対して不変であるという効果を有する。また、画像領域における回転の中心は、ｋ−空間が常にその中心で回転するので、わかる必要はない。

この方法の結果は、参照および補正された画像が空間的に配列されることである。さらに、凸集合への投影方法（ｔｈｅ ｍｅｔｈｏｄ Ｐｒｏｊｅｃｔｉｏｎ ｏｎｔｏ Ｃｏｎｖｅｘ Ｓｅｔｓ）とは異なり、参照およびデータ画像セットの、前もっての空間的整列は必要ない。
また、参照画像におけるすべての情報は、画像データに求められる動き補正を知らせるために使用されるが、補正された画像において使用される実際のデータ点のすべては、画像データから来る。
参照画像は、前に獲得された走査でもよい。進行性の障害を持つ患者はいくつかの走査をし、および履歴走査が使用されうる。実際の参照画像データは、補正された画像において使用されるものはまったくないので、該方法は、ゆっくり進行する障害の様々な段階において、走査の間に遭遇するかもしれない画像の間のわずかな相違とともに使用されうる。
あるいは、例えば、造影剤を増加して複数の画像が取られる場合、特定の画像が参照画像として取られうる。さらに、いくつかの画像技術は、２の画像をとり、およびより良い信号対ノイズ比を達成するために平均される。これらのいずれの場合にも、動きが誘発する画像アーチファクトからの解放という点で、最善の画像が使用されうる。これは、システム・オペレータによって選択されることができ、あるいは最善の画像が、適切な画像処理技術によって選択されうる。例えば、画像の乱れの量を測定するために、画像に関してエントロピに基づいた計算が実施されうる。より低いエントロピを有する画像はより多くの秩序とより少ないアーチファクトを持つ。当業者は、適切なエントロピ計算技術に精通しているだろう。

使用可能な参照データがない場合、参照画像は、使用可能なデータから作られうる。
図４は、被写体の２Ｄスライスを示す。静止した被写体の画像は、左に示されている。この画像は、参照画像として使用された。動きで崩れた画像が中央に示されているように取られた（被写体は、走査中にうなずいた）。本発明の方法を用いて、補正された画像が図のように得られた。補正された画像は、動きで崩れた画像に対して大幅に改善されていることがわかる。
図５は、被写体を通した３直交スライスを示す。ここでも、参照画像は左に示されており、動きで崩れた画像は中央に、および本発明に従って補正された画像は右に示されている。ここでも、補正された画像において改善が見られる。

さらなる実施態様において、上述のクロス・スペクトルに基づいた方法は、動きが決定されたことについての信頼性の表示を与えるために使用されうる。参照データ・セットおよび画像データ・セットが、並進動きによってのみ異なるという理想的な条件下では、Ｋ_cのフーリエ変換は、すべての他の値が０であるピークである。実際に、両方のデータ・セットは、ノイズにより、およびそれらが異なる時間に、おそらく異なる条件下で獲得されるために異なる。よって実際には、フーリエ変換は、非０である他の値を有するピークを生じる。ピークの高さ対他の値の平均との比は、信頼性の測定値として使用されることができ、比が高くなるほど、決定された動きはより信頼性がある。これは、単純に、動きの決定における信頼性の表示を与えるために使用されることができ、または異なる参照データ・セットが試されるべきであることを示すために使用されうる。該方法は自動的に比を決定することができ、および比が一定の閾値よりも低くなる場合には信号を生成する。
本発明の他の実施態様および効果は、ここに説明された本発明から逸脱することなく、当業者に明らかになるだろう。

図１は、磁気共鳴画像システムの概略図を示す。 図２は、図１のシステムの操作の機能図を示す。 図３は、本発明に従って実行される画像補正ルーチンの流れ図である。 図４は、本発明の方法を用いた、ＭＲＩ画像、参照画像、および動き補正画像を示す。 図５は、本発明を用いた、獲得された画像、参照画像、および動き補正画像に関する直交ＭＲＩを示す。

Claims (13)

1. 走査中に、対象の望ましくない動きによって取り込まれるアーチファクトを補正された、走査された対象の画像を作る画像処理方法であって；
   該対象の走査から導出された走査ｋ−空間データ・セットを受取り、
   該対象の参照ｋ−空間データ・セットを受取り、
   該走査ｋ−空間データ・セットにおいて走査中に起こる何らかの動きの程度を決定するために、該走査ｋ−空間データ・セットと参照ｋ−空間データ・セットとを比較し、
   動きを補正されたｋ−空間データ・セットを作るために、該走査ｋ−空間データ・セットおよび該走査ｋ−空間データ・セットの動きの程度に関する情報を使用し、および
   画像を作るために、該補正された走査ｋ−空間データ・セットを使用するステップ
   を含み、
   何らかの動きの程度を決定するために、該走査ｋ−空間データ・セットと参照ｋ−空間データ・セットとを比較する前記ステップが、回転動きの程度を決定するために、該走査ｋ−空間データ・セットおよび参照ｋ−空間データ・セットにおける選択された点の絶対値を比較するステップを含み、前記選択された点が、ｋ−空間を分割する複数のブロックの内の選択されたブロック内の複数の点の内の最大絶対値を伴う点であることを特徴とする方法。
2. 該方法は、該走査ｋ−空間データ・セットおよび参照ｋ−空間データ・セットをセクションに分割するステップを含み、および該走査ｋ−空間データ・セットと参照ｋ−空間データ・セットを比較するステップは、該走査ｋ−空間データ・セットの動きの程度を決定するために、該走査ｋ−空間データ・セットと参照ｋ−空間データ・セットをセクションごとに比較するステップを含むことを特徴とする、請求項１に記載の方法。
3. 該走査ｋ−空間データ・セットの各セクションは、位相符号化データの線であることを特徴とする請求項２に記載の方法。
4. 該走査ｋ−空間データ・セットの各セクションは、ｋ−空間データの平面であることを特徴とする請求項２に記載の方法。
5. 該方法は、何らかの回転の程度を決定するために、該走査ｋ−空間データ・セットと参照ｋ−空間データ・セットとを比較し、および何らかの並進の程度を決定するために、該走査ｋ−空間データ・セットと参照ｋ−空間データ・セットを別個に比較するステップを含むことを特徴とする、請求項１から４のいずれかに記載の方法。
6. 選択された点の該絶対値を比較するステップは、様々な可能性のある回転変換を用いて該参照ｋ−空間データ・セットまたは走査ｋ−空間データ・セットのいずれかにおける該選択された点を繰り返し変換し、および他のｋ−空間データ・セットと比較するステップを含むことを特徴とする、請求項１に記載の方法。
7. クロス・スペクトルに基づく方法は、並進動きの程度を決定するため、該走査ｋ−空間データ・セットと参照ｋ−空間データ・セットとを比較するために使用されることを特徴とする請求項１から６のいずれかに記載の方法。
8. 該クロス・スペクトルに基づいた方法は、組み合わされたｋ−空間データ・セットを形成するために、該参照ｋ−空間データ・セットと走査ｋ−空間データ・セットとを組み合わせ、該組み合わせたｋ−空間データ・セットのフーリエ変換を実行し、および並進動きの程度を決定するために、フーリエ変換された組み合わせたｋ−空間データ・セットのピークの位置および位相を決定するステップを伴うことを特徴とする、請求項７に記載の方法。
9. 該方法は、決定された該並進動きに関して該走査ｋ−空間データ・セットを補正し、およびサブピクセルの動きを決定するために、動きを補正された走査ｋ−空間データ・セットの位相と、該参照ｋ−空間データ・セットの位相とを比較するステップを含むことを特徴とする、請求項７または８に記載の方法。
10. 該参照ｋ−空間データ・セットおよび走査ｋ−空間データ・セットによって表される該画像は、該ｋ−空間データ・セットの比較の前に、ｋ−空間にぼけを取り込むためにフィルタにかけられることを特徴とする、請求項１から９のいずれかに記載の方法。
11. 該走査ｋ−空間データ・セットは、磁気共鳴画像走査手段から得られることを特徴とする、請求項１から１０のいずれかに記載の方法。
12. 該走査ｋ−空間データ・セットによって表される画像の一部は、比較の前にマスクされることを特徴とする、請求項１から１１のいずれかに記載の方法。
13. 対象を走査し、および該対象から受信された信号から走査ｋ−空間データ・セットを導出する走査手段と、および該走査中に前記対象の望ましくない動きによって取り込まれるアーチファクトに関して補正された、該走査された対象の画像を作るために、該走査ｋ−空間データ・セットを操作する処理手段を含む画像システムであって、該処理手段は、該対象走査から導出された走査ｋ−空間データ・セットを取り、該対象のｋ−空間参照データ・セットを取り、該走査ｋ−空間データ・セットにおける何らかの動きの程度を決定するために、該走査ｋ−空間データ・セットと参照ｋ−空間データ・セットとを比較し、動きを補正された走査ｋ−空間データ・セットを作るために、該走査ｋ−空間データ・セットの動きの程度に関する情報を使用し、および画像を作るために、該動き補正された走査ｋ−空間データ・セットを使用するように適応しており、前記比較が、回転動きの程度を決定するために、該走査ｋ−空間データ・セットおよび参照ｋ−空間データ・セットにおける選択された点の絶対値を比較することを含み、前記選択された点が、ｋ−空間を分割する複数のブロックの内の選択されたブロック内の複数の点の内の最大絶対値を伴う点であることを特徴とする画像システム。

Images