(57)【特許請求の範囲】
【請求項１】
被検体が置かれる空間に静磁場を発生する静磁場発生手段と、前記被検体に高周波磁場パルスを照射する照射手段と、高周波磁場パルスの照射により前記被検体から生じる核磁気共鳴信号を受信する受信手段と、前記核磁気共鳴信号に位置情報を付加するための傾斜磁場を発生する傾斜磁場発生手段と、前記照射手段、前記傾斜磁場発生手段及び前記受信手段を所定のパルスシーケンスに従って動作させる制御手段と、前記核磁気共鳴信号から画像を再構成する画像再構成手段と、前記傾斜磁場の強度分布の歪みによる画像の歪みを補正する歪み補正手段とを備える磁気共鳴イメージング装置であって、
前記歪み補正手段による補正後の画像から、当該補正により発生する虚像を除去する虚像除去手段を備え、
前記虚像除去手段は、前記補正後の画像の、前記傾斜磁場の強度分布が極値を取る領域に対応する画素より外側の画素の画素値を所定の値に置き換えることにより、前記虚像を除去すること
を特徴とする磁気共鳴イメージング装置。
【請求項２】
請求項１記載の磁気共鳴イメージング装置であって、
前記虚像除去手段は、前記補正後の画像の、ヤコビアンが所定の閾値以下となる画素を、前記極値を取る領域に対応する画素とすること
を特徴とする磁気共鳴イメージング装置。
【請求項３】
請求項１記載の磁気共鳴イメージング装置であって、
前記虚像除去手段は、前記歪み補正による各画素の微小領域の大きさの変化が所定の閾値以下の画素を、前記極値を取る領域に対応する画素とすること
を特徴とする磁気共鳴イメージング装置。
【請求項４】
請求項１から３いずれか１項記載の磁気共鳴イメージング装置であって、
前記静磁場発生手段は、前記被検体が置かれる空間側に窪み部を備え、
前記傾斜磁場発生手段は、前記窪み部に少なくとも一部が収納されていること
を特徴とする磁気共鳴イメージング装置。
【請求項５】
被検体が置かれる空間に静磁場を発生する静磁場発生手段と、前記被検体に高周波磁場パルスを照射する照射手段と、高周波磁場パルスの照射により前記被検体から生じる核磁気共鳴信号を受信する受信手段と、前記核磁気共鳴信号に位置情報を付加するための傾斜磁場を発生する傾斜磁場発生手段と、前記照射手段、前記傾斜磁場発生手段及び前記受信手段を所定のパルスシーケンスに従って動作させる制御手段とを備える磁気共鳴イメージング装置における画像再構成方法であって、
前記核磁気共鳴信号から画像を再構成する画像再構成ステップと、
前記再構成された画像について、前記傾斜磁場の強度分布の歪みによる画像の歪みを補正する歪み補正ステップと、
前記歪み補正手段による補正後の画像から、当該補正により発生する虚像を除去する虚像除去ステップと、を備え、
前記虚像除去ステップは、前記補正後の画像の、前記傾斜磁場の強度分布が極値を取る領域に対応する画素より外側の画素の画素値を所定の値に置き換えることにより、前記虚像を除去すること
を特徴とする画像再構成方法。
【請求項６】
被検体が置かれる空間に静磁場を発生する静磁場発生手段と、前記被検体に高周波磁場パルスを照射する照射手段と、高周波磁場パルスの照射により前記被検体から生じる核磁気共鳴信号を受信する受信手段と、前記核磁気共鳴信号に位置情報を付加するための傾斜磁場を発生する傾斜磁場発生手段と、前記照射手段、前記傾斜磁場発生手段及び前記受信手段を所定のパルスシーケンスに従って動作させる制御手段とを有する磁気共鳴イメージング装置が備えるコンピュータを、
前記核磁気共鳴信号から画像を再構成する画像再構成手段と、
前記傾斜磁場の強度分布の歪みによる画像の歪みを補正する歪み補正手段と、
前記歪み補正手段による補正後の画像から、当該補正により発生する虚像を、前記補正後の画像の、前記傾斜磁場の強度分布が極値を取る領域に対応する画素より外側の画素の画素値を所定の値に置き換えることにより除去する虚像除去手段、として機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、磁気共鳴イメージング（ＭＲＩ）装置における取得画像の虚像除去技術に関する。特に、傾斜磁場強度の分布（傾斜磁場分布）の非線形性を要因とする虚像の除去技術に関する。
【背景技術】
【０００２】
ＭＲＩ装置は、被検体、特に人体の組織を構成する原子核スピンが発生する核磁気共鳴（ＮＭＲ）信号を計測し、その頭部、腹部、四肢等の形態や機能を２次元的に或いは３次元的に画像化する装置である。撮影においては、ＮＭＲ信号には、傾斜磁場によって異なる位相エンコードが付与されるとともに周波数エンコードされて、時系列データとして計測される。計測されたＮＭＲ信号は、２次元又は３次元フーリエ変換され画像に再構成される。
【０００３】
正しい再構成画像を得るためには、ＮＭＲ信号に位相エンコードを付与する傾斜磁場分布が線形でなければならない。しかし、現在のＭＲＩ装置では、磁石のボアを短くしたり傾斜磁場の速度を改善したりするため、傾斜磁場分布の線形性は犠牲にされている。このような傾斜磁場分布の非線形性は、再構成画像に歪みを生じさせる。従って、正しい画像を得ることができない。このため、傾斜磁場分布の歪みから、傾斜磁場分布の非線形性による再構成画像の歪み量を推定し、画像の歪みを補正する手法が提案されている（例えば、非特許文献１参照）。
【０００４】
【非特許文献１】Ａｎｄｒｅｗ  Ｊａｎｋｅ，  Ｈｕａｗｅｉ  Ｚｈａｏ，  Ｇａｒｙ  Ｊ．  Ｃｏｗｉｎ，  Ｇｒａｈａｍ  Ｊ．  Ｇａｌｌｏｗａｙ  ａｎｄ  Ｄａｖｉｄ  Ｍ．  Ｄｏｄｄｒｅｌｌ  “Ｕｓｅ  ｏｆ  Ｓｐｈｅｒｉｃａｌ  Ｈａｒｍｏｎｉｃ  Ｄｅｃｏｎｖｏｌｕｔｉｏｎ  Ｍｅｔｈｏｄｓ  ｔｏ  Ｃｏｍｐｅｎｓａｔｅ  ｆｏｒ  Ｎｏｎｌｉｎｅａｒ  Ｇｒａｄｉｅｎｔ  Ｅｆｆｅｃｔｓ  ｏｎ  ＭＲＩ  Ｉｍａｇｅｓ”  Ｍａｇｎｅｔｉｃ  Ｒｅｓｏｎａｎｃｅ  ｉｎ  Ｍｅｄｉｃｉｎｅ  ５２：１１５−１２２（２００４）
【発明の開示】
【発明が解決しようとする課題】
【０００５】
非特許文献１に開示されている歪み補正手法では、傾斜磁場強度より歪んだ傾斜磁場の位置座標を算出し、歪みを補正した画像を得ている。ここで、撮影視野内の傾斜磁場強度分布に極値点が存在すると、極値点の両側に同じ傾斜磁場強度となる点が存在し、算出される位置座標も同じとなる。従って、非特許文献１に開示の歪み補正によれば、原画像上の１つの点が極値点の両側の２つの点に移動する。すなわち、傾斜磁場分布の非線形性による歪みを、非特許文献１に開示の歪み補正手法を用いて補正すると、傾斜磁場強度が極値点を有する場合、その外側に本来なら存在しない画像である虚像が発生する。
【０００６】
本発明は、上記事情に鑑みてなされたもので、傾斜磁場分布が極値点を持つために歪み補正により発生する虚像を除去し、虚像のない歪み補正画像を得ることを目的とする。
【課題を解決するための手段】
【０００７】
本発明は、画像の歪み分布量から、傾斜磁場の極値点に対応する再構成画像上の点を特定し、当該点の、磁場中心より外側の画像の画素値を所定の値にすることにより、歪み補正により傾斜磁場の極値点内の信号が作り出す虚像を除去する。
【０００８】
具体的には、被検体が置かれる空間に静磁場を発生する静磁場発生手段と、前記被検体に高周波磁場パルスを照射する照射手段と、高周波磁場パルスの照射により前記被検体から生じる核磁気共鳴信号を受信する受信手段と、前記核磁気共鳴信号に位置情報を付加するための傾斜磁場を発生する傾斜磁場発生手段と、前記照射手段、前記傾斜磁場発生手段及び前記受信手段を所定のパルスシーケンスに従って動作させる制御手段と、前記核磁気共鳴信号から画像を再構成する信号処理手段とを備える磁気共鳴イメージング装置であって、前記傾斜磁場の強度分布の歪みによる画像の歪みを補正する歪み補正手段と、前記歪み補正手段による補正後の画像から、当該補正により発生する虚像を除去する虚像除去手段と、を備え、前記虚像除去手段は、前記補正後の画像の、前記傾斜磁場の強度分布が極値を取る領域に対応する画素より外側の画素の画素値を所定の値に置き換えることにより、前記虚像を除去することを特徴とする磁気共鳴イメージング装置を提供する。
【発明の効果】
【０００９】
本発明によれば、傾斜磁場分布が極値点を持つために歪み補正により発生する虚像を除去することができ、虚像のない歪み補正画像を得ることができる。
【発明を実施するための最良の形態】
【００１０】
＜＜第一の実施形態＞＞
以下、本発明を適用した第一の実施形態を図面を用いて説明する。なお、発明の実施形態を説明するための全図において、同一機能を有するものは同一符号を付け、その繰り返しの説明は省略する。
【００１１】
まず、本実施形態のＭＲＩ装置の全体構成について説明する。図１は、本実施形態のＭＲＩ装置１００の全体構成を示すブロック図である。本実施形態のＭＲＩ装置１００は、ＮＭＲ現象を利用して被検体の断層画像を得るもので、静磁場発生系２と、傾斜磁場発生系３と、シーケンサ４と、送信系５と、受信系６と、信号処理系７とを備える。
【００１２】
静磁場発生系２は、垂直磁場方式であれば、被検体１の周りの空間にその体軸と直交する方向に、水平磁場方式であれば、被検体１の体軸方向に均一な静磁場を発生させる。被検体１の周りに配置される永久磁石方式、常電導方式あるいは超電導方式の静磁場発生源により実現される。
【００１３】
傾斜磁場発生系３は、ＭＲＩ装置１００の座標系（静止座標系）であるＸ、Ｙ、Ｚの３軸方向に巻かれた傾斜磁場コイル９と、それぞれの傾斜磁場コイルを駆動する傾斜磁場電源１０とを備える。後述のシ−ケンサ４からの命令に従ってそれぞれの傾斜磁場コイルの傾斜磁場電源１０を駆動することにより、Ｘ、Ｙ、Ｚの３軸方向に傾斜磁場Ｇｘ，Ｇｙ，Ｇｚを印加する。撮影時には、スライス面（撮影断面）に直交する方向にスライス方向傾斜磁場パルス（Ｇｓ）を印加して被検体１に対するスライス面を設定し、そのスライス面に直交して且つ互いに直交する残りの２つの方向に位相エンコード方向傾斜磁場パルス（Ｇｐ）と周波数エンコード方向傾斜磁場パルス（Ｇｒ）とを印加して、エコー信号にそれぞれの方向の位置情報をエンコードする。
【００１４】
送信系５は、被検体１の生体組織を構成する原子の原子核スピンに核磁気共鳴を起こさせるために、被検体１に高周波磁場（ＲＦ）パルスを照射するもので、高周波発振器１１と変調器１２と高周波増幅器１３と送信側の高周波コイル（送信コイル）１４ａとを備える。高周波発振器１１から出力されたＲＦパルスは後述するシーケンサ４からの指令によるタイミングで変調器１２により振幅変調され、高周波増幅器１３で増幅され、被検体１に近接して配置された送信コイル１４ａから被検体１に照射される。
【００１５】
受信系６は、被検体１の生体組織を構成する原子の原子核スピンの核磁気共鳴により放出されるエコー信号（ＮＭＲ信号）を検出するもので、受信側の高周波コイル（受信コイル）１４ｂと信号増幅器１５と直交位相検波器１６と、Ａ／Ｄ変換器１７とを備える。送信コイル１４ａから照射されたＲＦパルスによって誘起された応答のＮＭＲ信号は、被検体１に近接して配置された受信コイル１４ｂで検出され、信号増幅器１５で増幅され、後述するシーケンサ４からの指令によるタイミングで直交位相検波器１６により直交する二系統の信号に分割され、それぞれがＡ／Ｄ変換器１７でディジタル量に変換されて、計測データとして信号処理系７に送られる。
【００１６】
シーケンサ４は、所定のパルスシーケンスに従って、ＲＦパルスと傾斜磁場パルスとを繰り返し印加するよう制御するもので、後述する信号処理系７のディジタル信号処理装置８の制御で動作し、被検体１の断層画像を再構成するためのデータ収集に必要な種々の命令を送信系５、傾斜磁場発生系３、および受信系６に送る。
【００１７】
信号処理系７は、各種データ処理と処理結果の表示及び保存等を行うもので、ディジタル信号処理装置８と、ＲＯＭ、ＲＡＭ等の記憶装置１８と、光ディスク、磁気ディスク等の外部記憶装置１９と、表示装置２０とを備える。受信系６からの計測データがディジタル信号処理装置（ＣＰＵ）８に入力されると、ディジタル信号処理装置８が信号処理、画像再構成等の処理を実行し、その結果である被検体１の断層画像を表示装置２０に表示すると共に、記憶装置１８または外部記憶装置１９に記録する。
【００１８】
操作部２５は、ＭＲＩ装置１００自体の各種制御情報および信号処理系７で行う処理の各種制御情報の入力を受け付けるもので、トラックボール又はマウス２３、および、キーボード２４を備える。操作部２５は表示装置２０に近接して配置され、オペレータは、表示装置２０を見ながら操作部２５を介してインタラクティブにＭＲＩ装置１００の各種処理に必要な情報を入力する。
【００１９】
なお、図１において、送信側の高周波コイル１４ａと傾斜磁場コイル９は、被検体１が挿入される静磁場発生系２の静磁場空間内に、垂直磁場方式であれば被検体１に対向して、水平磁場方式であれば被検体１を取り囲むようにして設置されている。また、受信側の高周波コイル１４ｂは、被検体１に対向して、或いは取り囲むように設置されている。
【００２０】
現在ＭＲＩ装置の撮像対象核種は、臨床で普及しているものとしては、被検体の主たる構成物質である水素原子核（プロトン）である。プロトン密度の空間分布や、励起状態の緩和時間の空間分布に関する情報を画像化することで、人体頭部、腹部、四肢等の形態または、機能を２次元もしくは３次元的に撮像する。
【００２１】
なお、上記のＭＲＩ装置１００で得られた画像は、傾斜磁場コイル９に生じる非線形性の影響を受け、歪みが生じる。本実施形態では、非特許文献１に記載の手法を用い、再構成後の画像の座標から、傾斜磁場分布の歪みを推定し、歪み補正を行う。さらに、本実施形態では、歪み補正後の画像から、傾斜磁場分布の極値点に対応する領域から得られた画像上の座標（画像上の極値点と呼ぶ。）を特定し、傾斜磁場の傾きが極値を有するために発生する虚像を除去する。これを実現するため、本実施形態のディジタル信号処理装置８は、画像再構成部に加え、歪み補正部と、虚像除去部とを備える。
【００２２】
本実施形態のディジタル信号処理装置８は、ＣＰＵとメモリとを備え、ＣＰＵ８が、予め記憶装置１８等に格納されたプログラムをメモリにロードし、実行することにより、上記各機能を実現する。
【００２３】
本実施形態の歪み補正部は、上記非特許文献１に開示の歪み補正処理を行う。具体的には、以下の処理を行う。半径Ｒ０の球面の内部で、Ｘ軸、Ｙ軸、Ｚ軸の傾斜磁場コイルが発生させる、
【数１】
0005186698.tif 000002
はそれぞれ、Ｘ軸の傾斜磁場、Ｙ軸の傾斜磁場、Ｚ軸の傾斜磁場をルジャンドル関数で展開したときの、展開係数である。
【００２４】
各軸に印加した傾斜磁場強度(ＧＸ、ＧＹ、ＧＺ)で、傾斜磁場を割ると、各軸の位置座標が得られる。それぞれの傾斜磁場が、線形成分ＧＸＸ、ＧＹＹ、ＧＺＺを含んでいることに注意すると、傾斜磁場の非線形性による画像歪みδＸ、δＹ、δＺは、次の式（２）で得られる。
【数２】
0005186698.tif 000003
従って、画像上の各点を、(−δＸ、−δＹ、−δＺ)だけ移動することで、傾斜磁場の非線形による画像歪みを補正することができる。
【００２５】
次に、本実施形態の虚像除去部による処理の詳細を説明する。虚像除去部は、上述のように、傾斜磁場分布が極値を持つために、非特許文献１に開示の手法による歪み補正（以後、単に歪み補正と呼ぶ。）後の画像に発生する虚像を除去する。ここで、本実施形態では、傾斜磁場は、上述のように、スライス方向傾斜磁場パルスＧｓ、位相エンコード方向傾斜磁場パルスＧｐ、周波数エンコード方向傾斜磁場パルスＧｒの３方向に与えられる。ここでは、これらの傾斜磁場パルスによる傾斜磁場強度の分布を合成したものを傾斜磁場分布と呼ぶ。
【００２６】
歪み補正では、傾斜磁場強度により傾斜磁場の位置座標を算出している。従って、上述のように傾斜磁場分布に極値点があると、その両側で傾斜磁場強度が同じになる点が存在し、歪み補正時に原画像上の１つの点が２箇所に配置される。一般に、傾斜磁場分布の非線形性は、磁場中心から遠ざかるにつれて高くなる。このため、極値を持つほど傾斜磁場分布の非線形性が高くなるのは、撮影範囲の端部に近い領域である可能性が高い。すなわち、歪み補正により、極値点の内側の信号を基に、極値点の外側の、傾斜磁場強度が同じ領域に画像が作り出される。これが、歪み補正による虚像である。
【００２７】
本実施形態では、歪み補正後の再構成画像上の極値点を探し出し、その外側の領域の画素値を所定の値（例えば、０等）にし、虚像を除去する。すなわち、本実施形態の虚像除去部は、画像上の極値点を探索し、外側の領域を外側領域として特定し、外側領域の画素値を所定の値とする虚像除去処理を行う。以下、本実施形態の虚像除去処理の具体的な実現手法について説明する。
【００２８】
まず、画像上の極値点の探索手法について説明する。本実施形態では、歪み補正時の変換処理の特性を利用し、画像上の極値点を特定する。非線形な傾斜磁場分布環境で得られた計測データから再構成された画像は、歪んだ格子状でない座標系（ｘ，ｙ，ｚ）に配置される。以後、この座標系を第一の座標系と呼ぶ。また、理想的な傾斜磁場分布環境で得られる計測データが再構成されて配置される座標系（Ｘ，Ｙ，Ｚ）を、第二の座標系と呼ぶ。歪み補正では、第一の座標系を第二の座標系に変換する。このとき、歪み補正の性質上、画像上の極値点では、座標変換に伴う面積素（３次元画像の場合、体積素、以下、面積素で代表して記載する。）が０となる。すなわち、非特許文献１に記載されている歪み補正による座標変換では、画像上の極値点Ａを含む四角形(Ａ、Ｂ、Ｃ、Ｄ)のＢ、Ｃ、Ｄが、補正時に極値点Ａに移動する。従って、Ａ（ｘａ、ｙａ）、Ｄ（ｘｄ、ｙｄ）とすると、式（３）で表される四角形の面積Ｓは０となる。
【数３】
0005186698.tif 000004
よって、傾斜磁場強度の極値点に対応する画素の変換後の面積は０となる。本実施形態は、これを利用して、画像上の極値点を特定する。
【００２９】
座標変換に伴う面積素の変化率は、ヤコビアンで表される。本実施形態の第一の座標系から第二の座標系への変換におけるヤコビアンは、以下の式（４）で求められる。
【数４】
0005186698.tif 000005
【００３０】
従って、本実施形態では、上記式（４）で計算されるヤコビアンＪが０になる座標を計算し、それを画像上の極値点とする。なお、実際の処理においては、ヤコビアンＪが０となる座標を画像上の極値点とするのではなく、所定の閾値を設定し、ヤコビアンＪが当該閾値以下となる座標を、画像上の極値点とする。
【００３１】
次に、求めた画像上の極値点から、画素値を所定の値にする領域、すなわち、外側の領域を決定する手法について説明する。ここでは、求めた画像上の極値点から上下方向および左右方向の画像端までの距離を算出し、短い側を外側とする。具体的には、算出された画像上の極値点を中心として、極値点と補正画像の領域の端までの距離を、上下方向および左右方向について求める。このとき、算出する画像は、歪み補正後の画像を用い、領域の端とは、歪み補正処理前の画像の端のことである。
【００３２】
以下、本実施形態の虚像除去部による虚像除去処理全体について、その流れを説明する。図２は、本実施形態の虚像除去部による虚像除去処理の処理フローである。ここでは、ｎ×ｍ（ｎ、ｍは共に自然数）の画素からなる２次元画像について処理を行うものとする。また、各画素の座標を、（ｘｉ、ｙｊ）（ｉ＝１、２、・・・ｎ、ｊ＝１、２、・・・・ｍ）とする。また、画像上の極値点の外側と判定された領域の画素値を０とする場合を例にあげて説明する。
【００３３】
ｉ、ｊをそれぞれ初期化する。ここでは、１とする（ステップＳ２０１）。次に、座標（ｘｉ、ｙｊ）のヤコビアンＪを計算する（ステップＳ２０２）。計算したヤコビアンＪの値が予め定めた閾値（例えば、０．０２）未満であるか否か、判別する（ステップＳ２０３）。ヤコビアンＪが閾値以上であれば、ステップＳ２１３へ移行する。
【００３４】
ヤコビアンＪが閾値未満であれば、その座標（ｘｉ、ｙｉ）は、傾斜磁場分布の極値点に対応する座標と判断し、当該座標から、画像の左端までの距離Ｌｌ（＝｜ｘ１−ｘｉ｜）および右端までの距離Ｌｒ（＝｜ｘｎ−ｘｉ｜）を計算する（ステップＳ２０４）。次に、当該座標から画像の上端までの距離Ｌｕ（＝｜ｙ１−ｙｊ｜）および下端までの距離Ｌｄ（＝｜ｙｍ−ｙｊ｜）を計算する（ステップＳ２０５）。
【００３５】
そして、ＬｌとＬｒとの大きさ、および、ＬｕとＬｄとの大きさを比較し（ステップＳ２０６、Ｓ２０７，Ｓ２０８）、ＬｌがＬｒ以上、かつ、ＬｕがＬｄ以上の場合、座標（ｘｉ、ｙｉ）より右下の領域の画素値を０とする（ステップＳ２０９）。また、ＬｌがＬｒ以上、かつ、ＬｕがＬｄ未満の場合、座標（ｘｉ、ｙｉ）より右上の領域の画素値を０とする（ステップＳ２１０）。ＬｌがＬｒ未満、かつ、ＬｕがＬｄ以上の場合、座標（ｘｉ、ｙｉ）より左下の領域の画素値を０とする（ステップＳ２１１）。ＬｌがＬｒ未満、かつ、ＬｕがＬｄ未満の場合、座標（ｘｉ、ｙｉ）より左上の領域の画素値を０とする（ステップＳ２１２）。
【００３６】
次に、画像上の全座標の処理が終了したか否かを判別する（ステップＳ２１３）。ここでは、ｉ＝ｎであるか否かを判別し（ステップＳ２１４）、ｉ＜ｎであれば、ｉを１インクリメントし（ステップＳ２１５）、ステップＳ２０２に戻る。一方、ｉ＝ｎであれば、ｊ＝ｍであるか否かを判別する（ステップＳ２１６）。ｊ＜ｍであれば、ｊを１インクリメントし（ステップＳ２１７）、ステップＳ２０２に戻る。ｊ＝ｍであれば、全座標の処理が終了しているとし、虚像除去処理を終了する。
【００３７】
以上の手順で虚像除去処理を行うことにより、本実施形態の虚像除去部は、歪み補正により発生した虚像を除去することができる。
【００３８】
なお、上記実施形態では、画像上の各座標のヤコビアンを計算して画像上の極値点を特定している。しかし、極値点を特定する処理は、これに限られない。例えば、実際に座標変換後の微小領域の面積（３次元画像の場合体積）を計算するよう構成してもよい。
【００３９】
座標（ｘ、ｙ）における微小領域（単位面積）Ａ’の座標変換後の面積Ａは、以下の式（５）により求められる。
【数５】
0005186698.tif 000006
なお、計算する点の座標（ｘ、ｙ）、各点の歪みを（δｘｉ、δｙｉ）とすると、式（５）の座標は以下の通りとなる。
ｘ０＝ｘ−δｘ０、      ｙ０＝ｙ-δｙ０
ｘ１＝ｘ＋１−δｘ１、  ｙ１＝ｙ−δｙ１
ｘ２＝ｘ＋１−δｘ２、  ｙ２＝ｙ＋１−δｙ２
ｘ３＝ｘ−δｘ３、      ｙ３＝ｙ＋１−δｙ３
ｘ４＝ｘ０、            ｙ４＝ｙ０
【００４０】
すなわち、再構成画像上の各座標について式（５）に従ってＡを計算し、結果が０（または、所定の閾値以下）になる座標を画像上の極値点とする。
【００４１】
なお、３次元画像の場合、座標（Ｘ、Ｙ、Ｚ）における微小領域（単位体積）Ｖ’の座標変換後の体積Ｖは、以下の式（６）により求められる。
【数６】
0005186698.tif 000007
なお、Ａ１は、計算する座標（Ｘ、Ｙ、Ｚ）のボクセルの上面の面積、Ａ２は、下面の面積である。それぞれ、式（２）を用いて求める。また、ｔは、上面の座標を（ｘ、ｙ、ｚ）、下面の座標を（ｘ’、ｙ’、ｚ’）とすると、以下の式（７）で算出される。
【数７】
0005186698.tif 000008
このように、３次元画像の場合、各座標について、式（６）に従ってＶを計算し、結果が０（または、所定の閾値以下）になる座標を、画像上の極値点とする。
【００４２】
以上説明したように、本実施形態によれば、傾斜磁場分布が極値点を有するために歪み補正により発生する虚像を、効果的に除去することができる。従って、傾斜磁場分布が非線形であり、撮影領域内に極値点を持つ場合であっても、虚像が除去されたひずみ補正画像を得ることができる。
【００４３】
例えば、特許文献１に開示されているような、被検体の配置スペースを大きくとるとともにオープン性を高めるために傾斜磁場コイル９を静磁場発生系２内の磁石の窪み内に配置するＭＲＩ装置の場合、傾斜磁場分布の線形性を保つことができる範囲が狭くなっている。このような構成のＭＲＩ装置には、特に有効である。
【００４４】
【特許文献１】国際公開２００６／０５７３９５号パンフレット
【００４５】
また、本実施形態では、歪み補正、画像再構成、虚像除去の各処理を、ＭＲＩ装置１００内のデジタル信号処理装置８で行うよう構成しているが、本構成に限られない。例えば、ＭＲＩ装置１００とデータ送受信可能な独立した情報処理装置を備え、当該情報処理装置でこれらの処理を行うよう構成してもよい。
【図面の簡単な説明】
【００４６】
【図１】本発明の実施形態のＭＲＩ装置の全体構成を示すブロック図である。
【図２】本発明の実施形態の虚像除去処理のフローチャートである。
【符号の説明】
【００４７】
１：被検体、２：静磁場発生系、３：傾斜磁場発生系、４：シーケンサ、５：送信系、６：受信系、７：信号処理系、８：ディジタル信号処理装置（ＣＰＵ）、９：傾斜磁場コイル、１０：傾斜磁場電源、１１：高周波発信器、１２：変調器、１３：高周波増幅器、１４ａ：高周波コイル（送信コイル）、１４ｂ：高周波コイル（受信コイル）、１５：信号増幅器、１６：直交位相検波器、１７：Ａ／Ｄ変換器、１８:記憶装置、１９：外部記憶装置、２０：表示装置、２３：トラックボール又はマウス、２４：キーボード、２５：操作部、１００：ＭＲＩ装置
【図１】
0005186698.tif 000009
【図２】
0005186698.tif 000010