(57)【要約】【課題】２-Ｐｏｉｎｔ  Ｄｉｘｏｎ法において、静磁場不均一が反映された位相画像を正しく算出する確度を高めることで、水画像と脂肪画像とが入れ替わることを防ぎ、操作者に判断する負担を負わせることなく、水と脂肪とを分離する技術を提供する。【解決手段】静磁場不均一に起因した位相変化である位相分布を、正しく算出する可能性を高める。そして、その影響を取り除いた信号値を用いて、水と脂肪とを分離する。これを実現するため、位相分布の算出において、領域拡大法を用いた位相アンラップ処理の開始点として、水信号が主となる画素が選択される確率が高まるよう、開始点候補を特定する。さらに、算出した位相分布の正誤を、開始点候補の特性を用いて高い確度で判別する。【選択図】  図３
2014090949.tif 000002
【特許請求の範囲】
【請求項１】
被検体が置かれる空間に静磁場を生成する静磁生成磁石と、前記被検体に高周波磁場パルスを送信する送信部と、高周波磁場パルスの送信により前記被検体から生じる核磁気共鳴信号を受信する受信部と、前記核磁気共鳴信号に位置情報を付加するための傾斜磁場を印加する傾斜磁場印加部と、前記送信部、前記傾斜磁場印加部及び前記受信部の動作を制御するとともに前記受信した核磁気共鳴信号に対して演算処理を行う計算機と、を備える磁気共鳴撮影装置であって、
前記計算機は、
予め定めたパルスシーケンスに従って、異なる２つのエコー時間の前記核磁気共鳴信号を得る計測部と、
前記２つの核磁気共鳴信号から、前記異なる２つのエコー時間の原画像を再構成する再構成部と、
前記異なる２つのエコー時間の原画像から静磁場不均一によって変化する位相の分布である位相分布を算出する位相分布算出部と、
前記異なる２つのエコー時間の原画像と前記位相分布とを用いて水信号と脂肪信号とを分離する水脂肪分離部と、を備え、
前記位相分布算出部は、
前記異なる２つのエコー時間の原画像を用い、ノイズの影響を受けた画素群をシード点候補画素群として特定するシード点候補特定部と、
前記異なる２つのエコー時間の原画像と前記シード点候補画素群の情報とを用い、仮の前記位相分布である仮位相分布を算出する仮位相分布算出部と、
前記仮位相分布の正誤を判定する正誤判定部と、
前記正誤判定部による判定結果に基づいて真の位相分布を決定する真位相分布決定部と、を備えること
を特徴とする磁気共鳴撮影装置。
【請求項２】
請求項１記載の磁気共鳴撮影装置であって、
前記シード点候補特定部は、
前記異なる２つのエコー時間の原画像の比として得た比画像の絶対値を絶対値比画像として算出する絶対値比画像算出部と、
前記シード点候補画素群を特定するための閾値を少なくとも１つ設定する閾値設定部と、
前記閾値を用いて前記絶対値比画像における前記シード点候補画素群を特定する画素群特定部と、を備えること
を特徴とする磁気共鳴撮影装置。
【請求項３】
請求項２記載の磁気共鳴撮影装置であって、
前記閾値設定部は、前記絶対値比画像の画素値の分布を用いて前記閾値を設定すること
を特徴とする磁気共鳴撮影装置。
【請求項４】
請求項２記載の磁気共鳴撮影装置であって、
前記閾値設定部は、前記異なる２つのエコー時間の原画像の信号雑音比に基づいて、信号比ヒストグラムを算出し、当該モデルに基づいて前記閾値を設定すること
を特徴とする磁気共鳴撮影装置。
【請求項５】
請求項３記載の磁気共鳴撮影装置であって、
前記閾値は、前記絶対値比画像の画素値の最頻値に当該画素値の標準偏差を加算した値以上の値に設定されること
を特徴とする磁気共鳴撮影装置。
【請求項６】
請求項３記載の磁気共鳴撮影装置であって、
前記閾値は、前記絶対値比画像の画素値の最頻値より大きい値に設定されること
を特徴とする磁気共鳴撮影装置。
【請求項７】
請求項２記載の磁気共鳴撮影装置であって、
前記閾値は、１．０以上の値に設定されること
を特徴とする磁気共鳴撮影装置。
【請求項８】
請求項２記載の磁気共鳴撮影装置であって、
前記仮位相分布算出部は、
前記比画像の位相値に対し、予め定めた基準画素から領域拡大法を用いた位相アンラップ処理を実行し、前記仮位相分布を算出する位相アンラップ部を備え、
前記基準画素は、前記シード点候補画素群の中から選択されること
を特徴とする磁気共鳴撮影装置。
【請求項９】
請求項２記載の磁気共鳴撮影装置であって、
前記正誤判定部は、
前記比画像の前記シード点候補画素群の各信号から前記仮位相分布として得た位相値をそれぞれ除去し、前記各シード点候補画素群の各画素の除去後信号値を得る位相除去部と、
前記各シード点候補画素群の中で、前記除去後信号値が正の値である画素数と負の値である画素数とをそれぞれ算出する画素数算出部と、
前記正の値の画素数と前記負の値の画素数とを比較し、前記仮位相分布の正誤を判定する判定部と、を備えること
を特徴とする磁気共鳴撮像装置。
【請求項１０】
請求項９記載の磁気共鳴撮影装置であって、
前記判定部は、前記正の値の画素数が前記負の値の画素数より多い場合、正と判定すること
を特徴とする磁気共鳴撮影装置。
【請求項１１】
請求項９記載の磁気共鳴撮影装置であって、
前記真位相分布決定部は、
前記判定部の判定結果が誤の場合、前記仮位相分布を補正し、補正後の位相分布を、前記真の位相分布と決定し、当該判定結果被核結果が正の場合、仮の位相分布を前記真の位相分布と決定すること
を特徴とする磁気共鳴撮影装置。
【請求項１２】
予め定めたパルスシーケンスに従って計測した異なる２つのエコー時間の核磁気共鳴信号からそれぞれ再構成された、前記異なる２つのエコー時間の原画像から、水信号と脂肪信号とを分離する水脂肪分離方法であって、
前記異なる２つのエコー時間の原画像を用い、ノイズの影響を受けた画素群をシード点候補画素群として特定し、
前記異なる２つのエコー時間の原画像と前記シード点候補画素群の情報とを用い、仮の前記位相分布である仮位相分布を算出し、
前記仮位相分布の正誤を判定し、
前記正誤判定部による判定結果に基づいて真の位相分布を決定し、
前記異なる２つのエコー時間の原画像と前記真の位相分布とを用いて、水信号と脂肪信号とを分離すること
を特徴とする水脂肪分離方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、核磁気共鳴撮影（以下、ＭＲＩという）技術に関し、特に水・脂肪の分離撮影技術に関する。
【背景技術】
【０００２】
ＭＲＩは、プロトン密度や核磁気共鳴信号の緩和情報の空間分布を画像化することで人体頭部、腹部、四肢等の形態または機能を、２次元もしくは３次元で撮影する。現在のＭＲＩにおいて、臨床で普及している撮影対象核種は、被検体の主たる構成物質である水素原子核（プロトン）である。
【０００３】
人体を計測対象とした場合、ＭＲＩで検出できる主なプロトンには水と脂肪がある。ところが、人体に存在する脂肪信号は、被検体の腹部、脊椎、四肢などにおける画像コントラストを低下させる。そのため、臨床において脂肪信号を抑制する方法が提案されている。
【０００４】
その一つとして水と脂肪の位相差を用いて脂肪を抑制するＤｉｘｏｎ法が知られている（例えば、非特許文献１参照）。水と脂肪には、分子構造の違いによって共鳴周波数が−３．５ｐｐｍだけ異なる「化学シフト」が生じる。この化学シフトによる周波数差は、磁場強度に比例し、磁場強度が１．５テスラのときはおよそ−２２４Ｈｚとなる。Ｄｉｘｏｎ法は、水と脂肪の周波数の違いによって生じる水と脂肪の核磁化ベクトルの位相差を利用し、取得した画像の中の、水信号と脂肪信号とを分離する。これにより、脂肪が抑制された、水信号からなる画像（水画像）を得る。なお、分離後の脂肪信号からなる画像を脂肪画像と呼ぶ。
【０００５】
Ｄｉｘｏｎ法では、原子核スピンを励起してから信号を取得するまでの時間（以下、エコー時間という）が異なる２つの画像データを取得する。これら２つの画像データは、水と脂肪が同位相（ｉｎ  ｐｈａｓｅ）となるエコー時間ｔinで取得した原画像Ｉinと、水と脂肪が逆位相（ｏｕｔ  ｏｆ  ｐｈａｓｅ）となるエコー時間ｔoutで取得した原画像Ｉoutとである。これらを加算および減算することによって、水画像の信号強度Ｗと脂肪画像の信号強度Ｆとをそれぞれ算出し、水信号と脂肪信号とを分離する。
【０００６】
ただし、Ｄｉｘｏｎ法では、静磁場空間に被検体を挿入した際に生じる静磁場不均一による位相変化が考慮されていない。空間的な静磁場不均一が存在すると、位置によって、化学シフトとは異なる位相変化が生じるため、単純な加算、減算処理で水・脂肪の完全な分離ができない。
【０００７】
そこで、静磁場不均一を考慮した水・脂肪分離方法として、３-Ｐｏｉｎｔ  Ｄｉｘｏｎ法と呼ばれる手法が知られている（例えば、非特許文献２参照）。３-Ｐｏｉｎｔ  Ｄｉｘｏｎ法は、水画像の信号強度Ｗ、脂肪画像の信号強度Ｆ、および静磁場不均一による周波数差という３つの変数に対して、エコー時間の異なる３つの画像から、最小二乗推定処理を繰り返し行うことで、３変数を確定し、水と脂肪を分離する方法である。ところが、３-Ｐｏｉｎｔ  Ｄｉｘｏｎ法は、３回の計測を実施する必要があるため、計測時間が長引く。
【０００８】
静磁場不均一を考慮し、かつ、計測時間を長引かせない手法として、２-Ｐｏｉｎｔ  Ｄｉｘｏｎ法と呼ばれる手法が知られている（例えば、非特許文献３参照）。２-Ｐｏｉｎｔ  Ｄｉｘｏｎ法では、Ｄｉｘｏｎ法と同様に、水と脂肪が同位相（ｉｎ  ｐｈａｓｅ）となるエコー時間ｔin、および、水と脂肪が逆位相（ｏｕｔ  ｏｆ  ｐｈａｓｅ）となるエコー時間ｔoutで、それぞれ原画像Ｉinおよび原画像Ｉoutを取得する。その後、静磁場不均一に起因した位相変化（以下、単に位相分布と呼ぶ）を、領域拡大法を用いて位相のアンラップ処理を行うことにより推定する。そして、この位相変化を除去後の原画像ＩinとＩoutとを加算、減算することで水画像と脂肪画像とに分離する。
【０００９】
なお、２-Ｐｏｉｎｔ  Ｄｉｘｏｎ法では、水と脂肪が同位相になるエコー時間ｔinで原画像Ｉinを計測したのち、傾斜磁場パルスを反転して、エコー時間ｔoutで水と脂肪が逆位相となる原画像Ｉoutを取得する。つまり１回の計測で同位相と逆位相の画像を取得する。
【先行技術文献】
【非特許文献】
【００１０】
【非特許文献１】Ｔｈｏｍａｓ  Ｄｉｘｏｎ，ｅｔ  ａｌ．  ”Ｓｉｍｐｌｅ  Ｐｒｏｔｏｎ  Ｓｐｅｃｔｒｏｓｃｏｐｉｃ  Ｉｍａｇｉｎｇ”  Ｒａｄｉｏｌｏｇｙ、１５３巻、１９８４年、１８９−１９４頁
【非特許文献２】Ｓｃｏｔｔ  Ｂ. Ｒｅｅｄｅｒ，ｅｔ  ａｌ．  ”Ｍｕｌｔｉｃｏｉｌ  Ｄｉｘｏｎ  Ｃｈｅｍｉｃａｌ  Ｓｐｅｃｉｅｓ  Ｓｅｐａｒａｔｉｏｎ  Ｗｉｔｈ  ａｎ  Ｉｔｅｒａｔｉｖｅ  Ｌｅａｓｔ−Ｓｑｕａｒｅｓ  Ｅｓｔｉｍａｔｉｏｎ  Ｍｅｔｈｏｄ”  Ｍａｇｎｅｔｉｃ  Ｒｅｓｏｎａｎｃｅ  ｉｎ  Ｍｅｄｉｃｉｎｅ、５１巻、２００４年、３５−４５頁
【非特許文献３】Ｊｉｎｇｆｅｉ  Ｍａ  “Ｂｒｅａｔｈ−Ｈｏｌｄ  Ｗａｔｅｒ  ａｎｄ  Ｆａｔ  Ｉｍａｇｉｎｇ  Ｕｓｉｎｇ  ａ  Ｄｕａｌ−Ｅｃｈｏ  Ｔｗｏ−Ｐｏｉｎｔ  Ｄｉｘｏｎ  Ｔｅｃｈｎｉｑｕｅ  Ｗｉｔｈａｎ Ｅｆｆｉｃｉｅｎｔ  ａｎｄ Ｒｏｂｕｓｔ  Ｐｈａｓｅ−Ｃｏｒｒｅｃｔｉｏｎ  Ａｌｇｏｒｉｔｈｍ”  Ｍａｇｎｅｔｉｃ  Ｒｅｓｏｎａｎｃｅ  ｉｎ  Ｍｅｄｉｃｉｎｅ、５２巻、２００４年、４１５−４１９頁
【発明の概要】
【発明が解決しようとする課題】
【００１１】
非特許文献３に記載の２-Ｐｏｉｎｔ  Ｄｉｘｏｎ法は、傾斜磁場パルスを反転することで、水と脂肪が同位相の画像と逆位相の画像とを１回の計測で取得するため、計測時間が長引くことはない。しかし、位相分布を算出する際、アンラップ処理の開始点（以下、シード点と呼ぶ。）が、脂肪が主成分である画素に設定されると、得られた位相分布は、脂肪の化学シフトによって位相値πのオフセットが加算されたものとなる。この位相分布を用いて、水画像と脂肪画像とを算出すると、両者は入れ替わったものとなる。従って、非特許文献３に記載の２-Ｐｏｉｎｔ  Ｄｉｘｏｎ法によれば、ユーザーは、分離した画像のどちらが水画像で、どちらが脂肪画像であるかを撮影のたびに判断しなければならない。
【００１２】
本発明は、上記事情に鑑みてなされたもので、２-Ｐｏｉｎｔ  Ｄｉｘｏｎ法において、静磁場不均一が反映された位相画像を正しく算出する確度を高めることで、水画像と脂肪画像とが入れ替わることを防ぎ、操作者に判断する負担を負わせることなく、水と脂肪とを分離する技術を提供する。
【課題を解決するための手段】
【００１３】
本発明は、静磁場不均一に起因した位相変化である位相分布を、正しく算出する可能性を高める。そして、静磁場不均一による位相分布の影響を取り除いた信号値を用いて、水と脂肪とを分離する。これを実現するため、位相分布の算出において、領域拡大法を用いた位相アンラップ処理の開始点として、水信号が主となる画素が選択される確率が高まるよう、開始点候補を特定する。さらに、算出した位相分布の正誤を、開始点候補の特性を用いて高い確度で判別する。
【００１４】
すなわち、被検体が置かれる空間に静磁場を生成する静磁生成磁石と、前記被検体に高周波磁場パルスを送信する送信部と、高周波磁場パルスの送信により前記被検体から生じる核磁気共鳴信号を受信する受信部と、前記核磁気共鳴信号に位置情報を付加するための傾斜磁場を印加する傾斜磁場印加部と、前記送信部、前記傾斜磁場印加部及び前記受信部の動作を制御するとともに前記受信した核磁気共鳴信号に対して演算処理を行う計算機と、を備える磁気共鳴撮影装置であって、前記計算機は、予め定めたパルスシーケンスに従って、異なる２つのエコー時間の前記核磁気共鳴信号を得る計測部と、前記２つの核磁気共鳴信号から、前記異なる２つのエコー時間の原画像を再構成する再構成部と、前記異なる２つのエコー時間の原画像から静磁場不均一によって変化する位相の分布である位相分布を算出する位相分布算出部と、前記異なる２つのエコー時間の原画像と前記位相分布とを用いて水信号と脂肪信号とを分離する水脂肪分離部と、を備え、前記位相分布算出手段は、前記異なる２つのエコー時間の原画像を用い、ノイズの影響を受けた画素群をシード点候補画素群として特定するシード点候補特定部と、前記異なる２つのエコー時間の原画像と前記シード点候補画素群の情報とを用い、仮の前記位相分布である仮位相分布を算出する仮位相分布算出部と、前記仮位相分布の正誤を判定する正誤判定部と、前記正誤判定部による判定結果に基づいて真の位相分布を決定する真位相分布決定部と、を備えることを特徴とする磁気共鳴撮影装置を提供する。
【００１５】
また、予め定めたパルスシーケンスに従って計測した異なる２つのエコー時間の核磁気共鳴信号からそれぞれ再構成された、前記異なる２つのエコー時間の原画像から、水信号と脂肪信号とを分離する水脂肪分離方法であって、前記異なる２つのエコー時間の原画像を用い、ノイズの影響を受けた画素群をシード点候補画素群として特定し、前記異なる２つのエコー時間の原画像と前記シード点候補画素群の情報とを用い、仮の前記位相分布である仮位相分布を算出し、前記仮位相分布の正誤を判定し、前記正誤判定部による判定結果に基づいて真の位相分布を決定し、前記異なる２つのエコー時間の原画像と前記真の位相分布とを用いて、水信号と脂肪信号とを分離することを特徴とする水脂肪分離方法を提供する。
【発明の効果】
【００１６】
本発明によれば、静磁場不均一が反映された位相画像を正しく算出する確度を高めることで、水画像と脂肪画像とが入れ替わることを防ぎ、操作者に判断する負担を負わせることなく、水と脂肪とを分離できる。
【図面の簡単な説明】
【００１７】
【図１】（ａ）は、本発明の実施形態のＭＲＩ装置の中の、水平磁場方式のＭＲＩ装置の外観図である。（ｂ）は、本発明の実施形態のＭＲＩ装置の中の、垂直磁場方式のＭＲＩ装置の外観図である。（ｃ）は、本発明の実施形態のＭＲＩ装置の中の、トンネル型磁石を斜めに傾けたＭＲＩ装置の外観図である。
【図２】本発明の実施形態のＭＲＩ装置の機能構成図である。
【図３】本発明の実施形態の計算機の機能ブロック図である。
【図４】本発明の実施形態の計測処理のフローチャートである。
【図５】本発明の実施形態のグラディエントエコー型パルスシーケンスの一例を説明するための説明図である。
【図６】本発明の実施形態のスピンエコー型パルスシーケンスの一例を説明するための説明図である。
【図７】本発明の実施形態の位相分布算出処理のフローチャートである。
【図８】本発明の位相分布算出部の機能ブロック図である。
【図９】本発明の実施形態のシード点候補画素群特定処理のフローチャートである。
【図１０】本発明の実施形態の信号比分布のグラフである。
【図１１】本発明の実施形態の仮位相分布算出処理のフローチャートである。
【図１２】本発明の実施形態の真位相分布決定処理のフローチャートである。
【発明を実施するための形態】
【００１８】
以下、本発明を適用する実施形態について説明する。以下、実施形態を説明するための全図において、特に言及しない限り、同一機能を有するものは同一符号を付し、その繰り返しの説明は省略する。
【００１９】
まず、本実施形態の磁気共鳴撮影装置（ＭＲＩ装置）について説明する。図１は、本実施形態のＭＲＩ装置の外観図である。図１（ａ）は、ソレノイドコイルで静磁場を生成するトンネル型磁石を用いた水平磁場方式のＭＲＩ装置１００である。図１（ｂ）は、開放感を高めるために磁石を上下に分離したハンバーガー型（オープン型）の垂直磁場方式のＭＲＩ装置１２０である。また、図１（ｃ）は、図１（ａ）と同じトンネル型磁石を用い、磁石の奥行を短くし、かつ、斜めに傾けることによって、開放感を高めたＭＲＩ装置１３０である。本実施形態では、これらの外観を有するＭＲＩ装置のいずれを用いることもできる。なお、これらは一例であり、本実施形態のＭＲＩ装置はこれらの形態に限定されるものではない。本実施形態では、装置の形態やタイプを問わず、公知の各種のＭＲＩ装置を用いることができる。以下、特に区別する必要がない場合は、ＭＲＩ装置１００で代表する。
【００２０】
図２は、本実施形態のＭＲＩ装置１００の機能構成図である。本図に示すように、本実施形態のＭＲＩ装置１００は、被検体１０１が置かれる空間に、静磁場を生成する、例えば、静磁場コイルなどの静磁場生成磁石（以下、静磁場コイルという）１０２と、静磁場分布を調整するシムコイル１０４と、被検体１０１の計測領域に対し高周波磁場を送信する送信用高周波コイル１０５（以下、単に送信コイルという）と、被検体１０１から生じる核磁気共鳴信号を受信する受信用高周波コイル１０６（以下、単に受信コイルという）と、被検体１０１から生じる核磁気共鳴信号に位置情報を付加するため、ｘ方向、ｙ方向、ｚ方向それぞれに傾斜磁場を印加する傾斜磁場コイル１０３と、送信機１０７と、受信機１０８と、計算機１０９と、傾斜磁場用電源部１１２と、シム用電源部１１３と、シーケンス制御装置１１４と、を備える。
【００２１】
静磁場コイル１０２は、図１（ａ）、図１（ｂ）、図１（ｃ）にそれぞれ示した各ＭＲＩ装置１００、１２０、１３０の構造に応じて、種々の形態のものが採用される。傾斜磁場コイル１０３及びシムコイル１０４は、それぞれ傾斜磁場用電源部１１２及びシム用電源部１１３により駆動される。なお、本実施形態では、送信コイル１０５と受信コイル１０６とに別個のものを用いる場合を例にあげて説明するが、送信コイル１０５と受信コイル１０６との機能を兼用する１つのコイルで構成してもよい。送信コイル１０５が照射する高周波磁場は、送信機１０７により生成される。受信コイル１０６が検出した核磁気共鳴信号は、受信機１０８を通して計算機１０９に送られる。
【００２２】
シーケンス制御装置１１４は、傾斜磁場コイル１０３の駆動用電源である傾斜磁場用電源部１１２、シムコイル１０４の駆動用電源であるシム用電源部１１３、送信機１０７及び受信機１０８の動作を制御し、傾斜磁場、高周波磁場の印加および核磁気共鳴信号の受信のタイミングを制御する。制御のタイムチャートはパルスシーケンスと呼ばれ、計測に応じて予め設定され、後述する計算機１０９が備える記憶装置等に格納される。
【００２３】
計算機１０９は、ＭＲＩ装置１００全体の動作を制御するとともに、受信した核磁気共鳴信号に対して様々な演算処理を行う。本実施形態では、異なるエコー時間の原画像や分離後の水画像、脂肪画像、静磁場不均一による位相変化を示す位相分布を生成する。計算機１０９は、ＣＰＵ、メモリ、記憶装置などを備える情報処理装置であり、計算機１０９にはディスプレイ１１０、外部記憶装置１１１、入力装置１１５などが接続される。
【００２４】
ディスプレイ１１０は、演算処理で得られた結果等をオペレータに表示するインタフェースである。入力装置１１５は、本実施形態で行われる演算処理に必要な条件、パラメータ等をオペレータが入力するためのインタフェースである。外部記憶装置１１１は、記憶装置とともに、計算機１０９が実行する各種の演算処理に用いられるデータ、演算処理により得られるデータ、入力された条件、パラメータ等を保持する。
【００２５】
以下、本実施形態の計算機１０９が実現する機能について説明する。図３は、本実施形態の計算機１０９の機能ブロック図である。本実施形態の計算機１０９は、計測部３１０と、再構成部３２０と、位相分布算出部３３０と、水−脂肪分離部３４０とを備える。
【００２６】
計測部３１０は、パルスシーケンスにしたがってシーケンス制御装置１１４を動作させるとともに、ＭＲＩ装置１００の各部を制御して計測を行い、予め定めたエコー時間の核磁気共鳴信号（エコー信号）を得る。本実施形態では、異なる２つのエコー時間のエコー信号を得る。異なる２つのエコー時間として、水と脂肪とが同位相となるエコー時間ｔinと、水と脂肪とが逆位相となるエコー時間ｔoutとを用いる。
【００２７】
再構成部３２０は、取得したエコー信号をｋ空間上に配置し、フーリエ変換により画像を再構成する。本実施形態では、異なる２つのエコー時間ｔinおよびｔoutの原画像IinおよびＩoutをそれぞれ再構成する。
【００２８】
位相分布算出部３３０は、算出した原画像から静磁場不均一によって変化する位相の分布（位相分布）を算出する。本実施形態では、異なる２つのエコー時間の原画像を複素除算して得た比画像を用いて位相分布を算出する。
【００２９】
水−脂肪分離部３４０は、位相分布算出部３３０で算出した位相分布と異なる２つのエコー時間の原画像とから水信号と脂肪信号とを分離する。そして、水信号からなる水画像と脂肪信号からなる脂肪画像とを生成する。
【００３０】
なお、計算機１０９が実現する各種の機能は、記憶装置が保持するプログラムを、ＣＰＵがメモリにロードして実行することにより実現される。また、計算機１０９が実現する各種の機能のうち、少なくとも一つの機能は、ＭＲＩ装置１００とは独立した、情報処理装置であって、ＭＲＩ装置１００とデータの送受信が可能な情報処理装置により実現されていてもよい。
【００３１】
本実施形態の水、脂肪分離処理の説明に先立ち、上記Ｄｉｘｏｎ法による水、脂肪分離処理の原理、および、２ｐｏｉｎｔ  Ｄｉｘｏｎ法による処理の概要を説明する。
【００３２】
上述のように、Ｄｉｘｏｎ法では、エコー時間ｔinで取得した原画像Ｉinと、エコー時間ｔoutで取得した原画像Ｉoutとを加算および減算することによって、水画像の信号強度Ｗと脂肪画像の信号強度Ｆとをそれぞれ算出し、水信号と脂肪信号とを分離する。
【００３３】
ここで、水と脂肪の周波数差をｄｆとしたとき、水と脂肪が同位相になる時間ｔinは、ｎ／ｄｆ、水と脂肪が逆位相になる時間ｔoutは、（２ｎ＋１）／（２ｄｆ）となる。なお、ｎは整数である。したがって、エコー時間ｔinとｔoutとの差をｄｔとすると、ｄｔ＝１／（２ｄｆ）となる。
【００３４】
水画像の信号強度をＷ、脂肪画像の信号強度をＦとすると、静磁場不均一に起因する位相変化がない場合、ＩinとＩoutはそれぞれ以下の式（１−１）および式（１−２）で表現できる。
【数１】
2014090949.tif 000003
Ｄｉｘｏｎ法では、この２つの信号を、以下の式（２−１）および式（２−２）のように加算および減算することによって、水画像の信号強度Ｗと脂肪画像の信号強度Ｆとをそれぞれ算出し、水信号と脂肪信号とを分離する。
【数２】
2014090949.tif 000004
【００３５】
一方、２-Ｐｏｉｎｔ  Ｄｉｘｏｎ法では、上述のように、さらに、静磁場不均一に起因した位相変化（以下、単に位相分布と呼ぶ）を、領域拡大法を用いて、位相のアンラップ処理を行うことにより推定する。領域拡大法を用いた位相のアンラップ処理では、まず、開始点（シード点）とする基準画素を設定し、そのシード点と隣り合う画素の位相を、シード点との位相差に応じて位相変化が滑らかになるよう決定する。位相を決定した画素を次のシード点とし、画像の全画素の位相が決定するまで同様の処理を繰り返し、位相分布を算出する。そして、得られた位相分布の位相値を、原画像ＩinまたはＩoutから除去したのち、原画像ＩinとＩoutとを加算、減算することで水画像と脂肪画像とに分離する。
【００３６】
以下、上記本実施形態の計算機１０９の各機能による、計測全体の処理（計測処理）の流れについて簡単に説明する。図４は、本実施形態の計測処理の流れを説明するためのフローチャートである。
【００３７】
本実施形態では、予め定めた異なる２つのエコー時間ｔinおよびｔoutの原画像IinおよびＩoutを撮影し、原画像IinおよびＩoutから静磁場不均一による位相分布を算出し、位相分布と原画像IinおよびＩoutから、水と脂肪とを分離した、水画像および脂肪画像を算出する。本実施形態では、静磁場不均一による位相分布を算出する際、水が主成分である画素がシード点として設定される可能性を高める。また、得られた位相分布の正誤も判定する。
【００３８】
まず、計測部３１０は、予め定められたパルスシーケンスにしたがってシーケンス制御装置１１４を制御し、異なる２つのエコー時間ｔinおよびｔoutのエコー信号を計測する（ステップＳ１００１）。
【００３９】
その後、再構成部３２０は、得られたエコー信号を、それぞれ、ｋ空間上に配置してフーリエ変換することにより、異なる２つのエコー時間ｔinおよびｔoutの原画像IinおよびＩoutをそれぞれ再構成する（ステップＳ１００２）。
【００４０】
位相分布算出部３３０は、２つの原画像IinおよびＩoutの比をとり、その比画像の位相情報から静磁場不均一に起因した位相変化の分布(位相分布)を算出する（ステップＳ１００３）。
【００４１】
その後、水−脂肪分離部３４０は、算出した位相分布と２つの原画像IinおよびＩoutとから、水信号と脂肪信号とを分離し（ステップＳ１００４）、分離後の水信号からなる水画像と脂肪信号からなる脂肪画像とをディスプレイ１１０に表示する（ステップＳ１００５）。なお、ステップＳ１００４の処理を、水−脂肪分離処理と呼ぶ。また、分離後の水画像および脂肪画像を、分離画像と呼ぶ。
【００４２】
水画像および脂肪画像をディスプレイ１１０に表示する際、両者を区別可能とする情報を併せて表示してもよい。区別可能とする情報は、例えば、「水画像」、「脂肪画像」といったキャプションなどとする。また、必要に応じて、ステップＳ１００３で算出した位相分布をディスプレイ１１０に表示するよう構成してもよい。
【００４３】
ここで、計測部３１０がステップＳ１００１の計測において実行するパルスシーケンスの一例について説明する。図５に、本実施形態において採用されるパルスシーケンス５１０のタイムチャートの例を示す。このパルスシーケンス５１０は、グラディエントエコー型のパルスシーケンスである。なお、パルスシーケンス５１０において、ＲＦはＲＦパルスの、Ｇｓはスライス選択傾斜磁場の、Ｇｐは位相エンコード傾斜磁場の、Ｇｒは読み出し傾斜磁場の、それぞれ印加タイミングを、エコーは、エコー信号の取得タイミングをそれぞれ示す。以下、本明細書の各パルスシーケンスにおいて、同様である。
【００４４】
パルスシーケンス５１０では、１回の繰り返し時間内に以下の手順でエコー信号の計測を行う。まずＲＦパルス５１１を照射し、被検体１０１の水素原子核スピンを励起する。この際、被検体１０１の特定のスライスを選択するためにスライス選択傾斜磁場（Ｇｓ）５１２をＲＦパルス５１１と同時に印加する。続いてエコー信号に位相エンコードするための位相エンコード傾斜磁場（Ｇｐ）５１３を印加する。その後、最初のＲＦパルス５１１照射から時間ＴＥ1後に、読み出し傾斜磁場（Ｇｒ）５１４を印加してエコー信号（第一エコー）５１５を計測する。更に、第１エコー５１５の計測から時間ｄｔ後に、極性の反転した読み出し傾斜磁場（Ｇｒ）５１６を印加してエコー信号（第二エコー）５１７を計測する。第二エコー５１７を計測するエコー時間ＴＥ2は、ＴＥ2＝ＴＥ1＋ｄｔを満たす。
【００４５】
第一エコー５１５を計測する際のエコー時間ＴＥ1、および、第二エコー５１７を計測する際のエコー時間ＴＥ2は、上述したように，それぞれ、水と脂肪が同位相になる時間ｔinおよび水と脂肪が逆位相になる時間ｔoutとなるように設定する。ここで、水と脂肪の周波数差をｄｆとしたとき、水と脂肪が同位相になる時間ｔinは、ｎ／ｄｆ、水と脂肪が逆位相になる時間ｔoutは、（２ｎ＋１）／（２ｄｆ）である。なお、ｎは整数である。
【００４６】
本実施形態のパルスシーケンス５１０では、上述の条件を満たすエコー時間ＴＥ1およびエコー時間ＴＥ2が選択される。なお、エコー時間ｔinとｔoutの計測順序は問わない。
【００４７】
本実施形態の計測では、パルスシーケンス５１０を、位相エンコード傾斜磁場５１３の強度を変化させながら所定回数繰り返す。例えば１２８回、２５６回等、画像再構成に必要な数のエコー信号を繰り返し取得する。繰り返し回数分の第一エコーにより、１つの原画像（第一原画像）が形成され、繰り返し回数分の第二エコーにより、他の１つの原画像（第二原画像）が形成される。これらは後述する水画像および脂肪画像を算出するための演算用の原画像として使用される。
【００４８】
図６に、本実施形態に採用されるパルスシーケンスの他の例のタイムチャートを示す。このパルスシーケンス５２０は、スピンエコー型のパルスシーケンスである。
【００４９】
パルスシーケンス５２０では、1回の繰り返し時間内に以下の手順でエコー信号の計測を行う。まず、ＲＦパルス５２１を照射し、被検体１０１のスピンを励起する。この際、被検体１０１の特定のスライスを選択するためにスライス選択傾斜磁場（Ｇｓ）５２２をＲＦパルス５２１と同時に印加する。続いてエコー信号に位相エンコードするための位相エンコード傾斜磁場（Ｇｐ）５２３を印加し、更にスピンを反転させるためのＲＦパルス５２８をスライス選択傾斜磁場（Ｇｓ）５２９とともに照射する。その後、最初のＲＦパルス５２１照射から時間ＴＥ1後に、読み出し傾斜磁場（Ｇｒ）５２４印加してエコー信号（第一エコー）５２５を計測する。更に、第１エコー５２５の計測から時間ｄｔ後に、極性の反転した読み出し傾斜磁場（Ｇｒ）５２６を印加してエコー信号（第二エコー）５２７を計測する。
【００５０】
第一エコー５２５を計測する際のエコー時間ＴＥ1および、第２エコー５２７を計測する際のエコー時間ＴＥ2は、グラディエント型パルスシーケンスと同様の条件を満たすように選択される。
【００５１】
なお、図５および図６に示す例では、グラディエントエコー型のパルスシーケンス５１０の場合も、スピンエコー型のパルスシーケンス５２０の場合も、１回の繰り返し時間内にエコー時間の異なる２つの信号（ここでは、第一エコーおよび第二エコー）を計測している。しかし、１回の繰り返し時間には、１つのエコー信号を計測するよう構成してもよい。この場合、２回の繰り返し時間を用いて、第一エコーおよび第二エコーを計測する。
【００５２】
再構成部３２０は、上記ステップＳ１００２において、ステップＳ１００１で計測したエコー時間ｔinおよびｔoutのエコー信号をｋ空間上に各々配置し、フーリエ変換する。これにより、再構成部３２０は、異なる２つのエコー時間（ｔinおよびｔout）の原画像Ｉin（第一原画像）およびＩout（第二原画像）をそれぞれ算出する。
【００５３】
次に、上記ステップＳ１００３の、位相分布算出部３３０による位相分布の算出について説明する。本実施形態の位相分布算出部３３０は、まず、異なる２つのエコー時間の原画像ＩinおよびＩoutを複素除算して得た比画像において、位相分布の算出を開始する画素（以下、シード点）を特定する。そして、比画像の位相値を用いて、シード点から拡大領域法により各画素の位相値を決定し、位相分布を算出する。
【００５４】
これを実現するため、本実施形態の位相分布算出部３３０は、図３に示すように、シード点候補となる画素群を特定するシード点候補特定部３３１と、シード点候補として特定した画素群の中の１の画素をシード点として仮の位相分布を算出する仮位相分布算出部３３２と、仮位相分布算出部３３２が算出した仮の位相分布の正誤を判定する正誤判定部３３３と、正誤判定部３３３の正誤判定結果に基づいて、真の位相分布を決定する真位相分布決定部３３４と、を備える。
【００５５】
上記各機能による、本実施形態の位相分布算出処理の流れの概略を説明する。図７は、本実施形態の位相分布算出処理の流れを説明するためのフローチャートである。
【００５６】
シード点候補特定部３３１は、２つの原画像ＩinおよびＩoutを用いて、シード点の候補となる画素群（シード点候補画素群）を特定する（ステップＳ１１０１）。次に、仮位相分布算出部３３２は、特定したシード点候補画素群の中の任意の１点をシード点として、拡大領域法を用いて位相アンラップ処理を行い、仮の位相分布を算出する（ステップＳ１１０２）。次に、正誤判定部３３３は、仮の位相分布の正誤を判定する（ステップＳ１１０３）。そして、真位相分布決定部３３４は、正誤判定結果に基づいて真の位相分布を決定する（ステップＳ１１０４）。以下、各処理の詳細を説明する。
【００５７】
まず、ステップＳ１１０１のシード点候補特定部３３１によるシード点候補特定処理を説明する。シード点候補特定部３３１は、２つのエコー時間の原画像を用いて、シード点候補画素群を特定する。このとき、シード点として任意の画素を選択する際、水が主成分の画素が選択される確率が高まるよう、シード点候補画素群を決定する。本実施形態では、ノイズの影響を強く受けた画素を、水が主成分の画素である確率が高いシード点候補画素群として特定する。
【００５８】
これを実現するため、本実施形態のシード点候補特定部３３１は、図８に示すように、異なる２つのエコー時間の原画像から比画像を算出し、その絶対値を算出する絶対値比画像算出部４１１と、絶対値比画像からシード点候補の画素群を特定するための閾値を設定する閾値設定部４１２と、閾値を用いて絶対値比画像におけるシード点候補画素群を特定する画素群特定部４１３と、を備える。
【００５９】
このシード点候補特定部３３１の各部によるシード点候補画素群特定処理の流れの概略を説明する。図９は、本実施形態のシード点候補画素群特定処理の流れを説明するためのフローチャートである。
【００６０】
まず、絶対値比画像算出部４１１は、エコー時間ｔoutで計測したエコー信号から得た原画像Ｉoutを、エコー時間ｔinで計測したエコー信号から得た原画像Ｉinで複素除算し、比画像Ｉを算出する（ステップＳ１２０１）。そして、その絶対値をとり、絶対値比画像｜Ｉ｜を算出する（ステップＳ１２０２）。なお、絶対値比画像｜Ｉ｜は、以下の式（３）で算出される。
｜Ｉ｜＝｜Ｉout／Ｉin｜・・・（３）
【００６１】
次に、閾値設定部４１２は、少なくとも一つの閾値Ｔｈを設定する（ステップＳ１２０３）。画素群特定部４１３は、絶対値比画像｜Ｉ｜の画素値のうち、閾値Ｔｈ以上の画素値を持つ画素群を、シード点候補画素群として特定する（ステップＳ１２０４）。
【００６２】
ここで、閾値設定部４１２によるステップＳ１２０３の閾値Ｔｈの設定手法について説明する。閾値Ｔｈは、シード点候補画素群として、原画像または比画像において、水が主成分の画素の比率が高い画素群が特定されるよう設定される。
【００６３】
原画像Ｉinにおいて、水信号が１００％の画素（水画素）の信号値（単に水信号と呼ぶ）をＳwin、脂肪信号が１００％の画素（脂肪画素）の信号値（単に脂肪信号と呼ぶ）をＳfin、原画像Ｉoutにおける水信号をＳwout、脂肪信号をＳfout、ガウス分布に従う最頻値０のノイズの標準偏差をσとする。このとき、水信号ＳwoutをＳwinで複素除算したときの信号比（水信号比）をＳmw、脂肪信号ＳfoutをＳfinで複素除算したときの信号比（脂肪信号比）をＳmfとすると、水信号比Ｓmwおよび脂肪信号比Ｓmfはそれぞれ以下の式（４−１）、式（４−２）で表される。
Ｓmw＝Ｓwout／Ｓwin・・・（４−１）
Ｓmf＝Ｓfout／Ｓfin・・・（４−２）
なお、水信号比Ｓmwおよび脂肪信号比Ｓmfは、それぞれ、比画像における水画素および脂肪画素の信号値に相当する。
【００６４】
このとき、水信号比Ｓmwおよび脂肪信号比Ｓmfの標準偏差をそれぞれσmwおよびσmfとすると、誤差伝播の法則より、σmwおよびσmfはそれぞれ、以下の式（５−１）、式（５−２）で表される。
【数５】
2014090949.tif 000005
【００６５】
エコー時間ｔinおよびｔoutで計測したときの信号強度が変化しないと仮定して、Ｓwin＝Ｓwout＝Ｓw、Ｓfin＝Ｓfout＝Ｓfとすると、式（５−１）および式（５−２）は、以下の式（６−１）および式（６−２）で表される。
【数６】
2014090949.tif 000006
なお、ＳＮＲwおよびＳＮＲfは、水信号Ｓwと脂肪信号Ｓfの信号対雑音比（ＳＮＲ：Ｓｉｇｎａｌ  ｔｏ  Ｎｏｉｓｅ  Ｒａｔｉｏ）をそれぞれ表し、以下の式（７−１）および式（７−２）で定義される。
ＳＮＲw＝Ｓw／σ・・・（７−１）
ＳＮＲf＝Ｓf／σ・・・（７−２）
【００６６】
通常、２-Ｐｏｉｎｔ  Ｄｉｘｏｎ法における原画像ＩinおよびＩoutは、グラディエントエコー法もしくはスピンエコー法を用いたＴ1強調画像、または、高速スピンエコー法を用いたＴ2強調画像として得る。このとき、脂肪信号Ｓfは、水信号Ｓwより常に大きくなることが知られている。従って、脂肪信号ＳfのＳＮＲfは、水信号ＳwのＳＮＲwよりも高くなる。すなわち、脂肪信号ＳfのＳＮＲfと、水信号ＳwのＳＮＲwとは、以下の式（８）の関係がある。
ＳＮＲw＜ＳＮＲf・・・（８）
したがって、式（６−１）および式（６−２）より、水信号比Ｓmwの標準偏差σmwは脂肪信号比Ｓmfの標準偏差σmfよりも大きくなる。すなわち、標準偏差σmwと標準偏差σmfとは、以下の式（９）の関係がある。
σmw＞σmf・・・（９）
【００６７】
たとえば、水信号ＳwのＳＮＲwを１０、脂肪信号ＳfのＳＮＲfをＳＮＲwの２倍の２０とすると、式（６−１）および式（６−２）より、脂肪信号比Ｓmfの標準偏差σmfと水信号比Ｓwfの標準偏差σmwとの関係は、以下の式（１０）となる。
σmf＝σmw／２・・・（１０）
【００６８】
このときの、水信号比Ｓmwおよび脂肪信号比Ｓmfの、信号比（比画像Ｉの画素値）に応じた画素数（ボクセル数）のヒストグラム（信号比の分布）を図１０に示す。図１０の横軸は信号比、縦軸は画素数（ボクセル数）をそれぞれ表す。また、６０１は、水信号比Ｓmwの分布、６０２は、脂肪信号比Ｓmfの分布をそれぞれ表す。なお、ここでは、水信号Ｓwと脂肪信号Ｓfの画素数をそれぞれ１０4個とした。また、このときの水信号比Ｓmwおよび脂肪信号比Ｓmfの最頻値は、１である。
【００６９】
図１０に示すように、ＳＮＲwがＳＮＲfより小さいとき、水信号比Ｓmwの標準偏差σmwは脂肪信号比Ｓmfの標準偏差σmfよりも大きくなるため、水信号比Ｓmwのヒストグラムは、脂肪信号比Ｓmfのヒストグラよりも広い分布を持つ。従って、分布の両側の裾野部分の信号比においては、水信号比Ｓmwの画素数の方が、脂肪信号比Ｓmfの画素数より多いことが分かる。なお、信号比分布の裾野部分は、ノイズの影響を強く受けた画素群である。
【００７０】
本実施形態では、このような水信号比Ｓmwおよび脂肪信号比Ｓmfの分布の特性を利用し、水が主成分である画素の比率の高い画素群が特定される閾値Ｔｈを設定する。すなわち、絶対値比画像の画素値の分布を用いて、閾値Ｔｈを設定する。例えば、図１０に示すように、閾値Ｔｈ６０３を、信号比１．２に設定すると、信号比が１．２以上の画素は、水画素である確率が高い。
【００７１】
なお、ここでは、水または脂肪が画素内に１００％存在する場合を仮定した。しかし、実際の生体組織では、組織によって水と脂肪とが異なる割合で混在するため、エコー時間ｔinおよびｔoutで計測したときの信号強度が変化しないという仮定は、多くの場合成り立たない。通常はＳwin≧Ｓwout、Ｓfin≧Ｓfoutとなり、水信号比Ｓmwおよび脂肪信号比Ｓmfの最頻値は、１よりも小さい値をとる。従って、これを考慮して閾値Ｔｈを設定する。例えば、閾値Ｔｈとして、１．０以上に設定する。なお、組織によって水と脂肪とが異なる割合で混在する場合であっても、水信号比Ｓmwの画素数と脂肪信号比Ｓmfの画素数との関係は、上述と同様である。
【００７２】
例えば、閾値Ｔｈは、予め定めて外部記憶装置１１１あるいは計算機１０９の記憶装置に保持しておくよう構成してもよい。また、閾値設定部４１２は、計測毎に、絶対値比画像｜Ｉ｜について、信号比の分布（ヒストグラム）の最頻値および標準偏差を算出し、その最頻値と標準偏差の和、または、それ以上の任意の値を、閾値Ｔｈと設定してもよい。また、最頻値より大きい任意の値を閾値Ｔｈと設定してもよい。また、水信号ＳwのＳＮＲwと脂肪信号ＳfのＳＮＲfによって、信号比の分布の標準偏差は変化する。従って、本実施形態の閾値設定部４１２は、原画像ＩinまたはＩoutのＳＮＲに基づいて、ガウス分布を模擬した計算機シミュレーションを行い、信号比ヒストグラムを算出し、当該モデルに基づいて閾値Ｔｈを算出し、設定してもよい。
【００７３】
次に、上記ステップＳ１１０２の、仮位相分布算出部３３２による仮位相算出処理について説明する。本実施形態の仮位相分布算出部３３２は、シード点候補画素群から任意のシード点を選択し、比画像Ｉの位相値に対し、選択したシード点から領域拡大法により位相アンラップ処理を行い、各画素の位相値を決定し、仮の位相分布を算出する。これを実現するため、本実施形態の仮位相分布算出部３３２は、図８に示すように、シード点を基準画素として領域拡大法を用いた位相アンラップ処理を実行し、各画素の位相値を算出する位相アンラップ部４２１を備える。
【００７４】
位相アンラップ処理では、隣接する画素間での、水と脂肪との位相差による±πの位相とび、および、折り返しによる２πを超える位相とびを判定する。そして、それぞれ、シード点に隣接する画素の位相に、＋πあるいはーπを加算、＋２πあるいは−２πを加算し、滑らかに補正する。本実施形態の位相アンラップ部４２１は、シード点と隣接する画素との間の位相差Δｐを算出し、位相差Δｐの値に応じ、＋π、−π、＋２πあるいは−２πを、隣接する画素の位相に加算することにより、位相アンラップ処理を行う。
【００７５】
なお、本実施形態の位相アンラップ処理は上述の方法に限らない。既に公知の方法で位相をアンラップしてもよい。
【００７６】
上記ステップＳ１１０２の仮位相分布算出部３３２による仮の位相分布算出処理の流れの概略を説明する。図１１は、本実施形態の仮位相分布算出処理の流れを説明するためのフローチャートである。
【００７７】
まず、仮位相分布算出部３３２は、比画像Ｉにおいて、上述のシード点候補画素群の任意の１画素を、シード点と設定する（ステップＳ１３０１）。ここでは、シード点候補画素群の中からランダムに選択した画素を、シード点として設定する。なお、設定するシード点は、シード点候補画素群のうち、絶対比画像｜Ｉ｜の値が最大となる画素としてもよい。
【００７８】
次に、位相アンラップ部４２１は、シード点を基準画素として位相アンラップ処理を行い、仮の位相分布を算出する（ステップＳ１３０２）。具体的には、シード点の位相と、シード点と隣り合う画素の位相との位相差を算出する。そして、この位相差が小さくなるように上述の位相アンラップ処理を行う。隣接する画素が複数ある場合、全ての隣接する画素について、位相アンラップ処理を行う。
【００７９】
次に、仮位相分布算出部３３２は、位相アンラップ処理後の画素を新たなシード点として設定し（ステップＳ１３０３）、画像内の全画素の処理が終了するまで、処理対象の領域を拡大し、位相アンラップ処理を繰り返し（ステップＳ１３０４）、仮位相分布を得る。
【００８０】
なお、本実施形態では、仮の位相分布を、シード点候補画素群の中からシード点を選択し、領域拡大法を用いて位相アンラップ処理を行うことにより算出したが、シード点を用いずに位相分布を算出する公知の方法を用いて、仮の位相分布を算出してもよい。
【００８１】
次に、上記ステップＳ１１０３の正誤判定部３３３による正誤判定処理の詳細について説明する。本実施形態の正誤判定部３３３は、仮位相分布算出部３３２が算出した仮位相分布が、正しい位相分布であるか、位相値がπだけずれた誤った位相分布であるかを判定する。判定は、比画像Ｉの、シード点候補特定部３３１が特定したシード点候補画素群の信号を用いて行う。
【００８２】
シード点候補画素群の任意の画素の位置をｒ、位置ｒにおける水の信号強度をＷ（ｒ）、脂肪の信号強度をＦ（ｒ）、エコー時間ｔinとｔoutとの時間差ｄｔの間に生じた静磁場不均一による位相変化をφ（ｒ）とすると、比画像の信号Ｉ（ｒ）は、以下の式（１１）で表される。
【数１１】
2014090949.tif 000007
仮の位相分布における位置ｒの画素の位相値をθ（ｒ）、φ（ｒ）とθ（ｒ）との位相差をδ（ｒ）、比画像の信号Ｉ（ｒ）の位相からθ（ｒ）を除去した後の比画像の信号（除去後信号値）をＩ’（ｒ）とすると、Ｉ’（ｒ）は、以下の式（１２）で表される。
【数１２】
2014090949.tif 000008
仮の位相分布θが正しく算出されているとき、すなわち、仮の位相分布の正誤判定が真のときは、δ（ｒ）は０となる。従って、式（１２）は、以下の式（１３）で表される。
【数１３】
2014090949.tif 000009
一方、仮の位相分布が誤って算出されているとき、すなわち、仮の位相分布の正誤判定が偽のときは、水と脂肪の位相差によってδ（ｒ）はπとなる。従って、式（１２）は、以下の式（１４）で表される。
【数１４】
2014090949.tif 000010
【００８３】
したがって、式（１３）および式（１４）より、Ｗ（ｒ）＞Ｆ（ｒ）のとき、すなわち、位置ｒの画素の信号値の主成分が水のとき、Ｉ'true（ｒ）＞０、Ｉ’false（ｒ）＜０となる。また、Ｗ（ｒ）＜Ｆ（ｒ）のとき、すなわち、位置ｒの画素の信号値の主成分が脂肪信号のとき、Ｉ’true（ｒ）＜０、Ｉ’false（ｒ）＞０となる。なお、本実施形態では、Ｗ（ｒ）＝Ｆ（ｒ）のときは演算対象から除外する。
【００８４】
上述したとおり、シード点候補画素群に含まれる水画素の数は、脂肪画素の数よりも多い。従って、仮の位相分布が正しく算出されていれば、Ｉ’（ｒ）＞０を満たす画素数が、Ｉ’（ｒ）＜０を満たす画素数より多くなるはずである。本実施形態の正誤判定部３３３は、この特性を利用し、仮の位相分布の正誤を判定する。
【００８５】
正誤判定部３３３は、比画像Ｉのシード点候補画素群の中で、Ｉ’（ｒ）＞０を満たす画素数Ｎpと、Ｉ’（ｒ）＜０を満たす画素数Ｎnと、を算出する。そして、Ｎp＞Ｎnのとき、正誤判定結果は真、すなわち、仮の位相分布が正しく算出されているとする。一方、Ｎp≦Ｎnのとき、正誤判定結果は偽、すなわち、仮の位相分布が誤って出されているとする。
【００８６】
上記の手順により正誤判定を行うため、本実施形態の正誤判定部３３３は、図８に示すように、比画像Ｉのシード点候補画素群の各画素について、その信号の位相値から、仮位相分布のとして算出された位相値を除去する位相除去部４３１と、比画像Ｉのシード点候補画素群の各画素について、位相除去後の信号値（除去後信号値）が正の値である画素数Ｎpと、負の値である画素数Ｎnとを算出する画素数算出部４３２と、正の値である画素数Ｎpと負の値である画素数Ｎnとの大小を比較し、仮位相分布の正誤（真偽）を判定する判定部４３３と、を備える。
【００８７】
真位相分布決定部３３４は、正誤判定部３３３の判定結果が正（真）のときは、仮の位相分布を真の位相分布として決定する。また、正誤判定結果が誤（偽）のときは、仮の位相分布に位相値πを加算したものを、真の位相分布として決定する。
【００８８】
仮位相分布の正誤判定を行い、真の位相分布を決定する処理（真位相分布決定処理）の流れを、図１２のフローチャートを用いて説明する。図１２は、本実施形態の真位相分布決定処理のフローチャートである。
【００８９】
まず、位相除去部４３１は、比画像Ｉのシード点候補画素群の信号（比画像信号Ｉ）から、仮の位相分布の位相値を除去して（ステップＳ１４０１）、比画像信号Ｉ’を算出する。これにより、静磁場不均一の影響による位相変化が取り除かれる。
【００９０】
次に、画素数算出部４３２は、シード点候補画素群のうち、I’＞０を満たす画素の数Ｎpと、I’＜０を満たす画素の数Ｎnとを算出する（ステップＳ１４０２）。なお、I’＝０を満たす画素は演算対象から除外する。その後、判定部４３３は、画素数ＮpとＮnとの大きさを比較する（ステップＳ１４０３）。
【００９１】
判定結果がＮp＞Ｎnのとき、真位相分布決定部３３４は、仮の位相分布を真の位相分布として決定する（ステップＳ１４０５）。一方、判定結果が、Ｎp≦Ｎnのとき、真位相分布決定部３３４は、仮の位相分布に位相値πを加算し（ステップＳ１４０４）、加算後の位相分布を、真の位相分布として決定する（ステップＳ１４０５）。
【００９２】
以上の手順により、静磁場不均一に起因した位相分布を正確に算出することができる。すなわち、仮に、シード点が、脂肪が主成分の画素に設定され、算出した位相分布に脂肪の化学シフトによる位相値πのオフセットが加算されている場合であっても、それを判定することができ、オフセットを除いた位相分布を得ることができる。
【００９３】
次に、上記ステップＳ１００４の、水−脂肪分離部３４０による水−脂肪画像算出方法について説明する。水−脂肪分離部３４０は、上記ステップＳ１００２で得たエコー時間の異なる２つの原画像ＩinおよびＩoutと、上記ステップＳ１００３で算出した真の位相分布と、を用いて水画像と脂肪画像を算出する。具体的には、原画像Ｉoutの各画素の位相値を真の位相分布の各画素の位相値でそれぞれ除算（ｅｘｐ（−ｉθ））した後、または、原画像Ｉinの各画素の位相値に真の位相分布の各画素の位相値をそれぞれ乗算（ｅｘｐ（＋ｉθ））した後、式（２−１）および式（２−２）により、水画像と脂肪画像とを算出する。
【００９４】
以上説明したように、本実施形態のＭＲＩ装置１００は、被検体１０１が置かれる空間に静磁場を生成する静磁生成磁石１０２と、前記被検体１０１に高周波磁場パルスを送信する送信部と、高周波磁場パルスの送信により前記被検体から生じる核磁気共鳴信号を受信する受信部と、前記核磁気共鳴信号に位置情報を付加するための傾斜磁場を印加する傾斜磁場印加部と、前記送信部、前記傾斜磁場印加部及び前記受信部の動作を制御するとともに前記受信した核磁気共鳴信号に対して演算処理を行う計算機１０９と、を備えるＭＲＩ装置１００であって、前記計算機１０９は、予め定めたパルスシーケンスに従って、異なる２つのエコー時間の前記核磁気共鳴信号を得る計測部３１０と、前記２つの核磁気共鳴信号から、前記異なる２つのエコー時間の原画像を再構成する再構成部３２０と、前記異なる２つのエコー時間の原画像から静磁場不均一によって変化する位相の分布である位相分布を算出する位相分布算出部３３０と、前記異なる２つのエコー時間の原画像と前記位相分布とを用いて水信号と脂肪信号とを分離する水脂肪分離部３４０と、を備え、前記位相分布算出部３３０は、前記異なる２つのエコー時間の原画像を用い、ノイズの影響を受けた画素群をシード点候補画素群として特定するシード点候補特定部３３１と、前記異なる２つのエコー時間の原画像と前記シード点候補画素群の情報とを用い、仮の前記位相分布である仮位相分布を算出する仮位相分布算出部３３２と、前記仮位相分布の正誤を判定する正誤判定部３３３と、前記正誤判定部３３３による判定結果に基づいて真の位相分布を決定する真位相分布決定部３３４と、を備える。
【００９５】
なお、前記シード点候補特定部３３１は、前記異なる２つのエコー時間の原画像の比として得た比画像の絶対値を絶対値比画像として算出する絶対値比画像算出部４１１と、前記シード点候補画素群を特定するための閾値を少なくとも１つ設定する閾値設定部４１２と、前記閾値を用いて前記絶対値比画像において前記シード点候補画素群を特定する画素群特定部４１３と、を備えてもよい。
【００９６】
前記閾値設定部４１２は、前記絶対値比画像の画素値の分布を用いて前記閾値を設定してもよい。また、前記閾値設定部４１２は、前記異なる２つのエコー時間の原画像の信号雑音比に基づいて、信号比ヒストグラムモデルを算出し、当該モデルに基づいて前記閾値を設定してもよい。また、前記閾値は、前記絶対値比画像の画素値の最頻値に当該画素値の標準偏差を加算した値以上の値に設定されてもよい。また、前記閾値は、前記絶対値比画像の画素値の最頻値より大きい値に設定されてもよい。また、前記閾値は、１．０以上の値に設定されてもよい。
【００９７】
前記仮位相分布算出部３３２は、前記比画像の位相値に対し、予め定めた基準画素から領域拡大法を用いた位相アンラップ処理を実行し、前記仮位相分布を算出する位相アンラップ部４２１を備え、前記基準画素は、前記シード点候補画素群の中から選択されてもよい。
【００９８】
前記正誤判定部３３３は、前記比画像の前記シード点候補画素群の各信号から前記仮位相分布として得た位相値をそれぞれ除去し、前記各シード点候補画素群の各画素の除去後信号値を得る位相除去部４３１と、前記各シード点候補画素群の中で、前記除去後信号値が正の値である画素数と負の値である画素数とをそれぞれ算出する画素数算出部４３２と、前記正の値の画素数と前記負の値の画素数とを比較し、前記仮位相分布の正誤を判定する判定部４３３と、を備えてもよい。そして、前記判定部４３３は、前記正の値の画素数が前記負の値の画素数より多い場合、正と判定してもよい。
さらに、前記真位相分布決定部３３４は、前記判定部４３３の判定結果が誤の場合、前記仮位相分布を補正し、補正後の位相分布を、前記真の位相分布と決定し、当該判定結果被核結果が正の場合、仮の位相分布を前記真の位相分布と決定してもよい。
【００９９】
また、本実施形態の水脂肪分離方法は、予め定めたパルスシーケンスに従って計測した異なる２つのエコー時間の核磁気共鳴信号からそれぞれ再構成された、前記異なる２つのエコー時間の原画像から、水信号と脂肪信号とを分離する水脂肪分離方法であって、前記異なる２つのエコー時間の原画像を用い、ノイズの影響を受けた画素群をシード点候補画素群として特定し、前記異なる２つのエコー時間の原画像と前記シード点候補画素群の情報とを用い、仮の前記位相分布である仮位相分布を算出し、前記仮位相分布の正誤を判定し、前記正誤判定部による判定結果に基づいて真の位相分布を決定し、前記異なる２つのエコー時間の原画像と前記真の位相分布とを用いて、水信号と脂肪信号とを分離することを特徴とする。
【０１００】
このように、本実施形態によれば、静磁場不均一に起因した位相変化である位相分布を、異なる２つのエコー時間の原画像の比画像の位相値に対して領域拡大法を用いた位相アンラップ処理を行うことにより算出する。算出に先立ち、位相アンラップ処理の開始点であるシード点が、水信号が主である画素となる確率が高まるよう、シード点候補画素群を特定する。従って、位相アンラップ処理により算出される位相分布に、位相値πのオフセットが加算される可能性が低減する。すなわち、正確な位相分布が算出される可能性が高まる。
【０１０１】
さらに、算出された位相分布の正誤判定を行う。このとき、水信号が主である画素となる確率の高い、シード点候補画素群を用いて、画素数により判定を行う。従って、簡易に高精度な判定を行うことができる。高精度な判定結果に基づいて、最終的に位相分布を決定するため、正確な位相分布が算出される可能性は、さらに高まる。そして、本実施形態によれば、この正確な位相分布である可能性が極めて高い位相分布を用いて水信号と脂肪信号とを分離し、水画像、脂肪画像とを得る。
【０１０２】
従って、本実施形態によれば、得られる水画像、脂肪画像が、正しいものである可能性も極めて高いものとなる。すなわち、本実施形態によれば、２-Ｐｏｉｎｔ  Ｄｉｘｏｎ法を用いた水と脂肪の分離撮影において、静磁場不均一が反映された位相画像を正しく算出する確度を高めることで、水画像と脂肪画像とが入れ替わることを防ぎ、操作者に判断する負担を負わせることなく、高精度に水と脂肪とを分離できる。
【０１０３】
なお、上記実施形態では、エコー信号を取得する異なる２つのエコー時間として、水と脂肪とが同位相となるエコー時間ｔinと、水と脂肪とが逆位相となるエコー時間ｔoutとを用いる場合を例にあげて説明したが、エコー信号の取得時間は、これに限られない。例えば、２つのエコー時間が、以下の式（１５）を満たすＴＥ1およびＴＥ2であってもよい。
ｃｏｓ（２πｆFWＴＥ1）≠ｃｏｓ（２πｆFWＴＥ2）・・・（１５）
ｆFWは、水と脂肪との周波数差である。
【０１０４】
この場合、上記Ｓ１１０２の仮の位相分布の算出時に、比画像の代わりに、δW、δFを用いる。なお、δWは、生体組織の信号が全て水と仮定した場合に算出される位相分布であり、δFは、生体組織の信号が全て脂肪と仮定した場合に算出される位相分布である。そして、拡大領域法により位相アンラップ処理を行う際、シード点に隣接する画素の位相として、当該画素のδWの位相とδFの位相とのうち、シード点の画素の位相に近い方を採用する。
【符号の説明】
【０１０５】
１００：ＭＲＩ装置、１０１：被検体、１０２：静磁場コイル、１０３：傾斜磁場コイル、１０４：シムコイル、１０５：送信コイル、１０６：受信コイル、１０７：送信機、１０８：受信機、１０９：計算機、１１０：ディスプレイ、１１１：外部記憶装置、１１２：傾斜磁場用電源部、１１３：シム用電源部、１１４：シーケンス制御装置、１１５：入力装置、１２０：ＭＲＩ装置、１３０：ＭＲＩ装置、３１０：計測部、３２０：再構成部、３３０：位相分布算出部、３３１：シード点候補特定部、３３２：仮位相分布算出部、３３３：正誤判定部、３３４：真位相分布決定部、３４０：水−脂肪分離部、４１１：絶対値比画像算出部、４１２：閾値設定部、４１３：画素群特定部、４２１：位相アンラップ部、４３１：位相除去部、４３２：画素数算出部、４３３：判定部、５１０：パルスシーケンス、５１１：ＲＦパルス、５１２：スライス選択傾斜磁場、５１３：位相エンコード傾斜磁場、５１４：読み出し傾斜磁場、５１５：第一エコー、５１６：読み出し傾斜磁場、５１７：第二エコー、５２０：パルスシーケンス、５２１：ＲＦパルス、５２２：スライス選択傾斜磁場、５２３：位相エンコード傾斜磁場、５２４：読み出し傾斜磁場、５２５：第一エコー、５２６：読み出し傾斜磁場、５２７：第二エコー、５２８：ＲＦパルス、５２９：スライス選択傾斜磁場、６０１：水信号比分布、６０２：脂肪信号比分布、６０３：閾値Ｔｈ
【図１】
2014090949.tif 000011
【図２】
2014090949.tif 000012
【図３】
2014090949.tif 000013
【図４】
2014090949.tif 000014
【図５】
2014090949.tif 000015
【図６】
2014090949.tif 000016
【図７】
2014090949.tif 000017
【図８】
2014090949.tif 000018
【図９】
2014090949.tif 000019
【図１０】
2014090949.tif 000020
【図１１】
2014090949.tif 000021
【図１２】
2014090949.tif 000022