JP2007159718A - 磁気共鳴イメージング装置および画像補正評価方法 - Google Patents

Abstract

【課題】撮影領域に形成される静磁場の誤差分布が適切に推定されて静磁場の不均一性による画像歪みの補正が適切に行われたか否かを評価することが可能な磁気共鳴イメージング装置および画像補正評価方法である。
【解決手段】磁気共鳴イメージング装置２０は、異なる複数の撮影条件を設定する撮影条件設定手段４６と、複数の撮影条件に従って撮影領域から磁気共鳴信号を受信する受信手段と、磁気共鳴信号に基づいて複数の撮影条件に対応する複数の画像データを再構成する画像再構成手段４２と、撮影領域の磁場分布に基づいて複数の画像データの歪みを補正する画像歪み補正手段４５と、補正後の複数の画像データに基づいて少なくとも１つの画像データの補正が適切に行われたか否かを評価する画像補正評価手段４７、４８とを備える。
【選択図】 図２

Description

本発明は、被検体の原子核スピンをラーモア周波数のＲＦ信号で磁気的に励起し、この励起に伴って発生する磁気共鳴信号から画像を再構成する磁気共鳴イメージング装置および画像補正評価方法に係り、特に磁場の不均一性による画像歪みの補正が適切に行われたか否かを評価することが可能な磁気共鳴イメージング装置および画像補正評価方法に関する。

従来、医療現場におけるモニタリング装置として、磁気共鳴イメージング（ＭＲＩ：Magnetic Resonance Imaging）装置が利用される。

ＭＲＩ装置は、静磁場を形成する筒状の静磁場用磁石内部にセットされた被検体の撮像領域に傾斜磁場コイルでＸ軸、Ｙ軸、Ｚ軸方向の傾斜磁場を形成するとともにＲＦ（Ｒａｄｉｏ Ｆｒｅｑｕｅｎｃｙ）コイルから高周波（ＲＦ）信号を送信することにより被検体内の原子核スピンを磁気的に共鳴させ、励起により生じた核磁気共鳴（ＮＭＲ：Nuclear Magnetic Resonance）信号を利用して被検体の画像を再構成する装置である。

このようなＭＲＩ装置では、静磁場を形成する静磁場用磁石の性能限界や磁場中に被検体が入ることによる磁化率の変動によって、局所的に静磁場強度の均一性が乱されるという問題がある。静磁場強度の均一性が乱されると、再構成された画像に歪みが生じることが知られている。

静磁場強度の不均一性の影響を特に強く受けるケースとしては、ＥＰＩ（echo planar imaging）法による撮像を行う場合が挙げられる。ＥＰＩ法は、一回の励起パルスで画像再構成のためのすべてのデータを取得する超高速スキャン法であり、シングルショットＥＰＩ法とも呼ばれる。ＥＰＩ法による撮像は、例えば被検体の頭部の検査における脳梗塞などの描出の際に、被検体の体動や拍動の影響を避けるために行われるが、ＥＰＩ法によって撮像された拡散強調画像は静磁場強度の不均一性の影響を強く受け、静磁場強度が十分に均一でないと再構成された画像に歪みが生じる。

そこで、静磁場強度の不均一性による問題を解決する方法として、従来、静磁場強度の分布を推定し、推定して得られた静磁場強度の分布をもとに再構成された画像の位置および画像値を補正するという方法がとられている（例えば非特許文献１参照）。

静磁場強度の均一性を乱し、再構成された画像に歪みをもたらす要因としては、傾斜磁場の誤差、静磁場用磁石単体による静磁場誤差および被検体が静磁場内に入ることによる静磁場の誤差の３つの要因がある。このため、静磁場の誤差は、傾斜磁場の誤差、静磁場用磁石単体による静磁場誤差および被検体が静磁場内に入ることによる静磁場の誤差の総和で定義することができる。

さらに、静磁場の誤差から画像データの歪み量が画像空間の位置ごとに求められ、求められた画像データの歪み量に基づいて画像空間の位置ごとに画像値が補正される。すなわち、まず画像空間において、画像値が得られている３次元状の位置が、歪みがなかった場合の位置にシフトされる。しかし、シフトされた位置は、必ずしも本来あるべき格子点上とならないため、シフト後の各位置における複数の画像値から補間によって格子点上の画像値が求められる。
J.MICHIELS ET AL. "Geometric distortion in MRI for stereotactic neurosurgery" Magnetic Resonance Imaging, Volume 12,Nimber 5,1994,p.749-765

従来行われる静磁場強度の不均一性による画像歪みの補正においては、静磁場の各誤差要因に応じて静磁場の歪み量の分布を推定することが精度上重要である。特に静磁場内に入ることによる静磁場の誤差については、被検体が動くことによって変動するので、静磁場の誤差分布の推定値が誤差を持つ可能性がある。

そして、静磁場の誤差分布が精度よく推定されないまま、静磁場の誤差分布に基づいて画像の補正が行われると、適切に画像歪みが改善されない恐れがある。このため、静磁場の誤差分布が精度よく推定されているか否か、あるいは静磁場の不均一性による画像歪みの補正が適切に行われたか否かを評価する技術の開発が望まれる。

本発明はかかる従来の事情に対処するためになされたものであり、撮影領域に形成される静磁場の誤差分布が適切に推定されて静磁場の不均一性による画像歪みの補正が適切に行われたか否かを評価することが可能な磁気共鳴イメージング装置および画像補正評価方法を提供することを目的とする。

本発明に係る磁気共鳴イメージング装置は、上述の目的を達成するために、請求項１に記載したように、異なる複数の撮影条件を設定する撮影条件設定手段と、前記複数の撮影条件に従って撮影領域から磁気共鳴信号を受信する受信手段と、前記磁気共鳴信号に基づいて前記複数の撮影条件に対応する複数の画像データを再構成する画像再構成手段と、前記撮影領域の磁場分布に基づいて前記複数の画像データの歪みを補正する画像歪み補正手段と、補正後の複数の画像データに基づいて少なくとも１つの前記画像データの補正が適切に行われたか否かを評価する画像補正評価手段とを備えることを特徴とするものである。

また、本発明に係る画像補正評価方法は、上述の目的を達成するために、請求項１４に記載したように、異なる複数の撮影条件を設定するステップと、前記複数の撮影条件に従って撮影領域から磁気共鳴信号を受信するステップと、前記磁気共鳴信号に基づいて前記複数の撮影条件に対応する複数の画像データを再構成するステップと、前記撮影領域の磁場分布に基づいて前記複数の画像データの歪みを補正するステップと、補正後の複数の画像データに基づいて少なくとも１つの前記画像データの補正が適切に行われたか否かを評価するステップとを有することを特徴とするものである。

本発明に係る磁気共鳴イメージング装置および画像補正評価方法においては、撮影領域に形成される静磁場の誤差分布が適切に推定されて静磁場の不均一性による画像歪みの補正が適切に行われたか否かを評価することができる。

本発明に係る磁気共鳴イメージング装置および画像補正評価方法の実施の形態について添付図面を参照して説明する。

図１は本発明に係る磁気共鳴イメージング装置の実施の形態を示す構成図である。

磁気共鳴イメージング装置２０は、静磁場を形成する筒状の静磁場用磁石２１、この静磁場用磁石２１の内部に設けられたシムコイル２２、傾斜磁場コイルユニット２３およびＲＦコイル２４を図示しないガントリに内蔵した構成である。

また、磁気共鳴イメージング装置２０には、制御系２５が備えられる。制御系２５は、静磁場電源２６、傾斜磁場電源２７、シムコイル電源２８、送信器２９、受信器３０、シーケンスコントローラ３１およびコンピュータ３２を具備している。制御系２５の傾斜磁場電源２７は、Ｘ軸傾斜磁場電源２７ｘ、Ｙ軸傾斜磁場電源２７ｙおよびＺ軸傾斜磁場電源２７ｚで構成される。また、コンピュータ３２には、入力装置３３、表示装置３４、演算装置３５および記憶装置３６が備えられる。

静磁場用磁石２１は静磁場電源２６と接続され、静磁場電源２６から供給された電流により撮像領域に静磁場を形成させる機能を有する。尚、静磁場用磁石２１は超伝導コイルで構成される場合が多く、励磁の際に静磁場電源２６と接続されて電流が供給されるが、一旦励磁された後は非接続状態とされるのが一般的である。また、静磁場用磁石２１を永久磁石で構成し、静磁場電源２６が設けられない場合もある。

また、静磁場用磁石２１の内側には、同軸上に筒状のシムコイル２２が設けられる。シムコイル２２はシムコイル電源２８と接続され、シムコイル電源２８からシムコイル２２に電流が供給されて静磁場が均一化されるように構成される。

傾斜磁場コイルユニット２３は、Ｘ軸傾斜磁場コイル２３ｘ、Ｙ軸傾斜磁場コイル２３ｙおよびＺ軸傾斜磁場コイル２３ｚで構成され、静磁場用磁石２１の内部において筒状に形成される。傾斜磁場コイルユニット２３の内側には寝台３７が設けられて撮像領域とされ、寝台３７には被検体Ｐがセットされる。ＲＦコイル２４はガントリに内蔵されず、寝台３７や被検体Ｐ近傍に設けられる場合もある。

また、傾斜磁場コイルユニット２３は、傾斜磁場電源２７と接続される。傾斜磁場コイルユニット２３のＸ軸傾斜磁場コイル２３ｘ、Ｙ軸傾斜磁場コイル２３ｙおよびＺ軸傾斜磁場コイル２３ｚはそれぞれ、傾斜磁場電源２７のＸ軸傾斜磁場電源２７ｘ、Ｙ軸傾斜磁場電源２７ｙおよびＺ軸傾斜磁場電源２７ｚと接続される。

そして、Ｘ軸傾斜磁場電源２７ｘ、Ｙ軸傾斜磁場電源２７ｙおよびＺ軸傾斜磁場電源２７ｚからそれぞれＸ軸傾斜磁場コイル２３ｘ、Ｙ軸傾斜磁場コイル２３ｙおよびＺ軸傾斜磁場コイル２３ｚに供給された電流により、撮像領域にそれぞれＸ軸方向の傾斜磁場Ｇｘ、Ｙ軸方向の傾斜磁場Ｇｙ、Ｚ軸方向の傾斜磁場Ｇｚを形成することができるように構成される。

ＲＦコイル２４は、送信器２９および受信器３０と接続される。ＲＦコイル２４は、送信器２９からＲＦ信号を受けて被検体Ｐに送信する機能と、被検体Ｐ内部の原子核スピンのＲＦ信号による励起に伴って発生したＭＲ信号を受信して受信器３０に与える機能を有する。

一方、制御系２５のシーケンスコントローラ３１は、傾斜磁場電源２７、送信器２９および受信器３０と接続される。シーケンスコントローラ３１は傾斜磁場電源２７、送信器２９および受信器３０を駆動させるために必要な制御情報、例えば傾斜磁場電源２７に印加すべきパルス電流の強度や印加時間、印加タイミング等の動作制御情報を記述したパルスシーケンスを記憶する機能と、記憶した所定のパルスシーケンスで規定される撮影条件に従って傾斜磁場電源２７、送信器２９および受信器３０を駆動させることによりＸ軸傾斜磁場Ｇｘ、Ｙ軸傾斜磁場Ｇｙ，Ｚ軸傾斜磁場ＧｚおよびＲＦ信号を発生させる機能を有する。

また、シーケンスコントローラ３１は、受信器３０におけるＭＲ信号の検波およびＡ／Ｄ変換により得られた複素データである生データ（raw data）を受けてコンピュータ３２に与えるように構成される。

このため、送信器２９には、シーケンスコントローラ３１から受けた制御情報に基づいてＲＦ信号をＲＦコイル２４に与える機能が備えられる一方、受信器３０には、ＲＦコイル２４から受けたＭＲ信号を検波して所要の信号処理を実行するとともにＡ／Ｄ変換することにより、デジタル化された複素データである生データを生成する機能と生成した生データをシーケンスコントローラ３１に与える機能とが備えられる。

すなわち、静磁場用磁石２１、シムコイル２２、傾斜磁場コイルユニット２３、ＲＦコイル２４並びに制御系２５の各構成要素により、磁気共鳴イメージング装置２０は、パルスシーケンスとして設定された各撮影条件に従って、静磁場中の被検体Ｐに対して傾斜磁場の印加およびＲＦ信号の送信を行なう一方、被検体Ｐ内部におけるＲＦ信号による原子核の核磁気共鳴に伴って発生したＭＲ信号を受信し、デジタル化することにより生データを生成する。

また、コンピュータ３２の記憶装置３６に保存されたプログラムを演算装置３５で実行することにより、コンピュータ３２には各種機能が備えられる。ただし、プログラムによらず、特定の回路を設けてコンピュータ３２を構成してもよい。

図２は、図１に示す磁気共鳴イメージング装置２０におけるコンピュータ３２の機能ブロック図である。

コンピュータ３２は、画像補正評価プログラムによりシーケンスコントローラ制御部４０、生データデータベース４１、画像再構成部４２、画像データデータベース４３、静磁場分布データベース４４、画像歪み補正部４５、撮影条件設定部４６、画像補正評価手段の一例としての相互相関計算部４７、画像補正評価手段の一例としての信号強度比較部４８、磁場分布計算部４９およびフィードバック部５０として機能する。フィードバック部５０は、近似曲線次数変更部５０Ａ、再補正指示部５０Ｂおよびシミング実行部５０Ｃを備えている。

シーケンスコントローラ制御部４０は、入力装置３３またはその他の構成要素からの情報に基づいてシーケンスコントローラ３１に所要のパルスシーケンスを与えることにより駆動制御させる機能と、シーケンスコントローラ３１から生データを受けて生データデータベース４１に形成されたｋ空間（フーリエ空間）に配置する機能とを有する。

撮影条件設定部４６は、パルスシーケンスを生成することにより撮影条件を設定する機能と、生成したパルスシーケンスをシーケンスコントローラ制御部４０に与える機能とを有する。ここで、撮影条件設定部４６には、静磁場分布測定用のプリスキャンの実行のためのパルスシーケンスを生成する他、位相エンコード（ＰＥ）方向、リードアウト（ＲＯ）方向およびスライスエンコード（ＳＥ）方向の少なくとも１方向の傾斜磁場の極性を反転させて、複数通りの撮影条件によりスキャンが実行されるように複数通りのパルスシーケンスを生成する機能が備えられる。

つまり、撮影条件設定部４６は、画素値（または信号強度。以下同様。）の値の差が十分に得られるような複数の画像を撮影するための画像収集（イメージング）用の撮影条件を設定することができるように構成される。ここで、画素値（または信号強度）の差とは、複数の画像間で対応する位置における差である。十分な画素値の差を有する複数の画像を撮影するための撮影条件は、傾斜磁場の極性を反転させることにより容易に設定することができる。

この他、T1（縦）緩和時間やT2（横）緩和時間等の緩和時間Tが異なる複数の撮影条件や、FE(field echo)法やSE(spin echo)法等のイメージング法が異なる複数の撮影条件を設定することによっても、対応する位置において十分な画素値や信号値の差を得ることができる。

従って、十分な画素値の差を有する複数の画像を撮影するための撮影条件が設定されれば、傾斜磁場の極性を反転させる方法以外の方法により撮影条件を設定することも可能であるが、以下では、傾斜磁場の極性を反転させて複数の撮影条件を設定する場合について説明する。

尚、撮影条件として設定される画像収集用のパルスシーケンスとしては、例えばＥＰＩシーケンスが挙げられる。ＥＰＩシーケンスによる撮影では、静磁場強度の不均一性の影響を特に強く受けることが知られている。

このため、生データデータベース４１には、撮影条件設定部４６により設定された複数通りの撮影条件のスキャンの実行によって受信器３０において生成された各生データが保存され、生データデータベース４１に形成されたｋ空間に各生データがそれぞれ撮影条件ごとに配置される。また、静磁場分布測定用のプリスキャンが実行された場合には、静磁場分布測定用の生データが保存される。

画像再構成部４２は、生データデータベース４１から生データを取り込んでフーリエ変換処理や最大値投影処理等の所定の画像再構成処理を施すことにより、被検体Ｐの画像データを再構成する機能と、得られた画像データを画像データデータベース４３に書き込む機能とを有する。ただし、診断時以外には、被検体Ｐ以外の例えばファントムの画像データを再構成するようにしてもよい。

このため、画像データデータベース４３には、撮影条件ごとの被検体Ｐやファントムの画像データが保存される。

静磁場分布データベース４４には、予め静磁場分布が計測あるいは求められて保存される。静磁場分布データベース４４には、磁場分布計算部４９によって求められた静磁場分布を保存することもできる。

画像歪み補正部４５は、画像データデータベース４３に保存された撮影条件ごとの被検体Ｐやファントムの画像データの歪みを静磁場分布データベース４４から読み込んだ静磁場分布に基づいてそれぞれ補正する機能と、補正後の各画像データを相互相関計算部４７および信号強度比較部４８にそれぞれ与える機能とを有する。

相互相関計算部４７は、複数通りの撮影条件によって得られた補正後の画像データを画像歪み補正部４５から受けて、複数の画像データ間の相互相関係数を求める機能と、求められた相互相関係数に基づいて画像データの補正が適切に行われたか否かを評価する画像補正評価手段としての機能とを有する。また、必要に応じて相互相関計算部４７には、求められた複数の画像データ間の相互相関係数や画像データの補正が適切に行われたか否かの評価結果を表示装置３４に与えることにより表示させる機能が備えられる。

補正が適切であるか否かの評価方法としては、例えば、画像データ間の相互相関係数を予め設定された相互相関係数の閾値と比較し、相互相関係数が閾値以上または閾値を超える値である場合には補正が適切に行われたと評価する一方、相互相関係数が閾値未満または閾値以下である場合には補正が適切に行われていないと評価する方法が挙げられる。

そして、相互相関計算部４７は、補正が適切でないと判定した場合には、フィードバック部５０に補正を適正化するようにフィードバック指示を与えるように構成される。このとき、補正が適切でないと判定される領域が局所的である場合には、補正が適切でないと判定される領域も相互相関計算部４７からフィードバック部５０に通知されるように構成される。また、必要に応じて相互相関係数の閾値からの乖離量も相互相関計算部４７からフィードバック部５０に通知されるように構成される。

信号強度比較部４８は、複数通りの撮影条件によって得られた補正後の画像データを画像歪み補正部４５から受けて、複数の画像データ間でピクセルごとに信号強度を比較する機能と、信号強度の比較結果にに基づいて画像データの補正が適切に行われたか否かを評価する画像補正評価手段としての機能とを有する。また、必要に応じて信号強度比較部４８には、複数の画像データ間における信号強度の比較結果や画像データの補正が適切に行われたか否かの評価結果を表示装置３４に与えることにより表示させる機能が備えられる。

補正が適切であるか否かの評価方法としては、例えば、画像データ間の信号強度の差分値を予め設定された信号強度の差分値の閾値と比較し、信号強度の差分値が閾値以上または閾値を超える値である場合には補正が適切に行われていないと評価する一方、信号強度の差分値が閾値未満または閾値以下である場合には補正が適切に行われたと評価する方法が挙げられる。

そして、信号強度比較部４８は、補正が適切でないと判定した場合には、フィードバック部５０に補正を適正化するようにフィードバック指示を与えるように構成される。このとき、補正が適切でないと判定される領域が局所的である場合には、補正が適切でないと判定される領域も信号強度比較部４８からフィードバック部５０に通知されるように構成される。また、必要に応じて信号強度の差分値の閾値からの乖離量も信号強度比較部４８からフィードバック部５０に通知されるように構成される。

磁場分布計算部４９は、任意の方法により静磁場の分布を推定する機能と、推定して得られた静磁場分布データを静磁場分布データベース４４に書き込む機能とを有する。静磁場の分布の推定には、必要に応じて生データデータベース４１に保存された静磁場分布測定用の生データが参照される。さらに、傾斜磁場を測定するための傾斜磁場測定コイル５１を撮影領域近傍に設けて、傾斜磁場測定コイル５１による傾斜磁場の測定結果を読み込んで参照するように構成することもできる。

磁場分布計算部４９により推定される静磁場分布データは、望ましくは、位置ごとの静磁場の値に基づいて位置と静磁場の値を関連付けた任意の次数の近似曲線として求められる。近似曲線を求める際には、測定あるいは設計シミュレーションにより求められた静磁場の誤差データが参照される。このように、静磁場分布データを近似曲線で表せば、補正すべき点がいかなる位置であってもその点の誤差を求めることができる。

ただし、磁場分布計算部４９により推定される静磁場分布データを近似曲線とせずに、位置ごとの静磁場の値としてもよい。

フィードバック部５０は、相互相関計算部４７および信号強度比較部４８による補正の評価結果を任意の方法でフィードバックさせる機能を有する。すなわち、フィードバック部５０は、相互相関計算部４７または信号強度比較部４８からフィードバック指示を受けた場合に、近似曲線次数変更部５０Ａ、再補正指示部５０Ｂおよびシミング実行部５０Ｃの少なくとも１つを用いて、補正を適正化させるように構成される。

近似曲線次数変更部５０Ａは、相互相関計算部４７または信号強度比較部４８からフィードバック指示を受けた場合に、静磁場分布データとして静磁場分布データベース４４に保存された近似曲線の次数を増加させる指示を磁場分布計算部４９に与える機能と、磁場分布計算部４９により新たに得られた静磁場分布データを用いて画像データの補正を実行するように画像歪み補正部４５に指示を与える機能とを有する。つまり、静磁場分布データの近似曲線の次数を増加させると、静磁場分布の近似精度が向上し、その結果、静磁場分布データに基づく画像の歪み補正の精度が向上すると考えられる。

そこで、近似曲線次数変更部５０Ａは、磁場分布計算部４９に近似曲線の次数を増加させる指示を与えることにより、静磁場分布データの再計算を行わせ、静磁場分布データベース４４に保存される静磁場分布データの精度を向上させる。そして、近似曲線次数変更部５０Ａは、精度が向上された静磁場分布を用いて再度、画像データの歪み補正を実施するように画像歪み補正部４５に指示を与えるように構成される。

ただし、画像データの歪み補正が適切でないと判定される領域が局所的である場合には、処理量を低減させるために、画像データの歪み補正が適切でないと判定された領域に限定して静磁場分布データの精度を向上させるようにすることもできる。

再補正指示部５０Ｂは、相互相関計算部４７または信号強度比較部４８からフィードバック指示を受けた場合に、画像歪み補正部４５に補正に用いた静磁場分布データの値をある方向に変化させて再度画像データの補正を行うように指示する機能と、静磁場分布データを変化させる方向が、補正が改善される方向であるか否かを判定する機能を有する。

静磁場分布データを変化させる方向が、補正が改善される方向であるか否かの判定は、画像歪み補正部４５による再補正により、相互相関係数または信号強度の差分値のそれぞれの閾値に対する乖離量が低減されたか否かを判定することにより行うことができる。そして、乖離量が低減されていれば、再度静磁場分布データを同じ方向に変化させる一方、乖離量が増加していれば静磁場分布データを逆方向に変化させて再度画像データの補正を行うように画像歪み補正部４５に指示することができる。また、再度静磁場分布データを変化させる量を次第に小さくしていくことで、静磁場分布データを最適化することができる。

つまり、静磁場分布データをある方向に変化させて再度画像データの補正を実施し、補正の評価結果が改善されるようであれば、補正が適切であると判定されるまで繰り返し、静磁場分布データの変化および画像データの補正を行うという収束計算により、画像データの補正適正化を図ることができる。また、再補正の評価結果が改善されない場合には、静磁場分布データの変化方向を逆方向に変えて、再補正の評価が試みられる。そして、このような集束方向の判定および集束計算を繰返し行えば、相互相関係数または信号強度の差分値が閾値内となり、補正が適切に行われたと判定されるものと考えられる。

ただし、画像データの歪み補正が適切でないと判定される領域が局所的である場合には、処理量を低減させるために、画像データの歪み補正が適切でないと判定された領域に限定して静磁場分布データを変化させるようにすることもできる。

シミング実行部５０Ｃは、相互相関計算部４７または信号強度比較部４８からフィードバック指示を受けた場合に、シムコイル電源２８に制御信号を与えることにより、シミングを実施させる機能を有する。すなわち、シムコイル電源２８からシムコイル２２に出力される電流または電圧を制御することにより、シミング条件を変化させ、静磁場分布の均一性を向上させることができる。

静磁場分布の均一性が向上すると、画像データの補正精度も向上するため、補正の適切化を図ることができる。尚、シミング条件の調整の際には、予めシミング条件と静磁場分布の均一性を関連付けたシミング条件テーブルを参照することができる。このシミング条件テーブルは、必要に応じてシミング実行部５０Ｃに保存されるか参照できるように構成される。

次に磁気共鳴イメージング装置２０の動作について説明する。

図３は、図１に示す磁気共鳴イメージング装置２０により、静磁場の不均一性による画像歪みの補正が適切に行われたか否かを評価する際の流れの一例を示すフローチャートであり、図中Ｓに数字を付した符号はフローチャートの各ステップを示す。

まずステップＳ１において、予め任意の方法により静磁場の誤差分布が推定されて静磁場分布データベース４４に保存される。静磁場の誤差分布は計測してもよい。

静磁場の誤差分布の推定手法の一例としては、静磁場の誤差の発生要因を傾斜磁場の誤差、静磁場用磁石２１単体による静磁場誤差および被検体Ｐが静磁場内に入ることによる静磁場の誤差の３つの要因に分類し、各要因による静磁場の誤差をそれぞれ推定した後、加算する方法がある。

第１の要因である傾斜磁場の誤差は、理想的には線形であるべき傾斜磁場の線形からの誤差により生じるものである。傾斜磁場の誤差に基づく静磁場の歪み量は、傾斜磁場の強度とは独立して決定される。

撮影領域におけるガントリ座標系の任意の点の座標をｘ、ｙ、ｚとすると、傾斜磁場の誤差は、ｘ，ｙ，ｚ方向の各チャンネルにおいて勾配を有する３次元状の分布となる。そして、理想的な線形の傾斜磁場の強度をＧｘ，Ｇｙ，Ｇｚとして、ｘ，ｙ，ｚ方向の各チャンネルにおける傾斜磁場の誤差の割合ＧｘＥｒｒｏｒＲａｔｉｏ（ｘ，ｙ，ｚ）、ＧｙＥｒｒｏｒＲａｔｉｏ（ｘ，ｙ，ｚ）、ＧｚＥｒｒｏｒＲａｔｉｏ（ｘ，ｙ，ｚ）を、式（１−１）、式（１−２）、式（１−３）のように表すと、ΔＧｘ（ｘ，ｙ，ｚ）、ΔＧｙ（ｘ，ｙ，ｚ）、ΔＧｚ（ｘ，ｙ，ｚ）は、ｘ，ｙ，ｚ方向の各チャンネルにおける傾斜磁場の誤差量を意味することになる。

［数１］
ＧｘＥｒｒｏｒＲａｔｉｏ（ｘ，ｙ，ｚ）＝ΔＧｘ（ｘ，ｙ，ｚ）／Ｇｘ…（１−１）
ＧｙＥｒｒｏｒＲａｔｉｏ（ｘ，ｙ，ｚ）＝ΔＧｙ（ｘ，ｙ，ｚ）／Ｇｙ…（１−２）
ＧｚＥｒｒｏｒＲａｔｉｏ（ｘ，ｙ，ｚ）＝ΔＧｚ（ｘ，ｙ，ｚ）／Ｇｚ…（１−３）

さらに、傾斜磁場の誤差量ΔＧｘ（ｘ，ｙ，ｚ）、ΔＧｙ（ｘ，ｙ，ｚ）、ΔＧｚ（ｘ，ｙ，ｚ）と、ガントリ座標系の座標ｘ、ｙ、ｚとからｘ，ｙ，ｚ方向の各チャンネルにおける傾斜磁場の誤差に起因する静磁場の誤差分布ΔＢｇｒａｄｘ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄｙ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄｚ（ｘ，ｙ，ｚ）を式（２−１）、式（２−２）、式（２−３）のように推定することができる。

［数２］
ΔＢｇｒａｄｘ（ｘ，ｙ，ｚ）＝ΔＧｘ（ｘ，ｙ，ｚ）×ｘ …（２−１）
ΔＢｇｒａｄｙ（ｘ，ｙ，ｚ）＝ΔＧｙ（ｘ，ｙ，ｚ）×ｙ …（２−２）
ΔＢｇｒａｄｚ（ｘ，ｙ，ｚ）＝ΔＧｚ（ｘ，ｙ，ｚ）×ｚ …（２−３）

また、第２の要因である静磁場用磁石２１単体による静磁場誤差は、理想的な静磁場の強度からの実際の静磁場の強度の誤差により生じるものである。静磁場用磁石２１単体による静磁場誤差は、撮像視野（ＦＯＶ：ｆｉｅｌｄ ｏｆ ｖｉｅｗ）が大きくなると急激に悪化するという性質がある。また、静磁場用磁石２１単体による静磁場誤差は、傾斜磁場の強度に依存する一方、静磁場中における被検体Ｐの形状や有無からは独立している。

また、第３の要因である被検体Ｐが静磁場内に入ることによる静磁場の誤差は、被検体Ｐが静磁場内に入ることで磁化率に違いが生じた結果、生じる静磁場の誤差分である。被検体Ｐが静磁場内に入ることによる静磁場の誤差は、傾斜磁場の強度に依存し、かつ静磁場中における被検体Ｐの形状や有無にも依存する。

そして、静磁場の誤差分布は、第１、第２、第３の各要因による静磁場の誤差成分の和として定義することができる。すなわち、ガントリ座標系の座標ｘ、ｙ、ｚにおける静磁場の誤差をΔＢｏｘ（ｘ，ｙ，ｚ）、ΔＢｏｙ（ｘ，ｙ，ｚ）、ΔＢｏｚ（ｘ，ｙ，ｚ）、傾斜磁場に起因する誤差成分をΔＢｇｒａｄｘ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄｙ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄｚ（ｘ，ｙ，ｚ）、静磁場用磁石２１単体による静磁場誤差をΔＢｍａｇ（ｘ，ｙ，ｚ）、被検体Ｐが静磁場内に入ることによる静磁場の誤差をΔＢｏｂｊ（ｘ，ｙ，ｚ）とすると、式（３−１）、式（３−２）、式（３−３）のように表すことができる。

［数３］
ΔＢｏｘ（ｘ，ｙ，ｚ）＝ΔＢｇｒａｄｘ（ｘ，ｙ，ｚ）＋ΔＢｍａｇ（ｘ，ｙ，ｚ）＋ΔＢｏｂｊ（ｘ，ｙ，ｚ） …（３−１）
ΔＢｏｙ（ｘ，ｙ，ｚ）＝ΔＢｇｒａｄｙ（ｘ，ｙ，ｚ）＋ΔＢｍａｇ（ｘ，ｙ，ｚ）＋ΔＢｏｂｊ（ｘ，ｙ，ｚ） …（３−２）
ΔＢｏｚ（ｘ，ｙ，ｚ）＝ΔＢｇｒａｄｚ（ｘ，ｙ，ｚ）＋ΔＢｍａｇ（ｘ，ｙ，ｚ）＋ΔＢｏｂｊ（ｘ，ｙ，ｚ） …（３−３）

従って、式（３−１）、式（３−２）、式（３−３）からガントリ座標系のｘ、ｙ、ｚ方向の静磁場の歪み量ΔＸ、ΔＹ，ΔＺは、式（４−１）、式（４−２）、式（４−３）のように求めることができる。

［数４］
ΔＸ＝ΔＢｏｘ（ｘ，ｙ，ｚ）／Ｇｘ
＝ＧｘＥｒｒｏｒＲａｔｉｏ（ｘ，ｙ，ｚ）×ｘ＋｛ΔＢｍａｇ（ｘ，ｙ，ｚ）＋ΔＢｏｂｊ（ｘ，ｙ，ｚ）｝／Ｇｘ …（４−１）
ΔＹ＝ΔＢｏｙ（ｘ，ｙ，ｚ）／Ｇｙ
＝ＧｙＥｒｒｏｒＲａｔｉｏ（ｘ，ｙ，ｚ）×ｙ＋｛ΔＢｍａｇ（ｘ，ｙ，ｚ）＋ΔＢｏｂｊ（ｘ，ｙ，ｚ）｝／Ｇｙ …（４−２）
ΔＺ＝ΔＢｏｚ（ｘ，ｙ，ｚ）／Ｇｚ
＝ＧｚＥｒｒｏｒＲａｔｉｏ（ｘ，ｙ，ｚ）×ｚ＋｛ΔＢｍａｇ（ｘ，ｙ，ｚ）＋ΔＢｏｂｊ（ｘ，ｙ，ｚ）｝／Ｇｚ …（４−３）

式（４−１）、式（４−２）、式（４−３）において、｛ΔＢｍａｇ（ｘ，ｙ，ｚ）＋ΔＢｏｂｊ（ｘ，ｙ，ｚ）｝の値は、イメージング用のスキャンに先立って静磁場分布測定用のプリスキャンを実行することにより求めることができる。イメージング用のスキャンがＥＰＩシーケンスによるスキャンである場合には、静磁場分布測定用のプリスキャンをＥＰＩシーケンスの受信ゲイン調整や傾斜磁場の遅れ時間調整のためのプリスキャンとともに実行することができる。

すなわち、静磁場分布測定用のプリスキャンの実行により設定されたスライス面内の静磁場分布（中心周波数からのずれを位相で表した分布図）を測定するためのデータを収集することができる。この場合、静磁場分布測定用のプリスキャンの実行の際における面内分解能はイメージング用のスキャンの実行の際における面内分解能よりも粗くすることができる。

静磁場分布測定用のプリスキャンの一例としては、ＦｉｅｌｄＥｃｈｏ法によりエコー時間（ＴＥ：ｅｃｈｏ ｔｉｍｅ）を変えて２つの画像データを収集するものが挙げられる。そして、式（５）に示すように、静磁場分布測定用のプリスキャンにより収集された各画像データのピクセルごとの位相差Δθ（ｘ，ｙ，ｚ）から｛ΔＢｍａｇ（ｘ，ｙ，ｚ）＋ΔＢｏｂｊ（ｘ，ｙ，ｚ）｝の値を求めることができる。

［数５］
Δθ（ｘ，ｙ，ｚ）＝（ＴＥ１−ＴＥ２）×｛ΔＢｍａｇ（ｘ，ｙ，ｚ）＋ΔＢｏｂｊ（ｘ，ｙ，ｚ）｝ …（５）
但し、式（５）において、ＴＥ１、ＴＥ２は、それぞれ互いに異なるＴＥである。

つまり、ＴＥを変えて収集された画像データを差分することで、静磁場の誤差ΔＢｏｘ（ｘ，ｙ，ｚ）、ΔＢｏｙ（ｘ，ｙ，ｚ）、ΔＢｏｚ（ｘ，ｙ，ｚ）のうちΔＢｇｒａｄｘ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄｙ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄｚ（ｘ，ｙ，ｚ）の項がキャンセルされる。

従って、傾斜磁場の誤差に起因する静磁場の誤差分布ΔＢｇｒａｄ．ｘ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄ．ｙ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄ．ｚ（ｘ，ｙ，ｚ）を求めることができれば、式（３−１），式（３−２）、式（３−３）および式（４−１），式（４−２）、式（４−３）から静磁場の誤差ΔＢｏｘ（ｘ，ｙ，ｚ）、ΔＢｏｙ（ｘ，ｙ，ｚ）、ΔＢｏｚ（ｘ，ｙ，ｚ）および静磁場の歪み量ΔＸ、ΔＹ，ΔＺを求めることができる。

傾斜磁場の誤差に起因する静磁場の誤差分布ΔＢｇｒａｄ．ｘ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄ．ｙ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄ．ｚ（ｘ，ｙ，ｚ）は、例えば傾斜磁場測定コイル５１により測定することができる。

そこで、傾斜磁場の誤差に起因する静磁場の誤差分布ΔＢｇｒａｄ．ｘ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄ．ｙ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄ．ｚ（ｘ，ｙ，ｚ）が傾斜磁場測定コイル５１により測定されて磁場分布計算部４９に与えられる。

一方、撮影条件設定部４６により静磁場分布測定用のプリスキャンのためのパルスシーケンスが生成されてシーケンスコントローラ制御部４０に与えられる。このため、シーケンスコントローラ制御部４０は、入力装置３３またはその他の構成要素からの情報に基づいてシーケンスコントローラ３１に静磁場分布測定用のパルスシーケンスを与えることにより静磁場分布測定用のプリスキャンを実行させる。

すなわち、予め寝台３７に被検体Ｐあるいはファントムがセットされ、静磁場電源２６により励磁された静磁場用磁石２１（超伝導磁石）の撮像領域に静磁場が形成される。また、シムコイル電源２８からシムコイル２２に電流が供給されて撮像領域に形成された静磁場が均一化される。

さらに、シーケンスコントローラ３１は、シーケンスコントローラ制御部４０から受けた静磁場分布測定用のパルスシーケンスに従って傾斜磁場電源２７、送信器２９および受信器３０を駆動させることにより被検体Ｐあるいはファントムがセットされた撮像領域にＸ軸傾斜磁場Ｇｘ、Ｙ軸傾斜磁場Ｇｙ，Ｚ軸傾斜磁場Ｇｚを形成させるとともに、ＲＦ信号を発生させる。

このため、送信器２９からパルスシーケンスに応じてＲＦコイル２４にＲＦ信号が与えられ、ＲＦコイル２４から被検体ＰにＲＦ信号が送信される。ＲＦ信号が被検体Ｐに送信されると、被検体Ｐの内部において原子核の核磁気共鳴により生じたＭＲ信号が、ＲＦコイル２４により受信されて受信器３０に与えられる。受信器３０は、ＲＦコイル２４からＭＲ信号を受けて、前置増幅、中間周波変換、位相検波、低周波増幅、フィルタリング等の各種信号処理を実行する。さらに、受信器３０は、ＭＲ信号をＡ／Ｄ変換することにより、デジタルデータのＭＲ信号である生データを生成する。そして、受信器３０は、生成した生データをシーケンスコントローラ３１に与える。

シーケンスコントローラ３１は、受信器３０から受けた生データをシーケンスコントローラ制御部４０に与え、シーケンスコントローラ制御部４０は生データデータベース４１に生データを入力する。このようにして、例えばＴＥを変えたＦｉｅｌｄＥｃｈｏ法による２回に亘る静磁場分布測定用のプリスキャンが実行され、生データデータベース４１には、静磁場分布測定用のプリスキャンにより収集された生データが蓄積される。

そして、磁場分布計算部４９は、静磁場分布測定用のプリスキャンの実行により収集された生データを生データデータベース４１から読み込む一方、傾斜磁場測定コイル５１から傾斜磁場の誤差に起因する静磁場の誤差分布ΔＢｇｒａｄ．ｘ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄ．ｙ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄ．ｚ（ｘ，ｙ，ｚ）を受けて、例えば式（３−１），式（３−２）、式（３−３）、式（４−１），式（４−２）、式（４−３）および式（５）により静磁場の誤差ΔＢｏｘ（ｘ，ｙ，ｚ）、ΔＢｏｙ（ｘ，ｙ，ｚ）、ΔＢｏｚ（ｘ，ｙ，ｚ）および静磁場の歪み量ΔＸ、ΔＹ，ΔＺを求める。

さらに、磁場分布計算部４９は、求められた静磁場の歪み量ΔＸ、ΔＹ，ΔＺを静磁場分布データとして静磁場分布データベース４４に書き込む。

尚、傾斜磁場の誤差に起因する静磁場の誤差分布ΔＢｇｒａｄ．ｘ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄ．ｙ（ｘ，ｙ，ｚ）、ΔＢｇｒａｄ．ｚ（ｘ，ｙ，ｚ）は、傾斜磁場測定コイル５１による測定の他、例えば系を模擬した磁場モデルを用いてシミュレーションによる計算で求めることもできる。

次に、ステップＳ２において、撮影条件設定部４６により被検体Ｐの画像収集用のパルスシーケンスが生成されることにより、撮影条件が設定される。このとき、生成される画像収集用のパルスシーケンスは、ＰＥ方向、ＲＯ方向およびＳＥ方向の少なくとも１方向の傾斜磁場の極性を反転させて、例えば２通りの撮影条件によりスキャンが実行されるような２通りのパルスシーケンスとされる。

生成される画像収集用のパルスシーケンスの例としては、前述のようにＥＰＩシーケンスが挙げられる。

図４は、図１に示す磁気共鳴イメージング装置２０により静磁場の不均一性による画像歪みの補正が適切に行われたか否かを評価するために生成されるパルスシーケンスの一例を示す図である。

図４において、ＳＥ１、ＲＯ１、ＰＥ１は、それぞれＳＥ方向、ＲＯ方向およびＰＥ方向の傾斜磁場を規定する第１のパルスシーケンスであり、ＳＥ２、ＲＯ２、ＰＥ２は、それぞれＳＥ方向、ＲＯ方向およびＰＥ方向の傾斜磁場を規定する第２のパルスシーケンスである。図４に示すように、例えばＰＥ方向およびＲＯ方向の２方向の傾斜磁場の極性を互いに反転させて、２通りの撮影条件によりスキャンが実行されるような第１、第２の２通りのパルスシーケンスが撮影条件設定部４６により生成される。

そして、撮影条件設定部４６により生成された２通りのパルスシーケンスは、シーケンスコントローラ制御部４０に与える。

次に、ステップＳ３において、撮影条件設定部４６により設定された２つの撮影条件に従ってスキャンが実行され、画像データが再構成される。すなわち、シーケンスコントローラ制御部４０は、入力装置３３またはその他の構成要素からの情報に基づいてシーケンスコントローラ３１に画像収集用の２つのパルスシーケンスを与えることにより画像収集用のスキャンを実行させる。

この結果、傾斜磁場の極性を互いに反転させて得られた各生データが、生データデータベース４１に形成されたｋ空間に撮影条件ごとにそれぞれ配置される。

さらに、画像再構成部４２は、生データデータベース４１から各生データを取り込んでフーリエ変換処理等の所定の画像再構成処理をそれぞれ施すことにより、２つの画像データを再構成する。このようにして得られた各画像データは、傾斜磁場を反転させて収集された生データから再構成されたため、互いに歪み方が異なる画像データとなる。そして、各画像データは画像データデータベース４３に書き込まれて保存される。

次に、ステップＳ４において、画像データデータベース４３に保存された歪み方の異なる各画像データが画像歪み補正部４５に読み込まれて、それぞれ静磁場分布データベース４４に保存された静磁場分布に基づく各画像データの歪みの補正が行われる。

図５は、図１に示す磁気共鳴イメージング装置２０により静磁場分布に基づいて画像データを補正する際の方法の一例を示す説明図である。

まず、静磁場分布から各画像データの歪み量が画像空間の位置ごとに求められる。例えばｘ方向のみに無視できない歪が存在する場合には、ｘ方向の歪み量Δｘが求められる。そして、歪み量Δｘに基づいて、例えば歪み後の２点の座標（Ｘ_１，Ｙ_１，Ｚ_１）、（Ｘ_２，Ｙ_２，Ｚ_２）におけるそれぞれの画像値Ｐ_１（Ｘ_１，Ｙ_１，Ｚ_１）、Ｐ_２（Ｘ_２，Ｙ_２，Ｚ_２）が歪み前における各座標（Ｘ_１’，Ｙ_１’，Ｚ_１’）、（Ｘ_２’，Ｙ_２’，Ｚ_２’）における画像値Ｐ_１（Ｘ_１’，Ｙ_１’，Ｚ_１’）、Ｐ_２（Ｘ_２’，Ｙ_２’，Ｚ_２’）とされる。つまり、ｘ方向に歪み量Δｘだけ画像値Ｐ_１（Ｘ_１，Ｙ_１，Ｚ_１）、Ｐ_２（Ｘ_２，Ｙ_２，Ｚ_２）の位置（Ｘ_１，Ｙ_１，Ｚ_１）、（Ｘ_２，Ｙ_２，Ｚ_２）がシフトされる。

さらに、本来あるべき格子点の座標（ｘ，ｙ，ｚ）における画像値Ｐ（ｘ，ｙ，ｚ）が、格子点の座標（ｘ，ｙ，ｚ）近傍の位置（Ｘ_１’，Ｙ_１’，Ｚ_１’）、（Ｘ_２’，Ｙ_２’，Ｚ_２’）における画像値Ｐ_１（Ｘ_１’，Ｙ_１’，Ｚ_１’）、Ｐ_２（Ｘ_２’，Ｙ_２’，Ｚ_２’）から補間によって画像空間の位置ごとに求められる。このとき通常、補間にはｃｕｂｉｃ ｓｐｌｉｎｅ以上の関数が用いられ、Ｌｉｎｅａｒではボケが大きく不十分であると考えられる。

そして、このようにして画像歪み補正部４５により得られた補正後の２つの各画像データは相互相関計算部４７および信号強度比較部４８にそれぞれ与えられる。

次に、ステップＳ５において、相互相関計算部４７は、傾斜磁場を反転させた２つの撮影条件によって得られた補正後の画像データ間の相互相関係数を求め、求められた相互相関係数に基づいて画像データの補正が適切に行われたか否かを評価する。

つまり、２つの画像データの収集条件は、傾斜磁場の極性の向きが互いに反転していることを除けば実質的に同一であるため、画像データの歪みの補正が適切に行われている程、より各画像データ上の物体の形状が一致し、２つの画像データ間における相互相関が強くなって相互相関係数が大きくなると考えられる。逆に、画像データの歪みの補正が適切に行われていない場合には、２つの画像データ間における相互相関が弱くなって相互相関係数が小さくなると考えられる。

そこで、相互相関計算部４７は、２つの画像データ間の相互相関係数を予め設定された閾値と比較し、相互相関係数が閾値以上または閾値を超える値である場合には補正が適切に行われたと評価する一方、相互相関係数が閾値未満または閾値以下である場合には補正が適切に行われていないと評価する。

また、必要に応じて相互相関計算部４７は、求められた２つの画像データ間の相互相関係数や画像データの補正が適切に行われたか否かの評価結果を表示装置３４に与えて表示させる。このため、ユーザは、画像データの補正が適切に行われたか否かを容易に把握することができる。

また、ステップＳ６において、信号強度比較部４８は、傾斜磁場を反転させた２つの撮影条件によって得られた補正後の画像データ間でピクセルごとに信号強度を比較し、信号強度の比較結果に基づいて画像データの補正が適切に行われたか否かを評価する。

つまり、前述のように２つの画像データの収集条件は、傾斜磁場の極性の向きが互いに反転していることを除けば実質的に同一であるため、画像データの歪みの補正が適切に行われている程、より各画像データ上の物体の形状が一致し、２つの画像データ間における信号強度の差分が小さくなると考えられる。逆に、画像データの歪みの補正が適切に行われていない場合には、２つの画像データ間における信号強度の差分が大きくなると考えられる。

そこで、信号強度比較部４８は、２つの画像データ間の信号強度の差分値を予め設定された閾値と比較し、信号強度の差分値が閾値以上または閾値を超える値である場合には補正が適切に行われていないと評価する一方、信号強度の差分値が閾値未満または閾値以下である場合には補正が適切に行われたと評価する。

また、必要に応じて信号強度比較部４８は、求められた２つの画像データ間の信号強度の差分値や画像データの補正が適切に行われたか否かの評価結果を表示装置３４に与えて表示させる。このため、ユーザは、画像データの補正が適切に行われたか否かを容易に把握することができる。

次に、ステップＳ７において、画像データの補正が適切に行われていないと相互相関計算部４７および信号強度比較部４８の少なくとも一方により判定されると、相互相関計算部４７および／または信号強度比較部４８から、フィードバック部５０に補正を適正化するようにフィードバック指示が与えられる。

このとき、補正が適切でないと判定される領域が局所的である場合には、補正が適切でないと判定される領域も相互相関計算部４７および／または信号強度比較部４８からフィードバック部５０に通知される。また、必要に応じて相互相関係数の閾値からの乖離量が相互相関計算部４７からフィードバック部５０に通知され、信号強度の差分値の閾値からの乖離量が信号強度比較部４８からフィードバック部５０に通知される。

そうすると、ステップＳ８において、予め指定された任意の方法で補正のフィードバック処理が行われる。ここでは、フィードバック部５０において、静磁場分布データベース４４に静磁場分布データとして保存された近似曲線の次数を増加させる方法、静磁場分布データの値をある方向に変化させて再度画像データの補正を行う方法およびシミング条件を変化させ、静磁場分布の均一性を向上させた後に再度データ収集を行う方法のいずれかがフィードバック処理として指定されている場合について説明する。

静磁場分布データを示す近似曲線の次数を増加させる方法が選択されている場合には、相互相関計算部４７または信号強度比較部４８からフィードバック指示が近似曲線次数変更部５０Ａに与えられる。

そうすると、ステップＳ８Ａにおいて、近似曲線次数変更部５０Ａは、静磁場分布データとして静磁場分布データベース４４に保存された近似曲線の次数を増加させる指示を磁場分布計算部４９に与える一方、磁場分布計算部４９により新たに得られた静磁場分布データを用いて画像データの補正を実行するように画像歪み補正部４５に指示を与える。このため、磁場分布計算部４９は、近似曲線の次数を増加させて再度、静磁場分布データを計算する。

そして、再びステップＳ４において、画像歪み補正部４５は、新たに得られた静磁場分布データを用いて画像データの補正を実行する。このとき、静磁場分布データの近似曲線の次数が増加されたため、静磁場分布の近似精度が向上している。従って、静磁場分布データに基づく画像の歪み補正の精度が向上する。

ただし、画像データの歪み補正が適切でないと判定される領域が局所的である場合には、近似曲線次数変更部５０Ａから画像データの歪み補正が適切でないと判定された領域が磁場分布計算部４９および／または画像歪み補正部４５に通知されるようにし、磁場分布計算部４９および／または画像歪み補正部４５が当該領域に限定して静磁場分布データの計算および／または画像データの補正処理を行うようにしてもよい。このような領域限定により処理量を低減させることができる。

次に、同様にステップＳ５およびステップＳ６において、相互相関計算部４７および信号強度比較部４８により画像データの補正が適切であるか否かが判定される。そして、ステップＳ７において、画像データの補正が適切であると判定されるまで、同様なフィードバック処理が繰返し実行される。

一方、フィードバック処理として、静磁場分布データの値をある方向に変化させて再度画像データの補正を行う方法が選択されている場合には、相互相関計算部４７または信号強度比較部４８からフィードバック指示が再補正指示部５０Ｂに与えられる。このとき、相互相関係数および／または信号強度の差分値のそれぞれの閾値に対する乖離量が相互相関計算部４７および／または信号強度比較部４８から再補正指示部５０Ｂに通知される。

そうすると、ステップＳ８Ｂにおいて、再補正指示部５０Ｂは、画像歪み補正部４５に補正に用いた静磁場分布データの値をある方向に変化させて再度画像データの補正を行うように指示する。

このため、再びステップＳ４において再補正指示部５０Ｂは、静磁場分布データの値をある方向に変化させて再度画像データの補正を行う。そして、ステップＳ５およびステップＳ６において、相互相関計算部４７および信号強度比較部４８により画像データの補正が適切であるか否かが判定される。

ここで、静磁場分布データを変化させる方向が適切でない場合には、相互相関係数または信号強度の差分値のそれぞれの閾値に対する乖離量が増加するため、ステップＳ７において、相互相関計算部４７および信号強度比較部４８により画像データの補正が適切でないと判定されることとなる。

このような場合には、ステップＳ８Ｂにおいて、再補正指示部５０Ｂが、相互相関係数または信号強度の差分値のそれぞれの閾値に対する乖離量が増加し、静磁場分布データを変化させる方向が、補正が改善される方向でないと判定する。そして、再補正指示部５０Ｂは、静磁場分布データを逆方向に変化させて再度画像データの補正を行うように画像歪み補正部４５に指示する。

逆に、静磁場分布データを変化させる方向が適切であり、相互相関係数または信号強度の差分値のそれぞれの閾値に対する乖離量が減少した場合であっても、未だ相互相関係数または信号強度の差分値がそれぞれの閾値の範囲内とならない場合がある。このような場合にも、ステップＳ７において、相互相関計算部４７および信号強度比較部４８により画像データの補正が適切でないと判定されることとなる。

このような場合には、ステップＳ８Ｂにおいて、再補正指示部５０Ｂが、相互相関係数または信号強度の差分値のそれぞれの閾値に対する乖離量が減少し、静磁場分布データを変化させる方向が、補正が改善される方向であると判定する。そして、再補正指示部５０Ｂは、再度静磁場分布データを同じ方向に変化させて画像データの補正を行うように画像歪み補正部４５に指示する。

そして、ステップＳ７において、画像データの補正が適切であると判定されるまで、同様なフィードバック処理が繰返し実行される。ここで、静磁場分布データの変化を繰り返した結果、静磁場分布データを変化させる方向が、補正が改善される方向から補正が改善されない方向に変わった場合には、静磁場分布データの最適値が変化前後の値の間にあると考えられる。従って、このような場合には、再補正指示部５０Ｂにより静磁場分布データの変化方向が逆方向とされるとともに変化量が小さく設定される。

また、画像データの歪み補正が適切でないと判定される領域が局所的である場合には、再補正指示部５０Ｂから画像データの歪み補正が適切でないと判定された領域が画像歪み補正部４５に通知されるようにし、画像歪み補正部４５が当該領域に限定して画像データの補正処理を行うようにしてもよい。このような領域限定により処理量を低減させることができる。

一方、フィードバック処理として、シミング条件を変化させ、静磁場分布の均一性を向上させた後に再度データ収集を行う方法が選択されている場合には、相互相関計算部４７または信号強度比較部４８からフィードバック指示がシミング実行部５０Ｃに与えられる。

そうすると、ステップＳ８Ｃにおいて、シミング実行部５０Ｃは、予めシミング条件と静磁場分布の均一性を関連付けたシミング条件テーブルを参照し、静磁場分布の均一性が向上するように、シミング条件を設定する。そして、シミング実行部５０Ｃは、設定したシミング条件に従ってシムコイル電源２８に制御信号を与えることにより、シムコイル電源２８からシムコイル２２に出力される電流または電圧を制御する。これにより、再設定されたシミング条件に従ったシミングが実施され、静磁場分布の均一性が向上する。

このため、再度スキャンを実行すれば、より均一な静磁場分布下において、画像データを収集することが可能となるのみならず、補正に用いる静磁場分布データもより均一となる。この結果、画像データの補正精度を向上させることができる。

そして、ステップＳ７において、画像データの補正が適切であると判定されるまで、ステップＳ１からステップＳ６までの処理やスキャン、さらには同様なフィードバック処理が繰返し実行される。

このようにして、ステップＳ７において、画像データの補正が適切に行われたと相互相関計算部４７および信号強度比較部４８により判定された場合には、画像データ補正のフィードバックは行われず、良好に歪みが補正された画像データが表示装置３４に表示される。このため、ユーザは、適切に歪みが補正された画像データを参照して診断を行うことができる。

尚、上述の例では、単一の方法によるフィードバック処理が指定された場合について説明したが、複数の方法によるフィードバック処理を指定された優先順位で併用してもよい。例えば、静磁場分布データにおける近似曲線の次数の増加および静磁場分布データのある方向への変化を一定数回繰り返しても依然画像データが適切に補正されない場合に、シミングを実施して再びデータ収集するようにしてもよい。

つまり以上のような磁気共鳴イメージング装置２０は、傾斜磁場を反転させて２つの画像データを作成し、歪み補正後の２つの画像データの相関係数や信号強度の比較結果に基づいて補正の適否を評価するものである。

また、補正が適切でないとの評価が得られた場合には、適切な補正が実施されるようにフィードバック処理を行うものである。

このため、磁気共鳴イメージング装置２０によれば、撮影領域に形成される静磁場の誤差分布が適切に推定されて静磁場の不均一性による画像歪みの補正が適切に行われたか否かを評価することができる。

特に、静磁場強度の不均一性の影響を特に強く受けるＥＰＩシーケンスによる撮影では、画像歪みの補正が適切に行われていない恐れがあるが、磁気共鳴イメージング装置２０によれば、画像歪みの補正が適切に行われたか否かを評価することができる。

また、画像歪みの補正が適切に行われない場合でも、フィードバック処理により最終的により適切に補正された画像データを得ることができる。

尚、磁気共鳴イメージング装置２０の一部の機能や処理を省略してもよく、処理順序を変更してもよい。また、撮影条件は２通りに限らず、例えば、基準となる撮影条件、ＳＥ方向のみの傾斜磁場を反転させた撮影条件、ＲＯ方向のみの傾斜磁場を反転させた撮影条件の３通りの撮影条件とする場合のように、３つ以上の撮影条件としてもよい。

また、上述の実施形態では、静磁場強度の不均一性に起因する画像歪みの補正について評価を行う場合について説明したが、静磁場に限らず傾斜磁場強度の不均一性に起因する画像歪みの補正についても同様な原理に基づいて画像歪みの補正結果を評価することができる。すなわち、異なる撮影条件により得られ、かつ傾斜磁場分布の不均一性に応じて歪みが補正された複数の画像間の相互相関や信号強度の相異に基づいて、傾斜磁場強度分布に応じた画像歪みの補正結果を評価することができる。

さらに、画像歪みの補正が適切でないと評価された場合には、フィードバック処理により、適切に補正された画像データを得ることができる。

本発明に係る磁気共鳴イメージング装置の実施の形態を示す構成図。 図１に示す磁気共鳴イメージング装置におけるコンピュータの機能ブロック図。 図１に示す磁気共鳴イメージング装置により、静磁場の不均一性による画像歪みの補正が適切に行われたか否かを評価する際の流れの一例を示すフローチャート。 図１に示す磁気共鳴イメージング装置により静磁場の不均一性による画像歪みの補正が適切に行われたか否かを評価するために生成されるパルスシーケンスの一例を示す図。 図１に示す磁気共鳴イメージング装置により静磁場分布に基づいて画像データを補正する際の方法の一例を示す説明図。

符号の説明

２０ 磁気共鳴イメージング装置
２１ 静磁場用磁石
２２ シムコイル
２３ 傾斜磁場コイルユニット
２４ ＲＦコイル
２５ 制御系
２６ 静磁場電源
２７ 傾斜磁場電源
２８ シムコイル電源
２９ 送信器
３０ 受信器
３１ シーケンスコントローラ
３２ コンピュータ
３３ 入力装置
３４ 表示装置
３５ 演算装置
３６ 記憶装置
３７ 寝台
４０ シーケンスコントローラ制御部
４１ 生データデータベース
４２ 画像再構成部
４３ 画像データデータベース
４４ 静磁場分布データベース
４５ 画像歪み補正部
４６ 撮影条件設定部
４７ 相互相関計算部
４８ 信号強度比較部
４９ 磁場分布計算部
５０ フィードバック部
５０Ａ 近似曲線次数変更部
５０Ｂ 再補正指示部
５０Ｃ シミング実行部
５１ 傾斜磁場測定コイル
Ｐ 被検体

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