"被検体の撮像部位の実空間における動きの大きさの空間分布を示す動き情報を取得する動き情報取得部と、前記動き情報に基づいて、磁気共鳴イメージングのスキャンにより収集された前記被検体の撮像部位の画像データの第1の領域において前記動きに因る位置シフトおよび位相シフトのいずれかを第2の領域と異なる強度によるk空間上における線形の補正処理によって補正する補正部と、前記補正部により補正された前記第1の領域および前記第2の領域の各画像データをk空間上で合成する合成部と、前記k空間上で合成された画像データを実空間上の画像に再構成する再構成部と、を備える画像データ補正装置。前記補正部は、前記撮像部位の動きに伴って生じる前記画像データの空間的に不均一な劣化の程度に応じた線形の補正処理を行うように構成される請求項1に記載の画像データ補正装置。前記合成部は、前記第1の領域における第1の補正後における画像データと前記第2の領域における前記第1の補正と異なる第2の補正後における画像データとを合成するように構成される請求項1に記載の画像データ補正装置。前記合成部は、前記第1の領域における補正後における画像データと前記第2の領域における非補正の画像データとを合成するように構成される請求項1に記載の画像データ補正装置。前記動き情報取得部は、前記撮像部位からイメージング用のデータ収集を行うときにナビゲーター用のエコー信号を収集するナビゲーターエコー収集部と、前記エコー信号を処理して前記動き情報を生成する動き情報生成部とを備える請求項1に記載の画像データ補正装置。前記動き情報取得部は、前記撮像部位の動きを外部から光学的に又は空気圧で検知するセンサと、前記センサが検知した信号を処理して前記動き情報を生成する動き情報生成部と、を備える請求項1に記載の画像データ補正装置。前記補正部は、前記画像データの3つ以上の領域についてそれぞれ異なる補正を行うように構成される請求項1に記載の画像データ補正装置。前記動き情報取得部は、前記撮像部位からイメージング用のデータ収集を行うときにナビゲーター用のエコー信号を収集するナビゲーターエコー収集部と、前記エコー信号から各部位の動きの振幅および位相の平均値を含む処理情報を得る処理情報取得部と、前記空間分布の全体のプロファイルを予め測定又は仮定したモデルを用いることにより取得するプロファイル取得部と、前記プロファイルと前記処理情報とに基づき前記動き情報を生成する動き情報生成部と、を備える請求項1に記載の画像データ補正装置。前記動き情報取得部は、前記撮像部位からイメージング用のデータ収集を行うときに当該イメージングに伴うデータ読出し方向および位相エンコード方向のいずれかにおけるナビゲーター用のエコー信号を収集するナビゲーターエコー収集部と、前記エコー信号を処理してプロジェクションデータを生成するプロジェクションデータ生成部と、前記プロジェクションデータを前記動き情報として用いる動き情報生成部と、を備える請求項1に記載の画像データ補正装置。前記動き情報取得部は、前記撮像部位からイメージング用のデータ収集を行うときにナビゲーター用のエコー信号を収集するナビゲーターエコー収集部と、前記エコー信号のk空間上の位相分布およびr空間上の少なくとも一方向の位置のいずれかのシフトを前記動き情報として求める動き情報生成部と、を備える請求項1に記載の画像データ補正装置。前記画像データを単一の受信RFコイルを用いて収集する画像データ収集部を備える請求項1に記載の画像データ補正装置。前記補正部は、前記撮像部位における前記画像データに互いに異なる重み分布を有する複数のウィンドウ関数を乗じることにより、複数の画像データを生成するように構成される請求項1に記載の画像データ補正装置。前記補正部は、位相および位置のいずれかの平均シフト処理、位相および位置のいずれかのアフィン変換、または前記平均シフト処理と前記アフィン組み合わせで表される線形の補正処理を行うように構成される請求項1記載の画像データ補正装置。前記スキャンを、マルチショットタイプおよびシングルタイプのいずれかのパルス列であって、スピンワープ法、スパイラル法およびラジアル法の何れかに基づくパルス列で構成されたパルスシーケンスを用いて実行する画像データ収集部を備える請求項1記載の画像データ補正装置。前記動き情報取得部は、3次元の空間分布を示す動き情報を取得するように構成される一方、前記補正部は、前記3次元の空間分布を示す動き情報に基づいて、3次元的に補正を行うように構成される請求項1記載の画像データ補正装置。前記情報取得部は、非剛体の動きの大きさの空間分布を示す動き情報を取得するように構成される請求項1記載の画像データ補正装置。前記動き情報取得部は、前記撮像部位における予め仮定された動きの大きさの空間分布を用いるように構成される請求項1に記載の画像データ補正装置。前記補正部は、前記撮像部位の前記画像データにおける少なくとも3つ以上の領域の一部を互いに合成して得られる2つ以上の領域に対して非補正を含む互いに異なる補正を行うように構成される請求項1記載の画像データ補正装置。前記画像データを、予め前記撮像部位の動きの大きさに応じた感度分布を持って収集する画像データ収集部を備える請求項1に記載の画像データ補正装置。前記動き情報取得部は、腹部を前記撮像部位とし、前記腹部の背中側から腹壁側へ向かう方向において各部位の動きの振幅が前記背中側から前記腹壁側に向かって線形に増加する分布を前記空間分布として用いるように構成される請求項17に記載の画像データ補正装置。前記補正部は、複数のエレメントコイルを用いて収集した複数の画像データの一部を互いに合成することにより前記2つ以上の領域を得るように構成される請求項18に記載の画像データ補正装置。前記補正部は、前記撮像部位における前記画像データに互いに異なる重み分布を有する複数のウィンドウ関数を乗じることにより、合成対象となる前記3つ以上の画像データを生成するように構成される請求項18に記載の画像データ補正装置。前記画像データ収集部は、互いに異なる感度分布の複数のエレメントコイルを有するマルチコイルを用いて前記画像データを収集するように構成される請求項19に記載の画像データ補正装置。前記画像データ収集部は、腹部を前記撮像部位とし、前記複数のエレメントコイルとして背中側及び腹壁側にそれぞれ配置された2つの表面コイルを備える請求項23に記載の画像データ補正装置。前記画像データ収集部は、腹部を前記撮像部位とし、前記複数のエレメントコイルとして背中側及び腹壁側にそれぞれ配置され、且つ前記背中側及び前記腹壁側の少なくとも一方にて頭部から脚部に向かう方向に沿って複数個所に配置された複数の表面コイルを備える請求項23に記載の画像データ補正装置。前記補正部は、前記複数のエレメントコイルを用いて収集した複数の画像データの一部または全部に、互いに異なる重み分布を有する複数のウィンドウ関数を乗じることにより、複数の画像データを生成するように構成される請求項23に記載の画像データ補正装置。被検体の撮像部位の実空間における動きの大きさの空間分布を示す動き情報に基づいて、磁気共鳴イメージングのスキャンにより収集された前記撮像部位の画像データの第1の領域において前記動きに因る位置シフトおよび位相シフトのいずれかを第2の領域と異なる強度によるk空間上における線形の補正処理によって補正する補正部と、前記補正部により補正された前記第1の領域および前記第2の領域の各画像データをk空間上で合成する合成部と、前記k空間上で合成された画像データを実空間上の画像に再構成する再構成部と、を備える画像データ補正装置。前記補正部は、前記撮像部位の動きに伴って生じる前記画像データの空間的に不均一な劣化の程度に応じた線形の補正処理を行うように構成される請求項27に記載の画像データ補正装置。実空間における被検体の撮像部位の動きの大きさの空間分布を示す動き情報を取得する動き情報取得部と、前記被検体の撮像部位の第1の領域における実空間上のデータを、実空間とk空間のハイブリッド空間である第3の領域のデータに変換すると共に、前記被検体の撮像部位の第2の領域における実空間上のデータを、実空間とk空間のハイブリッド空間である第4の領域のデータに変換する第1のデータ変換部と、前記第3の領域における前記k空間上のデータの動きに因る位置シフトおよび位相シフトのいずれかに対し、前記動き情報に基づいて、前記第4の領域と異なる強度で線形の補正処理を施す補正部と、補正後における前記第3の領域における前記ハイブリッド空間のデータと、前記第4の領域における前記ハイブリッド空間のデータとを合成する合成部と、合成後における前記ハイブリッド空間のデータを前記実空間上のデータに変換する第2の変換部と、を備える画像データ補正装置。前記合成部は、補正後における前記第3の領域における前記ハイブリッド空間のデータおよび前記補正と異なる補正後における前記第4の領域における前記ハイブリッド空間のデータを合成するように構成される請求項29に記載の画像データ補正装置。前記合成部は、補正後における前記第3の領域における前記ハイブリッド空間のデータおよび非補正の前記第4の領域における前記ハイブリッド空間のデータを合成するように構成される請求項29に記載の画像データ補正装置。前記補正部は、前記撮像部位の動きに伴って生じる前記実空間上のデータの空間的に不均一な劣化の程度に応じた線形の補正処理を行うように構成される請求項29に記載の画像データ補正装置。被検体の撮像部位の実空間における動きの大きさの空間分布を示す動き情報を取得する動き情報取得部と、磁気共鳴イメージングのスキャンにより前記被検体の撮像部位の画像データを収集する画像データ収集部と、前記動き情報に基づいて、収集された前記画像データの第1の領域において前記動きに因る位置シフトおよび位相シフトのいずれかを第2の領域と異なる強度によるk空間上における線形の補正処理によって補正する補正部と、前記補正部により補正された前記第1の領域および前記第2の領域の各画像データをk空間上で合成する合成部と、前記k空間上で合成された画像データを実空間上の画像に再構成する再構成部と、を備える磁気共鳴イメージング装置。磁気共鳴イメージングのスキャンにより被検体の撮像部位の画像データを収集する画像データ収集部と、前記撮像部位の実空間における動きの大きさの空間分布を示す動き情報に基づいて、収集された前記画像データの第1の領域において前記動きに因る位置シフトおよび位相シフトのいずれかを第2の領域と異なる強度によるk空間上における線形の補正処理によって補正する補正部と、前記補正部により補正された前記第1の領域および前記第2の領域の各画像データをk空間上で合成する合成部と、前記k空間上で合成された画像データを実空間上の画像に再構成する再構成部と、を備える画像データ補正装置。磁気共鳴イメージングのスキャンにより被検体の撮像部位における実空間上のデータおよび前記被検体の撮像部位の実空間における動きの大きさの空間分布を示す動き情報を収集するデータ収集部と、第1の領域における前記実空間上のデータを、実空間とk空間のハイブリッド空間である第3の領域のデータに変換すると共に、第2の領域における前記実空間上のデータを、実空間とk空間のハイブリッド空間である第4の領域のデータに変換する第1のデータ変換部と、前記第3の領域における前記k空間上のデータの動きに因る位置シフトおよび位相シフトのいずれかに対し、前記動き情報に基づいて、前記第4の領域と異なる強度で線形の補正処理を施す補正部と、補正後における前記第3の領域における前記ハイブリッド空間のデータと、前記第4の領域における前記ハイブリッド空間のデータとを合成する合成部と、合成後における前記ハイブリッド空間のデータを前記実空間上のデータに変換する第2の変換部と、を備える画像データ補正装置。本発明は、上述した従来技術が抱える状況に鑑みてなされたもので、磁気共鳴イメージングにおいて、特にマルチショット型のパルスシーケンスを用いて収集充填したk空間のデータから画像を得る場合に、撮像部位の非線形の動きに因って空間的に不均一に劣化して収集されたエコーデータを比較的簡便にかつ高速に補正してMR画像を得ることが可能な磁気共鳴イメージング装置および画像データ補正装置を提供することを目的とする。画像データ補正装置は、上述した目的を達成するために、画像データ補正装置は、被検体の撮像部位の実空間における動きの大きさの空間分布を示す動き情報を取得する動き情報取得部と、前記動き情報に基づいて、磁気共鳴イメージングのスキャンにより収集された前記被検体の撮像部位の画像データの第1の領域において前記動きに因る位置シフトおよび位相シフトのいずれかを第2の領域と異なる強度によるk空間上における線形の補正処理によって補正する補正部と、前記補正部により補正された前記第1の領域および前記第2の領域の各画像データをk空間上で合成する合成部と、前記k空間上で合成された画像データを実空間上の画像に再構成する再構成部と、を備える。このような本発明に係る磁気共鳴イメージング装置および画像データ補正装置によれば、磁気共鳴イメージングにおいて、特にマルチショット型のパルスシーケンスを用いて収集充填したk空間のデータから画像を得る場合に、撮像部位の非線形の動きに因って空間的に不均一に劣化して収集されたエコーデータを、動きの程度に応じて線形補正したデータの和合成として補正することで、比較的簡便にかつ高速に補正されたMR画像が得られる。以下、図面を参照して本発明の実施の形態を説明する。第1の実施の形態以下、本発明の磁気共鳴イメージング装置に係る一実施形態を図1~5を参照して説明する。この実施形態にかかる磁気共鳴イメージングMRIMagneticResonanceImaging装置の概略構成を図1に示す。この磁気共鳴イメージング装置は、被検体Pを載せる寝台部と、静磁場を発生させる静磁場発生部と、静磁場に位置情報を付加するための傾斜磁場発生部と、高周波信号を送受信する送受信部と、システム全体のコントロール及び画像再構成を担う制御演算部とを備えている。静磁場発生部は、例えば超電導方式の磁石1と、この磁石1に電流を供給する静磁場電源2とを備え、被検体Pが遊挿される円筒状の開口部診断用空間の軸方向Z軸方向に静磁場H0を発生させる。なお、この磁石部にはシムコイル図示せずが設けられている。寝台部は、被検体Pを載せた天板Tを磁石1の開口部に退避可能に挿入できる。傾斜磁場発生部は、磁石1に組み込まれた傾斜磁場コイルユニット3を備える。この傾斜磁場コイルユニット3は、互いに直交するX軸方向、Y軸方向及びZ軸方向の傾斜磁場を発生させるための3組種類のx,y,zコイル3x~3zを備える。傾斜磁場部はまた、x,y,zコイル3x~3zに電流を供給する傾斜磁場電源4を備える。この傾斜磁場電源4は、後述するシーケンサ5の制御のもとで、x,y,zコイル3x~3zに傾斜磁場を発生させるためのパルス電流を供給する。傾斜磁場電源4からx,y,zコイル3x~3zに供給されるパルス電流を制御することにより、物理軸である3軸X軸,Y軸,Z軸方向の傾斜磁場を合成して、互いに直交するスライス方向傾斜磁場GS、位相エンコード方向傾斜磁場GE、および読出し方向周波数エンコード方向傾斜磁場GRから成る論理軸方向を任意に設定変更することができる。スライス方向、位相エンコード方向、および読出し方向の各傾斜磁場は静磁場H0に重畳される。送受信部は、磁石1内の撮影空間にて被検体Pの近傍に配設される送信RF高周波コイル7T及び受信RFコイル7Rと、このRFコイル7T,7Rにそれぞれ接続された送信器8T及び受信器8Rとを備える。受信RFコイル7Rは、コイルエレメントを1つであるシングルコイル又はコイルエレメントが複数のマルチコイルであり、何れも表面コイルとして形成されている。受信RFコイルは被検体P被験者の撮像部位としての腹部などの体表面に沿って配置される。なお、送信RFコイルは例えば全身用コイルとして形成されている。この送信器8T及び受信器8Rは、後述するシーケンサ5の制御のもとで動作する。この動作により、送信器8Tは、核磁気共鳴NMRnuclearmagneticresonanceを励起させるためのラーモア周波数のRF電流パルスを送信RFコイル7Tに供給する。受信器8Rは、受信RFコイル7Rが受信した磁気共鳴MRmagneticresonance信号高周波信号を取り込み、これに前置増幅、中間周波変換、位相検波、低周波増幅、フィルタリングなどの各種の信号処理を施した後、A/D変換してMR信号のデジタルデータ原データを生成する。さらに、制御演算部は、シーケンサシーケンスコントローラとも呼ばれる5、ホスト計算機6、演算ユニット10、記憶ユニット11、表示器12、及び入力器13を備える。この内、ホスト計算機6は、記憶したソフトウエア手順図示せずにより、シーケンサ5にパルスシーケンス情報を指令するとともに、装置全体の動作を統括する機能を有する。スキャンに用いるパルスシーケンスとしては、高速SE法、高速FE法、EPI法、FASEspinecho法、GRASE法などに拠る2次元又は3次元のマルチショットタイプのパルス列が使用される。補正をk空間で実行するとしても、r空間の情報に基づいてデータを分割するために、k空間を如何なる方法でサンプリングしても良い。このため、かかるパルス列は、具体的には、スピンワープspinwarp法、スパイラルspiral法、又は、ラジアルradial法の何れに基づくパルス列であってもよい。これにより、後述する補正によって、アーチファクトが発生源付近に局所化する場合に、特に大きな補正効果が発揮される。なお、発生するエコー間の動きが無視できないような磁気共鳴イメージングの場合には、シングルショットタイプのスキャンであっても、本発明に係る補正を適用することができる。シーケンサ5は、CPUおよびメモリを備えており、ホスト計算機6から送られてきたパルスシーケンス情報を記憶し、この情報にしたがって傾斜磁場電源4、送信器8T、受信器8Rの動作を制御するとともに、受信器8Rが出力した磁気共鳴信号のデジタルデータを一旦入力し、これを演算ユニット10に転送するように構成されている。ここで、パルスシーケンス情報とは、一連のパルスシーケンスにしたがって傾斜磁場電源4、送信器8Tおよび受信器8Rを動作させるために必要な全ての情報であり、例えばx,y,zコイル3x~3zに印加するパルス電流の強度、印加時間、印加タイミングなどに関する情報を含む。また、演算ユニット10は、受信器8Rが出力したデジタルデータ原データ又は生データとも呼ばれるを、シーケンサ5を通して入力し、その内部メモリによるk空間フーリエ空間または周波数空間とも呼ばれるにそのデジタルデータを配置し、このデータを1組毎に2次元または3次元のフーリエ変換に付して実空間の画像データに再構成する。また演算ユニット10は、本発明に係る補正処理を所定のアルゴリズムにしたがって実行できるようになっている。この補正処理は、現時点で再構成した或いは既に取得している画像データの、被検体の体動に因る空間的に不均一非線形な劣化を後処理として補正する処理である。被検体Pの例えば腹部を撮像する場合、腹部は呼吸性の動き体動に因って、その各サンプリング位置が空間的に非線形に動く。これにより、取得した画像データは空間的に不均一に劣化しているのである。このことは、本願発明の中枢を成す事項であるので、後で詳述する。さらに演算ユニット10は、必要に応じて、画像に関するデータの合成処理や差分演算処理も実行可能になっている。この合成処理には、画素毎に加算する処理、最大値投影MIPMaximumIntensityProjection処理などが含まれる。記憶ユニット11は、再構成された画像データのみならず、上述の合成処理や差分処理が施された画像データを保管することができる。表示器12は例えば再構成画像を表示するのに使用される。また入力器13を介して、術者が希望するパラメータ情報、スキャン条件、パルスシーケンス、画像合成や差分の演算に関する情報などをホスト計算機6に入力できる。ここで、上述した画像データに施される後処理としての補正処理について説明する。非線形な動きの性質被検体である生体の空間的な動きには、線形な動きと非線形の動きとがあることは既述した。このうち、空間的に非線形な動きとは、その動きの振幅や位相に空間的な分布があることを意味している。周期性に関して言えば、呼吸性や心拍などに因る動きは周期性である場合が多いが、不意の突発的な動きinvoluntarymotionは非周期的である。腹部臓器を画像化する場合、腹部の撮像領域の各サンプリング位置は呼吸性の動き体動の影響を強く受ける。この呼吸性の動きは、空間的に非線形不均一になる。そこで、本実施形態は、腹部臓器などを診断するときに、かかる呼吸性の非線形の動きに因って空間的に不均一に分布した動き成分が混じったエコー信号を収集し、このエコー信号の空間的に不均一に劣化した信号成分を補正する。尚、動きは、平行移動と回転のみを伴う剛体変形と、さらに拡大縮小やせん断を含めた線形変形や非線形変形を伴う非剛体変形とに分類することもできる。本発明によれば、より正確には非剛体変形の動きに因って空間的に不均一に分布する動き成分が混在するエコー信号に対して効果的に補正を行うことができる。従って、非剛体変形のうち線形変形により劣化したデータを対象として補正を行うことも可能であるが、以下では非線形変形に対する補正を対象として説明する。この補正を行うために、生体の動きに関してあるモデルを仮定する図2参照。このモデルは、空間的に非線形な動きを最も良く示す断面像として、腹部のアキシャル断面像を対象とするが、必ずしもアキシャル断面像でなくてもよい。腹部のアキシャル断面像の場合、図2に示すモデルのように、呼吸性の動きの振幅は背中側では無視できるが、その反対の腹壁側前側が最も大きくなり、その間の振幅は背中側から腹壁側へ漸増していると仮定できる。いま、ある2次元対象物の動きの振幅が、x方向では一様であり、y方向では背中側から腹壁側へその距離yに比例して大きくなっているとすれば、腹部のアキシャル断面形状がほぼ楕円であること及びコントラストもほぼ前後左右対称であることを考慮すると、あるタイミングで取得した画像全体の平均シフト量は被写体中心部の動きを反映していると考えられる。このとき、y方向での動きの振幅は、背中側端ではゼロとなり、腹壁側端では最大かつ中心部の2倍の変化を示すものと見做すことができる。そこで、本実施形態では、r空間rspaceで計測又は予測した非線形な動きの情報、すなわち、その動きの大きさオーダーを用いてk空間kspaceを分割し、その複数の分割k空間で互いに異なる0次の位相補正を行い、その補正されたデータを相互に重み付け加算するという手法を採用する。非線形な動きの検出とその補正次いで、上述した非線形な動きの検出及びその補正の具体的な手法を説明する。図3~5に、かかる手法の概要を示す。いま、撮像部位は被検体の腹部であり、その腹部の周囲にRF受信コイル7Rとしてシングルコイルから成る表面コイルが配置されているものとする図5a参照。前述したように、腹部の各位置は呼吸に伴って空間的に非線形に動くので、腹部を診断する場合、動きの大小に応じて空間的に分割して処理をする必要がある。一方で、処理の簡便さからk空間で位相補正を実施したいので、k空間に充填されたデータを空間的な動きの大きさに依存して分割したい。k空間に充填されたデータには、動きの分布に関する位置情報が既に均等に散逸してしまっているので、k空間の空間情報自体を用いた分割は困難である。そこで、動きの空間分布情報を用いてr空間のデータを分割する方法をとる。ここでは腹部のアキシャル断面を対象として、分布の指標として、a0次位相のモニタと動きの空間依存性は、腹側が背側より大きいという知識推定を用いるか、又は、b実際に被検体の前後方向背中側から腹壁方向に向かう方向を位相エンコード方向として、その方向の腹部の動きをナビゲーション用の傾斜磁場パルスに応答するナビゲーターエコー以下、ナビエコーと呼ぶなどでモニタする方法を採用する。なお、r空間における平均0次位相のシフトはk空間では1次の位相シフトになる。また、r空間における位置のシフトはk空間では1次の位相シフトになる。したがって、位相誤差がexpiΦ0の場合、その位相を補正するにはr空間及びk空間共に一様に各画素データとexp-iΦ0との積を演算すればよい。また、r空間でのシフト距離がy方向にΔYであるならば、k空間でexp[-2πKyΔY]との積を演算すればよい。上記a項に記載の、動きの平均振幅だけをモニタして位相エンコード方向に動きの振幅変化を線形に外装し補正する方法をとる場合、ナビエコーで位相エンコード方向に一定値の位相シフトをもたらす傾斜磁場を加えたパルスシーケンスを用いればよい。このパルスシーケンスの例を図3に示す。また、上述したb項に記載の位相エンコード方向のプロジェクションデータを収集する場合、ナビエコーの収集前に位相エンコード方向と読出し方向との両方に同様の傾斜磁場を印加すれよい。この例を図4に示す。なお、図3,4にはパルスシーケンスとしてスピンワープ法を適用した高速SE法に基づくパルス列を示すが、紙面の関係上、1ショットに対する1エコーのみを図示している。このような図3又は図4に示すパルスシーケンスを用いてエコー信号が収集され、そのエコーデータが2次元k空間上に配置される図5b参照。このエコーデータは、演算ユニット10により、2次元フーリエ変換2DFTtwodimensionalFouriertransformされた実空間の画像データに再構成される図5c参照。ここでは、腹部の空間的動きに応じてk空間のデータ分布を分割するに際して、最も単純に、k空間をその位相エンコード方向に2分割する例を示す。この場合、実空間の画像に位相エンコード方向に相当するy軸方向にウィンドウを掛け詳しくはウィンドウ関数を乗じる、動きの程度が腹壁側腹部の前側と背中側腹部の後ろ側とで各々支配的な2領域にデータ上で分割する図5c、d参照。このデータ上で分割された2つの実空間を2次元の逆フーリエ変換2DIFTtwodimensionalinverseFouriertransformにより、それぞれ2つの2次元k空間に変換する図5ef参照。一方、前記パルスシーケンスの実行によってショットn毎にナビエコーSnaviKy,nが収集されているので図5g参照、このナビエコーSnaviKy,nから腹部の動き情報としての位相誤差ΔΦ0Ky,nを測定する図5h参照。この測定は具体的には、[数1]Φ0Ky,n=2πYnKyΦ0Ky,nbase=2πYnbaseKyΔΦ0Ky,n=Φ0Ky,n-Φ0Ky,nbaseとしてショットn毎に実行される。但しYnは、ショットnにおける動きのシフト量である。次いで、腹壁側のデータを有するk空間のデータにのみ位相補正を行い図5i参照、背中側のデータを有するk空間のデータには位相補正を行わない。かかる位相補正は、ショットn毎に、各ショットのデータに[数2]exp[-iΔΦ0Ky,n]を乗じることで行なわれる。この2つのk空間を夫々2次元フーリエ変換2DFTして個々の実空間データに戻し、相互に画素値を加算して最終的な実空間データを得る図5j、k参照。ここで、上述したエコーデータナビエコーのデータを含むの収集及び補正処理を数式化して説明する。この処理は、本実施形態では、演算ユニット10により実行されるようになっているが、その処理の一部をホスト計算機6に任せるようにしてもよい。ここで、[数3]Simgx,y:イメージングする対象のエコーデータ、Snavix、Ky:1次元フーリエ変換FTされた、x軸読出し方向のナビエコー、Axnavi:n番目のショットでサンプリングされたx軸方向の動きの振幅、Aynavi:n番目のショットでサンプリングされたy軸方向の動きの振幅、ΦxnaviKx,n:n番目のショットでサンプリングされたx軸方向のナビエコーの位相、ΦynaviKy,n:n番目のショットでサンプリングされたy軸方向のナビエコーの位相、WL,M,Sy:人体における動きの支配領域を選択するためのウィンドウ関数、n:ショット番号1=1~Nmaxとする。最初に、前述した図3又は図4に示すマルチショットタイプの最も単純な1エコー/ショットである、通常のスピンワープ法によるFSE法に基づくパルスシーケンスを用いてエコーデータを収集して、その画像データが得られる処理1。つまり、エコーデータSimgKx,Kyとすると、その実空間の画像データPimgx,yは、[数4]Pimgx,y=FT2DKx→x、Ky→y[SimgKx,Ky]として得られる。このとき、ナビエコーSnaviKy,nも同時に収集される。次いで、演算ユニット10は、動きのオーダ大きさに応じて空間領域を選択する、すなわちk空間の領域を分割するためのウィンドウを決める処理2。具体的には、最初に位相エンコード方向の被検体の腹部の両縁位置ymin及びymaxを検出する。仮にナビエコーを収集して1次元の動きのプロファイルを測定し、これを使用する場合には、動きの空間分布のインデックスΔYyが演算される。これにより、y軸方向、すなわち位相エンコード方向のプロファイルWyが[数5]Wy=ΔYy/{Δymax-Δymin}又は、ΔYは位相Φに置換してもよいとして得られる。必要であれば、このプロファイルWyをスムージングしてプロファイルWyとしてもよい。したがって、k空間データを位相エンコード方向ここではy軸方向:図2参照に例えば2分割する場合、その2分割を行なうウィンドウ関数WL,M,Syは、[数6]WLy=WyWSy=1-WLyとなり、仮にk空間データを3分割する場合図6参照、ウィンドウ関数WL,M,Syは、[数7]WLy=2max[Wy-0.5,0]WMy={1-WLy}Wy>=0.5のとき=WyotherwiseWSy=1-WLy-WMyとなる。なお、図6は、k空間データを3分割するときの、ウィンドウ関数WL,M,Syの設定を説明するものである。同図は、r空間のy方向での動き分布を表す指標として、その縦軸を位相Φで示している。この縦軸は位置シフトΔYで表すようにしてもよく、何れであってもよい。次に、演算ユニット10は、k空間は3分割するものとすると、求めたウィンドウ関数WL,M,Syと実空間の画像データPimgx,yを使って、k空間のデータを動きの大きいlarger領域、中程度mediumの領域、及び小さいsmaller領域のデータSLKx、Ky、SMKx、Ky、及びSSKx、Kyにそれぞれ分割する。すなわち、k空間のデータが複数領域のデータに分割される処理3。それには、[数8]PLx、y=WLyPimgx,y動きの大きい領域の画像データに対してPMx、y=WMyPimgx,y動きの中程度の領域の画像データに対してPSx、y=WSyPimgx,y=Pimgx,y-PLx、y-PMx、y動きの小さい領域の画像データに対してSLKx、Ky=IFTx→Kx,y→Ky[PLx、y]SMKx、Ky=IFTx→Kx,y→Ky[PMx、y]SSKx、Ky=IFTx→Kx,y→Ky[PSx、y]=SKx、Ky-SLKx、Ky-SMKx、Kyの演算が行なわれる。次いで、演算ユニット10により、ショットn毎に、すなわちショットn=1~Nmaxまで順次繰り返して以下の処理4~6が実行される。まず、ナビエコーSnaviKy,nから、動きに因って生じたk空間ΦKy,nにおける位相分布が求められる処理4。[数9]Pnaviy,n=FTKy→y[SnaviKy,n]ΔYy,n=Yy,n-Yy,nbaseΔYmeann=∫ΔYy,ndyΔYmaxn=maxof[ΔYy,n]ここで、ベースのエコーデータはk空間の中心データn=Nmax/2又は全てのΔYy,nの平均に基づいて選択される。次いで、k空間における位相分布の最大値ΦmaxKy,nが演算される処理5。具体的には、線形モデル及びゼロ次位相を用いる場合、[数10]ΔΦ0Ky,n=2πKyΔYmeannΔΦmaxKy,n=2maxof[ΔΦ0Ky,n]となり、y軸方向の1次元投影のプロファイルΔYy,nを用いる場合、[数11]ΔΦmaxKy,n=2πKyΔYmaxy,nとなる。次いで、同一ショットにおけるKyに対して位相補正が行なわれる処理6。具体的には、[数12]SL.corKx,Ky=SLKx,Kyexp[-iΔΦmaxKy,n]SM.corKx,Ky=SMKx,Kyexp[-iΔΦmaxKy,n/2]となる。これで、ショットn毎のループは終る。次いで、演算ユニット10は全ての分割k空間のデータを下記式によって合成する処理7。[数13]Simg.corKx,Ky=SL.corSM.corSSこのようにして求めたk空間全体の補正データSimg.corKx,Kyを2次元フーリエ変換により実空間の画像データに再構成する処理8。これを数式化すれば、[数14]Pimg.corx,y=FT2DKx→x,Ky→y[Simg.corKx,Ky]となる。このように、本実施形態によれば、r空間で計測又は予測した腹部の動きの情報を用いてk空間データを動きの程度に鑑みて分割し、かつ、かかる動きの程度に応じた互いに異なる線形の動き補正を分割データ毎に実施し、補正された分割データを最終的に合成するようにしている。ここで、線形の補正とはアフィン変換で定義される補正の全てが含まれる。つまり、従来では複雑で大量の演算が必要であって実質的に実施困難であった、撮像部位の非線形の動きに因って劣化したデータを、動きの程度に応じて画像データを分割し、各分割データに線形の動き補正を行い、この補正後の分割データを合成することで代替的に非線形のデータ補正を実現している。この代替的な非線形補正により、マルチショット型のMRイメージングを行なう場合でも、腹部などの空間的な非線形は動きに因って空間的に劣化したデータを容易に且つ高速に位相補正することができる。この結果、かかる非線形な動きが影響して空間的な不均一な非線形の位置シフト及び位相シフトが生じることに因るアーチファクトの補正を簡便にかつ高速に行なうことができ、画質を向上させることができる。ここで、前述した図3のパルスシーケンスについて補足しておく。このパルスシーケンスは、ナビエコーでショット毎にイメージングの位相エンコード方向に沿った動きに因るr空間で0次位相シフトのみを測定する場合のシーケンスである。ナビエコーSnaviKy,nの前に位相エンコード方向Geに強度Gm、時間Tmの傾斜磁場ナビゲータを加えて、r空間での平均位相k空間では中心の位相に相当を測定し、ベースショットとの差を求める。[数15]Φnavi0n=arg[SnaviKy=0,n]ΔΦnavi0y,n=Φnavi0n-Φnavi0nbaseなお、Φnavi0n=2πGmTmΔY0nの関係を用いてr空間での平均シフト量ΔY0nにも変換可能である。また、ナビエコーSnaviKy,nに対するA/Dサンプリングは、k空間の中心DCの位相のみを測定できればよいので、狭くてもよい。ナビエコー用のパルス列部分はSE法ではなく、時間短縮できるFE法でも信号さえ十分確保できればよい。それは、ベースショットによるエコーとの差分をとるため、FE法で生じるサスセプタビリティ効果susceptibilityeffectsや静磁場由来の位相シフトはキャンセルするためである。第2の実施の形態次に、本発明の磁気共鳴イメージング装置に係る一実施形態を、図1,7を参照して説明する。なお、本実施形態に係る磁気共鳴イメージング装置のハード的な構成は前述した図1に記載のものと同一であるので、その説明は省略する。本実施形態は、ハード的には受信RFコイル7Rにマルチコイルを使用し、このマルチコイルを成す2つのコイルエレメントそれぞれの感度領域を利用したデータ収集と本願発明の補正処理とをリンクさせて、効率的な補正を実現するものである。とくに、腹部を撮像する場合を考えると、腹壁側からのゴースト動きに因る非線形成分が背中側に及んでいる場合、動きに関する図2の仮定が正しいとしても、その影響は最終的に無視できるかどうかという懸念が残る。本実施形態は、そのような懸念を確実に払拭できるようにする。図7aに示すように、この受信RFコイルとしてのマルチコイル7Rは、腹部の背中側と腹壁側とにそれぞれ相対して配置する2つの表面コイル7RL、7RUをエレメントコイルとする。このようなマルチコイルと専用の画像再構成を組合せた撮像は、撮像時間の短縮とSNRを向上させるために最近の磁気共鳴イメージングでは広く実用化されている。腹部を撮像する場合には、図7aに示すように、背中側と腹壁側とに2個配置する場合が多い。このように表面コイル7RU、7RLを配置することで、一方の腹壁側の表面コイル7RUはその腹壁側により強い検出感度を有する。このため、表面コイル7RUはからの信号は腹壁側に重みが大きくかかっており、その腹壁側の動きを支配的な成分として検出する。他方の背中側の表面コイル7RLはその背中側により強い検出感度を有する。このため、この表面コイル7RLからの信号は背中側の動きが小さい成分を支配的に検出する。つまり、この2つの表面コイル7RU、7RLを腹部の腹壁側及び背中側に沿ってそれぞれ配置することで、都合良く、最初から、腹部の非線形の動きの程度を反映して分割された2組のk空間データがそれぞれ収集される図7b、c参照。この2組のk空間データは、前述した第1の実施形態における図5e、fにそれぞれ対応するものであるから、図5c、dに対応するフーリエ変換及びウィンドウ処理は不要になる。図7d~hの処理は、前述した図5g~kの処理と同じである。このように配置された2つの表面コイル7RU、7RLによって、空間的に異なる重みがつけられたエコーデータを被検体Pの腹部から直接に収集できる。とくに、腹部の診断の場合の受信RFコイルは従来から腹壁側と背中側とに表面コイルを配置するようになっていることが多い。このため、2つの表面コイル7RU、7RLから腹壁側の大きい動きの影響を支配的に受けたエコーデータ及び背中側の小さい動きの影響を支配的に受けたエコーデータが最初から収集される。そこで、腹壁側からのエコーデータのみについて前述した補正処理を行い、それぞれのk空間データを再構成して相互に加算し、最終的な実空間の画像データPimg.corx,yを得ることができる図7d~h参照。したがって、この第2の実施形態によれば、前述した第1の実施形態とは異なり、腹壁側の動きに因り背中側までゴーストが及んでいたとしても、ゴーストの発生元信号は腹壁側の表面コイル7RUからのエコーデータのみに重点的に含まれている。つまり、反対側の表面コイル7RLからのエコーデータに含まれるゴーストの影響は最初から抑制されている。このため、腹壁側の表面コイル7RUからのエコーデータをそのまま補正すれば、腹壁側から遠い背中側に及んだゴーストを含めて補正したことと等価になり、第1の実施形態で採用した単一コイルの場合に比べて、処理時間的にも補正精度的にもより向上した画質の再構成像を最終的に得ることができる。なお、上述したコイル配置において、図8に示すように、エレメントコイル7A~7Dを腹部の背中側及び腹壁側それぞれにて上下方向頭部から脚部に向かう方向に分割して配置した4個の配置構成とすることもできる。この場合には、横隔膜付近にあたる腹部の上側頭部側の動きが大きいので、4個をそれぞれ別々の補正にするとよい。それにより、コロナル面やサジタル面の画像に対して補正効果が大きい。さらに、第1の実施形態に係る単独の表面コイルと、上述したエレメントコイル4個とを組み合わせて配置する14個の配置構成を採ることもできる。この場合も、横隔膜付近にあたる腹部の上側頭部側の動きが大きいので、分割配置の4個のエレメントコイルそれぞれに別々の補正を施すようにしてもよく、これにより、コロナル面やサジタル面に対して大きな補正効果が得られる。さらに、撮像部位は腹部以外でもよい。マルチコイルを成すコイルエレメントから出力される信号の相互間で撮像部位の動き成分が、最初からある程度分離できている場合には腹部以外の撮像部位であってもよい。例えば、撮像部位が心臓の場合には、同一心拍位相においても呼吸性の動きに因り、心臓の下部と上部では動きが非線形であることは知られている。この場合に、本第2の実施形態に係るマルチコイルによる撮像を実施するには、体壁から心臓を囲むように上下左右の複数の表面コイルを配置する。これにより、各表面コイルはそのコイル感度に応じて動きの程度を分離した状態で検出できので、本第2の実施形態に係る手法を実施することができる。腹部の撮像において、腹壁側からのゴーストが背中側に及んでいる場合の効果は、本発明者が実施したシミュレーションによっても既に確認されている。なお、マルチコイルを用いた動きの空間分布の測定は、各コイルエレメントが出力したナビエコーから位置の平均シフト量を演算し、その平均シフト量を補正に用いてもよい。そのようにすれば、コイルの位置関係に関する情報を同定しなくても補正可能であり、部位やコイル形態を問わずに一般的な応用が可能になる。第3の実施形態第3の実施形態に係る磁気共鳴イメージング装置は、マルチコイルに基づく処理つまり、上述した第2の実施形態で説明される構成と、画像処理によるウィンドウ処理とによる動き成分の分割つまり、前述した第1の実施形態で説明される構成とを組み合わせて実施するものである。すなわち、第3の実施形態に係る磁気共鳴イメージング装置では、動き成分を細分化するために、各コイルエレメントからの画像それぞれにウィンドウ処理を施してデータを分割する。図9は、本発明の第3の実施形態に係る磁気共鳴イメージング装置において受信RFコイルとして用いられるシングルコイルおよびマルチコイルとその感度プロファイル並びに被検体の非剛体モーションの大きさを示す図である。図9aに示すように、第3の実施形態では、図7aに示す第2の実施形態と同様に、受信RFコイルとしてのマルチコイル7RMが用いられる。図9aに示すようにマルチコイル7RMは、2つの表面コイル7RU、7RDをエレメントコイルとしており、各表面コイル7RU、7RDは、それぞれ被検体Pの腹部の背中側と腹壁側とにそれぞれ相対して配置される。また、必要に応じて、受信RFコイルとして筒状のシングルコイル7RSも被検体Pを覆うように設けられる。そして、シングルコイル7RSおよびマルチコイル7RMの任意の一方または双方を受信用RFコイルとして使用することができる。図9bには、被検体Pの腹部における非剛体モーションの大きさを示す。図9bの縦軸は図9aに示す画像データの位置y[pixel]を示し、横軸は、ピクセルyにおける動きの大きさAy[pixel]を示す。図9中の実線Ayに示すように、動きはyuとydとの間においてy方向の1次分布をとるものとする。すなわち、動きの大きさAyは、yuにおいて最大値Amaxをとり、ydにおいて0であるものとする。図9cには、各表面コイル7RU、7RDおよびシングルコイル7RSの正規化されたy方向の各感度分布Wcuy,Wcdy,Wcsyを示す。図9cの縦軸は図9aに示す画像データの位置y[pixel]を示し、横軸は、各表面コイル7RU、7RDおよびシングルコイル7RSの正規化されたy方向の各感度分布Wcuy,Wcdy,Wcsyを示す。図9cに示すように、各表面コイル7RU、7RDの感度分布Wcuy,Wcdyはそれぞれ動き分布と同様にy方向に1次分布するものとし、シングルコイル7RSの感度分布Wcsyはy方向に一定であるものとする。そして、このような各表面コイル7RU、7RDまたはシングルコイル7RSにより収集されたk空間データを、ウィンドウ関数により複数の領域に分割することができる。特に、2つの表面コイル7RU、7RDにより収集されたk空間データは、ハード的に2つの領域に分割されているため、各領域のうち一方または双方を更にウィンドウ関数を用いてソフト的に複数の領域に分割することができる。また、シングルコイル7RSにより収集されたk空間データは、ウィンドウ関数を用いてソフト的に任意数の領域に分割することができる。例えば、シングルコイル7RSにより収集されたk空間データは、第1の実施形態と同様なウィンドウ関数を用いてソフト的に2分割または3分割することができる。図10は、図9に示すような1次元の動き分布と同様な感度分布を有する2つの表面コイル7RU、7RDを用いて収集された各k空間データをそれぞれ2つの領域に分割することにより等価的に3分割する場合におけるウィンドウ関数を示す図である。図10において、縦軸は、ウィンドウ関数WLWMWSの値を示し、横軸は図9に示すy方向のピクセル位置を示す。図10に示す3つのウィンドウ関数WLWMWSを用いて2つの表面コイル7RU、7RDにより収集された2つのk空間データを動きが大きいL領域、動きが中程度のM領域および動きが小さいS領域の3つの領域に分割することができる。そして、分割された各領域のk空間データに対して、動きの大きさに応じた強度の補正を施すことができる。つまり、2つの表面コイル7RU、7RDによりハード的に2つに分割されたk空間データを、ウィンドウ関数を用いてソフト的に更に2つの領域に分割するとともに、4つの領域のうち動きが中程度の2つの領域に含まれるk空間データが単一の領域に含まれるように領域の一部を合成することができる。このように、複数のエレメントコイルにより分割されて収集されたk空間データの一部を他の一部と合成して領域を形成してもよい。尚、複数のエレメントコイルにより、ハード的に分割された複数のk空間データに限らず、ウィンドウ関数により一旦分割されたk空間データの一部を同様に他の一部と合成した後、補正を施すようにしてもよい。補正前にk空間データの一部の合成処理を行うことにより、より標準的なウィンドウ関数を用いて、k空間データを所望の不均等な複数の領域に分割することが可能となるのみならず、処理を容易にすることができる。図10に示す各ウィンドウ関数は、それぞれ次式により正規化された動き分布Wおよび表面コイル7RU、7RDの感度分布Wcから求めることができる。[数16]WLWcW≒W2≒WWMWc1W1WcW≒2W1W≒WWS1WLWM≒1W2≒Wそして、第1の実施形態においてソフト的にk空間データを3分割する場合と同様に、上記のウィンドウ関数WLWMWSを2つの表面コイル7RU、7RDにより得られた実空間の各画像データPimguPimgdに乗じて、動きが大きい領域、中程度の領域および小さい領域の3つの画像データPLPMPSを生成する。そして、生成された各画像データを2次元逆フーリエ変換することにより、動きの大きさに応じて3つの領域に分割されたk空間データSLSMSSを得ることができる。さらに、各k空間データSLSMSSに対して動きの大きさに応じて異なる強度で位相補正が実施され、位相補正後の各k空間データSLcorSMcorSSを合成したk空間データSimgcorを2次元フーリエ変換することにより、補正後の実空間の画像データPimgcorを得ることができる。このように第3の実施形態によれば、複数のエレメントコイルおよびウィンドウ関数を用いて、k空間データをハード的かつソフト的に分割することができる。従って、第1の実施形態および第2の実施形態により得られる双方のメリットを享受することができる。また、ハード的またはソフト的に一旦、分割されたk空間データの一部を合成することにより、簡易な処理で容易にk空間データを所望の領域に分割した上で動きの大きさに応じた強度の補正を行うことが可能となる。以上のような各実施形態におけるデータの補正処理の共通原理について一般化して説明する。図11には、上述した各実施形態で説明した補正処理の原理を示す。まず、非剛体の動きの情報Input1を何らかの方法で取得する。例えばナビエコーから動きを測定しても良いし、モデル化することにより動きの情報Input1を得ても良い。次に、動きに起因して空間的に劣化した画像データInput2を、動きの情報Input1を用いて、動きの大きさ応じて空間的に複数の成分に分割する。そして、各画像データからそれぞれ分割された複数のk空間データを得る。このk空間データの分割は、ウィンドウ関数によりソフト的に行ってもよいし、マルチコイルを用いてハード的に行ってもよい。また、マルチコイルを用いてハード的に分割されたk空間データをウィンドウ関数によりソフト的に更に分割してもよい。次に、N個に分割されたk空間データの各成分i1≦i≦Nに対して、動きの大きさ応じた線形変換により線形補正を行う。この線形変換には、affine変換、回転移動、拡大縮小、せん断変形が含まれる。すなわち、より動きが大きい成分に対してより強度が強い線形補正を行う一方、より動きが小さい成分に対してより強度が弱い線形補正を行う。動きが小さい成分に対しては、強度をゼロとした補正、すなわち非補正としてもよい。次に、補正後におけるk空間データの各成分の合成処理を行った後、合成後のk空間データからFT処理により全体的な補正後の画像データOutputが得られる。また、図12には、上述した各実施形態で説明した補正処理の流れを一般化して示す。最初に、磁気共鳴イメージング装置により、撮像部位をスキャンしてイメージングを行い、エコーデータを収集するステップS1。このエコーデータの収集と一緒に又は単独にナビエコーデータを収集するステップS1 ́。なお、このエコーデータ及びナビエコーデータは既に収集済みのものを採用してもよい。次いで、かかる処理は磁気共鳴イメージング装置の演算ユニット10に任される。つまり、演算ユニット10により、種々のコントロール情報動きの分割領域数、ウィンドウ処理の有無、撮影部位、体位、コイルの種類シングルコイルかマルチコイルか、コイル配置位置、パルスシーケンスの種類などに基づいて、本発明に係る、腹部などの非線形な動きに応じた補正処理を伴った画像再構成が実施される。まず、演算ユニット10は、動きの程度に応じたソフト的なデータ分割が必要か否かを判断するステップS2。この判断は、使用する受信RFコイルの種類、コイル配置位置などを考慮してなされる。この判断がYESとなる場合には、前述した第1の実施形態で説明した2次元フーリエ変換ステップS3:図5の処理c参照、ウィンドウ処理ステップS4:図5の処理d参照、2次元逆フーリエ変換処理ステップS5:図5の処理e,f参照、r空間上での動き情報の検出ステップS6:図5の処理g、h参照、k空間上のデータの非線形な動きによる劣化の補正ステップS7:図5の処理i参照、及び、2次元フーリエ及び加算処理による画像生成ステップS8:図5の処理j、k参照が実行される。一方、ステップS2に判断がNOとなる場合、前述した第2の実施形態で説明したマルチコイルを受信RFコイルとして用いた処理が実行される。つまり、ステップS3~S5の処理を割愛して、直接に、r空間上での動き情報の検出ステップS6:図7の処理d、e参照、k空間上のデータの非線形な動きによる劣化の補正ステップS7:図7の処理f参照、及び、2次元フーリエ及び加算処理による画像生成ステップS8:図7の処理g、h参照が実行される。このような流れに沿って、空間的に非線形な動きの影響を受けた収集データ又は空間的に非線形の動きの直接反映させたデータを収集し、非線形成分を線形成分の和に置換して、比較的簡単な少ない処理でかつ高速に、撮像部位の非線形な動きの影響を確実に排除したMR画像を再構成することができる。ところで、近年、心臓に対する磁気共鳴イメージングの分野において、前述した非特許文献8に見られるように、直前にモニタしたナビエコーから横隔膜の3次元的な位置をリアルタイムに計算して傾斜磁場強度やRF周波数を制御することにより、スラブ面を心臓の動きに合わせて3次元的に追随させ、モーションアーチファクトやボケを低減する手法が報告されている。この手法は、線形変換であるaffine変換、拡大縮小、せん断変形の拡大縮小に限れば、後処理も不要であり画像劣化もないので優れた方法ではある。この手法を本発明の手法と比較した場合、本発明の手法は、まず、線形変換に限らず非線形の動きによるシフトを補正することができる。すなわち、動きによる位置シフトが線形変換では不可能な高次の空間分布にも適用できるという優位さがある。また、本発明の手法は、ウィンドウ処理が必要な場合でも、補正前の画像再構成後に分割するので、後処理ではあるが、演算は簡単になり、収集終了後、即座に補正された画像を出力できる。さらにマルチコイルを用いた場合、演算はさらに簡単化され、ショット単位で動きのナビゲーションを本イメージングの直前で実施すれば、リアルタイムに補正処理を行うことができるという優位性もある。なお、上述した各実施形態では1次元の補正処理、つまり動きの方向がY方向のみで振幅が1次元の分布である場合について説明したが、動き情報の測定検出及びその補正処理を3次元的にx、y、zの各方向で行うようにしてもよい。特に、動きの情報は、k空間における位相分布又はr空間における少なくとも1方向の位置のシフトに基づいて測定できる。ここで、3次元的に非剛体の動きを補正する場合について説明する。非剛体のモーションを一般化し、動きによるシフトの方向および大きさが共に場所毎および時間毎に異なる場合には、すなわちシフトの方向と変位が4次元的な分布である場合には、データを分割する重み関数を動きによるシフトの分布に応じて高次元にすればよい。具体的には、変位方向が3方向ΔX、ΔY、ΔZに混在している場合には、shot毎に測定または推定されたx、y、zの3方向の別々の独立した動きの大きさに応じて、3方向について別々にwindow関数を求めることにより動きの大きさに応じて分割し、補正を加法的に実施すればよい。例えば、window関数を動きの方向別に動きの大きさ毎に作成する。そして、各方向について動きの大きさにより分割された成分毎に補正した後、加算することによりデータ全体の補正を行うことができる。この場合のアルゴリズムの例を以下に示す。図13は、本発明により3次元的に非剛体の動きを補正する場合の手順を示すフローチャートであり、図中Sに数字を付した符号はフローチャートの各ステップを示す。また、図14は、被検体Pにおいて、非剛体の3次元的な動きがある場合にy方向からみた動き方向の例を示す図であり、図15は、図14に示す被検体Pにおける非剛体のx方向からみた動き方向を示す図である。まず、ステップS10において、動きによるある位置のシフトにおけるshotn毎のある基準からの変位ベクトル図14参照の3方向別の空間分布ΔX,ΔY,ΔZを得る。次にステップS11において、次式により3D空間における3方向での最大変位ΔXmax,ΔYmax,ΔZmaxを算出する。[数17]ΔXmaxmaxΔXΔYmaxmaxΔYΔZmaxmaxΔZ次にステップS12において、次式により3方向各々の3D分割ウィンドウ関数重み関数WxWyWzを算出する。[数18]WxΔXΔXmaxWyΔYΔYmaxWzΔZΔZmax次にステップS13において、次式により3方向成分について、各方向のウィンドウ関数で実空間の画像データPを分割する。尚、動きの大きさ別にL大きいM中程度S小さいの添え字をつけるものとする。[数19]PxLMSPWxLMSPyLMSPWyLMSPzLMSPWzLMS次にステップS14において、次式に示すように、分割された画像データPxLMSPyLMSPzLMSを逆フーリエ変換してk空間データSxLMSSyLMSSzLMSを得る。[数20]SxLMSIFTPxLMSSyLMSIFTPyLMSSxLMSIFTPyLMS次にステップS15において、次式に示すように、k空間データについて3成分毎に位相補正を行う。[数21]動き最大部の補正SxcorLSxLexpiΦxmaxSycorLSyLexpiΦymaxSzcorLSzLexpiΦzmax動き中間部の補正SxcorMSxMexpiΦxmax2SycorMSyMexpiΦymax2SzcorMSzMexpiΦzmax2但し、Φxmax2πkxΔXmaxKxΦymax2πkyΔYmaxKyΦzmax2πkzΔZmaxKzである。次にステップS16において、次式に示すように、k空間データの補正成分SxcorLSycorLSzcorLSxcorMSycorMSzcorMおよび非補正成分SxSSySSzSを合成する。具体的には、補正成分における各動き成分を3方向別に加算した後、加算後における3方向の成分SxcorSycorSzcorを更に加算する。[数22]SxcorSxcorLSxcorMSxSSycorSycorLSycorMSySSzcorSzcorLSzcorMSzSScorSxcorSycorSzcorこのようにして3次元的な補正を施したk空間データScorを得ることができる。さらに、ステップS17において、合成後におけるk空間データの3DFTにより3次元的な動きの補正が施された3D画像データPcorを得ることができる。ところで、以上のような各実施形態において、動き情報を検出するための手法は、傾斜磁場をナビゲーションパルスとして印加して収集するナビエコーを用いるものに限定されることなく、外部の光学、磁気、機械的センサを用いたものでもよい。その場合に、かかるセンサで動きの空間分布までは計測できなくても、腹壁のみの動きをモニタして、事前に計測又は推定した体動分布のモデリング結果と組み合わせて空間分布を推定するようにしてもよい。これにより、磁気共鳴イメージングのデータ収集と並行した計測が可能になるので、心臓などの高速イメージングにも対応可能になる。また、線形変換で補正可能な成分のみ、高周波や傾斜磁場などのオンザフライontheflyでの制御で行い、残った2次以上の非線形な成分のみを本発明に係る手法で補正するといった組合せで実施することも可能である。さらに、前述した補正処理において、測定した位相誤差が装置として対応し得る所定の限界値を超えている場合には、その限界値で閾値処理をするようにしてもよい。一方、上述した各実施形態において、補正すべき画像データがあるいは必要な画像データが分割された領域の一部であるような場合には、必ずしも分割された補正後の画像データを合成する必要はない。加えて、上述した各実施形態では、r空間上の画像データを分割した後、分割された各画像データを変換することにより複数のk空間データを生成したが、r空間上のデータの性質を引き継ぐことが可能な変換によりr空間上のデータから別の空間上のデータを生成し、この別の空間上のデータを変換することによりk空間データを生成するようにしてもよい。また、逆に、k空間上のデータの性質を引き継ぐことが可能な変換により得られる他の空間のデータにr空間上のデータから一旦変換し、この他の空間のデータからk空間データを生成するようにしてもよい。以上の実施形態及び変形例に共通する代表的な効果は、以下のようにまとめることができる。1次元又は2次元/3次元の動き位置シフト又は位相シフトの空間分布が測定できれば、少なくとも2つの異なる補正データの線形和を演算することで、非線形な動きに因る劣化を受けたデータの位相や信号値を比較的簡単に補正できる。位置のシフトの補正の場合は、r空間でボクセルを跨いで位置シフトの計測および補正処理する場合に比べて、k空間での同一波数での測定および位相補正であるから高速に処理が可能になる。ボクセル内の位相補正の場合はr空間での補正が適しているが、ボクセル毎に処理する場合に比べて、簡便な演算でかつ高速化が可能になる。また、動きの空間分布を厳密に測定しなくても、全体の平均や腹壁など典型的な部分の動きの測定データから、非線形な場合でも全体の空間分布をモデル的に推定するので測定系が簡便になるナビエコーで測定する場合はシーケンスも短く高速化が可能になる。例えば、腹部のアキシャル断面なら背中から動きに因る位置シフトが腹壁側へ漸増していると仮定できる。さらに、単一コイルを用いた場合でも、k空間データの分割数分だけのデータの再構成を行なえばよいので、r空間で補正するよりは高速化が可能になる。加えて、腹部をマルチコイルを用いて撮像する場合、エコーデータを直接取得できるので、補正精度が向上するのみならず、ウィンドウ分割が不要なので、より一層の高速化を図ることができる。ここで、本発明により補正を行って得られた被検体Pの腹部における断層画像を補正を行わずに得られた従来の断層画像と比較して示す。図16は、位相エンコードPEphaseencode方向と垂直な方向に1次元の分布を有する非剛体の動きがあった場合に、本発明によりシングルコイルを用いてソフト的にk空間データを分割して補正することにより得られた断層画像を示す図である。図16aは、補正を行わずに得られた画像、図16bは、従来の一様平均シフト補正により得られた画像、図16cは、本発明によりソフト的に2分割された各k空間データをそれぞれ線形補正して得られた画像、図16dは、本発明によりソフト的に3分割された各k空間データをそれぞれ線形補正して得られた画像である。図16aに示す補正前の画像に対し、従来の一様平均シフト補正を施すと、図16bに示すように補正量に一致した被検体Pの前後方向における中心部のゴーストは低減されるものの補正が不十分であると言える。一方、図16cに示すように、本発明によりデータをソフト的に2分割してそれぞれ異なる強度で線形補正を施すと、ウィンドウの中間被検体Pの前後方向中間部でSNRの低下とゴーストの補正不足が見られるものの、従来の一様平均シフト補正により得られた画像と比較して顕著にゴーストが低減されているのが確認できる。さらに、図16dに示すように、本発明によりデータをソフト的に3分割してそれぞれ異なる強度で線形補正を施すと、被検体Pの前後方向中間部におけるSNRの低下およびゴーストが十分に補正され、顕著な補正効果が確認できる。図17は、PE方向に1次元の分布を有する非剛体の動きがあった場合に、本発明によりシングルコイルを用いてソフト的にk空間データを分割して補正することにより得られた断層画像を示す図である。図17aは、補正を行わずに得られた画像、図17bは、従来の一様平均シフト補正により得られた画像、図17cは、本発明によりソフト的に2分割された各k空間データをそれぞれ線形補正して得られた画像、図17dは、本発明によりソフト的に3分割された各k空間データをそれぞれ線形補正して得られた画像、図17eは、動きの1次分布を示す図縦軸:ピクセル位置、横軸:正規化された動きの大きさである。図17aに示す補正前の画像に対し、従来の一様平均シフト補正を施すと、図17bに示すようにゴーストは新たに発生したものも含めて前後に同程度に生じているのが分かる。一方、図17cに示すように、本発明によりデータをソフト的に2分割してそれぞれ異なる強度で線形補正を施すとゴーストが低減されているのが確認できる。さらに図17dに示すように、本発明によりデータをソフト的に3分割してそれぞれ異なる強度で線形補正を施すと、ゴーストが更に低減されるとともに、腹壁と背中の中間部分における信号の低下が改善されているのが確認できる。図18は、PE方向に1次元の分布を有する非剛体の動きがあった場合に、本発明によりシングルコイルを用いてソフト的に空気部分を含めてk空間データを分割して補正することにより得られた断層画像を示す図である。図18aは、本発明によりソフト的に空気部分を含めて2分割された各k空間データをそれぞれ線形補正して得られた画像、図18bは、本発明によりソフト的に空気部分を含めて3分割された各k空間データをそれぞれ線形補正して得られた画像である。空気部分では動く部分から飛んだゴーストが支配的である。このため、空気部分を含めて動きが最大Lのデータとして画像データを補正すると、図18aおよび図18bに示すように、空気部分を含めない図17cおよび図17dと比較してゴーストが低減された。図19は、PE方向に1次元の分布を有する非剛体の動きがあった場合に、本発明により2つのエレメントコイルを有するマルチコイルを用いてハード的に2分割したk空間データをウィンドウ関数を用いてソフト的に更に2分割し、等価的に3分割されたk空間データを線形補正することにより得られた断層画像を示す図である。図19aは、2つのエレメントコイルによりハード的に2分割したk空間データをソフト的に分割せずに線形補正を行って得られた画像であり、図19bは、2つのエレメントコイルによりハード的に2分割したk空間データをソフト的に更に2分割し、等価的に3分割されたk空間データを補正することにより得られた画像である。図19aは、2つのエレメントコイルにより得られたデータのうち、動きが支配的なデータ側のみを補正して得られた画像であるが、シングルコイルを用いてソフト的に2分割されたデータを補正して得られた図17cの画像と比較して、空間的に遠くに飛んだ背側のゴーストが低減されている。また、図19bは、2つのエレメントコイルにより得られたデータをソフト的に2分割することにより等価的に3分割されたデータを補正して得られた画像である。図19bによれば、更に中間部のSNRの低下が改善され、図16a,b,c,d、図17a,b,c,d、図18a,b、図19aに示す他の画像と比較して最も良好な補正結果が得られているのが確認できる。尚、補正前の画像、一様平均シフト補正により得られた画像およびPE方向に垂直な読出し方向に動きが1次分布する場合の画像に対する補正では、シングルコイルを用いた場合とマルチコイルを用いた場合とで差異はない。なお、本発明は上記実施形態に記載のものに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。なお、心臓を画像化する場合、心拍に拠る動きと呼吸性の動きとの両方の影響を受けるが、心拍ゲーティングに加え、ナビエコーで呼吸周期をモニタし、スライス面の位置をダイナミックにコントロールする方法が採られている。"