(57)【要約】      （修正有）【課題】リンギングやブラーなどのアーチファクトを低減する磁気共鳴装置を提供する。【解決手段】イメージングシーケンスＩＳｎの他に、ナビゲータシーケンスＮＳｎを実行する。ナビゲータシーケンスＮＳｎは、肺と肝臓との境界部分を含むナビゲータ領域Ｒｎａｖ内のＳＩ方向の信号強度の変化を表す信号強度データを取得するためのシーケンスである。ナビゲータシーケンスＮＳｎにより得られた信号強度データＤ１〜Ｄｚごとに、肝臓のエッジの位置Ｐ０から肝臓側（Ｉ側）に距離Δｄだけ離れた位置Ｐ０＋Δｄの信号強度Ｎ１〜Ｎｚを検出する。そして、信号強度Ｎ１〜Ｎｚに基づいて、イメージングシーケンスにより得られた撮影領域のイメージングデータＳ１〜Ｓｚを補正する。【選択図】図９
2014023703.tif 000002
【特許請求の範囲】
【請求項１】
被検体の所定の部位を含む第１の領域からナビゲータ信号を収集するとともに、前記所定の部位を含む第２の領域から前記所定の部位の画像を再構成するときに使用されるイメージングデータを収集する磁気共鳴装置であって、
前記ナビゲータ信号と前記イメージングデータとを収集するためのスキャンを実行するスキャン手段と、
前記ナビゲータ信号に基づいて、前記イメージングデータの信号強度を補正する補正手段と、
信号強度が補正されたイメージングデータに基づいて画像再構成を行う再構成手段と、
を有する磁気共鳴装置。
【請求項２】
前記補正手段は、
前記ナビゲータ信号に基づいて前記第１の領域の信号強度データを作成する信号強度データ作成手段と、
前記信号強度データの中から前記所定の部位の信号強度を検出する検出手段とを有し、
前記信号強度に基づいて前記イメージングデータの信号強度を補正する、請求項１に記載の磁気共鳴装置。
【請求項３】
前記検出手段は、
前記信号強度データに基づいて前記所定の部位のエッジの位置を検出し、前記所定の部位のエッジの位置から所定の距離離れた位置の信号強度を、前記所定の部位の信号強度として検出する、請求項２に記載の磁気共鳴装置。
【請求項４】
前記スキャン手段は、
前記ナビゲータ信号と前記イメージングデータとを交互に収集し、
前記補正手段は、
ｉ番目のナビゲータ信号により得られた信号強度と、ｎ番目のナビゲータ信号により得られた信号強度とに基づいて、前記ｎ番目のナビゲータ信号の直前又は直後に収集されたイメージングデータを補正する、請求項２又は３に記載の磁気共鳴装置。
【請求項５】
前記スキャン手段は、
前記ナビゲータ信号と複数のイメージングデータとを交互に実行し、
前記補正手段は、
ｉ番目のナビゲータ信号により得られた信号強度と、ｎ番目のナビゲータ信号により得られた信号強度とに基づいて、前記ｎ番目のナビゲータ信号とｎ＋１番目のナビゲータ信号との間に収集されたイメージングデータ、又は前記ｎ番目のナビゲータ信号とｎ−１番目のナビゲータ信号との間に収集されたイメージングデータを補正する、請求項２又は３に記載の磁気共鳴装置。
【請求項６】
前記ｉ番目のナビゲータ信号は、１番目のナビゲータ信号である、請求項４又は５に記載の磁気共鳴装置。
【請求項７】
前記所定の部位の信号強度の時間変化を表すプロファイルを推定する推定手段を有する、請求項１〜６のうちのいずれか一項に記載の磁気共鳴装置。
【請求項８】
前記補正手段は、
前記ナビゲータ信号に基づいて前記第１の領域の信号強度を算出する信号強度算出手段を有し、
前記信号強度に基づいて前記イメージングデータの信号強度を補正する、請求項１に記載の磁気共鳴装置。
【請求項９】
前記スキャン手段は、
前記ナビゲータ信号と前記イメージングデータとを交互に収集し、
前記補正手段は、
ｉ番目のナビゲータ信号により得られた信号強度と、ｎ番目のナビゲータ信号により得られた信号強度とに基づいて、前記ｎ番目のナビゲータ信号の直前又は直後に収集されたイメージングデータを補正する、請求項８に記載の磁気共鳴装置。
【請求項１０】
前記スキャン手段は、
前記ナビゲータ信号と複数のイメージングデータとを交互に実行し、
前記補正手段は、
ｉ番目のナビゲータ信号により得られた信号強度と、ｎ番目のナビゲータ信号により得られた信号強度とに基づいて、前記ｎ番目のナビゲータ信号とｎ＋１番目のナビゲータ信号との間に収集されたイメージングデータ、又は前記ｎ番目のナビゲータ信号とｎ−１番目のナビゲータ信号との間に収集されたイメージングデータを補正する、請求項８に記載の磁気共鳴装置。
【請求項１１】
前記ｉ番目のナビゲータ信号は、１番目のナビゲータ信号である、請求項９又は１０に記載の磁気共鳴装置。
【請求項１２】
前記第１の領域は、前記第２の領域と同じである、請求項８〜１１のうちのいずれか一項に記載の磁気共鳴装置。
【請求項１３】
前記所定の部位は、体動する部位であり、
前記スキャン手段は、前記被検体が息止めをしている間に、前記ナビゲータ信号と前記イメージングデータとを収集するためのスキャンを実行する、請求項１〜１２のうちのいずれか一項に記載の磁気共鳴装置。
【請求項１４】
前記所定の部位は、体動する部位であり、
前記スキャン手段は、呼吸信号を取得するためのプリスキャンと、前記ナビゲータ信号と前記イメージングデータとを収集するための本スキャンとを実行し、
前記ナビゲータ信号は、前記体動する部位のエッジの位置と、前記体動する部位の信号強度とを求めるために使用される、請求項１〜３のうちのいずれか一項に記載の磁気共鳴装置。
【請求項１５】
前記呼吸信号に基づいて、前記イメージングデータを画像再構成用のデータとして採用するか否かを判断するときに使用される呼吸の範囲を設定し、
本スキャンでは、前記ナビゲータ信号と前記イメージングデータとを交互に収集し、前記ナビゲータ信号により求められた前記体動する部位のエッジの位置が、前記呼吸の範囲に含まれている場合、前記ナビゲータ信号の直前又は直後に収集されたイメージングデータを、画像再構成用のデータとして採用する、請求項１４に記載の磁気共鳴装置。
【請求項１６】
前記呼吸信号に基づいて前記体動する部位のエッジの基準位置を設定し、
前記本スキャンでは、前記体動する部位のエッジの位置が前記基準位置に到達したときに、前記イメージングデータを収集する、請求項１４に記載の磁気共鳴装置。
【請求項１７】
前記ナビゲータ信号および前記イメージングデータは別々のシーケンスで収集される、請求項１〜１６のうちのいずれか一項に記載の磁気共鳴装置。
【請求項１８】
前記ナビゲータ信号および前記イメージングデータは同じシーケンスで収集される、請求項１〜１６のうちのいずれか一項に記載の磁気共鳴装置。
【請求項１９】
被検体の所定の部位を含む第１の領域からナビゲータ信号を収集するとともに、前記所定の部位を含む第２の領域から前記所定の部位の画像を再構成するときに使用されるイメージングデータを収集する磁気共鳴装置であって、
前記ナビゲータ信号と前記イメージングデータとを収集するためのスキャンを実行するスキャン手段と、
前記ナビゲータ信号に基づいて、前記イメージングデータを収集するときに使用されるシーケンスのＲＦパルスのフリップ角を補正する補正手段と、
を有する、磁気共鳴装置。
【請求項２０】
前記補正手段は、
前記ナビゲータ信号に基づいて前記第１の領域の信号強度データを作成する信号強度データ作成手段と、
前記信号強度データの中から前記所定の部位の信号強度を検出する検出手段とを有し、
前記信号強度に基づいて前記フリップ角を補正する、請求項１９に記載の磁気共鳴装置。
【請求項２１】
前記検出手段は、
前記信号強度データに基づいて前記所定の部位のエッジの位置を検出し、前記所定の部位のエッジの位置から所定の距離離れた位置の信号強度を、前記所定の部位の信号強度として検出する、請求項２０に記載の磁気共鳴装置。
【請求項２２】
前記スキャン手段は、
前記ナビゲータ信号を収集するためのナビゲータシーケンスと、前記イメージングデータを収集するためのイメージングシーケンスとを交互に実行し、
前記補正手段は、
ｉ番目のナビゲータ信号により得られた信号強度と、ｎ番目のナビゲータ信号により得られた信号強度とに基づいて、前記ｎ番目のナビゲータ信号の直後に実行されるイメージングシーケンスのＲＦパルスのフリップ角を補正する、請求項２０又は２１に記載の磁気共鳴装置。
【請求項２３】
前記ｉ番目のナビゲータ信号は、１番目のナビゲータ信号である、請求項２２に記載の磁気共鳴装置。
【請求項２４】
前記補正手段は、
前記ナビゲータ信号に基づいて前記第１の領域の信号強度を算出する信号強度算出手段を有し、
前記信号強度に基づいて前記フリップ角を補正する、請求項１９に記載の磁気共鳴装置。
【請求項２５】
前記スキャン手段は、
前記ナビゲータ信号と前記イメージングデータとを収集するためのシーケンスを複数回実行し、
前記補正手段は、
ｎ番目のシーケンスのナビゲータ信号により得られた信号強度を補正し、補正された信号強度に基づいて、ｎ＋１番目のシーケンスのＲＦパルスのフリップ角を補正する、請求項２４に記載の磁気共鳴装置。
【請求項２６】
前記第１の領域は、前記第２の領域と同じである、請求項２４又は２５に記載の磁気共鳴装置。
【請求項２７】
前記所定の部位は、体動する部位であり、
前記スキャン手段は、前記被検体が息止めをしている間に、前記ナビゲータ信号と前記イメージングデータとを収集する、請求項１９〜２６のうちのいずれか一項に記載の磁気共鳴装置。
【請求項２８】
前記所定の部位は、体動する部位であり、
前記スキャン手段は、呼吸信号を取得するためのプリスキャンと、前記ナビゲータ信号と前記イメージングデータとを収集するための本スキャンとを実行し、
前記ナビゲータ信号は、前記体動する部位のエッジの位置と、前記体動する部位の信号強度とを求めるために使用される、請求項１９〜２１のうちのいずれか一項に記載の磁気共鳴装置。
【請求項２９】
前記呼吸信号に基づいて、前記イメージングデータを画像再構成用のデータとして採用するか否かを判断するときに使用される呼吸の範囲を設定し、
本スキャンでは、前記ナビゲータ信号と前記イメージングデータとを交互に収集し、前記ナビゲータ信号により求められた前記体動する部位のエッジの位置が、前記呼吸の範囲に含まれている場合、前記ナビゲータ信号の直前又は直後に収集されたイメージングデータを、画像再構成用のデータとして採用する、請求項２８に記載の磁気共鳴装置。
【請求項３０】
前記呼吸信号に基づいて前記体動する部位のエッジの基準位置を設定し、
前記本スキャンでは、前記体動する部位のエッジの位置が前記基準位置に到達したときに、前記イメージングデータを収集する、請求項２８に記載の磁気共鳴装置。
【請求項３１】
前記被検体に造影剤を注入し、前記ナビゲータ信号および前記イメージングデータとを収集する、請求項１〜３０のうちのいずれか一項に記載の磁気共鳴装置。
【請求項３２】
被検体の所定の部位を含む第１の領域からナビゲータ信号を収集するとともに、前記所定の部位を含む第２の領域から前記所定の部位の画像を再構成するときに使用されるイメージングデータを収集するためのスキャンを実行する磁気共鳴装置のプログラムであって、
前記ナビゲータ信号に基づいて、前記イメージングデータの信号強度を補正する補正処理と、
信号強度が補正されたイメージングデータに基づいて画像再構成を行う再構成処理と、
を計算機に実行させるためのプログラム。
【請求項３３】
被検体の所定の部位を含む第１の領域からナビゲータ信号を収集するとともに、前記所定の部位を含む第２の領域から前記所定の部位の画像を再構成するときに使用されるイメージングデータを収集するためのスキャンを実行する磁気共鳴装置のプログラムであって、
前記ナビゲータ信号に基づいて、前記イメージングデータを収集するときに使用されるシーケンスのＲＦパルスのフリップ角を補正する補正処理、
を計算機に実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、所定の部位を撮影する磁気共鳴装置、およびその磁気共鳴装置に適用されるプログラムに関する。
【背景技術】
【０００２】
磁気共鳴装置を用いて被検体の診断を行う方法として、被検体に造影剤を投与し、造影剤を投与した後に現れる各時相の画像データを取得する方法が知られている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１１−１６７４０８号
【発明の概要】
【発明が解決しようとする課題】
【０００４】
近年、肝臓のＭＲＩ検査ではガドキセト酸ナトリウム（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α｝
【０１３６】
したがって、補正後のフリップ角α４_ｃｏｒｒを求めることができる。補正後のフリップ角α４_ｃｏｒｒを求めたら、シーケンスＣ４のＲＦパルスＸ４のフリップ角αをα４_ｃｏｒｒに補正し、シーケンスＣ４を実行する。したがって、ナビゲータ信号Ａ４とイメージングデータＳ４とを得ることができる。信号強度算出手段９１１ｄは、撮影領域Ｒから得られたナビゲータ信号Ａ４の積分値を計算し、撮影領域Ｒの信号強度Ｎ４を算出する。
【０１３７】
以下同様に、シーケンスＣｎを実行するたびに、式（７）を用いて信号強度Ｎｎを補正してから、式（８）を用いて補正後のフリップ角αｎ＋１_ｃｏｒｒを求め、シーケンスを実行する。最後のシーケンスＣｚが終了したら、画像再構成を行う。これにより、撮影領域Ｒの画像データを得ることができる。
【０１３８】
第６の形態では、一つのシーケンスＣｎで、ナビゲータ信号ＡｎとイメージングデータＳｎとの両方を収集している。この場合は、式（７）を用いて信号強度を補正した後で、フリップ角を補正することにより、息止めスキャンＢＨにより得られるイメージングデータの信号強度の差を小さくすることができるので、リンギングやブラーなどのアーチファクトを低減することができる。
【０１３９】
第５および第６の形態では息止めスキャンの場合について説明したが、第３および第４の形態と同様に、被検体に自由に呼吸をさせた状態でスキャンを実行し、フリップ角を補正してもよい。
【０１４０】
尚、第１〜第６の形態では、造影剤を用いてスキャンを行う場合について説明されているが、本発明は、イメージングデータの信号強度が時間とともに変化する場合であれば、造影剤を用いないスキャンを実行する場合にも適用できる。
【０１４１】
更に、第１〜第６の形態では、体動する部位の画像再構成を行う場合について説明されているが、本発明は、体動しない部位の画像再構成を行う場合にも適用することができる。
【符号の説明】
【０１４２】
２  マグネット
３  テーブル
３ａ  クレードル
４  受信コイル
５  造影剤注入装置
６  送信器
７  勾配磁場電源
８  受信器
９  制御部
１０  操作部
１１  表示部
１２  被検体
２１  ボア
２２  超伝導コイル
２３  勾配コイル
２４  ＲＦコイル
９１ａ、９１１ａ  信号強度データ作成手段
９１ｂ、９１１ｂ  検出手段
９１ｃ、９１１ｃ  演算手段
９２  再構成手段
１００、２００、３００、４００  ＭＲ装置
【図１】
2014023703.tif 000003
【図２】
2014023703.tif 000004
【図３】
2014023703.tif 000005
【図４】
2014023703.tif 000006
【図５】
2014023703.tif 000007
【図６】
2014023703.tif 000008
【図７】
2014023703.tif 000009
【図８】
2014023703.tif 000010
【図９】
2014023703.tif 000011
【図１０】
2014023703.tif 000012
【図１１】
2014023703.tif 000013
【図１２】
2014023703.tif 000014
【図１３】
2014023703.tif 000015
【図１４】
2014023703.tif 000016
【図１５】
2014023703.tif 000017
【図１６】
2014023703.tif 000018
【図１７】
2014023703.tif 000019
【図１８】
2014023703.tif 000020
【図１９】
2014023703.tif 000021
【図２０】
2014023703.tif 000022
【図２１】
2014023703.tif 000023
【図２２】
2014023703.tif 000024
【図２３】
2014023703.tif 000025
【図２４】
2014023703.tif 000026
【図２５】
2014023703.tif 000027
【図２６】
2014023703.tif 000028
【図２７】
2014023703.tif 000029
【図２８】
2014023703.tif 000030
【図２９】
2014023703.tif 000031
【図３０】
2014023703.tif 000032
【図３１】
2014023703.tif 000033
【図３２】
2014023703.tif 000034
【図３３】
2014023703.tif 000035
【図３４】
2014023703.tif 000036
【図３５】
2014023703.tif 000037
【図３６】
2014023703.tif 000038
【図３７】
2014023703.tif 000039
【図３８】
2014023703.tif 000040
【図３９】
2014023703.tif 000041
【図４０】
2014023703.tif 000042
【図４１】
2014023703.tif 000043
【図４２】
2014023703.tif 000044
【図４３】
2014023703.tif 000045
【図４４】
2014023703.tif 000046