"印加軸が同一で異なるb値に対応するMPGパルスの印加を伴って拡散強調イメージングを実行することによって前記異なるb値に対応する複数の拡散強調画像データを収集するイメージング手段と、前記複数の拡散強調画像データのうち歪補正の対象となる拡散強調画像データに対応するb値と異なるb値に対応する拡散強調画像データに基づいて前記歪補正の対象となる拡散強調画像データに対する歪補正係数を算出し、算出した前記歪補正係数を用いて前記歪補正の対象となる拡散強調画像データの歪補正を実行することによって前記歪補正が施された画像データを生成する歪補正手段と、を備える磁気共鳴イメージング装置。前記イメージング手段は、絶対値が同一で符号が互いに逆のb値に対応する2つの拡散強調画像データを収集するように構成され、前記歪補正手段は、前記2つの拡散強調画像データに基づいて前記2つの拡散強調画像データの一方に対する歪補正係数を算出するように構成される請求項1記載の磁気共鳴イメージング装置。前記イメージング手段は、絶対値が異なるb値に対応する複数の拡散強調画像データを収集するように構成され、前記歪補正手段は、基準となるb値に対応する歪補正係数に基づいて前記基準となるb値よりも大きいb値に対応する歪補正係数を算出するように構成される請求項1又は2記載の磁気共鳴イメージング装置。前記歪補正手段により算出された前記歪補正係数を記憶する記憶手段を更に備え、前記歪補正手段は、前記記憶手段を参照して前記歪補正係数を取得できるように構成される請求項1乃至3のいずれか1項に記載の磁気共鳴イメージング装置。前記イメージング手段は、b値をゼロとして更に画像データを収集するように構成され、前記歪補正手段は、前記b値をゼロとして収集された画像データにも基づいて前記歪補正係数を算出するように構成される請求項2記載の磁気共鳴イメージング装置。前記歪補正手段は、前記2つの拡散強調画像データ又は前記2つの拡散強調画像データに基づいて生成される2つの画像データに対する各座標変換係数をパラメータとして、前記2つの拡散強調画像データ又は前記2つの画像データの、b値をゼロとして収集された画像データに対する位置ずれの指標を最小にする最適化処理によって前記各座標変換係数を計算し、より小さい指標に対応する座標変換係数を前記歪補正係数とするように構成される請求項2記載の磁気共鳴イメージング装置。前記歪補正手段は、前記2つの拡散強調画像データの一方又は前記2つの拡散強調画像データに基づいて生成される2つの画像データの一方に対する座標変換係数をパラメータとして、前記一方の、b値をゼロとして収集された画像データに対する位置ずれの指標及び他方に対する位置ずれの指標を最小にする最適化処理によって計算された座標変換係数に基づいて前記歪補正係数を決定するように構成される請求項2記載の磁気共鳴イメージング装置。前記歪補正手段は、前記2つの拡散強調画像データの一方又は前記2つの拡散強調画像データに基づいて生成される2つの画像データの一方に対する座標変換係数をパラメータとして、前記一方の他方に対する位置ずれの指標を最小にする最適化処理によって計算された座標変換係数に基づいて前記歪補正係数を決定するように構成される請求項2記載の磁気共鳴イメージング装置。前記歪補正手段は、前記2つの拡散強調画像データに対するエッジ検出を含む処理によって前記歪補正係数を算出するように構成される請求項2記載の磁気共鳴イメージング装置。前記歪補正前及び前記歪補正後における所定の複数の拡散強調画像データの少なくとも一方の平均をとるアベレージング処理を行う画像処理手段を更に備える請求項1記載の磁気共鳴イメージング装置。前記画像処理手段は、前記歪補正後において互いに異なるb値に対応する複数の拡散強調画像データに対して前記アベレージング処理を行うように構成される請求項10記載の磁気共鳴イメージング装置。DWI及びDTIでは、画像の歪を一層良好に補正することが課題である。本発明は、より良好に歪補正を施したDWI又はDTIを得ることが可能な磁気共鳴イメージング装置を提供することを目的とする。本発明の実施形態に係る磁気共鳴イメージング装置は、イメージング手段及び歪補正手段を備える。イメージング手段は、印加軸が同一で異なるb値に対応するMPGパルスの印加を伴って拡散強調イメージングを実行することによって前記異なるb値に対応する複数の拡散強調画像データを収集する。歪補正手段は、前記複数の拡散強調画像データのうち歪補正の対象となる拡散強調画像データに対応するb値と異なるb値に対応する拡散強調画像データに基づいて前記歪補正の対象となる拡散強調画像データに対する歪補正係数を算出し、算出した前記歪補正係数を用いて前記歪補正の対象となる拡散強調画像データの歪補正を実行することによって前記歪補正が施された画像データを生成する。本発明の実施形態に係る磁気共鳴イメージング装置について添付図面を参照して説明する。図1は本発明の実施形態に係る磁気共鳴イメージング装置の構成図である。磁気共鳴イメージング装置20は、静磁場を形成する筒状の静磁場用磁石21、この静磁場用磁石21の内部に設けられたシムコイル22、傾斜磁場コイル23及びRFコイル24を備えている。また、磁気共鳴イメージング装置20には、制御系25が備えられる。制御系25は、静磁場電源26、傾斜磁場電源27、シムコイル電源28、送信器29、受信器30、シーケンスコントローラ31及びコンピュータ32を具備している。制御系25の傾斜磁場電源27は、X軸傾斜磁場電源27x、Y軸傾斜磁場電源27y及びZ軸傾斜磁場電源27zで構成される。また、コンピュータ32には、入力装置33、表示装置34、演算装置35及び記憶装置36が備えられる。静磁場用磁石21は静磁場電源26と接続され、静磁場電源26から供給された電流により撮像領域に静磁場を形成させる機能を有する。尚、静磁場用磁石21は超伝導コイルで構成される場合が多く、励磁の際に静磁場電源26と接続されて電流が供給されるが、一旦励磁された後は非接続状態とされるのが一般的である。また、静磁場用磁石21を永久磁石で構成し、静磁場電源26が設けられない場合もある。また、静磁場用磁石21の内側には、同軸上に筒状のシムコイル22が設けられる。シムコイル22はシムコイル電源28と接続され、シムコイル電源28からシムコイル22に電流が供給されて静磁場が均一化されるように構成される。傾斜磁場コイル23は、X軸傾斜磁場コイル23x、Y軸傾斜磁場コイル23y及びZ軸傾斜磁場コイル23zで構成され、静磁場用磁石21の内部において筒状に形成される。傾斜磁場コイル23の内側には寝台37が設けられて撮像領域とされ、寝台37には被検体Pがセットされる。RFコイル24にはガントリに内蔵されたRF信号の送受信用の全身用コイルや寝台37や被検体P近傍に設けられるRF信号の受信用の局所コイルなどがある。また、傾斜磁場コイル23は、傾斜磁場電源27と接続される。傾斜磁場コイル23のX軸傾斜磁場コイル23x、Y軸傾斜磁場コイル23y及びZ軸傾斜磁場コイル23zはそれぞれ、傾斜磁場電源27のX軸傾斜磁場電源27x、Y軸傾斜磁場電源27y及びZ軸傾斜磁場電源27zと接続される。そして、X軸傾斜磁場電源27x、Y軸傾斜磁場電源27y及びZ軸傾斜磁場電源27zからそれぞれX軸傾斜磁場コイル23x、Y軸傾斜磁場コイル23y及びZ軸傾斜磁場コイル23zに供給された電流により、撮像領域にそれぞれX軸方向の傾斜磁場Gx、Y軸方向の傾斜磁場Gy、Z軸方向の傾斜磁場Gzを形成することができるように構成される。RFコイル24は、送信器29及び受信器30の少なくとも一方と接続される。送信用のRFコイル24は、送信器29からRF信号を受けて被検体Pに送信する機能を有し、受信用のRFコイル24は、被検体P内部の原子核スピンのRF信号による励起に伴って発生したMR信号を受信して受信器30に与える機能を有する。一方、制御系25のシーケンスコントローラ31は、傾斜磁場電源27、送信器29及び受信器30と接続される。シーケンスコントローラ31は傾斜磁場電源27、送信器29及び受信器30を駆動させるために必要な制御情報、例えば傾斜磁場電源27に印加すべきパルス電流の強度や印加時間、印加タイミング等の動作制御情報を記述したシーケンス情報を記憶する機能と、記憶した所定のシーケンスに従って傾斜磁場電源27、送信器29及び受信器30を駆動させることによりX軸傾斜磁場Gx、Y軸傾斜磁場Gy,Z軸傾斜磁場Gz及びRF信号を発生させる機能を有する。また、シーケンスコントローラ31は、受信器30におけるMR信号の検波及びAD変換により得られた複素データである生データを受けてコンピュータ32に与えるように構成される。このため、送信器29には、シーケンスコントローラ31から受けた制御情報に基づいてRF信号をRFコイル24に与える機能が備えられる一方、受信器30には、RFコイル24から受けたMR信号を検波して所要の信号処理を実行するとともにA/D変換することにより、デジタル化された複素データである生データを生成する機能と生成した生データをシーケンスコントローラ31に与える機能とが備えられる。また、コンピュータ32の記憶装置36に保存されたプログラムを演算装置35で実行することにより、コンピュータ32には各種機能が備えられる。ただし、プログラムの少なくとも一部に代えて、各種機能を有する特定の回路を磁気共鳴イメージング装置20に設けてもよい。図2は、図1に示すコンピュータ32の機能ブロック図である。コンピュータ32の演算装置35は、記憶装置36に保存されたプログラムを実行することにより撮像条件設定部40及びデータ処理部41として機能する。データ処理部41は、画像再構成部41A、歪補正部41B及び画像処理部41Cを有する。また、記憶装置36は、k空間データ記憶部42、画像データ記憶部43及び歪係数記憶部44として機能する。撮像条件設定部40は、入力装置33から入力された情報に基づいてEPIシーケンスを含む撮像条件を設定し、設定した撮像条件をシーケンスコントローラ31に与えることにより駆動制御させる機能を有する。図3は、図2に示す撮像条件設定部40において設定されるEPIシーケンスの一例を示す図である。図3において横軸は時間を、RFはRF送信パルスを、GSSはスライス選択傾斜磁場パルスを、GPEは位相エンコード傾斜磁場パルスを、GROは読み出し傾斜磁場パルスを、ECHOはMR受信エコー信号を、それぞれ示す。図3に示すように、EPIシーケンスは、90°RFパルス及び180°RFパルスの印加によってSS傾斜磁場パルスの印加により選択されたスライスを励起した後、PE傾斜磁場パルス及びRO傾斜磁場パルスを繰り返し印加することによって、選択されたスライスから連続的に複数のMR受信エコー信号を収集するパルスシーケンスである。90°RFパルスの印加時刻からk空間中心におけるエコー信号の収集時刻までの時間は、実効エコー時間と呼ばれる。また、撮像条件設定部40は、MPGパルスの効果の指標であるb値を様々な値に設定してDWIを行うためのパルスシーケンスを設定する機能を有する。すなわち、撮像条件設定部40では、互いに異なるb値に対応する複数のDWIデータを収集するための撮像条件を設定することができる。特に、撮像条件設定部40は、印加軸及びb値の絶対値が互いに同一でかつ極性が互いに異なるMPGパルスの印加を伴う2つのDWIシーケンスを設定する機能も備えている。図4は、図2に示す撮像条件設定部40において設定される2つのDWIシーケンスの一例を示す図である。図4において横軸は時間を、RFはRF送信パルスを、GSSはSS傾斜磁場パルスを、GPEはPE傾斜磁場パルスを、GROはRO傾斜磁場パルスを、ECHOはMR受信エコー信号を、それぞれ示す。図4は、PE方向を印加軸とする極性が正のMPGパルスの印加を伴うDWIシーケンスの例を示す。また、図4は、PE方向を印加軸とする極性が負のMPGパルスの印加を伴うDWIシーケンスの例を示す。図4に示すように、MR受信エコー信号の収集、つまりPE傾斜磁場パルスの印加に先だって印加軸及び面積が同一で極性がそれぞれ正及び負のMPGパルスの印加を伴う2つのDWIシーケンスを撮像条件として撮像条件設定部40において設定することができる。尚、図4は、MPGパルスの印加軸がPE方向である例を示しているが、診断目的に応じて他の方向が印加軸とされる場合もある。また、DWIデータは、複数の画像データの平均をとるアベレージング処理によりAVERAGE画像データとして収集することができる。アベレージング処理を行えば、SNRを向上させたDWIデータを得ることができる。従って、平均をとるために加算される画像データの数であるAVERAGE数、加算対象となるDWIデータに対応するMPGパルスの印加軸及び加算対象となるDWIデータのb値が撮像条件として撮像条件設定部40により設定される。データ処理部41は、シーケンスコントローラ31からMR信号を受けてk空間データ記憶部42に形成されたk空間にk空間データとして配置する機能、k空間データ記憶部42に保存されたk空間データに対する画像再構成処理及び歪補正処理を含むデータ処理によって、歪補正を施した画像データを生成する機能、生成した画像データを画像データ記憶部43に書き込む機能、DWIデータの歪補正に用いられる歪補正係数を求めて歪係数記憶部44に書き込む機能、画像データを表示装置34に表示させる機能を有する。画像再構成部41Aは、k空間データ記憶部42からk空間データを取り込んでフーリエ変換を含む画像再構成処理を施すことにより画像データを再構成する機能を有する。MPGパルスを印加して収集された画像データはDWIデータとなり、各DWIデータは、b値及びMPGパルスの印加軸に対応する画像データとなる。尚、DWIデータではない画像データは、b値がゼロの画像データに相当する。また、画像再構成部41Aは、必要に応じてMPGパルスの印加軸及びb値を同一に設定して収集された複数のDWIデータの平均をとる第1のアベレージング処理によってAVERAGE画像データを生成する機能を有する。すなわち、同一のDWIの条件で収集されたDWIデータを加算することによって特定の印加軸及びb値に対応するDWIデータのSNRを向上させることができる。歪補正部41Bは、MPGパルスの印加を伴って収集されたDWIデータに対して歪補正を行うための歪補正係数を求める機能と、求めた歪補正係数又は歪係数記憶部44を参照して取得した歪補正係数を用いて歪補正の対象となるDWIデータに対する歪補正を行うことによって、歪補正が施された画像データを生成する機能を有する。特に、歪補正部41Bは、DWIによって収集された複数のDWIデータのうち歪補正の対象となるDWIデータに対応するb値と異なるb値に対応するDWIデータに基づいて、歪補正の対象となるDWIデータにに対する歪補正係数を算出するように構成される。また、歪補正係数の算出は、必要に応じてMPGパルスを印加せずにb値をゼロとして収集されたDWIでない画像データにも基づいて行うことができる。より具体的には、印加軸及びb値の絶対値が互いに同一でかつ異なる極性のMPGパルスに対応する2つのDWIデータに基づいて、より高精度に2つのDWIデータの一方に対する歪補正係数を求めることができる。図5は、図2に示す歪補正部41Bにおいて歪補正係数を求めて歪補正を行う第1の方法を説明する図である。図5において各縦軸はx軸を示す。例えば、b値を0500500に設定してそれぞれ画像データを収集すると、図5に示すように、b値をゼロとして収集された基準画像データIb0は、MPGパルスの影響がないため、理想的には歪のない画像データとなる。一方、b値が500に対応するDWIデータIb500と、b値が500に対応するDWIデータIb500は、それぞれMPGパルスの影響を受けて歪を有する画像データとなる。また、b値が500に対応するDWIデータIb500と、b値が500に対応するDWIデータIb500は、それぞれ互いに逆の極性のMPGパルスの影響を受けているため、他の影響が無視できれば、歪方向の極性が原点に関して互いに対称となる。尚、図5は、説明簡易化のためにDWIデータIb500Ib500にx軸方向の1次元の歪みが生じた場合の例を示している。b値が±500である場合に限らず、絶対値が同じで符号が逆のb値に対応する2つのDWIデータは、他の影響が無視できれば、歪方向が原点に関して対称で、歪量が同一となる。このような、0500500のb値に対応する各画像データIb0Ib500Ib500に対して、閾値を用いた2値化処理を行うことによって、図5に示すように、それぞれ0500500のb値に対応する2値化画像データIb0Ib500Ib500が生成される。次に、500のb値に対応する2値化DWIデータIb500の各画像点を座標変換して移動させながら、0のb値に対応する2値化基準画像データIb0の各画像点における画素値と500のb値に対応する2値化DWIデータIb500の各画像点における画素値との間の2乗差の総和が計算される。すなわち、2値化DWIデータIb500の座標変換を表す各係数をパラメータとして、2値化基準画像データIb0に対する2値化DWIデータIb500の画素値の2乗誤差を最小にする最適化処理が実行される。座標変換は、平行移動、回転移動及び伸縮変形の組み合わせて表されるアフィン変換とすることが処理の簡易化に繋がる。アフィン変換による座標変換は、式により表される。但し、式において、fは座標変換を表す関数、は座標変換前における画像点の空間座標、は座標変換後における画像点の空間座標、p0p1p2p8TxTyTzは座標変換の係数である。式の各座標変換係数をパラメータとして、2値化DWIデータIb500の画素値の2乗誤差を最小とする最適化処理を行うと、図5に示すように2値化基準画像データIb0に対する位置ずれが最小となる2値化補正DWIデータIb500corが生成されるとともに行列として複数の座標変換係数が得られる。すなわち、基準画像データIb0と正のb値に対応するDWIデータIb500間の2値化処理及び座標変換処理後における各画素値の差分値の総和が最小となる座標変換係数の組合せK1を得ることができる。尚、最適化処理は、座標変換係数の逐次探索処理及び画素値の補間処理を含む処理となる。同様な方法で、基準画像データIb0と負のb値に対応するDWIデータIb500間の2値化処理及び座標変換処理後における各画素値の差分値の総和が最小となる座標変換係数の組合せK2とともに2値化補正DWIデータIb500corを得ることができる。そして、2値化基準画像データIb0に対する画素値の2乗差の総和がより小さくなる2値化補正DWIデータの生成に用いられた座標変換係数の組合せを選択して歪補正係数として採用することができる。更に、2値化基準画像データIb0に対する画素値の2乗差の総和がより小さくなる2値化補正DWIデータの生成に用いられたDWIデータを、歪補正係数を用いて座標変換することにより、図5に示すように歪補正後のDWIデータIb500corを得ることができる。図6は、図2に示す歪補正部41Bにおいて歪補正係数を求めて歪補正を行う第2の方法を説明する図である。図6において各縦軸はx軸を示す。図5と同様に、図6に示すような歪のない基準画像データIb0及びb値が±500に対応する歪方向が互いに逆の2つのDWIデータIb500Ib500を収集して2値化処理を行うことによって、図6に示すように、それぞれ0500500のb値に対応する2値化画像データIb0Ib500Ib500を生成することができる。ここで、仮に、500のb値に対応する2値化DWIデータIb500の2値化基準画像データIb0に対する各画素値の2乗差の総和を最小にする最適な各座標変換係数を2倍にした座標変換係数の組合せ2K1を用いて500のb値に対応する2値化DWIデータIb500の座標変換を実行すると、座標変換によって生成される2値化変換DWIデータIb500doubleは、理論的には500のb値に対応する2値化DWIデータIb500と等しくなる。そこで、500のb値に対応する2値化DWIデータIb500の各座標変換係数をパラメータとして、500のb値に対応する2値化DWIデータIb500の2値化基準画像データIb0に対する各画素値の2乗差の総和SumError1と、2倍の座標変換係数の組合せで座標変換された2値化変換DWIデータIb500doubleと500のb値に対応する2値化DWIデータIb500との間における各画素値の2乗差の総和SumError2の和SumError1SumError2を最小とする最適化処理が実行される。この最適化処理の結果、図6に示すような、2値化補正DWIデータIb500cor2及び2値化変換DWIデータIb500doubleが生成される。また、500のb値に対応する2値化DWIデータIb500とゼロのb値に対応する2値化基準画像データIb0との間における位置ずれ及び500のb値に対応する2値化DWIデータIb500と500のb値に対応する2値化DWIデータIb500との間における実質的な位置ずれの双方を最小とする座標変換係数の組合せK12を得ることができる。そして得られた座標変換係数の組合せK12を500のb値に対応するDWIデータIb500に対する歪補正係数として設定することができる。更に、歪補正係数を用いて500のb値に対応するDWIデータIb500を座標変換することにより、図6に示すような歪補正後のDWIデータIb500corを得ることができる。尚、上述の例では、500のb値に対応する2値化DWIデータIb500の各座標変換係数をパラメータとして最適化処理を行う場合について説明したが、同様に500のb値に対応する2値化DWIデータIb500の各座標変換係数をパラメータとする最適化処理を行うことによって歪補正係数を求めることもできる。また、b値をゼロとした画像データを用いずに歪補正係数を求めることもできる。図7は、図2に示す歪補正部41Bにおいて歪補正係数を求めて歪補正を行う第3の方法を説明する図である。図7において各縦軸はx軸を示す。図5と同様に、図7に示すようなb値が±500に対応する歪方向が互いに逆の2つのDWIデータIb500Ib500を収集して2値化処理を行うことによって、図7に示すように、それぞれ500500のb値に対応する2値化画像データIb500Ib500を生成することができる。次に、500のb値に対応する2値化画像データIb500の各座標変換係数をパラメータとして、500のb値に対応する2値化DWIデータIb500と500のb値に対応する2値化画像データIb500との間における各画素値の2乗差の総和を最小とする最適化処理が実行される。この最適化処理の結果、図7に示すような、500のb値に対応する2値化画像データIb500の座標変換後における2値化変換DWIデータIb500doubleが生成される。また、500のb値に対応する2値化DWIデータIb500と500のb値に対応する2値化DWIデータIb500との間における実質的な位置ずれを最小とする座標変換係数の組合せK21を得ることができる。そうすると、2値化変換DWIデータIb500doubleの生成に用いられた各座標変換係数の組合せK21は、500のb値に対応する2値化DWIデータIb500の2値化基準画像データIb0に対する各画素値の2乗差の総和を最小にする最適な各座標変換係数の2倍に相当する。そこで、2値化変換DWIデータIb500doubleの生成に用いられた各座標変換係数の12の座標変換係数の組合せ2を500のb値に対応するDWIデータIb500に対する歪補正係数として設定することができる。更に、歪補正係数を用いて500のb値に対応するDWIデータIb500を座標変換することにより、図7に示すような歪補正後のDWIデータIb500corを得ることができる。尚、図5、図6及び図7に示す第1、第2及び第3の方法において、最適化処理の対象を画素値の2乗差の総和としたが、2値化画像データ間における位置ずれを表す相関係数等の他の指標を最適化処理の対象としてもよい。また、エッジ検出処理としての2値化処理を行わずにDWIデータIb500Ib500に対する座標変換係数をパラメータとして画像データ間の位置ずれの指標を最小化する最適化処理を実行するようにしてもよい。或いは、2値化処理を含む又は2値化処理を伴わない任意のエッジ検出処理によってDWIデータを含む複数の画像データの各エッジを検出し、検出したエッジ間の位置ずれを表す指標を最小とする最適化処理を行ってもよい。そして、2つのDWIデータに対するエッジ検出を含む処理によって歪補正係数を算出することができる。従って、図5に示す第1の方法は、2つのDWIデータ又は2つのDWIデータに基づいて生成される2つの画像データに対する各座標変換係数をパラメータとして、2つのDWIデータ又は2つの画像データの、b値をゼロとして収集された画像データに対する位置ずれの指標を最小にする最適化処理によって各座標変換係数を計算し、より小さい指標に対応する座標変換係数を歪補正係数とする方法であると言うことができる。また、図6に示す第2の方法は、2つのDWIデータの一方又は2つのDWIデータに基づいて生成される2つの画像データの一方に対する座標変換係数をパラメータとして、一方の、b値をゼロとして収集された画像データに対する位置ずれの指標及び他方に対する位置ずれの指標を最小にする最適化処理によって計算された座標変換係数に基づいて歪補正係数を決定する方法であると言うことができる。更に、図7に示す第3の方法は、2つのDWIデータの一方又は2つのDWIデータに基づいて生成される2つの画像データの一方に対する座標変換係数をパラメータとして、一方の他方に対する位置ずれの指標を最小にする最適化処理によって計算された座標変換係数に基づいて歪補正係数を決定する方法であると言うことができる。ここまでは、絶対値が同一で符号が互いに逆のb値に対応する2つのDWIデータを収集し、歪補正部41Bが、2つのDWIデータに基づいて2つのDWIデータの一方に対する歪補正係数を算出する例について説明したが、歪補正部41Bでは、補正対象となるDWIデータに対応するb値と絶対値が異なるb値に対応し、かつMPGパルスの印加軸が同一のDWIデータに基づいて歪補正係数を求めることもできる。以下、絶対値が異なるb値に対応する複数のDWIデータを収集し、歪補正部41Bが、基準となるb値に対応する歪補正係数に基づいて基準となるb値よりも大きいb値に対応する歪補正係数を算出する例について説明する。具体的には、b値と各歪補正係数Kとの関係式Kfに基づいて、各b値に対応する歪補正係数Kを求めることができる。関係式の係数は、係数の数に対応する数の基準となるb値及び基準となるb値に対応する歪補正係数Krefから求めることができる。b値と各歪補正係数Kとの関係式は、高次式、指数関数、対数関数等の任意の関数で近似することができる。近似に用いる適切な関数の種類及び高次式で近似する場合に2次以上の各項がそれぞれ無視できるか否かは、MPGパルスの振幅を制御することによってb値を調整するのかMPGパルスの印加時間を制御することによってb値を調整するのか或いはMPGパルスの振幅及び印加時間の双方を制御することによってb値を調整するのかといったようなMPGパルスの設計方法に依存する。そこで、MPGパルスの条件を用いたシミュレーション、装置の据付時又はMPGパルスの設計時におけるファントムを用いた試験撮像或いはMPGパルスの印加を伴うプレスキャンによって収集された画像データの画像処理によって予めb値と各歪補正係数Kとの関係を表す近似式を決定しておくことができる。b値と各歪補正係数が正比例するとみなせる場合には、式に示すように、1組の基準となるb値に対応する歪補正係数Krefに基づいて各b値に対応する歪補正係数Kを求めることができる。Kb基準となるb値は、ゼロでなく、かつDWIの撮像条件として設定され得る小さい値とすることが、歪の小さいDWIに対応する歪補正係数Krefを基準として、各b値に対応する歪補正係数Kを精度良く求める観点から好適である。すなわち、より小さなb値に対応する歪補正係数に基づいて、より大きなb値に対応する歪補正係数を推定することができる。例えば、b0の画像データ及びb500b1000のDWIデータを収集する場合には、b500として収集されたDWIデータに対応する歪補正係数に基づいて、式によりb1000として収集されたDWIデータに対応する歪補正係数を計算することができる。基準となるb値に対応する歪補正係数Krefは、図5、図6及び図7を参照して上述したように符号の異なるb値に対応する2つのDWIデータIbrefIbref又は符号の異なるb値に対応する2つのDWIデータIbrefIbrefとb値がゼロの基準画像データIb0に基づいて求めることができる。或いは、正又は負の単一の符号のb値に対応するDWIデータIbrefとb値がゼロの基準画像データIb0に基づいて基準となるb値に対応する歪補正係数Krefを求めるようにしてもよい。この場合、DWIデータIbrefを2値化した2値化DWIデータIbrefの座標変換係数をパラメータとして、2値化DWIデータIbrefの2値化基準画像データIb0に対する各画素値の2乗差の総和を最小にする最適化計算によって座標変換係数として基準となるb値に対応する歪補正係数Krefを求めることができる。このようにして歪補正部41Bにおいて算出された歪補正係数は、歪係数記憶部44に記憶させることができる。すなわち、MPGパルスの印加軸及びb値ごとに対応する歪補正係数を歪係数記憶部44に保存し、歪補正部41BがDWIの撮像条件に対応する歪補正係数を歪係数記憶部44から読み込んで取得することができる。また、b値から歪補正係数を求める関数自体又は関数の係数を歪係数記憶部44に保存し、b値に応じて歪補正係数を計算できるようにしてもよい。このため、ファントムを用いた試験撮像や過去の撮像において得られたDWIに基づいて歪補正係数を求め、求めた歪補正係数を歪係数記憶部44に保存すれば、歪補正部41Bが歪係数記憶部44を参照するのみで歪補正係数を取得することができる。この結果、各撮像の際に歪補正部41Bによる歪補正係数の計算処理を省略することができる。画像処理部41Cは、データ処理部41において生成された画像データ又は画像データ記憶部43から読み込んだ画像データに対して画像表示のための必要な画像処理を施す機能と、画像処理後における表示用の画像データを表示装置34に表示させる機能を有する。例えば、必要に応じて歪補正後における複数のDWIデータに第2のアベレージング処理を施すことができる。これにより、診断に必要なSNRを有する歪補正後のDWIデータを得ることができる。画像再構成部41Aにおいて実行される第1のアベレージング処理は、歪補正前におけるDWIデータに対して施される。従って、第1のアベレージング処理の対象は、歪量が同等となる、MPGパルスの印加軸及びb値が同一の複数のDWIデータを対象とすることが前提条件になる。これに対して、第2のアベレージング処理は、歪補正後における複数のDWIデータに対して施される。このため、第2のアベレージング処理は、互いに異なるb値に対応する複数のDWIデータに対して行うことができる。そこで、MR信号の絶対値を信号値とする画像再構成処理によってDWIデータが生成される場合には、絶対値が同一で極性が異なるb値等の互いに異なるb値に対応する複数のDWIデータを対象として第2のアベレージング処理を実行することが好適である。これにより、歪補正係数を求めるために追加的に収集されたDWIデータをSNRの向上のために利用することができる。すなわち、診断用に必要でないb値に対応するDWIデータが歪補正係数の算出のために収集された場合であっても、診断用のDWIデータに必要とされるSNRを得るために用いることができる。この結果、必要なSNRを得るために加算対象となる、診断用に必要なb値に対応するDWIデータの数を減らすことができる。一方、診断用に必要となるb値に対応するDWIデータが歪補正係数の算出のために収集された場合には、そのまま診断用のDWIデータの生成用に用いることができる。換言すれば、診断用に本来収集されるDWIデータを歪補正係数の算出用のDWIデータとして利用することができる。従って、上述したような歪補正部41Bにおける歪補正係数の算出に伴う実質的なデータ収集時間及びデータ収集量の増加を回避することができる。換言すれば、単位データ収集時間当たりのSNR向上率を一定にすることができる。第2のアベレージング処理は、MPGパルスの印加軸及びb値が同一の複数のDWIデータを対象として実行されるようにしてもよい。また、第1及び第2のアベレージング処理の一方又は双方を実行しないようにしても良い。従って、歪補正前及び前記歪補正後における所定の複数のDWIデータの少なくとも一方の平均をとるアベレージング処理を行うようにしてもよい。そして、このような構成を有する磁気共鳴イメージング装置20により、印加軸が同一で異なるb値に対応するMPGパルスの印加を伴ってDWIを実行し、異なるb値に対応する複数のDWIデータを収集することができる。また、b値をゼロとして画像データを収集することもできる。次に磁気共鳴イメージング装置20の動作及び作用について説明する。まず、磁気共鳴イメージング装置20により、歪補正係数の算出処理を伴ってDWIを行う場合について説明する。図8は、図1に示す磁気共鳴イメージング装置20により歪補正係数の算出処理を伴ってDWIを行う際の流れを示すフローチャートである。まずステップS1において、撮像条件設定部40は、入力装置33から入力された情報に基づいて図3に示すようにb値をゼロとしたEPIシーケンス及び診断目的に応じて必要なb値に対応するDWIシーケンスを設定する。更に、撮像条件設定部40は、最も小さいb値と絶対値が同じで符号が逆のb値に対応するDWIシーケンスを設定する。すなわち、図4に示すように印加方向が同一で極性が互いに逆のMPGパルスの印加を伴う2つのDWIシーケンスを少なくとも含む撮像条件が撮像条件設定部40において設定される。また、必要に応じて、DWIデータの収集及び加算の回数を決定するために必要となる第1及び第2のアベレージング処理のための各AVERAGE数が撮像条件として設定される。ここでは、例えば、b500及びb1000に対応する2つのDWIデータが診断に必要であるものとする。この場合には、b0b500b500b1000に対応するEPIシーケンス及び複数のDWIシーケンスが設定される。また、同一のb値を有するDWIデータを対象とする第1のアベレージング処理のAVERAGE数が2であれば、例えばb0に対応するEPIシーケンス、b500の1AVERAGE目に対応するDWIシーケンス、b500の1AVERAGE目に対応するDWIシーケンス、b1000の1AVERAGE目に対応するDWIシーケンス、b500の2AVERAGE目に対応するDWIシーケンス、b500の2AVERAGE目に対応するDWIシーケンス及びb1000の2AVERAGE目に対応するDWIシーケンスの順にシーケンスの順序が設定される。次にステップS2において、撮像条件設定部40において設定された撮像条件に従ってDWIスキャンを含むイメージングスキャンが実行される。そのために予め寝台37に被検体Pがセットされ、静磁場電源26により励磁された静磁場用磁石21の撮像領域に静磁場が形成される。また、シムコイル電源28からシムコイル22に電流が供給されて撮像領域に形成された静磁場が均一化される。そして、撮像条件設定部40から撮像条件として設定された複数のEPIシーケンスがシーケンスコントローラ31に順次与えられる。そうすると、シーケンスコントローラ31は、パルスシーケンスに従って傾斜磁場電源27、送信器29及び受信器30を駆動させることにより被検体Pがセットされた撮像領域に傾斜磁場を形成させるとともに、RFコイル24からRF信号を発生させる。このため、被検体Pの内部における核磁気共鳴により生じたMR信号が、RFコイル24により受信されて受信器30に与えられる。受信器30は、RFコイル24からMR信号を受けて、所要の信号処理を実行した後、AD変換することにより、デジタルデータのMR信号である生データを生成する。そして、受信器30は、生成したMR信号をシーケンスコントローラ31に与える。更にシーケンスコントローラ31は、MR信号をコンピュータ32のデータ処理部41に出力する。そうすると、データ処理部41は、MR信号をk空間データとしてk空間データ記憶部42に形成されたk空間に配置する。この結果、k空間データ記憶部42には、b0に対応するk空間データ、b500の1AVERAGE目に対応するk空間データ、b500の1AVERAGE目に対応するk空間データ、b1000の1AVERAGE目に対応するk空間データ、b500の2AVERAGE目に対応するk空間データ、b500の2AVERAGE目に対応するk空間データ及びb1000の2AVERAGE目に対応するk空間データが順次書き込まれて保存される。次にステップS3において、画像再構成部41Aは、k空間データ記憶部42からk空間データを取り込んで画像再構成処理を施すことにより画像データを再構成する。この結果、b0に対応する非DWIデータである基準画像データ、b500の1AVERAGE目及び2AVERAGE目に対応する2フレーム分のDWIデータ、b500の1AVERAGE目及び2AVERAGE目に対応する2フレーム分のDWIデータ、b1000の1AVERAGE目及び2AVERAGE目に対応する2フレーム分のDWIデータがそれぞれ生成される。また、画像再構成部41Aは、MPGパルスの印加軸及びb値を同一に設定して収集された複数のDWIデータを加算する第1のアベレージング処理によってAVERAGE画像データを生成する。すなわち、b500の1AVERAGE目及び2AVERAGE目に対応する2フレーム分のDWIデータ同士の加算、b500の1AVERAGE目及び2AVERAGE目に対応する2フレーム分のDWIデータ同士の加算並びにb1000の1AVERAGE目及び2AVERAGE目に対応する2フレーム分のDWIデータ同士の加算が画像再構成部41Aにおいて行われる。この結果、b500b500及びb1000にそれぞれ対応する3フレーム分のDWIAVERAGEデータが生成される。次にステップS4において、歪補正部41Bは、b値が小さく、かつ符号が逆の2フレーム分のDWIAVERAGEデータに基づいて小さいb値に対応する歪補正係数を算出する。すなわち、歪補正部41Bは、b±500に対応する2フレーム分のDWIデータに基づいてb500に対応する歪補正係数を算出する。歪補正係数は、図5、図6又は図7に示すように、2値化処理後のb500に対応する2値化DWIAVERAGEデータ及び2値化処理後のb500に対応する2値化DWIAVERAGEデータの一方又は双方に対する座標変換係数をパラメータとする最適化処理によって計算することができる。図5又は図6で示される方法で歪補正係数が算出される場合には、2値化処理後のb0に対応する2値化基準画像データも用いられる。そして、図5に示す方法で歪補正係数が算出される場合には、最適化処理において最小とする対象が2値化基準画像データに対するb500及びb500に対応する2値化DWIAVERAGEデータの各画素値の2乗差の総和に設定される。また、図6に示す方法で歪補正係数が算出される場合には、最適化処理において最小とする対象が2値化基準画像データ及びb500に対応する2値化DWIAVERAGEデータに対するb500に対応する2値化DWIAVERAGEデータの各画素値の2乗差の総和に設定される。また、図7に示す方法で歪補正係数が算出される場合には、最適化処理において最小とする対象がb500に対応する2値化DWIAVERAGEデータに対するb500に対応する2値化DWIAVERAGEデータの各画素値の2乗差の総和に設定される。このような最適化処理により算出される歪補正係数は、歪方向が互いに逆のb±500に対応するDWIデータに基づいて算出されているため、2値化基準画像データに対するb500に対応する2値化DWIAVERAGEデータの各画素値の2乗差の総和を最小とする最適化処理によって算出される歪補正係数に比べて高精度となる。次にステップS5において、歪補正部41Bは、小さいb値に対応する歪補正係数に基づいて大きいb値に対応する歪補正係数を求める。すなわち、歪補正部41Bは、b値と歪補正係数が比例するものとして、式によりb500に対応する歪補正係数からb1000に対応する歪補正係数を算出する。これにより、撮像条件として設定された全てのb値に対応する歪補正係数を簡易に取得することができる。取得したb値ごとの歪補正係数は、必要に応じて歪係数記憶部44に書き込んで保存することができる。次にステップS6において、歪補正部41Bは、各b値に対応する歪補正係数を用いてそれぞれDWIAVERAGEデータの歪補正を行う。すなわち、b500及びb1000にそれぞれ対応する各歪補正係数を用いてb500及びb1000に対応するDWIAVERAGEデータの歪が補正される。次にステップS7において、画像処理部41Cは、b0に対応する画像データ及び歪補正後のDWIAVERAGEデータに対して、画像表示のための必要な画像処理を施す。例えば、異なるb値に対応する歪補正後の複数のDWIAVERAGEデータ間において第2のアベレージング処理を施すことができる。本例では、b500及びb1000に対応する2つのDWIデータが診断に必要であり、b500に対応するDWIAVERAGEデータは診断には不要である。そこで、b±500に対応するDWIAVERAGEデータに対して第2のアベレージング処理を実行することができる。これにより、SNRを向上させたb500に対応する診断用のDWIAVERAGEデータを生成することができる。次にステップS8において、画像処理部41Cは、画像処理後における表示用のb0に対応する画像データ及びDWIAVERAGEデータを表示装置34に表示させる。これにより、ユーザは、表示装置34に表示されたb0に対応する画像、b500に対応するDWI及びb1000に対応するDWIを確認することができる。このとき表示装置34に表示されるDWIは、精度良く算出された歪補正係数に基づいて歪補正されているため、良好に歪が補正された画像となる。加えて、表示装置34に表示されるDWIは、第1及び第2のアベレージング処理によってSNRが向上した画像となる。次に、磁気共鳴イメージング装置20により、予め取得した歪補正係数を用いた歪補正処理を伴ってDWIを行う場合について説明する。図9は、図1に示す磁気共鳴イメージング装置20により予め取得した歪補正係数を用いた歪補正処理を伴ってDWIを行う際の流れを示すフローチャートである。尚、図8に示すフローチャートのステップと同様なステップには同符号を付して説明を省略する。まず予めファントムを用いた試験撮像や過去の撮像において求められた歪補正係数が歪係数記憶部44に保存される。但し、歪補正係数は、図5、図6又は図7に示すように、絶対値が同一でかつ符号が互いに逆のb値に対応する2フレーム分のDWIAVERAGEデータに基づく座標変換係数をパラメータとする最適化処理によって計算される。また、歪補正係数の算出には、DWIの撮像条件として設定され得る小さいb値に対応するDWIAVERAGEデータが用いられる。例えば、b500に対応する歪補正係数が歪係数記憶部44に保存されるものとする。次に、ステップS11において、撮像条件設定部40は、入力装置33から入力された情報に基づいてb値をゼロとしたEPIシーケンス及び診断目的に応じて必要なb値に対応するDWIシーケンスを設定する。また、必要に応じてDWIのAVERAGE数が設定される。尚、図8のステップS1とは異なり、あるb値と絶対値が同じで符号が逆のb値に対応するDWIシーケンスを設定する必要はない。このため、例えば、b0b500及びb1000に対応するEPIシーケンスが設定されるものとする。次に、ステップS12において、撮像条件設定部40において設定された撮像条件に従ってDWIスキャンを含むイメージングスキャンが実行される。すなわち、図8のステップS1と同様な流れで、絶対値が互いに異なる複数のb値にそれぞれ対応するk空間データが収集される。次に、ステップS13において、画像再構成部41Aは、絶対値が互いに異なる複数のb値にそれぞれ対応する画像データ及びDWIデータを生成する。そのために、k空間データに対する画像再構成処理が実行される他、撮像条件として設定されたAVERAGE数に応じてDWIデータに対する第1のアベレージング処理が実行される。第1のアベレージング処理が実行された場合には、DWIAVERAGEデータが生成される。次に、ステップS14において、歪補正部41Bは、歪係数記憶部44からb値が小さく、かつ符号が逆の2フレーム分のDWIAVERAGEデータ又はDWIデータに基づいて算出された小さいb値に対応する歪補正係数を取得する。すなわち、歪補正部41Bは、b500に対応する歪補正係数を歪係数記憶部44から取得する。これにより、ステップS5において、歪補正部41Bは、b500に対応する歪補正係数に基づいてb1000に対応する歪補正係数を算出することが可能となる。そして図8のステップS6からステップS8と同様に歪補正係数を用いた歪補正後のDWIに必要な画像処理を施した画像を表示装置34に表示させることができる。つまり以上のような磁気共鳴イメージング装置20は、印加軸及びb値の絶対値が同一で極性が互いに逆のMPGパルスを印加して拡散イメージングすることによって歪み極性が逆の2つのDWIデータを収集し、歪み極性が逆の2つのDWIデータに基づいて精度良く歪補正係数を求めるようにしたものである。このため、磁気共鳴イメージング装置20によれば、従来よりも精度よく歪み補正が施されたDWIを得ることができる。特に、被検体Pの腹部には限られた方向にのみ拡散する神経の束がない。このため、撮像部位を腹部としてDWIを収集する場合には、MPGパルスの印加軸を1軸方向として、より精度良く歪補正係数を求めて良好に歪が補正されたDWIを得ることができる。一方、MPGパルスの印加軸が複数方向であっても、印加軸ごとに互いに極性の異なるMPGパルスを印加することによって、より精度良く歪補正係数を計算することができる。また、磁気共鳴イメージング装置20によれば、あるb値に対応する歪補正係数に基づいて、他のb値に対応する歪補正係数を計算することができる。このため、歪補正係数の算出に必要なデータ処理量を低減させることができる。以上、特定の実施形態について記載したが、記載された実施形態は一例に過ぎず、発明の範囲を限定するものではない。ここに記載された新規な方法及び装置は、様々な他の様式で具現化することができる。また、ここに記載された方法及び装置の様式において、発明の要旨から逸脱しない範囲で、種々の省略、置換及び変更を行うことができる。添付された請求の範囲及びその均等物は、発明の範囲及び要旨に包含されているものとして、そのような種々の様式及び変形例を含んでいる。"