JP2014023703A - 磁気共鳴装置およびプログラム - Google Patents

Abstract

【課題】リンギングやブラーなどのアーチファクトを低減する磁気共鳴装置を提供する。
【解決手段】イメージングシーケンスＩＳ_ｎの他に、ナビゲータシーケンスＮＳ_ｎを実行する。ナビゲータシーケンスＮＳ_ｎは、肺と肝臓との境界部分を含むナビゲータ領域Ｒ_ｎａｖ内のＳＩ方向の信号強度の変化を表す信号強度データを取得するためのシーケンスである。ナビゲータシーケンスＮＳ_ｎにより得られた信号強度データＤ_１〜Ｄ_ｚごとに、肝臓のエッジの位置Ｐ_０から肝臓側（Ｉ側）に距離Δｄだけ離れた位置Ｐ_０＋Δｄの信号強度Ｎ_１〜Ｎ_ｚを検出する。そして、信号強度Ｎ_１〜Ｎ_ｚに基づいて、イメージングシーケンスにより得られた撮影領域のイメージングデータＳ_１〜Ｓ_ｚを補正する。
【選択図】図９

Description

本発明は、所定の部位を撮影する磁気共鳴装置、およびその磁気共鳴装置に適用されるプログラムに関する。

磁気共鳴装置を用いて被検体の診断を行う方法として、被検体に造影剤を投与し、造影剤を投与した後に現れる各時相の画像データを取得する方法が知られている。

特開２０１１−１６７４０８号

近年、肝臓のＭＲＩ検査ではガドキセト酸ナトリウム（Gd-EOB-DTPA：以下、EOBと呼ぶ）を用いた造影ダイナミック撮像検査が頻繁に行われるようになったが、従前より使用されていたガドリニウムキレート（Gd-DTPA）と比較して、投与量が半分である。したがって、EOBを用いた撮像では造影剤の濃度が一定になる時間が短く、リンギングやブラーなどのアーチファクトが現れやすいという問題がある。このため、リンギングやブラーなどのアーチファクトを低減することが望まれている。

本発明の第１の態様は、被検体の所定の部位を含む第１の領域からナビゲータ信号を収集するとともに、前記所定の部位を含む第２の領域から前記所定の部位の画像を再構成するときに使用されるイメージングデータを収集する磁気共鳴装置であって、
前記ナビゲータ信号と前記イメージングデータとを収集するためのスキャンを実行するスキャン手段と、
前記ナビゲータ信号に基づいて、前記イメージングデータの信号強度を補正する補正手段と、
信号強度が補正されたイメージングデータに基づいて画像再構成を行う再構成手段と、
を有する磁気共鳴装置である。

本発明の第２の態様は、被検体の所定の部位を含む第１の領域からナビゲータ信号を収集するとともに、前記所定の部位を含む第２の領域から前記所定の部位の画像を再構成するときに使用されるイメージングデータを収集する磁気共鳴装置であって、
前記ナビゲータ信号と前記イメージングデータとを収集するためのスキャンを実行するスキャン手段と、
前記ナビゲータ信号に基づいて、前記イメージングデータを収集するときに使用されるシーケンスのＲＦパルスのフリップ角を補正する補正手段と、
を有する、磁気共鳴装置である。

本発明の第３の態様は、被検体の所定の部位を含む第１の領域からナビゲータ信号を収集するとともに、前記所定の部位を含む第２の領域から前記所定の部位の画像を再構成するときに使用されるイメージングデータを収集するためのスキャンを実行する磁気共鳴装置のプログラムであって、
前記ナビゲータ信号に基づいて、前記イメージングデータの信号強度を補正する補正処理と、
信号強度が補正されたイメージングデータに基づいて画像再構成を行う再構成処理と、
を計算機に実行させるためのプログラムである。

本発明の第４の態様は、被検体の所定の部位を含む第１の領域からナビゲータ信号を収集するとともに、前記所定の部位を含む第２の領域から前記所定の部位の画像を再構成するときに使用されるイメージングデータを収集するためのスキャンを実行する磁気共鳴装置のプログラムであって、
前記ナビゲータ信号に基づいて、前記イメージングデータを収集するときに使用されるシーケンスのＲＦパルスのフリップ角を補正する補正処理、
を計算機に実行させるためのプログラムである。

信号強度に基づいてイメージングデータ又はフリップ角を補正するので、リンギングやブラーなどのアーチファクトを低減することができる。

本発明の第１の形態の磁気共鳴装置の概略図である。 第１の形態で実行されるスキャンを示す図である。 撮影領域を概略的に示す図である。 第１の形態において被検体を撮影するときのフローを示す図である。 息止めスキャンＢＨの説明図である。 息止めスキャンＢＨの間における肝臓のエッジのＳＩ方向の位置Ｐ_０を検出するときの説明図である。 位置Ｐ_０＋Δｄの信号強度Ｎ_１〜Ｎ_ｚを示す図である。 プロファイルＦを概略的に示す図である。 式（１）′を用いてイメージングデータＳ_１〜Ｓ_ｚを補正するときの説明図である。 ナビゲータシーケンスＮＳ_ｎとＮＳ_ｎ＋１との間に複数のイメージングシーケンスが実行される場合の一例を示す図である。 第２の形態のＭＲ装置の概略図である。 第２の形態で実行されるスキャンを示す図である。 撮影領域を概略的に示す図である。 第２の形態のフローを示す図である。 息止めスキャンＢＨの説明図である。 撮影領域Ｒの信号強度を算出するときの説明図である。 式（２）′を用いてイメージングデータＳ_１〜Ｓ_ｚを補正するときの説明図である。 第３の形態で実行されるスキャンの説明図である。 第３の形態において被検体を撮影するときのフローを示す図である。 プリスキャンＡの説明図である。 エッジ位置プロファイルを概略的に示す図である。 吐き終りの範囲ＡＷの一例を示す図である。 本スキャンＢの説明図である。 検出された信号強度を示す図である。 式（３）を用いてイメージングデータを補正するときの説明図である。 第４の形態で実行されるスキャンの説明図である。 第４の形態において被検体を撮影するときのフローを示す図である。 トリガレベルＴＬの一例を示す図である。 本スキャンＢの説明図である。 検出された信号強度を示す図である。 式（４）を用いてイメージングデータを補正するときの説明図である。 第５の形態のＭＲ装置の概略図である。 第５の形態で実行されるスキャンを示す図である。 ナビゲータシーケンスＮＳ_１により収集されたナビゲータ信号Ａ_１と、イメージングシーケンスＩＳ_１により収集されたイメージングデータＳ_１とを示す図である。 ナビゲータシーケンスＮＳ_２により収集されたナビゲータ信号Ａ_２を示す図である。 フリップ角αの補正方法の説明図である。 息止めスキャンＢＨにより得られたイメージングデータＳ_１〜Ｓ_ｚを示す図である。 第６の形態のＭＲ装置を示す概略図である。 第６の形態で実行されるスキャンを示す図である。 シーケンスＣ_１により収集されたナビゲータ信号Ａ_１およびイメージングデータＳ_１を示す図である。 シーケンスＣ_２により収集されたナビゲータ信号Ａ_２およびイメージングデータＳ_２を示す図である。 シーケンスＣ_３のＲＦパルスＸ_３のフリップ角αの補正方法の説明図である。 補正後の信号強度Ｎ_{３＿ｃｏｒｒ}が算出された様子を示す図である。 シーケンスＣ_４のフリップ角を補正するときの説明図である。

以下、発明を実施するための形態について説明するが、本発明は、以下の形態に限定されることはない。

（１）第１の形態
図１は、本発明の第１の形態の磁気共鳴装置の概略図である。
磁気共鳴装置（以下、「ＭＲ装置」と呼ぶ。ＭＲ：Magnetic Resonance）１００は、マグネット２、テーブル３、受信コイル４などを有している。

マグネット２は、被検体１２が収容されるボア２１と、超伝導コイル２２と、勾配コイル２３と、ＲＦコイル２４とを有している。超伝導コイル２２は静磁場を印加し、勾配コイル２３は勾配磁場を印加し、ＲＦコイル２４はＲＦパルスを送信する。尚、超伝導コイル２２の代わりに、永久磁石を用いてもよい。

テーブル３は、クレードル３ａを有している。クレードル３ａは、ボア２１内に移動できるように構成されている。クレードル３ａによって、被検体１２はボア２１に搬送される。

受信コイル４は、被検体１２の腹部に取り付けられている。受信コイル４は、被検体１２からの磁気共鳴信号を受信する。

ＭＲ装置１００は、更に、造影剤注入装置５、送信器６、勾配磁場電源７、受信器８、制御部９、操作部１０、および表示部１１などを有している。

造影剤注入装置５は、被検体に造影剤を注入する。
送信器６はＲＦコイル２４に電流を供給し、勾配磁場電源７は勾配コイル２３に電流を供給する。
受信器８は、受信コイル４から受け取った信号に対して、検波などの信号処理を実行する。

制御部９は、表示部１１に必要な情報を伝送したり、受信器８から受け取ったデータに基づいて画像を再構成するなど、ＭＲ装置１００の各種の動作を実現するように、ＭＲ装置１００の各部の動作を制御する。制御部９は、例えばコンピュータ（computer）によって構成される。制御部９は、補正手段９１および再構成手段９２などを有している。

補正手段９１は、イメージングデータの信号強度を補正する。補正手段９１は、信号強度データ作成手段９１ａ、検出手段９１ｂ、演算手段９１ｃなどを有している。

信号強度データ作成手段９１ａは、ナビゲータ信号Ａ_ｎに基づいて、ナビゲータ領域Ｒ_ｎａｖ内のＳＩ方向の信号強度の変化を表す信号強度データＤ_ｎを作成する（例えば、図５参照）。
検出手段９１ｂは、信号強度データに基づいて肝臓の位置における信号強度を検出する。
演算手段９１ｃは、後述する式を用いて、補正後のイメージングデータを求める。

再構成手段９２は、補正後のイメージングデータを用いて画像を再構成する。

制御部９は、補正手段９１および再構成手段９２を構成する一例であり、所定のプログラムを実行することにより、これらの手段として機能する。

操作部１０は、オペレータにより操作され、種々の情報を制御部９に入力する。表示部１１は種々の情報を表示する。
ＭＲ装置１００は、上記のように構成されている。

図２は第１の形態で実行されるスキャンを示す図、図３は撮影領域を概略的に示す図である。
第１の形態では、息止めスキャンＢＨが実行される。
息止めスキャンＢＨは、被検体に造影剤を注入し、被検体が息止めをした状態で実行されるスキャンである。息止めスキャンＢＨでは、ナビゲータシーケンスＮＳ_ｎ（ｎ＝１〜ｚ）と、イメージングシーケンスＩＳ_ｎ（ｎ＝１〜ｚ）とが交互に実行される。

ナビゲータシーケンスＮＳ_ｎは、肺と肝臓との境界部分を含むナビゲータ領域Ｒ_ｎａｖ内のＳＩ方向の信号強度の変化を表す信号強度データを取得するためのシーケンスである。イメージングシーケンスＩＳ_ｎは肝臓を含む撮影領域Ｒの画像を再構成するときに使用されるイメージングデータを取得するためのシーケンスである。
以下に、息止めスキャンＢＨを実行するときのフローについて説明する。

図４は、第１の形態において被検体を撮影するときのフローを示す図である。
ステップＳＴ１では、被検体に造影剤を注入し、息止めスキャンＢＨを実行する。息止めスキャンＢＨでは、図５に示すように、ナビゲータシーケンスＮＳ_ｎとイメージングシーケンスＩＳ_ｎとが交互に実行される。ナビゲータシーケンスＮＳ_ｎを実行することにより、ナビゲータ領域Ｒ_ｎａｖからナビゲータ信号Ａ_ｎ（ｎ＝１〜ｚ）が収集され、イメージングシーケンスＩＳ_ｎを実行することにより、撮影領域Ｒの画像を再構成するときに使用されるイメージングデータＳ_ｎ（ｎ＝１〜ｚ）が収集される。信号強度データ作成手段９１ａ（図１参照）は、ナビゲータ信号Ａ_ｎを一次元フーリエ変換（１Ｄ−ＦＴ）し、ナビゲータ領域Ｒ_ｎａｖ内のＳＩ方向の信号強度の変化を表す信号強度データＤ_１〜Ｄ_ｚを作成する。信号強度データＤ_１〜Ｄ_ｚを得た後、ステップＳＴ２に進む。

ステップＳＴ２では、検出手段９１ｂ（図１参照）が、息止めスキャンＢＨの間における肝臓のエッジのＳＩ方向の位置Ｐ_０を検出する。

図６は、息止めスキャンＢＨの間における肝臓のエッジのＳＩ方向の位置Ｐ_０を検出するときの説明図である。

本形態では、肺側が低信号となり肝臓側が高信号となるようなナビゲータシーケンスを用いている。したがって、肺と肝臓との境界には信号強度の段差が現れるので、この段差を検出することにより、肝臓のエッジの位置Ｐ_０を検出することができる。尚、本形態では、被検体は息を止めているので、スキャンの間、肝臓のエッジの位置は動かないと考えることができる。したがって、信号強度データＤ_１〜Ｄ_ｚのうちのいずれかの信号強度データを用いて検出された肝臓のエッジの位置を、息止めスキャンの間における肝臓のエッジの位置Ｐ_０とすることができる。ここでは、信号強度データＤ_１を用いて検出された肝臓のエッジの位置を、息止めスキャンの間における肝臓のエッジの位置Ｐ_０とする。

肝臓のエッジの位置Ｐ_０を検出した後、ステップＳＴ３に進む。
ステップＳＴ３では、検出手段９１ｂが、信号強度データＤ_１〜Ｄ_ｚごとに、肝臓のエッジの位置Ｐ_０から肝臓側（Ｉ側）に距離Δｄだけ離れた位置Ｐ_０＋Δｄの信号強度Ｎ_１〜Ｎ_ｚを検出する。図７に、位置Ｐ_０＋Δｄの信号強度Ｎ_１〜Ｎ_ｚを示す。信号強度Ｎ_１〜Ｎ_ｚを検出することにより、肝臓における信号強度の時間変化を表すプロファイルＦを得ることができる。図８に、そのプロファイルＦを概略的に示す。肝臓に流入する造影剤の濃度が増加するに伴って、信号強度Ｎ_ｎの値は増加する。したがって、プロファイルＦから、時間の経過とともに、肝臓に流入する造影剤の濃度が増加していることがわかる。尚、Δｄの値が小さすぎると、肝臓と肺との境界に現れる段差の信号強度を検出してしまう可能性が高くなるので、肺側の領域の信号強度を検出してしまう恐れがある。一方、Δｄの値が大きすぎると、腸の領域の信号強度を検出してしまう恐れがある。したがって、Δｄは、肝臓の内側の信号強度を検出することができるように設定する必要がある。一例としては、Δｄ＝５ｃｍとすることにより、肝臓とは別の部位の信号強度を検出する危険を回避することができる。
プロファイルＦを得た後、ステップＳＴ４に進む。

ステップＳＴ４では、信号強度Ｎ_１〜Ｎ_ｚに基づいて、イメージングシーケンスにより得られた撮影領域のイメージングデータＳ_１〜Ｓ_ｚを補正する。この補正は、以下の式を用いて行われる。
Ｓ_{ｎ，ｃｏｒｒ}＝Ｓ_ｎ＊Ｎ_ｉ／Ｎ_ｎ ・・・（１）

Ｓ_ｎ：ｎ番目のナビゲータ信号の直後に得られたイメージングデータ
Ｓ_{ｎ，ｃｏｒｒ}：Ｓ_ｎを補正した後のイメージングデータ
Ｎ_ｉ：信号強度Ｎ_１〜Ｎ_ｚの中から選択されたｉ番目の信号強度
Ｎ_ｎ：ｎ番目の信号強度

演算手段９１１ｃ（図１参照）は、式（１）を用いて、補正後のイメージングデータＳ_{ｎ，ｃｏｒｒ}を計算する。式（１）では、イメージングデータＳ_ｎをＮ_ｉ／Ｎ_ｎ倍している。したがって、肝臓に流入する造影剤の濃度の違いにより生じるイメージングデータＳ_１〜Ｓ_ｚの信号強度の差を補正することができる。式（１）のＮ_ｉは、信号強度Ｎ_１〜Ｎ_ｚの中から選択されたｉ番目の信号強度であるが、Ｎ_ｉは、例えば、ｉ＝１、すなわち、Ｎ_１とすることができる。ｉ＝１の場合、式（１）は以下のように書き換えられる。
Ｓ_{ｎ，ｃｏｒｒ}＝Ｓ_ｎ＊Ｎ_１／Ｎ_ｎ ・・・（１）′

以下に、式（１）′を用いてイメージングデータＳ_１〜Ｓ_ｚを補正する方法について、図９を参照しながら説明する。

図９は、式（１）′を用いてイメージングデータＳ_１〜Ｓ_ｚを補正するときの説明図である。

イメージングデータＳ_１を補正する場合、式（１）′のｎに、ｎ＝１を代入すればよい。
Ｓ_{１，ｃｏｒｒ}＝Ｓ_１＊Ｎ_１／Ｎ_１
＝Ｓ_１
したがって、補正後のイメージングデータＳ_{１，ｃｏｒｒ}は、補正前のイメージングデータＳ_１と同じである。

しかし、イメージングデータＳ_２を補正する場合、式（１）′のｎに、ｎ＝２が代入されるので、補正後のイメージングデータＳ_{２，ｃｏｒｒ}は、以下の式で表される。
Ｓ_{２，ｃｏｒｒ}＝Ｓ_２＊Ｎ_１／Ｎ_２

したがって、補正後のイメージングデータＳ_{２，ｃｏｒｒ}は、補正前のイメージングデータＳ_２のＮ_１／Ｎ_２倍になる。

以下同様に、式（１）′のｎの値を変更し、Ｓ_ｎをＮ_１／Ｎ_ｎ倍することにより、補正後のイメージングデータＳ_{ｎ，ｃｏｒｒ}を得ることができる。イメージングデータを補正した後、ステップＳＴ５に進む。

ステップＳＴ５では、補正後のイメージングデータＳ_{１，ｃｏｒｒ}〜Ｓ_{ｚ，ｃｏｒｒ}を用いて、画像再構成を行い、フローを終了する。

第１の形態では、ナビゲータシーケンスで得られた肝臓の位置における信号強度Ｎ_１〜Ｎ_ｚに基づいて、イメージングシーケンスで取得されたイメージングデータＳ_１〜Ｓ_ｚを補正している。したがって、肝臓に流入する造影剤の濃度の違いにより生じるイメージングデータＳ_１〜Ｓ_ｚの信号強度の差を補正することができるので、リンギングやブラーなどのアーチファクトを低減することができる。

尚、上記の例では、信号強度Ｎ_１〜Ｎ_ｚは、それぞれ信号強度データＤ_１〜Ｄ_ｚから検出されている。しかし、信号強度Ｎ_１〜Ｎ_ｚの中の一部の信号強度を、他の信号強度から推定してもよい。第１の形態では、信号強度Ｎ_１〜Ｎ_ｚは、時間とともに増加する傾向があるので（図８参照）、信号強度Ｎ_１とＮ_ｚとが分かれば、信号強度が時間とともに増加するプロファイルＦを推定することができる。したがって、信号強度Ｎ_１とＮ_ｚとに基づいてプロファイルＦを推定する推定手段を設けておき、推定されたプロファイルＦから、他の信号強度Ｎ_２〜Ｎ_ｚ−１を推定してもよい。この場合、ナビゲータシーケンスＮＳ_１およびＮＳ_ｚのみを実行するだけで信号強度Ｎ_１〜Ｎ_ｚを求めることができるので、息止めスキャンＢＨの時間を更に短くすることができる。

第１の形態では、ｎ番目の信号強度Ｎ_ｎを用いて、ｎ番目のナビゲータ信号の直後に得られたｎ番目のイメージングデータＳ_ｎを補正している。しかし、ｎ番目の信号強度Ｎ_ｎを用いて、ｎ番目のナビゲータ信号の直前に収集されたｎ−１番目のイメージングデータＳ_ｎ−１を補正してもよい。

尚、上記の説明では、ナビゲータシーケンスＮＳ_ｎとＮＳ_ｎ＋１との間に１つのイメージングシーケンスＩＳ_ｎが実行される場合について述べられている。しかし、ナビゲータシーケンスＮＳ_ｎとＮＳ_ｎ＋１との間に複数のイメージングシーケンスを実行してもよい（図１０参照）。

図１０は、ナビゲータシーケンスＮＳ_ｎとＮＳ_ｎ＋１との間に複数のイメージングシーケンスが実行される場合の一例を示す図である。

図１０では、ナビゲータシーケンスＮＳ_ｎとＮＳ_ｎ＋１との間に、複数のイメージングデータが収集される。このような場合は、ナビゲータシーケンスＮＳ_ｎとＮＳ_ｎ＋１との間に収集された複数のイメージングデータの各々をＮ_１／Ｎ_ｎ倍すればよい。例えば、ナビゲータシーケンスＮＳ_２とＮＳ_３との間では、イメージングデータＳ_２１〜Ｓ_２ｑが収集されるので、これらのイメージングデータＳ_２１〜Ｓ_２ｑの各々を、Ｎ_１／Ｎ_２倍すればよい。これにより、補正後のイメージングデータＳ_{２１，ｃｏｒｒ}〜Ｓ_{２ｑ，ｃｏｒｒ}を得ることができる。

尚、図１０では、ｎ番目の信号強度Ｎ_ｎを用いて、ｎ番目のナビゲータ信号とｎ＋１番目のナビゲータ信号との間に収集されたイメージングデータを補正している。しかし、ｎ番目の信号強度Ｎ_ｎを用いて、ｎ−１番目のナビゲータ信号とｎ番目のナビゲータ信号との間に収集されたイメージングデータを補正してもよい。

また、第１の形態では、息止めスキャンＢＨを終了した後でイメージングデータを補正しているが、息止めスキャンＢＨを実行しながら、イメージングデータを補正してもよい。

（２）第２の形態
第１の形態では、ナビゲータ信号とイメージングデータとを別々のシーケンスで収集していたが、第２の形態では、ナビゲータ信号とイメージングデータとを同じシーケンスで収集する場合について説明する。

図１１は、第２の形態のＭＲ装置の概略図である。
第２の形態のＭＲ装置２００は、第１の形態のＭＲ装置１００の信号強度データ作成手段９１ａおよび検出手段９１ｂの代わりに、信号強度算出手段９１ｄを有している。尚、その他の構成は、第１の形態のＭＲ装置１００と同じである。

図１２は第２の形態で実行されるスキャンを示す図、図１３は撮影領域を概略的に示す図である。
第２の形態では、息止めスキャンＢＨが実行される。第２の形態の息止めスキャンＢＨでは、シーケンスＣ_ｎ（ｎ＝１〜ｚ）が実行される。シーケンスＣ_ｎは、撮影領域Ｒの画像を再構成するときに使用されるイメージングデータＳ_ｎを取得するためのシーケンスである。また、シーケンスＣ_ｎには、勾配磁場ＧｘおよびＧｙが印加される直前に待ち時間Ｔ_ｗａｉｔが設けられており、待ち時間Ｔ_ｗａｉｔのときにナビゲータ信号Ａ_ｎが取得される。ナビゲータ信号Ａ_ｎは、撮影領域Ｒのｋ空間の中心のデータを表している。待ち時間Ｔ_ｗａｉｔは、例えば２０μｍである。

したがって、シーケンスＣ_ｎを実行するたびに、イメージングデータＳ_ｎの他に、ナビゲータ信号Ａ_ｎを得ることができる。第２の形態では、シーケンスＣ_ｎにより得られたナビゲータ信号Ａ_ｎを用いてイメージングデータＳ_ｎを補正する。以下に、第２の形態において被検体を撮影するときのフローについて説明する。

図１４は、第２の形態のフローを示す図である。
ステップＳＴ１では、被検体に造影剤を注入し、息止めスキャンＢＨを実行する。息止めスキャンＢＨでは、図１５に示すように、シーケンスＣ_ｎが実行される。シーケンスＣ_ｎを実行することにより、ナビゲータ信号Ａ_ｎ（ｎ＝１〜ｚ）とイメージングデータＳ_ｎ（ｎ＝１〜ｚ）が収集される。息止めスキャンＢＨを実行した後、ステップＳＴ２に進む。

ステップＳＴ２では、信号強度算出手段９１ｄ（図１１参照）が、ナビゲータ信号Ａ_１〜Ａ_ｚごとに撮影領域Ｒの信号強度を算出する（図１６参照）。

図１６は、撮影領域Ｒの信号強度を算出するときの説明図である。
第２の形態では、信号強度算出手段９１ｄは、撮影領域Ｒから得られたナビゲータ信号Ａ_１〜Ａ_ｚの積分値を計算する。これにより、撮影領域Ｒの信号強度Ｎ_１〜Ｎ_ｚを算出することができる。信号強度Ｎ_１〜Ｎ_ｚを算出した後、ステップＳＴ３に進む。

ステップＳＴ３では、信号強度Ｎ_１〜Ｎ_ｚに基づいて、シーケンスＣ_１〜Ｃ_ｚにより得られた撮影領域ＲのイメージングデータＳ_１〜Ｓ_ｚを補正する。この補正は、以下の式を用いて行われる。
Ｓ_{ｎ，ｃｏｒｒ}＝Ｓ_ｎ＊Ｎ_ｉ／Ｎ_ｎ ・・・（２）

Ｓ_ｎ：ｎ番目のイメージングデータ
Ｓ_{ｎ，ｃｏｒｒ}：Ｓ_ｎを補正した後のイメージングデータ
Ｎ_ｉ：信号強度Ｎ_１〜Ｎ_ｚの中から選択されたｉ番目の信号強度
Ｎ_ｎ：ｎ番目の信号強度

演算手段９１ｃ（図１１参照）は、式（２）を用いて、補正後のイメージングデータＳ_{ｎ，ｃｏｒｒ}を計算する。式（２）では、イメージングデータＳ_ｎをＮ_ｉ／Ｎ_ｎ倍している。したがって、肝臓に流入する造影剤の濃度の違いにより生じるイメージングデータＳ_１〜Ｓ_ｚの信号強度の差を補正することができる。式（２）のＮ_ｉは、信号強度Ｎ_１〜Ｎ_ｚの中から選択されたｉ番目の信号強度であるが、Ｎ_ｉは、例えば、ｉ＝１、すなわち、Ｎ_１とすることができる。ｉ＝１の場合、式（２）は以下のように書き換えられる。
Ｓ_{ｎ，ｃｏｒｒ}＝Ｓ_ｎ＊Ｎ_１／Ｎ_ｎ ・・・（２）′

以下に、式（２）′を用いてイメージングデータＳ_１〜Ｓ_ｚを補正する方法について、図１７を参照しながら説明する。

図１７は、式（２）′を用いてイメージングデータＳ_１〜Ｓ_ｚを補正するときの説明図である。
イメージングデータＳ_１を補正する場合、式（２）′のｎに、ｎ＝１を代入すればよい。
Ｓ_{１，ｃｏｒｒ}＝Ｓ_１＊Ｎ_１／Ｎ_１
＝Ｓ_１

したがって、補正後のイメージングデータＳ_{１，ｃｏｒｒ}は、補正前のイメージングデータＳ_１と同じである。

しかし、イメージングデータＳ_２を補正する場合、式（２）′のｎに、ｎ＝２が代入されるので、補正後のイメージングデータＳ_{２，ｃｏｒｒ}は、以下の式で表される。
Ｓ_{２，ｃｏｒｒ}＝Ｓ_２＊Ｎ_１／Ｎ_２

したがって、補正後のイメージングデータＳ_{２，ｃｏｒｒ}は、補正前のイメージングデータＳ_２のＮ_１／Ｎ_２倍になる。

以下同様に、式（２）′のｎの値を変更し、Ｓ_ｎをＮ_１／Ｎ_ｎ倍することにより、補正後のイメージングデータＳ_{ｎ，ｃｏｒｒ}を得ることができる。イメージングデータを補正した後、ステップＳＴ４に進む。

ステップＳＴ４では、補正後のイメージングデータＳ_{１，ｃｏｒｒ}〜Ｓ_{ｚ，ｃｏｒｒ}を用いて、画像再構成を行い、フローを終了する。

第２の形態では、シーケンスＣ_ｎで得られた肝臓の位置における信号強度Ｎ_１〜Ｎ_ｚに基づいて、イメージングデータＳ_１〜Ｓ_ｚを補正している。したがって、肝臓に流入する造影剤の濃度の違いにより生じるイメージングデータＳ_１〜Ｓ_ｚの信号強度の差を補正することができるので、リンギングやブラーなどのアーチファクトを低減することができる。

また、第２の形態では、シーケンスＣ_ｎを実行するときに励起される撮影領域Ｒから、ナビゲータ信号とイメージングデータとを取得することができる。したがって、ナビゲータ信号を取得するときに励起される領域と、イメージングデータを取得するときに励起される領域とが同じであるので、イメージングデータの補正の精度を向上させることができる。

尚、第２の形態では、ｎ番目のシーケンスＣ_ｎにより収集された信号強度Ｎ_ｎを用いて、ｎ番目のシーケンスＣ_ｎにより収集されたイメージングデータＳ_ｎを補正している。しかし、ｎ番目のシーケンスＣ_ｎにより収集された信号強度Ｎ_ｎを用いて、ｎ−１番目のシーケンスＣ_ｎ−１により収集されたイメージングデータＳ_ｎ−１を補正してもよい。

（３）第３の形態
第１および第２の形態では、被検体に息止めさせた状態でスキャンを実行する例について説明したが、第３の形態では、被検体に自由に呼吸をさせた状態でスキャンを実行する例について説明する。尚、ＭＲ装置のハードウェア構成は、第１の形態と同じである。

図１８は、第３の形態で実行されるスキャンの説明図である。
第３の形態では、プリスキャンＡと、本スキャンＢとが実行される。

プリスキャンＡは、呼吸信号を検出するためのスキャンである。プリスキャンＡでは、ナビゲータシーケンスＮＳ_１〜ＮＳ_ｘが実行される。ナビゲータシーケンスＮＳ_１〜ＮＳ_ｘは、第１の形態のナビゲータシーケンスＮＳ_ｎ（ｎ＝１〜ｚ）と同様に、ナビゲータ領域Ｒ_ｎａｖ内のＳＩ方向の信号強度の変化を表す信号強度データを取得するためのシーケンスである。

本スキャンＢは、肝臓を含む撮影領域Ｒの画像を再構成するときに使用されるイメージングデータを取得するためのスキャンである。本スキャンＢでは、第１の形態の息止めスキャンＢＨと同様に、ナビゲータシーケンスＮＳ_ｎ（ｎ＝１〜ｚ）と、イメージングシーケンスＩＳ_ｎ（ｎ＝１〜ｚ）とが交互に実行される。

以下に、プリスキャンＡおよび本スキャンＢにより被検体を撮影するときのフローについて説明する。

図１９は、第３の形態において被検体を撮影するときのフローを示す図である。
ステップＳＴ１では、プリスキャンＡが実行される（図２０参照）。

図２０は、プリスキャンＡの説明図である。
プリスキャンＡでは、ナビゲータシーケンスＮＳ_１〜ＮＳ_ｘが実行される。ナビゲータシーケンスＮＳ_１〜ＮＳ_ｘを実行することにより、ナビゲータ領域Ｒ_ｎａｖからナビゲータ信号Ａ_ｎ（ｎ＝１〜ｘ）が収集される。信号強度データ作成手段９１ａ（図１参照）は、ナビゲータ信号Ａ_ｎを一次元フーリエ変換（１Ｄ−ＦＴ）し、ナビゲータ領域Ｒ_ｎａｖ内のＳＩ方向の信号強度の変化を表す信号強度データＤ_１〜Ｄ_ｘを作成する。肺と肝臓との境界には、信号強度の段差が現れるので、この信号強度の段差を検出することにより、肝臓のＳＩ方向のエッジの位置Ｐ_１〜Ｐ_ｘを検出することができる。信号強度データＤ_１〜Ｄ_ｘごとに肝臓のエッジ位置を検出することにより、肝臓のエッジ位置が時間とともにどのように変化するかを知ることができる。したがって、肝臓のエッジ位置の時間変化を表すエッジ位置プロファイルを得ることができる。図２１に、エッジ位置プロファイルを概略的に示す。プリスキャンＡを実行した後、ステップＳＴ２に進む。

ステップＳＴ２では、エッジ位置プロファイルに基づいて、被検体１１の呼吸の吐き終わりの範囲ＡＷを求める。図２２に吐き終りの範囲ＡＷの一例を示す。呼吸の吐き終りの範囲ＡＷは、例えば、肝臓のエッジ位置の最小値から一定の範囲として求めることができる。呼吸の吐き終りの範囲ＡＷを求めた後、ステップＳＴ３に進む。

ステップＳＴ３では、本スキャンＢが実行される。
図２３は、本スキャンＢの説明図である。
本スキャンＢでは、ナビゲータシーケンスＮＳ_ｎ（ｎ＝１〜ｚ）と、イメージングシーケンスＩＳ_ｎ（ｎ＝１〜ｚ）とが交互に実行される。

本スキャンＢでは、ナビゲータシーケンスで検出された肝臓のエッジ位置が、呼吸の吐き終わりの範囲ＡＷに含まれている場合、その直後のイメージングシーケンスで取得されたイメージングデータを、画像再構成のデータとして採用する。一方、ナビゲータシーケンスで検出された肝臓のエッジ位置が、呼吸の吐き終わりの範囲ＡＷに含まれていない場合、その直後のイメージングシーケンスで取得されたイメージングデータは、画像再構成のデータとしては採用されない。例えば、図２３では、ナビゲータシーケンスＮＳ_１により検出された肝臓のエッジ位置Ｐ_１は、呼吸の吐き終わりの範囲ＡＷから外れているので、その直後のイメージングシーケンスＩＳ_１で取得されたイメージングデータＳ_１は、画像再構成のデータとしては採用されず、破棄される。

一方、ナビゲータシーケンスＮＳ_２により検出された肝臓のエッジ位置Ｐ_２は、呼吸の吐き終わりの範囲ＡＷに含まれているので、その直後のイメージングシーケンスＩＳ_２で取得されたイメージングデータＳ_２は、画像再構成のデータとして採用される。同様に、ナビゲータシーケンスＮＳ_３により検出された肝臓のエッジ位置Ｐ_３も、呼吸の吐き終わりの範囲ＡＷに含まれているので、その直後のイメージングシーケンスＩＳ_３で取得されたイメージングデータＳ_３も、画像再構成のデータとして採用される。
本スキャンＢを実行した後、ステップＳＴ４に進む。

ステップＳＴ４では、信号強度データＤ_１〜Ｄ_ｚの中から、肝臓のエッジの位置が呼吸の吐き終りの範囲ＡＷに含まれている信号強度データを特定する。信号強度データＤ_１〜Ｄ_ｚにおける肝臓のエッジの位置Ｐ_１〜Ｐ_ｚは既知であり、呼吸の吐き終りの範囲ＡＷも既知であるので、肝臓のエッジの位置Ｐ_１〜Ｐ_ｚが呼吸の吐き終りの範囲ＡＷに含まれているか否かを判断することができる。したがって、信号強度データＤ_１〜Ｄ_ｚの中から、肝臓のエッジの位置が呼吸の吐き終りの範囲ＡＷに含まれている信号強度データを特定することができる。図２３では、例えば、信号強度データＤ_１は、肝臓のエッジの位置Ｐ_１が呼吸の吐き終わりの範囲ＡＷから外れているが、信号強度データＤ_２、Ｄ_３、Ｄ_ｚなどは、肝臓のエッジの位置が呼吸の吐き終わりの範囲ＡＷに含まれている。信号強度データを特定した後、ステップＳＴ５に進む。

ステップＳＴ５では、検出手段９１ｂ（図１参照）が、ステップＳＴ４で特定された信号強度データに対して、肝臓のエッジの位置から肝臓側（Ｉ側）に距離Δｄだけ離れた位置における信号強度を検出する。図２４に、検出された信号強度を示す。図２４では、代表して、信号強度データＤ_２、Ｄ_３、Ｄ_ｚに対して検出された信号強度Ｎ_２、Ｎ_３、Ｎ_ｚが示されている。信号強度Ｎ_２、Ｎ_３、Ｎ_ｚは、それぞれ、肝臓のエッジの位置Ｐ_２、Ｐ_３、Ｐ_ｚから肝臓側（Ｉ側）に距離Δｄだけ離れた位置における信号強度を表している。信号強度を検出した後、ステップＳＴ６に進む。

ステップＳＴ６では、画像再構成用のデータとして採用されたイメージングデータを補正する。この補正は、以下の式を用いて行われる。
Ｓ_{ｎ，ｃｏｒｒ}＝Ｓ_ｎ＊Ｎ_２／Ｎ_ｎ ・・・（３）

Ｓ_ｎ：画像再構成用のデータとして採用されたｎ番目のイメージングデータ
Ｓ_{ｎ，ｃｏｒｒ}：Ｓ_ｎを補正した後のイメージングデータ
Ｎ_２：データＤ_２の位置Ｐ_２＋Δｄにおける信号強度
Ｎ_ｎ：データＤ_ｎの位置Ｐ_ｎ＋Δｄにおける信号強度

演算手段９１ｃ（図１１参照）は、式（３）を用いて、補正後のイメージングデータＳ_{ｎ，ｃｏｒｒ}を計算する。以下に、式（３）を用いてイメージングデータを補正する方法について説明する。

図２５は、式（３）を用いてイメージングデータを補正するときの説明図である。
イメージングデータＳ_２を補正する場合、式（３）のｎに、ｎ＝２を代入すればよい。
Ｓ_{２，ｃｏｒｒ}＝Ｓ_２＊Ｎ_２／Ｎ_２
＝Ｓ_２

したがって、補正後のイメージングデータＳ_{２，ｃｏｒｒ}は、補正前のイメージングデータＳ_２と同じである。

しかし、イメージングデータＳ_３を補正する場合、式（３）のｎに、ｎ＝３が代入されるので、補正後のイメージングデータＳ_{３，ｃｏｒｒ}は、以下の式で表される。
Ｓ_{３，ｃｏｒｒ}＝Ｓ_３＊Ｎ_２／Ｎ_３
したがって、補正後のイメージングデータＳ_{３，ｃｏｒｒ}は、補正前のイメージングデータＳ_３のＮ_２／Ｎ_３倍になる。

以下同様に、式（３）のｎの値を変更し、画像再構成用のデータとして採用されたイメージングデータＳ_ｎをＮ_２／Ｎ_ｎ倍することにより、補正後のイメージングデータＳ_{ｎ，ｃｏｒｒ}を得ることができる。イメージングデータを補正した後、ステップＳＴ７に進む。

ステップＳＴ７では、補正後のイメージングデータを用いて画像再構成を行い、フローを終了する。

第３の形態では、画像再構成用のデータとして採用されたイメージングデータＳ_ｎをＮ_２／Ｎ_ｎ倍している。したがって、肝臓に流入する造影剤の濃度の違いにより生じるイメージングデータＳ_ｎの信号強度の差を補正することができるので、リンギングやブラーなどのアーチファクトを低減することができる。

また、第３の形態では、肝臓のエッジの位置が呼吸の吐き終りの範囲ＡＷに含まれているときに取得されたイメージングデータのみを用いて画像再構成が行われるので、自由呼吸下で被検体を撮影しても、ゴーストなどの体動アーチファクトを低減することができる。

尚、式（３）の信号強度Ｎ_２の代わりに、別の信号強度Ｎ_ｉを用いてイメージングデータを補正してもよい。

（４）第４の形態
第４の形態では、被検体に自由に呼吸をさせた状態でスキャンを実行する別の例について説明する。尚、ＭＲ装置のハードウェア構成は、第１の形態と同じである。

図２６は、第４の形態で実行されるスキャンの説明図である。
第４の形態では、プリスキャンＡと、本スキャンＢとが実行される。

プリスキャンＡでは、第３の形態と同様に、ナビゲータシーケンスＮＳ_１〜ＮＳ_ｘが実行される。

本スキャンＢでは、肝臓を含む撮影領域の画像を再構成するときに使用されるイメージングデータを収集するためのイメージングシーケンスＩＳ_１〜ＩＳ_ｗが実行される。また、イメージングシーケンスＩＳ_１〜ＩＳ_ｗの他に、ナビゲータシーケンスＮＳが実行される。ナビゲータシーケンスＮＳは、各イメージングシーケンスＩＳ_１〜ＩＳ_ｗの直前に１回以上実行される。

以下に、プリスキャンＡおよび本スキャンＢにより被検体を撮影するときのフローについて説明する。

図２７は、第４の形態において被検体を撮影するときのフローを示す図である。
ステップＳＴ１では、プリスキャンＡが実行される。プリスキャンＡを実行することにより、第３の形態と同様に、肝臓のエッジ位置の時間変化を表すエッジ位置プロファイル（図２１参照）を得ることができる。プリスキャンＡを実行した後、ステップＳＴ２に進む。

ステップＳＴ２では、エッジ位置プロファイルに基づいて、トリガレベルＴＬが設定される。図２８に、設定されたトリガレベルＴＬの一例を示す。トリガレベルＴＬは、本スキャンＢを実行するときに使用されるものであり、具体的には、本スキャンＢにおいてイメージングシーケンスＩＳ_１〜ＩＳ_ｗ（図２６参照）を実行するときの肝臓のエッジの基準位置を表している。トリガレベルＴＬが本スキャンＢを実行するときにどのように使用されるかについては、後述する。トリガレベルＴＬは、例えば、肝臓のエッジの位置の最大値と最小値の中間に設定することができる。トリガレベルＴＬを設定した後、ステップＳＴ３に進む。

ステップＳＴ３では、本スキャンＢが実行される。
図２９は、本スキャンＢの説明図である。
本スキャンＢでは、先ず、ナビゲータシーケンスＮＳ_１１〜ＮＳ_１Ｐが実行される。ナビゲータシーケンスＮＳ_１１〜ＮＳ_１Ｐにより収集されたナビゲータ信号Ａ_１１〜Ａ_１Ｐに基づいてナビゲータ領域Ｒ_ｎａｖの信号強度データＤ_１１〜Ｄ_１Ｐを求め、信号強度データＤ_１１〜Ｄ_１Ｐごとに肝臓のエッジの位置Ｐ_１１〜Ｐ_１Ｐを検出する。そして、肝臓のエッジの位置がトリガレベルＴＬの下側（Ｉ側）から上側（Ｓ側）に移動したときに、イメージングシーケンスＩＳ_１を実行する。イメージングシーケンスＩＳ_１を実行することにより、撮影領域ＲのイメージングデータＳ_１が収集される。

イメージングシーケンスＩＳ_１を実行した後、ナビゲータシーケンスＮＳ_２１〜ＮＳ_２ｑが実行される。そして、肝臓のエッジの位置がトリガレベルＴＬの下側（Ｉ側）から上側（Ｓ側）に移動したときに、次のイメージングシーケンスＩＳ_２を実行する。イメージングシーケンスＩＳ_２を実行することにより、撮影領域ＲのイメージングデータＳ_２が収集される。

以下同様にナビゲータシーケンスとイメージングシーケンスを実行する。最後にナビゲータシーケンスＮＳ_ｚ１〜ＮＳ_ｚｒを実行し、肝臓のエッジの位置がトリガレベルＴＬの下側（Ｉ側）から上側（Ｓ側）に移動したときに、イメージングシーケンスＩＳ_ｗを実行する。イメージングシーケンスＩＳ_ｗを実行することにより、撮影領域ＲのイメージングデータＳ_ｗが収集され、本スキャンＢが終了する。本スキャンＢが終了した後、ステップＳＴ４に進む。

ステップＳＴ４では、検出手段９１ｂ（図１参照）は、イメージングシーケンスＩＳ_１〜ＩＳ_ｗの直前に取得された信号強度データＤ_１ｐ〜Ｄ_ｗｒに対して、肝臓のエッジの位置から肝臓側（Ｉ側）に距離Δｄだけ離れた位置における信号強度を検出する。図３０に、検出された信号強度を示す。図３０では、代表して、信号強度データＤ_１ｐ、Ｄ_２ｑ、Ｄ_ｗｒに対して検出された信号強度Ｎ_１、Ｎ_２、Ｎ_ｗが示されている。信号強度Ｎ_１、Ｎ_２、Ｎ_ｗは、それぞれ、肝臓のエッジの位置Ｐ_１ｐ、Ｐ_２ｑ、Ｐ_ｗｒから肝臓側（Ｉ側）に距離Δｄだけ離れた位置における信号強度を表している。信号強度を検出した後、ステップＳＴ５に進む。

ステップＳＴ５では、イメージングデータＳ_１〜Ｓ_ｗを補正する。この補正は、以下の式を用いて行われる。
Ｓ_{ｎ，ｃｏｒｒ}＝Ｓ_ｎ＊Ｎ_１／Ｎ_ｎ ・・・（４）

Ｓ_ｎ：ｎ番目のイメージングデータ
Ｓ_{ｎ，ｃｏｒｒ}：Ｓ_ｎを補正した後のイメージングデータ
Ｎ_１：１番目のイメージングデータＳ_１の直前の信号強度データから検出された信号強度
Ｎ_ｎ：ｎ番目のイメージングデータＳ_ｎの直前の信号強度データから検出された信号強度

以下に、式（４）を用いてイメージングデータを補正する方法について説明する。
図３１は、式（４）を用いてイメージングデータを補正するときの説明図である。
イメージングデータＳ_１を補正する場合、式（４）のｎに、ｎ＝１を代入すればよい。
Ｓ_{１，ｃｏｒｒ}＝Ｓ_１＊Ｎ_１／Ｎ_１
＝Ｓ_１

したがって、補正後のイメージングデータＳ_{１，ｃｏｒｒ}は、補正前のイメージングデータＳ_１と同じである。

しかし、イメージングデータＳ_２を補正する場合、式（４）のｎに、ｎ＝２が代入されるので、補正後のイメージングデータＳ_{２，ｃｏｒｒ}は、以下の式で表される。
Ｓ_{２，ｃｏｒｒ}＝Ｓ_２＊Ｎ_１／Ｎ_２

したがって、補正後のイメージングデータＳ_{２，ｃｏｒｒ}は、補正前のイメージングデータＳ_２のＮ_１／Ｎ_２倍になる。

以下同様に、式（４）のｎの値を変更し、イメージングデータＳ_ｎをＮ_１／Ｎ_ｎ倍することにより、補正後のイメージングデータＳ_{ｎ，ｃｏｒｒ}を得ることができる。イメージングデータを補正した後、ステップＳＴ６に進む。

ステップＳＴ６では、補正後のイメージングデータを用いて画像再構成を行い、フローを終了する。

第４の形態では、イメージングデータＳ_ｎをＮ_１／Ｎ_ｎ倍しているので、肝臓に流入する造影剤の濃度の違いにより生じるイメージングデータＳ_ｎの信号強度の差を補正することができるので、リンギングやブラーなどのアーチファクトを低減することができる。

また、第４の形態では、肝臓のエッジの位置がトリガレベルＴＬに到達したときにイメージングデータを取得するので、イメージングデータＳ_１〜Ｓ_Ｗを取得するときの肝臓のエッジの位置のばらつきを小さくすることができる。したがって、ゴーストなどの体動アーチファクトを低減することができる。

尚、式（４）の信号強度Ｎ_１の代わりに、別の信号強度Ｎ_ｉを用いてイメージングデータを補正してもよい。

（５）第５の形態
第１〜第４の形態では、イメージングデータを補正することによりリンギングやブラーなどのアーチファクトを低減する例について説明したが、第５の形態では、ＲＦパルスのフリップ角を補正することによりリンギングやブラーなどのアーチファクトを低減する例について説明する。

図３２は、第５の形態のＭＲ装置の概略図である。
第５の形態のＭＲ装置３００は、第１の形態のＭＲ装置１００の補正手段９１（図１参照）の代わりに、ＲＦパルスのフリップ角を補正する補正手段９１１を有している。尚、その他の構成は、第１の形態のＭＲ装置１００と同じである。制御部９は、補正手段９１１および再構成手段９２を構成する一例であり、所定のプログラムを実行することにより、これらの手段として機能する。

図３３は、第５の形態で実行されるスキャンを示す図である。
第５の形態では、息止めスキャンＢＨが実行される。息止めスキャンＢＨでは、ナビゲータシーケンスＮＳ_ｎ（ｎ＝１〜ｚ）と、イメージングシーケンスＩＳ_ｎ（ｎ＝１〜ｚ）とが交互に実行される。図３３の下側に、イメージングシーケンスＩＳ_ｎの一例が示されている。ここでは、イメージングシーケンスＩＳ_ｎの一例として、ＳＰＧＲ（Spoiled Gradient Recalled）シーケンスが示されている。ＲＦパルスＸ_ｎのフリップ角αは、例えば９０°である。
以下に、息止めスキャンＢＨを実行する方法について説明する。

息止めスキャンＢＨでは、先ず、ナビゲータシーケンスＮＳ_１とイメージングシーケンスＩＳ_１とが実行される。図３４に、ナビゲータシーケンスＮＳ_１により収集されたナビゲータ信号Ａ_１と、イメージングシーケンスＩＳ_１により収集されたイメージングデータＳ_１とを示す。信号強度データ作成手段９１１ａ（図３２参照）は、ナビゲータ信号Ａ_１を一次元フーリエ変換し、これにより、ナビゲータ領域Ｒ_ｎａｖ内のＳＩ方向の信号強度の変化を表す信号強度データＤ_１が得られる。

検出手段９１１ｂ（図３２参照）は、信号強度データＤ_１の信号強度の段差を検出し、肝臓のエッジのＳＩ方向の位置Ｐ_１を検出する。そして、肝臓のエッジの位置Ｐ_１から肝臓側（Ｉ側）に距離Δｄだけ離れた位置Ｐ_１＋Δｄの信号強度Ｎ_１を検出する。

イメージングシーケンスＩＳ_１を実行した後、ナビゲータシーケンスＮＳ_２が実行される。図３５に、ナビゲータシーケンスＮＳ_２により収集されたナビゲータ信号Ａ_２を示す。ナビゲータ信号Ａ_２は一次元フーリエ変換され、ナビゲータ領域Ｒ_ｎａｖ内のＳＩ方向の信号強度の変化を表す信号強度データＤ_２が得られる。そして、信号強度データＤ_２の位置Ｐ_１＋Δｄにおける信号強度Ｎ_２を検出する。

次に、イメージングシーケンスＩＳ_２のＲＦパルスＸ_２のフリップ角αを補正する。以下に、図３６を参照しながら、フリップ角αの補正方法について説明する。フリップ角αは、以下の式を用いて補正される。
α_{ｎ_ｃｏｒｒ}＝arcsin｛（Ｎ_１／Ｎ_ｎ）＊sinα｝ ・・・（５）

α：イメージングシーケンスＩＳ_ｎの補正前のフリップ角
α_{ｎ_ｃｏｒｒ}：イメージングシーケンスＩＳ_ｎの補正後のフリップ角
Ｎ_１：信号強度データＤ_１の位置Ｐ_１＋Δｄにおける信号強度
Ｎ_ｎ：信号強度データＤ_ｎの位置Ｐ_１＋Δｄにおける信号強度

演算手段９１１ｃ（図３２参照）は、式（５）を用いて、補正後のフリップ角α_{ｎ，ｃｏｒｒ}を計算する。ここでは、ｎ＝２、つまり、イメージングシーケンスＩＳ_２のＲＦパルスＸ_２のフリップ角αを補正するので、式（５）のｎに、ｎ＝２を代入すればよい。
α_{２_ｃｏｒｒ}＝arcsin｛（Ｎ_１／Ｎ_２）＊sinα｝

したがって、補正後のフリップ角α_{２_ｃｏｒｒ}を求めることができる。補正後のフリップ角α_{２_ｃｏｒｒ}を求めたら、イメージングシーケンスＩＳ_２のＲＦパルスＸ_２のフリップ角αをα_{２_ｃｏｒｒ}に補正し、撮影領域Ｒのスキャンを実行する。したがって、フリップ角α＝α_{２_ｃｏｒｒ}のときのイメージングデータＳ_２を得ることができる。

以下同様に、ナビゲータシーケンスＮＳ_ｎを実行するたびに、式（５）を用いて補正後のフリップ角α_{ｎ_ｃｏｒｒ}を計算し、イメージングシーケンスＩＳ_ｎを実行する。したがって、フリップ角α＝α_{ｎ_ｃｏｒｒ}のときのイメージングデータＳ_ｎを得ることができる。図３７に、息止めスキャンＢＨにより得られたイメージングデータＳ_１〜Ｓ_ｚを示す。イメージングデータＳ_１〜Ｓ_ｚを収集した後、これらのイメージングデータＳ_１〜Ｓ_ｚを用いて画像再構成を行う。これにより、撮影領域Ｒの画像データを得ることができる。

ＳＰＧＲシーケンスでは、（Ｎ_１／Ｎ_ｎ）＊sinαのarcsineを求めることにより、造影剤の濃度の違いが補正されたイメージングデータＳ_１〜Ｓ_ｚを収集することができるので、リンギングやブラーなどのアーチファクトを低減することができる。

（６）第６の形態
第５の形態では、ナビゲータ信号とイメージングデータとを別々のシーケンスで収集する例について示されているが、第６の形態では、ナビゲータ信号とイメージングデータとを一つのシーケンスＣ_ｎで収集する例について説明する。

図３８は、第６の形態のＭＲ装置を示す概略図である。
第６の形態のＭＲ装置４００は、第５の形態のＭＲ装置３００（図３２参照）の信号強度データ作成手段９１１ａおよび検出手段９１１ｂの代わりに、信号強度算出手段９１１ｄを備えているが、その他の構成は、第５の形態と同じである。

図３９は、第６の形態で実行されるスキャンを示す図である。
第６の形態では、息止めスキャンＢＨが実行される。息止めスキャンＢＨでは、第２の形態と同様に、シーケンスＣ_ｎ（ｎ＝１〜ｚ）が実行される。

先ず、シーケンスＣ_１が実行される。図４０に、シーケンスＣ_１により収集されたナビゲータ信号Ａ_１およびイメージングデータＳ_１を示す。シーケンスＣ_１を実行することによりナビゲータ信号Ａ_１およびイメージングデータＳ_１が収集される。信号強度算出手段９１１ｄ（図３８参照）は、撮影領域Ｒ（図１３参照）から得られたナビゲータ信号Ａ_１の積分値を計算する。これにより、撮影領域Ｒの信号強度Ｎ_１を算出することができる。

シーケンスＣ_１を実行した後、次のシーケンスＣ_２が実行される。図４１に、シーケンスＣ_２により収集されたナビゲータ信号Ａ_２およびイメージングデータＳ_２を示す。信号強度算出手段９１１ｄは、撮影領域Ｒから得られたナビゲータ信号Ａ_２の積分値を計算する。これにより、撮影領域Ｒの信号強度Ｎ_２を算出することができる。

シーケンスＣ_２を実行した後、シーケンスＣ_３が実行される。ただし、シーケンスＣ_３を実行する前に、シーケンスＣ_３のＲＦパルスＸ_３のフリップ角αが補正される。以下に、シーケンスＣ_３のＲＦパルスＸ_３のフリップ角αの補正方法について述べる。

図４２は、シーケンスＣ_３のＲＦパルスＸ_３のフリップ角αの補正方法の説明図である。
シーケンスＣ_３のＲＦパルスＸ_３のフリップ角αは、以下の式を用いて補正される。
α_{３_ｃｏｒｒ}＝arcsin｛（Ｎ_１／Ｎ_２）＊sinα｝ ・・・（６）

α：シーケンスＣ_３のＲＦパルスＸ_３の補正前のフリップ角
α_{３_ｃｏｒｒ}：シーケンスＣ_３のＲＦパルスＸ_３の補正後のフリップ角
Ｎ_１：ナビゲータ信号Ａ_１により得られた撮影領域Ｒの信号強度
Ｎ_２：ナビゲータ信号Ａ_２により得られた撮影領域Ｒの信号強度

演算手段９１１ｃ（図３８参照）は、式（６）を用いて、補正後のフリップ角α_{３_ｃｏｒｒ}を計算する。式（６）の信号強度Ｎ_１およびＮ_２は既に算出されており、フリップ角αは事前に決まっている値である。したがって、式（６）にＮ_１、Ｎ_２、およびα_３を代入することにより、補正後のフリップ角α_{３_ｃｏｒｒ}を求めることができる。補正後のフリップ角α_{３_ｃｏｒｒ}を求めたら、シーケンスＣ_３のＲＦパルスＸ_３のフリップ角αをα_{３_ｃｏｒｒ}に補正し、シーケンスＣ_３を実行する。したがって、ナビゲータ信号Ａ_３とイメージングデータＳ_３とを得ることができる。信号強度算出手段９１１ｄは、撮影領域Ｒから得られたナビゲータ信号Ａ_３の積分値を計算し、撮影領域Ｒの信号強度Ｎ_３を算出する。

シーケンスＣ_３を実行した後、シーケンスＣ_４が実行される。ただし、シーケンスＣ_４を実行する前に、シーケンスＣ_４のＲＦパルスＸ_４のフリップ角αが補正される。以下に、シーケンスＣ_４のＲＦパルスＸ_４のフリップ角αの補正方法について述べる。

シーケンスＣ_３のフリップ角は、αからα_{３＿ｃｏｒｒ}に補正されているので、シーケンスＣ_３の信号強度Ｎ_３は、フリップ角α_{３＿ｃｏｒｒ}により得られた値である。一方、シーケンスＣ_１の信号強度Ｎ_１はフリップ角αにより得られた値である。したがって、信号強度Ｎ_１およびＮ_３は異なるフリップ角により得られている。そこで、信号強度Ｎ_３を、信号強度Ｎ_１と同様に、フリップ角αにより得られた値に補正する。信号強度Ｎ_３は、以下の式（７）を用いて補正することができる。
Ｎ_{ｎ＿ｃｏｒｒ}＝Ｎ_ｎ（sinα／sinα_{ｎ＿ｃｏｒｒ}） ・・・（７）

α：シーケンスＣ_ｎのＲＦパルスＸ_ｎの補正前のフリップ角
α_{ｎ＿ｃｏｒｒ}：シーケンスＣ_ｎのＲＦパルスＸ_ｎの補正後のフリップ角
Ｎ_ｎ：ナビゲータ信号Ａ_ｎにより得られた撮影領域Ｒの信号強度
Ｎ_{ｎ＿ｃｏｒｒ}：Ｎ_ｎの補正後の値

ここでは、ｎ＝３、つまり、信号強度Ｎ_３を補正するので、式（７）のｎに、ｎ＝３を代入すればよい。
Ｎ_{３＿ｃｏｒｒ}＝Ｎ_ｎ（sinα／sinα_{３＿ｃｏｒｒ}）

したがって、フリップ角α_{３＿ｃｏｒｒ}により得られた信号強度Ｎ_３を、フリップ角αにより得られた信号強度Ｎ_{３＿ｃｏｒｒ}に補正することができる。図４３に、補正後の信号強度Ｎ_{３＿ｃｏｒｒ}が算出された様子を示す。補正後の信号強度Ｎ_{３＿ｃｏｒｒ}を求めた後、次のシーケンスＣ_４のフリップ角を補正する。

図４４は、シーケンスＣ_４のフリップ角を補正するときの説明図である。
シーケンスＣ_４のＲＦパルスＸ_４のフリップ角αは、以下の式を用いて補正される。
α_{ｎ＋１_ｃｏｒｒ}＝arcsin｛（Ｎ_１／Ｎ_{ｎ＿ｃｏｒｒ}）＊sinα｝ ・・・（８）

演算手段９１１ｃは、式（８）を用いて、補正後のフリップ角α_{ｎ_ｃｏｒｒ}を計算する。ここでは、シーケンスＣ_４のＲＦパルスＸ_４のフリップ角αを補正するので、式（８）のｎに、ｎ＝３を代入すればよい。
α_{４_ｃｏｒｒ}＝arcsin｛（Ｎ_１／Ｎ_{３＿ｃｏｒｒ}）＊sinα｝

したがって、補正後のフリップ角α_{４_ｃｏｒｒ}を求めることができる。補正後のフリップ角α_{４_ｃｏｒｒ}を求めたら、シーケンスＣ_４のＲＦパルスＸ_４のフリップ角αをα_{４_ｃｏｒｒ}に補正し、シーケンスＣ_４を実行する。したがって、ナビゲータ信号Ａ_４とイメージングデータＳ_４とを得ることができる。信号強度算出手段９１１ｄは、撮影領域Ｒから得られたナビゲータ信号Ａ_４の積分値を計算し、撮影領域Ｒの信号強度Ｎ_４を算出する。

以下同様に、シーケンスＣ_ｎを実行するたびに、式（７）を用いて信号強度Ｎ_ｎを補正してから、式（８）を用いて補正後のフリップ角α_{ｎ＋１_ｃｏｒｒ}を求め、シーケンスを実行する。最後のシーケンスＣ_ｚが終了したら、画像再構成を行う。これにより、撮影領域Ｒの画像データを得ることができる。

第６の形態では、一つのシーケンスＣ_ｎで、ナビゲータ信号Ａ_ｎとイメージングデータＳ_ｎとの両方を収集している。この場合は、式（７）を用いて信号強度を補正した後で、フリップ角を補正することにより、息止めスキャンＢＨにより得られるイメージングデータの信号強度の差を小さくすることができるので、リンギングやブラーなどのアーチファクトを低減することができる。

第５および第６の形態では息止めスキャンの場合について説明したが、第３および第４の形態と同様に、被検体に自由に呼吸をさせた状態でスキャンを実行し、フリップ角を補正してもよい。

尚、第１〜第６の形態では、造影剤を用いてスキャンを行う場合について説明されているが、本発明は、イメージングデータの信号強度が時間とともに変化する場合であれば、造影剤を用いないスキャンを実行する場合にも適用できる。

更に、第１〜第６の形態では、体動する部位の画像再構成を行う場合について説明されているが、本発明は、体動しない部位の画像再構成を行う場合にも適用することができる。

２ マグネット
３ テーブル
３ａ クレードル
４ 受信コイル
５ 造影剤注入装置
６ 送信器
７ 勾配磁場電源
８ 受信器
９ 制御部
１０ 操作部
１１ 表示部
１２ 被検体
２１ ボア
２２ 超伝導コイル
２３ 勾配コイル
２４ ＲＦコイル
９１ａ、９１１ａ 信号強度データ作成手段
９１ｂ、９１１ｂ 検出手段
９１ｃ、９１１ｃ 演算手段
９２ 再構成手段
１００、２００、３００、４００ ＭＲ装置

Claims (33)

1. 被検体の所定の部位を含む第１の領域からナビゲータ信号を収集するとともに、前記所定の部位を含む第２の領域から前記所定の部位の画像を再構成するときに使用されるイメージングデータを収集する磁気共鳴装置であって、
   前記ナビゲータ信号と前記イメージングデータとを収集するためのスキャンを実行するスキャン手段と、
   前記ナビゲータ信号に基づいて、前記イメージングデータの信号強度を補正する補正手段と、
   信号強度が補正されたイメージングデータに基づいて画像再構成を行う再構成手段と、
   を有する磁気共鳴装置。
2. 前記補正手段は、
   前記ナビゲータ信号に基づいて前記第１の領域の信号強度データを作成する信号強度データ作成手段と、
   前記信号強度データの中から前記所定の部位の信号強度を検出する検出手段とを有し、
   前記信号強度に基づいて前記イメージングデータの信号強度を補正する、請求項１に記載の磁気共鳴装置。
3. 前記検出手段は、
   前記信号強度データに基づいて前記所定の部位のエッジの位置を検出し、前記所定の部位のエッジの位置から所定の距離離れた位置の信号強度を、前記所定の部位の信号強度として検出する、請求項２に記載の磁気共鳴装置。
4. 前記スキャン手段は、
   前記ナビゲータ信号と前記イメージングデータとを交互に収集し、
   前記補正手段は、
   ｉ番目のナビゲータ信号により得られた信号強度と、ｎ番目のナビゲータ信号により得られた信号強度とに基づいて、前記ｎ番目のナビゲータ信号の直前又は直後に収集されたイメージングデータを補正する、請求項２又は３に記載の磁気共鳴装置。
5. 前記スキャン手段は、
   前記ナビゲータ信号と複数のイメージングデータとを交互に実行し、
   前記補正手段は、
   ｉ番目のナビゲータ信号により得られた信号強度と、ｎ番目のナビゲータ信号により得られた信号強度とに基づいて、前記ｎ番目のナビゲータ信号とｎ＋１番目のナビゲータ信号との間に収集されたイメージングデータ、又は前記ｎ番目のナビゲータ信号とｎ−１番目のナビゲータ信号との間に収集されたイメージングデータを補正する、請求項２又は３に記載の磁気共鳴装置。
6. 前記ｉ番目のナビゲータ信号は、１番目のナビゲータ信号である、請求項４又は５に記載の磁気共鳴装置。
7. 前記所定の部位の信号強度の時間変化を表すプロファイルを推定する推定手段を有する、請求項１〜６のうちのいずれか一項に記載の磁気共鳴装置。
8. 前記補正手段は、
   前記ナビゲータ信号に基づいて前記第１の領域の信号強度を算出する信号強度算出手段を有し、
   前記信号強度に基づいて前記イメージングデータの信号強度を補正する、請求項１に記載の磁気共鳴装置。
9. 前記スキャン手段は、
   前記ナビゲータ信号と前記イメージングデータとを交互に収集し、
   前記補正手段は、
   ｉ番目のナビゲータ信号により得られた信号強度と、ｎ番目のナビゲータ信号により得られた信号強度とに基づいて、前記ｎ番目のナビゲータ信号の直前又は直後に収集されたイメージングデータを補正する、請求項８に記載の磁気共鳴装置。
10. 前記スキャン手段は、
    前記ナビゲータ信号と複数のイメージングデータとを交互に実行し、
    前記補正手段は、
    ｉ番目のナビゲータ信号により得られた信号強度と、ｎ番目のナビゲータ信号により得られた信号強度とに基づいて、前記ｎ番目のナビゲータ信号とｎ＋１番目のナビゲータ信号との間に収集されたイメージングデータ、又は前記ｎ番目のナビゲータ信号とｎ−１番目のナビゲータ信号との間に収集されたイメージングデータを補正する、請求項８に記載の磁気共鳴装置。
11. 前記ｉ番目のナビゲータ信号は、１番目のナビゲータ信号である、請求項９又は１０に記載の磁気共鳴装置。
12. 前記第１の領域は、前記第２の領域と同じである、請求項８〜１１のうちのいずれか一項に記載の磁気共鳴装置。
13. 前記所定の部位は、体動する部位であり、
    前記スキャン手段は、前記被検体が息止めをしている間に、前記ナビゲータ信号と前記イメージングデータとを収集するためのスキャンを実行する、請求項１〜１２のうちのいずれか一項に記載の磁気共鳴装置。
14. 前記所定の部位は、体動する部位であり、
    前記スキャン手段は、呼吸信号を取得するためのプリスキャンと、前記ナビゲータ信号と前記イメージングデータとを収集するための本スキャンとを実行し、
    前記ナビゲータ信号は、前記体動する部位のエッジの位置と、前記体動する部位の信号強度とを求めるために使用される、請求項１〜３のうちのいずれか一項に記載の磁気共鳴装置。
15. 前記呼吸信号に基づいて、前記イメージングデータを画像再構成用のデータとして採用するか否かを判断するときに使用される呼吸の範囲を設定し、
    本スキャンでは、前記ナビゲータ信号と前記イメージングデータとを交互に収集し、前記ナビゲータ信号により求められた前記体動する部位のエッジの位置が、前記呼吸の範囲に含まれている場合、前記ナビゲータ信号の直前又は直後に収集されたイメージングデータを、画像再構成用のデータとして採用する、請求項１４に記載の磁気共鳴装置。
16. 前記呼吸信号に基づいて前記体動する部位のエッジの基準位置を設定し、
    前記本スキャンでは、前記体動する部位のエッジの位置が前記基準位置に到達したときに、前記イメージングデータを収集する、請求項１４に記載の磁気共鳴装置。
17. 前記ナビゲータ信号および前記イメージングデータは別々のシーケンスで収集される、請求項１〜１６のうちのいずれか一項に記載の磁気共鳴装置。
18. 前記ナビゲータ信号および前記イメージングデータは同じシーケンスで収集される、請求項１〜１６のうちのいずれか一項に記載の磁気共鳴装置。
19. 被検体の所定の部位を含む第１の領域からナビゲータ信号を収集するとともに、前記所定の部位を含む第２の領域から前記所定の部位の画像を再構成するときに使用されるイメージングデータを収集する磁気共鳴装置であって、
    前記ナビゲータ信号と前記イメージングデータとを収集するためのスキャンを実行するスキャン手段と、
    前記ナビゲータ信号に基づいて、前記イメージングデータを収集するときに使用されるシーケンスのＲＦパルスのフリップ角を補正する補正手段と、
    を有する、磁気共鳴装置。
20. 前記補正手段は、
    前記ナビゲータ信号に基づいて前記第１の領域の信号強度データを作成する信号強度データ作成手段と、
    前記信号強度データの中から前記所定の部位の信号強度を検出する検出手段とを有し、
    前記信号強度に基づいて前記フリップ角を補正する、請求項１９に記載の磁気共鳴装置。
21. 前記検出手段は、
    前記信号強度データに基づいて前記所定の部位のエッジの位置を検出し、前記所定の部位のエッジの位置から所定の距離離れた位置の信号強度を、前記所定の部位の信号強度として検出する、請求項２０に記載の磁気共鳴装置。
22. 前記スキャン手段は、
    前記ナビゲータ信号を収集するためのナビゲータシーケンスと、前記イメージングデータを収集するためのイメージングシーケンスとを交互に実行し、
    前記補正手段は、
    ｉ番目のナビゲータ信号により得られた信号強度と、ｎ番目のナビゲータ信号により得られた信号強度とに基づいて、前記ｎ番目のナビゲータ信号の直後に実行されるイメージングシーケンスのＲＦパルスのフリップ角を補正する、請求項２０又は２１に記載の磁気共鳴装置。
23. 前記ｉ番目のナビゲータ信号は、１番目のナビゲータ信号である、請求項２２に記載の磁気共鳴装置。
24. 前記補正手段は、
    前記ナビゲータ信号に基づいて前記第１の領域の信号強度を算出する信号強度算出手段を有し、
    前記信号強度に基づいて前記フリップ角を補正する、請求項１９に記載の磁気共鳴装置。
25. 前記スキャン手段は、
    前記ナビゲータ信号と前記イメージングデータとを収集するためのシーケンスを複数回実行し、
    前記補正手段は、
    ｎ番目のシーケンスのナビゲータ信号により得られた信号強度を補正し、補正された信号強度に基づいて、ｎ＋１番目のシーケンスのＲＦパルスのフリップ角を補正する、請求項２４に記載の磁気共鳴装置。
26. 前記第１の領域は、前記第２の領域と同じである、請求項２４又は２５に記載の磁気共鳴装置。
27. 前記所定の部位は、体動する部位であり、
    前記スキャン手段は、前記被検体が息止めをしている間に、前記ナビゲータ信号と前記イメージングデータとを収集する、請求項１９〜２６のうちのいずれか一項に記載の磁気共鳴装置。
28. 前記所定の部位は、体動する部位であり、
    前記スキャン手段は、呼吸信号を取得するためのプリスキャンと、前記ナビゲータ信号と前記イメージングデータとを収集するための本スキャンとを実行し、
    前記ナビゲータ信号は、前記体動する部位のエッジの位置と、前記体動する部位の信号強度とを求めるために使用される、請求項１９〜２１のうちのいずれか一項に記載の磁気共鳴装置。
29. 前記呼吸信号に基づいて、前記イメージングデータを画像再構成用のデータとして採用するか否かを判断するときに使用される呼吸の範囲を設定し、
    本スキャンでは、前記ナビゲータ信号と前記イメージングデータとを交互に収集し、前記ナビゲータ信号により求められた前記体動する部位のエッジの位置が、前記呼吸の範囲に含まれている場合、前記ナビゲータ信号の直前又は直後に収集されたイメージングデータを、画像再構成用のデータとして採用する、請求項２８に記載の磁気共鳴装置。
30. 前記呼吸信号に基づいて前記体動する部位のエッジの基準位置を設定し、
    前記本スキャンでは、前記体動する部位のエッジの位置が前記基準位置に到達したときに、前記イメージングデータを収集する、請求項２８に記載の磁気共鳴装置。
31. 前記被検体に造影剤を注入し、前記ナビゲータ信号および前記イメージングデータとを収集する、請求項１〜３０のうちのいずれか一項に記載の磁気共鳴装置。
32. 被検体の所定の部位を含む第１の領域からナビゲータ信号を収集するとともに、前記所定の部位を含む第２の領域から前記所定の部位の画像を再構成するときに使用されるイメージングデータを収集するためのスキャンを実行する磁気共鳴装置のプログラムであって、
    前記ナビゲータ信号に基づいて、前記イメージングデータの信号強度を補正する補正処理と、
    信号強度が補正されたイメージングデータに基づいて画像再構成を行う再構成処理と、
    を計算機に実行させるためのプログラム。
33. 被検体の所定の部位を含む第１の領域からナビゲータ信号を収集するとともに、前記所定の部位を含む第２の領域から前記所定の部位の画像を再構成するときに使用されるイメージングデータを収集するためのスキャンを実行する磁気共鳴装置のプログラムであって、
    前記ナビゲータ信号に基づいて、前記イメージングデータを収集するときに使用されるシーケンスのＲＦパルスのフリップ角を補正する補正処理、
    を計算機に実行させるためのプログラム。
    
    

Images