(57)【要約】      （修正有）【課題】磁気共鳴装置による呼吸同期法を用いた腹部撮影において、呼吸信号に対するサイドローブの影響を低減する。【解決手段】｜Ｗ´（Ｋ´）｜Ｓ´（Ｋ´）＝｜Ｗ（Ｋ）｜Ｓ（Ｋ）およびθ´＝−θを満たすように設定されたナビゲータシーケンスＮＳ１およびＮＳ２を実行するＳＴ１。ナビゲータシーケンスＮＳ１により収集されたエコーからナビゲータデータを得るとともに、ナビゲータシーケンスＮＳ２により得られたエコーからナビゲータデータを得る。そして、ナビゲータデータを加算しＳＴ２、加算ナビゲータデータに基づいて、肝臓のエッジの位置を検出するＳＴ３。【選択図】図１１
2016049298.tif 000002
【特許請求の範囲】
【請求項１】
第１の励起ｋ空間において座標点（ｋｘ，ｋｙ）を含む第１のｋ空間軌跡を辿るように被検体の第１の部位を励起し、前記第１の部位の情報を含む第１のデータを取得するための第１のシーケンスであって、前記第１の部位を励起するための第１の励起パルスと、第１の軸に印加される第１の傾斜磁場と、第２の軸に印加される第２の傾斜磁場とを有し、前記第１の傾斜磁場および前記第２の傾斜磁場は、前記第１のｋ空間軌跡が得られるように構成されている第１のシーケンス、並びに
第２の励起ｋ空間において座標点（ｋｘ´，ｋｙ´）を含む第２のｋ空間軌跡を辿るように前記第１の部位を励起し、前記第１の部位の情報を含む第２のデータを取得するための第２のシーケンスであって、前記第１の部位を励起するための第２の励起パルスと、前記第１の軸に印加される第３の傾斜磁場と、前記第２の軸に印加される第４の傾斜磁場とを有し、前記第３の傾斜磁場と前記第４の傾斜磁場は、前記座標点（ｋｘ´，ｋｙ´）が前記座標点（ｋｘ，ｋｙ）とは符号が反対の座標点（−ｋｘ，−ｋｙ）に設定された前記第２のｋ空間軌跡が得られるように構成されている第２のシーケンス、を実行するためのスキャン手段と、
前記第１のデータと前記第２のデータとを合成する合成手段と、
を備えており、
前記第１のシーケンスおよび前記第２のシーケンスは、以下の式を満たすように設定されている、磁気共鳴装置。
【数２０】
2016049298.tif 000022
ここで、Ｂ１：前記（ｋｘ，ｋｙ）における前記第１の励起パルスの値
｜ｇ｜：前記（ｋｘ，ｋｙ）における前記第１の傾斜磁場と前記第２の傾斜磁場との合成傾斜磁場の値の絶対値
Ｓ（Ｋ）：前記（ｋｘ，ｋｙ）におけるサンプリング密度関数
Ｂ１´：前記（ｋｘ´，ｋｙ´）における前記第２の励起パルスの値
｜ｇ´｜：前記（ｋｘ´，ｋｙ´）における前記第３の傾斜磁場と前記第４の傾斜磁場との合成傾斜磁場の値の絶対値
Ｓ´（Ｋ´）：前記（ｋｘ´，ｋｙ´）におけるサンプリング密度関数
γ：磁気回転比
θ：前記（ｋｘ，ｋｙ）における前記第１の励起パルスの位相
θ´：前記（ｋｘ´，ｋｙ´）における前記第２の励起パルスの位相
【請求項２】
｜ｇ｜および｜ｇ´｜は、以下の式で表される、請求項１に記載の磁気共鳴装置。
【数２１】
2016049298.tif 000023
ここで、ｇｘ：前記（ｋｘ，ｋｙ）における前記第１の傾斜磁場の値
ｇｙ：前記（ｋｘ，ｋｙ）における前記第２の傾斜磁場の値
ｇｘ´：前記（ｋｘ´，ｋｙ´）における前記第３の傾斜磁場の値
ｇｙ´：前記（ｋｘ´，ｋｙ´）における前記第４の傾斜磁場の値
【請求項３】
｜ｇ｜および｜ｇ´｜は、以下の関係を有している、請求項２に記載の磁気共鳴装置。
｜ｇ´｜＝｜ｇ｜
【請求項４】
Ｂ１およびＢ１´は、以下の式で表される、請求項１〜３のうちのいずれか一項に記載の磁気共鳴装置。
Ｂ１＝ｂ・ｅｉθ
Ｂ１´＝ｂ´・ｅｉθ´
ここで、ｂ：前記（ｋｘ，ｋｙ）における前記第１の励起パルスの大きさ
ｂ´：前記（ｋｘ´，ｋｙ´）における前記第２の励起パルスの大きさ
【請求項５】
ｂ＝ｂ´である、請求項４に記載の磁気共鳴装置。
【請求項６】
前記第２の励起パルスの励起時間は、前記第１の励起パルスの励起時間とは異なっている、請求項１〜５のうちのいずれか一項に記載の磁気共鳴装置。
【請求項７】
前記第１のｋ空間軌跡および前記第２のｋ空間軌跡はスパイラル軌跡である、請求項１〜６のうちのいずれか一項に記載の磁気共鳴装置。
【請求項８】
前記スパイラル軌跡の始点は、ｋ空間の原点とは異なる座標点に設定されており、前記スパイラル軌跡の終点は、ｋ空間の原点に設定されている、請求項７に記載の磁気共鳴装置。
【請求項９】
前記スパイラル軌跡の始点は、ｋ空間の原点に設定されており、前記スパイラル軌跡の終点は、ｋ空間の原点とは異なる座標点に設定されている、請求項７に記載の磁気共鳴装置。
【請求項１０】
前記サンプリング密度関数Ｓ（Ｋ）およびＳ´（Ｋ´）は、以下の関係を有している、請求項１〜９のうちのいずれか一項に記載の磁気共鳴装置。
Ｓ´（Ｋ´）＝Ｓ（Ｋ）
【請求項１１】
前記スキャン手段は、
前記被検体の呼吸信号を取得するためのプリスキャンと、前記被検体の第２の部位の画像を取得するための本スキャンとを実行し、
前記プリスキャンでは、前記第１のシーケンスおよび前記第２のシーケンスのセットが実行され、
前記本スキャンでは、前記第１のシーケンスおよび前記第２のシーケンスのセットと、前記第２の部位の画像を取得するための第３のシーケンスとが実行される、請求項１〜１０のうちのいずれか一項に記載の磁気共鳴装置。
【請求項１２】
前記第１のデータおよび前記第２のデータは、前記第１の部位の位置情報を含むナビゲータデータである、請求項１〜１１のうちのいずれか一項に記載の磁気共鳴装置。
【請求項１３】
前記第１のデータおよび前記第２のデータは、前記第１の部位のイメージングデータである、請求項１〜１１のうちのいずれか一項に記載の磁気共鳴装置。
【請求項１４】
第１の励起ｋ空間において座標点（ｋｘ，ｋｙ）を含む第１のｋ空間軌跡を辿るように被検体の第１の部位を励起し、前記第１の部位の情報を含む第１のデータを取得するための第１のシーケンスであって、前記第１の部位を励起するための第１の励起パルスと、第１の軸に印加される第１の傾斜磁場と、第２の軸に印加される第２の傾斜磁場とを有し、前記第１の傾斜磁場および前記第２の傾斜磁場は、前記第１のｋ空間軌跡が得られるように構成されている第１のシーケンス、並びに
第２の励起ｋ空間において座標点（ｋｘ´，ｋｙ´）を含む第２のｋ空間軌跡を辿るように前記第１の部位を励起し、前記第１の部位の情報を含む第２のデータを取得するための第２のシーケンスであって、前記第１の部位を励起するための第２の励起パルスと、前記第１の軸に印加される第３の傾斜磁場と、前記第２の軸に印加される第４の傾斜磁場とを有し、前記第３の傾斜磁場と前記第４の傾斜磁場は、前記座標点（ｋｘ´，ｋｙ´）が前記座標点（ｋｘ，ｋｙ）とは符号が反対の座標点（−ｋｘ，−ｋｙ）に設定された前記第２のｋ空間軌跡が得られるように構成されている第２のシーケンス、を実行するためのスキャン手段であって、以下の式を満たすように設定された前記第１のシーケンスおよび前記第２のシーケンスを実行するスキャン手段を有する磁気共鳴装置に適用されるプログラムであって、
前記第１のデータと前記第２のデータとを合成する合成処理をコンピュータに実行させるためのプログラム。
【数２２】
2016049298.tif 000024
ここで、Ｂ１：前記（ｋｘ，ｋｙ）における前記第１の励起パルスの値
｜ｇ｜：前記（ｋｘ，ｋｙ）における前記第１の傾斜磁場と前記第２の傾斜磁場との合成傾斜磁場の値の絶対値
Ｓ（Ｋ）：前記（ｋｘ，ｋｙ）におけるサンプリング密度関数
Ｂ１´：前記（ｋｘ´，ｋｙ´）における前記第２の励起パルスの値
｜ｇ´｜：前記（ｋｘ´，ｋｙ´）における前記第３の傾斜磁場と前記第４の傾斜磁場との合成傾斜磁場の値の絶対値
Ｓ´（Ｋ´）：前記（ｋｘ´，ｋｙ´）におけるサンプリング密度関数
γ：磁気回転比
θ：前記（ｋｘ，ｋｙ）における前記第１の励起パルスの位相
θ´：前記（ｋｘ´，ｋｙ´）における前記第２の励起パルスの位相
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被検体をスキャンする磁気共鳴装置、およびその磁気共鳴装置に適用されるプログラムに関する。
【背景技術】
【０００２】
腹部の撮影では、動きアーチファクトを低減するために呼吸同期法を用いた撮影方法が実行される。呼吸同期法の撮影方法として、肺と肝臓との境界部分をペンシルビーム型に選択励起し、肝臓のエッジの動きを検出する方法が知られている（特許文献１参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００９−２６１５７４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかし、ペンシルビーム型の選択励起では、肺と肝臓との境界部分で、大きい横磁化を持つメインローブが現れるが、メインローブから離れた位置に、小さい横磁化を持つサイドローブが現れる。したがって、このサイドローブによって、肺と肝臓との境界から離れた部位（例えば、体表面）でも励起が行われる。このように、サイドローブによる励起が行われると、不要な信号がナビゲータエコーに混入するので、肝臓のエッジの検出誤差が大きくなるという問題がある。
【０００５】
したがって、サイドローブの影響を低減することができる技術が望まれている。
【課題を解決するための手段】
【０００６】
本発明の第１の観点は、第１の励起ｋ空間において座標点（ｋｘ，ｋｙ）を含む第１のｋ空間軌跡を辿るように被検体の第１の部位を励起し、前記第１の部位の情報を含む第１のデータを取得するための第１のシーケンスであって、前記第１の部位を励起するための第１の励起パルスと、第１の軸に印加される第１の傾斜磁場と、第２の軸に印加される第２の傾斜磁場とを有し、前記第１の傾斜磁場および前記第２の傾斜磁場は、前記第１のｋ空間軌跡が得られるように構成されている第１のシーケンス、並びに
第２の励起ｋ空間において座標点（ｋｘ´，ｋｙ´）を含む第２のｋ空間軌跡を辿るように前記第１の部位を励起し、前記第１の部位の情報を含む第２のデータを取得するための第２のシーケンスであって、前記第１の部位を励起するための第２の励起パルスと、前記第１の軸に印加される第３の傾斜磁場と、前記第２の軸に印加される第４の傾斜磁場とを有し、前記第３の傾斜磁場と前記第４の傾斜磁場は、前記座標点（ｋｘ´，ｋｙ´）が前記座標点（ｋｘ，ｋｙ）とは符号が反対の座標点（−ｋｘ，−ｋｙ）に設定された前記第２のｋ空間軌跡が得られるように構成されている第２のシーケンス、を実行するためのスキャン手段と、
前記第１のデータと前記第２のデータとを合成する合成手段と、
を備えており、
前記第１のシーケンスおよび前記第２のシーケンスは、以下の式を満たすように設定されている、磁気共鳴装置。
【数１】
2016049298.tif 000003
ここで、Ｂ１：前記（ｋｘ，ｋｙ）における前記第１の励起パルスの値
｜ｇ｜：前記（ｋｘ，ｋｙ）における前記第１の傾斜磁場と前記第２の傾斜磁場との合成傾斜磁場の値の絶対値
Ｓ（Ｋ）：前記（ｋｘ，ｋｙ）におけるサンプリング密度関数
Ｂ１´：前記（ｋｘ´，ｋｙ´）における前記第２の励起パルスの値
｜ｇ´｜：前記（ｋｘ´，ｋｙ´）における前記第３の傾斜磁場と前記第４の傾斜磁場との合成傾斜磁場の値の絶対値
Ｓ´（Ｋ´）：前記（ｋｘ´，ｋｙ´）におけるサンプリング密度関数
γ：磁気回転比
θ：前記（ｋｘ，ｋｙ）における前記第１の励起パルスの位相
θ´：前記（ｋｘ´，ｋｙ´）における前記第２の励起パルスの位相
【０００７】
本発明の第２の観点は、第１の励起ｋ空間において座標点（ｋｘ，ｋｙ）を含む第１のｋ空間軌跡を辿るように被検体の第１の部位を励起し、前記第１の部位の情報を含む第１のデータを取得するための第１のシーケンスであって、前記第１の部位を励起するための第１の励起パルスと、第１の軸に印加される第１の傾斜磁場と、第２の軸に印加される第２の傾斜磁場とを有し、前記第１の傾斜磁場および前記第２の傾斜磁場は、前記第１のｋ空間軌跡が得られるように構成されている第１のシーケンス、並びに
第２の励起ｋ空間において座標点（ｋｘ´，ｋｙ´）を含む第２のｋ空間軌跡を辿るように前記第１の部位を励起し、前記第１の部位の情報を含む第２のデータを取得するための第２のシーケンスであって、前記第１の部位を励起するための第２の励起パルスと、前記第１の軸に印加される第３の傾斜磁場と、前記第２の軸に印加される第４の傾斜磁場とを有し、前記第３の傾斜磁場と前記第４の傾斜磁場は、前記座標点（ｋｘ´，ｋｙ´）が前記座標点（ｋｘ，ｋｙ）とは符号が反対の座標点（−ｋｘ，−ｋｙ）に設定された前記第２のｋ空間軌跡が得られるように構成されている第２のシーケンス、を実行するためのスキャン手段であって、以下の式を満たすように設定された前記第１のシーケンスおよび前記第２のシーケンスを実行するスキャン手段を有する磁気共鳴装置に適用されるプログラムであって、
前記第１のデータと前記第２のデータとを合成する合成処理をコンピュータに実行させるためのプログラム。
【数２】
2016049298.tif 000004
ここで、Ｂ１：前記（ｋｘ，ｋｙ）における前記第１の励起パルスの値
｜ｇ｜：前記（ｋｘ，ｋｙ）における前記第１の傾斜磁場と前記第２の傾斜磁場との合成傾斜磁場の値の絶対値
Ｓ（Ｋ）：前記（ｋｘ，ｋｙ）におけるサンプリング密度関数
Ｂ１´：前記（ｋｘ´，ｋｙ´）における前記第２の励起パルスの値
｜ｇ´｜：前記（ｋｘ´，ｋｙ´）における前記第３の傾斜磁場と前記第４の傾斜磁場との合成傾斜磁場の値の絶対値
Ｓ´（Ｋ´）：前記（ｋｘ´，ｋｙ´）におけるサンプリング密度関数
γ：磁気回転比
θ：前記（ｋｘ，ｋｙ）における前記第１の励起パルスの位相
θ´：前記（ｋｘ´，ｋｙ´）における前記第２の励起パルスの位相
【発明の効果】
【０００８】
｜Ｗ´（Ｋ´）｜ Ｓ´（Ｋ´）＝｜Ｗ（Ｋ）｜ Ｓ（Ｋ）およびθ´＝−θを満たす第１のシーケンスおよび第２のシーケンスを実行し、第１のシーケンスにより得られた第１のデータと、第２のシーケンスにより得られた第２のデータとを合成することにより、サイドローブの影響を低減することができる。
【図面の簡単な説明】
【０００９】
【図１】本発明の一形態の磁気共鳴装置の概略図である。
【図２】プロセッサ９が実行する処理を示す図である。
【図３】肝臓のエッジ位置を検出するときに使用される通常のナビゲータシーケンスの一例であるペンシルビーム型ナビゲータシーケンスＮＳ１を示す図である。
【図４】ナビゲータ領域Ｒｎａｖを概略的に示す図である。
【図５】ナビゲータシーケンスＮＳ１により得られたエコーＥ１に基づいて肝臓のエッジＱの位置Ｐを検出するときの説明図である。
【図６】シミュレーション結果の説明図である。
【図７】本形態で使用されるナビゲータシーケンスを示す図である。
【図８】ナビゲータシーケンスＮＳ１およびＮＳ２の説明図である。
【図９】ナビゲータシーケンスＮＳ１のｋ空間軌跡Ｊ１と、ナビゲータシーケンスＮＳ２のｋ空間軌跡Ｊ２とを示す図である。
【図１０】シミュレーション結果の説明図である。
【図１１】本形態において肝臓のエッジの位置を検出するときのフローの一例を示す図である。
【図１２】ナビゲータシーケンスＮＳ１およびＮＳ２を実行するときの説明図である。
【図１３】加算ナビゲータデータＳ１２を示す図である。
【図１４】信号強度プロファイルＩＳを概略的に示す図である。
【図１５】本形態で実行されるスキャンを示す図である。
【図１６】図１５に示すスキャンを実行するときのフローを示す図である。
【図１７】ローカライザスキャンＬＸにより取得された画像ＬＤを概略的に示す図である。
【図１８】撮影部位を概略的に示す図である。
【図１９】プレスキャンＡで実行されるシーケンスの説明図である。
【図２０】ナビゲータシーケンスのセットＮ１を実行したときの肝臓のエッジの位置Ｐ１を示す図である。
【図２１】ナビゲータシーケンスＮ１〜ＮＳｆの各々を実行したときの肝臓のエッジの位置の時間変化を概略的に示す図である。
【図２２】ウインドウＷの一例を示す図である。
【図２３】本スキャンＢの説明図である。
【図２４】ナビゲータシーケンスを用いた呼吸同期法により撮影された画像を示す図である。
【発明を実施するための形態】
【００１０】
以下、発明を実施するための形態について説明するが、本発明は、以下の形態に限定されることはない。
【００１１】
図１は、本発明の一形態の磁気共鳴装置の概略図である。
磁気共鳴装置（以下、「ＭＲ装置」と呼ぶ。ＭＲ：Magnetic Resonance）１００は、マグネット２、テーブル３、受信ＲＦコイル（以下、「受信コイル」と呼ぶ）４などを有している。
【００１２】
マグネット２は、被検体１３が収容されるボア２１を有している。また、マグネット２は、静磁場を発生させるための超伝導コイル、傾斜磁場を印加するための傾斜磁場コイル、およびＲＦパルスを送信するためのＲＦコイルなどを有している。超伝導コイルの代わりに、永久磁石を用いてもよい。
【００１３】
テーブル３は、クレードル３ａを有している。クレードル３ａは、ボア２１内に移動できるように構成されている。クレードル３ａによって、被検体１３はボア２１に搬送される。
【００１４】
受信コイル４は、被検体１３の胴部に取り付けられている。受信コイル４は、被検体１３からの磁気共鳴信号を受信する。
【００１５】
ＭＲ装置１００は、更に、送信器５、傾斜磁場電源６、受信器７、コンピュータ８、操作部１１、および表示部１２などを有している。
【００１６】
送信器５はＲＦコイルに電流を供給し、傾斜磁場電源６は傾斜磁場コイルに電流を供給する。受信器７は、受信コイル４から受け取った信号に対して、検波などの信号処理を実行する。尚、マグネット２、送信器５、傾斜磁場電源６、および受信器７を合わせたものが、スキャン手段に相当する。
【００１７】
コンピュータ８は、表示部１２に必要な情報を伝送したり、画像を再構成するなど、ＭＲ装置１００の各種の動作を実現するように、ＭＲ装置１００の各部の動作を制御する。コンピュータ８は、プロセッサ９およびメモリ１０などを有している。
【００１８】
図２に、プロセッサ９が実行する処理を示す。メモリ１０には、プロセッサ９により実行されるプログラムが記憶されている。プロセッサ９は、メモリ１０に記憶されているプログラムを読み出し、プログラムに記述されている処理を実行する。プロセッサ９は、メモリ１０に記憶されているプログラムを読み出すことにより、合成手段９１〜判断手段９５などを構成する。
【００１９】
合成手段９１は、後述するナビゲータデータＳ１およびＳ２を合成する。
変換手段９２は、合成手段９１により得られた加算ナビゲータデータＳ１２（後述する図１３参照）をフーリエ変換する。
検出手段９３は、フーリエ変換により得られたプロファイルに基づいて、肝臓のエッジの位置を検出する。
ウインドウ設定手段９４は、後述するウインドウＷ（図２２参照）を設定する。
判断手段９５は、イメージングシーケンスを実行するか否かを判断する。
【００２０】
プロセッサ９は、合成手段９１〜判断手段９５を構成する一例であり、メモリ１０に記憶されたプログラムを実行することによりこれらの手段として機能する。
【００２１】
操作部１１は、オペレータにより操作され、種々の情報をコンピュータ８に入力する。表示部１２は種々の情報を表示する。
ＭＲ装置１００は、上記のように構成されている。
【００２２】
本形態では、肝臓のエッジの位置を検出するためのナビゲータシーケンスを実行する。以下に、本形態におけるナビゲータシーケンスについて説明する。尚、以下では、本形態の効果を明確にするために、先ず、肝臓のエッジ位置の検出するための通常のナビゲータシーケンスについて説明し、通常のナビゲータシーケンスを説明した後、本形態のナビゲータシーケンスについて説明する。
【００２３】
図３は、肝臓のエッジ位置を検出するときに使用される通常のナビゲータシーケンスの一例であるペンシルビーム型ナビゲータシーケンスＮＳ１を示す図、図４は、ペンシルビーム型ナビゲータシーケンスＮＳ１を実行するときに設定されるナビゲータ領域Ｒｎａｖを概略的に示す図である。
【００２４】
ナビゲータシーケンスＮＳ１は、肝臓のエッジＱのｚ軸方向の位置Ｐを検出するために実行されるシーケンスであり、励起パルスα、傾斜磁場Ｇｘ、傾斜磁場Ｇｙ、および読出し傾斜磁場Ｇｚを有している。
【００２５】
励起パルスαはナビゲータ領域Ｒｎａｖを励起するためのＲＦパルスである。励起パルスαは、時点ｔ０において印加が開始され、時点ｔｚにおいて印加が終了するように設定されている。したがって、時点ｔ０から時点ｔｚまでが、励起パルスαの印加時間Ｔとして設定されている。
【００２６】
傾斜磁場ＧｘおよびＧｙは、ナビゲータ領域Ｒｎａｖを励起するときのｋ空間（以下、「励起ｋ空間」と呼ぶ）におけるｋ空間軌跡を決めるための傾斜磁場である。傾斜磁場ＧｘおよびＧｙは、負の極性と正の極性が交互に現れる波形を有している。傾斜磁場ＧｘおよびＧｙは、励起パルスαと同様に、時点ｔ０において印加が開始され、時点ｔｚにおいて印加が終了するように設定されている。
【００２７】
図３の右側には、励起ｋ空間内において傾斜磁場ＧｘおよびＧｙにより設定されるｋ空間軌跡Ｊ１が概略的に示されている。傾斜磁場ＧｘおよびＧｙは、スパイラル状のｋ空間軌跡Ｊ１を描きながらナビゲータ領域Ｒｎａｖの励起が行われるように設計されている。
【００２８】
傾斜磁場ＧｘおよびＧｙが印加された後、ナビゲータ領域ＲｎａｖからエコーＥ１を読み出すための読出し傾斜磁場Ｇｚが印加される。
ナビゲータシーケンスＮＳ１は、上記のように構成されている。
【００２９】
図５は、ナビゲータシーケンスＮＳ１により得られたエコーＥ１に基づいて肝臓のエッジＱの位置Ｐを検出するときの説明図である。
エコーＥ１は受信コイル４で受信され、受信器７に送信される。受信器７は、受信コイル４から受け取った信号に対して検波などの信号処理を行う。これにより、肝臓のエッジＱの位置情報を含むナビゲータデータＳ１が得られる。そして、ナビゲータデータＳ１をフーリエ変換することにより、ｚ軸方向の位置と信号強度との関係を表す信号強度プロファイルＩが得られる。信号強度プロファイルＩの肺側の信号強度は肝臓側の信号強度に比べて小さいので、肺と肝臓との境界部分には信号強度が急減に変化するする位置Ｐが現れる。この位置Ｐが、肝臓のエッジＱのｚ方向の位置Ｐとして検出される。
【００３０】
しかし、ナビゲータシーケンスＮＳ１を用いて肝臓のエッジ位置を検出した場合、エッジ位置の検出誤差が大きくなることがある。この理由について説明するため、ナビゲータシーケンスＮＳ１の励起直後の時点ｔｚにおける横磁化Ｍｘｙを求めるためのシミュレーションを行った。以下に、シミュレーション結果について説明する。
【００３１】
図６は、シミュレーション結果の説明図である。
図６には、ナビゲータシーケンスＮＳ１の励起直後の時点ｔｚの横磁化を表すプロファイルが示されている。図６（ａ）は、ｙ軸方向の位置と横磁化のＭｙ成分との関係を表すプロファイルＦ１１を示しており、図６（ｂ）は、ｙ軸方向の位置と横磁化のＭｘ成分との関係を表すプロファイルＦ１２を示している。プロファイルＦ１１およびＦ１２では、ナビゲータ領域Ｒｎａｖのｙ軸方向の位置を、ｙ＝０で示してある。
【００３２】
ナビゲータシーケンスＮＳ１では、図６（ａ）に示すように、励起の中心位置（ｙ＝０）に、大きい横磁化（Ｍｙ）を持つメインローブが現れる。しかし、図６（ｂ）に示すように、メインローブから離れた位置に、小さい横磁化（Ｍｘ）を持つサイドローブが現れる。したがって、メインローブはナビゲータ領域Ｒｎａｖ（図４参照）に現れるが、サイドローブはナビゲータ領域Ｒｎａｖから離れた位置に現れるので、ナビゲータ領域Ｒｎａｖだけでなく、ナビゲータ領域Ｒｎａｖから離れた部位でも励起が行われる。このため、不要な信号がエコーＥ１に混入してしまい、肝臓のエッジの検出誤差が大きくなるという問題がある。
【００３３】
そこで、本願発明者は、鋭意研究を重ね、肝臓のエッジの検出精度を向上させることができるシーケンスを考え出した。以下に、本形態におけるシーケンスについて説明する。
【００３４】
図７は、本形態で使用されるナビゲータシーケンスを示す図である。
本形態では、２つのナビゲータシーケンスＮＳ１およびＮＳ２が実行される。ナビゲータシーケンスＮＳ１は、励起パルスα、傾斜磁場Ｇｘ、傾斜磁場Ｇｙ、および読出し傾斜磁場Ｇｚを有している。一方、ナビゲータシーケンスＮＳ２は、励起パルスα´、傾斜磁場Ｇｘ´、傾斜磁場Ｇｙ´、および読出し傾斜磁場Ｇｚ´を有している。
【００３５】
以下、ナビゲータシーケンスＮＳ１とナビゲータシーケンスＮＳ２について具体的に説明する。
【００３６】
図８は、ナビゲータシーケンスＮＳ１およびＮＳ２の説明図である。
先ず、ナビゲータシーケンスＮＳ１について説明する。
図８に示すナビゲータシーケンスＮＳ１は、図３に示すナビゲータシーケンスＮＳ１と同じシーケンスである。したがって、図３に示すナビゲータシーケンスＮＳ１を実行することにより、エコーＥ１を収集することができる。
【００３７】
また、図８には、ナビゲータシーケンスＮＳ１の中に、励起パルスαの印加開始時点ｔ０から時間Δｔ離れた時点ｔが示されている。時点ｔにおいて、励起パルスαの値は「Ｂ１」である。Ｂ１は、以下の式で規定される。
Ｂ１＝ｂ・ｅｉθ  ・・・（１）
ここで、ｂ：時点ｔにおける励起パルスαの大きさ
θ：時点ｔにおける励起パルスαの位相
【００３８】
また、時点ｔにおいて、傾斜磁場Ｇｘの値は「ｇｘ」であり、傾斜磁場Ｇｙの値は「ｇｙ」である。
【００３９】
次に、ナビゲータシーケンスＮＳ２について説明する。
ナビゲータシーケンスＮＳ２は、励起パルスα´、傾斜磁場Ｇｘ´、傾斜磁場Ｇｙ´、および読出し傾斜磁場Ｇｚ´を有している。以下に、励起パルスα´、傾斜磁場Ｇｘ´、傾斜磁場Ｇｙ´、および読出し傾斜磁場Ｇｚ´について順に説明する。
【００４０】
（１）励起パルスα´について
また、励起パルスα´の印加開始時点ｔ０から時間Δｔ離れた時点ｔにおける励起パルスα´の値は「Ｂ１´」である。Ｂ１´は、以下の式で規定される。
Ｂ１´＝ｂ´・ｅｉθ´  ・・・（２）
ここで、ｂ´：時点ｔにおける励起パルスα´の大きさ
θ´：時点ｔにおける励起パルスα´の位相
【００４１】
（２）傾斜磁場Ｇｘ´について
図８では、傾斜磁場Ｇｘ´と傾斜磁場Ｇｘとの違いを明確にするために、ナビゲータシーケンスＮＳ２のｘ軸上に、傾斜磁場Ｇｘ´とＧｘの両方が示されている。傾斜磁場Ｇｘ´は実線で示してあり、傾斜磁場Ｇｘは破線で示してある。傾斜磁場Ｇｘ´は傾斜磁場Ｇｘの極性（正負）を逆にすることにより得られている。したがって、時点ｔにおける傾斜磁場Ｇｘ´の値を「ｇｘ´」で表すと、ｇｘ´とｇｘとの間には、以下の関係式が成り立つ。
ｇｘ´＝−ｇｘ
【００４２】
（３）傾斜磁場Ｇｙ´について
図８では、傾斜磁場Ｇｙ´と傾斜磁場Ｇｙとの違いを明確にするために、ナビゲータシーケンスＮＳ２のｙ軸上に、傾斜磁場Ｇｙ´とＧｙの両方が示されている。傾斜磁場Ｇｙ´は実線で示してあり、傾斜磁場Ｇｙは破線で示してある。傾斜磁場Ｇｙ´は傾斜磁場Ｇｙの極性（正負）を逆にすることにより得られている。したがって、時点ｔにおける傾斜磁場Ｇｙの値を「ｇｙ´」で表すと、ｇｙ´とｇｙとの間には、以下の関係式が成り立つ。
ｇｙ´＝−ｇｙ
【００４３】
（４）読出し傾斜磁場Ｇｚ´について
読出し傾斜磁場Ｇｚ´は、読出し傾斜磁場Ｇｚと同じであるので、説明は省略する。
【００４４】
次に、ナビゲータシーケンスＮＳ１およびＮＳ２のｋ空間軌跡について説明する（図９参照）。
図９は、ナビゲータシーケンスＮＳ１のｋ空間軌跡Ｊ１と、ナビゲータシーケンスＮＳ２のｋ空間軌跡Ｊ２とを示す図である。
尚、図９では、説明の便宜上、ナビゲータシーケンスＮＳ１およびＮＳ２のｚ軸は図示省略されている。
【００４５】
先ず、ナビゲータシーケンスＮＳ１のｋ空間軌跡Ｊ１について説明する。時点ｔ０では、ｋ空間軌跡Ｊ１の始点は、座標点（０，ｋｙ０）に設定されており、ｋ空間軌跡Ｊ１の終点は、原点Ｃ＝（０，０）に設定されている。傾斜磁場ＧｘおよびＧｙの印加が開始されると、ｋ空間軌跡Ｊ１は座標点（０，ｋｙ０）から螺旋を描きながら反時計回りに移動し、時点ｔｚにおいて、原点Ｃ＝（０，０）に到達する。ナビゲータシーケンスＮＳ１では、ｋ空間軌跡Ｊ１を辿るように励起を行う。尚、図９には、時点ｔ０から時間Δｔ離れた時点ｔにおけるｋ空間軌跡Ｊ１上の座標点を、（ｋｘ，ｋｙ）で示してある。座標点（ｋｘ，ｋｙ）では、ナビゲータシーケンスＮＳ１の励起パルスαの値はＢ１＝ｂ・ｅｉθ（式（１）参照）であり、傾斜磁場Ｇｘの値は「ｇｘ」であり、傾斜磁場Ｇｙの値は「ｇｙ」であることがわかる。
【００４６】
ここで、座標点（ｋｘ，ｋｙ）における傾斜磁場Ｇｘと傾斜磁場Ｇｙとの合成傾斜磁場Ｇについて考える。座標点（ｋｘ，ｋｙ）における合成傾斜磁場Ｇの値を「ｇ」とすると、ｇの絶対値｜ｇ｜は、座標点（ｋｘ，ｋｙ）における傾斜磁場Ｇｘの値ｇｘと傾斜磁場Ｇｙの値ｇｙを用いて、以下の式で表される。
【数３】
2016049298.tif 000005
【００４７】
次に、ナビゲータシーケンスＮＳ２のｋ空間軌跡Ｊ２について説明する。本形態では、上記のように、ナビゲータシーケンスＮＳ２の傾斜磁場Ｇｘ´およびＧｙ´の極性は、ナビゲータシーケンスＮＳ１の傾斜磁場ＧｘおよびＧｙの極性とは反対の極性を有するように設定されている（図８参照）。したがって、ナビゲータシーケンスＮＳ２のｋ空間軌跡Ｊ２は、ｋ空間軌跡Ｊ１に対して、原点Ｃを中心として１８０°ずれている。このため、ｋ空間軌跡Ｊ２の始点は、座標点（０，−ｋｙ０）に設定されており、ｋ空間軌跡Ｊ２の終点は、原点Ｃ＝（０，０）に設定されている。傾斜磁場Ｇｘ´およびＧｙ´の印加が開始されると、ｋ空間軌跡Ｊ２は座標点（０，−ｋｙ０）から螺旋を描きながら反時計回りに移動し、時点ｔｚにおいて、原点Ｃ＝（０，０）に到達する。ナビゲータシーケンスＮＳ２では、ｋ空間軌跡Ｊ２を辿るように励起を行う。尚、図９には、時点ｔ０から時間Δｔ離れた時点ｔにおけるｋ空間軌跡Ｊ２上の座標点を、（ｋｘ´，ｋｙ´）で示してある。
【００４８】
ここで、ｋ空間軌跡Ｊ２の座標点（ｋｘ´，ｋｙ´）と、ｋ空間軌跡Ｊ１の座標点（ｋｘ，ｋｙ）との関係について考える。ｋ空間軌跡Ｊ２は、ｋ空間軌跡Ｊ１に対して原点Ｃを中心として１８０°ずれている。したがって、（ｋｘ´，ｋｙ´）は、以下の式で表すことができる。
（ｋｘ´，ｋｙ´）＝（−ｋｘ，−ｋｙ）  ・・・（４）
【００４９】
つまり、時点ｔにおけるｋ空間軌跡Ｊ１の座標点（ｋｘ，ｋｙ）の符号を逆にすることにより、時点ｔにおけるｋ空間軌跡Ｊ２の座標点（ｋｘ´，ｋｙ´）を求めることができる。
【００５０】
座標点（ｋｘ´，ｋｙ´）では、ナビゲータシーケンスＮＳ２の励起パルスα´の値はＢ１´＝ｂ´・ｅｉθ´（式（２）参照）である。尚、本形態では、座標点（ｋｘ´，ｋｙ´）におけるｂ´およびθ´と、座標点（ｋｘ，ｋｙ）におけるｂおよびθは、以下の関係を満たすように設定されているとする。
ｂ´＝ｂ      ・・・（５ａ）
θ´＝−θ    ・・・（５ｂ）
【００５１】
式（５ａ）および（５ｂ）を式（２）に代入すると、以下の式が得られる。
Ｂ１´＝ｂ・ｅｉ（−θ）  ・・・（６）
したがって、座標点（ｋｘ´，ｋｙ´）におけるＢ１´は、式（６）で表される。
【００５２】
また、座標点（ｋｘ´，ｋｙ´）では、傾斜磁場Ｇｘの値は「ｇｘ´」であり、傾斜磁場Ｇｙの値は「ｇｙ´」である。
【００５３】
ここで、座標点（ｋｘ´，ｋｙ´）における傾斜磁場Ｇｘ´と傾斜磁場Ｇｙ´との合成傾斜磁場Ｇ´について考える。座標点（ｋｘ´，ｋｙ´）における合成傾斜磁場Ｇ´の値を「ｇ´」とすると、ｇ´の絶対値｜ｇ´｜は、座標点（ｋｘ´，ｋｙ´）における傾斜磁場Ｇｘ´の値ｇｘ´と傾斜磁場Ｇｙ´の値ｇｙ´を用いて、以下の式で表される。
【数４】
2016049298.tif 000006
【００５４】
式（３）と式（７）とを比較すると、以下の関係式が成り立つことがわかる。
｜ｇ´｜＝｜ｇ｜  ・・・（８）
【００５５】
本形態では、ナビゲータシーケンスＮＳ１だけでなく、ナビゲータシーケンスＮＳ２も実行し、２つのエコーＥ１およびＥ２を得る。２つのエコーＥ１およびＥ２を用いることにより、サイドローブを低減することができる。以下に、２つのエコーＥ１およびＥ２を用いることによりサイドローブを低減できる理由について説明する。
【００５６】
先ず、ナビゲータシーケンスＮＳ１により励起した直後の時点ｔｚにおける横磁化Ｍｘｙについて考察する。Ｍｘｙは、例えば以下の式で表すことができる（参考文献：J.Pauly, D.Nishimura, and A.Macovski, “A k-Space Analysis of Small-Tip-Angle Excitation”, J. Magn. Reson. 81, 43 (1989)）。
【数５】
2016049298.tif 000007
ここで、Ｍ０：熱平衡状態での縦磁化の大きさ
γ：磁気回転比
Ｋ：時点ｔにおける座標点（ｋｘ，ｋｙ）
Ｗ（Ｋ）：座標点Ｋ＝（ｋｘ，ｋｙ）における重み付け関数
Ｓ（Ｋ）：座標点Ｋ＝（ｋｘ，ｋｙ）におけるサンプリング密度関数
【００５７】
式（９）のＷ（Ｋ）は、座標点（ｋｘ，ｋｙ）における励起パルスαの値Ｂ１、磁気回転比γ、および座標点（ｋｘ，ｋｙ）における合成傾斜磁場Ｇの絶対値｜ｇ｜を用いて表される関数であり、以下の式で規定される。
【数６】
2016049298.tif 000008
【００５８】
式（１０）のＢ１は式（１）で表されるので、式（１）を式（１０）に代入すると、以下の式が得られる。
【数７】
2016049298.tif 000009
【００５９】
式（１１）は以下のように変形することができる。
【数８】
2016049298.tif 000010
【００６０】
式（１２）を式（９）に代入すると、以下の式が得られる。
【数９】
2016049298.tif 000011
【００６１】
次に、ナビゲータシーケンスＮＳ２により励起した直後の時点ｔｚにおける横磁化Ｍｘｙ´について考察する。Ｍｘｙ´は以下の式で表すことができる。
【数１０】
2016049298.tif 000012
ここで、Ｍ０：熱平衡状態での縦磁化の大きさ
γ：磁気回転比
Ｋ´：時点ｔにおける座標点（ｋｘ´，ｋｙ´）
Ｗ´（Ｋ）：座標点Ｋ´＝（ｋｘ´，ｋｙ´）における重み付け関数
Ｓ´（Ｋ´）：座標点Ｋ´＝（ｋｘ´，ｋｙ´）におけるサンプリング密度関数
【００６２】
式（１５）のＷ´（Ｋ´）は、座標点（ｋｘ´，ｋｙ´）における励起パルスα´の値Ｂ１´、磁気回転比γ、および座標点（ｋｘ´，ｋｙ´）における合成傾斜磁場Ｇ´の絶対値｜ｇ´｜を用いて表される関数であり、以下の式で規定される。
【数１１】
2016049298.tif 000013
【００６３】
式（１６）のＢ１´は式（２）で表されるので、式（２）を式（１６）に代入すると、以下の式が得られる。
【数１２】
2016049298.tif 000014
【００６４】
式（１７）は以下のように変形することができる。
【数１３】
2016049298.tif 000015
【００６５】
式（１８）を式（１５）に代入すると、以下の式が得られる。
【数１４】
2016049298.tif 000016
【００６６】
また、式（２０）の｜Ｗ´（Ｋ´）｜は、式（１９）で表すことができる。ここで、式（１９）の右辺のｂ´および｜ｇ´｜は、それぞれｂ´＝ｂ（式（５ａ）参照）、｜ｇ´｜＝｜ｇ｜（式（８）参照）に設定されているので、｜Ｗ´（Ｋ´）｜は、以下のように表すことができる。
【数１５】
2016049298.tif 000017
【００６７】
式（２１）の右辺は、式（１３）の右辺と同じである。したがって、式（１３）を式（２１）に代入すると、以下の式が得られる。
｜Ｗ´（Ｋ´）｜＝｜Ｗ（Ｋ）｜  ・・・（２２）
【００６８】
ここで、Ｋは、Ｋ＝（ｋｘ，ｋｙ）である。一方、Ｋ´は、式（４）から、Ｋ´＝（ｋｘ´，ｋｙ´）＝（−ｋｘ，−ｋｙ）で表される。したがって、Ｋは、Ｋ´を用いて以下の式で表すことができる。
Ｋ＝−Ｋ´  ・・・（２３）
【００６９】
式（２３）を式（２２）に代入すると、以下の式が得られる。
｜Ｗ´（Ｋ´）｜＝｜Ｗ（−Ｋ´）｜  ・・・（２４）
【００７０】
また、本形態では、サンプリング密度関数Ｓ´（Ｋ´）（式（１５）参照）は、サンプリング密度関数Ｓ（Ｋ）（式（９）参照）と同じ値に設定されている。したがって、Ｓ´（Ｋ´）とＳ（Ｋ）との間には、以下の関係式が成り立つ。
Ｓ´（Ｋ´）＝Ｓ（Ｋ）  ・・・（２５）
【００７１】
式（２３）を式（２５）に代入すると、以下の式が得られる。
Ｓ´（Ｋ´）＝Ｓ（−Ｋ´）  ・・・（２６）
【００７２】
式（２４）および式（２６）から、以下の関係式が得られる。
｜Ｗ´（Ｋ´）｜Ｓ´（Ｋ´）＝｜Ｗ（−Ｋ´）｜Ｓ（−Ｋ´）  ・・・（２７）
【００７３】
式（２７）を式（２０）に代入すると、以下の式が得られる。
【数１６】
2016049298.tif 000018
【００７４】
ここで、式（２３）および式（５ｂ）を式（２８）に代入し、変数変換すると、以下の式が得られる。
【数１７】
2016049298.tif 000019
【００７５】
式（２９）を整理すると、以下の式が得られる。
【数１８】
2016049298.tif 000020
【００７６】
式（１４）に含まれている二重積分と式（３０）に含まれている二重積分とを比較すると、複素共役の関係であることがわかる。したがって、式（１４）と式（３０）とを加算すると、以下の式が得られる。
【数１９】
2016049298.tif 000021
ここで、Ｒｅ｛｝：二重積分の実部
【００７７】
本形態では、ナビゲータシーケンスＮＳ１およびＮＳ２は、式（２７）を満たすように設定されている。また、ナビゲータシーケンスＮＳ１の励起パルスαの位相θと、ナビゲータシーケンスＮＳ２の励起パルスα´の位相θ´は、式（５ｂ）を満たすように設定されている。式（２７）および式（５ｂ）が成り立つ場合、ＭｘｙとＭｘｙ´とを加算することにより、式（３１）を得ることができる。式（３１）では、二重積分の実部が残り、虚部が除去されている。したがって、サイドローブを低減できることがわかる。
【００７８】
ここで、横磁化（Ｍｘｙ＋Ｍｘｙ´）によりサイドローブを低減できることを検証するため、横磁化Ｍｘｙ、横磁化Ｍｘｙ´、および横磁化（Ｍｘｙ＋Ｍｘｙ´）それぞれの値を求めるためのシミュレーションを行った。以下に、シミュレーション結果について説明する。
【００７９】
図１０はシミュレーション結果の説明図である。
図１０（ａ）は、ナビゲータシーケンスＮＳ１の横磁化Ｍｘｙのシミュレーション結果を示している。図１０（ａ１）は、ｙ軸方向の位置と横磁化のＭｙ成分との関係を表すプロファイルＦ１１を示しており、図１０（ａ２）は、ｙ軸方向の位置と横磁化のＭｘ成分との関係を表すプロファイルＦ１２を示している。
【００８０】
図１０（ｂ）は、ナビゲータシーケンスＮＳ２の横磁化Ｍｘｙ´のシミュレーション結果を示している。図１０（ｂ１）は、ｙ軸方向の位置と横磁化のＭｙ成分との関係を表すプロファイルＦ２１を示しており、図１０（ｂ２）は、ｙ軸方向の位置と横磁化のＭｘ成分との関係を表すプロファイルＦ２２を示している。
【００８１】
図１０（ａ）および図１０（ｂ）を参照すると、どちらのナビゲータシーケンスであっても、励起直後の横磁化には、サイドローブが現れることがわかる（図１０（ａ２）および（ｂ２）参照）。
【００８２】
これに対し、図１０（ｃ）は、横磁化（Ｍｘｙ＋Ｍｘｙ´）のシミュレーション結果を示している。図１０（ｃ１）は、ｙ軸方向の位置と横磁化のＭｙ成分との関係を表すプロファイルＦ３１を示しており、図１０（ｃ２）は、ｙ軸方向の位置と横磁化のＭｘ成分との関係を表すプロファイルＦ３２を示している。
【００８３】
図１０（ｃ２）のプロファイルＦ３２を参照すると、サイドローブの横磁化は十分に小さくなっていることがわかる。したがって、横磁化（Ｍｘｙ＋Ｍｘｙ´）により、サイドローブの影響を低減できることがわかる。
【００８４】
次に、図７のナビゲータシーケンスＮＳ１およびＮＳ２を用いて肝臓のエッジの位置を検出するときのフローの一例について説明する。
図１１は、本形態において肝臓のエッジの位置を検出するときのフローの一例を示す図である。
ステップＳＴ１では、ナビゲータシーケンスＮＳ１およびＮＳ２を順に実行する。図１２は、ナビゲータシーケンスＮＳ１およびＮＳ２を実行するときの説明図である。ナビゲータシーケンスＮＳ１を実行することにより、エコーＥ１を得ることができる。このエコーＥ１は受信コイル４で受信され、受信器７に送信される。受信器７は、受信コイル４から受け取った信号に対して検波などの信号処理を行う。これにより、肝臓のエッジの位置情報を含むナビゲータデータＳ１が得られる。
【００８５】
一方、ナビゲータシーケンスＮＳ２を実行することにより、エコーＥ２を得ることができる。このエコーＥ１は受信コイル４で受信され、受信器７に送信される。受信器７は、受信コイル４から受け取った信号に対して検波などの信号処理を行う。これにより、肝臓のエッジの位置情報を含むナビゲータデータＳ２が得られる。
ナビゲータシーケンスＮＳ１およびＮＳ２を実行した後、ステップＳＴ２に進む。
【００８６】
ステップＳＴ２では、合成手段９１（図２参照）が、ナビゲータデータＳ１とナビゲータデータＳ２とを合成する。本形態では、合成手段９１は、ナビゲータデータＳ１とナビゲータデータＳ２とを加算することによって、ナビゲータデータＳ１とナビゲータデータＳ２とを合成する。図１３に、ナビゲータデータＳ１とナビゲータデータＳ２とを加算することにより得られた加算ナビゲータデータＳ１２を示す。
【００８７】
先に説明したように、ナビゲータシーケンスＮＳ１およびＮＳ２は、式（２７）および式（５ｂ）を満たすように設定されている。式（２７）および式（５ｂ）が成り立つ場合、サイドローブを低減することができる（式（３１）参照）。したがって、ナビゲータデータＳ１とＳ２とを加算することにより、サイドローブの影響が低減された加算ナビゲータデータＳ１２を得ることができる。加算ナビゲータデータＳ１２を得た後、ステップＳＴ３に進む。
【００８８】
ステップＳＴ３では、変換手段９２（図２参照）が、加算ナビゲータデータＳ１２をフーリエ変換する。図１４に、加算ナビゲータデータＳ１２をフーリエ変換することにより得られた信号強度プロファイルＩＳを概略的に示す。次に、検出手段９３（図２参照）は、信号強度プロファイルＩＳの中から信号強度が急減に変化する位置を、肝臓のエッジＱの位置Ｐとして検出する。このようにして、図１１のフローを終了する。
【００８９】
本形態では、ナビゲータシーケンスＮＳ１およびＮＳ２は、式（２７）および式（５ｂ）を満たすように設定されている、式（２７）および式（５ｂ）が成り立つ場合、サイドローブを低減することができるので（式（３１）参照）、肝臓のエッジの位置の検出精度を高めることができる。
【００９０】
尚、本形態では、ナビゲータデータＳ１およびＳ２を加算することにより、ナビゲータデータＳ１とナビゲータデータＳ２とを合成しているが、ナビゲータデータの合成は、ナビゲータデータＳ１とナビゲータデータＳ２との加算に限定されることはない。例えば、ナビゲータデータＳ１およびＳ２を重み付け係数で重み付けし、重み付けされたナビゲータデータＳ１およびＳ２を加算してもよい。
【００９１】
本形態では、ナビゲータシーケンスＮＳ１およびＮＳ２は、｜Ｗ´（Ｋ´）｜＝｜Ｗ（Ｋ）｜およびＳ´（Ｋ´）＝Ｓ（Ｋ）を満たしている（式（２２）および式（２５）参照）。しかし、式（２７）が成り立つのであれば、｜Ｗ´（Ｋ´）｜≠｜Ｗ（Ｋ）｜でもよいし、Ｓ´（Ｋ´）≠Ｓ（Ｋ）であってもよい。
【００９２】
次に、ナビゲータシーケンスＮＳ１およびＮＳ２を用いた呼吸同期法で被検体をスキャンするときの一例について説明する。
図１５は本形態で実行されるスキャンを示す図である。
本形態では、ローカライザスキャンＬＸ、プレスキャンＡ、および本スキャンＢが実行される。
【００９３】
ローカライザスキャンＬＸは、撮影部位などを設定するために使用される画像を取得するためのスキャンである。
プレスキャンＡは、後述するウインドウＷ（図２２参照）求めるためのデータを取得するためのスキャンである。
本スキャンＢは、肝臓の画像を取得するためのスキャンである。
以下に、ローカライザスキャンＬＸ、プレスキャンＡ、および本スキャンＢを実行するときのフローについて説明する。
【００９４】
図１６は、図１５に示すスキャンを実行するときのフローを示す図である。
ステップＳＴ１では、ローカライザスキャンＬＸを実行する。図１７に、ローカライザスキャンＬＸにより取得された画像ＬＤを概略的に示す。ローカライザスキャンＬＸを実行した後、ステップＳＴ２に進む。
【００９５】
ステップＳＴ２では、画像ＬＤに基づいて撮影部位Ｒが設定される。図１８に、撮影部位を概略的に示す。オペレータは、画像ＬＤを参考にして、撮影部位Ｒと、ナビゲータ領域Ｒｎａｖとを設定する。撮影部位Ｒおよびナビゲータ領域Ｒｎａｖを設定した後、ステップＳＴ３に進む。
【００９６】
ステップＳＴ３では、プレスキャンＡが実行される（図１９参照）。
図１９は、プレスキャンＡで実行されるシーケンスの説明図である。
プレスキャンＡでは、２つのナビゲータシーケンスＮＳ１およびＮＳ２のセットが繰り返し実行される。図１９では、ナビゲータシーケンスのセットを符号「Ｎ１」、「Ｎ２」、・・・「Ｎｆ」で示してある。以下、ナビゲータシーケンスのセットＮ１〜Ｎｆについて説明する。
【００９７】
プレスキャンＡでは、先ず、ナビゲータシーケンスのセットＮ１が実行される。ナビゲータシーケンスのセットＮ１は、２つのナビゲータシーケンスＮＳ１およびＮＳ２を有しているので、ナビゲータシーケンスＮＳ１およびＮＳ２が順に実行される。ナビゲータシーケンスＮＳ１およびＮＳ２から収集したエコーに基づいて、肝臓のエッジの位置を検出することができる。肝臓のエッジの位置を検出するときの手順は、図１１に示すフローと同じであるので、説明は省略する。したがって、ナビゲータシーケンスのセットＮ１を実行したときの肝臓のエッジの位置を検出することができる。図２０に、ナビゲータシーケンスのセットＮ１を実行したときの肝臓のエッジの位置を符号「Ｐ１」で概略的に示してある。
【００９８】
以下同様に、残りのナビゲータシーケンスのセットＮ２〜Ｎｆの各々を実行するときも、図１１に示すフローが実行される。したがって、ナビゲータシーケンスのセットＮ１〜Ｎｆの各々を実行したときの肝臓のエッジの位置を検出することができる。図２１に、ナビゲータシーケンスのセットＮ１〜Ｎｆの各々を実行したときの肝臓のエッジの位置の時間変化を概略的に示す。
【００９９】
したがって、プレスキャンＡを実行している間の被検体の呼吸信号Ｓｒｅｓを得ることができる。
プレスキャンＡを実行した後、ステップＳＴ４（図１６参照）に進む。
【０１００】
ステップＳＴ４では、ウインドウ設定手段９４（図３参照）が、呼吸信号Ｓｒｅｓに基づいてウインドウＷを設定する。図２２に、ウインドウＷの一例を示す。ウインドウＷは、本スキャンＢにおいてイメージングシーケンスＤＡＱ１〜ＤＡＱｚ（図２３参照）を実行するときの基準となる肝臓のエッジの範囲を表している。以下に、ウインドウＷを設定する方法の一例について簡単に説明する。
【０１０１】
ウインドウ設定手段９４は、呼吸信号Ｓｒｅｓから、被検体が息を吐き終わったときの肝臓のエッジの位置Ｐｂを特定する。被検体が息を吐いている間、エッジの位置は増加するが、被検体が息を吸い始めると、エッジの位置は減少し始める。したがって、エッジの位置が増加してから減少に転じるときの極大値を検出することにより、被検体が息を吐き終わったときのエッジの位置Ｐｂを特定することができる。位置Ｐｂを特定した後、ウインドウ設定手段９４は、この位置Ｐｂに対して±ｘｍｍ（例えば、ｘ＝２０）の範囲を、ウインドウＷとして設定する。ウインドウＷを設定した後、ステップＳＴ５に進む。
【０１０２】
ステップＳＴ５では、本スキャンＢが実行される。
図２３は、本スキャンＢの説明図である。
本スキャンＢでは、先ず、ナビゲータシーケンスＮＳ１およびＮＳ２のセットＮｆ＋１〜Ｎｐが実行される。ナビゲータシーケンスのセットＮｆ＋１〜Ｎｐの各々が実行されるたびに、図１１に示すフローに従って肝臓のエッジの位置が検出される。図２３では、ナビゲータシーケンスのセットＮｆ＋１〜Ｎｐを実行することにより検出された肝臓のエッジの位置が、符号「Ｐｆ＋１」、「Ｐｆ＋２」、・・・「Ｐｐ」で示されている。
【０１０３】
判断手段９５（図２参照）は、ナビゲータシーケンスのセットＮｆ＋１〜Ｎｐを実行するたびに、検出された肝臓のエッジの位置がウインドウＷに入っているか否かを判断する。そして、肝臓のエッジの位置がウインドウＷの外側からウインドウＷの内側に入り込んだときに、イメージングシーケンスＤＡＱ１を実行する。
【０１０４】
ナビゲータシーケンスのセットＮｆ＋１〜Ｎｐ−１では、肝臓のエッジの位置はウインドウＷの外側であるが、ナビゲータシーケンスのセットＮｐを実行したときに肝臓のエッジの位置はウインドウＷの内側に入り込んでいる。したがって、ナビゲータシーケンスのセットＮｐを実行した直後に、イメージングシーケンスＤＡＱ１が実行される。
【０１０５】
イメージングシーケンスＤＡＱ１を実行した後、ナビゲータシーケンスのセットＮｐ＋１〜Ｎｑを実行し、肝臓のエッジの位置を検出する。そして、肝臓のエッジの位置がウインドウＷの外側からウインドウＷの内側に入り込んだときに、次のイメージングシーケンスＤＡＱ２を実行する。ナビゲータシーケンスのセットＮｐ＋１〜Ｎｑ−１では、肝臓のエッジの位置はウインドウＷの外側であるが、ナビゲータシーケンスのセットＮｑを実行したときに肝臓のエッジの位置はウインドウＷの内側に入り込んでいる。したがって、ナビゲータシーケンスのセットＮｑを実行した直後に、イメージングシーケンスＤＡＱ２が実行される。
【０１０６】
以下同様にナビゲータシーケンスのセットを実行し、肝臓のエッジの位置がウインドウＷに入ったときに、イメージングシーケンスを実行する。そして、最後のイメージングシーケンスＤＡＱｚを実行したら本スキャンＢを終了し、フローを終了する。
【０１０７】
本形態では、位相プロファイルの位相情報に基づいて、強度プロファイルに現れた不要信号の強度を除去することができるので、肝臓のエッジの位置の検出精度を高めることができる。したがって、信頼性の高い呼吸信号を得ることができるので、体動アーチファクトが低減された高品質な画像を得ることができる。図２４に、ナビゲータシーケンスを用いた呼吸同期法により撮影された画像を示す。図２４（ａ）は、ナビゲータシーケンスＮＳ１のみを用いた呼吸同期法により撮影された画像ａを示し、図２４（ｂ）は、ナビゲータシーケンスＮＳ１およびＮＳ２の両方を用いた呼吸同期法により撮影された画像ｂを示す。図２４（ａ）および（ｂ）を比較すると、画像ｂは、画像ａよりも体動アーチファクトが低減されていることがわかる。したがって、ナビゲータシーケンスＮＳ１およびＮＳ２の両方を用いた呼吸同期法が有効であることがわかる。
【０１０８】
尚、本形態では、励起パルスα´の励起時間Ｔ´は、励起パルスαの励起時間Ｔと同じ値に設定されている。しかし、ｋ空間軌跡Ｊ１上の座標点（ｋｘ，ｋｙ）と、ｋ空間軌跡Ｊ２上の座標点（ｋｘ´，ｋｙ´）において、式（２７）および式（５ｂ）の関係を満たしているのであれば、Ｔ´≠Ｔに設定してもよい。
【０１０９】
本形態では、信号強度プロファイルから肝臓のエッジの位置を検出することにより、被検体の呼吸信号を得ているが、肝臓のエッジの位置を検出する代わりに、ＬＳＱ（Least Square）法などを用いて肝臓のエッジの移動量を求め、肝臓のエッジの移動量から呼吸信号を作成してもよい。
【０１１０】
本形態では、ナビゲータシーケンスＮＳ１およびＮＳ２を用いて呼吸信号を求めているが、本発明は、呼吸信号以外の生体信号（例えば、心拍信号）を求める場合や、イメージングデータを収集する場合にも適用することができる。例えば、式（２７）および式（３ｂ）を満たすように構成された２つのイメージングシーケンスを実行することにより、サイドローブが低減されたイメージングデータを収集することができるので、高品質な画像を得ることができる。
【０１１１】
尚、本形態では、動きのある部位をスキャンする例について説明しているが、本発明は、動きのある部位のスキャンに限定されることはなく、静止部位をスキャンする場合にも適用することができる。
【符号の説明】
【０１１２】
２  マグネット
３  テーブル
３ａ  クレードル
４  受信コイル
５  送信器
６  傾斜磁場電源
７  受信器
８  コンピュータ
９  プロセッサ
１０  メモリ
１１  操作部
１２  表示部
１３  被検体
２１  ボア
９１  合成手段
９２  変換手段
９３  検出手段
９４  ウインドウ設定手段
９５  判断手段
１００  ＭＲ装置
【図１】
2016049298.tif 000025
【図２】
2016049298.tif 000026
【図３】
2016049298.tif 000027
【図４】
2016049298.tif 000028
【図５】
2016049298.tif 000029
【図６】
2016049298.tif 000030
【図７】
2016049298.tif 000031
【図８】
2016049298.tif 000032
【図９】
2016049298.tif 000033
【図１０】
2016049298.tif 000034
【図１１】
2016049298.tif 000035
【図１２】
2016049298.tif 000036
【図１３】
2016049298.tif 000037
【図１４】
2016049298.tif 000038
【図１５】
2016049298.tif 000039
【図１６】
2016049298.tif 000040
【図１７】
2016049298.tif 000041
【図１８】
2016049298.tif 000042
【図１９】
2016049298.tif 000043
【図２０】
2016049298.tif 000044
【図２１】
2016049298.tif 000045
【図２２】
2016049298.tif 000046
【図２３】
2016049298.tif 000047
【図２４】
2016049298.tif 000048