【手続補正書】
【提出日】平成30年11月12日(2018.11.12)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】請求項２
【補正方法】変更
【補正の内容】
【請求項２】
請求項１に記載の磁気共鳴イメージング装置であって、
前記位相分布分離部は、
複数の前記低解像度画像について、画像間の位相差を算出する位相差算出部と、
前記位相差算出部が算出した位相差の位相折り返しを除去する折り返し除去部と、
前記折り返し除去部により折り返しが除去された位相差を、周波数に換算して周波数分布を算出する周波数変換部と、
前記周波数変換部が算出した周波数分布を、大域的周波数分布と局所的周波数分布とに分離する周波数分布分離部と、
前記低解像度画像と前記周波数変換部が算出した周波数分布とを用いて、オフセット位相分布を算出するオフセット位相分布算出部と、
を備えることを特徴とする磁気共鳴イメージング装置。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】００５９
【補正方法】変更
【補正の内容】
【００５９】
［低解像度化：Ｓ１１０１］
低解像度変換部３３１は、計測したマルチエコー複素画像Ｉ（ｎ）を線形補間やＣｕｂｉｃ補間など公知の手法を用いて、計測した複素画像Ｉ（ｎ）の分解能以下となる所望の解像度に変換し、低解像複素画像ｉ（ｎ）を算出する。なお、低解像度変換部３３１は、計測した被検体の構造に併せて、被検体内部の分解能を粗く、その辺縁を細かくするなど、領域によって任意の解像度に変換した低解像複素画像ｉ（ｎ）を算出してもよい。あるいは、複数の解像度の複素画像ｉ（ｎ）を算出してもよい。
【手続補正３】
【補正対象書類名】明細書
【補正対象項目名】００８１
【補正方法】変更
【補正の内容】
【００８１】
［ＳＷＩ法を用いた磁化率強調画像の算出法］
以下に、磁化率強調画像算出部７２０で行われる、ＳＷＩ法による磁化率強調画像の算出法について説明する。
本実施形態では、まず、マスク生成部７２１が、局所的周波数分布Ｆｌｏｃａｌから磁化率を強調する磁化率強調マスクを作成する。具体的には、任意の周波数ｆｍ以下の値を０、周波数ｆｍより多い周波数を１、その間を線形に結んだ値に設定する強調マスクを作成する。
【手続補正４】
【補正対象書類名】明細書
【補正対象項目名】００８２
【補正方法】変更
【補正の内容】
【００８２】
その後、磁化率強調画像作成部７２２で、強調マスクを任意の回数乗算したのち、任意のＴＥの複素画像の絶対値成分（絶対値画像）に乗算することで、磁化率強調画像を算出する。
なお、本実施形態の磁化率強調画像の算出法は上記の方法に限られない。公知のさまざまな手法を適用してもよい。
【手続補正５】
【補正対象書類名】明細書
【補正対象項目名】００８３
【補正方法】変更
【補正の内容】
【００８３】
［画像表示：Ｓ１００５］
分布画像算出部３４０が算出した定量的磁化率分布や磁化率強調画像は、ディスプレイ１１０（図２）に表示することができる。或いは外部記憶装置１１１に画像データとして格納し、所望の表示装置で表示されてもよい。
【手続補正６】
【補正対象書類名】明細書
【補正対象項目名】００８８
【補正方法】変更
【補正の内容】
【００８８】
また必須ではないが、本実施形態の計算機は、低解像オフセット位相分布及び／又は低解像大域的周波数分布を高解像度に変換する際に、平滑化する手段や辺縁のデータを外挿する手段を備えることができる。
【手続補正７】
【補正対象書類名】明細書
【補正対象項目名】００９４
【補正方法】変更
【補正の内容】
【００９４】
まず、位相分布分離部３３２は、低解像度変換部３３１にて複素画像Ｉ（ｎ）を低解像度化した低解像複素画像ｉ（ｎ）において、隣り合うエコー間の位相差分布Δｐ（ｎ）を（Ｎ−１）セット算出する（Ｓ１２０１）。そして、算出した位相差分布Δｐ（ｎ）を加重平均して一つの位相差分布Δｐを算出する（Ｓ１２０２）。位相差分布Δｐに対して位相アンラップ処理を実施し（Ｓ１２０３）、周波数成分に換算し静磁場不均一周波数分布ｆｉｎｈｏｍｏを算出する（Ｓ１２０４）。周波数分布分離による大域的周波数分布算出処理により、静磁場不均一周波数分布ｆｉｎｈｏｍｏから大域的周波数分布ｆｇｌｏｂａｌを算出する（Ｓ１２０５）。低解像複素画像ｉ（ＴＥｎ）から、静磁場不均一周波数分布ｆｉｎｈｏｍｏによる位相変化を除去することで、各エコーのオフセット位相分布ｐｏｆｆｓｅｔ（ｎ）を算出する（Ｓ１２０６）。各エコーのオフセット位相分布ｐｏｆｆｓｅｔ（ｎ）を加重平均することでオフセット位相分布ｐｏｆｆｓｅｔを算出する（Ｓ１２０７）。以上の手順により、信号フィッティングを用いずに低解像大域的周波数分布ｆｇｌｏｂａｌと低解像オフセット位相分布ｐｏｆｆｓｅｔを算出することができる。
【手続補正８】
【補正対象書類名】明細書
【補正対象項目名】００９７
【補正方法】変更
【補正の内容】
【００９７】
［位相差分布算出：Ｓ１２０１］
位相差算出部９０１において、低解像マルチエコー複素画像ｉ（ｎ）における隣り合うエコー間の位相差分布Δｐ（ｎ）を（Ｎ−１）セット算出する。ここで、（ｎ＋１）番目のエコー時間ＴＥｎ＋１における低解像複素画像ｉ（ｎ＋１）と、ｎ番目のエコー時間ＴＥｎにおける低解像複素画像ｉ（ｎ）は、低解像度の静磁場不均一周波数分布ｆｉｎｈｏｍｏと、低解像オフセット位相ｐｏｆｆｓｅｔを用いて、それぞれ式（１４）と式（１５）で表される。
【数１４】
2017184935.tif 000001
【数１５】
2017184935.tif 000002
【手続補正９】
【補正対象書類名】明細書
【補正対象項目名】０１１５
【補正方法】変更
【補正の内容】
【０１１５】
＜第三の実施形態＞
上述した第二の実施形態では、図８及び図９に示したように、位相分布分離部３３２は、まず低解像複素画像ｉ（ｎ）の位相差を求め、オフセット位相分布を除去してから位相アンラップ及び周波数変換を行って低解像大域的周波数分布を算出し、さらに低解像オフセット位相分布を算出した。
【手続補正１０】
【補正対象書類名】明細書
【補正対象項目名】０１２２
【補正方法】変更
【補正の内容】
【０１２２】
［フィッティング：Ｓ１３０１］
信号フィッティング部９１１において、低解像マルチエコー複素画像ｉ（ｎ）をＧｒＥ法の信号モデルに信号フィッティングさせる。ＧｒＥ法で計測した画素内の信号モデルＳ（ｎ）は、画素内のプロトン密度分布をＭ０、見かけの横緩和速度分布をＲ２＊とすると、式（２５）で表される。
【数２５】
2017184935.tif 000003
【手続補正１１】
【補正対象書類名】明細書
【補正対象項目名】０１３３
【補正方法】変更
【補正の内容】
【０１３３】
位相分布分離部３３２は、低解像度変換部３３１にて複素画像Ｉ（ｎ）を低解像度化した低解像複素画像ｉ（ｎ）において、ＧｒＥ法の信号モデルを用いたフィッティングにより、水の信号強度Ｗと脂肪の信号強度Ｆ、静磁場不均一周波数分布ｆｉｎｈｏｍｏとゼロ次位相分布を算出する（Ｓ１４０１）。そして、信号フィッティングで算出したゼロ次位相成分をオフセット位相分布ｐｏｆｆｓｅｔとして設定する（Ｓ１４０２）。オフセット位相分布ｐｏｆｆｓｅｔと、水と脂肪の周波数差Ｆｆａｔに起因してＴＥｎによって変化する位相成分pｆａｔ（ｎ）を合成して合成オフセット位相分布ｐａｄｄ（ｎ）を算出する（Ｓ１４０３）。そして、信号フィッティングで算出した静磁場不均一周波数を周波数アンラップして静磁場不均一周波数分布ｆｉｎｈｏｍｏを算出する（Ｓ１４０４）。大域的周波数算出処理により、周波数アンラップした静磁場不均一周波数分布ｆｉｎｈｏｍｏから大域的周波数分布ｆｇｌｏｂａｌを算出する（Ｓ１４０５）。以上の手順により、低解像大域的周波数分布ｆｇｌｏｂａｌと、低解像オフセット位相分布ｐｏｆｆｓｅｔ（水と脂肪の周波数差に起因する位相差pｆａｔを含む）をそれぞれ算出することができる。
【手続補正１２】
【補正対象書類名】明細書
【補正対象項目名】０１４０
【補正方法】変更
【補正の内容】
【０１４０】
［周波数分布分離：Ｓ１４０５］
次に、大域的周波数分布算出部（周波数分布分離部）９２５において、静磁場不均一周波数分布ｆｉｎｈｏｍｏから、生体形状などに起因して生じる局所的周波数変化と大域的周波数変化とを分離し、大域的周波数分布ｆｇｌｏｂａｌを算出する。この算出手法は、第二の実施形態の大域的周波数分布算出部（周波数分布分離部）９０４の処理（図８：Ｓ１２０５）と同様であり、例えば、ＳＭＶフィルタ処理や、その他の手法が採用できる。
【手続補正１３】
【補正対象書類名】明細書
【補正対象項目名】０１４２
【補正方法】変更
【補正の内容】
【０１４２】
その後、図４に示したように、図１２に示す位相分布分離処理において算出した低解像大域的周波数分布と低解像合成オフセット位相分布とをもとの解像度に変換した後（図６：Ｓ１１０３）、これらを計測した複素画像Ｉ（ｎ）から除去し（図６：Ｓ１１０４）、局所的周波数分布ｆｌｏｃａｌを算出する。さらに局所的周波数分布を用いて、磁場分布や、定量的磁化率分布、磁化率強調画像などの分布画像を取得することができる（図４：Ｓ１００４）。
【手続補正１４】
【補正対象書類名】明細書
【補正対象項目名】０１４４
【補正方法】変更
【補正の内容】
【０１４４】
１００：ＭＲＩ装置、１０１：被検体、１０２：静磁場コイル、１０３：傾斜磁場コイル、１０４：シムコイル、１０５：送信コイル、１０６：受信コイル、１０７：送信機、１０８：受信機、１０９：計算機、１１０：ディスプレイ、１１１：外部記憶装置、１１２：傾斜磁場用電源部、１１３：シム用電源部、１１４：シーケンス制御装置、１１５：入力装置、１２０：ＭＲＩ装置、１３０：ＭＲＩ装置、３１０：計測制御部、３２０：画像再構成部、３３０：局所的周波数分布算出部、３３１：低解像度変換部（第一の解像度変換部）、３３２：位相分布分離部、３３３：高解像度変換部（第二の解像度変換部）、３３４：位相除去部、３３５：加重平均部、３４０：分布画像算出部、９０１：位相差算出部、９０２：位相アンラップ部、９０３：周波数変換部、９０４：大域的周波数分布算出部（周波数分布分離部）、９０５：オフセット位相算出部、９１１：信号フィッティング部、９１２：オフセット位相分布設定部、９１３：周波数アンラップ部、９１４：大域的周波数分布算出部（周波数分布分離部）、９２１：信号フィッティング部、９２２：オフセット位相分布設定部、９２３：合成オフセット位相分布算出部、９２４：周波数アンラップ部、９２５：大域的周波数分布算出部（周波数分布分離部）
【手続補正１５】
【補正対象書類名】図面
【補正対象項目名】図３
【補正方法】変更
【補正の内容】
【図３】
2017184935.tif 000004