(57)【要約】【課題】、静磁場不均一による位相誤差成分を補正することにより、画像歪みのない良好な画像を得ることが可能な磁気共鳴イメージング装置を実現する。【解決手段】読出し方向のみならず位相エンコード方向もフーリエ変換を行った計測データに予め各エコー時間で計測した前計測データに含む位相が０となるように位相回転し静磁場不均一による画像歪を補正する。静磁場不均一成分Sn(x,y)を算出するリファレンス計測で得た画像にFFTを実行し静磁場不均一成分を算出する（ステップ201,202）。本計測データのうち同一エコー計測時のデータのみ用い2DFFTを行い各エコー計測時の画像Sn(x,y)を作成し各エコー計測時での静磁場不均一成分Bn(x,y)を補正する（ステップ301,302）。補正データS'n(x,y)は逆FFTで計測空間へ戻し本来の位相エンコード位置に再配列しデータIF[S'n(x,y)]を加算し１計測空間データにする（ステップ303,304）。【選択図面】      図２【特許請求の範囲】
【請求項１】
少なくとも一回のＲＦ励起パルスの印加後に、位相エンコード傾斜磁場パルスを離散的に印加しながら読み出し方向傾斜磁場を反転させることにより画像形成のための全エコー信号もしくは一部のエコー信号を計測する第１パルスシーケンスを行う磁気共鳴イメージング装置において、
上記第１パルスシーケンスにおける各エコー信号の計測時間と同一の時間において画像再構成に必要なエコー信号を取得する第２パルスシーケンスを実行し、
上記第２パルスシーケンスにより計測したエコー信号を用いて、このエコー信号の計測時間毎の画像を作成し、
該画像から静磁場不均一成分を得て、
上記第１パルスシーケンスで取得したエコー信号を、そのエコー計測時間に対応する静磁場不均一成分により位相補正を行い、
上記位相補正を行った各エコー信号を計測空間に配列し、二次元フーリエ変換を行うことにより最終的な画像を得ることを特徴とする核磁気共鳴イメージング装置。
【請求項２】
請求項１記載の磁気共鳴イメージング装置において、第２パルスシーケンスでは、第１パルスシーケンスに対して位相エンコード傾斜磁場パルスを印加しない状態で第１パルスシーケンスと同一のエコー計測時間で各エコー信号を計測し、少なくとも一回のＲＦ励起パルスの印加後に実行されるパルスシーケンスを位相エンコードのオフセットを与える傾斜磁場パルスの強度を変化させつつ複数回繰り返すことにより各エコー計測時間において画像再構成に必要なエコー信号を取得することを特徴とする磁気共鳴イメージング装置。
【請求項３】
請求項１記載の磁気共鳴イメージング装置において、第１パルスシーケンスで計測されたエコー信号の位相補正は、第１パルスシーケンスで計測された全エコー信号のうち、同一エコー計測時間に計測されたエコー信号のみを用いて作成した画像の位相を、上記静磁場不均一成分を用いて回転させて行うことを特徴とする磁気共鳴イメージング装置。
【請求項４】
請求項３記載の磁気共鳴イメージング装置において、同一エコー計測時間に計測されたエコー信号を計測空間に配列する際、本来の位相エンコード位置又は低周波域のデータとみなして配列することを特徴とする磁気共鳴イメージング装置。
【請求項５】
請求項３記載の磁気共鳴イメージング装置において、位相補正を行った各エコー信号を計測空間に配列する際、位相補正を行った画像を２次元逆フーリエ変換することにより得られたエコー計測時間毎の計測空間データで、同一位相エンコード位置のデータ同士を単純加算することを特徴とする磁気共鳴イメージング装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、核磁気共鳴現象を利用する磁気共鳴イメージング装置に関する。
【背景技術】
【０００２】
磁気共鳴イメージング装置（ＭＲＩ装置）は、静磁場中に置かれた被検体に高周波磁場を印加することにより被検体に生じるエコー信号を検出し、これを信号処理し画像化する装置である。
【０００３】
このＭＲＩ装置では、Spin Echo法やGradient Echo法に代表されるようなエコー信号計測方法が多数存在する。このエコー信号計測方法の中で特に、少なくとも一回のＲＦ励起パルス照射後に、読み出し方向傾斜磁場を何度も反転させることにより、複数のエコー信号を計測する超高速撮像法がある。
【０００４】
この超高速撮像法としては、例えば、シングルショットＥＰＩ（Echo Planar Imaging）、マルチショットＥＰＩ法、Gradient Spin Echo Imaging(ＧＲＡＳＥ法)等がある。
【０００５】
図４に、上記超高速撮像法の代表例として、シングルショットＥＰＩ法のパルスシーケンスの一例を示す。このパルスシーケンスでは、高周波パルス１０１とスライス選択傾斜磁場パルス１０２とを印加した後、位相エンコードのオフセットを与えるパルス１０３と読み出し傾斜磁場パルスのオフセットを与えるパルス１０４を印加する。
【０００６】
次に、位相エンコード傾斜磁場パルス１０５を離散的に印加しながら、反転する読み出し傾斜磁場１０６の各周期内で発生するエコー信号１０１を、各々の時間範囲１０７の間にサンプリングする。
【０００７】
マルチショットＥＰＩ法では、図４に示した一連のパルスシーケンス１０９を複数回繰り返し、画像再構成に必要な全てのエコー信号を取得する。
【０００８】
ＭＲＩ装置では、一般的に、静磁場不均一性や、傾斜磁場の極性を高速に反転させることにより発生する渦電流や不完全性等により、図４に示したパルスシーケンス１０９により計測したエコー信号に位相誤差が生じる。
【０００９】
従来技術においては、上記エコー信号の位相誤差を補正するために、被験者の画像形成のための撮影（以下本計測と呼ぶ）に先立って、あらかじめ補正用のデータを取得しておき（以下前計測と呼ぶ）、この補正用データを用いて本計測でのエコー信号を補正する。
【００１０】
この補正手段として、例えば補正用データである一組のエコー信号を取得し、これらエコー信号をフーリエ変換した後に、エコー信号間の位相差を求め、本計測データを補正する方法（特許文献１）や、本計測に先立って、位相エンコード傾斜磁場を加えない状態で補正用のデータを取得しておき、この補正用データ及び本計測データをそれぞれ読み出し方向に１次元フーリエ変換し、本計測データから補正用データの位相を減算する方法（特許文献２）等がある。
【００１１】
【特許文献１】特開平８−２１５１７４号公報
【００１２】
【特許文献２】特開平５−６８６７４号公報
【発明の開示】
【発明が解決しようとする課題】
【００１３】
ところで、ＭＲＩ装置では、静磁場不均一の影響により、ＲＦ励起パルス照射からの経過時間に比例して位相誤差が生じる。このような位相誤差の含まれたエコー信号を計測空間に配置すると、図５に示すように位相エンコード方向への位相傾斜が生じる。
【００１４】
上記位相傾斜は実空間上の静磁場不均一の強度に比例するため、画像中の各画素で異なる位相傾斜が生じることになる。このような静磁場不均一による位相誤差を補正するためには、画像中の画素毎に異なる補正値が必要となる。
【００１５】
従来技術における補正法で取り扱うデータは、読み出し方向にのみ１次元フーリエ変換されたものであるため、読み出し方向の各画素に含まれる静磁場不均一による位相誤差は補正することができる。
【００１６】
しかしながら、位相エンコード方向について見ると、各画素の位相誤差を位相エンコード方向に積分した値について補正することになるため、各画素の位相を正権に補正することができない。
【００１７】
このため、画像に歪みが生じ、画質の向上が困難であった。
【００１８】
本発明の目的は、少なくとも一回のＲＦ励起パルス照射後に、読み出し方向傾斜磁場を何度も反転させることにより複数のエコー信号を計測するパルスシーケンスで計測される各エコー信号に含まれる位相誤差のうち、静磁場不均一による位相誤差成分を補正することにより、画像歪みのない良好な画像を得ることが可能な磁気共鳴イメージング装置を実現することである。
【課題を解決するための手段】
【００１９】
上記目的を達成するため、本発明は次のように構成される。  
（１）少なくとも一回のＲＦ励起パルスの印加後に、位相エンコード傾斜磁場パルスを離散的に印加しながら読み出し方向傾斜磁場を反転させることにより画像形成のための全エコー信号もしくは一部のエコー信号を計測する第１パルスシーケンスを行う磁気共鳴イメージング装置において、上記第１パルスシーケンスにおける各エコー信号の計測時間と同一の時間において画像再構成に必要なエコー信号を取得する第２パルスシーケンスを実行し、上記第２パルスシーケンスにより計測したエコー信号を用いて、このエコー信号の計測時間毎の画像を作成し、
該画像から静磁場不均一成分を得て、上記第１パルスシーケンスで取得したエコー信号を、そのエコー計測時間に対応する静磁場不均一成分により位相補正を行い、上記位相補正を行った各エコー信号を計測空間に配列し、二次元フーリエ変換を行うことにより最終的な画像を得る。
【００２０】
（２）好ましくは、上記（１）において、第２パルスシーケンスでは、第１パルスシーケンスに対して位相エンコード傾斜磁場パルスを印加しない状態で第１パルスシーケンスと同一のエコー計測時間で各エコー信号を計測し、少なくとも一回のＲＦ励起パルスの印加後に実行されるパルスシーケンスを位相エンコードのオフセットを与える傾斜磁場パルスの強度を変化させつつ複数回繰り返すことにより各エコー計測時間において画像再構成に必要なエコー信号を取得する。
【００２１】
（３）また、好ましくは、上記（１）において、第１パルスシーケンスで計測されたエコー信号の位相補正は、第１パルスシーケンスで計測された全エコー信号のうち、同一エコー計測時間に計測されたエコー信号のみを用いて作成した画像の位相を、上記静磁場不均一成分を用いて回転させて行う。
【００２２】
（４）また、好ましくは、上記（３）において、同一エコー計測時間に計測されたエコー信号を計測空間に配列する際、本来の位相エンコード位置又は低周波域のデータとみなして配列する。
【００２３】
（５）また、好ましくは、上記（３）において、位相補正を行った各エコー信号を計測空間に配列する際、位相補正を行った画像を２次元逆フーリエ変換することにより得られたエコー計測時間毎の計測空間データで、同一位相エンコード位置のデータ同士を単純加算する。
【００２４】
本発明によれば、読み出し方向のみならず、位相エンコード方向にもフーリエ変換を行った計測データに対して、予め各エコー時間で計測しておいた前計測データに含まれる位相が０となるように位相回転することによって、静磁場不均一による画像歪みを補正したので、静磁場不均一による画像歪みを抑制した画像を得ることができる。
【発明の効果】
【００２５】
本発明によれば、少なくとも一回のＲＦ励起パルス照射後に、読み出し方向傾斜磁場を何度も反転させることにより複数のエコー信号を計測するパルスシーケンスで計測される各エコー信号に含まれる位相誤差のうち、静磁場不均一による位相誤差成分を補正し、画像歪みのない良好な画像を得ることが可能な磁気共鳴イメージング装置を実現することができる。
【発明を実施するための最良の形態】
【００２６】
まず、本発明の原理について、以下に説明する。
【００２７】
図４に示したパルスシーケンス１０９において計測される静磁場不均一Ｂ（ｘ，ｙ）の影響を含んだＭＲＩ信号ｓ（ｎｘ、ｎｙ）は、実空間上の読み出し方向をｘ、位相エンコード方向をｙ、エコー時間をｔeとすると、下記の式（１）で表わせる。
【００２８】
【数１】
2005028019.tif 000002
【００２９】
ただし、上記式（１）において、（ｎｘ、ｎｙ）は計測空間上の点を示し、ρ（ｘ，ｙ）はプロトン密度、ωは共鳴周波数、Ｇｙは読み出し傾斜磁場の印加強度、Ｇｙは位相エンコード傾斜磁場の印加強度、△ｔはデータサンプル時間、△τは位相エンコード傾斜磁場パルスの印加時間である。
【００３０】
ここで、上記式（１）は下式（２）のように示すことができる。
【００３１】
【数２】
2005028019.tif 000003
【００３２】
このとき、ｓ（ｎｘ、ｎｙ）を２次元フーリエ変換（Ｆ｛ｓ（x,y）｝）した画像Ｓ（ｘ，ｙ）は、下式（３）となる。
【００３３】
【数３】
2005028019.tif 000004
【００３４】
なお、上記式（３）は、静磁場不均一の影響が、画素中の位相に変換されることを示している。
【００３５】
これにより、図４に示したパルスシーケンス１０９での各エコー時間（ｔe1，ｔe2，ｔe3，…、ｔen）における静磁場不均一の位相情報は、それぞれ、e[φ1]、e[φ2]、・・・e[φn]であることが分かる。
【００３６】
静磁場不均一による画像歪みは、これらエコー毎の位相情報が異なる時に発生する。例えば、ある画素（ｘｐ，ｙｐ）での位相e[φ1]、e[φ2]、・・・e[φn]が次第に増加、もしくは減少している場合、計測空間に各エコーを配列すると、図５に示すような位相エンコード方向への位相傾斜となる。
【００３７】
このため、画像再構成を行うと、画素（ｘｐ，ｙｐ）は、位相エンコード方向に位相傾斜に応じた位置移動が生じる。
【００３８】
従来技術における補正方法で取り扱うデータは、読み出し方向にのみ１次元フーリエ変換されたものであるため、読み出し方向の各画素に含まれる静磁場不均一による位相誤差は補正することができるが、位相エンコード方向について見ると、各画素の位相誤差を位相エンコード方向に積分した値について補正することになるため、各画素の位相を正確に補正することができない。
【００３９】
以上のことから、読み出し方向のみならず、位相エンコード方向にもフーリエ変換を行った計測データに対して、予め各エコー時間で計測しておいた前計測データに含まれる位相e[φn]、（ｔen＝te1,te2, te3,・・・，ｔen）が０となるように位相回転することによって、静磁場不均一による画像歪みを補正できる。
【００４０】
次に、上記原理を用いるＭＲＩ装置について説明する。  
図１は、本発明が適用されるＭＲＩ装置の全体概略構成図である。
【００４１】
図１において、ＭＲＩ装置は、静磁場発生回路１と、傾斜磁場発生系２と、送信系３と、受信系４と、信号処理系５と、シーケンサ６と、中央処理装置（ＣＰＵ）７と、操作部８とを備える。
【００４２】
静磁場発生回路１は、被検体９の周りにその体軸方向または体軸と直交する方向に均一な静磁場を発生させる。また、静磁場発生回路１は、被検体９の周りのある広がりをもった空間に永久磁石方式又は常電導方式あるいは超電導方式の磁場発生手段が配置されている。
【００４３】
傾斜磁場発生系２は、Ｘ、Ｙ、Ｚの三軸方向に巻かれた傾斜磁場コイル１０と、それぞれのコイルを駆動する傾斜磁場電源１１とを備える。また、傾斜磁場発生系２は、シーケンサ６から命令に従って、それぞれのコイルの傾斜磁場電源１１を駆動することにより、Ｘ、Ｙ、Ｚの三軸方向の傾斜磁場Ｇｓ、Ｇｐ、Ｇｆを被検体９に印加する。
【００４４】
この傾斜磁場の印加方法により、被検体９に対するスライス面を設定することができる。
【００４５】
送信系３は、シーケンサ６から送出される高周波磁場パルスにより被検体９の生体組織を構成する原子の原子核に核磁気共鳴を起こさせるために高周波信号を照射するもので、高周波発振器１２と、変調器１３と、高周波増幅器１４と、送信側の高周波コイル１５とを備える。
【００４６】
高周波発振器１２から出力された高周波パルスを変調器１３を介して高周波増幅器１４で増幅した後に、被検体９に近接して配置された送信側の高周波コイル１５に供給することにより、電磁波が被検体９に照射される。
【００４７】
受信系４は、被検体９の生体組織の原子核の核磁気共鳴により放出されるエコー信号（ＮＭＲ信号）を検出する。この受信系４は、受信側の高周波コイル１６と、増幅器１７と、直交位相検波器１８と、Ａ／Ｄ変換器１９とを備える。
【００４８】
送信側の高周波コイル１５から被検体９に電磁波が照射されると、この電磁波の被検体９からの応答電磁波（ＮＭＲ信号）は、被検体９に近接して配置された受信側の高周波コイル１６で検出される。
【００４９】
そして、高周波コイル１６で検出された電磁波は、増幅器１７及び直交位相検波器１８を介してＡ／Ｄ変換器１９に入力され、ディジタル量に変換される。このＡ／Ｄ変換器１９において、シーケンサ６からの命令によるタイミングで直交位相検波器１８によりサンプリングされた二系列の収集データとされ、その信号が、Ａ／Ｄ変換器１９から信号処理系５に送られる。
【００５０】
信号処理系５は、上記受信系４で検出したエコー信号を用いて画像再構成演算を行うと共に画像表示をする。この信号処理系５は、ＣＰＵ７と、ＲＯＭ（読み出し専用メモリ）２０と、ＲＡＭ（随時書き込み読み出しメモリ）２１と、光磁気ディスク２２及び磁気ディスク２４と、ディスプレイ２３とを備える。
【００５１】
ＣＰＵ７は、上記エコー信号についてフーリエ変換、補正係数計算、画像再構成等の処理及びシーケンサ６の制御を行う。また、ＲＯＭ２０は、経時的な画像解析処理及び計測を行うプログラムやその実行において用いる不変のパラメータなどを記憶する。
【００５２】
また、ＲＡＭ２１は、前計測で得た計測パラメータや上記受信系４で検出したエコー信号、及び関心領域設定に用いる画像を一時保管すると共にその関心領域を設定するためのパラメータなどを記憶する。
【００５３】
また、光磁気ディスク２２及び磁気ディスク２４は、ＣＰＵ７で再構成された画像データを記録するデータ格納部である。そして、ディスプレイ２３は、光磁気ディスク２２又は磁気ディスク２４から読み出した画像データを映像化して断層像として表示する。
【００５４】
シーケンサ６は、被検体９の生体組織を構成する原子の原子核に核磁気共鳴を起こさせる高周波磁場パルスをある所定のパルスシーケンスで繰り返し印加する制御手段である。このシーケンサ６は、ＣＰＵ７の制御で動作し、被検体９の断層像のデータ収集に必要な種々の命令を送信系３及び傾斜磁場発生系２並びに受信系４に送る。
【００５５】
また、操作部８は、信号処理系５で行う処理の制御情報を入力するものであり、マウス２５と、キーボード２６とを備える。
【００５６】
次に、本発明の実施形態におけるパルスシーケンスについて説明する。なお、本発明は、少なくとも一回のＲＦ励起パルス信号の印加後に、位相エンコード傾斜磁場パルス信号を離散的に印加しながら読み出し方向傾斜磁場を何度も反転させることにより画像形成のための全エコー信号、もしくは一部のエコー信号を計測する全てのパルスシーケンスに適用可能である。ただし、以下の実施形態は、本発明をＧＥ-ＥＰＩシーケンスに適用した場合を例として説明する。
【００５７】
図２は、本発明の一実施形態における画像作成の動作フローチャートである。なお、本計測の一例として、図４を用いて説明したパルスシーケンスを使用する。
【００５８】
本計測（第１パルスシーケンス）により計測したデータの静磁場不均一の影響を除去するためには、図４中、各エコーの計測時（ｔe1，ｔe2，ｔe3，…、ｔen）の静磁場不均一成分Ｂn（ｘ，ｙ）を算出する必要がある。
上記各エコーの計測時（ｔe1，ｔe2，ｔe3，…、ｔen）での静磁場不均一成分Ｂn（ｘ，ｙ）を算出するために使用するパルスシーケンスの例を図３に示す。以下この計測をリファレンス計測（第２パルスシーケンス）と呼ぶ。
【００５９】
このパルスシーケンスでは、高周波パルス１０１とスライス選択傾斜磁場パルス１０２とを印加した後、位相エンコードのオフセットを与えるパルス１０３と読み出し傾斜磁場パルスのオフセットを与えるパルス１０４とを印加する。
【００６０】
次に、反転する読み出し傾斜磁場１０６の各周期内で発生するエコー信号１０１を、各々の時間範囲１０７の間サンプリングする。すなわち、リファレンス計測では、図３に示した一連のパルスシーケンス１０９において計測されるエコー信号１０８は、同一の位相エンコード量が印加されたデータである。
【００６１】
リファレンス計測では、図３に示した一連のパルスシーケンス１０９を位相エンコードのオフセットを与えるパルス１０３の強度を変化させつつ、複数回繰り返し、各エコーの計測時（ｔe1，ｔe2，ｔe3，…、ｔen）において画像再構成に必要なエコー信号を取得する。
【００６２】
この際、正確な静磁場不均一成分Ｂn（ｘ，ｙ）を求めるためには、本計測で印加される全ての位相エンコードを印加したエコー信号を、各エコーの計測時（ｔe1，ｔe2，ｔe3，…、ｔen）において計測しなければならない。
【００６３】
しかしながら、このような計測は、長い計測時間を必要とするため、実用上には間題がある。
【００６４】
そこで、本発明に使用される静磁場不均一成分Ｂn（ｘ，ｙ）の代わりに、本計測で印加される低周波域の位相エンコードのみ印加したエコー信号にて算出できる静磁場不均一の低周波成分Ｂｏn（ｘ，ｙ）を使用しても良い。
【００６５】
なお、リファレンス計測は、本計測前に実行しても、本計測後に実行しても良い。
【００６６】
リファレンス計測実行後、各エコーの計測時（ｔe1，ｔe2，ｔe3，…、ｔen）における静磁場不均一成分Ｂn（ｘ，ｙ）を算出する。
【００６７】
リファレンス計測で得られたデータは、図２に示す、二次元（２Ｄ）高速フーリエ変換（以下ＦＦＴ）を行うステップ２０１によって画像化する。
【００６８】
各エコーの計測時（ｔe1，ｔe2，ｔe3，…、ｔen）における静磁場不均一成分Ｂn（ｘ，ｙ）は、各エコーの計測時毎のリファレンス計測データに対して、２ＤＦＦＴを施すことにより作成した画像データＲn（ｘ，ｙ）を用いて、次式（４）で算出できる（ステップ２０２）。
【００６９】
Bn(x,y)=atan｛im(Rn(x,y))/re(Rn(x,y))｝      −−−（４）
なお、ｒｅ（）信号の実部を示し、ｉｍ（）は虚部を示す。また、ａｔａｎ（）は、ａｒｃｔａｎｇｅｎｔ関数である。
【００７０】
次に、本計測実行後、各エコーの計測時（ｔe1，ｔe2，ｔe3，…、ｔen）における静磁場不均一成分Ｂn（ｘ，ｙ）を補正する処理を行う（ステップ３０２に進む）。
【００７１】
ステップ３０２の処理前に、本計測データのうち、同一のエコー計測時に計測されたデータのみを用いて２ＤＦＦＴを行うことによって、各エコー計測時（ｔe1，ｔe2，ｔe3，…、ｔen）の画像Ｓn（ｘ，ｙ）｛（ｎ＝１，２，３、…、Ｎ）｝を作成する（ステップ３０１）。
【００７２】
なお、上記同一のエコー計測時に計測されたデータは、本来の位相エンコード位置に配列しても良いし、低周波域のデータとみなして配列しても良い。
【００７３】
上記エコー配列の際、他のエコー計測時に計測されたデータの領域には０詰めを行う。
【００７４】
以上の動作により作成された各エコー計測時（ｔe1，ｔe2，ｔe3，…、ｔen）の画像Ｓn（ｘ，ｙ）｛（ｎ＝１，２，３、…、Ｎ）｝は、リファレンス計測データにより作成した同一エコー計測時の静磁場不均一成分Ｂn（ｘ，ｙ）を用いて、次式（５−１）、（５−２）に従って位相回転することにより、静磁場不均一による位相誤差を補正することができる（ステップ３０２）。
【００７５】
re｛S'n(x,y)｝=re｛Sn(x,y)｝・cos｛Bn(x,y)｝+im｛Sn(x,y)｝・sin｛Bn(x,y)｝
−−−（５−１）
im｛S'n(x,y)｝=im｛Sn(x,y)｝・cos｛Bn(x,y)｝-re｛Sn(x,y)｝・sin｛Bn(x,y)｝
−−−（５−２）
ただし、上記（５−１）、（５−２）式において、sin（）はsine関数であり、cos（）はcosine関数である。
【００７６】
上記（５−１）、（５−２）式により補正したデータＳ’n（ｘ，ｙ）｛（ｎ＝１，２，３、…、Ｎ）｝は、逆ＦＦＴ（２ＤＩＦＦＴ）を行うステップ３０３により計測空間へ戻し、本来の位相エンコード位置に再配列する。
【００７７】
この際、ステップ３０３により逆ＦＦＴされたデータＩＦ[Ｓ’n（ｘ，ｙ）]では、ステップ３０１の処理時に計測空間に配列した領域以外の領域にもデータが発生するが、これは、（５−１）、（５−２）式に従い実空間上で位相回転を行ったことによりデータの位置移動が生じている為である。
【００７８】
このため、逆ＦＦＴを行ったデータＩＦ[Ｓ’n（ｘ，ｙ）]を単純加算して１つの計測空間データにするステップ３０４を行う。
【００７９】
以上の処理により作成した計測空間上のデータに対してＦＦＴを行い（ステップ３０５）、静磁場不均一による画像歪みを抑制した画像を得ることができる。
【００８０】
つまり、本発明の一実施形態によれば、読み出し方向のみならず、位相エンコード方向にもフーリエ変換を行った計測データに対して、予め各エコー時間で計測しておいた前計測データに含まれる位相e[φn]、（ｔen＝te1,te2, te3,・・・，ｔen）が０となるように位相回転することによって、静磁場不均一による画像歪みを補正したので、静磁場不均一による画像歪みを抑制した画像を得ることができる。
【００８１】
なお、上述した例においては、静磁場の不均一を上記式（４）により算出するように構成したが、静磁場不均一検出手段により静磁場不均一を検出してもよい。
【図面の簡単な説明】
【００８２】
【図１】本発明に適用されるＭＲＩ装置の概略構成図である。
【図２】本発明の一実施形態における信号処理のフローチャートである。
【図３】本発明の一実施形態におけるリファレンス計測のパルスシーケンスを示す図である。
【図４】シングルショットＥＰＩパルスシーケンスの一例を示す図である。
【図５】静磁場不均一の影響による計測空間上での位相誤差の説明図である。
【符号の説明】
【００８３】
１            静磁場発生回路
２            傾斜磁場発生系
３            送信系
４            受信系
５            信号処理系
６            シーケンサ
７            ＣＰＵ
８            操作部
９            被検体
１０          傾斜磁場コイル
１１          傾斜磁場電源
１２          高周波発振器
１３          変調器
１４、１７    高周波増幅器
１５          送信用高周波照射コイル
１６          受信用高周波照射コイル
１８          直交位相検波器
１９          Ａ／Ｄ変換器
２０          ＲＯＭ
２１          ＲＡＭ
２２          光磁気ディスク
２３          ディスプレイ
２４          磁気ディスク
２５          トラックボール又はマウス
２６          キーボード
【図１】
2005028019.tif 000005
【図２】
2005028019.tif 000006
【図３】
2005028019.tif 000007
【図４】
2005028019.tif 000008
【図５】
2005028019.tif 000009