JP5165791B2 - 磁気共鳴撮影装置 - Google Patents

Description

本発明は、スペクトロスコピックイメージングに関し、特に渦電流などによるスペクトル歪みを補正する渦電流補正の精度を高める技術に関する。

ＭＲＩで計測される核磁気共鳴信号には、分子構造の違いによって、共鳴周波数が僅かに異なる化学シフト現象が生じる。この現象を利用して、分子（代謝物質）ごとに核磁気共鳴信号を分離し、スペクトルを得るＭＲＳ（Ｍａｇｎｅｔｉｃ Ｒｅｓｏｎａｎｃｅ Ｓｐｅｃｔｒｏｓｃｏｐｙ）、代謝物質ごとの空間的な信号強度分布を画像化するＣＳＩ（Ｃｈｅｍｉｃａｌ Ｓｈｉｆｔ Ｉｍａｇｉｎｇ：化学シフトイメージング）、ＭＲＳＩ（Ｍａｇｎｅｔｉｃ Ｒｅｓｏｎａｎｃｅ Ｓｐｅｃｔｒｏｓｃｏｐｉｃ Ｉｍａｇｉｎｇ：核磁気共鳴スペクトロスコピックイメージング）が知られている。

ＭＲＳやＭＲＳＩで検出できる人体の主な代謝物質には、コリン（Ｃｈｏ）、クレアチン（Ｃｒ）、Ｎ−アセチルアスパラギン酸（ＮＡＡ）、乳酸（Ｌａｃ）等がある。これら代謝物質の量から、がん等の代謝異常疾患の進行度判定や早期診断が可能となる。また、腫瘍の悪性度診断を非侵襲的に行うことが可能になると考えられる。

ＭＲＳおよびＭＲＳＩでは、計測時に印加する傾斜磁場によって渦電流が生じる。渦電流は、空間的、時間的に静磁場不均一を生じさせ、計測によって得られたスペクトルの形状を歪ませる。このスペクトル歪みは、一般に参照信号の位相値を用いて補正される。例えば、参照信号として代謝物質よりも信号強度が大きい水信号を用いる。この水信号から、位相値を取得し、補正に使用する（例えば、非特許文献１参照。）。非特許文献１に開示の手法では、この水信号から空間的、時間的な位相値を算出し、代謝物質画像データに対して位相補正を行い、渦電流によるスペクトル歪みを補正している。

ＣＳＩあるいはＭＲＳＩは、計測時間およびＳＮＲ（信号対雑音比）の観点から、計測するマトリクス（ボクセル）数が８×８〜３２×３２程度と非常に少ない。このため、画像再構成で実施するフーリエ変換によって、トランケーションが生じ、遠くのボクセルの信号が混入する。その結果、静磁場不均一が存在する場合、対象とするボクセル内の水信号と異なる周波数の水信号が混入する。

異なる周波数の水信号が混入すると、渦電流補正に用いる位相値の時間変化に、単位時間当たりの位相値の変化量の変動が他の箇所に比べて突出して大きい、いわゆる、位相とび領域が発生する。この位相とび領域における位相変化量の大きさは、混入した水信号の濃度に比例する。ここで、図１３を用いて、この位相とびについて説明する。

図１３（ａ）に、計算機シミュレーションによる、水信号のスペクトルを示す。図１３（ａ）におけるＷａｔｅｒ１およびＷａｔｅｒ２は、周波数がそれぞれ２Ｈｚと５Ｈｚ、濃度比が１．０対０．９の水信号のスペクトルである。この２つの信号が混在する場合の、時間領域における位相値の時間変化を図１３（ｂ）に示す。図１３（ａ）、（ｂ）に示すように、周波数差Δｆに対して、１／Δｆの時間間隔で、濃度比に比例した位相とびが生じる。なお、ここでは、渦電流による静磁場変動を付加し、緩やかな位相変化を生じさせている。

ここで、時間方向に図１３（ｂ）に示す変化を示す位相データを用いて、代謝物質の渦電流によるスペクトル歪みを補正する場合の、補正前後の様子を、それぞれ図１３（ｃ）、（ｄ）に示す。図１３（ｃ）には、渦電流補正前の様子を、図１３（ｄ）は渦電流補正後の様子をそれぞれ示す。図１３（ｄ）に示すように、位相とびのある位相データを用いて補正すると、位相とびによって、リンギングアーチファクトが生じ、渦電流補正処理によって逆にスペクトルが劣化する。

これに対し、例えば、水信号のスペクトルに対してローパスフィルタをかけることによって、リンギングアーチファクトを低減する手法がある（例えば、非特許文献２参照。）。また、位相データに現れる位相とびを補正し、リンギングアーチファクトを低減する手法がある（例えば、非特許文献３参照。）。ここでは、時間領域における水信号の絶対値強度が極値を取るタイミングを位相とびの発生箇所とする。そして、時間領域における水信号の位相データに関して時間ｔについての１次導関数を算出し、上記の位相とび発生箇所の周辺をモデル関数でフィッティングする。このとき、フィッティングした領域を位相とび補正領域とし、上記位相値を線形補間することにより位相とびを補正する。

Ｕｗｅ Ｋｌｏｓｅ "Ｉｎ Ｖｉｖｏ Ｐｒｏｔｏｎ Ｓｐｅｃｔｒｏｓｃｏｐｙ ｉｎ Ｐｒｅｓｅｎｃｅ ｏｆ Ｅｄｄｙ Ｃｕｒｒｅｎｔｓ" Ｍａｇｎｅｔｉｃ Ｒｅｓｏｎａｎｃｅ Ｉｎ Ｍｅｄｉｃｉｎｅ、１４巻、２６−３０頁（１９９０年） Ｊ． Ｍ． Ｗｉｌｄ "Ａｒｔｉｆａｃｔｓ Ｉｎｔｒｏｄｕｃｅｄ ｂｙ Ｚｅｒｏ Ｏｒｄｅｒ Ｐｈａｓｅ Ｃｏｒｒｅｃｔｉｏｎ ｉｎ Ｐｒｏｔｏｎ ＮＭＲ Ｓｐｅｃｔｒｏｓｃｏｐｙ ａｎｄ ａ Ｍｅｔｈｏｄ ｏｆ Ｅｌｉｍｉｎａｔｉｏｎ ｂｙ Ｐｈａｓｅ Ｆｉｌｔｅｒｉｎｇ" Ｊｏｕｒｎａｌ ｏｆ Ｍａｇｎｅｔｉｃ Ｒｅｓｏｎａｎｃｅ、１３７巻、４３０−４３６頁（１９９９年） Ａ． Ｗ． Ｓｉｍｏｎｅｔｔｉ，ｅｔ ａｌ．"Ａｕｔｏｍａｔｅｄ ｃｏｒｒｅｃｔｉｏｎ ｏｆ ｕｎｗａｎｔｅｄ ｐｈａｓｅ ｊｕｍｐｓ ｉｎ ｒｅｆｅｒｅｎｃｅ ｓｉｇｎａｌｓ ｗｈｉｃｈ ｃｏｒｒｕｐｔ ＭＲＳＩ ｓｐｅｃｔｒａ ａｆｔｅｒ ｅｄｄｙ ｃｕｒｒｅｎｔ ｃｏｒｒｅｃｔｉｏｎ" Ｊｏｕｒｎａｌ ｏｆ Ｍａｇｎｅｔｉｃ Ｒｅｓｏｎａｎｃｅ、１５９巻、１５１−１５７頁（２００２年）

しかしながら、非特許文献２に記載の手法では、位相とびの低周波成分が残るため、リンギングアーチファクトを完全には除去できない。また、渦電流の高周波成分をカットしてしまうため、十分な渦電流の補正効果が得られない。

非特許文献３に記載の手法では、時間領域における水信号の絶対値強度に極値が存在しなくても位相とびが発生するときには、位相とびが補正されない。また、位相変化が急峻な箇所が複数ある場合、その中から位相とびの箇所を特定し、抽出することが難しい。さらに、フィッティング領域を指定していないため、位相変化量が異なる複数の位相とびがある場合、フィッティング精度が低下する。すなわち、フィッティング精度を高めるためには、位相とび発生箇所毎にフィッティングを行う必要があり、処理が煩雑となるとともに、自動的に位相とびを補正することが難しい。

本発明は、上記事情に鑑みてなされたものであり、ＭＲＩ装置における渦電流の影響を補正する処理を行うとき、補正処理によるアーチファクトを防ぎ、補正結果の精度を向上させることを目的とする。

本発明は、計測対象の代謝物質より信号強度の大きい物質の信号の位相データを用いて渦電流の影響を補正する渦電流補正処理において、時間領域における位相変化量の大きさにより、位相データ内の位相とびが発生している領域箇所を特定し、位相とびと特定された箇所のみ、位相とび補正を行う。変化量の大きさは、予め定められた領域分割に従って分割された各領域内で決定する。

具体的には、被検体に静磁場を印加する静磁場印加手段と、前記被検体に傾斜磁場を印加する傾斜磁場印加手段と、前記被検体に高周波磁場パルスを照射する高周波磁場パルス照射手段と、前記被検体から核磁気共鳴信号を受信する受信手段と、前記核磁気共鳴信号に渦電流補正を実施する渦電流補正手段と、を備える磁気共鳴撮影装置であって、前記渦電流補正手段は、計測対象の代謝物質より信号強度の大きい物質の核磁気共鳴信号の位相データを算出する位相データ算出手段と、前記位相データを、時間方向に複数の小領域に分割する位相データ分割手段と、前記小領域の中から位相とびが存在する位相とび発生領域を抽出する位相とび発生領域抽出手段と、前記位相とび発生領域内の位相とびを補正する位相とび補正手段と、を備えることを特徴とする磁気共鳴撮影装置を提供する。

本発明によれば、ＭＲＩ装置における渦電流の影響を補正する処理を行うとき、補正処理によるアーチファクトを防ぎ、補正結果の精度を向上させることができる。

（ａ）〜（ｃ）は、本発明の実施形態の磁気共鳴撮影装置の外観図である。 （ａ）は、本発明の実施形態の磁気共鳴撮影装置の機能構成図であり、（ｂ）は、本発明の実施形態の磁気共鳴撮影装置が備える計算機の機能ブロック図である。 本発明の実施形態のＭＲＳＩパルスシーケンスの一例を示す図である。 （ａ）〜（ｃ）は、本発明の実施形態のＭＲＳＩパルスシーケンスにより励起される領域を説明するための図である。 本発明の実施形態の位相折り返し接続処理の一例を説明するための図である。 本発明の実施形態の位相とび補正処理のフローチャートである。 （ａ）は、本発明の実施形態の領域分割処理のフローチャートであり、（ｂ）および（ｃ）は、本発明の実施形態の領域分割処理を説明するための図である。 本発明の実施形態の位相変化量算出処理を説明するための図である。 本発明の実施形態の位相とび補正領域決定処理および位相とび補正処理を説明するための図である。 （ａ）は、本発明の実施形態の位相とび補正処理を実施する前の位相データと、位相とび領域を１箇所補正した後の位相データである。（ｂ）は、本発明の実施形態の位相とび補正処理を実施する前の位相データと、全ての位相とび領域を補正した後の位相データを示す図である。 （ａ）〜（ｃ）は、計算機シミュレーション結果を用いて本発明の実施形態による効果を説明するための図である。 （ａ）〜（ｃ）は、ファントム実験結果を用いて本発明の実施形態による効果を説明するための図である。 （ａ）〜（ｄ）は、水信号の位相データを用いて渦電流補正を行う場合、リンギングアーチファクトが生じる原因を説明するための図である。

以下、本発明を適用する実施形態について説明する。以下、実施形態を説明するための全図において、同一機能を有するものは同一符号を付し、その繰り返しの説明は省略する。

まず、本実施形態の磁気共鳴撮影装置（ＭＲＩ装置）について説明する。図１は、本実施形態のＭＲＩ装置の外観図である。図１（ａ）は、ソレノイドコイルで静磁場を生成するトンネル型磁石を用いた水平磁場方式のＭＲＩ装置１００である。図１（ｂ）は、開放感を高めるために磁石を上下に分離したハンバーガー型（オープン型）の垂直磁場方式のＭＲＩ装置１１０である。また、図１（ｃ）は、図１（ａ）と同じトンネル型磁石を用い、磁石の奥行を短くし且つ斜めに傾けることによって、開放感を高めたＭＲＩ装置１２０である。本実施形態では、これらの外観を有するＭＲＩ装置のいずれを用いることもできる。なお、これらは一例であり、本実施形態のＭＲＩ装置はこれらの形態に限定されるものではない。本実施形態では、装置の形態やタイプを問わず、公知の各種のＭＲＩ装置を用いることができる。以下、特に区別する必要がない場合は、ＭＲＩ装置１００で代表する。

図２（ａ）は、本実施形態のＭＲＩ装置１００の機能構成図である。本図に示すように、本実施形態のＭＲＩ装置１００は、被検体１が置かれる空間に、静磁場を生成する静磁場コイル２と、ｘ方向、ｙ方向、ｚ方向にそれぞれ傾斜磁場を発生させる傾斜磁場コイル３と、静磁場分布を調整するシムコイル４と、被検体１の計測領域に対し高周波磁場を照射する計測用高周波コイル５（以下、単に送信コイルという）と、被検体１から発生する核磁気共鳴信号を受信する受信用高周波コイル６（以下、単に受信コイルという）とを備える。

静磁場コイル２は、図１に示した各ＭＲＩ装置１００、１１０、１２０の構造に応じて、種々の形態のものが採用される。傾斜磁場コイル３及びシムコイル４は、それぞれ傾斜磁場用電源部１２及びシム用電源部１３により駆動される。なお、本実施形態では、送信コイル５と受信コイル６とに別個のものを用いる場合を例にあげて説明するが、送信コイル５と受信コイル６との機能を兼用する１のコイルで構成してもよい。送信コイル５が照射する高周波磁場は、送信機７により生成される。受信コイル６が検出した核磁気共鳴信号は、受信機８を通して計算機９に送られる。

シーケンス制御装置１４は、傾斜磁場発生コイル３の駆動用電源である傾斜磁場用電源部１２、シムコイル４の駆動用電源であるシム用電源部１３、送信機７及び受信機８の動作を制御し、傾斜磁場、高周波磁場の印加および核磁気共鳴信号の受信のタイミングを制御する。制御のタイムチャートはパルスシーケンスと呼ばれ、計測に応じて予め設定され、後述する計算機９が備える記憶装置等に格納される。

計算機９は、受け取った核磁気共鳴信号に対して様々な演算処理を行い、画像情報やスペクトル情報を生成するとともに、ＭＲＩ装置１００全体の動作を制御する。計算機９は、ＣＰＵ、メモリ、記憶装置などを備える情報処理装置であり、計算機９にはディスプレイ１０、外部記憶装置１１、入力装置１５などが接続される。ディスプレイ１０は、演算処理で得られた結果等をオペレータに表示するインタフェースである。入力装置１５は、本実施形態で行われる演算処理に必要な条件、パラメータ等をオペレータが入力するためのインタフェースである。外部記憶装置１１は、計算機９が備える記憶装置とともに、計算機９が実行する各種の演算処理に用いられるデータ、演算処理により得られるデータ、入力された条件、パラメータ等を保持する。

以下、本実施形態の計算機９が実現する機能について説明する。図２（ｂ）は、本実施形態の計算機９の機能ブロック図である。本実施形態の計算機９は、計測制御部９１０と、表示情報生成部９２０と、渦電流補正部９３０とを備える。計測制御部９１０は、パルスシーケンスに従ってシーケンス制御装置１４を動作させるとともに各部を制御し、計測を行う。表示情報生成部９２０は、計測により得られた核磁気共鳴信号に対して様々な演算処理を行い画像情報やスペクトル情報を生成する。

渦電流補正部９３０は、渦電流によるスペクトル歪みを補正する渦電流補正処理を行う。本実施形態では、渦電流によるスペクトル歪みを、代謝物質よりも信号強度の大きい水信号の位相データにより補正する。これを実現するため、渦電流補正部９３０は、得られた核磁気共鳴信号から水信号の位相データを算出する位相データ算出部９４０と、位相データ内の全ての位相とびを繰り返し計算によって補正する位相データ補正部９５０とを備える。さらに、本実施形態の位相データ補正部９５０は、位相データを時間方向に複数の小領域に分割する領域分割部９５３と、分割された小領域の中から位相とびが発生している領域を抽出する位相とび発生領域抽出部９５１と、前記位相とび発生領域抽出部９５１が抽出した小領域内の位相とびを補正する位相とび補正部９５２とを備える。

なお、計算機９が実現する各種の機能は、記憶装置が保持するプログラムを、ＣＰＵがメモリにロードして実行することにより実現される。また、計算機９が実現する各種の機能のうち、少なくとも１の機能は、ＭＲＩ装置１００とは独立した、情報処理装置であって、ＭＲＩ装置１００とデータの送受信が可能な情報処理装置により実現されていてもよい。

以下、上記各機能による、本実施形態の計測全体の流れについて簡単に説明する。本実施形態では、渦電流補正用信号として水信号を用いる。従って、まず、計測制御部９１０は、予め定められたパルスシーケンスに従ってシーケンス制御装置１４を制御し、非水抑圧計測を行い、水信号Ｆ（ｔ）を得る。なお、ここで、ｔは、サンプリング時間であり、離散値をとる変数である。得られた水信号Ｆ（ｔ）から位相データ算出部９４０は、位相データφ（ｔ）を算出する。そして、位相データ補正部９５０は、位相データφ（ｔ）に存在する全ての位相とびが補正された位相データφｃ（ｔ）を得る。次に、計測制御部９１０は、予め定められたパルスシーケンスに従ってシーケンス制御装置１４を制御し、水抑圧計測を行い、代謝物質信号Ｓ（ｔ）を得る。渦電流補正部９３０は、得られた代謝物質信号Ｓ（ｔ）を位相データφｃ（ｔ）で補正し、渦電流補正後の代謝物質信号Ｓ_ｅｃｃ（ｔ）を得る。表示情報生成部９２０は、代謝物質信号Ｓ_ｅｃｃ（ｔ）をフーリエ変換し、代謝物質のスペクトルまたは分布画像を得る。

ここで、計測制御部９１０が上記計測に用いるパルスシーケンスの一例について説明する。ここでは、代謝物質を画像化する領域選択型ＭＲＳＩのパルスシーケンス（以降、ＭＲＳＩパルスシーケンスと呼ぶ）を例にあげて説明する。

図３は、ＭＲＳＩパルスシーケンス３００の一例である。図３において、ＲＦは高周波磁場パルスの印加タイミングを示す。Ｇｘ、Ｇｙ、Ｇｚは、それぞれ、ｘ、ｙ、ｚ方向の傾斜磁場パルスの印加タイミングを示す。Ａ／Ｄは、信号の計測期間Ｔｐ１を示す。図３に示すＭＲＳＩパルスシーケンス３００は、公知のＭＲＳＩパルスシーケンスと同じである。すなわち、１つの励起パルスＲＦ１と２つの反転パルスＲＦ２およびＲＦ３とを用いて、所定の関心領域を選択的に励起し、この関心領域からＦＩＤ信号（自由誘導減衰）ＦＩＤ１を得ることを、繰り返し時間ＴＲ間隔で繰り返す。

このＭＲＳＩパルスシーケンス３００に従って、励起される領域を図４に示す。図４は、本計測に先立って行われる計測により得る位置決め用スカウト画像であって、それぞれ、図４（ａ）はトランス像、図４（ｂ）はサジタル像、図４（ｃ）はコロナル像である。以下、各部の動作と励起される領域との関係を図３および図４を用いて説明する。

まず高周波磁場ＲＦ１とz方向の傾斜磁場パルスＧｓ１、Ｇｓ１’とを印加して、ｚ方向の断面４０１を励起する。ＴＥ／４（ここで、ＴＥはエコー時間）時間後に、高周波磁場ＲＦ２とｙ方向の傾斜磁場パルスＧｓ２とを印加する。その結果、ｚ方向の断面４０１とｙ方向の断面４０２とが交差する領域における核磁化の位相のみが反転する。続いて、高周波磁場ＲＦ２印加からＴＥ／２後に高周波磁場ＲＦ３とｘ方向の傾斜磁場パルスＧｓ３とを印加する。それによって、ｚ方向の断面４０１、ｙ方向の断面４０２、ｘ方向の断面４０３が交差する関心領域４０４における核磁化の位相のみが反転し、ここから自由誘導減衰信号ＦＩＤ１が生じる。この自由誘導減衰信号ＦＩＤ１を計測する。なお各方向の傾斜磁場パルスＧｄ１〜Ｇｄ３およびＧｄ１’〜Ｇｄ３’は、高周波磁場ＲＦ１で励起された核磁化の位相は乱さず、ＲＦ２、ＲＦ３で励起された核磁化の位相をディフェイズするための傾斜磁場である。また、高周波磁場ＲＦ３の後には、位相エンコード傾斜磁場Ｇｐ１、Ｇｐ２を印加する。以上により、関心領域４０４の核磁気共鳴信号を得る。

次に、渦電流補正部９３０による渦電流補正を説明する。渦電流補正部９３０は、代謝物質信号をＳ（ｔ）から、位相とび補正後の水信号の位相データφｃ（ｔ）を用いて、渦電流補正後の代謝物質信号Ｓ_ｅｃｃ（ｔ）を、以下の式（１）に従って算出する。
Ｓ_ｅｃｃ（ｔ）＝Ｓ（ｔ）・ｅｘｐ（−ｉ・φｃ（ｔ）） （１）
ここで、ｉは虚数単位である。
なお、今後、ｅｘｐ（−ｉ・φｃ（ｔ））を位相補正量と呼ぶ。

まず、位相データ算出部９４０による位相データの算出処理を説明する。位相データ算出部９４０は、計測された水信号Ｆ（ｔ）から水信号の位相データφ（ｔ）を、以下の式（２）に従って算出する。
φ（ｔ）＝ｔａｎ^−１（Ｉｍ（Ｆ（ｔ）））／（Ｒｅ（Ｆ（ｔ））） （２）
ここで、ｔａｎ^−１は、アークタンジェント関数、Ｉｍ（Ｆ（ｔ））は複素数Ｆ（ｔ）の虚部、Ｒｅ（Ｆ（ｔ））は複素数Ｆ（ｔ）の実部を表す。

次に、位相データ補正部９５０による、位相データφ（ｔ）に存在する全ての位相とびを繰り返し計算により補正する位相データ補正処理について説明する。なお、位相データ補正部９５０は、位相データ補正処理の実施に先立ち、上記式（２）で得られた位相データφ（ｔ）に対し、位相折り返し接続処理を行う。まず、この位相折返し接続処理について図５を用いて説明する。図５において、横軸は、水信号の計測開始からの時間（ｍｓ）、縦軸は位相値（ｒａｄ）である。また、破線は計測結果から算出された位相データφ（ｔ）、実線は位相折返し接続処理後の位相データφｚ（ｔ）を示す。

位相データφ（ｔ）は、−πから＋πの間の値として算出される。しかし、上記ＭＲＳＩパルスシーケンス３００により取得する水信号の位相には、−πから＋πの範囲を超えた値が生じる。位相にこのような値が発生するタイミングでは、折り返されて−πから＋πの間の値が得られる。本図に示すように、折り返された部分では、値に不連続な状態が生じる。そこで、このような不連続な位相値の時間的な変化を取り除く位相折り返し接続処理を行い、本来の位相変化の状態である位相データφｚ（ｔ）を得る。なお、位相折り返し接続処理は、既存の各種の位相折り返し接続処理を用いることができる。このとき、ノイズによる位相ばらつきを防ぐため、適宜スムージングを行ってもよい。

次に、位相データ補正部９５０による、位相折返し接続処理後の位相データφｚ（ｔ）に対する位相データ補正処理について説明する。一般に位相とびが発生している箇所は、その近傍で位相値に急激な変化が発生し、そこに変曲点が生じる。すなわち、位相とびが生じる領域は、位相データの１次時間微分値が上あるいは下に凸となり、２次時間微分値が０となる点が含まれる。以降、位相データの１次時間微分値の変化量を単に「位相変化量」、２次時間微分値が０となる点を単に「変曲点」と呼ぶ。ここで、位相値に不規則な変化が発生している箇所であっても、その位相変化量が所定以下の箇所は、位相とびによるものではなく、渦電流の影響など、他の理由による。本実施形態では、位相データ補正部９５０は、位相変化量の大きさにより、位相データφｚ（ｔ）内の位相とびが発生している領域を特定し、位相とびが発生していると特定された領域内でのみ、位相とび補正を行う。

具体的には、位相データ補正部９５０は、位相データφｚ（ｔ）を時間方向に複数の小領域に分割し、小領域内で位相変化量を算出する。そして、この変化量の大きい領域から順に位相とびを補正する。このとき、補正する毎に、位相補正量をフーリエ変換して得る値に与える影響を評価する。この値に与える影響が略なくなったところで、位相とびが発生している箇所の補正を終えたものと判断する。

なお、本実施形態では、分割した各小領域内で、当該領域内の位相データφｚ（ｔ）の１次微分φｚ’（ｔ）の最大値と最小値との差分の絶対値を位相変化量とする。

図６は、位相データ補正部９５０による本実施形態の位相データφｚ（ｔ）の位相とび補正処理の処理フローである。まず、領域分割部９５３は、位相データφｚ（ｔ）を時間方向に複数の小領域に分割する領域分割処理を行う（ステップＳ６０１）。

次に、位相とび発生領域抽出部９５１は、ステップＳ６０１で分割した小領域毎に、位相変化量を算出する（ステップＳ６０２）。そして、位相とび発生領域抽出部９５１は、各小領域と算出した位相変化量とを対応づけて、記憶装置等に保持する（ステップＳ６０３）。

次に、位相とび発生領域抽出部９５１は、記憶装置に保持された位相変化量の中から、最も大きい位相変化量に対応づけて保持された小領域を、位相とび発生領域として抽出する（ステップＳ６０４）。

次に、位相とび補正部９５２は、ステップＳ６０４で抽出した位相とび発生領域において、位相とび補正領域を決定する（ステップＳ６０５）。その後、位相とび補正部９５２は、位相とび補正領域において、位相とびを補正する（ステップＳ６０６）。

次に、位相とび発生領域抽出部９５１は、予め定められた評価関数に従って評価値を算出する（ステップＳ６０７）。位相とび発生領域抽出部９５１は、算出した評価値が予め定められた終了条件を満たすか否かを判別する（ステップＳ６０８）。終了条件を満たさない場合、位相とび発生領域がまだ残っているものとして、ステップＳ６０６で補正を行った小領域に対応して記憶装置に保持された位相変化量を０とし（ステップＳ６０９）、ステップＳ６０４に戻る。一方、終了条件を満たす場合、位相データ補正部９５０は、処理を終了し、最終的に位相とび補正後の位相データφｃ（ｔ）を得る。なお、ステップＳ６０９の処理は、上記に限られず、処理済みの小領域を、以降の処理から除くことができるようステップＳ６０３で記憶した情報を更新するものであればよい。

ここで、上記ステップＳ６０１の領域分割処理の詳細について説明する。上述したように、位相とびが生じる領域は変曲点が含まれる。したがって、本実施形態では、領域分割部９５３は、位相データφｚ（ｔ）の変曲点を算出し、分割領域に変曲点が１つ含まれるように分割し、各小領域を得る。図７は、領域分割部９５３による、本実施形態の領域分割処理を説明するための図であり、（ａ）は、処理フロー、（ｂ）および（ｃ）は、分割手法の説明図である。

まず、領域分割部９５３は、水信号の位相データφｚ（ｔ）の２次微分φｚ”（ｔ）を計算し、その値が０となる点を、図７（ｂ）に示すように、水信号の位相データφｚ（ｔ）の変曲点８１０として抽出する（ステップＳ８０１）。なお、図７（ｂ）のグラフにおいて、横軸は、水信号の計測開始からの時間（ｍｓ）、縦軸は２次微分値（ａ．ｕ）である。ここでは、計測期間（ｔ_０〜ｔ_Ｆ）中に５個の変曲点８１０を抽出した場合を例示する。

次に、領域分割部９５３は、各変曲点８１０の時間ｔ_ＩＰｎ（ｎは１、２、３、４、５のいずれか）を特定し（ステップＳ８０２）、時間方向に順に境界点８１１を決定する（ステップＳ８０３）。ここでは、領域分割部９５３は、各変曲点８１０を時間方向に順に取り上げ、最初の変曲点８１０および最後の変曲点８１０については、それぞれ最寄りのデータ端点（ｔ_０またはｔ_Ｆ）との間の中間の時間ｔ_ＩＰＭ１、ｔ_ＩＰＭ６を、他の変曲点８１０の場合は、変曲点８１０（ｔ_ＩＰｎとｔ_{ＩＰｎ−１}）間の中間の時間ｔ_ＩＰＭｎを、それぞれ境界点８１１の時間とする。そして、領域分割部９５３は、境界点８１１により、水信号の位相データφｚ（ｔ）の時間領域を、図７（ｃ）に示すように、例えば、小領域８２１、８２２、８２３、８２４、８２５、８２６、８２７に分割する。ここで、図７（ｃ）のグラフにおいて、横軸は、水信号の計測開始からの時間（ｍｓ）、縦軸は位相データの位相値（ｒａｄ）である。

本実施形態では、領域分割部９５３は、以上の手順で分割した各領域を特定する情報、例えば、各小領域８２１、８２２、８２３、８２４、８２５、８２６、８２７に対応づけて、それぞれの両端の境界点８１１の時間（ｔ_０、ｔ_ＩＰＭ１）、（ｔ_ＩＰＭ１、ｔ_ＩＰＭ２）、（ｔ_ＩＰＭ２、ｔ_ＩＰＭ３）、（ｔ_ＩＰＭ３、ｔ_ＩＰＭ４）、（ｔ_ＩＰＭ４、ｔ_ＩＰＭ５）、（ｔ_ＩＰＭ５、ｔ_ＩＰＭ６）、（ｔ_ＩＰＭ６、ｔ_Ｆ）を記憶する。

次に、上記ステップＳ６０２における、位相とび発生領域抽出部９５１による位相変化量を算出する位相変化量算出処理について、図８を用いて説明する。図８は、位相データφｚ（ｔ）の、１次微分φｚ’（ｔ）のグラフ８３１であり、横軸は、水信号の計測開始からの時間ｔ（ｍｓ）、縦軸は１次微分値（ａ．ｕ）である。位相変化量算出処理では、位相とび発生領域抽出部９５１は、まず、位相データφｚ（ｔ）の１次微分φｚ’（ｔ）を求め、ステップＳ６０１の分割で得た小領域毎に、当該領域内の１次微分φｚ’（ｔ）の最大値と最小値との差の絶対値を算出する。そして、位相とび発生領域抽出部９５１は、算出した絶対値を、位相変化量とする。

なお、位相変化量を算出する位相変化量算出処理はこれに限られない。位相データφｚ（ｔ）の１次微分φｚ’（ｔ）を求め、ステップＳ６０１の分割で得た小領域毎に、ガウシアンなどのモデル関数で１次微分φｚ’（ｔ）をフィッティングし、各領域のフィッティング関数の高さを位相変化量としてもよい。

次に、位相とび補正部９５２による、ステップＳ６０５の位相とび補正領域を決定する補正領域決定処理と、ステップＳ６０６の位相とび補正処理と、について、図９を用いて説明する。図９は、位相とび発生領域として決定された小領域近傍の１次微分グラフ８３１であり、横軸は時間（ｍｓ）、縦軸は１次微分値（ａ．ｕ．）である。ここでは、一例として小領域８２２が、位相とび発生領域８４１と決定されたものとして説明する。

補正領域決定処理では、位相とび補正部９５２は、位相とび発生領域８４１において、一次微分φｚ’（ｔ）のピーク半値幅を算出し、その２倍の幅を位相とび補正領域８４２とする。なお、位相とび補正領域１００２の決定手法はこれに限られない。当該小領域における位相変化量が、最小となるような領域であればよく、例えば、変曲点に対応する１次微分ピーク点の両側を１点ずつ除外したときに、位相変化量が最小となるように決定してもよい。また、上記のように位相変化量算出処理において、ガウシアンなどのモデル関数でフィッティングし、位相とび発生領域８４１を決定した場合は、フィッティング関数の半値幅の２倍を位相とび補正領域８４２とするよう構成してもよい。

また、位相とび補正処理では、位相とび補正部９５２は、例えば、図９に示すように、１次微分φｚ’（ｔ）のグラフ上で、位相とび補正領域８４２の両端を直線で結ぶ線形補間により補正後の１次微分φｃ’（ｔ）を決定する。そして、位相とび補正部９５２は、補正後の１次微分φｃ’（ｔ）を積分することで位相とび補正後の位相値φｃ（ｔ）を算出する。ここで、補正後の１次微分φｃ’（ｔ）を、積分ではなく累乗和から位相とび補正後の位相値φｃ（ｔ）を算出してもよい。

ここで、上記位相とび補正処理の前後の位相データのグラフの一例を図１０（ａ）に示す。破線が位相とび補正前の水信号の位相データφｚ（ｔ）のグラフ、実線が、位相とび発生領域８４１（小領域８２２）の位相とびを補正後の水信号位相データφｃ（ｔ）のグラフである。本図において、横軸は、時間（ｍｓ）、縦軸は、位相（ｒａｄ）を示す。本図に示すように、上記位相とび補正処理により、水信号の位相とびが解消されていることがわかる。

なお、位相とび補正の手法は、これに限られない。例えば、位相とび補正部９５２は、位相とび補正領域ｒの両側で一次微分φ’（ｔ）が滑らかな曲線になるよう、位相とび補正領域ｒを補間し、補正後の１次微分φｃ’（ｔ）を決定してもよい。

次に、ステップＳ６０７およびＳ６０８における、位相とび発生領域抽出部９５１による評価処理について説明する。本実施形態では、位相とび発生領域抽出部９５１は、ｎ回目の位相とび補正後の位相データφ_ｎ（ｔ）に関して、以下の式（３）を計算する。
ｇ_ｎ（ｔ）＝ｅｘｐ（ｉ・φ_ｎ（ｔ）） （３）
次に、位相とび発生領域抽出部９５１は、ｇ_ｎ（ｔ）をフーリエ変換したＧ_ｎ（ｆ）を算出する（式（４））。
Ｇ_ｎ（ｆ）＝Ｆ[ｇｎ（ｔ）] （４）
ここで、Ｆ[]は、フーリエ変換を示す。位相とび発生領域抽出部９５１は、算出したＧ_ｎ（ｆ）の絶対値｜Ｇ_ｎ（ｆ）｜の最大ピークＰ_ｎを算出する。そして、位相とび発生領域抽出部９５１は、最大ピークＰ_ｎと、ｎ−１回目の位相とび補正後の最大ピークＰ_ｎ−１との差分の絶対値｜Ｐ_ｎ−Ｐ_ｎ−１｜を評価関数とし、この評価関数を計算することにより得られるΔＰ_ｎを評価値とする。

また、本実施形態の評価処理の終了条件は、得られた評価値が閾値未満となることとする。閾値としては、限りなく０に近い値を設定する。位相とび補正後の｜Ｇ_ｎ（ｆ）｜の最大ピークＰ_ｎの変化が０に漸近するということは、これ以上、新たな領域を抽出して位相とび補正を行っても効果が変わらないことを意味する。なぜなら、渦電流やノイズなどに起因する位相とび以外の位相変化は、上記評価値を著しく変化させるほど大きな変化ではないからである。従って、本実施形態では、本終了条件を満たす場合、位相とびが発生している全ての箇所の位相とび補正を終えたものと判断し、処理を終了させる。

なお、評価関数として、｜Ｇ_ｎ（ｆ）｜の最大ピーク値の位相とび補正毎の変化ではなく、｜Ｇ_ｎ（ｆ）｜の標準偏差の位相とび補正毎の変化を用いてもよい。また、｜Ｇ_ｎ（ｆ）｜が閾値Ｔ以上の値をとる信号値ｆの数Ｎｎの、位相とび補正ごとの変化を用いてもよい。

図１０（ｂ）に、位相とび補正処理を実施する前の位相データのグラフと、全ての位相とび領域を補正した後の位相データのグラフの一例を示す。破線が本実施形態の位相とび補正処理前の水信号の位相データφｚ（ｔ）のグラフである。そして、実線が、全ての位相とびを補正した後の水信号の位相データφｃ（ｔ）のグラフである。本図において、横軸は、時間（ｍｓ）、縦軸は、位相（ｒａｄ）を示す。本図に示すように、上記繰り返し計算による位相とび補正処理により、水信号の全ての位相とびが解消されていることがわかる。

以上説明したように、本実施形態によれば、渦電流補正用信号の位相データを時間方向に複数の領域に分割し、その中で位相変化量が最も大きい領域を抽出する。その領域の位相とびを補正し、評価関数を基に評価値を算出する。この処理を繰り返し、評価値が０に近づいたとき、処理を終了する。このように構成することで、複数の位相変化がある位相とびを自動で判別し除去する。また、除去後の位相データを用い、渦電流補正を実施する。

本実施形態によれば、評価値が０に漸近するまで、位相とび発生領域を抽出し、位相とび補正を繰り返すため、位相データの値に不規則な変化が発生している箇所において、位相とび箇所のみ抽出して補正することができる。従って、効率よく、また、精度よく、位相データの位相とびを除去することができる。これにより、渦電流補正において、位相とびのない位相データを用いることができ、リンギングアーチファクトが生じない良好なスペクトルを取得することができる。

ここで、本実施形態による渦電流補正の効果を計算機シミュレーションおよびファントム実験にて確認する。

計算機シミュレーション結果を図１１に示す。ここで、図１１（ａ）は、渦電流補正を行わなかった場合の代謝物質スペクトル、（ｂ）は、本実施形態の位相とび補正を行わない位相データを用いて渦電流補正を行った場合の代謝物質スペクトル、（ｃ）は、本実施形態の位相とび補正を実施した位相データを用いて渦電流補正を行った場合の代謝物質スペクトルをそれぞれ示す。

これらの図に示すように、本実施形態の位相とび補正後に渦電流補正を行うと、位相とびによるリンギングアーチファクトが除去でき、かつ、渦電流によるスペクトル歪みも補正できることがわかる。

次に、ファントムを用い、実際に計測した結果を図１２に示す。ここで用いたファントムは、濃度１０ｍＭのアセチルアラニンファントムである。図１２（ａ）は、渦電流補正を行わなかった場合の代謝物質スペクトル、（ｂ）は、本実施形態の位相とび補正を行わない位相データを用いて渦電流補正を行った場合の代謝物質スペクトル、（ｃ）は、本実施形態の位相とび補正を実施した位相データを用いて渦電流補正を行った場合の代謝物質スペクトルをそれぞれ示す。

これらの図に示すように、ファントムを用いて実際に計測した結果においても、本実施形態の位相とび補正後に渦電流補正を行うと、位相とびによるリンギングアーチファクトが除去でき、かつ、渦電流によるスペクトル歪みも補正できることがわかる。

以上、説明したように、本実施形態によれば、渦電流歪みを補正するために用いる位相データに、位相変化量が大きく異なる等の不規則な変化をする箇所が複数ある場合であっても、リンギングアーチファクトの原因となる位相とびの箇所のみを抽出して位相とび補正を行うことができるため、位相データの位相とびを効率よく補正することができる。従って、渦電流によるスペクトル歪みの補正にあたり、リングングアーチファクトの発生を効果的に防ぐことができ、渦電流によるスペクトル歪みを良好に補正することができる。

このため、本実施形態によれば、渦電流の影響を補正する処理を行うとき、補正処理によるアーチファクトを防ぎ、補正結果の精度を高めることができる。例えば、ＭＲＳＩにおいて、リンギングアーチファクトのない、良好な代謝物質スペクトルを得ることができる。

なお、本実施形態では、渦電流の補正に用いる位相データを水信号から得る場合を例にあげているが、これに限られない。計測対象の代謝物質よりも信号強度の大きい物質の信号であれば、特に水信号に限られない。

また、上記位相とび補正処理では、領域分割部９５３が領域分割後、位相とび発生領域抽出部９５１が全小領域について、位相変化量を算出し、記憶装置に保持しているが、これに限られない。例えば、位相とび補正処理を、位相とび補正毎に、領域を分割し直し、位相変化量が最も大きい小領域を算出するよう構成してもよい。すなわち、ステップＳ６０２の後に、その結果を保持せず、ステップＳ６０４の位相とび発生領域を抽出し、ステップＳ６０５、Ｓ６０６、Ｓ６０７、Ｓ６０８と処理を進め、条件を満たさない場合、ステップＳ６０６で補正後の位相データについて、ステップＳ６０１から処理を繰り返す。

さらに、上記実施形態では、評価関数を用い、位相とび補正が渦電流補正係数に与える影響を評価し、位相とび発生領域を繰り返し抽出しているが、これに限られない。例えば、予め閾値を定めておき、位相変化量が閾値以上である領域を位相とび発生領域として抽出し、抽出した全領域に対し、位相とび補正処理を行うよう構成してもよい。

また、一般に、局所的な時間範囲で考えた場合、位相とびは、強度および周波数の異なる２つのベクトルのベクトル和によって生じる位相変化である。そのため、２つのベクトルの強度が近ければ近いほど、２つのベクトルが対向する短い時間幅で急激に位相が変化する。逆に、２つのベクトルの強度が異なるほど、強度が小さいベクトルの影響が小さくなるため、対向する時間幅で緩やかに位相が変化する。したがって、位相変化量の大小と、位相とび補正領域の時間幅（変化幅）とは、反比例の関係にある。これを利用し、位相変化量の大きい領域から順に位相とびを補正するのではなく、前記変化幅の小さい領域から順に、位相とび補正の影響がなくなるまで、位相とび補正を行うよう構成してもよい。この場合、位相折返し接続処理後、ノイズによる影響を軽減するためスムージング処理を行う。

また、位相とび補正を行う領域を決定するに当たり、位相変化量の大きさと、上述の変化幅とを併用するよう構成してもよい。例えば、予めそれぞれ閾値を定めておき、位相変化量が予め定めた第一の閾値以上であり、かつ、変化幅が、予め定めた第二の閾値以下である場合、位相とび発生領域と判別し、位相とび補正を行うよう構成してもよい。

１：被検体、２：静磁場コイル、３：傾斜磁場コイル、４：シムコイル、５：送信コイル、６：受信コイル、７：送信機、８：受信機、９：計算機、１０：ディスプレイ、１１：外部記憶装置、１２：傾斜磁場用電源部、１３：シム用電源部、１４：シーケンス制御装置、１５：入力装置、１００：ＭＲＩ装置、１１０：ＭＲＩ装置、１２０：ＭＲＩ装置、３００：ＭＲＳＩパルスシーケンス、４０１：断面、４０２：断面、４０３：断面、４０４：関心領域、８１０：変曲点、８１１：境界点、８２１：小領域、８２２：小領域、８２３：小領域、８２４：小領域、８２５：小領域、８２６：小領域、８２７：小領域、８３１：グラフ、８４１：位相とび発生領域、８４２：位相とび補正領域、９１０：計測制御部、９２０：表示情報生成部、９３０：渦電流補正部、９４０：位相データ算出部、９５０：位相データ補正部、９５１：位相とび発生領域抽出部、９５２：位相とび補正部、９５３：領域分割部

Claims (18)

1. 被検体に静磁場を印加する静磁場印加手段と、前記被検体に傾斜磁場を印加する傾斜磁場印加手段と、前記被検体に高周波磁場パルスを照射する高周波磁場パルス照射手段と、前記被検体から核磁気共鳴信号を受信する受信手段と、前記核磁気共鳴信号に渦電流補正を実施する渦電流補正手段と、を備える磁気共鳴撮影装置であって、
   前記渦電流補正手段は、
   計測対象の代謝物質より信号強度の大きい物質の核磁気共鳴信号の位相データを算出する位相データ算出手段と、
   前記位相データを、時間方向に複数の小領域に分割する位相データ分割手段と、
   前記小領域の中から位相とびが存在する位相とび発生領域を抽出する位相とび発生領域抽出手段と、
   前記位相とび発生領域内の位相とびを補正する位相とび補正手段と、を備えること
   を特徴とする磁気共鳴撮影装置。
2. 請求項１記載の磁気共鳴撮影装置であって、
   前記位相とび発生領域抽出手段は、前記位相データの位相変化量が、予め定められた条件を満たす領域を、前記位相とび発生領域として抽出すること
   を特徴とする磁気共鳴撮影装置。
3. 請求項１または２記載の磁気共鳴撮影装置であって、
   前記位相とび補正手段は、前記位相とび発生領域抽出手段により位相とび発生領域が抽出される毎に前記位相とびの補正を行い、
   前記位相とび発生領域抽出手段は、
   前記位相とび補正手段によって位相とびが補正される毎に、当該補正後の前記位相データの評価値を算出する評価値算出手段と、
   前記評価値が算出される毎に、前記評価値が所定の閾値以上であるか否かを判別し、閾値以上である場合、前記位相とび発生領域の抽出を継続するよう決定し、閾値未満である場合、前記抽出を終了するよう決定する判別手段と、を備え、
   前記判別手段が抽出を継続すると判別する間、前記位相データの位相変化量が最も大きい小領域から順に１の前記小領域を前記位相とび発生領域として抽出すること
   を特徴とする磁気共鳴撮影装置。
4. 請求項１から３いずれか１項記載の磁気共鳴撮影装置であって、
   前記位相データ分割手段は、
   前記位相データの変曲点を算出する変曲点算出手段を備え、
   前記各小領域に１の前記変曲点算出手段が算出した変曲点が含まれるよう分割すること
   を特徴とする磁気共鳴撮影装置。
5. 請求項３記載の磁気共鳴撮影装置であって、
   前記位相とび発生領域抽出手段は、
   前記位相データの１次微分値を算出する１次微分算出手段と、
   前記小領域ごとに、前記１次微分値の最大値と最小値との差分の絶対値を算出する位相変化量算出手段と、を備え、
   前記絶対値が最大の小領域を、前記位相変化量が最も大きい小領域とすること
   を特徴とする磁気共鳴撮影装置。
6. 請求項３記載の磁気共鳴撮影装置であって、
   前記位相とび発生領域抽出手段は、
   前記位相データの１次微分値を算出する変化量算出手段と、
   前記小領域毎に、前記一次微分値をモデル関数でフィッティングするフィッティング手段と、
   前記小領域毎に、前記フィッティング後のモデル関数の高さを算出する高さ算出手段と、を備え、
   前記高さが最も高い小領域を、前記位相変化量が最も大きい小領域とすること
   を特徴とする磁気共鳴撮影装置。
7. 請求項５記載の磁気共鳴撮影装置であって、
   前記位相とび補正手段は、
   前記位相とび発生領域として抽出された小領域における、前記１次微分値のピーク値の半値幅を、前記位相とびを補正する位相とび補正領域と設定する位相とび補正領域設定手段を備えること
   を特徴とする磁気共鳴撮影装置。
8. 請求項６記載の磁気共鳴撮影装置であって、
   前記位相とび補正手段は、
   前記フィッティング後のモデル関数の半値幅の２倍の領域を、前記位相とびを補正する位相とび補正領域と設定する位相とび補正領域設定手段を備えること
   を特徴とする磁気共鳴撮影装置。
9. 請求項７または８記載の磁気共鳴撮影装置であって、
   前記位相とび補正手段は、前記位相とび補正領域の前記１次微分値を補間により接続し、前記位相とびを補正すること
   を特徴とする磁気共鳴撮影装置。
10. 請求項９記載の磁気共鳴撮影装置であって、
    前記補間は、前記１次微分値の前記位相とび補正領域の両端を直線で接続する直線補間であること
    を特徴とする磁気共鳴撮影装置。
11. 請求項９または１０記載の磁気共鳴撮影装置であって、
    前記位相とび補正手段は、前記補間により接続後の前記１次微分値を積分し、前記位相とび補正後の位相データとすること
    を特徴とする磁気共鳴撮影装置。
12. 請求項９または１０記載の磁気共鳴撮影装置であって、
    前記位相とび補正手段は、前記補間により接続後の前記１次微分値の累乗和を、前記位相とび補正後の位相データとすること
    を特徴とする磁気共鳴撮影装置。
13. 請求項３から１２いずれか１項記載の磁気共鳴撮影装置であって、
    前記評価値算出手段は、前記補正後の位相データから算出した補正スペクトルのピークの高さと、補正前の位相データから算出した補正スペクトルのピークの高さとの差分を前記評価値とすること
    を特徴とする磁気共鳴撮影装置。
14. 請求項３から１２いずれか１項記載の磁気共鳴撮影装置であって、
    前記評価値算出手段は、前記補正後の位相データから算出した補正スペクトルの標準偏差と、補正前の位相データから算出した補正スペクトルの標準偏差との差分を前記評価値とすること、
    を特徴とする磁気共鳴撮影装置。
15. 請求項３から１２いずれか１項記載の磁気共鳴撮影装置であって、
    前記評価値算出手段は、前記補正後の位相データから算出した信号値であって、予め定めた値以上の値を有する信号値の数と、補正前の位相データから算出した信号値であって、予め定めた値以上の値を有する信号値の数との差を、前記評価値とすること
    を特徴とする磁気共鳴撮影装置。
16. 被検体に静磁場を印加する静磁場印加手段と、前記被検体に傾斜磁場を印加する傾斜磁場印加手段と、前記被検体に高周波磁場パルスを照射する高周波磁場パルス照射手段と、前記被検体から核磁気共鳴信号を受信する受信手段と、前記核磁気共鳴信号に渦電流補正を実施する渦電流補正演算手段と、を備える磁気共鳴撮影装置において、記受信手段で取得した核磁気共鳴信号の中で、計測対称の代謝物質より信号強度の大きい物質の核磁気共鳴信号である渦電流補正信号の位相データに発生する全ての位相とびを除去する位相データ補正方法であって、
    前記位相データを、時間方向に複数の小領域に分割する位相データ分割ステップと、
    前記小領域の中から位相とびが存在する位相とび発生領域を抽出する位相とび発生領域抽出ステップと、
    前記抽出した領域における位相とびを補正する位相とび補正ステップと、を備えること
    を特徴とする位相データ補正方法。
17. 請求項１６記載の位相データ補正方法であって、
    前記位相とび補正ステップによって補正後の位相データの評価値を算出する評価値算出ステップと、
    前記評価値が所定の閾値以上であれば、前記位相とび発生領域抽出ステップと、前記位相とび補正ステップと、前記評価値算出ステップとを実行する繰返しステップと、をさらに備えること
    を特徴とする位相データ補正方法。
18. 被検体に静磁場を印加する静磁場印加手段と、前記被検体に傾斜磁場を印加する傾斜磁場印加手段と、前記被検体に高周波磁場パルスを照射する高周波磁場パルス照射手段と、前記被検体から核磁気共鳴信号を受信する受信手段と、前記核磁気共鳴信号に渦電流補正を実施する渦電流補正手段と、を備える磁気共鳴撮影装置のコンピュータを、
    前記受信手段で取得した核磁気共鳴信号の中で、計測対称の代謝物質より信号強度の大きい物質の核磁気共鳴信号である渦電流補正信号の位相データを算出する位相データ算出手段と、
    前記位相データを、時間方向に複数の小領域に分割する位相データ分割手段と、
    前記小領域の中から位相とびが存在する位相とび発生領域を抽出する位相とび発生領域抽出手段と、
    前記位相とび発生領域内の位相とびを補正する位相とび補正手段と
    して機能させることを特徴とするプログラム。

Images