(57)【要約】      （修正有）【課題】強度プロファイルに現れる脂肪などの不要信号による強度を低減する磁気共鳴装置を提供する。【解決手段】ＭＲ装置１００は、ＭＲ信号を受信する受信コイル４を有する。また、ＭＲ装置１００は、受信コイル４により受信されたＭＲ信号に基づいて、ＭＲ信号の位相を表す位相情報を含む位相プロファイルと、ＭＲ信号の強度を表す強度情報を含む強度プロファイルとを生成する生成手段と、位相プロファイルに基づいて、強度プロファイルの中から範囲ｘを求める手段と、範囲ｘにおける強度をゼロに変更する強度変更手段とを有する。【選択図】図１
2016010582.tif 000002
【特許請求の範囲】
【請求項１】
被検体からＭＲ信号を発生させるためのシーケンスを実行する磁気共鳴装置であって、
前記ＭＲ信号を受信するコイルと、
前記コイルにより受信されたＭＲ信号に基づいて、前記ＭＲ信号の位相を表す位相情報を含む位相プロファイルと、前記ＭＲ信号の強度を表す強度情報を含む強度プロファイルとを生成する生成手段と、
前記位相プロファイルに基づいて、前記強度プロファイルの全範囲の中から一部の範囲を求める手段と、
前記一部の範囲における第１の強度を、前記第１の強度よりも小さい第２の強度に変更する強度変更手段と、を有する磁気共鳴装置。
【請求項２】
前記一部の範囲を求める手段は、
前記位相プロファイルに基づいて、前記強度プロファイルの中から、前記一部の範囲を求めるための基準位置を特定する位置特定手段と、
前記基準位置に基づいて、前記一部の範囲を特定する範囲特定手段と、
を有する請求項１に記載の磁気共鳴装置。
【請求項３】
前記シーケンスにより励起される第１の領域を設定する領域設定手段を有する、請求項２に記載の磁気共鳴装置。
【請求項４】
前記位相プロファイルは、前記第１の領域内の第１の方向の各位置における位相を表しており、
前記強度プロファイルは、前記第１の領域内の前記第１の方向の各位置における強度を表している、請求項３に記載の磁気共鳴装置。
【請求項５】
前記位置特定手段は、
前記位相プロファイルに基づいて、前記第１の方向において互いに隣接する位置における位相の位相差を表す位相差プロファイルを作成し、前記位相差プロファイルに基づいて、前記基準位置を特定する、請求項４に記載の磁気共鳴装置。
【請求項６】
前記位置特定手段は、
前記位相差を閾値と比較し、その比較結果に基づいて前記基準位置を特定する、請求項５に記載の磁気共鳴装置。
【請求項７】
前記位置特定手段は、
前記第１の領域の前記１の方向の中心位置から近い順に、前記位相差と前記閾値との比較を行い、その比較結果に基づいて前記基準位置を特定する、請求項６に記載の磁気共鳴装置。
【請求項８】
前記コイルは、
前記ＭＲ信号を受信する複数のチャネルを有しており、
前記作成手段は、
前記チャネルごとに、前記位相プロファイルおよび前記強度プロファイルを作成し、
前記強度変更手段は、
前記チャネルごとに、前記強度プロファイルの第１の強度を第２の強度に変更する、請求項３〜７のうちのいずれか一項に記載の磁気共鳴装置。
【請求項９】
強度が変更された後の各チャネルの強度プロファイルを合成することにより、合成プロファイルを作成する合成プロファイル作成手段を有する、請求項８に記載の磁気共鳴装置。
【請求項１０】
前記第１の領域は第１の部位および第２の部位を含んでおり、
前記範囲特定手段は、
前記強度プロファイルのうち、前記基準位置に対して前記第１の部位側の部分を、前記一部の範囲として特定する、請求項３〜９のうちのいずれか一項に記載の磁気共鳴装置。
【請求項１１】
前記第１の部位は肺を含んでおり、前記第２の部位は肝臓を含んでいる、請求項１０に記載の磁気共鳴装置。
【請求項１２】
前記合成プロファイルに基づいて、前記第１の部位のエッジの位置を検出する検出手段を有する、請求項９〜１１のうちのいずれか一項に記載の磁気共鳴装置。
【請求項１３】
前記磁気共鳴装置は、
前記第１の領域からＭＲ信号を収集するための第１のシーケンスと、前記第１の領域からＭＲ信号を収集するための第２のシーケンスとを実行し、
前記合成プロファイル作成手段は、
前記第１のシーケンスにより得られた第１の合成プロファイルと、前記第２のシーケンスにより得られた第２の合成プロファイルとを作成する、請求項９〜１１のうちのいずれか一項に記載の磁気共鳴装置。
【請求項１４】
前記第１の合成プロファイルと前記第２の合成プロファイルとに基づいて、前記第１の部位のエッジの移動量を算出する算出手段を有する、請求項１３に記載の磁気共鳴装置。
【請求項１５】
前記算出手段は、
前記第１の合成プロファイルと前記第２の合成プロファイルとの二乗誤差又は相関に基づいて前記移動量を算出する、請求項１４に記載の磁気共鳴装置。
【請求項１６】
前記第２の強度はゼロである、請求項１〜１５のうちのいずれか一項に記載の磁気共鳴装置。
【請求項１７】
前記一部の範囲には、第１の組織によるＭＲ信号の強度が含まれている、請求項１〜１６のうちのいずれか一項に記載の磁気共鳴装置。
【請求項１８】
前記第１の組織によるＭＲ信号は、前記第１の領域の内側又は外側で発生する、請求項１７に記載の磁気共鳴装置。
【請求項１９】
前記第１の組織は脂肪を含んでいる、請求項１７又は１８に記載の磁気共鳴装置。
【請求項２０】
被検体からＭＲ信号を発生させるためのシーケンスを実行する磁気共鳴装置であって、前記ＭＲ信号を受信するコイルを有する磁気共鳴装置に適用されるプログラムであって、
前記コイルにより受信されたＭＲ信号に基づいて、前記ＭＲ信号の位相を表す位相情報を含む位相プロファイルと、前記ＭＲ信号の強度を表す強度情報を含む強度プロファイルとを生成する生成処理と、
前記位相プロファイルに基づいて、前記強度プロファイルの全範囲の中から一部の範囲を求める処理と、
前記一部の範囲における第１の強度を、前記第１の強度よりも小さい第２の強度に変更する強度変更処理と、
をコンピュータに実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ＭＲ信号の位相情報および強度情報を取得する磁気共鳴装置、およびこの磁気共鳴装置に適用されるプログラムに関する。
【背景技術】
【０００２】
腹部の撮影では、動きアーチファクトを低減するために呼吸同期法を用いた撮影方法が実行される。呼吸同期法の撮影方法として、肺と肝臓との境界部分に肝臓のエッジの動きを検出するためのナビゲータ領域を設定し、ナビゲータエコーを収集する撮影方法が知られている（特許文献１参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１０−１５５０２１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
特許文献１では、ナビゲータ領域を励起するためのナビゲータシーケンスが実行され、ナビゲータ領域からナビゲータエコーを収集する。そして、ナビゲータエコーに基づいて、ナビゲータエコーの強度情報を表す強度プロファイルを作成し、肝臓のエッジの動きを検出している。
【０００５】
しかし、例えば、ペンシルビーム（pencil beam）励起型のナビゲータシーケンスを実行すると、脂肪などの不要信号がナビゲータエコーに混入してしまい、強度プロファイルに不要信号による強度が現れることがある。このような不要信号による強度は、肝臓のエッジの検出誤差の原因となる。
したがって、強度プロファイルに現れる脂肪などの不要信号による強度を低減する技術が望まれている。
【課題を解決するための手段】
【０００６】
本発明の第１の観点は、被検体からＭＲ信号を発生させるためのシーケンスを実行する磁気共鳴装置であって、
前記ＭＲ信号を受信するコイルと、
前記コイルにより受信されたＭＲ信号に基づいて、前記ＭＲ信号の位相を表す位相情報を含む位相プロファイルと、前記ＭＲ信号の強度を表す強度情報を含む強度プロファイルとを生成する生成手段と、
前記位相プロファイルに基づいて、前記強度プロファイルの全範囲の中から一部の範囲を求める手段と、
前記一部の範囲における第１の強度を、前記第１の強度よりも小さい第２の強度に変更する強度変更手段と、を有する磁気共鳴装置である。
【０００７】
本発明の第２の観点は、被検体からＭＲ信号を発生させるためのシーケンスを実行する磁気共鳴装置であって、前記ＭＲ信号を受信するコイルを有する磁気共鳴装置に適用されるプログラムであって、
前記コイルにより受信されたＭＲ信号に基づいて、前記ＭＲ信号の位相を表す位相情報を含む位相プロファイルと、前記ＭＲ信号の強度を表す強度情報を含む強度プロファイルとを生成する生成処理と、
前記位相プロファイルに基づいて、前記強度プロファイルの全範囲の中から一部の範囲を求める処理と、
前記一部の範囲における第１の強度を、前記第１の強度よりも小さい第２の強度に変更する強度変更処理と、をコンピュータに実行させるためのプログラムである。
【発明の効果】
【０００８】
前記位相プロファイルを用いることにより、強度プロファイルに現れる不要な信号の強度を低減することができる。
【図面の簡単な説明】
【０００９】
【図１】本発明の第１の形態の磁気共鳴装置の概略図である。
【図２】受信コイル４の概略拡大図である。
【図３】プロセッサ９が実行する処理を示す図である。
【図４】第１の形態で実行されるスキャンを示す図である。
【図５】図４に示すスキャンを実行するときのフローを示す図である。
【図６】ローカライザスキャンＬＸにより取得された画像ＬＤを概略的に示す図である。
【図７】撮影領域を概略的に示す図である。
【図８】プレスキャンＡで実行されるシーケンスの説明図である。
【図９】チャネルごとに得られたナビゲータデータＮＤ１〜ＮＤｎを概略的に示す図である。
【図１０】各チャネルのナビゲータデータＮＤ１〜ＮＤｎから得られた強度プロファイルＩＮ１〜ＩＮｎの一例を概略的に示す図である。
【図１１】強度プロファイルＩＮ１〜ＩＮｎを合成することにより得られた合成プロファイルＩＮを概略的に示す図である。
【図１２】不要信号の強度が現れた強度プロファイルの一例を概略的に示す図である。
【図１３】第１の形態におけるエッジの検出方法のフローを示す図である。
【図１４】強度プロファイルＩＮ１〜ＩＮｎの一例を示す図である。
【図１５】チャネルごとに作成された位相プロファイルＰＨ１〜ＰＨｎの一例を示す図である。
【図１６】位相プロファイルＰＨ１と、位相プロファイルＰＨ１から作成された位相差プロファイルＤＩＦ１とを示す図である。
【図１７】閾値ＴＨを用いて位相差Δθｉが大きい値を持つ位置Ｐｉを特定する方法の説明図である。
【図１８】強度プロファイルＩＮ１、位相差プロファイルＤＩＦ１、および位相プロファイルＰＨ１における位置Ｐｃ＋ｍを示す図である。
【図１９】範囲ｘを示す図である。
【図２０】範囲ｘの強度がゼロに変更された後の強度プロファイルＩＮ１´を示す図である。
【図２１】範囲ｘの強度がゼロに変更される前の強度プロファイルＩＮ１と、範囲ｘの強度がゼロに変更された後の強度プロファイルＩＮ１´とを示す図である。
【図２２】不要信号の強度がゼロに変更された後のチャネルＣＨ１〜ＣＨｎの強度プロファイルＩＮ１´〜ＩＮｎ´を示す図である。
【図２３】合成プロファイルＦ１を示す図である。
【図２４】検出された肝臓のエッジの位置Ｐａ１を示す図である。
【図２５】ナビゲータシーケンスＮＳ１〜ＮＳｆの各々を実行したときの肝臓のエッジの位置の時間変化を概略的に示す図である。
【図２６】ウインドウＷの一例を示す図である。
【図２７】本スキャンＢの説明図である。
【図２８】強度プロファイルＩＮ１の範囲ｘの強度を１／１０の値に変更したときの強度プロファイルＩＮ１´を概略的に示す図である。
【図２９】第２の形態におけるプロセッサ９が実行する処理を示す図である。
【図３０】ナビゲータシーケンスＮＳ１を実行し、合成プロファイルを作成するときのフローを示す図である。
【図３１】ナビゲータシーケンスＮＳ１を実行することにより得られた合成プロファイルＦ１を概略的に示す図である。
【図３２】ナビゲータシーケンスＮＳ２を実行するときのフローを示す図である。
【図３３】ナビゲータシーケンスＮＳ２を実行することにより得られた合成プロファイルＦ２を概略的に示す図である。
【図３４】合成プロファイルＦ１およびＦ２に基づいて肝臓のエッジの移動量を算出するときの説明図である。
【図３５】算出された二乗誤差ＳＥ＝ＳＥ０を概略的に示す図である。
【図３６】合成プロファイルＦ２をＩ方向に移動させ、二乗誤差を算出するときの説明図である。
【図３７】合成プロファイルＦ２をＳ方向に移動させ、二乗誤差を算出するときの説明図である。
【図３８】合成プロファイＦ１およびＦ２を用いて算出された肝臓のエッジの移動量ｄ＝Δｄ２を示す図である。
【図３９】ナビゲータシーケンスＮＳ３を実行することにより得られた合成プロファイルＦ３を概略的に示す図である。
【図４０】算出された移動量ｄ＝Δｄ３を示す図である。
【図４１】ナビゲータシーケンスＮＳ１〜ＮＳｆを実行することにより求められた肝臓のエッジの移動量ｄを概略的に示す図である。
【図４２】ウインドウＷの一例を示す図である。
【図４３】本スキャンＢで実行されるナビゲータシーケンスＮＳｆ＋１〜ＮＳｕおよびイメージングシーケンスＤＡＱ１〜ＤＡＱｚを概略的に示す図である。
【図４４】ナビゲータシーケンスＮＳｆ＋１を実行することにより得られた合成プロファイルＦｆ＋１を概略的に示す図である。
【図４５】肝臓のエッジの移動量を算出するときの方法を示す図である。
【図４６】ナビゲータシーケンスＮＳｆ＋１〜ＮＳｐを実行することにより算出された肝臓のエッジの移動量を概略的に示す図である。
【図４７】最後のイメージングシーケンスＤＡＱｚが実行されたときの様子を概略的に示す図である。
【図４８】ＬＳＱ法を用いて肝臓のエッジの移動量を求める利点の説明図である。
【発明を実施するための形態】
【００１０】
以下、発明を実施するための形態について説明するが、本発明は、以下の形態に限定されることはない。
【００１１】
（１）第１の形態
図１は、本発明の第１の形態の磁気共鳴装置の概略図である。
磁気共鳴装置（以下、「ＭＲ装置」と呼ぶ。ＭＲ：Magnetic Resonance）１００は、マグネット２、テーブル３、受信ＲＦコイル（以下、「受信コイル」と呼ぶ）４などを有している。
【００１２】
マグネット２は、被検体１３が収容されるボア２１を有している。また、マグネット２は、静磁場を発生させるための超伝導コイル、勾配磁場を印加するための勾配コイル、およびＲＦパルスを送信するためのＲＦコイルなどを有している。超伝導コイルの代わりに、永久磁石を用いてもよい。
【００１３】
テーブル３は、クレードル３ａを有している。クレードル３ａは、ボア２１内に移動できるように構成されている。クレードル３ａによって、被検体１３はボア２１に搬送される。
【００１４】
受信コイル４は、被検体１３の腹部に取り付けられている。受信コイル４は、被検体１３からの磁気共鳴信号を受信する。図２に受信コイル４の概略拡大図を示す。受信コイル４は、ＭＲ信号を受信するためのｎ個のチャネルＣＨ１〜ＣＨｎを有している。
図１に戻って説明を続ける。
【００１５】
ＭＲ装置１００は、更に、送信器５、勾配磁場電源６、受信器７、コンピュータ８、操作部１１、および表示部１２などを有している。
【００１６】
送信器５はＲＦコイルに電流を供給し、勾配磁場電源６は勾配コイルに電流を供給する。受信器７は、受信コイル４から受け取った信号に対して、検波などの信号処理を実行する。
【００１７】
コンピュータ８は、表示部１２に必要な情報を伝送したり、画像を再構成するなど、ＭＲ装置１００の各種の動作を実現するように、ＭＲ装置１００の各部の動作を制御する。コンピュータ８は、プロセッサ９およびメモリ１０などを有している。
【００１８】
図３に、プロセッサ９が実行する処理を示す。メモリ１０には、プロセッサ９により実行されるプログラムが記憶されている。プロセッサ９は、メモリ１０に記憶されているプログラムを読み出し、プログラムに記述されている処理を実行する。プロセッサ９は、メモリ１０に記憶されているプログラムを読み出すことにより、強度プロファイル作成手段９１〜判断手段９９などを構成する。
【００１９】
強度プロファイル作成手段９１は、受信コイルのチャネルごとに、ナビゲータ領域から得られたＭＲ信号の強度を表す強度情報を含む強度プロファイルを作成する。
位相プロファイル作成手段９２は、受信コイルのチャネルごとに、ナビゲータ領域から得られたＭＲ信号の位相を表す位相情報を含む位相プロファイルを作成する。
位置特定手段９３は、位相プロファイルに基づいて、強度プロファイルの中から、後述する範囲ｘ（図１９参照）を求めるときの基準となる基準位置を特定する。
範囲特定手段９４は、位置特定手段９３が特定した基準位置に基づいて、強度プロファイルの中から範囲ｘを特定する。
強度変更手段９５は、範囲特定手段により特定された範囲ｘ内の強度をゼロに変更する。
合成プロファイル作成手段９６は、各チャネルの強度プロファイルを合成する。
検出手段９７は、合成プロファイルに基づいて、肝臓のエッジを検出する。
ウインドウ設定手段９８は、後述するウインドウＷ（図２６参照）を設定する。
判断手段９９は、イメージングシーケンスを実行するか否かを判断する。
【００２０】
プロセッサ９は、強度プロファイル作成手段９１〜判断手段９９を構成する一例であり、メモリ１０に記憶されたプログラムを実行することによりこれらの手段として機能する。尚、受信器７、強度プロファイル作成手段９１、および位相プロファイル作成手段９２を合わせたものが、生成手段に相当する。
【００２１】
操作部１１は、オペレータにより操作され、種々の情報をコンピュータ８に入力する。表示部１２は種々の情報を表示する。
ＭＲ装置１００は、上記のように構成されている。
【００２２】
図４は第１の形態で実行されるスキャンを示す図である。
第１の形態では、ローカライザスキャンＬＸ、プレスキャンＡ、および本スキャンＢが実行される。
【００２３】
ローカライザスキャンＬＸは、撮影領域などを設定するために使用される画像を取得するためのスキャンである。
プレスキャンＡは、後述するウインドウＷ（図２６参照）求めるためのデータを取得するためのスキャンである。
本スキャンＢは、肝臓の画像を取得するためのスキャンである。
【００２４】
以下に、ローカライザスキャンＬＸ、プレスキャンＡ、および本スキャンＢを実行するときのフローについて説明する。
【００２５】
図５は、図４に示すスキャンを実行するときのフローを示す図である。
ステップＳＴ１では、ローカライザスキャンＬＸを実行する。図６に、ローカライザスキャンＬＸにより取得された画像ＬＤを概略的に示す。ローカライザスキャンＬＸを実行した後、ステップＳＴ２に進む。
【００２６】
ステップＳＴ２では、画像ＬＤに基づいて撮影領域Ｒが設定される。図７に、撮影領域を概略的に示す。オペレータは、操作部１１を操作し、画像ＬＤを参考にして、撮影領域Ｒを設定するために必要な情報を入力する。領域設定手段９０（図３参照）は、入力された情報に基づいて、撮影領域Ｒを設定する。第１の形態では、肝臓を撮影するので、肝臓を含む領域が撮影領域Ｒとして設定される。また、オペレータは、操作部１１を操作し、画像ＬＤを参考にして、肝臓のエッジの動きを検出するための領域（以下、「ナビゲータ領域」と呼ぶ）Ｒｎａｖを設定するために必要な情報も入力する。領域設定手段９０は、入力された情報に基づいて、ナビゲータ領域Ｒｎａｖを設定する。ナビゲータ領域ＲｎａｖのＳＩ方向の範囲は位置Ｐ１〜Ｐｚに設定されている。また、ナビゲータ領域Ｒｎａｖの中心位置Ｐｃは、肝臓の肺側のエッジ１３ａの近傍に位置決めされている。撮影領域Ｒおよびナビゲータ領域Ｒｎａｖを設定した後、ステップＳＴ３に進む。
【００２７】
ステップＳＴ３では、プレスキャンＡが実行される（図８参照）。
図８は、プレスキャンＡで実行されるシーケンスの説明図である。
プレスキャンＡでは、肝臓の動きを検出するためのナビゲータシーケンスＮＳ１〜ＮＳｆが順に実行される。以下、ナビゲータシーケンスＮＳ１〜ＮＳｆについて説明する。
【００２８】
プレスキャンＡでは、先ず、ナビゲータシーケンスＮＳ１が実行される。ナビゲータシーケンスＮＳ１が実行されることにより、ナビゲータ領域Ｒｎａｖ（図７参照）から、肝臓の動きを検出するためのＭＲ信号が発生する。このＭＲ信号は受信コイル４（図２参照）で受信される。受信コイル４は、図２に示すように、ｎ個のチャネルＣＨ１〜ＣＨｎを有しているので、ＭＲ信号は各チャネルで受信される。チャネルごとに受信されたＭＲ信号は受信器７（図１参照）に供給される。受信器７は、各チャネルから受け取った信号に対して検波などの信号処理を行い、信号処理により得られたデータを出力する。したがって、ナビゲータシーケンスＮＳ１を実行することにより、受信コイル４のチャネルごとに、肝臓の動きを検出するためのナビゲータデータを得ることができる。図９に、チャネルごとに得られたナビゲータデータＮＤ１〜ＮＤｎを概略的に示す。これらのナビゲータデータＮＤ１〜ＮＤｎは、コンピュータ８（図１参照）に供給される。
【００２９】
コンピュータ８では、受信器７から受け取った各チャネルのナビゲータデータＮＤ１〜ＮＤｎに基づいて、チャネルごとに、ナビゲータ領域Ｒｎａｖ内のＳＩ方向の信号強度の変化を表す強度プロファイルを作成する。図１０に、各チャネルのナビゲータデータＮＤ１〜ＮＤｎから得られた強度プロファイルＩＮ１〜ＩＮｎの一例を概略的に示す。強度プロファイルの横軸はＳＩ方向の位置を表し、縦軸は信号強度を表している。第１の形態では、ナビゲータシーケンスＮＳ１〜ＮＳｆは肝臓が高信号となるように設計されている。したがって、強度プロファイルＩＮ１〜ＩＮｎの肝臓側は高信号となる。一方、肺は空気を含んでいるので、肺側の信号はゼロに近い値にとなる。したがって、強度プロファイルＩＮ１〜ＩＮｎは、肺と肝臓との間に、信号強度の段差が現れる。
【００３０】
強度プロファイルＩＮ１〜ＩＮｎを作成した後、これらの強度プロファイルＩＮ１〜ＩＮｎを合成することにより、合成プロファイルを得ることができる。図１１に、強度プロファイルＩＮ１〜ＩＮｎを合成することにより得られた合成プロファイルＩＮを概略的に示す。合成方法としては、ＲＭＳ（Root Mean Square）などを用いることができる。
【００３１】
上記のように、強度プロファイルＩＮ１〜ＩＮｎは、肺と肝臓との間に信号強度の段差が現れる。したがって、合成プロファイルＩＮにも、肺と肝臓との間に信号強度の段差が現れる。この段差を検出することにより、肝臓のエッジ１３ａの位置Ｐｒを検出することができる。被検体が呼吸をすると、肝臓のエッジは動くので、ナビゲータシーケンスＮＳ１〜ＮＳｆを実行することにより、肝臓のエッジの位置の時間変化を表すデータを得ることができる。したがって、被検体の呼吸信号を得ることが可能となる。
【００３２】
しかし、実際には、ナビゲータシーケンスＮＳ１を実行した場合、Ｂ０不均一などが原因でナビゲータ領域Ｒｎａｖの外側も励起されてしまう。したがって、受信コイル４が受信するＭＲ信号には、ナビゲータ領域Ｒｎａｖに含まれる肝臓や肺からの信号だけでなく、ナビゲータ領域Ｒｎａｖの外側に存在する脂肪などの不要なＭＲ信号が混入する。また、ナビゲータ領域Ｒｎａｖに脂肪（内臓脂肪）が含まれている場合は、受信コイル４が受信するＭＲ信号には、ナビゲータ領域Ｒｎａｖの内側に存在する脂肪による不要なＭＲ信号も混入する。したがって、実際には、チャネルごとに得られた強度プロファイルＩＮ１〜ＩＮｎには、脂肪などの不要信号による強度が現れる（図１２参照）。
【００３３】
図１２は、不要信号の強度が現れた強度プロファイルの一例を概略的に示す図である。
強度プロファイルＩＮ１〜ＩＮｎの肝臓側は強度値が大きい。したがって、肝臓側に不要信号が混入しても、不要信号の強度は肝臓の信号値に埋もれてしまうので、肝臓側では不要信号の強度は現れにくい。
【００３４】
しかし、強度プロファイルＩＮ１〜ＩＮｎの肺側は空気の影響を受けるので、強度が小さくなる。したがって、肺側に不要信号が混入してしまうと、不要信号の強度が現れる。図１２では、強度プロファイルＩＮ１〜ＩＮｎの肺側に現れた不要信号の強度を、符号Ｃ１〜Ｃｎで示してある。したがって、これらの強度プロファイルＩＮ１〜ＩＮｎを合成すると、合成プロファイルＩＮにも、不要信号の強度Ｃが現れる。
【００３５】
図１２に示すように、肺側に不要信号の強度Ｃが現れてしまうと、合成プロファイルＩＮには、肝臓のエッジの実際の位置に現れる信号強度の段差Ｓ１と、不要信号の強度による信号強度の段差Ｓ２が現われる。したがって、信号強度の段差Ｓ２を、肝臓のエッジの位置に現れる信号強度の段差Ｓ１として誤検出してしまうことがある。そこで、本願発明者は、強度プロファイルに不要信号による強度が現れても、肝臓のエッジの実際の位置を検出することができるように、鋭意研究を重ね、強度プロファイルに現れる不要信号の強度Ｃ１〜Ｃｎを除去する方法を考え出した。以下に、不要信号の強度を除去する原理について簡単に説明する。
【００３６】
肺側に不要信号が混入すると、肺側の信号の位相が大きく変化する。本願発明者はこの点に着目し、位相が大きく変化する位置を特定することにより、不要信号が混入した位置を特定できることを見出した。したがって、位相情報を解析することにより、不要信号が混入した位置を特定することができるので、不要信号の強度Ｃ１〜Ｃｎを除去することが可能となり、エッジの検出精度を高めることができる。以下に、第１の形態におけるエッジの検出方法について具体的に説明する。尚、ナビゲータシーケンスＮＳ１〜ＮＳｆのうちのどのナビゲータシーケンスを実行した場合でも、エッジの検出方法は同じである。したがって、以下では、ナビゲータシーケンスＮＳ１〜ＮＳｆのうち、ナビゲータシーケンスＮＳ１を実行する場合を取り上げて、エッジの検出方法を説明する。
【００３７】
図１３は、第１の形態におけるエッジの検出方法のフローを示す図である。
ステップＳＴ１０では、プレスキャンＡにおけるナビゲータシーケンスＮＳ１を実行する。ナビゲータシーケンスＮＳ１を実行することにより、ナビゲータデータＮＤ１〜ＮＤｎが得られる（図９参照）。ナビゲータシーケンスＮＳ１を実行した後、ステップＳＴ２０に進む。
【００３８】
ステップＳＴ２０では、強度プロファイル作成手段９１（図３参照）が、各チャネルのナビゲータデータＮＤ１〜ＮＤｎに基づいて、チャネルごとに、ナビゲータ領域Ｒｎａｖ内のＳＩ方向の信号強度の変化を表す強度プロファイルＩＮ１〜ＩＮｎを作成する。図１４に、強度プロファイルＩＮ１〜ＩＮｎの一例を示す。強度プロファイルＩＮ１〜ＩＮｎを作成した後、ステップＳＴ３０に進む。
【００３９】
ステップＳＴ３０では、位相プロファイル作成手段９２（図３参照）が、各チャネルのナビゲータデータＮＤ１〜ＮＤｎに基づいて、チャネルごとに、ナビゲータ領域Ｒｎａｖ内のＳＩ方向の位相の変化を表す位相プロファイルを作成する。図１５に、チャネルごとに作成された位相プロファイルＰＨ１〜ＰＨｎの一例を示す。位相プロファイルの横軸はＳＩ方向の位置を表しており、縦軸は位相θｉを表している。また、位相プロファイルＰＨ１〜ＰＨｎの各々の下側には、任意の位置Ｐｉにおける位相θｉと、位置Ｐｉに隣接する位置Ｐｉ＋１における位相θｉ＋１が拡大して示してある。各位相プロファイルの位相θｉは、位相θｉとθｉ＋１との位相差Δθｉが−π≦Δθ≦πの範囲に収まるように、アンラップ処理されている。位相プロファイルＰＨ１〜ＰＨｎを作成した後、ステップＳＴ４０に進む。
【００４０】
ステップＳＴ４０では、位置特定手段９３（図３参照）が、位相プロファイルに基づいて、不要信号の位相が現れる位置を特定する。先に説明したように、肺側に不要信号が混入すると、肺側の信号の位相が大きく変化する。したがって、位相が大きく変化する位置を特定することにより、不要信号の位相が現れる位置を特定することができる。以下に、ステップＳＴ４０において、位相プロファイルから、不要信号の位相が現れる位置を特定する方法について説明する。尚、位相プロファイルＰＨ１〜ＰＨｎのうちのどの位相プロファイルであっても、不要信号の位相が現れる位置の特定方法は同じである。したがって、以下では、位相プロファイルＰＨ１〜ＰＨｎのうちの位相プロファイルＰＨ１を取り上げて、不要信号の位相が現れる位置の特定方法を説明する。
【００４１】
ステップＳＴ４０では、先ず、ステップＳＴ４１が実行される。
ステップＳＴ４１では、位置特定手段９３が、位相プロファイルＰＨ１を用いて、位置Ｐｉにおける位相差Δθｉを表す位相差プロファイルＤＩＦ１を作成する。図１６に、位相プロファイルＰＨ１と、位相プロファイルＰＨ１から作成された位相差プロファイルＤＩＦ１とを示す。図１６の右側には、各プロファイルを拡大して示してある。以下に、位相差プロファイルＤＩＦ１の作成方法について説明する。
【００４２】
位置特定手段９３は、以下の式を用いて位置Ｐｉにおける位相差Δθｉを算出する。
Δθｉ＝｜θｉ−θｉ＋１｜  ・・・（１）
ここで、θｉ：位相プロファイルＰＨ１の任意の位置Ｐｉにおける位相
θｉ＋１：位相プロファイルＰＨ１の位置Ｐｉに隣接する位置Ｐｉ＋１に
おける位相
ｉ：１〜ｚ−１の間の整数
【００４３】
式（１）によって、位置Ｐｉにおける位相差Δθｉを算出することができる。例えば、ｉ＝１、つまり、位置Ｐ１における位相差Δθ１を求める場合は、式（１）のｉに、ｉ＝１を代入すればよい。式（１）のｉに、ｉ＝１を代入すると、以下の式が得られる。
Δθ１＝｜θ１−θ２｜  ・・・（２）
【００４４】
したがって、位相プロファイルＰＨ１の位置Ｐ１における位相θ１から、位置Ｐ２における位相θ２を減算し、減算値の絶対値を求めることにより、位置Ｐ１における位相差Δθ１を算出することができる。
【００４５】
以下同様に、ｉの値を、ｉ＝２、３、４、・・・ｚ−１の順で変更することにより、各位置における位相差Δθｉを算出することができる。図１６では、ｉ＝ｋ、つまり、位置Ｐｋにおける位相差Δθｋと、ｉ＝ｚ−１、つまり、位置Ｐｚ−１における位相差Δθｚ−１も具体的に示してある。式（１）を用いることにより、位相差プロファイルＤＩＦ１を作成することができる。位相差プロファイルＤＩＦ１を作成した後、ステップＳＴ４２に進む。
【００４６】
ステップＳＴ４２では、位置特定手段９３が、位相差プロファイルＤＩＦ１の位置Ｐ１〜Ｐｚ−１の中から、位相差Δθｉが大きい値を持つ位置Ｐｉを特定する。位相差Δθｉが大きい値であるほど、位相が大きく変化していることを意味する。したがって、位相差Δθｉが大きい値を持つ位置Ｐｉを特定することにより、位相が大きく変化する位置Ｐｉを特定することができる。第１の形態では、位相差Δθｉが大きいか小さいかを判断するための閾値ＴＨを用いて、位相差Δθｉが大きい値を持つ位置Ｐｉを特定する。以下に、閾値ＴＨを用いて位相差Δθｉが大きい値を持つ位置Ｐｉを特定する方法について具体的に説明する。
【００４７】
図１７は、閾値ＴＨを用いて位相差Δθｉが大きい値を持つ位置Ｐｉを特定する方法の説明図である。
尚、閾値ＴＨは、脂肪などの不要信号の位相情報などを考慮して予め決められている値であり、例えば、ＴＨ＝１５°である。
【００４８】
位置特定手段９３は、先ず、位相差プロファイルＤＩＦ１のＳＩ方向の中心位置Ｐｃにおける位相差Δθｃを取り出す。そして、位置特定手段９３は、位相差Δθｃと閾値ＴＨとを比較し、位相差Δθｃが閾値ＴＨよりも大きいか否かを判断する。図１７には、位相差プロファイルＤＩＦ１の中心位置Ｐｃ付近の拡大図が示されている。
【００４９】
位相差Δθｃが閾値ＴＨよりも小さい場合（Δθｃ＜ＴＨ）、位置Ｐｃにおける位相差は小さいと判断する。一方、位相差Δθｃが閾値ＴＨ以上の場合（Δθｃ≧ＴＨ）、位置Ｐｃにおける位相差は大きいと判断する。図１７では、Δθｃ＜ＴＨであるので、位置特定手段９３は、位置Ｐｃにおける位相差Δθｃは小さいと判断する。位置Ｐｃにおける位相差Δθｃが小さいと判断されると、位置特定手段９３は、位相差プロファイルＤＩＦ１の位置Ｐ１〜Ｐｚ−１のうち、位置Ｐｃに近い位置から順に位相差を取り出し、位相差が大きいと判断されるまで、各位相差と閾値とを比較する。
【００５０】
具体的には、位置特定手段９３は、位置Ｐｃに対して肺側に隣接する位置Ｐｃ−１の位相差Δθｃ−１を取り出し、位相差Δθｃ−１と閾値ＴＨとを比較する。図１７に示すように、位相差Δθｃ−１は閾値ＴＨよりも小さいので、位置Ｐｃ−１における位相差Δθｃ−１は小さいと判断される。位相差Δθｃ−１が小さいと判断されると、位置特定手段９３は、位置Ｐｃに対して肝臓側に隣接する位置Ｐｃ＋１の位相差を取り出し、位相差Δθｃ＋１と閾値ＴＨとを比較する。図１７に示すように、位相差Δθｃ＋１は閾値ＴＨよりも小さいので、位置Ｐｃ＋１における位相差Δθｃ＋１は小さいと判断される。
【００５１】
以下同様に、位相差プロファイルＤＩＦｉの位置Ｐｉを、中心位置Ｐｃから近い順に、肺側、肝臓側に交互にずらしていき、位相差Δθｉが閾値ＴＨ以上になるまで、位相差Δθｉと閾値ＴＨとを比較する。図１７を参照すると、位置Ｐｃ−ｍ≦Ｐｉ≦Ｐｃ＋（ｍ−１）の範囲では、位相差Δθｉは、閾値ＴＨよりも小さい。しかし、位置Ｐｃ＋ｍにおける位相差θｃ＋ｍは閾値ＴＨを超えている。したがって、位置Ｐｃ＋ｍが、位相が大きく変化する位置として求められる。図１８に、強度プロファイルＩＮ１、位相差プロファイルＤＩＦ１、および位相プロファイルＰＨ１における位置Ｐｃ＋ｍを示す。
【００５２】
ステップＳＴ４０では、図１７を参照しながら説明したように、中心位置Ｐｃに近い位置から順に位相差Δθｉが閾値ＴＨを超えるか否かを判断している。中心位置Ｐｃは肝臓のエッジの近傍に位置決めされているので（図７参照）、中心位置Ｐｃに近い位置から順に位相差が閾値ＴＨを超えるか否かを判断することにより、肝臓のエッジ１３ａの近傍において、位相が大きく変化する位置Ｐｃ＋ｍを特定することができる。位置Ｐｃ＋ｍは、次のステップＳＴ５０において。後述する範囲ｘ（図１９参照）を求めるための基準位置として用いられる。位置Ｐｃ＋ｍを特定した後、ステップＳＴ５０に進む。
【００５３】
ステップＳＴ５０では、範囲特定手段９４（図３参照）が、位置Ｐｃ＋ｍに基づいて、強度プロファイルＩＮ１の全範囲の中から、不要信号による強度が現れている範囲を特定する。上記のように、位置Ｐｃ＋ｍは、不要信号により位相が大きく変化する位置を表している。したがって、強度プロファイルＩＮ１の位置Ｐｃ＋ｍに対して肺側の領域には、肺に混入した不要信号による強度が現れていると考えられる。そこで、範囲特定手段９４は、強度プロファイルＩＮ１の位置Ｐｃ＋ｍに対して肺側の範囲ｘ（位置Ｐ１〜Ｐｃ＋ｍの範囲）を、不要信号による強度が現れている範囲として特定する。図１９に範囲ｘを示す。範囲ｘを特定した後、ステップＳＴ６０に進む。
【００５４】
ステップＳＴ６０では、強度変更手段９５（図３参照）が、範囲ｘに含まれる強度をゼロに変更する。図２０に、範囲ｘの強度がゼロに変更された後の強度プロファイルＩＮ１´を示す。したがって、強度プロファイルＩＮ１から、肺側に混入した不要信号による強度を除去することができる。図２１に、範囲ｘの強度がゼロに変更される前の強度プロファイルＩＮ１と、範囲ｘの強度がゼロに変更された後の強度プロファイルＩＮ１´とを比較して示す。
【００５５】
尚、図１６〜図２０では、チャネルＣＨ１の強度プロファイルＩＮ１における不要信号の強度をゼロに変更する方法について示されているが、他のチャネルの強度プロファイルについても、同じ方法で不要信号の強度がゼロに変更される。したがって、全てのチャネルＣＨ１〜ＣＨｎの強度プロファイルＩＮ１〜ＩＮｎにおける不要信号の強度をゼロに変更することができる。図２２に、不要信号の強度がゼロに変更された後のチャネルＣＨ１〜ＣＨｎの強度プロファイルＩＮ１´〜ＩＮｎ´を示す。
不要信号の強度をゼロに変更した後、ステップＳＴ７０に進む。
【００５６】
ステップＳＴ７０では、合成プロファイル作成手段９６（図３参照）が、不要信号の強度がゼロに変更された強度プロファイルＩＮ１´〜ＩＮｎ´を合成する。図２３に、合成することにより得られた合成プロファイルＦ１を概略的に示す。合成プロファイルＦ１を得た後、ステップＳＴ８０に進む。
【００５７】
ステップＳＴ８０では、検出手段９７（図３参照）が、合成プロファイルＦ１の中から、肝臓のエッジの位置を検出する。図２４に、検出された肝臓のエッジの位置Ｐａ１を示す。検出手段９７は、例えば、合成プロファイルＦＩの強度が最大値Ｓ０の半分の値Ｓ０／２になるときの位置を、肝臓のエッジの位置Ｐａ１として検出することができる。
【００５８】
以下同様に、残りのナビゲータシーケンスＮＳ２〜ＮＳｆの各々を実行するときも、図１３に示すフローが実行される。したがって、ナビゲータシーケンスＮＳ１〜ＮＳｆの各々を実行したときの肝臓のエッジの位置を検出することができる。図２５に、ナビゲータシーケンスＮＳ１〜ＮＳｆの各々を実行したときの肝臓のエッジの位置の時間変化を概略的に示す。
【００５９】
したがって、プレスキャンＡを実行している間の被検体の呼吸信号Ｓｒｅｓを得ることができる。
プレスキャンＡを実行した後、ステップＳＴ４（図５参照）に進む。
【００６０】
ステップＳＴ４では、ウインドウ設定手段９８（図３参照）が、呼吸信号Ｓｒｅｓに基づいてウインドウＷを設定する。図２６に、ウインドウＷの一例を示す。ウインドウＷは、本スキャンＢにおいてイメージングシーケンスＤＡＱ１〜ＤＡＱｚ（図２７参照）を実行するときの基準となる肝臓のエッジ１３ａの範囲を表している。以下に、ウインドウＷを設定する方法の一例について簡単に説明する。
【００６１】
ウインドウ設定手段９８は、呼吸信号Ｓｒｅｓから、被検体が息を吐き終わったときの肝臓のエッジ１３ａの位置Ｐｂを特定する。被検体が息を吐いている間、エッジ１３ａの位置は増加するが、被検体が息を吸い始めると、エッジ１３ａの位置は減少し始める。したがって、エッジの位置が増加してから減少に転じるときの極大値を検出することにより、被検体が息を吐き終わったときのエッジ１３ａの位置Ｐｂを特定することができる。位置Ｐｂを特定した後、ウインドウ設定手段９８は、この位置Ｐｂに対して±ｘ％（例えば、ｘ＝２０）の範囲を、ウインドウＷとして設定する。ウインドウＷが本スキャンＢを実行するときにどのように使用されるかについては後述する。ウインドウＷを設定した後、ステップＳＴ５に進む。
【００６２】
ステップＳＴ５では、本スキャンＢが実行される。
図２７は、本スキャンＢの説明図である。
本スキャンＢでは、先ず、ナビゲータシーケンスＮＳｆ＋１〜ＮＳｐが実行される。ナビゲータシーケンスＮＳｆ＋１〜ＮＳｐの各々が実行されるたびに、図１３に示すフローに従って肝臓のエッジ１３ａの位置が検出される。図２７では、ナビゲータシーケンスＮＳｆ＋１〜ＮＳｐを実行することにより検出された肝臓のエッジ１３ａの位置が、符号「Ｐｆ＋１」、「Ｐｆ＋２」、・・・「Ｐｐ」で示されている。
【００６３】
判断手段９９（図３参照）は、ナビゲータシーケンスＮＳｆ＋１〜ＮＳｐを実行するたびに、検出された肝臓のエッジの位置がウインドウＷに入っているか否かを判断する。そして、肝臓のエッジ１３ａの位置がウインドウＷの外側からウインドウＷの内側に入り込んだときに、イメージングシーケンスＤＡＱ１を実行する。
【００６４】
図２７では、ナビゲータシーケンスＮＳｆ＋１〜ＮＳｐ−１では、肝臓のエッジ１３ａの位置はウインドウＷの外側であるが、ナビゲータシーケンスＮＳｐにおいて肝臓のエッジの位置はウインドウＷの内側に入り込む。したがって、ナビゲータシーケンスＮＳｐを実行した直後に、イメージングシーケンスＤＡＱ１が実行される。
【００６５】
イメージングシーケンスＤＡＱ１を実行した後、ナビゲータシーケンスＮＳｐ＋１〜ＮＳｑを実行し、肝臓のエッジの位置を検出する。そして、肝臓のエッジの位置がウインドウＷの外側からウインドウＷの内側に入り込んだときに、次のイメージングシーケンスＤＡＱ２を実行する。図２７では、ナビゲータシーケンスＮＳｐ＋１〜ＮＳｑ−１では、肝臓のエッジ１３ａの位置はウインドウＷの外側であるが、ナビゲータシーケンスＮＳｑにおいて肝臓のエッジ１３ａの位置はウインドウＷの内側に入り込む。したがって、ナビゲータシーケンスＮＳｑを実行した直後に、イメージングシーケンスＤＡＱ２が実行される。
【００６６】
以下同様にナビゲータシーケンスを実行し、肝臓のエッジ１３ａの位置がウインドウＷに入ったときに、イメージングシーケンスを実行する。そして、最後のイメージングシーケンスＤＡＱｚを実行したら本スキャンＢを終了し、フローを終了する。
【００６７】
第１の形態では、位相プロファイルの位相情報に基づいて、強度プロファイルに現れた不要信号の強度を除去することができるので、肝臓のエッジ１３ａの位置の検出精度を高めることができる。したがって、信頼性の高い呼吸信号を得ることができるので、体動アーチファクトが低減された高品質な画像を得ることができる。
【００６８】
第１の形態では、強度プロファイルの範囲ｘ（図２０参照）の強度をゼロに変更している。しかし、一般的には、範囲ｘの強度を十分に小さいに値に設定しておけば、肝臓のエッジの位置の検出精度には大きな問題はないと考えられる。したがって、強度プロファイルＩＮ１の範囲ｘの強度を、ゼロに変更する代わりに、例えば、１／Ｎの値に変更してもよい（Ｎは、例えばＮ＝１０）。図２８に、強度プロファイルＩＮ１の範囲ｘの強度を１／１０の値に変更したときの強度プロファイルＩＮ１´を概略的に示す。図２８では、範囲ｘの強度はゼロではないが、肝臓側の強度と比べると、十分に小さい値に設定されているので、肝臓のエッジの位置を十分な精度で検出することができる。
【００６９】
（２）第２の形態
第１の形態では、合成プロファイルに基づいて肝臓のエッジの位置を検出し、検出したエッジの位置を用いて呼吸信号を得る例について説明した。第２の形態では、肝臓のエッジの位置を検出する代わりに、肝臓のエッジの移動量を算出し、算出した移動量を用いて呼吸信号を得る例について説明する。
【００７０】
第２の形態のＭＲ装置は、第１の形態のＭＲ装置１００と比較すると、プロセッサが実行する処理に異なる点があるが、その他の構成は、第１の形態のＭＲ装置１００と同じである。したがって、第２の形態のＭＲ装置の説明に当たっては、プロセッサについて主に説明する。
【００７１】
図２９は、第２の形態におけるプロセッサ９が実行する処理を示す図である。
メモリ１０には、プロセッサ９により実行されるプログラムが記憶されている。プロセッサ９は、メモリ１０に記憶されているプログラムを読み出し、プログラムに記述されている処理を実行する。プロセッサ９は、メモリ１０に記憶されているプログラムを読み出すことにより、強度プロファイル作成手段９１〜判断手段９９などを構成する。
【００７２】
強度プロファイル作成手段９１〜合成プロファイル作成手段９６は、第１の形態と同じであるので、説明は省略する。
【００７３】
保存手段９６１は、後述する合成プロファイルＦ１（図３１参照）を保存する。
算出手段９６２は、合成プロファイルに基づいて肝臓のエッジの移動量を算出する。
ウインドウ設定手段９８は、後述するウインドウＷ（図４２参照）を設定する。
判断手段９９は、イメージングシーケンスを実行するか否かを判断する。
【００７４】
プロセッサ９は、強度プロファイル作成手段９１〜判断手段９９を構成する一例であり、メモリ１０に記憶されたプログラムを実行することによりこれらの手段として機能する。
【００７５】
次に、第２の形態において、図４のスキャン（ローカライザスキャンＬＸ、プレスキャンＡ、および本スキャンＢ）を実行するときのフローについて、図５を参照しながら説明する。
【００７６】
ステップＳＴ１およびＳＴ２は、第１の形態と同じであるので、説明は省略する。ステップＳＴ２において、撮影領域Ｒおよびナビゲータ領域Ｒｎａｖ（図７参照）を設定した後、ステップＳＴ３に進む。
【００７７】
ステップＳＴ３では、プレスキャンＡが実行される。プレスキャンＡでは、第１の形態と同様に、肝臓の動きを検出するためのナビゲータシーケンスＮＳ１〜ＮＳｆ（図８参照）が順に実行される。以下、ナビゲータシーケンスＮＳ１〜ＮＳｆについて説明する。
【００７８】
プレスキャンＡでは、先ず、ナビゲータシーケンスＮＳ１が実行される。第２の形態でも、第１の形態と同様に、ナビゲータシーケンスＮＳ１により得られたナビゲータデータに基づいて、合成プロファイルを作成する。図３０は、ナビゲータシーケンスＮＳ１を実行し、合成プロファイルを作成するときのフローを示す図である。
【００７９】
ステップＳＴ１０では、ナビゲータシーケンスＮＳ１が実行される。ナビゲータシーケンスＮＳ１を実行することにより、第１の形態で説明したように、受信コイルのチャネルごとにナビゲータデータＮＤ１〜ＮＤｎが得られる（図９参照）。ナビゲータシーケンスＮＳ１を実行した後、ステップＳＴ２０に進む。
【００８０】
ステップＳＴ２０〜ＳＴ７０は、第１の形態で説明した図１３のフローと同じであるので説明は省略する。ステップＳＴ２０〜ＳＴ７０を実行することにより、合成プロファイルが得られる。図３１に、ナビゲータシーケンスＮＳ１を実行することにより得られた合成プロファイルＦ１を概略的に示す。合成プロファイルＦ１を作成した後、ステップＳＴ７１に進む。
【００８１】
ステップＳＴ７１では、保存手段９６１（図２９参照）が合成プロファイルＦ１を保存する。合成プロファイルＦ１を保存したら、図３０のフローが終了する。したがって、第２の形態では、肝臓のエッジの位置を検出せずに図３０のフローを終了する。
【００８２】
ナビゲータシーケンスＮＳ１を実行した後、次のナビゲータシーケンスＮＳ２が実行される。図３２は、ナビゲータシーケンスＮＳ２を実行するときのフローを示す図である。以下に、図３２に示すフローについて説明する。
【００８３】
ステップＳＴ１０〜ＳＴ７０は、図３０と同じである。したがって、ナビゲータシーケンスＮＳ２を実行することにより、合成プロファイルＦ２が得られる。図３３に、ナビゲータシーケンスＮＳ２を実行することにより得られた合成プロファイルＦ２を概略的に示す。合成プロファイルＦ２を作成した後、ステップＳＴ７２に進む。
【００８４】
ステップＳＴ７２では、算出手段９６２（図２９参照）が、合成プロファイルＦ１およびＦ２に基づいて、肝臓のエッジの移動量を算出する。以下に肝臓のエッジの移動量の算出方法について説明する。
【００８５】
図３４は、合成プロファイルＦ１およびＦ２に基づいて肝臓のエッジの移動量を算出するときの説明図である。
【００８６】
図３４の上段には、ナビゲータシーケンスＮＳ１を実行したときの肝臓が破線で示してあり、ナビゲータシーケンスＮＳ２を実行したときの肝臓が実線で示してある。
図３４の中段には、ナビゲータシーケンスＮＳ１を実行することにより得られた合成プロファイルＦ１が示されており、図３４の下段には、ナビゲータシーケンスＮＳ２を実行することにより得られた合成プロファイルＦ２が示されている。
【００８７】
ナビゲータシーケンスＮＳ１を実行してからナビゲータシーケンスＮＳ２が実行されるまでの間、肝臓のエッジは移動する。したがって、合成プロファイルＦ１の立上り位置Ｐｒ１と、合成プロファイルＦ１の立上り位置Ｐｒ２には、肝臓のエッジが移動したことによる位置ずれ量ｄが生じる。この位置ずれ量ｄは、肝臓のエッジの移動量と考えることができるので、ｄを算出することができれば、肝臓のエッジの移動量を求めることができる。以下に、ｄを算出する方法について具体的に説明する。尚、第２の形態では、ｄを算出する方法の一例として、ＬＳＱ（Least Square）を用いた方法について説明する。
【００８８】
先ず、算出手段９６２は、合成プロファイルＦ１とＦ２との二乗誤差ＳＥ（Square Error）を算出する。図３５に、算出された二乗誤差ＳＥ＝ＳＥ０を概略的に示す。
【００８９】
次に、算出手段９６２は、合成プロファイルＦ２をＩ方向にΔｄだけ移動させ、二乗誤差を求める（図３６参照）。
【００９０】
図３６は、合成プロファイルＦ２をＩ方向に移動させ、二乗誤差を算出するときの説明図である。
算出手段９６２は、合成プロファイルＦ２をＩ方向にΔｄ＝Δｄ１、Δｄ２、・・・Δｄｈだけ移動させ、Δｄごとに、二乗誤差ＳＥ＝ＳＥ１、ＳＥ２、・・・ＳＥｈを算出する。
【００９１】
二乗誤差ＳＥ＝ＳＥ１、ＳＥ２、・・・ＳＥｈを算出した後、算出手段９６２は、合成プロファイルＦ２をＳ方向に移動させ、二乗誤差を算出する（図３７参照）。
【００９２】
図３７は、合成プロファイルＦ２をＳ方向に移動させ、二乗誤差を算出するときの説明図である。
尚、図３７に示す合成プロファイルＦ２の移動方向（Ｓ方向）は、図３６に示す合成プロファイルＦ１の移動方向（Ｉ方向）とは反対方向である。そこで、図３７に示す合成プロファイルＦ２の移動量Δｄは、負の値で表してある。
【００９３】
算出手段９６２は、合成プロファイルＦ２をＳ方向にΔｄ＝−Δｄ１、−Δｄ２、・・・−Δｄｈだけ移動させ、Δｄごとに、二乗誤差ＳＥ＝ＳＥ−１、ＳＥ−２、・・・ＳＥ−ｈを算出する。
【００９４】
したがって、Δｄの値に応じて二乗誤差ＳＥ０、ＳＥ１、ＳＥ２、・・・ＳＥｈ、ＳＥ−１、ＳＥ−２、・・・ＳＥ−ｈを算出することができる。
【００９５】
合成プロファイルＦ１の立上り位置Ｐｒ１と合成プロファイルＦ２の立上り位置Ｐｒ２との間のずれが大きい場合、二乗誤差ＳＥは大きくなる傾向がある。一方、合成プロファイルＦ１の立上り位置Ｐｒ１と合成プロファイルＦ２の立上り位置Ｐｒ２との間のずれが小さい場合、二乗誤差ＳＥは小さくなる傾向がある。したがって、二乗誤差ＳＥが最小値になるときのΔｄは、肝臓エッジ１３ａの移動量ｄに十分に近い値を有すると考えることができる。そこで、算出手段９６２は、二乗誤差ＳＥ０、ＳＥ１、ＳＥ２、・・・ＳＥｈ、ＳＥ−１、ＳＥ−２、・・・ＳＥ−ｈの中から最小値ＳＥｍｉｎを特定し、二乗誤差ＳＥ＝ＳＥｍｉｎになるときのΔｄを、肝臓のエッジ１３ａの移動量ｄとして算出する。ここでは、二乗誤差ＳＥ２（図３６参照）が最小値ＳＥｍｉｎであるとする。二乗誤差ＳＥ２はΔｄ＝Δｄ２のときに得られた値であるので、肝臓のエッジ１３ａの移動量ｄはｄ＝Δｄ２と求められる。図３８に、合成プロファイＦ１およびＦ２を用いて算出された肝臓のエッジの移動量ｄ＝Δｄ２を示す。
【００９６】
したがって、ナビゲータシーケンスＮＳ１における肝臓のエッジの位置をｄ＝０の位置で表すとすると、ナビゲータシーケンスＮＳ２における肝臓のエッジの位置は、ｄ＝０に対してΔｄ２だけずれていることがわかる。
【００９７】
ナビゲータシーケンスＮＳ２を実行した後、次のナビゲータシーケンスＮＳ３が実行される。ナビゲータシーケンスＮＳ３を実行した場合も、図３２に示したフローが実行される。したがって、ナビゲータシーケンスＮＳ３を実行することにより、合成プロファイルＦ３が作成される。図３９に、ナビゲータシーケンスＮＳ３を実行することにより得られた合成プロファイルＦ３を概略的に示す。合成プロファイルＦ３を求めた後、算出手段９６２は、合成プロファイルＦ１と合成プロファイルＦ３とを用いて、肝臓のエッジの移動量ｄを算出する。図４０に、算出された移動量ｄ＝Δｄ３を示す。
【００９８】
したがって、ナビゲータシーケンスＮＳ３における肝臓のエッジの位置は、ｄ＝０に対してΔｄ３だけずれていることがわかる。
【００９９】
以下同様に、残りのナビゲータシーケンスＮＳ４〜ＮＳｆを実行した場合も、図３２に示すフローが実行される。したがって、ナビゲータシーケンスを実行するたびに肝臓のエッジの移動量が求められる。図４１に、ナビゲータシーケンスＮＳ１〜ＮＳｆを実行することにより求められた肝臓のエッジの移動量ｄを概略的に示す。ｄを求めることにより、ナビゲータシーケンスＮＳ２〜ＮＳｆにおける肝臓のエッジの位置が、ｄ＝０に対してどれだけずれているかを知ることができる。したがって、プレスキャンＡを実行している間の被検体の呼吸信号Ｓｒｅｓを得ることができる。
プレスキャンＡを実行した後、ステップＳＴ４（図５参照）に進む。
【０１００】
ステップＳＴ４では、ウインドウ設定手段９８（図２９参照）が、呼吸信号Ｓｒｅｓに基づいてウインドウＷを設定する。図４２に、ウインドウＷの一例を示す。以下に、ウインドウＷを設定する方法の一例について簡単に説明する。
【０１０１】
ウインドウ設定手段９８は、呼吸信号Ｓｒｅｓから、被検体が息を吐き終わったときの肝臓のエッジの移動量Δｄｘを特定する。被検体が息を吐いている間、移動量ｄは増加するが、被検体が息を吸い始めると、移動量ｄは減少し始める。したがって、移動量ｄが増加してから減少に転じるときの極大値を検出することにより、被検体が息を吐き終わったときの肝臓のエッジの移動量Δｄ０を求めることができる。移動量Δｄｘを求めた後、ウインドウ設定手段９８は、移動量Δｄｘに基づいて、一定の範囲をウインドウＷとして設定する。
ウインドウＷを設定した後、ステップＳＴ５に進む。
【０１０２】
ステップＳＴ５では、本スキャンＢが実行される。図４３に、本スキャンＢで実行されるナビゲータシーケンスＮＳｆ＋１〜ＮＳｕおよびイメージングシーケンスＤＡＱ１〜ＤＡＱｚを概略的に示す。本スキャンＢでは、プレスキャンＡと同様に、図３２に示すフローに従って肝臓のエッジの移動量を算出する。以下に、本スキャンＢにおける肝臓のエッジの移動量の算出方法について、図３２のフローを参照しながら説明する。
【０１０３】
ステップＳＴ１０では、ナビゲータシーケンスＮＳｆ＋１が実行される。ナビゲータシーケンスＮＳｆ＋１を実行した後、ステップＳＴ２０〜ＳＴ７０を実行する。ステップＳＴ２０〜ＳＴ７０は第１の形態で説明したので、ステップＳＴ２０〜ＳＴ７０の説明は省略する。ステップＳＴ７０を実行することにより、合成プロファイルＦｆ＋１が得られる。図４４には、ナビゲータシーケンスＮＳｆ＋１を実行することにより得られた合成プロファイルＦｆ＋１を概略的に示してある。合成プロファイルＦｆ＋１を作成した後、ステップＳＴ７２に進む。
【０１０４】
ステップＳＴ７２では、算出手段９６２が、肝臓のエッジの移動量を算出する。図４５に、肝臓のエッジの移動量を算出するときの方法を示す。肝臓のエッジの移動量は、第１の形態と同様に、プレスキャンＡで求められた合成プロファイルＦ１（図３１参照）を用いて算出される。算出手段９６２は、合成プロファイルＦｆ＋１をＳＩ方向にΔｄだけ移動し、Δｄの値ごとに、合成プロファイルＦ１と合成プロファイルＦｆ＋１との二乗誤差ＳＥを算出する。そして、二乗誤差ＳＥが最小になるときのΔｄを求める。ここでは、Δｄ＝−Δｄｋであるとする。したがって、肝臓のエッジ１３ａの移動量ｄ＝−Δｄｋが求められる。
【０１０５】
以下同様に、ナビゲータシーケンスを実行し、肝臓のエッジの移動量を算出する。図４６に、ナビゲータシーケンスＮＳｆ＋１〜ＮＳｐを実行することにより算出された肝臓のエッジの移動量を概略的に示す。図４６では、ナビゲータシーケンスＮＳｆ＋１〜ＮＳｐ−１では、肝臓のエッジ１３ａの移動量はウインドウＷの外側であるが、ナビゲータシーケンスＮＳｐにおいて肝臓のエッジの移動量はウインドウＷの内側に入り込む。したがって、ナビゲータシーケンスＮＳｐを実行した直後に、イメージングシーケンスＤＡＱ１が実行される。
【０１０６】
イメージングシーケンスＤＡＱ１を実行した後、再びナビゲータシーケンスを実行し、肝臓のエッジの移動量がウインドウＷに入ったときに、イメージングシーケンスを実行する。図４７に、最後のイメージングシーケンスＤＡＱｚが実行されたときの様子を概略的に示す。最後のイメージングシーケンスＤＡＱｚが実行されたら本スキャンＢを終了する。
【０１０７】
第２の形態でも、第１の形態と同様に、位相プロファイルの位相情報に基づいて、強度プロファイルに現れた不要信号の強度を除去する。したがって、肝臓のエッジの移動量の計算精度を高めることができるので、信頼性の高い呼吸信号を得ることができ、体動アーチファクトが低減された高品質な画像を得ることができる。
【０１０８】
また、第１の形態では、肝臓のエッジの位置を検出しているが、第２の形態では、肝臓のエッジの位置を検出せずに、ＬＳＱ法を用いて肝臓のエッジの移動量を求めている。以下に、肝臓のエッジの位置を検出せずに、ＬＳＱ法を用いて肝臓のエッジの移動量を求める利点について説明する。
【０１０９】
図４８は、ＬＳＱ法を用いて肝臓のエッジの移動量を求める利点の説明図である。
図４８には、２つの合成プロファイルＦ１およびＦ２が示されている。合成プロファイルＦ１は、肺側の不要信号の強度が完全に除去されているが、合成プロファイルＦ２は、一部の不要信号の強度Ｃが除去されずに残った例が示されている。
【０１１０】
合成プロファイルの肺側の強度と肝臓側の強度とを比較すると、肝臓側の強度が大きい。したがって、合成プロファイルＦ１とＦ２との二乗誤差を求めた場合、二乗誤差の値は、肺側の強度よりも、肝臓側の強度の影響を強く受ける。したがって、肺側に不要信号の強度Ｃが残っていても、不要信号の強度Ｃが二乗誤差の値に与える影響は小さいと考えられる。したがって、二乗誤差を算出し、肝臓のエッジの移動量ｄを求めることにより、より信頼性の高い呼吸信号を得ることができる。
【０１１１】
尚、第２の形態では、二乗誤差を用いて肝臓のエッジの移動量を算出しているが、二乗誤差の代わりに、例えば相関（correlation）を用いて肝臓のエッジの移動量を算出してもよい。
【０１１２】
第１および第２の形態では、ナビゲータ領域Ｒｎａｖは肝臓と肺との境界に設定されているが、本発明は、ナビゲータ領域Ｒｎａｖを別の部位に設定した場合にも適用することができる。
【０１１３】
第１および第２の形態では、呼吸同期法を例にとって本発明について説明したが、本発明は、強度プロファイルに現れる不要信号の強度を除去したい場合であれば、呼吸同期法を用いない撮影方法にも適用することができる。
【０１１４】
また、第１および第２の形態では、マルチチャネルの受信コイルを用いた例について説明されているが、本発明は、シングルチャネルの受信コイルにも適用することができる。
【符号の説明】
【０１１５】
２  マグネット
３  テーブル
３ａ  クレードル
４  受信コイル
５  送信器
６  勾配磁場電源
７  受信器
８  コンピュータ
９  プロセッサ
１０  メモリ
１１  操作部
１２  表示部
１３  被検体
２１  ボア
９０  領域設定手段
９１  強度プロファイル作成手段
９２  位相プロファイル作成手段
９３  位置特定手段
９４  範囲特定手段
９５  強度変更手段
９６  合成プロファイル作成手段
９７  検出手段
９８  ウインドウ設定手段
９９  判断手段
１００  ＭＲ装置
９６１  保存手段
９６２  算出手段
【図１】
2016010582.tif 000003
【図２】
2016010582.tif 000004
【図３】
2016010582.tif 000005
【図４】
2016010582.tif 000006
【図５】
2016010582.tif 000007
【図６】
2016010582.tif 000008
【図７】
2016010582.tif 000009
【図８】
2016010582.tif 000010
【図９】
2016010582.tif 000011
【図１０】
2016010582.tif 000012
【図１１】
2016010582.tif 000013
【図１２】
2016010582.tif 000014
【図１３】
2016010582.tif 000015
【図１４】
2016010582.tif 000016
【図１５】
2016010582.tif 000017
【図１６】
2016010582.tif 000018
【図１７】
2016010582.tif 000019
【図１８】
2016010582.tif 000020
【図１９】
2016010582.tif 000021
【図２０】
2016010582.tif 000022
【図２１】
2016010582.tif 000023
【図２２】
2016010582.tif 000024
【図２３】
2016010582.tif 000025
【図２４】
2016010582.tif 000026
【図２５】
2016010582.tif 000027
【図２６】
2016010582.tif 000028
【図２７】
2016010582.tif 000029
【図２８】
2016010582.tif 000030
【図２９】
2016010582.tif 000031
【図３０】
2016010582.tif 000032
【図３１】
2016010582.tif 000033
【図３２】
2016010582.tif 000034
【図３３】
2016010582.tif 000035
【図３４】
2016010582.tif 000036
【図３５】
2016010582.tif 000037
【図３６】
2016010582.tif 000038
【図３７】
2016010582.tif 000039
【図３８】
2016010582.tif 000040
【図３９】
2016010582.tif 000041
【図４０】
2016010582.tif 000042
【図４１】
2016010582.tif 000043
【図４２】
2016010582.tif 000044
【図４３】
2016010582.tif 000045
【図４４】
2016010582.tif 000046
【図４５】
2016010582.tif 000047
【図４６】
2016010582.tif 000048
【図４７】
2016010582.tif 000049
【図４８】
2016010582.tif 000050