磁気共鳴イメージング装置、及び、位相アンラップ方法
【課題】エコープレナー法で計測したエコー信号を前計測で得た補正データで位相補正する際に、信号強度が弱い領域やデータ点があることによって生じる位相乱れを排除した補正データを取得し、これによりストリークアーチファクトのない画像を得ること。
【解決手段】エコー信号を奇数エコーと偶数エコーに分けて、それぞれエコー間の位相差分に対し周波数方向にフィッティングし、補正値となる位相変化量を求める。この位相変化量をエコートレイン方向の累積する際に基準となる奇数エコー及び偶数エコーについても、両者間の位相差分を求め、フィッティングする。フィッティングによって求めたフィッティングデータは、信号値が極端に弱い領域や点に限定的に用い、十分な信号強度が得られる領域は実際の位相値を用い、フィッティング誤差の発生を防止する。
【解決手段】エコー信号を奇数エコーと偶数エコーに分けて、それぞれエコー間の位相差分に対し周波数方向にフィッティングし、補正値となる位相変化量を求める。この位相変化量をエコートレイン方向の累積する際に基準となる奇数エコー及び偶数エコーについても、両者間の位相差分を求め、フィッティングする。フィッティングによって求めたフィッティングデータは、信号値が極端に弱い領域や点に限定的に用い、十分な信号強度が得られる領域は実際の位相値を用い、フィッティング誤差の発生を防止する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、磁気共鳴イメージング装置（以下、ＭＲＩ装置という）におけるエコープラナー（ＥＰＩ）計測時に生じるストリークアーチファクトを抑制する技術に関する。
【背景技術】
【０００２】
エコープラナー計測では、エコー信号を取得するときの読出し傾斜磁場の極性が異なることにより、エコートレインのうち奇数番目に取得されるエコー（以下、奇数エコーという）と偶数番目に取得されるエコー（以下、偶数エコーという）とにピークずれを生じる。このピークずれは、エコープラナー計測で得た画像上では、Ｎ/２アーチファクトとなって現れるが、前計測データを用いて奇数エコー及び／又は偶数エコーの位相を補正することにより、抑制することができる。
【０００３】
ストリークアーチファクトは、主として、位相エンコード方向の位相変化や二次的に混入する位相ノイズによって引き起こされる。特許文献１には、Ｎ/２アーチファクトの補正と同時に、ストリークアーチファクトを抑制する技術が開示されている。
【０００４】
特許文献１に記載されたストリークアーチファクト抑制方法では、位相エンコード傾斜磁場パルスを印加しない撮像で計測した前計測エコーを、奇数エコーと偶数エコーとに分離した後、奇・偶数エコーそれぞれのエコートレイン方向の位相変化を滑らかに接続することで、前計測データに含まれる位相乱れを抑制する。そして、この位相乱れが補正された前計測データを用いて本計測データを補正することにより、ストリークアーチファクトを低減する。
【０００５】
ＭＲＩで取得される信号や画像データの位相を滑らかに接続する処理には、一般にフィッティングが用いられるが、その前提技術として、隣接する点の位相が２πを超えた場合に、その整数倍（２π×Ｎ（Ｎ＝０，１，２・・・））を一方に加算又は減算する処理即ち主値回りした位相回転量を戻す処理（位相アンラップ処理と呼ばれる）があり、ＭＲＩの画像処理における必須技術となっている。位相アンラップ処理では、信号値が低いデータでは位相が乱れやすく信頼性のある処理が行えないという問題がある。
【０００６】
この問題に対し、特許文献２には、二次元画像データに対する位相アンラップ処理として、データを閾値処理し、所定の閾値以上のデータの領域（島領域）のデータについてのみ位相アンラップ処理し、島領域と島領域との間については、隣接する領域の位相との整合が取れるように所望の値を選択する処理が提案されている。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２００９−２７３５３０号公報
【特許文献２】国際特許公報ＷＯ０２/０８２９９４
【発明の概要】
【発明が解決しようとする課題】
【０００８】
特許文献１に記載された技術では、奇数エコー間又は偶数エコー間の位相差分及び奇・偶数エコー間の位相差分に対し、周波数方向の位相変化を直線でフィッティングし、それを位相エンコード方向に位相累積することにより、エコートレイン方向の位相変化を滑らかに接続している。
【０００９】
この技術では、周波数方向のフィッティングに単純な一次フィッティングを採用しているが、フィッティング処理には少なからずフィッティング誤差が生じ、Ｎ/２アーチファクトの原因となりうる。また奇・偶数エコーそれぞれについて、位相差分データを位相エンコード方向に位相累積する際には、開始点となるエコーを基準エコーとして、順次位相累積が行われるが、これら基準エコーの信号強度が極端に低く、両者の位相差に大きな乱れが生じる場合には、ストリークアーチファクトが生じてしまう可能性がある。
【００１０】
特許文献１の技術では、エコートレイン方向の中で最も信号強度が高いエコーを基準エコーとして採用することで極力高い信号強度が得られるような工夫がなされているが、周波数方向に対しては十分な信号強度が得られない場合がある。例えば、マルチプルアレイコイルのように空間的に不均一な感度分布を持った受信コイルを用いて撮像する場合には、局所的に信号強度が近い場所が存在するため、その領域における基準エコー間の位相差に大きな乱れが生じる可能性があり、それによりストリークアーチファクトが発生する。
【００１１】
一方、特許文献２に記載された技術は、信号強度が低く、位相情報の精度が低いデータを位相アンラップ処理から除外するものであり、精度の高いアンラップ処理を行うことができるが、位相アンラップ処理する際に、２π×ＮのＮを検索するステップや加算又は減算かを判定するステップ等が必須であり、位相アンラップ処理されない領域（島領域と島領域との間）については、評価関数を用いた評価を必要とし、複雑な処理を要する。
【課題を解決するための手段】
【００１２】
上記課題を解決するため、本発明のＭＲＩ装置は、エコープラナー計測によって取得したエコーデータを補正するための補正用データを、前計測データの奇エコーと偶数エコーとの位相差情報をフィッティンングして求める際に、フィッティングを信号強度の比較的低いデータ領域に限定的に行い、フィッティイング誤差の発生を抑制する。
【００１３】
また本発明のＭＲＩ装置は、奇数エコー及び偶数エコーそれぞれを位相エンコード方向に位相累積する際の基準エコーについて、奇・偶数エコー間の位相差についても上記限定されたフィッティングを行い、基準エコー間の位相乱れに起因するアーチファクトの発生を抑制する。
【００１４】
本発明の位相アンラップ処理方法は、隣り合う点の位相差を単位ベクトルに変換し、位相差として再計算する。これにより、もとの位相差に±２π×Ｎ（Ｎ＝０，１，２・・・）を加算して最小値を求める計算と同等な処理を簡易に行うことができる。
【００１５】
即ち、本発明のＭＲＩ装置は、検査対象の核磁気共鳴により発生するエコー信号を取得する撮像部と、前記撮像部が取得したエコー信号を処理し、核磁気共鳴画像を作成する画像処理部とを備え、前記画像処理部が処理するエコー信号は、第１の計測（前計測）で取得した第１計測データと、第２の計測(本計測)で取得した第２計測データとを含み、且つ第１計測データ及び第２計測データは、それぞれ、異なる極性の読出し傾斜磁場を用いて取得した第１群のエコーデータ（奇数エコーデータ）と第２群のエコーデータ（偶数エコーデータ）を含み、前記画像処理部は、前記第１計測データを用いて前記第２計測データを補正する補正手段を備え、前記補正手段は、前記第１計測データの第１群のエコーデータと第２群のエコーデータとの位相差分データを周波数方向に連続する位相差分データに変換して補正用データを作成し、その際、信号強度が閾値以下の周波数方向のデータ点または領域のみを予め求めたフィッティングデータで置換する。
【００１６】
本発明のＭＲＩ装置において、前記補正手段は、第１群及び第２群のエコーデータのうち、閾値以上のデータ点を選択し、選択されなかったデータ点または領域によって隔てられたデータ点の領域（島領域）毎に、位相差分を周波数方向に連続にする処理を行い、隣接する島領域と不連続な点或いは領域によって隔てられた島領域について、当該島領域の端部の位相値を、前記隣接する島領域から外挿によって予測される目標値と実際の位相値との差異が最小にするように決定する。
【００１７】
また本発明の位相アンラップ方法は、位相の主値回りを生じているデータに対し、主値回りを除去する方法であり、アンラップ方向の先端にあるデータ点の位相値（第１位相値）とそれに隣接するデータ点の位相値（第２位相値）の位相差が±２πを超えるときに、前記第１位相値と前記第２位相値との位相差を複素数に変換するステップと、前記複素数の実部と虚部の値から位相値を算出するステップと、前記算出された位相値を前記第１位相値に加算した値を、第２位相値に置き換えるステップと、を含む。
【発明の効果】
【００１８】
本発明のＭＲＩ装置によれば、ＥＰＩデータの補正用データを前計測によって得る際に、前計測エコーデータの周波数方向の位相フィッティングを信号値の弱い領域等に限定することにより、フィッティング誤差によるストリークアーチファクトが低減されたＭＲ画像を得ることができる。
【００１９】
また本発明のＭＲＩ装置によれば、奇数エコーと偶数エコーの基準となるエコー間の位相差に対してもフィッティング処理を加えることで、従来よりもアーチファクトを高精度に低減した画像を得ることができる。
【００２０】
さらに本発明によれば、ＭＲＩの基本技術である位相アンラップ方法として、精度が高く且つ従来の位相アンラップ方法に比べ簡便は位相アンラップ方法が提供される。
【図面の簡単な説明】
【００２１】
【図１】本発明が適用されるＭＲＩ装置の全体概要を示すブロック図
【図２】ＭＲＩ装置の動作の一実施形態を示すフローチャート
【図３】本計測のパルスシーケンスの一例を示す図
【図４】前計測のパルスシーケンスの一例を示す
【図５】図２のフローチャートのステップ230の詳細を示すフローチャート
【図６】実施形態における閾値処理の説明図で、（ａ）、（ｂ）はそれぞれ閾値処理前のエコーデータの振幅データと位相データを示し、（ｃ）、（ｄ）はそれぞれ閾値処理後のエコーデータの振幅データと位相データを示す。
【図７】実施形態における第１位相アンラップ処理の説明図で、（ａ）は差分後の位相データを示し、（ｂ）は位相アンラップ処理後の位相データを示す。
【図８】実施形態における最大島探索処理の説明図で、（ａ）は信号強度によって分けられた島の振幅データを示し、（ｂ）は各島の位相データを示す。
【図９】実施形態における最大島に対するフィッティング処理の説明図で、（ａ）は最大島の位相データ、（ｂ）は最大島の位相データを用いたフィッティング結果を示す図。
【図１０】実施形態における第２位相アンラップ処理の説明図で、（ａ）、（ｂ）は、それぞれ、目標値と実際の位相から求めた位相値に実際の位相を置換する様子を示す図、（ｃ）は第２位相アンラップ処理後の位相を示す図。
【図１１】実施形態におけるフィッティング処理の説明図で、（ａ）は第２位相アンラップ処理後の位相データを示す図、（ｂ）は（ａ）の位相データのうち、位相がゼロの領域或いは点のみをフィッティングデータで置換或いは線形補間する場合を示す図、（ｃ）は位相データ全体をフィッティングデータで置換する場合を示す図。
【図１２】実施形態による位相補正の効果の実例を示す図で、（ａ）は従来法による補正を行った場合、（ｂ）は本発明による補正を行った場合を示す。
【発明を実施するための形態】
【００２２】
以下、添付図面を参照して本発明のＭＲＩ装置の好ましい実施形態について説明する。
【００２３】
図１は、本発明を適用するＭＲＩ装置の全体構成を示すブロック図である。このＭＲＩ装置は、静磁場発生系１１、傾斜磁場発生系１２、送信系１３、受信系１４、シーケンサ１５、再構成演算部（画像処理部）１６、制御系１７および表示部（ディスプレイ）１８を備え、また図示しない操作部を備えている。主として、静磁場発生系１１、傾斜磁場発生系１２、送信系１３、受信系１４及びシーケンサ１５が、本発明における撮像部の機能を担っている。
【００２４】
静磁場発生系１１は、永久磁石、常電導磁石または超電導磁石のいずれかからなり、被検体１０が置かれる空間に均一な静磁場を発生させる。静磁場の方向は、通常、被検体１０の体軸方向または被検体１０の体軸と直交する方向である。
【００２５】
傾斜磁場発生系１２は、Ｘ、Ｙ、Ｚの３軸方向に巻かれた傾斜磁場コイル１２１と、これらの各々のコイルを磁化させる傾斜磁場電源１２２とからなり、シーケンサ１５からの命令に従って傾斜磁場電源１２２の各々のコイルを磁化させることによりＸ、Ｙ、Ｚの３軸方向の傾斜磁場Ｇｓ、Ｇｐ、Ｇｆを被検体１０に印加する。この傾斜磁場の加え方により、被検体１０の撮影断面が設定される。
【００２６】
送信系１３は、高周波発振器１３１、変調器１３２、高周波増幅器１３３及び高周波照射コイル１３４とからなり、傾斜磁場発生系１２で設定された被検体１０の断面を構成する原子の原子核に核磁気共鳴を起こさせるために、高周波発振器１３１から出力されたＲＦパルスを高周波増幅器１３３で増幅した後に、被検体１０に近接して設置された高周波照射コイル１３４に供給して被検体１０に照射する。
【００２７】
受信系１４は、高周波受信コイル１４１、受信回路１４２及びアナログ／ディジタル（以下「Ａ／Ｄ」という）変換器１４３からなり、送信系１３の高周波照射コイル１３４から照射された電磁波による被検体１０の原子核の核磁気共鳴によるＮＭＲ信号であるエコー信号を、被検体１０に近接して配置された高周波受信コイル１４１で検出し、受信回路１４２を介してＡ／Ｄ変換器１４３に入力し、ディジタル信号に変換して、さらにシーケンサ１５からの命令によるタイミングでサンプリングされたデータ（以下「サンプリングデータ」という）として、再構成演算部１６に送る。
【００２８】
制御系１７は、シーケンサ１５に指令を送り、シーケンサに設定された所定のパルスシーケンスに従い、上述した傾斜磁場発生系１２、送信系１３および受信系１４を制御するとともに、再構成演算部１６の動作を制御する。パルスシーケンスは、撮像方法によって種々のものがあり、それらはプログラム化されて制御系１７に格納されている。本実施形態のＭＲＩ装置では、ＥＰＩ系のパルスシーケンスを備えている。
【００２９】
再構成演算部１６は、サンプリングデータにフーリエ変換などの各種演算を行なうとともに、後述する位相変化や補正用データの算出、補正用データを用いた補正などを行なう。再構成演算部１６で再構成された画像は、表示部（ディスプレイ）１８に表示される。
【００３０】
次に上記構成におけるＭＲＩ装置の動作を説明する。図２に、動作の手順を示す。
本実施の形態のＭＲＩ装置では、画像化するエコー信号を計測する本計測（ステップ220）と、本計測とは別に実行される前計測（ステップ210）を実行する。ここでは、ステップ210の計測を「前計測」と呼んでいるが、前計測を実行するタイミングは、本計測の前と後のいずれでもよい。
【００３１】
本計測では、マルチエコー撮像によるパルスシーケンスを実行し、本計測によって得られたエコー信号（本計測エコー）３０を、前計測によって得られたエコー信号（前計測エコー）２０を用いて補正し（ステップ222）、画像化する（ステップ223）。本計測エコー３０の補正は、本計測エコー３０および前計測エコー２０をそれぞれ周波数方向にフーリエ変換（ステップ211、221）し、ｘ−ｋｙデータとし、差分することにより行なう。
【００３２】
前計測エコー２０を補正用データとして補正に供するに際し、奇数エコー２１と偶数エコー２２に分けて、それぞれの差分（位相差分）を周波数方向フィッティングし、エコートレイン方向に位相累積することにより、前計測エコーに２次的に混入する位相ノイズを除去するとともに時間変動を含む位相変化を反映した補正用データを得る。以下、各ステップの詳細を説明する。
【００３３】
＜計測ステップ＞
本計測（ステップ220）及び前計測（ステップ210）では、シングルショット又はマルチショットＥＰＩシーケンスにより、1回の励起で複数のエコー信号を計測する。図３及び図４に、本計測および前計測で用いられるパルスシーケンスの一例を示す。これらパルスシーケンスは公知のSpin Echo系ＥＰＩシーケンスであり、簡単に説明する。
【００３４】
まず検知する磁化を含む被検体にＲＦパルス301を照射すると同時にスライスを選択する傾斜磁場パルス302を印加し、画像化するスライスを選択後、前記傾斜磁場パルスによる位相分散を収束させるためのリフォーカス傾斜磁場パルス303を印加する。次いで、スピンエコーを生成するために１８０度ＲＦパルス305と２番目のスライス選択用の傾斜磁場パルス304、およびその前後に磁化の位相を分散させるためのクラッシャー傾斜磁場パルス（304、307）を印加する。
【００３５】
その後、周波数エンコード傾斜磁場のオフセットを与えるパルス309を印加した後、連続して反転する周波数エンコード傾斜磁場パルス311を印加し、反転する周波数エンコード傾斜磁場311の各周期内でエコー信号313が時系列的に発生するので、これをA/Dの時間範囲(312)の間おのおのサンプリングし時系列データを得る。反転する周波数エンコード傾斜磁場311の強度、印加タイミング、A/Dの時間範囲(312)は、本計測シーケンスと前計測と同じである。しかし、図３に示す本計測シーケンスでは、各エコー信号313に位相エンコードするための傾斜磁場パルス308、310が印加されるのに対し、図４に示す前計測シーケンスでは、位相エンコードは用いない。
【００３６】
シングルショットシーケンスの場合には、図示する繰り返し単位314を１回で、マルチショットシーケンスであれば、繰り返し単位314を複数回繰り返して、それぞれｋ空間全体を埋めるデータを取得する。1回の繰り返し単位で取得する複数のエコーをエコートレインと言い、奇数番目に取得されるエコーを奇数エコー、偶数番目に取得されるエコーを偶数エコーという。
【００３７】
このような計測をスライス傾斜磁場（302,306）を異ならせて行い、複数のスライスのデータを取得する。スライス数は特に限定されるものではないが、本実施の形態では、少なくとも前計測において、複数のスライスのデータを取得することが望ましい。複数のスライスのデータを利用することにより、より精度の高い補正用データを作成することができる。
【００３８】
このようなＥＰＩを基本とする計測シーケンスでは、周波数エンコード傾斜磁場の印加タイミングおよびサンプリング時間の誤差があると、エコーのピークがｋ空間中央（ｋｙ軸）に対してずれる。ずれる方向が、奇数エコーと偶数エコーとでは反対になるため、Ｎ/２アーチファクトを生じる。また位相エンコード方向のオフセットに誤差があると、エコーのピークのｋｙ軸からのずれが、エコー番号が大きくなるほど拡大する。以下の処理では、前計測エコーから、エコーのピークのずれに対応する位相情報（補正用データ）を得て、本計測エコーのピークずれを解消する。
【００３９】
＜周波数方向フーリエ変換（ステップ211、221）＞
本計測および前計測で得られたエコー信号２０、３０をそれぞれ周波数方向にフーリエ変換し、ｘ−ｋｙデータとする。周波数方向のフーリエ変換は、計測空間（ｋｘ−ｋｙ空間）におけるエコー信号のピークずれを、ｘ−ｋｙ空間上で位相補正するために行われる。即ち、計測空間におけるエコー信号のピークずれは、フーリエ変換後の空間（ｘ−ｋｙ空間）での位相変化に対応するので、前計測のｘ−ｋｙデータから求めた位相変化を用いて、ｘ−ｋｙ空間上で本計測データの位相を補正することにより、本計測データのピークずれが補正される。
【００４０】
＜補正用データの作成＞
前計測エコーをフーリエ変換した後のｘ−ｋｙデータ（以下、エコーデータという）を、奇数エコーと偶数エコーに分けて、それぞれについて補正用データを得る。以下、奇数エコーと偶数エコー各処理の詳細を説明する。
【００４１】
＜奇数エコーデータの処理＞
奇数エコー間の位相差分を求める。奇数エコー間或いは偶数エコー間で位相差分した場合、周波数エンコード傾斜磁場の印加タイミングとサンプリング時間の誤差はキャンセルされるので、位相エンコード方向のオフセット誤差のみが位相差分として求められる。この位相差分は、周波数方向に直線となるものであり、位相差分後の前計測データをフィッティングすることにより求められる。
【００４２】
具体的には、エコートレインを構成する複数のエコーのうち、２つのエコーデータ、例えばエコートレインの１番目のエコーと３番目のエコーとの位相差分を取る。位相差分処理を行うエコーの一方は、エコートレインのうち最も信号強度が高いＴＥの位置のエコーを用いることが好ましい。このようなエコーは、後述する位相累積処理を開始する基準エコーとなる。
【００４３】
＜＜閾値処理（ステップ212）＞＞
２つの奇数エコー間で位相差分する際に、信号強度が低く位相値の信頼性が低いデータをフィッティングに用いないために、エコーデータから信号強度が低いデータを除外する。除外するデータは、所定の閾値以下、例えば対象とするエコーデータの最大振幅値に対して所定の比率（例えば２０％）以下のデータとする。閾値Thresholdは次式（１）で表わされる。
（数１）
Threshold=Signalmax×０．２ （１）
ここで、Signalmaxはエコーデータ内の最大振幅値を意味する。
【００４４】
＜＜周波数方向のフィッティング（ステップ213）＞＞
次に２つのエコーデータを位相差分したデータに対し、周波数方向のフィッティングを行う。フィッティングは、最小二乗法等の公知の手法を用いることができ、式（２）で表わされる一次関数と実際の各データ点における位相差との差が最小となるように、式（３）及び（４）により、式（２）の係数ａ、ｂを求める。
【００４５】
【数２】
ここで、Difｏ(x)は奇数エコー間の位相差分をフィッティングした一次関数であり、Xiは周波数方向の座標、Ｙiは位相値で、nは周波数点数（データ点の点数）を意味する（以下、同様）。
【００４６】
なお、フィッティングに先立ち、位相の主値回りを除去するための位相アンラップ処理を行う。具体的には、隣合うデータ点の位相差を順次検索し、一つのデータ点の位相と次のデータ点の位相との差が２πを超えたときに、次のデータ点の位相に±２π×Ｎを加算する処理を行う。位相アンラップ処理は、一般的な手法を用いてもよいし、後述する簡易アンラップ処理を行ってもよい。
【００４７】
奇数エコー間の位相差分は、２つのエコーから求めることができるが、２つのエコーのみを用いるのではなく、複数の奇数エコー間で式（２）の係数を求めて、これを平均することが望ましい。
またマルチスライス計測の場合には、スライス毎に奇数エコー間の位相差分を求めて、これを平均してもよいし、特許文献１に記載されるように、スライスの信号を合成した後、エコー間の位相差分を求めてもよい。
【００４８】
＜＜位相累積（ステップ217）＞＞
２つのエコーデータ間の位相差は、位相エンコード方向のオフセット誤差に対応し、位相エンコード方向のオフセット誤差はエコートレイン方向に蓄積されるものである。そこで、上記フィッティングステップで求めた位相差を、エコートレイン方向に位相累積する。具体的には、位相累積を開始する基準となるエコーをエコートレインの（２×Ｂ+１）番目のエコーであるとすると、位相累積後の（２×ｉ+１）番目のエコーの位相変化Ｐは、次式（５）で表わされる。
【００４９】
（数３）
P(x,2×i+1)=P(x,2×B+1)+ Difｏ(x)×(i-B) （５）
【００５０】
以上の各ステップ211〜213、217により、エコートレインを構成する各奇数エコーについて、ピークを揃えるための補正用データが得られたことになる。
【００５１】
＜偶数エコーデータの処理＞
偶数エコーデータ２２については、周波数エンコード傾斜磁場の印加タイミングおよびサンプリング時間の誤差と、位相エンコード方向のオフセット誤差とを同時に補正するために、まず奇数エコーとの位相差分を行い（ステップ214）、位相差分データについて周波数方向のフィッティングを行う（ステップ230）。奇数エコーと偶数エコーとの位相差分も、理論的には一次直線でフィッティングできるものであるが、奇数エコー間或いは偶数エコー間の位相差分が小さい値の範囲に留まるのに対し、スケール（値が変化する範囲）が大きくなるため、適切な位相アンラップ処理を行えない可能性があり、単純な一次フィッティングを行った場合、フィッティング誤差が生じる可能性が高まる。
【００５２】
そこで偶数エコーデータ（奇数エコーとの位相差分データ）の処理ステップ230では、奇数エコー間の処理とは異なる位相アンラップ処理と限定されたフィッティングを行う。以下、ステップ230の処理の詳細を、図５を参照して説明する。図５は、ステップ230内の処理を示すフローである。
【００５３】
＜＜閾値処理（ステップ231）＞＞
まず、奇数エコーデータに対する閾値処理212と同様に、エコーデータの振幅値に対して閾値（例えば最大振幅値の２０％）を設定し、それ以下の周波数方向の点は、これ以降のステップから除外する（ステップ231）。
【００５４】
図６に、閾値処理前後のエコーデータを示す。図６（ａ）、（ｃ）は閾値処理前のエコーデータの振幅及び位相、図６（ｂ）、（ｄ）は閾値処理後のエコーデータの振幅及び位相である。図示するように、最大振幅の２０％以下のデータは、位相がゼロであるとして扱う。図６（ｄ）からもわかるように、エコー信号位相値は主値（−π〜＋πの間）の範囲で回転するため、周波数方向の位相変化を一次の線形変化にするために、主値回りを取り除く処理（位相アンラップ処理）が必要となる。
【００５５】
＜＜第１位相アンラップ処理（ステップ232）＞＞
位相アンラップは、主値回りが生じる点で本来の位相変化よりも±２πの変化が生じることを利用して、隣り合う点同士で位相差が十分に小さくなるよう、±２π加算する処理を行う。しかし、本実施形態では、±２π加算するのではなく、位相差を単位ベクトルに変換することで、位相差を−π〜＋πの間の主値内に納める簡易な位相アンラップ処理を行う。
【００５６】
具体的には、隣り合うデータ点の位相θn、θn-1の差分が−π〜＋πの間に入らない場合には、この差分を、式（６）、（７）に示すように、複素数（実部・虚部）に変換する。この実部と虚部のデータを用いて、式（８）にて位相値Δθを算出する。この位相値Δθを、一方のデータ点の位相θｎに加算し（式（９））、それに隣接するデータ点の位相θｎ+1を求める。これにより、θｎとθｎ+1間の±２πの変化を除去できる。
【数４】
【００５７】
この位相アンラップ処理は、±２π×Ｎの位相をＮを変えながら適当に加えて、目標位相との誤差が最も小さくなるＮを探索する処理を必要とする技術（特許文献２に記載された処理）に比べ、簡便な演算処理といえる。
【００５８】
ステップ232の位相アンラップ処理前後の位相データを図７に示す。図７（a）に示した主値内で回転する位相は、位相アンラップ処理によって、図７（b）に示すような本来の１次の線形な変化に戻すことができる。ただし、これら位相データは、ステップ231の閾値処理を行うことによって、処理から除外されたデータ点或いはデータ点の領域（除外領域）が存在するため、位相が連続しているデータ点の領域が除外領域によって隔てられる可能性がある。「位相が連続している」とは位相飛びがないことを意味し、この位相が連続しているデータ点の領域を「島領域」と呼ぶ。これら島領域が離れて存在する場合、位相アンラップ処理の対象となる隣り合う点の距離が離れることになり、正しく位相アンラップ処理を行えない場合がある。
【００５９】
例えば、図８に示す４つの島（島１〜島４）で構成されたデータの場合、小さな島々で構成される左から２番目の島２の右端の点における位相値と、左から３番目の島３の左端の点における位相値の差は、−π〜＋πの範囲の間に収まるため、上述した簡易アンラップ処理では、そのままの位相になるが、本来は−２π加算しなければならない。このように、本来隣接する点の位相データがないと連続性が確認できない。そこで、位相の連続性を持たせるための第２の位相アンラップ処理を行う。
【００６０】
＜＜第２の位相アンラップ処理（ステップ233〜235）＞＞
第２の位相アンラップ処理を行うために、まず複数の島のうち、連続するデータ点の数が最も多い最大島の探索する（ステップ233）。次いで、最大島の位相データに対して、最小二乗法を用いたフィッティング処理を行って、近似関数Dif(x)を算出し、その１次の傾きを求める（ステップ234）。最小二乗法を用いたフィッティングは、奇数エコーのフィッティング（ステップ214）で説明した式（２）〜（４）と同様の式を用いることができる。
【００６１】
上述した最大島の検索と、最大島の位相データを用いたフィッティングの様子を図９に示す。図９（ａ）は、連続するデータ点の点数が最も多い島の位相データを取り出した状態を示し、（ｂ）は（ａ）の位相データを用いて式（３）、（４）によりフィッティングして求めた一次関数（点線）を示す。
【００６２】
次にステップ234で求めた近似関数の傾きを用いて、島と島との間の位相値に連続性を持たせるための位相アンラップ処理（第２アンラップ処理）を行う（ステップ235）。具体的には、位相データがない（ゼロである）部分によって隔てられた２つの島の隣り合う点となる端部の位相値の一方を、他方の位相値と近似関数の傾きを用いて、位相アンラップ処理する。位相アンラップ処理は、周波数方向に配列するデータ点の左側から開始しても右側開始してもよいが、以下、左側から開始する場合を例に説明する。
【００６３】
まず２つの島のうち左側の島の右端の位相値と上記近似関数の傾きを用いて、右側の島の左端の位相値の目標値を算出する。目標値とはステップ２３４で求めた１次の位相変化に従えば、隣接する島の左端の点は本来この位相値になるであろう値を意味する。この目標値と実際の位相は大きくずれることはないため、この位相差が−π〜＋πに収まるよう、ステップ２３２の位相アンラプ処理と同様の演算処理を行う。即ち、近似関数から推定した位相目標値θtargetと、島の左端のデータ点の位相値θnとの差分を以下の式（１０）、（１１）で示すように、複素数に変換し、その実部と虚部のデータを用いて式（１２）により位相値Δθを算出する。この位相値Δθを位相目標値に加算し（式（１３））、そのデータ点の位相値とする。
【００６４】
【数５】
【００６５】
右側の島が、位相値が連続する複数のデータ点からなる場合には、島の各データ点全てに、式（１２）で算出した位相値Δθを、第１位相アンラップ処理後の位相値に加算する。
【００６６】
第２アンラップ処理の様子を図１０に示す。図示する例では、検索を左から行い、小さい島からなる島２と、その右側の島３との間でアンラップ処理する場合（ａ）及び島３と島４との間でアンラップ処理する場合（ｂ）を示している。島２と島３との間には、比較的大きいデータ欠落部分があり、最初の位相アンラップ処理（第１アンラップ処理）（ステップ232）では、この間のアンラップ処理が適切に行われていない。そこで島３の左端のデータ点について、ステップで求めたフィッティング関数の傾きから、この点の目標値を求め、実際の位相値との差分を求め、上記式（１０）〜（１３）から算出したΔθを島３の左端の実際の位相値に加算する。島３の全てのデータ点（図では２点）の位相値にも同じΔθを加算する。島３と島４との間でも、同様にして、島３の右端のデータ点の位相値とフィッティング関数の係数から、島４の左端のデータ点の位相値の目標値を算出し、目標値と実際の位相値との差が−π〜＋πの範囲の間に収まるように、Δθを求め加算する（図１０（ｂ））。図１０（ｃ）は、この第２アンラップ処理後の状態を示す図であり、閾値処理によって除去されたデータ点以外の全てのデータ点を通して、位相飛びのない（連続性のある）位相値になっている。
【００６７】
＜＜全データに対するフィッティング（ステップ236）＞＞
第２アンラップ処理を終えた後、除去されたデータ点以外の全てのエコーデータの位相値を用いて、再度最小二乗法を用いてフィッティング処理を行う（ステップ236）。フィッティング処理は、奇数エコーのフィッティングで用いた式と同様の式（１４）〜（１６）を用いて、フィッティング関数を求める。
【数６】
式中、DifOEは奇数エコー（基準エコー）と偶数エコー（基準エコー）との位相差分を表し、ｃはその傾き、ｄは切片を表す。
【００６８】
このフィッティング関数から算出される各データ点の位相値（フィッティングデータ）を、そのまま、補正用に用いることも可能であるが、本実施形態では、フィッティングデータを信号強度に応じて制限して使用し、フィッティングに伴うフィッティング誤差を極力小さくする。即ち、フィッティングには必ずフィッティング誤差が生じ、この誤差は直接Ｎ／２アーチファクトの発生に繋がる。そのため、位相変化に乱れが生じていない場合には、今回の抑制ターゲットとなるストリークアーチファクトも発生しないため、フィッティングデータを敢えて利用する必要はない。一方、信号値が弱い領域では、フィッティング誤差によって生じるＮ／２アーチファクトの信号強度も相対的に低い。このような理由から、本実施形態では、信号値が低く位相が乱れやすい領域・点に対してのみ、フィッティングデータを限定的に用い置換する。
【００６９】
＜＜フィッティングデータを用いたデータ置換（ステップ237〜239）＞＞
フィッティングデータの限定的な使用は、マルチスライス計測の場合に、エコー全体の信号が弱いスライスのみを対象とし、十分な信号値が得られたスライス内ではエコーデータのうち信号値の弱い領域のみを対象とする。
【００７０】
まず、スライス毎にエコー全体の信号が弱いか否かの判定を行う（ステップ237）。このため、マルチスライス間で各最大信号値の比較を行い、スライス間で閾値を設定する。閾値は、例えばマルチスライスの最大信号値に対して２０％とし、最大信号値が閾値以下のスライスは、エコー全体の信号が弱いと判定し、そのスライスのエコーデータ全体をフィッティングデータに置換する（ステップ239）。
【００７１】
十分に信号値が得られているスライスに対しては、既に閾値処理ステップ231によって信号値が弱い領域は除外されているので、この除外した領域に対してデータの置換を行う（ステップ238）。この置換処理は、全島の両端に対しては、ステップ゜236のフィッティング処理で求めた１次の傾きｃを用いて、外挿処理を行う。島と島との間については、既に第２位相アンラップ処理（ステップ235）において位相の連続性が保たれているため、線形補間処理を行えばよい。
【００７２】
図１１は、フィッティング処理の様子を示す図であり、図１１（ａ）は処理前の位相データ（第２アンラップ処理後の位相データ：図１０（ｃ））を示し、図１１（ｂ）は、十分に信号値が得られているスライスに対する処理、図１１（ｃ）は、エコー全体の信号値が低いスライスに対する処理を示す。図１１（ｂ）で点線で囲った部分がフィッティングを行った部分であり、両端はフィッティングデータを用いて外挿し、中間部分は線形補間している。また図１１（ｃ）に示すスライスでは、位相データ全体がフィッティングデータにより置換されている。このようなフィッティングを行うことにより、フィッティング誤差を最小に留め且つ信号値の弱い領域についても誤差の少ない位相値を得ることができる。
【００７３】
以上、図５を参照して、図２のステップ230の詳細を説明したが、このステップ230は、複数の奇数エコーと偶数エコーの組（例えば第１エコーと第２エコー、及び、第３エコーと第４エコー）について行ってもよいし、一組について行ってもよいが、少なくとも、次の位相累積処理（ステップ217）で基準となるエコーに対して行うことが好ましい。複数組について位相差分のフィッティングデータを求めた場合、その平均値を基準エコー間の位相差分データとしてもよい。
再び図２を参照して、ステップ230以降の処理を説明する。
【００７４】
＜＜偶数エコー間の位相差分（ステップ215）＞＞
偶数エコーについても、エコートレイン方向のピークずれを補正するために偶数エコー間の位相差分を行う必要がある。偶数エコー間の位相差分の処理は、奇数エコーと同様に、２つの偶数エコー、例えばエコートレインの２番目のエコーと４番目のエコーについて、エコーデータ（ｘ−ｋｙデータ）の位相差分を取り、周波数方向にフィッティングすることにより求めることもできるが、本実施形態では、ステップ214或いはステップ230で求めた奇数エコーと偶数エコーとの位相差分について、エコートレイン方向の位相変化（差分）を求める。これにより、Ｎ/２アーチファクトの原因である奇数エコーと偶数エコーのピークずれを併せて補正できるデータを得ることができる。具体的には、ステップ230で求めた奇数エコーの基準エコー（奇数基準エコー）と偶数エコーの基準エコー（偶数基準エコー）との位相差分に対する、他の隣接する奇数エコーと偶数エコーとの位相差分との差分を取ることになる。
【００７５】
この差分データに対して、閾値処理すること、最小二乗法を用いたフィッティングを行う（ステップ216）。この処理は、奇数エコーに対し行った処理（ステップ212、213）と同じであり、説明を省略する。ステップ216では、次式（１７）で示す一次関数が得られる。
（数７）
Difoe(x)=[DifOE(x)]+ex+g （１７）
式中、[DifOE(x)]は、ステップ230で最終的に得られた、偶数基準エコーと奇数基準エコーとの位相差分を表し、ｅは一次関数の傾き、ｇは切片を示す。
【００７６】
＜＜エコートレイン方向位相累積＞＞
位相累積は、奇数エコーの処理と同様に、基準エコーを開始点として、偶数エコー間の位相差分を順次エコートレイン方向に累積する。
【００７７】
奇数エコーから偶数エコーへの位相差分について位相累積を行った場合、位相累積された位相データには、奇数エコーと偶数エコーとのピークずれの情報と基準エコーに対するピークずれの情報とを含んでいる。そこでステップ218において、位相差分データに対し、位相累積後の奇数エコーデータの位相を加算し、偶数エコーの補正用データを得る。具体的には、式（１７）を用いて、次式（１８）の計算を行う。
（数８）
P’(x,2*i+2)=P’(x,2*i+1)+Difｏe(x) （１８）
式中「P’(x,ky)」は位相累積後の前計測データの位相変化である。
【００７８】
以上の処理（ステップ211〜218）は、マルチスライスの場合には、スライス毎に行い、スライス毎に奇数エコー及び偶数エコーの補正用データを得る。
【００７９】
＜本計測データの補正（ステップ222、223）＞
周波数方向にフーリエ変換した本計測データ（ｘ−ｋｙデータ）の位相に対し、上記ステップ217及び218で得た、前計測の奇数エコーの位相と偶数エコーの位相を加算することにより、ｘ−ｋｙ空間上で本計測データの位相が補正される（ステップ222）。ｘ−ｋｙ空間上での位相補正は、ｋ空間上でのピークずれの補正を意味し、この補正によってｋ空間上でのピークずれが解消される。マルチスライス計測の場合には、補正用データはスライス毎に用意されるので、本計測の各スライスに対し、対応するスライスの補正用データを用いて位相補正する。
【００８０】
位相補正後のｘ−ｋｙデータを位相方向にフーリエ変換することにより（ステップ223）、最終的に画像31が再構成される。再構成画像では、ピークずれに起因するＮ/２アーチファクトやストリークアーチファクトが大幅に低減される。
【００８１】
図１２は、前後の腹部ＤＷ−ＥＰＩ画像を示す図であり、（ａ）は、奇数エコー及び偶数エコーについて、それぞれ位相データに対する単純な一次フィッティングを行って位相補正した場合の画像、（ｂ）は、図５に示す処理（ステップ230）を行って位相補正した場合の画像である。図１２（ａ）ではストリークアーチファクトが発生していたが、（ｂ）では発生していないことが確認できる。
【００８２】
以上、本発明の一実施形態を説明し、その実施形態では、エコートレイン方向の位相累積を行う際の基準となるエコーデータの、奇数・偶数基準エコー間位相差分について、限定的なフィッティング（図４のステップ230）を行う場合を説明したが、奇数エコー間位相差分に対する周波数方向フィッティング（ステップ213）や偶数エコー間位相差分に対する周波数方向フィッティング（ステップ216）についても、同様の限定的なフィッティングを採用することも可能である。
【００８３】
また上記実施形態では、奇数エコーについて奇数エコー間位相差分を求め、偶数エコーについては奇数エコーとの差分に対しエコー間位相差分を行う場合を説明したが、逆に奇数エコーと偶数エコーの処理を逆にしても同様の結果を得ることができる。
【００８４】
さらに本発明の簡易位相アンラップ処理を含む周波数フィッティングの技術は、上述したＥＰＩ計測におけるエコーデータ間の位相差分に対してのみならず、他のエコー信号の位相補正や画像データの位相補正にも拡張して適用することが可能である。例えば、水脂肪分離計測におけるピークずれの補正、real IR シーケンスやラジアルシーケンスの補正、および、エコーピーク検出による傾斜磁場のディレイ測定にも適用することが可能である。
【産業上の利用可能性】
【００８５】
本発明によれば、ＥＰＩ系のパルスシーケンスを用いた撮像において、静磁場の不均一や傾斜磁場のオフセット誤差等の影響を低減し、ストリークアーチファクトのない画像を得ることができる。
本発明によれば、ＭＲＩ信号処理の基本的な位相アンラップ技術として、簡便で且つフィッティング誤差の少ない位相アンラップ技術が提供される。
【符号の説明】
【００８６】
１１・・・静磁場発生系、１２・・・傾斜磁場発生系、１３・・・送信系、１４・・・受信系、１５・・・シーケンサ、１６・・・再構成演算部（画像処理部）、１７・・・制御系、１８・・・表示部（ディスプレイ）。
【技術分野】
【０００１】
本発明は、磁気共鳴イメージング装置（以下、ＭＲＩ装置という）におけるエコープラナー（ＥＰＩ）計測時に生じるストリークアーチファクトを抑制する技術に関する。
【背景技術】
【０００２】
エコープラナー計測では、エコー信号を取得するときの読出し傾斜磁場の極性が異なることにより、エコートレインのうち奇数番目に取得されるエコー（以下、奇数エコーという）と偶数番目に取得されるエコー（以下、偶数エコーという）とにピークずれを生じる。このピークずれは、エコープラナー計測で得た画像上では、Ｎ/２アーチファクトとなって現れるが、前計測データを用いて奇数エコー及び／又は偶数エコーの位相を補正することにより、抑制することができる。
【０００３】
ストリークアーチファクトは、主として、位相エンコード方向の位相変化や二次的に混入する位相ノイズによって引き起こされる。特許文献１には、Ｎ/２アーチファクトの補正と同時に、ストリークアーチファクトを抑制する技術が開示されている。
【０００４】
特許文献１に記載されたストリークアーチファクト抑制方法では、位相エンコード傾斜磁場パルスを印加しない撮像で計測した前計測エコーを、奇数エコーと偶数エコーとに分離した後、奇・偶数エコーそれぞれのエコートレイン方向の位相変化を滑らかに接続することで、前計測データに含まれる位相乱れを抑制する。そして、この位相乱れが補正された前計測データを用いて本計測データを補正することにより、ストリークアーチファクトを低減する。
【０００５】
ＭＲＩで取得される信号や画像データの位相を滑らかに接続する処理には、一般にフィッティングが用いられるが、その前提技術として、隣接する点の位相が２πを超えた場合に、その整数倍（２π×Ｎ（Ｎ＝０，１，２・・・））を一方に加算又は減算する処理即ち主値回りした位相回転量を戻す処理（位相アンラップ処理と呼ばれる）があり、ＭＲＩの画像処理における必須技術となっている。位相アンラップ処理では、信号値が低いデータでは位相が乱れやすく信頼性のある処理が行えないという問題がある。
【０００６】
この問題に対し、特許文献２には、二次元画像データに対する位相アンラップ処理として、データを閾値処理し、所定の閾値以上のデータの領域（島領域）のデータについてのみ位相アンラップ処理し、島領域と島領域との間については、隣接する領域の位相との整合が取れるように所望の値を選択する処理が提案されている。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２００９−２７３５３０号公報
【特許文献２】国際特許公報ＷＯ０２/０８２９９４
【発明の概要】
【発明が解決しようとする課題】
【０００８】
特許文献１に記載された技術では、奇数エコー間又は偶数エコー間の位相差分及び奇・偶数エコー間の位相差分に対し、周波数方向の位相変化を直線でフィッティングし、それを位相エンコード方向に位相累積することにより、エコートレイン方向の位相変化を滑らかに接続している。
【０００９】
この技術では、周波数方向のフィッティングに単純な一次フィッティングを採用しているが、フィッティング処理には少なからずフィッティング誤差が生じ、Ｎ/２アーチファクトの原因となりうる。また奇・偶数エコーそれぞれについて、位相差分データを位相エンコード方向に位相累積する際には、開始点となるエコーを基準エコーとして、順次位相累積が行われるが、これら基準エコーの信号強度が極端に低く、両者の位相差に大きな乱れが生じる場合には、ストリークアーチファクトが生じてしまう可能性がある。
【００１０】
特許文献１の技術では、エコートレイン方向の中で最も信号強度が高いエコーを基準エコーとして採用することで極力高い信号強度が得られるような工夫がなされているが、周波数方向に対しては十分な信号強度が得られない場合がある。例えば、マルチプルアレイコイルのように空間的に不均一な感度分布を持った受信コイルを用いて撮像する場合には、局所的に信号強度が近い場所が存在するため、その領域における基準エコー間の位相差に大きな乱れが生じる可能性があり、それによりストリークアーチファクトが発生する。
【００１１】
一方、特許文献２に記載された技術は、信号強度が低く、位相情報の精度が低いデータを位相アンラップ処理から除外するものであり、精度の高いアンラップ処理を行うことができるが、位相アンラップ処理する際に、２π×ＮのＮを検索するステップや加算又は減算かを判定するステップ等が必須であり、位相アンラップ処理されない領域（島領域と島領域との間）については、評価関数を用いた評価を必要とし、複雑な処理を要する。
【課題を解決するための手段】
【００１２】
上記課題を解決するため、本発明のＭＲＩ装置は、エコープラナー計測によって取得したエコーデータを補正するための補正用データを、前計測データの奇エコーと偶数エコーとの位相差情報をフィッティンングして求める際に、フィッティングを信号強度の比較的低いデータ領域に限定的に行い、フィッティイング誤差の発生を抑制する。
【００１３】
また本発明のＭＲＩ装置は、奇数エコー及び偶数エコーそれぞれを位相エンコード方向に位相累積する際の基準エコーについて、奇・偶数エコー間の位相差についても上記限定されたフィッティングを行い、基準エコー間の位相乱れに起因するアーチファクトの発生を抑制する。
【００１４】
本発明の位相アンラップ処理方法は、隣り合う点の位相差を単位ベクトルに変換し、位相差として再計算する。これにより、もとの位相差に±２π×Ｎ（Ｎ＝０，１，２・・・）を加算して最小値を求める計算と同等な処理を簡易に行うことができる。
【００１５】
即ち、本発明のＭＲＩ装置は、検査対象の核磁気共鳴により発生するエコー信号を取得する撮像部と、前記撮像部が取得したエコー信号を処理し、核磁気共鳴画像を作成する画像処理部とを備え、前記画像処理部が処理するエコー信号は、第１の計測（前計測）で取得した第１計測データと、第２の計測(本計測)で取得した第２計測データとを含み、且つ第１計測データ及び第２計測データは、それぞれ、異なる極性の読出し傾斜磁場を用いて取得した第１群のエコーデータ（奇数エコーデータ）と第２群のエコーデータ（偶数エコーデータ）を含み、前記画像処理部は、前記第１計測データを用いて前記第２計測データを補正する補正手段を備え、前記補正手段は、前記第１計測データの第１群のエコーデータと第２群のエコーデータとの位相差分データを周波数方向に連続する位相差分データに変換して補正用データを作成し、その際、信号強度が閾値以下の周波数方向のデータ点または領域のみを予め求めたフィッティングデータで置換する。
【００１６】
本発明のＭＲＩ装置において、前記補正手段は、第１群及び第２群のエコーデータのうち、閾値以上のデータ点を選択し、選択されなかったデータ点または領域によって隔てられたデータ点の領域（島領域）毎に、位相差分を周波数方向に連続にする処理を行い、隣接する島領域と不連続な点或いは領域によって隔てられた島領域について、当該島領域の端部の位相値を、前記隣接する島領域から外挿によって予測される目標値と実際の位相値との差異が最小にするように決定する。
【００１７】
また本発明の位相アンラップ方法は、位相の主値回りを生じているデータに対し、主値回りを除去する方法であり、アンラップ方向の先端にあるデータ点の位相値（第１位相値）とそれに隣接するデータ点の位相値（第２位相値）の位相差が±２πを超えるときに、前記第１位相値と前記第２位相値との位相差を複素数に変換するステップと、前記複素数の実部と虚部の値から位相値を算出するステップと、前記算出された位相値を前記第１位相値に加算した値を、第２位相値に置き換えるステップと、を含む。
【発明の効果】
【００１８】
本発明のＭＲＩ装置によれば、ＥＰＩデータの補正用データを前計測によって得る際に、前計測エコーデータの周波数方向の位相フィッティングを信号値の弱い領域等に限定することにより、フィッティング誤差によるストリークアーチファクトが低減されたＭＲ画像を得ることができる。
【００１９】
また本発明のＭＲＩ装置によれば、奇数エコーと偶数エコーの基準となるエコー間の位相差に対してもフィッティング処理を加えることで、従来よりもアーチファクトを高精度に低減した画像を得ることができる。
【００２０】
さらに本発明によれば、ＭＲＩの基本技術である位相アンラップ方法として、精度が高く且つ従来の位相アンラップ方法に比べ簡便は位相アンラップ方法が提供される。
【図面の簡単な説明】
【００２１】
【図１】本発明が適用されるＭＲＩ装置の全体概要を示すブロック図
【図２】ＭＲＩ装置の動作の一実施形態を示すフローチャート
【図３】本計測のパルスシーケンスの一例を示す図
【図４】前計測のパルスシーケンスの一例を示す
【図５】図２のフローチャートのステップ230の詳細を示すフローチャート
【図６】実施形態における閾値処理の説明図で、（ａ）、（ｂ）はそれぞれ閾値処理前のエコーデータの振幅データと位相データを示し、（ｃ）、（ｄ）はそれぞれ閾値処理後のエコーデータの振幅データと位相データを示す。
【図７】実施形態における第１位相アンラップ処理の説明図で、（ａ）は差分後の位相データを示し、（ｂ）は位相アンラップ処理後の位相データを示す。
【図８】実施形態における最大島探索処理の説明図で、（ａ）は信号強度によって分けられた島の振幅データを示し、（ｂ）は各島の位相データを示す。
【図９】実施形態における最大島に対するフィッティング処理の説明図で、（ａ）は最大島の位相データ、（ｂ）は最大島の位相データを用いたフィッティング結果を示す図。
【図１０】実施形態における第２位相アンラップ処理の説明図で、（ａ）、（ｂ）は、それぞれ、目標値と実際の位相から求めた位相値に実際の位相を置換する様子を示す図、（ｃ）は第２位相アンラップ処理後の位相を示す図。
【図１１】実施形態におけるフィッティング処理の説明図で、（ａ）は第２位相アンラップ処理後の位相データを示す図、（ｂ）は（ａ）の位相データのうち、位相がゼロの領域或いは点のみをフィッティングデータで置換或いは線形補間する場合を示す図、（ｃ）は位相データ全体をフィッティングデータで置換する場合を示す図。
【図１２】実施形態による位相補正の効果の実例を示す図で、（ａ）は従来法による補正を行った場合、（ｂ）は本発明による補正を行った場合を示す。
【発明を実施するための形態】
【００２２】
以下、添付図面を参照して本発明のＭＲＩ装置の好ましい実施形態について説明する。
【００２３】
図１は、本発明を適用するＭＲＩ装置の全体構成を示すブロック図である。このＭＲＩ装置は、静磁場発生系１１、傾斜磁場発生系１２、送信系１３、受信系１４、シーケンサ１５、再構成演算部（画像処理部）１６、制御系１７および表示部（ディスプレイ）１８を備え、また図示しない操作部を備えている。主として、静磁場発生系１１、傾斜磁場発生系１２、送信系１３、受信系１４及びシーケンサ１５が、本発明における撮像部の機能を担っている。
【００２４】
静磁場発生系１１は、永久磁石、常電導磁石または超電導磁石のいずれかからなり、被検体１０が置かれる空間に均一な静磁場を発生させる。静磁場の方向は、通常、被検体１０の体軸方向または被検体１０の体軸と直交する方向である。
【００２５】
傾斜磁場発生系１２は、Ｘ、Ｙ、Ｚの３軸方向に巻かれた傾斜磁場コイル１２１と、これらの各々のコイルを磁化させる傾斜磁場電源１２２とからなり、シーケンサ１５からの命令に従って傾斜磁場電源１２２の各々のコイルを磁化させることによりＸ、Ｙ、Ｚの３軸方向の傾斜磁場Ｇｓ、Ｇｐ、Ｇｆを被検体１０に印加する。この傾斜磁場の加え方により、被検体１０の撮影断面が設定される。
【００２６】
送信系１３は、高周波発振器１３１、変調器１３２、高周波増幅器１３３及び高周波照射コイル１３４とからなり、傾斜磁場発生系１２で設定された被検体１０の断面を構成する原子の原子核に核磁気共鳴を起こさせるために、高周波発振器１３１から出力されたＲＦパルスを高周波増幅器１３３で増幅した後に、被検体１０に近接して設置された高周波照射コイル１３４に供給して被検体１０に照射する。
【００２７】
受信系１４は、高周波受信コイル１４１、受信回路１４２及びアナログ／ディジタル（以下「Ａ／Ｄ」という）変換器１４３からなり、送信系１３の高周波照射コイル１３４から照射された電磁波による被検体１０の原子核の核磁気共鳴によるＮＭＲ信号であるエコー信号を、被検体１０に近接して配置された高周波受信コイル１４１で検出し、受信回路１４２を介してＡ／Ｄ変換器１４３に入力し、ディジタル信号に変換して、さらにシーケンサ１５からの命令によるタイミングでサンプリングされたデータ（以下「サンプリングデータ」という）として、再構成演算部１６に送る。
【００２８】
制御系１７は、シーケンサ１５に指令を送り、シーケンサに設定された所定のパルスシーケンスに従い、上述した傾斜磁場発生系１２、送信系１３および受信系１４を制御するとともに、再構成演算部１６の動作を制御する。パルスシーケンスは、撮像方法によって種々のものがあり、それらはプログラム化されて制御系１７に格納されている。本実施形態のＭＲＩ装置では、ＥＰＩ系のパルスシーケンスを備えている。
【００２９】
再構成演算部１６は、サンプリングデータにフーリエ変換などの各種演算を行なうとともに、後述する位相変化や補正用データの算出、補正用データを用いた補正などを行なう。再構成演算部１６で再構成された画像は、表示部（ディスプレイ）１８に表示される。
【００３０】
次に上記構成におけるＭＲＩ装置の動作を説明する。図２に、動作の手順を示す。
本実施の形態のＭＲＩ装置では、画像化するエコー信号を計測する本計測（ステップ220）と、本計測とは別に実行される前計測（ステップ210）を実行する。ここでは、ステップ210の計測を「前計測」と呼んでいるが、前計測を実行するタイミングは、本計測の前と後のいずれでもよい。
【００３１】
本計測では、マルチエコー撮像によるパルスシーケンスを実行し、本計測によって得られたエコー信号（本計測エコー）３０を、前計測によって得られたエコー信号（前計測エコー）２０を用いて補正し（ステップ222）、画像化する（ステップ223）。本計測エコー３０の補正は、本計測エコー３０および前計測エコー２０をそれぞれ周波数方向にフーリエ変換（ステップ211、221）し、ｘ−ｋｙデータとし、差分することにより行なう。
【００３２】
前計測エコー２０を補正用データとして補正に供するに際し、奇数エコー２１と偶数エコー２２に分けて、それぞれの差分（位相差分）を周波数方向フィッティングし、エコートレイン方向に位相累積することにより、前計測エコーに２次的に混入する位相ノイズを除去するとともに時間変動を含む位相変化を反映した補正用データを得る。以下、各ステップの詳細を説明する。
【００３３】
＜計測ステップ＞
本計測（ステップ220）及び前計測（ステップ210）では、シングルショット又はマルチショットＥＰＩシーケンスにより、1回の励起で複数のエコー信号を計測する。図３及び図４に、本計測および前計測で用いられるパルスシーケンスの一例を示す。これらパルスシーケンスは公知のSpin Echo系ＥＰＩシーケンスであり、簡単に説明する。
【００３４】
まず検知する磁化を含む被検体にＲＦパルス301を照射すると同時にスライスを選択する傾斜磁場パルス302を印加し、画像化するスライスを選択後、前記傾斜磁場パルスによる位相分散を収束させるためのリフォーカス傾斜磁場パルス303を印加する。次いで、スピンエコーを生成するために１８０度ＲＦパルス305と２番目のスライス選択用の傾斜磁場パルス304、およびその前後に磁化の位相を分散させるためのクラッシャー傾斜磁場パルス（304、307）を印加する。
【００３５】
その後、周波数エンコード傾斜磁場のオフセットを与えるパルス309を印加した後、連続して反転する周波数エンコード傾斜磁場パルス311を印加し、反転する周波数エンコード傾斜磁場311の各周期内でエコー信号313が時系列的に発生するので、これをA/Dの時間範囲(312)の間おのおのサンプリングし時系列データを得る。反転する周波数エンコード傾斜磁場311の強度、印加タイミング、A/Dの時間範囲(312)は、本計測シーケンスと前計測と同じである。しかし、図３に示す本計測シーケンスでは、各エコー信号313に位相エンコードするための傾斜磁場パルス308、310が印加されるのに対し、図４に示す前計測シーケンスでは、位相エンコードは用いない。
【００３６】
シングルショットシーケンスの場合には、図示する繰り返し単位314を１回で、マルチショットシーケンスであれば、繰り返し単位314を複数回繰り返して、それぞれｋ空間全体を埋めるデータを取得する。1回の繰り返し単位で取得する複数のエコーをエコートレインと言い、奇数番目に取得されるエコーを奇数エコー、偶数番目に取得されるエコーを偶数エコーという。
【００３７】
このような計測をスライス傾斜磁場（302,306）を異ならせて行い、複数のスライスのデータを取得する。スライス数は特に限定されるものではないが、本実施の形態では、少なくとも前計測において、複数のスライスのデータを取得することが望ましい。複数のスライスのデータを利用することにより、より精度の高い補正用データを作成することができる。
【００３８】
このようなＥＰＩを基本とする計測シーケンスでは、周波数エンコード傾斜磁場の印加タイミングおよびサンプリング時間の誤差があると、エコーのピークがｋ空間中央（ｋｙ軸）に対してずれる。ずれる方向が、奇数エコーと偶数エコーとでは反対になるため、Ｎ/２アーチファクトを生じる。また位相エンコード方向のオフセットに誤差があると、エコーのピークのｋｙ軸からのずれが、エコー番号が大きくなるほど拡大する。以下の処理では、前計測エコーから、エコーのピークのずれに対応する位相情報（補正用データ）を得て、本計測エコーのピークずれを解消する。
【００３９】
＜周波数方向フーリエ変換（ステップ211、221）＞
本計測および前計測で得られたエコー信号２０、３０をそれぞれ周波数方向にフーリエ変換し、ｘ−ｋｙデータとする。周波数方向のフーリエ変換は、計測空間（ｋｘ−ｋｙ空間）におけるエコー信号のピークずれを、ｘ−ｋｙ空間上で位相補正するために行われる。即ち、計測空間におけるエコー信号のピークずれは、フーリエ変換後の空間（ｘ−ｋｙ空間）での位相変化に対応するので、前計測のｘ−ｋｙデータから求めた位相変化を用いて、ｘ−ｋｙ空間上で本計測データの位相を補正することにより、本計測データのピークずれが補正される。
【００４０】
＜補正用データの作成＞
前計測エコーをフーリエ変換した後のｘ−ｋｙデータ（以下、エコーデータという）を、奇数エコーと偶数エコーに分けて、それぞれについて補正用データを得る。以下、奇数エコーと偶数エコー各処理の詳細を説明する。
【００４１】
＜奇数エコーデータの処理＞
奇数エコー間の位相差分を求める。奇数エコー間或いは偶数エコー間で位相差分した場合、周波数エンコード傾斜磁場の印加タイミングとサンプリング時間の誤差はキャンセルされるので、位相エンコード方向のオフセット誤差のみが位相差分として求められる。この位相差分は、周波数方向に直線となるものであり、位相差分後の前計測データをフィッティングすることにより求められる。
【００４２】
具体的には、エコートレインを構成する複数のエコーのうち、２つのエコーデータ、例えばエコートレインの１番目のエコーと３番目のエコーとの位相差分を取る。位相差分処理を行うエコーの一方は、エコートレインのうち最も信号強度が高いＴＥの位置のエコーを用いることが好ましい。このようなエコーは、後述する位相累積処理を開始する基準エコーとなる。
【００４３】
＜＜閾値処理（ステップ212）＞＞
２つの奇数エコー間で位相差分する際に、信号強度が低く位相値の信頼性が低いデータをフィッティングに用いないために、エコーデータから信号強度が低いデータを除外する。除外するデータは、所定の閾値以下、例えば対象とするエコーデータの最大振幅値に対して所定の比率（例えば２０％）以下のデータとする。閾値Thresholdは次式（１）で表わされる。
（数１）
Threshold=Signalmax×０．２ （１）
ここで、Signalmaxはエコーデータ内の最大振幅値を意味する。
【００４４】
＜＜周波数方向のフィッティング（ステップ213）＞＞
次に２つのエコーデータを位相差分したデータに対し、周波数方向のフィッティングを行う。フィッティングは、最小二乗法等の公知の手法を用いることができ、式（２）で表わされる一次関数と実際の各データ点における位相差との差が最小となるように、式（３）及び（４）により、式（２）の係数ａ、ｂを求める。
【００４５】
【数２】
ここで、Difｏ(x)は奇数エコー間の位相差分をフィッティングした一次関数であり、Xiは周波数方向の座標、Ｙiは位相値で、nは周波数点数（データ点の点数）を意味する（以下、同様）。
【００４６】
なお、フィッティングに先立ち、位相の主値回りを除去するための位相アンラップ処理を行う。具体的には、隣合うデータ点の位相差を順次検索し、一つのデータ点の位相と次のデータ点の位相との差が２πを超えたときに、次のデータ点の位相に±２π×Ｎを加算する処理を行う。位相アンラップ処理は、一般的な手法を用いてもよいし、後述する簡易アンラップ処理を行ってもよい。
【００４７】
奇数エコー間の位相差分は、２つのエコーから求めることができるが、２つのエコーのみを用いるのではなく、複数の奇数エコー間で式（２）の係数を求めて、これを平均することが望ましい。
またマルチスライス計測の場合には、スライス毎に奇数エコー間の位相差分を求めて、これを平均してもよいし、特許文献１に記載されるように、スライスの信号を合成した後、エコー間の位相差分を求めてもよい。
【００４８】
＜＜位相累積（ステップ217）＞＞
２つのエコーデータ間の位相差は、位相エンコード方向のオフセット誤差に対応し、位相エンコード方向のオフセット誤差はエコートレイン方向に蓄積されるものである。そこで、上記フィッティングステップで求めた位相差を、エコートレイン方向に位相累積する。具体的には、位相累積を開始する基準となるエコーをエコートレインの（２×Ｂ+１）番目のエコーであるとすると、位相累積後の（２×ｉ+１）番目のエコーの位相変化Ｐは、次式（５）で表わされる。
【００４９】
（数３）
P(x,2×i+1)=P(x,2×B+1)+ Difｏ(x)×(i-B) （５）
【００５０】
以上の各ステップ211〜213、217により、エコートレインを構成する各奇数エコーについて、ピークを揃えるための補正用データが得られたことになる。
【００５１】
＜偶数エコーデータの処理＞
偶数エコーデータ２２については、周波数エンコード傾斜磁場の印加タイミングおよびサンプリング時間の誤差と、位相エンコード方向のオフセット誤差とを同時に補正するために、まず奇数エコーとの位相差分を行い（ステップ214）、位相差分データについて周波数方向のフィッティングを行う（ステップ230）。奇数エコーと偶数エコーとの位相差分も、理論的には一次直線でフィッティングできるものであるが、奇数エコー間或いは偶数エコー間の位相差分が小さい値の範囲に留まるのに対し、スケール（値が変化する範囲）が大きくなるため、適切な位相アンラップ処理を行えない可能性があり、単純な一次フィッティングを行った場合、フィッティング誤差が生じる可能性が高まる。
【００５２】
そこで偶数エコーデータ（奇数エコーとの位相差分データ）の処理ステップ230では、奇数エコー間の処理とは異なる位相アンラップ処理と限定されたフィッティングを行う。以下、ステップ230の処理の詳細を、図５を参照して説明する。図５は、ステップ230内の処理を示すフローである。
【００５３】
＜＜閾値処理（ステップ231）＞＞
まず、奇数エコーデータに対する閾値処理212と同様に、エコーデータの振幅値に対して閾値（例えば最大振幅値の２０％）を設定し、それ以下の周波数方向の点は、これ以降のステップから除外する（ステップ231）。
【００５４】
図６に、閾値処理前後のエコーデータを示す。図６（ａ）、（ｃ）は閾値処理前のエコーデータの振幅及び位相、図６（ｂ）、（ｄ）は閾値処理後のエコーデータの振幅及び位相である。図示するように、最大振幅の２０％以下のデータは、位相がゼロであるとして扱う。図６（ｄ）からもわかるように、エコー信号位相値は主値（−π〜＋πの間）の範囲で回転するため、周波数方向の位相変化を一次の線形変化にするために、主値回りを取り除く処理（位相アンラップ処理）が必要となる。
【００５５】
＜＜第１位相アンラップ処理（ステップ232）＞＞
位相アンラップは、主値回りが生じる点で本来の位相変化よりも±２πの変化が生じることを利用して、隣り合う点同士で位相差が十分に小さくなるよう、±２π加算する処理を行う。しかし、本実施形態では、±２π加算するのではなく、位相差を単位ベクトルに変換することで、位相差を−π〜＋πの間の主値内に納める簡易な位相アンラップ処理を行う。
【００５６】
具体的には、隣り合うデータ点の位相θn、θn-1の差分が−π〜＋πの間に入らない場合には、この差分を、式（６）、（７）に示すように、複素数（実部・虚部）に変換する。この実部と虚部のデータを用いて、式（８）にて位相値Δθを算出する。この位相値Δθを、一方のデータ点の位相θｎに加算し（式（９））、それに隣接するデータ点の位相θｎ+1を求める。これにより、θｎとθｎ+1間の±２πの変化を除去できる。
【数４】
【００５７】
この位相アンラップ処理は、±２π×Ｎの位相をＮを変えながら適当に加えて、目標位相との誤差が最も小さくなるＮを探索する処理を必要とする技術（特許文献２に記載された処理）に比べ、簡便な演算処理といえる。
【００５８】
ステップ232の位相アンラップ処理前後の位相データを図７に示す。図７（a）に示した主値内で回転する位相は、位相アンラップ処理によって、図７（b）に示すような本来の１次の線形な変化に戻すことができる。ただし、これら位相データは、ステップ231の閾値処理を行うことによって、処理から除外されたデータ点或いはデータ点の領域（除外領域）が存在するため、位相が連続しているデータ点の領域が除外領域によって隔てられる可能性がある。「位相が連続している」とは位相飛びがないことを意味し、この位相が連続しているデータ点の領域を「島領域」と呼ぶ。これら島領域が離れて存在する場合、位相アンラップ処理の対象となる隣り合う点の距離が離れることになり、正しく位相アンラップ処理を行えない場合がある。
【００５９】
例えば、図８に示す４つの島（島１〜島４）で構成されたデータの場合、小さな島々で構成される左から２番目の島２の右端の点における位相値と、左から３番目の島３の左端の点における位相値の差は、−π〜＋πの範囲の間に収まるため、上述した簡易アンラップ処理では、そのままの位相になるが、本来は−２π加算しなければならない。このように、本来隣接する点の位相データがないと連続性が確認できない。そこで、位相の連続性を持たせるための第２の位相アンラップ処理を行う。
【００６０】
＜＜第２の位相アンラップ処理（ステップ233〜235）＞＞
第２の位相アンラップ処理を行うために、まず複数の島のうち、連続するデータ点の数が最も多い最大島の探索する（ステップ233）。次いで、最大島の位相データに対して、最小二乗法を用いたフィッティング処理を行って、近似関数Dif(x)を算出し、その１次の傾きを求める（ステップ234）。最小二乗法を用いたフィッティングは、奇数エコーのフィッティング（ステップ214）で説明した式（２）〜（４）と同様の式を用いることができる。
【００６１】
上述した最大島の検索と、最大島の位相データを用いたフィッティングの様子を図９に示す。図９（ａ）は、連続するデータ点の点数が最も多い島の位相データを取り出した状態を示し、（ｂ）は（ａ）の位相データを用いて式（３）、（４）によりフィッティングして求めた一次関数（点線）を示す。
【００６２】
次にステップ234で求めた近似関数の傾きを用いて、島と島との間の位相値に連続性を持たせるための位相アンラップ処理（第２アンラップ処理）を行う（ステップ235）。具体的には、位相データがない（ゼロである）部分によって隔てられた２つの島の隣り合う点となる端部の位相値の一方を、他方の位相値と近似関数の傾きを用いて、位相アンラップ処理する。位相アンラップ処理は、周波数方向に配列するデータ点の左側から開始しても右側開始してもよいが、以下、左側から開始する場合を例に説明する。
【００６３】
まず２つの島のうち左側の島の右端の位相値と上記近似関数の傾きを用いて、右側の島の左端の位相値の目標値を算出する。目標値とはステップ２３４で求めた１次の位相変化に従えば、隣接する島の左端の点は本来この位相値になるであろう値を意味する。この目標値と実際の位相は大きくずれることはないため、この位相差が−π〜＋πに収まるよう、ステップ２３２の位相アンラプ処理と同様の演算処理を行う。即ち、近似関数から推定した位相目標値θtargetと、島の左端のデータ点の位相値θnとの差分を以下の式（１０）、（１１）で示すように、複素数に変換し、その実部と虚部のデータを用いて式（１２）により位相値Δθを算出する。この位相値Δθを位相目標値に加算し（式（１３））、そのデータ点の位相値とする。
【００６４】
【数５】
【００６５】
右側の島が、位相値が連続する複数のデータ点からなる場合には、島の各データ点全てに、式（１２）で算出した位相値Δθを、第１位相アンラップ処理後の位相値に加算する。
【００６６】
第２アンラップ処理の様子を図１０に示す。図示する例では、検索を左から行い、小さい島からなる島２と、その右側の島３との間でアンラップ処理する場合（ａ）及び島３と島４との間でアンラップ処理する場合（ｂ）を示している。島２と島３との間には、比較的大きいデータ欠落部分があり、最初の位相アンラップ処理（第１アンラップ処理）（ステップ232）では、この間のアンラップ処理が適切に行われていない。そこで島３の左端のデータ点について、ステップで求めたフィッティング関数の傾きから、この点の目標値を求め、実際の位相値との差分を求め、上記式（１０）〜（１３）から算出したΔθを島３の左端の実際の位相値に加算する。島３の全てのデータ点（図では２点）の位相値にも同じΔθを加算する。島３と島４との間でも、同様にして、島３の右端のデータ点の位相値とフィッティング関数の係数から、島４の左端のデータ点の位相値の目標値を算出し、目標値と実際の位相値との差が−π〜＋πの範囲の間に収まるように、Δθを求め加算する（図１０（ｂ））。図１０（ｃ）は、この第２アンラップ処理後の状態を示す図であり、閾値処理によって除去されたデータ点以外の全てのデータ点を通して、位相飛びのない（連続性のある）位相値になっている。
【００６７】
＜＜全データに対するフィッティング（ステップ236）＞＞
第２アンラップ処理を終えた後、除去されたデータ点以外の全てのエコーデータの位相値を用いて、再度最小二乗法を用いてフィッティング処理を行う（ステップ236）。フィッティング処理は、奇数エコーのフィッティングで用いた式と同様の式（１４）〜（１６）を用いて、フィッティング関数を求める。
【数６】
式中、DifOEは奇数エコー（基準エコー）と偶数エコー（基準エコー）との位相差分を表し、ｃはその傾き、ｄは切片を表す。
【００６８】
このフィッティング関数から算出される各データ点の位相値（フィッティングデータ）を、そのまま、補正用に用いることも可能であるが、本実施形態では、フィッティングデータを信号強度に応じて制限して使用し、フィッティングに伴うフィッティング誤差を極力小さくする。即ち、フィッティングには必ずフィッティング誤差が生じ、この誤差は直接Ｎ／２アーチファクトの発生に繋がる。そのため、位相変化に乱れが生じていない場合には、今回の抑制ターゲットとなるストリークアーチファクトも発生しないため、フィッティングデータを敢えて利用する必要はない。一方、信号値が弱い領域では、フィッティング誤差によって生じるＮ／２アーチファクトの信号強度も相対的に低い。このような理由から、本実施形態では、信号値が低く位相が乱れやすい領域・点に対してのみ、フィッティングデータを限定的に用い置換する。
【００６９】
＜＜フィッティングデータを用いたデータ置換（ステップ237〜239）＞＞
フィッティングデータの限定的な使用は、マルチスライス計測の場合に、エコー全体の信号が弱いスライスのみを対象とし、十分な信号値が得られたスライス内ではエコーデータのうち信号値の弱い領域のみを対象とする。
【００７０】
まず、スライス毎にエコー全体の信号が弱いか否かの判定を行う（ステップ237）。このため、マルチスライス間で各最大信号値の比較を行い、スライス間で閾値を設定する。閾値は、例えばマルチスライスの最大信号値に対して２０％とし、最大信号値が閾値以下のスライスは、エコー全体の信号が弱いと判定し、そのスライスのエコーデータ全体をフィッティングデータに置換する（ステップ239）。
【００７１】
十分に信号値が得られているスライスに対しては、既に閾値処理ステップ231によって信号値が弱い領域は除外されているので、この除外した領域に対してデータの置換を行う（ステップ238）。この置換処理は、全島の両端に対しては、ステップ゜236のフィッティング処理で求めた１次の傾きｃを用いて、外挿処理を行う。島と島との間については、既に第２位相アンラップ処理（ステップ235）において位相の連続性が保たれているため、線形補間処理を行えばよい。
【００７２】
図１１は、フィッティング処理の様子を示す図であり、図１１（ａ）は処理前の位相データ（第２アンラップ処理後の位相データ：図１０（ｃ））を示し、図１１（ｂ）は、十分に信号値が得られているスライスに対する処理、図１１（ｃ）は、エコー全体の信号値が低いスライスに対する処理を示す。図１１（ｂ）で点線で囲った部分がフィッティングを行った部分であり、両端はフィッティングデータを用いて外挿し、中間部分は線形補間している。また図１１（ｃ）に示すスライスでは、位相データ全体がフィッティングデータにより置換されている。このようなフィッティングを行うことにより、フィッティング誤差を最小に留め且つ信号値の弱い領域についても誤差の少ない位相値を得ることができる。
【００７３】
以上、図５を参照して、図２のステップ230の詳細を説明したが、このステップ230は、複数の奇数エコーと偶数エコーの組（例えば第１エコーと第２エコー、及び、第３エコーと第４エコー）について行ってもよいし、一組について行ってもよいが、少なくとも、次の位相累積処理（ステップ217）で基準となるエコーに対して行うことが好ましい。複数組について位相差分のフィッティングデータを求めた場合、その平均値を基準エコー間の位相差分データとしてもよい。
再び図２を参照して、ステップ230以降の処理を説明する。
【００７４】
＜＜偶数エコー間の位相差分（ステップ215）＞＞
偶数エコーについても、エコートレイン方向のピークずれを補正するために偶数エコー間の位相差分を行う必要がある。偶数エコー間の位相差分の処理は、奇数エコーと同様に、２つの偶数エコー、例えばエコートレインの２番目のエコーと４番目のエコーについて、エコーデータ（ｘ−ｋｙデータ）の位相差分を取り、周波数方向にフィッティングすることにより求めることもできるが、本実施形態では、ステップ214或いはステップ230で求めた奇数エコーと偶数エコーとの位相差分について、エコートレイン方向の位相変化（差分）を求める。これにより、Ｎ/２アーチファクトの原因である奇数エコーと偶数エコーのピークずれを併せて補正できるデータを得ることができる。具体的には、ステップ230で求めた奇数エコーの基準エコー（奇数基準エコー）と偶数エコーの基準エコー（偶数基準エコー）との位相差分に対する、他の隣接する奇数エコーと偶数エコーとの位相差分との差分を取ることになる。
【００７５】
この差分データに対して、閾値処理すること、最小二乗法を用いたフィッティングを行う（ステップ216）。この処理は、奇数エコーに対し行った処理（ステップ212、213）と同じであり、説明を省略する。ステップ216では、次式（１７）で示す一次関数が得られる。
（数７）
Difoe(x)=[DifOE(x)]+ex+g （１７）
式中、[DifOE(x)]は、ステップ230で最終的に得られた、偶数基準エコーと奇数基準エコーとの位相差分を表し、ｅは一次関数の傾き、ｇは切片を示す。
【００７６】
＜＜エコートレイン方向位相累積＞＞
位相累積は、奇数エコーの処理と同様に、基準エコーを開始点として、偶数エコー間の位相差分を順次エコートレイン方向に累積する。
【００７７】
奇数エコーから偶数エコーへの位相差分について位相累積を行った場合、位相累積された位相データには、奇数エコーと偶数エコーとのピークずれの情報と基準エコーに対するピークずれの情報とを含んでいる。そこでステップ218において、位相差分データに対し、位相累積後の奇数エコーデータの位相を加算し、偶数エコーの補正用データを得る。具体的には、式（１７）を用いて、次式（１８）の計算を行う。
（数８）
P’(x,2*i+2)=P’(x,2*i+1)+Difｏe(x) （１８）
式中「P’(x,ky)」は位相累積後の前計測データの位相変化である。
【００７８】
以上の処理（ステップ211〜218）は、マルチスライスの場合には、スライス毎に行い、スライス毎に奇数エコー及び偶数エコーの補正用データを得る。
【００７９】
＜本計測データの補正（ステップ222、223）＞
周波数方向にフーリエ変換した本計測データ（ｘ−ｋｙデータ）の位相に対し、上記ステップ217及び218で得た、前計測の奇数エコーの位相と偶数エコーの位相を加算することにより、ｘ−ｋｙ空間上で本計測データの位相が補正される（ステップ222）。ｘ−ｋｙ空間上での位相補正は、ｋ空間上でのピークずれの補正を意味し、この補正によってｋ空間上でのピークずれが解消される。マルチスライス計測の場合には、補正用データはスライス毎に用意されるので、本計測の各スライスに対し、対応するスライスの補正用データを用いて位相補正する。
【００８０】
位相補正後のｘ−ｋｙデータを位相方向にフーリエ変換することにより（ステップ223）、最終的に画像31が再構成される。再構成画像では、ピークずれに起因するＮ/２アーチファクトやストリークアーチファクトが大幅に低減される。
【００８１】
図１２は、前後の腹部ＤＷ−ＥＰＩ画像を示す図であり、（ａ）は、奇数エコー及び偶数エコーについて、それぞれ位相データに対する単純な一次フィッティングを行って位相補正した場合の画像、（ｂ）は、図５に示す処理（ステップ230）を行って位相補正した場合の画像である。図１２（ａ）ではストリークアーチファクトが発生していたが、（ｂ）では発生していないことが確認できる。
【００８２】
以上、本発明の一実施形態を説明し、その実施形態では、エコートレイン方向の位相累積を行う際の基準となるエコーデータの、奇数・偶数基準エコー間位相差分について、限定的なフィッティング（図４のステップ230）を行う場合を説明したが、奇数エコー間位相差分に対する周波数方向フィッティング（ステップ213）や偶数エコー間位相差分に対する周波数方向フィッティング（ステップ216）についても、同様の限定的なフィッティングを採用することも可能である。
【００８３】
また上記実施形態では、奇数エコーについて奇数エコー間位相差分を求め、偶数エコーについては奇数エコーとの差分に対しエコー間位相差分を行う場合を説明したが、逆に奇数エコーと偶数エコーの処理を逆にしても同様の結果を得ることができる。
【００８４】
さらに本発明の簡易位相アンラップ処理を含む周波数フィッティングの技術は、上述したＥＰＩ計測におけるエコーデータ間の位相差分に対してのみならず、他のエコー信号の位相補正や画像データの位相補正にも拡張して適用することが可能である。例えば、水脂肪分離計測におけるピークずれの補正、real IR シーケンスやラジアルシーケンスの補正、および、エコーピーク検出による傾斜磁場のディレイ測定にも適用することが可能である。
【産業上の利用可能性】
【００８５】
本発明によれば、ＥＰＩ系のパルスシーケンスを用いた撮像において、静磁場の不均一や傾斜磁場のオフセット誤差等の影響を低減し、ストリークアーチファクトのない画像を得ることができる。
本発明によれば、ＭＲＩ信号処理の基本的な位相アンラップ技術として、簡便で且つフィッティング誤差の少ない位相アンラップ技術が提供される。
【符号の説明】
【００８６】
１１・・・静磁場発生系、１２・・・傾斜磁場発生系、１３・・・送信系、１４・・・受信系、１５・・・シーケンサ、１６・・・再構成演算部（画像処理部）、１７・・・制御系、１８・・・表示部（ディスプレイ）。
【特許請求の範囲】
【請求項１】
検査対象の核磁気共鳴により発生するエコー信号を取得する撮像部と、前記撮像部が取得したエコー信号を処理し、核磁気共鳴画像を作成する画像処理部とを備えた磁気共鳴イメージング装置において、
前記画像処理部が処理するエコー信号は、第１の計測で取得した第１計測データと、第２の計測で取得した第２計測データとを含み、且つ第１計測データ及び第２計測データは、それぞれ、異なる極性の読出し傾斜磁場を用いて取得した第１群のエコーデータと第２群のエコーデータを含み、
前記画像処理部は、前記第１計測データを用いて前記第２計測データを補正する補正手段を備え、前記補正手段は、前記第１計測データの第１群のエコーデータと第２群のエコーデータとの位相差分データを周波数方向に連続する位相差分データに変換して補正用データを作成し、その際、信号強度が閾値以下の周波数方向のデータ点または領域のみを予め求めたフィッティングデータで置換することを特徴とする磁気共鳴イメージング装置。
【請求項２】
請求項１に記載の磁気共鳴イメージング装置において、
前記補正手段は、前記第１群のエコーデータと第２群のエコーデータとの位相差分データに加えて、前記第１群のエコーデータ間の位相差分と、前記第２群のエコーデータ間の位相差分と、を用いて補正用データを作成するものであり、
前記第１群のエコーデータと第２群のエコーデータとの位相差分は、前記第１群のエコーデータ間の位相差分を算出する際の基準となる第１群中のエコーデータと、前記第２群のエコーデータ間の位相差分を算出する際の基準となる第２群中のエコーデータと、の位相差分であることを特徴とする磁気共鳴イメージング装置。
【請求項３】
請求項１に記載の磁気共鳴イメージング装置において、
前記補正手段は、第１群及び第２群のエコーデータのうち、閾値以上のデータ点を選択し、選択されなかったデータ点または領域によって隔てられたデータ点の領域（島領域）毎に、位相差分を周波数方向に連続にする処理を行い、
隣接する島領域と不連続な点或いは領域によって隔てられた島領域について、当該島領域の端部の位相値を、前記隣接する島領域から外挿によって予測される目標値と実際の位相値との差異が最小にするように決定することを特徴とする磁気共鳴イメージング装置。
【請求項４】
請求項３に記載の磁気共鳴イメージング装置において、
前記補正手段は、前記島領域が複数存在するとき、データ点数が最大である島領域を検索し、前記最大の島領域の位相差情報をフィッティングすることにより近似関数を求め、当該近似関数の傾きを用いて前記目標値を算出することを特徴とする磁気共鳴イメージング装置。
【請求項５】
請求項１に記載の磁気共鳴イメージング装置において、
前記画像処理部が処理する第１計測データ及び第２計測データは、それぞれ、複数のスライスの計測データを含み、
前記補正手段は、複数のスライスの計測データのうち、信号強度がスライス間で設定された閾値以下のスライスについては、第１群のエコーデータと第２群のエコーデータとの位相差分を、フィッティングデータで置換することにより補正用データを作成することを特徴とする磁気共鳴イメージング装置。
【請求項６】
請求項１に記載の磁気共鳴イメージング装置において、
前記補正手段が行う位相差を周波数方向に連続する処理は、隣接する２つのデータ点の位相差を複素数に変換するステップと、前記複素数の実部と虚部の値から位相値を算出するステップと、前記隣接する２つのデータ点の一方のデータ点の位相に、前記算出された位相値を加算し、他方のデータ点の位相とするステップとを含むことを特徴とする磁気共鳴イメージング装置。
【請求項７】
位相の主値回りを生じているデータに対し、主値回りを除去する位相アンラップ方法であって、
アンラップ方向の先端にあるデータ点の位相値（第１位相値）とそれに隣接するデータ点の位相値（第２位相値）の位相差が±２πを超えるときに、
前記第１位相値と前記第２位相値との位相差を複素数に変換するステップと
前記複素数の実部と虚部の値から位相値を算出するステップと、
前記算出された位相値を前記第１位相値に加算した値を、第２位相値に置き換えるステップと、を含む位相アンラップ方法。
【請求項８】
請求項７に記載の位相アンラップ方法であって、
位相アンラップに先立ち、位相値が閾値以下であるデータ点の位相値をゼロとするステップを含み、
位相値がゼロであるデータ点からなる区間によって隔てられた、位相値が閾値を超えるデータ点からなる島領域のデータについて、それぞれ独立して位相アンラップ処理を行うステップと、
島領域のうち最もデータ点数の大きい島領域について、位相アンラップ処理後の位相値の傾斜をフィッティングするステップと、
隣り合う２つの島領域の端部のデータ点の一方について、他方のデータ点の位相値と前記フィッティングするステップで求めたフィッティング関数の傾斜を用いて位相目標値を算出するステップと、
前記一方のデータ点の実際の位相値とフィッティング関数を用いて算出された位相目標値との位相差を求め、当該位相差を前記位相目標値に加算した値を実際の位相値に置換するステップとを含むことを特徴とする位相アンラップ方法。
【請求項９】
請求項８に記載の位相アンラップ方法であって、
前記位相値がゼロである区間の位相値を、当該区間に隣接する島領域の端部の位相値を用いて線形補間して求めることを特徴とする位相アンラップ方法。
【請求項１０】
請求項８に記載の位相アンラップ方法であって、
前記位相値がゼロである区間の位相値を、前記フィッティングステップで求めたフィッティング関数を用いて決定することを特徴とする位相アンラップ方法。
【請求項１】
検査対象の核磁気共鳴により発生するエコー信号を取得する撮像部と、前記撮像部が取得したエコー信号を処理し、核磁気共鳴画像を作成する画像処理部とを備えた磁気共鳴イメージング装置において、
前記画像処理部が処理するエコー信号は、第１の計測で取得した第１計測データと、第２の計測で取得した第２計測データとを含み、且つ第１計測データ及び第２計測データは、それぞれ、異なる極性の読出し傾斜磁場を用いて取得した第１群のエコーデータと第２群のエコーデータを含み、
前記画像処理部は、前記第１計測データを用いて前記第２計測データを補正する補正手段を備え、前記補正手段は、前記第１計測データの第１群のエコーデータと第２群のエコーデータとの位相差分データを周波数方向に連続する位相差分データに変換して補正用データを作成し、その際、信号強度が閾値以下の周波数方向のデータ点または領域のみを予め求めたフィッティングデータで置換することを特徴とする磁気共鳴イメージング装置。
【請求項２】
請求項１に記載の磁気共鳴イメージング装置において、
前記補正手段は、前記第１群のエコーデータと第２群のエコーデータとの位相差分データに加えて、前記第１群のエコーデータ間の位相差分と、前記第２群のエコーデータ間の位相差分と、を用いて補正用データを作成するものであり、
前記第１群のエコーデータと第２群のエコーデータとの位相差分は、前記第１群のエコーデータ間の位相差分を算出する際の基準となる第１群中のエコーデータと、前記第２群のエコーデータ間の位相差分を算出する際の基準となる第２群中のエコーデータと、の位相差分であることを特徴とする磁気共鳴イメージング装置。
【請求項３】
請求項１に記載の磁気共鳴イメージング装置において、
前記補正手段は、第１群及び第２群のエコーデータのうち、閾値以上のデータ点を選択し、選択されなかったデータ点または領域によって隔てられたデータ点の領域（島領域）毎に、位相差分を周波数方向に連続にする処理を行い、
隣接する島領域と不連続な点或いは領域によって隔てられた島領域について、当該島領域の端部の位相値を、前記隣接する島領域から外挿によって予測される目標値と実際の位相値との差異が最小にするように決定することを特徴とする磁気共鳴イメージング装置。
【請求項４】
請求項３に記載の磁気共鳴イメージング装置において、
前記補正手段は、前記島領域が複数存在するとき、データ点数が最大である島領域を検索し、前記最大の島領域の位相差情報をフィッティングすることにより近似関数を求め、当該近似関数の傾きを用いて前記目標値を算出することを特徴とする磁気共鳴イメージング装置。
【請求項５】
請求項１に記載の磁気共鳴イメージング装置において、
前記画像処理部が処理する第１計測データ及び第２計測データは、それぞれ、複数のスライスの計測データを含み、
前記補正手段は、複数のスライスの計測データのうち、信号強度がスライス間で設定された閾値以下のスライスについては、第１群のエコーデータと第２群のエコーデータとの位相差分を、フィッティングデータで置換することにより補正用データを作成することを特徴とする磁気共鳴イメージング装置。
【請求項６】
請求項１に記載の磁気共鳴イメージング装置において、
前記補正手段が行う位相差を周波数方向に連続する処理は、隣接する２つのデータ点の位相差を複素数に変換するステップと、前記複素数の実部と虚部の値から位相値を算出するステップと、前記隣接する２つのデータ点の一方のデータ点の位相に、前記算出された位相値を加算し、他方のデータ点の位相とするステップとを含むことを特徴とする磁気共鳴イメージング装置。
【請求項７】
位相の主値回りを生じているデータに対し、主値回りを除去する位相アンラップ方法であって、
アンラップ方向の先端にあるデータ点の位相値（第１位相値）とそれに隣接するデータ点の位相値（第２位相値）の位相差が±２πを超えるときに、
前記第１位相値と前記第２位相値との位相差を複素数に変換するステップと
前記複素数の実部と虚部の値から位相値を算出するステップと、
前記算出された位相値を前記第１位相値に加算した値を、第２位相値に置き換えるステップと、を含む位相アンラップ方法。
【請求項８】
請求項７に記載の位相アンラップ方法であって、
位相アンラップに先立ち、位相値が閾値以下であるデータ点の位相値をゼロとするステップを含み、
位相値がゼロであるデータ点からなる区間によって隔てられた、位相値が閾値を超えるデータ点からなる島領域のデータについて、それぞれ独立して位相アンラップ処理を行うステップと、
島領域のうち最もデータ点数の大きい島領域について、位相アンラップ処理後の位相値の傾斜をフィッティングするステップと、
隣り合う２つの島領域の端部のデータ点の一方について、他方のデータ点の位相値と前記フィッティングするステップで求めたフィッティング関数の傾斜を用いて位相目標値を算出するステップと、
前記一方のデータ点の実際の位相値とフィッティング関数を用いて算出された位相目標値との位相差を求め、当該位相差を前記位相目標値に加算した値を実際の位相値に置換するステップとを含むことを特徴とする位相アンラップ方法。
【請求項９】
請求項８に記載の位相アンラップ方法であって、
前記位相値がゼロである区間の位相値を、当該区間に隣接する島領域の端部の位相値を用いて線形補間して求めることを特徴とする位相アンラップ方法。
【請求項１０】
請求項８に記載の位相アンラップ方法であって、
前記位相値がゼロである区間の位相値を、前記フィッティングステップで求めたフィッティング関数を用いて決定することを特徴とする位相アンラップ方法。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【公開番号】特開２０１３−７０８７３（Ｐ２０１３−７０８７３Ａ）
【公開日】平成２５年４月２２日（２０１３．４．２２）
【国際特許分類】
【出願番号】特願２０１１−２１３０８２（Ｐ２０１１−２１３０８２）
【出願日】平成２３年９月２８日（２０１１．９．２８）
【出願人】（０００１５３４９８）株式会社日立メディコ (1,613)
【Ｆターム（参考）】
【公開日】平成２５年４月２２日（２０１３．４．２２）
【国際特許分類】
【出願日】平成２３年９月２８日（２０１１．９．２８）
【出願人】（０００１５３４９８）株式会社日立メディコ (1,613)
【Ｆターム（参考）】
[ Back to top ]