"検査対象の核磁気共鳴により発生するエコー信号を取得する撮像部と、前記撮像部が取得したエコー信号を処理し、核磁気共鳴画像を作成する画像処理部とを備えた磁気共鳴イメージング装置において、前記画像処理部が処理するエコー信号は、第1の計測で取得した第1計測データと、第2の計測で取得した第2計測データとを含み、且つ第1計測データ及び第2計測データは、それぞれ、異なる極性の読出し傾斜磁場を用いて取得した第1群のエコーデータと第2群のエコーデータを含み、前記画像処理部は、前記第1計測データを用いて前記第2計測データを補正する補正手段を備え、前記補正手段は、前記第1計測データの第1群のエコーデータと第2群のエコーデータとの位相差分データを周波数方向に連続する位相差分データに変換して補正用データを作成し、その際、信号強度が閾値以下の周波数方向のデータ点または領域のみを予め求めたフィッティングデータで置換することを特徴とする磁気共鳴イメージング装置。請求項1に記載の磁気共鳴イメージング装置において、前記補正手段は、前記第1群のエコーデータと第2群のエコーデータとの位相差分データに加えて、前記第1群のエコーデータ間の位相差分と、前記第2群のエコーデータ間の位相差分と、を用いて補正用データを作成するものであり、前記第1群のエコーデータと第2群のエコーデータとの位相差分は、前記第1群のエコーデータ間の位相差分を算出する際の基準となる第1群中のエコーデータと、前記第2群のエコーデータ間の位相差分を算出する際の基準となる第2群中のエコーデータと、の位相差分であることを特徴とする磁気共鳴イメージング装置。請求項1に記載の磁気共鳴イメージング装置において、前記補正手段は、第1群及び第2群のエコーデータのうち、閾値以上のデータ点を選択し、選択されなかったデータ点または領域によって隔てられたデータ点の領域島領域毎に、位相差分を周波数方向に連続にする処理を行い、隣接する島領域と不連続な点或いは領域によって隔てられた島領域について、当該島領域の端部の位相値を、前記隣接する島領域から外挿によって予測される目標値と実際の位相値との差異が最小にするように決定することを特徴とする磁気共鳴イメージング装置。請求項3に記載の磁気共鳴イメージング装置において、前記補正手段は、前記島領域が複数存在するとき、データ点数が最大である島領域を検索し、前記最大の島領域の位相差情報をフィッティングすることにより近似関数を求め、当該近似関数の傾きを用いて前記目標値を算出することを特徴とする磁気共鳴イメージング装置。請求項1に記載の磁気共鳴イメージング装置において、前記画像処理部が処理する第1計測データ及び第2計測データは、それぞれ、複数のスライスの計測データを含み、前記補正手段は、複数のスライスの計測データのうち、信号強度がスライス間で設定された閾値以下のスライスについては、第1群のエコーデータと第2群のエコーデータとの位相差分を、フィッティングデータで置換することにより補正用データを作成することを特徴とする磁気共鳴イメージング装置。請求項1に記載の磁気共鳴イメージング装置において、前記補正手段が行う位相差を周波数方向に連続する処理は、隣接する2つのデータ点の位相差を複素数に変換するステップと、前記複素数の実部と虚部の値から位相値を算出するステップと、前記隣接する2つのデータ点の一方のデータ点の位相に、前記算出された位相値を加算し、他方のデータ点の位相とするステップとを含むことを特徴とする磁気共鳴イメージング装置。位相の主値回りを生じているデータに対し、主値回りを除去する位相アンラップ方法であって、アンラップ方向の先端にあるデータ点の位相値第1位相値とそれに隣接するデータ点の位相値第2位相値の位相差が±2πを超えるときに、前記第1位相値と前記第2位相値との位相差を複素数に変換するステップと前記複素数の実部と虚部の値から位相値を算出するステップと、前記算出された位相値を前記第1位相値に加算した値を、第2位相値に置き換えるステップと、を含む位相アンラップ方法。請求項7に記載の位相アンラップ方法であって、位相アンラップに先立ち、位相値が閾値以下であるデータ点の位相値をゼロとするステップを含み、位相値がゼロであるデータ点からなる区間によって隔てられた、位相値が閾値を超えるデータ点からなる島領域のデータについて、それぞれ独立して位相アンラップ処理を行うステップと、島領域のうち最もデータ点数の大きい島領域について、位相アンラップ処理後の位相値の傾斜をフィッティングするステップと、隣り合う2つの島領域の端部のデータ点の一方について、他方のデータ点の位相値と前記フィッティングするステップで求めたフィッティング関数の傾斜を用いて位相目標値を算出するステップと、前記一方のデータ点の実際の位相値とフィッティング関数を用いて算出された位相目標値との位相差を求め、当該位相差を前記位相目標値に加算した値を実際の位相値に置換するステップとを含むことを特徴とする位相アンラップ方法。請求項8に記載の位相アンラップ方法であって、前記位相値がゼロである区間の位相値を、当該区間に隣接する島領域の端部の位相値を用いて線形補間して求めることを特徴とする位相アンラップ方法。請求項8に記載の位相アンラップ方法であって、前記位相値がゼロである区間の位相値を、前記フィッティングステップで求めたフィッティング関数を用いて決定することを特徴とする位相アンラップ方法。特許文献1に記載された技術では、奇数エコー間又は偶数エコー間の位相差分及び奇偶数エコー間の位相差分に対し、周波数方向の位相変化を直線でフィッティングし、それを位相エンコード方向に位相累積することにより、エコートレイン方向の位相変化を滑らかに接続している。この技術では、周波数方向のフィッティングに単純な一次フィッティングを採用しているが、フィッティング処理には少なからずフィッティング誤差が生じ、N2アーチファクトの原因となりうる。また奇偶数エコーそれぞれについて、位相差分データを位相エンコード方向に位相累積する際には、開始点となるエコーを基準エコーとして、順次位相累積が行われるが、これら基準エコーの信号強度が極端に低く、両者の位相差に大きな乱れが生じる場合には、ストリークアーチファクトが生じてしまう可能性がある。特許文献1の技術では、エコートレイン方向の中で最も信号強度が高いエコーを基準エコーとして採用することで極力高い信号強度が得られるような工夫がなされているが、周波数方向に対しては十分な信号強度が得られない場合がある。例えば、マルチプルアレイコイルのように空間的に不均一な感度分布を持った受信コイルを用いて撮像する場合には、局所的に信号強度が近い場所が存在するため、その領域における基準エコー間の位相差に大きな乱れが生じる可能性があり、それによりストリークアーチファクトが発生する。一方、特許文献2に記載された技術は、信号強度が低く、位相情報の精度が低いデータを位相アンラップ処理から除外するものであり、精度の高いアンラップ処理を行うことができるが、位相アンラップ処理する際に、2π×NのNを検索するステップや加算又は減算かを判定するステップ等が必須であり、位相アンラップ処理されない領域島領域と島領域との間については、評価関数を用いた評価を必要とし、複雑な処理を要する。上記課題を解決するため、本発明のMRI装置は、エコープラナー計測によって取得したエコーデータを補正するための補正用データを、前計測データの奇エコーと偶数エコーとの位相差情報をフィッティンングして求める際に、フィッティングを信号強度の比較的低いデータ領域に限定的に行い、フィッティイング誤差の発生を抑制する。また本発明のMRI装置は、奇数エコー及び偶数エコーそれぞれを位相エンコード方向に位相累積する際の基準エコーについて、奇偶数エコー間の位相差についても上記限定されたフィッティングを行い、基準エコー間の位相乱れに起因するアーチファクトの発生を抑制する。本発明の位相アンラップ処理方法は、隣り合う点の位相差を単位ベクトルに変換し、位相差として再計算する。これにより、もとの位相差に±2π×NN=0,1,2を加算して最小値を求める計算と同等な処理を簡易に行うことができる。即ち、本発明のMRI装置は、検査対象の核磁気共鳴により発生するエコー信号を取得する撮像部と、前記撮像部が取得したエコー信号を処理し、核磁気共鳴画像を作成する画像処理部とを備え、前記画像処理部が処理するエコー信号は、第1の計測前計測で取得した第1計測データと、第2の計測で取得した第2計測データとを含み、且つ第1計測データ及び第2計測データは、それぞれ、異なる極性の読出し傾斜磁場を用いて取得した第1群のエコーデータ奇数エコーデータと第2群のエコーデータ偶数エコーデータを含み、前記画像処理部は、前記第1計測データを用いて前記第2計測データを補正する補正手段を備え、前記補正手段は、前記第1計測データの第1群のエコーデータと第2群のエコーデータとの位相差分データを周波数方向に連続する位相差分データに変換して補正用データを作成し、その際、信号強度が閾値以下の周波数方向のデータ点または領域のみを予め求めたフィッティングデータで置換する。本発明のMRI装置において、前記補正手段は、第1群及び第2群のエコーデータのうち、閾値以上のデータ点を選択し、選択されなかったデータ点または領域によって隔てられたデータ点の領域島領域毎に、位相差分を周波数方向に連続にする処理を行い、隣接する島領域と不連続な点或いは領域によって隔てられた島領域について、当該島領域の端部の位相値を、前記隣接する島領域から外挿によって予測される目標値と実際の位相値との差異が最小にするように決定する。また本発明の位相アンラップ方法は、位相の主値回りを生じているデータに対し、主値回りを除去する方法であり、アンラップ方向の先端にあるデータ点の位相値第1位相値とそれに隣接するデータ点の位相値第2位相値の位相差が±2πを超えるときに、前記第1位相値と前記第2位相値との位相差を複素数に変換するステップと、前記複素数の実部と虚部の値から位相値を算出するステップと、前記算出された位相値を前記第1位相値に加算した値を、第2位相値に置き換えるステップと、を含む。本発明のMRI装置によれば、EPIデータの補正用データを前計測によって得る際に、前計測エコーデータの周波数方向の位相フィッティングを信号値の弱い領域等に限定することにより、フィッティング誤差によるストリークアーチファクトが低減されたMR画像を得ることができる。また本発明のMRI装置によれば、奇数エコーと偶数エコーの基準となるエコー間の位相差に対してもフィッティング処理を加えることで、従来よりもアーチファクトを高精度に低減した画像を得ることができる。さらに本発明によれば、MRIの基本技術である位相アンラップ方法として、精度が高く且つ従来の位相アンラップ方法に比べ簡便は位相アンラップ方法が提供される。以下、添付図面を参照して本発明のMRI装置の好ましい実施形態について説明する。図1は、本発明を適用するMRI装置の全体構成を示すブロック図である。このMRI装置は、静磁場発生系11、傾斜磁場発生系12、送信系13、受信系14、シーケンサ15、再構成演算部画像処理部16、制御系17および表示部ディスプレイ18を備え、また図示しない操作部を備えている。主として、静磁場発生系11、傾斜磁場発生系12、送信系13、受信系14及びシーケンサ15が、本発明における撮像部の機能を担っている。静磁場発生系11は、永久磁石、常電導磁石または超電導磁石のいずれかからなり、被検体10が置かれる空間に均一な静磁場を発生させる。静磁場の方向は、通常、被検体10の体軸方向または被検体10の体軸と直交する方向である。傾斜磁場発生系12は、X、Y、Zの3軸方向に巻かれた傾斜磁場コイル121と、これらの各々のコイルを磁化させる傾斜磁場電源122とからなり、シーケンサ15からの命令に従って傾斜磁場電源122の各々のコイルを磁化させることによりX、Y、Zの3軸方向の傾斜磁場Gs、Gp、Gfを被検体10に印加する。この傾斜磁場の加え方により、被検体10の撮影断面が設定される。送信系13は、高周波発振器131、変調器132、高周波増幅器133及び高周波照射コイル134とからなり、傾斜磁場発生系12で設定された被検体10の断面を構成する原子の原子核に核磁気共鳴を起こさせるために、高周波発振器131から出力されたRFパルスを高周波増幅器133で増幅した後に、被検体10に近接して設置された高周波照射コイル134に供給して被検体10に照射する。受信系14は、高周波受信コイル141、受信回路142及びアナログ/ディジタル以下A/Dという変換器143からなり、送信系13の高周波照射コイル134から照射された電磁波による被検体10の原子核の核磁気共鳴によるNMR信号であるエコー信号を、被検体10に近接して配置された高周波受信コイル141で検出し、受信回路142を介してA/D変換器143に入力し、ディジタル信号に変換して、さらにシーケンサ15からの命令によるタイミングでサンプリングされたデータ以下サンプリングデータというとして、再構成演算部16に送る。制御系17は、シーケンサ15に指令を送り、シーケンサに設定された所定のパルスシーケンスに従い、上述した傾斜磁場発生系12、送信系13および受信系14を制御するとともに、再構成演算部16の動作を制御する。パルスシーケンスは、撮像方法によって種々のものがあり、それらはプログラム化されて制御系17に格納されている。本実施形態のMRI装置では、EPI系のパルスシーケンスを備えている。再構成演算部16は、サンプリングデータにフーリエ変換などの各種演算を行なうとともに、後述する位相変化や補正用データの算出、補正用データを用いた補正などを行なう。再構成演算部16で再構成された画像は、表示部ディスプレイ18に表示される。次に上記構成におけるMRI装置の動作を説明する。図2に、動作の手順を示す。本実施の形態のMRI装置では、画像化するエコー信号を計測する本計測ステップ220と、本計測とは別に実行される前計測ステップ210を実行する。ここでは、ステップ210の計測を前計測と呼んでいるが、前計測を実行するタイミングは、本計測の前と後のいずれでもよい。本計測では、マルチエコー撮像によるパルスシーケンスを実行し、本計測によって得られたエコー信号本計測エコー30を、前計測によって得られたエコー信号前計測エコー20を用いて補正しステップ222、画像化するステップ223。本計測エコー30の補正は、本計測エコー30および前計測エコー20をそれぞれ周波数方向にフーリエ変換ステップ211、221し、x-kyデータとし、差分することにより行なう。前計測エコー20を補正用データとして補正に供するに際し、奇数エコー21と偶数エコー22に分けて、それぞれの差分位相差分を周波数方向フィッティングし、エコートレイン方向に位相累積することにより、前計測エコーに2次的に混入する位相ノイズを除去するとともに時間変動を含む位相変化を反映した補正用データを得る。以下、各ステップの詳細を説明する。<計測ステップ>本計測ステップ220及び前計測ステップ210では、シングルショット又はマルチショットEPIシーケンスにより、1回の励起で複数のエコー信号を計測する。図3及び図4に、本計測および前計測で用いられるパルスシーケンスの一例を示す。これらパルスシーケンスは公知のSpinEcho系EPIシーケンスであり、簡単に説明する。まず検知する磁化を含む被検体にRFパルス301を照射すると同時にスライスを選択する傾斜磁場パルス302を印加し、画像化するスライスを選択後、前記傾斜磁場パルスによる位相分散を収束させるためのリフォーカス傾斜磁場パルス303を印加する。次いで、スピンエコーを生成するために180度RFパルス305と2番目のスライス選択用の傾斜磁場パルス304、およびその前後に磁化の位相を分散させるためのクラッシャー傾斜磁場パルス304、307を印加する。その後、周波数エンコード傾斜磁場のオフセットを与えるパルス309を印加した後、連続して反転する周波数エンコード傾斜磁場パルス311を印加し、反転する周波数エンコード傾斜磁場311の各周期内でエコー信号313が時系列的に発生するので、これをADの時間範囲の間おのおのサンプリングし時系列データを得る。反転する周波数エンコード傾斜磁場311の強度、印加タイミング、ADの時間範囲は、本計測シーケンスと前計測と同じである。しかし、図3に示す本計測シーケンスでは、各エコー信号313に位相エンコードするための傾斜磁場パルス308、310が印加されるのに対し、図4に示す前計測シーケンスでは、位相エンコードは用いない。シングルショットシーケンスの場合には、図示する繰り返し単位314を1回で、マルチショットシーケンスであれば、繰り返し単位314を複数回繰り返して、それぞれk空間全体を埋めるデータを取得する。1回の繰り返し単位で取得する複数のエコーをエコートレインと言い、奇数番目に取得されるエコーを奇数エコー、偶数番目に取得されるエコーを偶数エコーという。このような計測をスライス傾斜磁場302306を異ならせて行い、複数のスライスのデータを取得する。スライス数は特に限定されるものではないが、本実施の形態では、少なくとも前計測において、複数のスライスのデータを取得することが望ましい。複数のスライスのデータを利用することにより、より精度の高い補正用データを作成することができる。このようなEPIを基本とする計測シーケンスでは、周波数エンコード傾斜磁場の印加タイミングおよびサンプリング時間の誤差があると、エコーのピークがk空間中央ky軸に対してずれる。ずれる方向が、奇数エコーと偶数エコーとでは反対になるため、N2アーチファクトを生じる。また位相エンコード方向のオフセットに誤差があると、エコーのピークのky軸からのずれが、エコー番号が大きくなるほど拡大する。以下の処理では、前計測エコーから、エコーのピークのずれに対応する位相情報補正用データを得て、本計測エコーのピークずれを解消する。<周波数方向フーリエ変換ステップ211、221>本計測および前計測で得られたエコー信号20、30をそれぞれ周波数方向にフーリエ変換し、x-kyデータとする。周波数方向のフーリエ変換は、計測空間kx-ky空間におけるエコー信号のピークずれを、x-ky空間上で位相補正するために行われる。即ち、計測空間におけるエコー信号のピークずれは、フーリエ変換後の空間x-ky空間での位相変化に対応するので、前計測のx-kyデータから求めた位相変化を用いて、x-ky空間上で本計測データの位相を補正することにより、本計測データのピークずれが補正される。<補正用データの作成>前計測エコーをフーリエ変換した後のx-kyデータ以下、エコーデータというを、奇数エコーと偶数エコーに分けて、それぞれについて補正用データを得る。以下、奇数エコーと偶数エコー各処理の詳細を説明する。<奇数エコーデータの処理>奇数エコー間の位相差分を求める。奇数エコー間或いは偶数エコー間で位相差分した場合、周波数エンコード傾斜磁場の印加タイミングとサンプリング時間の誤差はキャンセルされるので、位相エンコード方向のオフセット誤差のみが位相差分として求められる。この位相差分は、周波数方向に直線となるものであり、位相差分後の前計測データをフィッティングすることにより求められる。具体的には、エコートレインを構成する複数のエコーのうち、2つのエコーデータ、例えばエコートレインの1番目のエコーと3番目のエコーとの位相差分を取る。位相差分処理を行うエコーの一方は、エコートレインのうち最も信号強度が高いTEの位置のエコーを用いることが好ましい。このようなエコーは、後述する位相累積処理を開始する基準エコーとなる。<<閾値処理ステップ212>>2つの奇数エコー間で位相差分する際に、信号強度が低く位相値の信頼性が低いデータをフィッティングに用いないために、エコーデータから信号強度が低いデータを除外する。除外するデータは、所定の閾値以下、例えば対象とするエコーデータの最大振幅値に対して所定の比率例えば20以下のデータとする。閾値Thresholdは次式1で表わされる。数1ThresholdSignalmax×0.21ここで、Signalmaxはエコーデータ内の最大振幅値を意味する。<<周波数方向のフィッティングステップ213>>次に2つのエコーデータを位相差分したデータに対し、周波数方向のフィッティングを行う。フィッティングは、最小二乗法等の公知の手法を用いることができ、式2で表わされる一次関数と実際の各データ点における位相差との差が最小となるように、式3及び4により、式2の係数a、bを求める。ここで、Difoは奇数エコー間の位相差分をフィッティングした一次関数であり、Xiは周波数方向の座標、Yiは位相値で、nは周波数点数データ点の点数を意味する以下、同様。なお、フィッティングに先立ち、位相の主値回りを除去するための位相アンラップ処理を行う。具体的には、隣合うデータ点の位相差を順次検索し、一つのデータ点の位相と次のデータ点の位相との差が2πを超えたときに、次のデータ点の位相に±2π×Nを加算する処理を行う。位相アンラップ処理は、一般的な手法を用いてもよいし、後述する簡易アンラップ処理を行ってもよい。奇数エコー間の位相差分は、2つのエコーから求めることができるが、2つのエコーのみを用いるのではなく、複数の奇数エコー間で式2の係数を求めて、これを平均することが望ましい。またマルチスライス計測の場合には、スライス毎に奇数エコー間の位相差分を求めて、これを平均してもよいし、特許文献1に記載されるように、スライスの信号を合成した後、エコー間の位相差分を求めてもよい。<<位相累積ステップ217>>2つのエコーデータ間の位相差は、位相エンコード方向のオフセット誤差に対応し、位相エンコード方向のオフセット誤差はエコートレイン方向に蓄積されるものである。そこで、上記フィッティングステップで求めた位相差を、エコートレイン方向に位相累積する。具体的には、位相累積を開始する基準となるエコーをエコートレインの2×B1番目のエコーであるとすると、位相累積後の2×i1番目のエコーの位相変化Pは、次式5で表わされる。数3PPDifo×5以上の各ステップ211~213、217により、エコートレインを構成する各奇数エコーについて、ピークを揃えるための補正用データが得られたことになる。<偶数エコーデータの処理>偶数エコーデータ22については、周波数エンコード傾斜磁場の印加タイミングおよびサンプリング時間の誤差と、位相エンコード方向のオフセット誤差とを同時に補正するために、まず奇数エコーとの位相差分を行いステップ214、位相差分データについて周波数方向のフィッティングを行うステップ230。奇数エコーと偶数エコーとの位相差分も、理論的には一次直線でフィッティングできるものであるが、奇数エコー間或いは偶数エコー間の位相差分が小さい値の範囲に留まるのに対し、スケール値が変化する範囲が大きくなるため、適切な位相アンラップ処理を行えない可能性があり、単純な一次フィッティングを行った場合、フィッティング誤差が生じる可能性が高まる。そこで偶数エコーデータ奇数エコーとの位相差分データの処理ステップ230では、奇数エコー間の処理とは異なる位相アンラップ処理と限定されたフィッティングを行う。以下、ステップ230の処理の詳細を、図5を参照して説明する。図5は、ステップ230内の処理を示すフローである。<<閾値処理ステップ231>>まず、奇数エコーデータに対する閾値処理212と同様に、エコーデータの振幅値に対して閾値例えば最大振幅値の20を設定し、それ以下の周波数方向の点は、これ以降のステップから除外するステップ231。図6に、閾値処理前後のエコーデータを示す。図6a、cは閾値処理前のエコーデータの振幅及び位相、図6b、dは閾値処理後のエコーデータの振幅及び位相である。図示するように、最大振幅の20以下のデータは、位相がゼロであるとして扱う。図6dからもわかるように、エコー信号位相値は主値-π~πの間の範囲で回転するため、周波数方向の位相変化を一次の線形変化にするために、主値回りを取り除く処理位相アンラップ処理が必要となる。<<第1位相アンラップ処理ステップ232>>位相アンラップは、主値回りが生じる点で本来の位相変化よりも±2πの変化が生じることを利用して、隣り合う点同士で位相差が十分に小さくなるよう、±2π加算する処理を行う。しかし、本実施形態では、±2π加算するのではなく、位相差を単位ベクトルに変換することで、位相差を-π~πの間の主値内に納める簡易な位相アンラップ処理を行う。具体的には、隣り合うデータ点の位相θn、θn1の差分が-π~πの間に入らない場合には、この差分を、式6、7に示すように、複素数実部虚部に変換する。この実部と虚部のデータを用いて、式8にて位相値Δθを算出する。この位相値Δθを、一方のデータ点の位相θnに加算し式9、それに隣接するデータ点の位相θn1を求める。これにより、θnとθn1間の±2πの変化を除去できる。この位相アンラップ処理は、±2π×Nの位相をNを変えながら適当に加えて、目標位相との誤差が最も小さくなるNを探索する処理を必要とする技術特許文献2に記載された処理に比べ、簡便な演算処理といえる。ステップ232の位相アンラップ処理前後の位相データを図7に示す。図7aに示した主値内で回転する位相は、位相アンラップ処理によって、図7bに示すような本来の1次の線形な変化に戻すことができる。ただし、これら位相データは、ステップ231の閾値処理を行うことによって、処理から除外されたデータ点或いはデータ点の領域除外領域が存在するため、位相が連続しているデータ点の領域が除外領域によって隔てられる可能性がある。位相が連続しているとは位相飛びがないことを意味し、この位相が連続しているデータ点の領域を島領域と呼ぶ。これら島領域が離れて存在する場合、位相アンラップ処理の対象となる隣り合う点の距離が離れることになり、正しく位相アンラップ処理を行えない場合がある。例えば、図8に示す4つの島島1~島4で構成されたデータの場合、小さな島々で構成される左から2番目の島2の右端の点における位相値と、左から3番目の島3の左端の点における位相値の差は、-π~πの範囲の間に収まるため、上述した簡易アンラップ処理では、そのままの位相になるが、本来は-2π加算しなければならない。このように、本来隣接する点の位相データがないと連続性が確認できない。そこで、位相の連続性を持たせるための第2の位相アンラップ処理を行う。<<第2の位相アンラップ処理ステップ233~235>>第2の位相アンラップ処理を行うために、まず複数の島のうち、連続するデータ点の数が最も多い最大島の探索するステップ233。次いで、最大島の位相データに対して、最小二乗法を用いたフィッティング処理を行って、近似関数Difを算出し、その1次の傾きを求めるステップ234。最小二乗法を用いたフィッティングは、奇数エコーのフィッティングステップ214で説明した式2~4と同様の式を用いることができる。上述した最大島の検索と、最大島の位相データを用いたフィッティングの様子を図9に示す。図9aは、連続するデータ点の点数が最も多い島の位相データを取り出した状態を示し、bはaの位相データを用いて式3、4によりフィッティングして求めた一次関数点線を示す。次にステップ234で求めた近似関数の傾きを用いて、島と島との間の位相値に連続性を持たせるための位相アンラップ処理第2アンラップ処理を行うステップ235。具体的には、位相データがないゼロである部分によって隔てられた2つの島の隣り合う点となる端部の位相値の一方を、他方の位相値と近似関数の傾きを用いて、位相アンラップ処理する。位相アンラップ処理は、周波数方向に配列するデータ点の左側から開始しても右側開始してもよいが、以下、左側から開始する場合を例に説明する。まず2つの島のうち左側の島の右端の位相値と上記近似関数の傾きを用いて、右側の島の左端の位相値の目標値を算出する。目標値とはステップ234で求めた1次の位相変化に従えば、隣接する島の左端の点は本来この位相値になるであろう値を意味する。この目標値と実際の位相は大きくずれることはないため、この位相差が-π~πに収まるよう、ステップ232の位相アンラプ処理と同様の演算処理を行う。即ち、近似関数から推定した位相目標値θtargetと、島の左端のデータ点の位相値θnとの差分を以下の式10、11で示すように、複素数に変換し、その実部と虚部のデータを用いて式12により位相値Δθを算出する。この位相値Δθを位相目標値に加算し式13、そのデータ点の位相値とする。右側の島が、位相値が連続する複数のデータ点からなる場合には、島の各データ点全てに、式12で算出した位相値Δθを、第1位相アンラップ処理後の位相値に加算する。第2アンラップ処理の様子を図10に示す。図示する例では、検索を左から行い、小さい島からなる島2と、その右側の島3との間でアンラップ処理する場合a及び島3と島4との間でアンラップ処理する場合bを示している。島2と島3との間には、比較的大きいデータ欠落部分があり、最初の位相アンラップ処理第1アンラップ処理ステップ232では、この間のアンラップ処理が適切に行われていない。そこで島3の左端のデータ点について、ステップで求めたフィッティング関数の傾きから、この点の目標値を求め、実際の位相値との差分を求め、上記式10~13から算出したΔθを島3の左端の実際の位相値に加算する。島3の全てのデータ点図では2点の位相値にも同じΔθを加算する。島3と島4との間でも、同様にして、島3の右端のデータ点の位相値とフィッティング関数の係数から、島4の左端のデータ点の位相値の目標値を算出し、目標値と実際の位相値との差が-π~πの範囲の間に収まるように、Δθを求め加算する図10b。図10cは、この第2アンラップ処理後の状態を示す図であり、閾値処理によって除去されたデータ点以外の全てのデータ点を通して、位相飛びのない連続性のある位相値になっている。<<全データに対するフィッティングステップ236>>第2アンラップ処理を終えた後、除去されたデータ点以外の全てのエコーデータの位相値を用いて、再度最小二乗法を用いてフィッティング処理を行うステップ236。フィッティング処理は、奇数エコーのフィッティングで用いた式と同様の式14~16を用いて、フィッティング関数を求める。式中、DifOEは奇数エコー基準エコーと偶数エコー基準エコーとの位相差分を表し、cはその傾き、dは切片を表す。このフィッティング関数から算出される各データ点の位相値フィッティングデータを、そのまま、補正用に用いることも可能であるが、本実施形態では、フィッティングデータを信号強度に応じて制限して使用し、フィッティングに伴うフィッティング誤差を極力小さくする。即ち、フィッティングには必ずフィッティング誤差が生じ、この誤差は直接N/2アーチファクトの発生に繋がる。そのため、位相変化に乱れが生じていない場合には、今回の抑制ターゲットとなるストリークアーチファクトも発生しないため、フィッティングデータを敢えて利用する必要はない。一方、信号値が弱い領域では、フィッティング誤差によって生じるN/2アーチファクトの信号強度も相対的に低い。このような理由から、本実施形態では、信号値が低く位相が乱れやすい領域点に対してのみ、フィッティングデータを限定的に用い置換する。<<フィッティングデータを用いたデータ置換ステップ237~239>>フィッティングデータの限定的な使用は、マルチスライス計測の場合に、エコー全体の信号が弱いスライスのみを対象とし、十分な信号値が得られたスライス内ではエコーデータのうち信号値の弱い領域のみを対象とする。まず、スライス毎にエコー全体の信号が弱いか否かの判定を行うステップ237。このため、マルチスライス間で各最大信号値の比較を行い、スライス間で閾値を設定する。閾値は、例えばマルチスライスの最大信号値に対して20とし、最大信号値が閾値以下のスライスは、エコー全体の信号が弱いと判定し、そのスライスのエコーデータ全体をフィッティングデータに置換するステップ239。十分に信号値が得られているスライスに対しては、既に閾値処理ステップ231によって信号値が弱い領域は除外されているので、この除外した領域に対してデータの置換を行うステップ238。この置換処理は、全島の両端に対しては、ステップ ゚236のフィッティング処理で求めた1次の傾きcを用いて、外挿処理を行う。島と島との間については、既に第2位相アンラップ処理ステップ235において位相の連続性が保たれているため、線形補間処理を行えばよい。図11は、フィッティング処理の様子を示す図であり、図11aは処理前の位相データ第2アンラップ処理後の位相データ:図10cを示し、図11bは、十分に信号値が得られているスライスに対する処理、図11cは、エコー全体の信号値が低いスライスに対する処理を示す。図11bで点線で囲った部分がフィッティングを行った部分であり、両端はフィッティングデータを用いて外挿し、中間部分は線形補間している。また図11cに示すスライスでは、位相データ全体がフィッティングデータにより置換されている。このようなフィッティングを行うことにより、フィッティング誤差を最小に留め且つ信号値の弱い領域についても誤差の少ない位相値を得ることができる。以上、図5を参照して、図2のステップ230の詳細を説明したが、このステップ230は、複数の奇数エコーと偶数エコーの組例えば第1エコーと第2エコー、及び、第3エコーと第4エコーについて行ってもよいし、一組について行ってもよいが、少なくとも、次の位相累積処理ステップ217で基準となるエコーに対して行うことが好ましい。複数組について位相差分のフィッティングデータを求めた場合、その平均値を基準エコー間の位相差分データとしてもよい。再び図2を参照して、ステップ230以降の処理を説明する。<<偶数エコー間の位相差分ステップ215>>偶数エコーについても、エコートレイン方向のピークずれを補正するために偶数エコー間の位相差分を行う必要がある。偶数エコー間の位相差分の処理は、奇数エコーと同様に、2つの偶数エコー、例えばエコートレインの2番目のエコーと4番目のエコーについて、エコーデータx-kyデータの位相差分を取り、周波数方向にフィッティングすることにより求めることもできるが、本実施形態では、ステップ214或いはステップ230で求めた奇数エコーと偶数エコーとの位相差分について、エコートレイン方向の位相変化差分を求める。これにより、N2アーチファクトの原因である奇数エコーと偶数エコーのピークずれを併せて補正できるデータを得ることができる。具体的には、ステップ230で求めた奇数エコーの基準エコー奇数基準エコーと偶数エコーの基準エコー偶数基準エコーとの位相差分に対する、他の隣接する奇数エコーと偶数エコーとの位相差分との差分を取ることになる。この差分データに対して、閾値処理すること、最小二乗法を用いたフィッティングを行うステップ216。この処理は、奇数エコーに対し行った処理ステップ212、213と同じであり、説明を省略する。ステップ216では、次式17で示す一次関数が得られる。数7DifoeDifOEexg17式中、DifOEは、ステップ230で最終的に得られた、偶数基準エコーと奇数基準エコーとの位相差分を表し、eは一次関数の傾き、gは切片を示す。<<エコートレイン方向位相累積>>位相累積は、奇数エコーの処理と同様に、基準エコーを開始点として、偶数エコー間の位相差分を順次エコートレイン方向に累積する。奇数エコーから偶数エコーへの位相差分について位相累積を行った場合、位相累積された位相データには、奇数エコーと偶数エコーとのピークずれの情報と基準エコーに対するピークずれの情報とを含んでいる。そこでステップ218において、位相差分データに対し、位相累積後の奇数エコーデータの位相を加算し、偶数エコーの補正用データを得る。具体的には、式17を用いて、次式18の計算を行う。数8P’P’Difoe18式中P’は位相累積後の前計測データの位相変化である。以上の処理ステップ211~218は、マルチスライスの場合には、スライス毎に行い、スライス毎に奇数エコー及び偶数エコーの補正用データを得る。<本計測データの補正ステップ222、223>周波数方向にフーリエ変換した本計測データx-kyデータの位相に対し、上記ステップ217及び218で得た、前計測の奇数エコーの位相と偶数エコーの位相を加算することにより、x-ky空間上で本計測データの位相が補正されるステップ222。x-ky空間上での位相補正は、k空間上でのピークずれの補正を意味し、この補正によってk空間上でのピークずれが解消される。マルチスライス計測の場合には、補正用データはスライス毎に用意されるので、本計測の各スライスに対し、対応するスライスの補正用データを用いて位相補正する。位相補正後のx-kyデータを位相方向にフーリエ変換することによりステップ223、最終的に画像31が再構成される。再構成画像では、ピークずれに起因するN2アーチファクトやストリークアーチファクトが大幅に低減される。図12は、前後の腹部DW-EPI画像を示す図であり、aは、奇数エコー及び偶数エコーについて、それぞれ位相データに対する単純な一次フィッティングを行って位相補正した場合の画像、bは、図5に示す処理ステップ230を行って位相補正した場合の画像である。図12aではストリークアーチファクトが発生していたが、bでは発生していないことが確認できる。以上、本発明の一実施形態を説明し、その実施形態では、エコートレイン方向の位相累積を行う際の基準となるエコーデータの、奇数偶数基準エコー間位相差分について、限定的なフィッティング図4のステップ230を行う場合を説明したが、奇数エコー間位相差分に対する周波数方向フィッティングステップ213や偶数エコー間位相差分に対する周波数方向フィッティングステップ216についても、同様の限定的なフィッティングを採用することも可能である。また上記実施形態では、奇数エコーについて奇数エコー間位相差分を求め、偶数エコーについては奇数エコーとの差分に対しエコー間位相差分を行う場合を説明したが、逆に奇数エコーと偶数エコーの処理を逆にしても同様の結果を得ることができる。さらに本発明の簡易位相アンラップ処理を含む周波数フィッティングの技術は、上述したEPI計測におけるエコーデータ間の位相差分に対してのみならず、他のエコー信号の位相補正や画像データの位相補正にも拡張して適用することが可能である。例えば、水脂肪分離計測におけるピークずれの補正、realIRシーケンスやラジアルシーケンスの補正、および、エコーピーク検出による傾斜磁場のディレイ測定にも適用することが可能である。"