(57)【要約】      （修正有）【課題】ＲＰＧ法に適した歪み量マップを生成する技術を提供する。【解決手段】ＭＲＩ装置は、シーケンスＶＰ１を実行することにより得られたデータに基づいて、撮影部位の画像Ｐ１を生成するとともに、第２のシーケンスＶＮ１を実行することにより得られたデータに基づいて、撮影部位の画像Ｑ１を生成する画像生成手段と、画像Ｐ１と画像Ｑ１との位置合わせを行う位置合わせ手段と、位置合わせされた後の画像Ｐ１と画像Ｑ１とに基づいて、画像Ｐ１および画像Ｑ１の歪み量を表す歪み量マップＤ１を生成する歪み量マップ生成手段とを有する。【選択図】図８
2019050937.tif 000002
【特許請求の範囲】
【請求項１】
第１の極性を有する第１の複数の位相エンコード勾配パルスを含みシングルショットで撮影部位の画像を得るための第１のシーケンスと、前記第１の極性とは反対の第２の極性を有する第２の複数の位相エンコード勾配パルスを含みシングルショットで前記撮影部位の画像を得るための第２のシーケンスとを実行するＭＲＩ装置で使用される画像処理装置であって、
前記第１のシーケンスを実行することにより得られたデータに基づいて、前記撮影部位の第１の画像を生成するとともに、前記第２のシーケンスを実行することにより得られたデータに基づいて、前記撮影部位の第２の画像を生成する画像生成手段と、
前記第１の画像と第２の画像との位置合わせを行う位置合わせ手段と、
位置合わせされた後の前記第１の画像と前記第２の画像に基づいて、前記第１の画像および前記第２の画像の歪み量を補正するための歪み量マップを生成する歪み量マップ生成手段と、
を有する画像処理装置。
【請求項２】
前記位置合わせ手段は、前記第１の画像と第２の画像との位置合わせを行うための基準画像として、前記第１の画像を用いる、請求項１に記載の画像処理装置。
【請求項３】
前記歪み量マップに基づいて、前記第１の画像又は前記第２の画像の歪みを補正する補正手段を有する請求項２に記載の画像処理装置。
【請求項４】
前記ＭＲＩ装置は、ｂ値が異なる複数の第１のシーケンスを実行し、
前記画像生成手段は、第１のシーケンスごとに、前記撮影部位の第１の画像を生成し、
前記位置合わせ手段は、前記複数の第１の画像のうちの一つの第１の画像を前記基準画像として、前記一つの第１の画像と、他の第１の画像との位置合わせを行い、
前記補正手段は、前記歪み量マップに基づいて、前記複数の第１の画像のうちの少なくとも一つの画像の歪みを補正する、請求項３に記載の画像処理装置。
【請求項５】
前記位置合わせ手段は、前記複数の第１のシーケンスのｂ値の違いに起因した歪み量の違いが低減されるように、前記一つの第１の画像と、前記他の第１の画像との位置合わせを行う、請求項４に記載の画像処理装置。
【請求項６】
前記ＭＲＩ装置は、ｂ値が異なる複数の第２のシーケンスを実行し、
前記画像生成手段は、第２のシーケンスごとに、前記撮影部位の第２の画像を生成し、
前記位置合わせ手段により位置合わせされる第１の画像および第２の画像は、同じｂ値を有する第１のシーケンスおよび第２のシーケンスにより得られた画像である、請求項１〜５のうちのいずれか一項に記載の画像処理装置。
【請求項７】
前記第１の画像を得るために実行される第１のシーケンスのｂ値は、前記第２の画像を得るために実行される第２のシーケンスのｂ値と同じである、請求項１に記載の画像処理装置。
【請求項８】
第１の極性を有する第１の複数の位相エンコード勾配パルスを含みマルチショットで撮影部位の画像を得るための複数の第１のシーケンスを含む第１のシーケンス群と、前記第１の極性とは反対の第２の極性を有する第２の複数の位相エンコード勾配パルスを含みマルチショットで前記撮影部位の画像を得るための複数の第２のシーケンスを含む第２のシーケンス群とを実行するＭＲＩ装置で使用される画像処理装置であって、
前記第１のシーケンス群を実行することにより得られたデータに基づいて、前記撮影部位の第１の画像を生成するとともに、前記第２のシーケンス群を実行することにより得られたデータに基づいて、前記撮影部位の第２の画像を生成する画像生成手段と、
前記第１の画像と第２の画像との位置合わせを行う位置合わせ手段と、
位置合わせされた後の前記第１の画像と前記第２の画像に基づいて、前記第１の画像および前記第２の画像の歪み量を補正するための歪み量マップを生成する歪み量マップ生成手段と、
を有する画像処理装置。
【請求項９】
請求項１〜８のうちのいずれか一項に記載の画像処理装置を有するＭＲＩ装置。
【請求項１０】
第１の極性を有する第１の複数の位相エンコード勾配パルスを含みシングルショットで撮影部位の画像を得るための第１のシーケンスと、前記第１の極性とは反対の第２の極性を有する第２の複数の位相エンコード勾配パルスを含みシングルショットで前記撮影部位の画像を得るための第２のシーケンスとを実行するＭＲＩ装置で使用される画像処理装置に適用されるプログラムであって、
前記第１のシーケンスを実行することにより得られたデータに基づいて、前記撮影部位の第１の画像を生成するとともに、前記第２のシーケンスを実行することにより得られたデータに基づいて、前記撮影部位の第２の画像を生成する画像生成処理と、
前記第１の画像と第２の画像との位置合わせを行う位置合わせ処理と、
位置合わせされた後の前記第１の画像と前記第２の画像に基づいて、前記第１の画像および前記第２の画像の歪み量を補正するための歪み量マップを生成する歪み量マップ生成処理と、
をコンピュータに実行させるためのプログラム。
【請求項１１】
第１の極性を有する第１の複数の位相エンコード勾配パルスを含みマルチショットで撮影部位の画像を得るための複数の第１のシーケンスを含む第１のシーケンス群と、前記第１の極性とは反対の第２の極性を有する第２の複数の位相エンコード勾配パルスを含みマルチショットで前記撮影部位の画像を得るための複数の第２のシーケンスを含む第２のシーケンス群とを実行するＭＲＩ装置に適用されるプログラムであって、
前記第１のシーケンス群を実行することにより得られたデータに基づいて、前記撮影部位の第１の画像を生成するとともに、前記第２のシーケンス群を実行することにより得られたデータに基づいて、前記撮影部位の第２の画像を生成する画像生成処理と、
前記第１の画像と第２の画像との位置合わせを行う位置合わせ処理と、
位置合わせされた後の前記第１の画像と前記第２の画像に基づいて、前記第１の画像および前記第２の画像の歪み量を補正するための歪み量マップを生成する歪み量マップ生成処理と、
をコンピュータに実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、歪み量マップを生成するための処理を実行する画像処理装置、および当該画像処理装置を有する磁気共鳴イメージング装置に関する。
【背景技術】
【０００２】
磁気共鳴イメージング装置（ＭＲＩ：magnetic resonance imaging）を用いて撮像する方法として、ＥＰＩ（Echo Planar Imaging）法を用いた撮像方法が知られている。ＥＰＩは、高速で被検体の画像を生成することができるので、スキャン時間の短縮に有効な手法である。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１５−１１２４７４号公報
【非特許文献】
【０００４】
【非特許文献１】Holland et al. Neuroimage. 2010 March ; 50(1): 175Morgan PS et al., JMRI 2004 19; 499-507
【発明の概要】
【発明が解決しようとする課題】
【０００５】
ＥＰＩは、１回の励起で複数のエコーを収集することなどが原因で画像が歪むことが知られている。そこで、画像の歪みを補正するために、歪みマップを生成し、歪みマップに基づいて画像の歪みを補正する方法が開示されている（特許文献１参照）。
【０００６】
また、ＥＰＩ法による画像歪みを補正する方法として、位相エンコード勾配パルスの極性を反転させる方法が知られている（非特許文献１参照）。この方法は、正極性の位相エンコード勾配パルスを用いたシーケンスと、負極性の位相エンコード勾配パルスを用いたシーケンスとを実行する。そして、正極性の画像と負極性の画像を用いて、歪み量マップを生成し、この歪み量マップを用いて画像歪みを補正している。
【０００７】
しかし、位相エンコード勾配パルスの極性を反転させる当該方法（以下、「ＲＰＧ法」と呼ぶ。ＲＰＧ：Reversed Polarity Gradient）は、静磁場不均一、スキャン中の被検体の体動、中心周波数のずれなどにより、歪み量が過大評価された歪みマップが生成されることがある。この場合、画像の歪み補正がうまくいかないという問題がある。
【０００８】
したがって、ＲＰＧ法において画像の歪み量が適切に評価された歪み量マップを生成する技術が望まれている。
【課題を解決するための手段】
【０００９】
本発明の第１の観点は、第１の極性を有する第１の複数の位相エンコード勾配パルスを含みシングルショットで撮影部位の画像を得るための第１のシーケンスと、前記第１の極性とは反対の第２の極性を有する第２の複数の位相エンコード勾配パルスを含みシングルショットで前記撮影部位の画像を得るための第２のシーケンスとを実行するＭＲＩ装置で使用される画像処理装置であって、
前記第１のシーケンスを実行することにより得られたデータに基づいて、前記撮影部位の第１の画像を生成するとともに、前記第２のシーケンスを実行することにより得られたデータに基づいて、前記撮影部位の第２の画像を生成する画像生成手段と、
前記第１の画像と第２の画像との位置合わせを行う位置合わせ手段と、
位置合わせされた後の前記第１の画像と前記第２の画像に基づいて、前記第１の画像および前記第２の画像の歪み量を補正するための歪み量マップを生成する歪み量マップ生成手段と、
を有する画像処理装置である。
【００１０】
本発明の第２の観点は、第１の極性を有する第１の複数の位相エンコード勾配パルスを含みマルチショットで撮影部位の画像を得るための複数の第１のシーケンスを含む第１のシーケンス群と、前記第１の極性とは反対の第２の極性を有する第２の複数の位相エンコード勾配パルスを含みマルチショットで前記撮影部位の画像を得るための複数の第２のシーケンスを含む第２のシーケンス群とを実行するＭＲＩ装置で使用される画像処理装置であって、
前記第１のシーケンス群を実行することにより得られたデータに基づいて、前記撮影部位の第１の画像を生成するとともに、前記第２のシーケンス群を実行することにより得られたデータに基づいて、前記撮影部位の第２の画像を生成する画像生成手段と、
前記第１の画像と第２の画像との位置合わせを行う位置合わせ手段と、
位置合わせされた後の前記第１の画像と前記第２の画像に基づいて、前記第１の画像および前記第２の画像の歪み量を補正するための歪み量マップを生成する歪み量マップ生成手段と、
を有する画像処理装置である。
【００１１】
本発明の第３の観点は、第１又は第２の観点の画像処理装置を有するＭＲＩ装置である。
【００１２】
本発明の第４の観点は、第１の極性を有する第１の複数の位相エンコード勾配パルスを含みシングルショットで撮影部位の画像を得るための第１のシーケンスと、前記第１の極性とは反対の第２の極性を有する第２の複数の位相エンコード勾配パルスを含みシングルショットで前記撮影部位の画像を得るための第２のシーケンスとを実行するＭＲＩ装置で使用される画像処理装置に適用されるプログラムであって、
前記第１のシーケンスを実行することにより得られたデータに基づいて、前記撮影部位の第１の画像を生成するとともに、前記第２のシーケンスを実行することにより得られたデータに基づいて、前記撮影部位の第２の画像を生成する画像生成処理と、
前記第１の画像と第２の画像との位置合わせを行う位置合わせ処理と、
位置合わせされた後の前記第１の画像と前記第２の画像に基づいて、前記第１の画像および前記第２の画像の歪み量を補正するための歪み量マップを生成する歪み量マップ生成処理と、
をコンピュータに実行させるためのプログラムである。
【００１３】
本発明の第５の観点は、第１の極性を有する第１の複数の位相エンコード勾配パルスを含みマルチショットで撮影部位の画像を得るための複数の第１のシーケンスを含む第１のシーケンス群と、前記第１の極性とは反対の第２の極性を有する第２の複数の位相エンコード勾配パルスを含みマルチショットで前記撮影部位の画像を得るための複数の第２のシーケンスを含む第２のシーケンス群とを実行するＭＲＩ装置に適用されるプログラムであって、
前記第１のシーケンス群を実行することにより得られたデータに基づいて、前記撮影部位の第１の画像を生成するとともに、前記第２のシーケンス群を実行することにより得られたデータに基づいて、前記撮影部位の第２の画像を生成する画像生成処理と、
前記第１の画像と第２の画像との位置合わせを行う位置合わせ処理と、
位置合わせされた後の前記第１の画像と前記第２の画像に基づいて、前記第１の画像および前記第２の画像の歪み量を補正するための歪み量マップを生成する歪み量マップ生成処理と、
をコンピュータに実行させるためのプログラムである。
【発明の効果】
【００１４】
画像の位置合わせを行った後で歪み量マップを生成するので、歪み量の過大評価を回避することができる。
【図面の簡単な説明】
【００１５】
【図１】本発明の第１の形態の磁気共鳴イメージング装置の概略図である。
【図２】制御部５の機能ブロック図の説明図である。
【図３】撮影部位を概略的に示す図である。
【図４】ＲＰＧ法を用いて撮影部位の画像を得るために実行されるスキャンＳＣの説明図である。
【図５】２個のシーケンスＶＰ１およびＶＰ２並びに２個のシーケンスＶＮ１およびＶＮ２が実行されるスキャンを示す図である。
【図６】ＭＲＩ装置の動作フローの一例を示す図である。
【図７】シーケンスの各々を実行することにより得られたスライスＳＬの画像Ｐ１およびＰ２並びに画像Ｑ1およびＱ２を概略的に示す図である。
【図８】画像Ｐ１およびＱ１の歪みを補正するために実行される処理を示す図である。
【図９】画像Ｐ２およびＱ２の歪みを補正するために実行される処理を示す図である。
【図１０】ｍ≠ｎの場合に実行されるスキャンＳＣの説明図である。
【図１１】第２の形態におけるＭＲＩ装置の動作フローの説明図である。
【図１２】シーケンスの各々を実行することにより得られたスライスＳＬの画像Ｐ１、Ｐ２、およびＱ１を概略的に示す図である。
【図１３】画像Ｐ１およびＱ１の歪みを補正するために実行される処理を示す図である。
【図１４】画像Ｐ２の歪みを補正するために実行される処理を示す図である。
【図１５】マルチショットで画像を取得するためのスキャンＳＣの一例を示す図である。
【図１６】第３の形態における歪み補正の説明図である。
【図１７】位置合わせが行われないファントム画像に基づいて歪み量マップを生成し、当該歪み量マップを用いて歪み補正されたファントム画像を示す図である。
【図１８】位置合わせが行われたファントム画像に基づいて歪み量マップを生成し、当該歪み量マップを用いて歪み補正されたファントム画像を示す図である。
【発明を実施するための形態】
【００１６】
以下、発明を実施するための形態について説明するが、本発明は、以下の形態に限定されることはない。
【００１７】
図１は、本発明の第１の形態の磁気共鳴イメージング装置の概略図である。
【００１８】
磁気共鳴イメージング装置（以下、「ＭＲＩ装置」と呼ぶ。ＭＲＩ：Magnetic Resonance Imaging）１は、マグネット２、テーブル３、受信コイル４などを有している。
【００１９】
マグネット２は、被検体１２が収容される収容空間２１を有している。また、マグネット２は、超伝導コイル２２と、勾配コイル２３と、ＲＦコイル２４とを有している。超伝導コイル２２は静磁場を印加し、勾配コイル２３は勾配パルスを印加し、ＲＦコイル２４はＲＦパルスを印加する。尚、超伝導コイル２２の代わりに、永久磁石を用いてもよい。
【００２０】
テーブル３は、被検体１２を搬送するためのクレードル３ａを有している。クレードル３ａによって、被検体１２は収容空間２１に搬送される。
【００２１】
受信コイル４は、被検体１２に取り付けられている。受信コイル４は、被検体１２からの磁気共鳴信号を受信する。
【００２２】
ＭＲＩ装置１は、更に、制御部５、送信器６、勾配磁場電源７、受信器８、記憶部９、操作部１０、および表示部１１などを有している。
【００２３】
制御部５は、ＲＦパルスのデータに基づいて送信器６を制御し、勾配パルスのデータに基づいて勾配磁場電源７を制御する。また、制御部５は、クレードル３ａの移動の制御なども行う。尚、図１では、制御部５が、送信器６、勾配磁場電源７、クレードル３ａなどの制御を行っているが、送信器６、勾配磁場電源７、クレードル３ａなどの制御を複数の制御部で行ってもよい。例えば、送信器６および勾配磁場電源７を制御する制御部と、クレードル３ａを制御する制御部とを別々に設けてもよい。
【００２４】
送信器６は、制御部５から受け取ったデータに基づいて、ＲＦコイル２４に電流を供給する。
【００２５】
勾配磁場電源７は、制御部５から受け取ったデータに基づいて、勾配コイル２３に電流を供給する。
【００２６】
受信器８は、受信コイル４で受信された磁気共鳴信号に対して、検波などの処理を行い制御部５に出力する。
【００２７】
記憶部９には、制御部５により実行されるプログラムなどが記憶されている。尚、記憶部９は、ハードディスク、ＣＤ−ＲＯＭなどの非一過性の記憶媒体であってもよい。制御部５は、記憶部９に記憶されているプログラムを読み出し、プログラムに記述されている処理を実行するプロセッサとして動作する。制御部５は、プログラムに記述されている処理を実行することにより、種々の手段を実現する。図２は、制御部５の機能ブロック図の説明図である。
【００２８】
画像生成手段１０１は、被検体１２の撮影部位の画像を生成する。
【００２９】
画像選択手段１０２は、歪み量マップの生成に使用される画像を選択する。
【００３０】
位置合わせ手段１０３は、画像の位置あわせを行う。
【００３１】
歪み量マップ生成手段１０４は、位置合わせされた後の画像に基づいて、画像の歪み量を表す歪み量マップを生成する。
【００３２】
補正手段１０５は、歪み量マップを用いて、位置合わせされた後の画像の歪み補正を行う。
【００３３】
これらの手段１０１〜１０５の詳細については後述する。
【００３４】
制御部５は、例えば、コンピュータの少なくとも一部として構成することができる。制御部５は、記憶部９に記憶されているプログラムを読み出すことにより、画像生成手段１０１〜補正手段１０５などを実現する。尚、制御部５は、一つのプロセッサで画像生成手段１０１〜補正手段１０５を実現してもよいし、２つ以上のプロセッサで、画像生成手段１０１〜補正手段１０５を実現してもよい。また、制御部５が実行するプログラムは、一つの記憶部に記憶させておいてもよいし、複数の記憶部に分けて記憶させておいてもよい。制御部５は画像処理装置の一例に相当する。
【００３５】
図１に戻って説明を続ける。
【００３６】
操作部１０は、オペレータにより操作され、種々の情報を入力する。表示部１１は種々の情報を表示する。
【００３７】
ＭＲＩ装置１は、上記のように構成されている。
【００３８】
図３は撮影部位を概略的に示す図である。
【００３９】
第１の形態では、撮影部位は頭部であり、頭部を横切るように、スライスＳＬが設定される。尚、第１の形態では、説明の便宜上、スライスＳＬが１枚のみ設定された例について説明されているが、撮影部位に複数枚のスライスを設定してもよい。また、撮影部位は頭部に限定されることはなく、頭部とは別の部位（例えば、腹部、脚部）を撮影部位としてもよい。
【００４０】
図４は、ＲＰＧ法を用いて撮影部位の画像を得るために実行されるスキャンＳＣの説明図である。
【００４１】
第１の形態では、スキャンＳＣにおいて、正極性の位相エンコード勾配パルスを有するｍ個のシーケンスＶＰ１〜ＶＰｍと、負極性の位相エンコード勾配パルスを有するｎ個のシーケンスＶＮ１〜ＶＮｎが実行される。
【００４２】
先ず、ｍ個のシーケンスＶＰ１〜ＶＰｍについて説明する。
【００４３】
図４には、シーケンスＶＰ１の具体的な一例が示されている。シーケンスＶＰ１は、ＥＰＩ法を用いてシングルショットにより撮影部位の画像を得るためのシーケンスである。シーケンスＶＰ１は９０°パルスおよび１８０°パルスを有している。１８０°パルスは９０°パルスから所定の時間間隔をあけて印加される。
【００４４】
また、シーケンスＶＰ１は、スライス選択方向ＳＳに印加されるスライス選択勾配パルスＧｓ１およびＧｓ２を有している。スライス選択勾配パルスＧｓ１およびＧｓ２は、それぞれ、９０°パルスおよび１８０°パルスの印加時に印加される。また、シーケンスＶＰ１は、スライス選択方向ＳＳに印加されるＭＰＧ（Motion Probing Gradient）を有している。
【００４５】
また、シーケンスＶＰ１は、読出し方向ＲＥに印加されるリードアウト用の勾配パルスＧｒを有している。リーとアウト用の勾配パルスＧｒは、正極性のリードアウト用の勾配パルスと負極性のリードアウト用の勾配パルスとを有している。正極性のリードアウト用の勾配パルスと負極性のリードアウト用の勾配パルスとは、交互に印加される。
【００４６】
更に、シーケンスＶＰ１は、位相エンコード方向ＰＥに印加される補正勾配パルスＧｐ０と、位相エンコード用勾配パルスＧｐ１を有している。補正勾配パルスＧｐ０は負極性であり、位相エンコード用勾配パルスＧｐ１は正極性である。
【００４７】
シーケンスＶＰ１は、上記のようなパルスを有している。
【００４８】
他の（ｍ−１）個のシーケンスＶＰ２〜ＶＰｍも、シーケンスＶＰ１と同じシーケンスである。しかし、ｍ個のシーケンスＶＰ１〜ＶＰｍは、ＭＰＧのｂ値が異なる値を有するように設計されている。図４では、シーケンスＶＰ１、ＶＰ２、・・・ＶＰｍのｂ値は、それぞれ、ｂ＝ｂ１、ｂ２、・・・ｂｍで表されている。尚、ｍ個のｂ値（ｂ１〜ｂｍ）には、ｂ＝０の値が含まれていてもよい。ＭＰＧの磁場強度をゼロに設定する（即ち、ＭＰＧを印加しない）ことにより、ｂ＝０に設定することができる。
【００４９】
次に、ｎ個のシーケンスＶＮ１〜ＶＮｎについて説明する。
【００５０】
図４には、シーケンスＶＮ１の具体的な一例が示されている。シーケンスＶＮ１は、補正勾配パルスＧｐ０が正極性を有し、位相エンコード勾配パルスＧｐ１が負極性を有している点を除いて、シーケンスＶＰ１と同じである。したがって、シーケンスＶＮ１は、シーケンスＶＰ１と同じｂ値（ｂ＝ｂ１）を有している。
【００５１】
他の（ｎ−１）個のシーケンスＶＮ２〜ＶＮｎも、シーケンスＶＮ１と同じシーケンスである。しかし、ｎ個のシーケンスＶＮ１〜ＶＮｎは、ＭＰＧのｂ値が異なる値を有するように設計されている。図４では、シーケンスＶＮ１、ＶＮ２、・・・ＶＮｎのｂ値は、それぞれ、ｂ＝ｂ１、ｂ２・・・ｂｎで表されている。尚、ｎ個のｂ値（ｂ１〜ｂｎ）には、ｂ＝０の値が含まれていてもよい。ＭＰＧの磁場強度をゼロに設定する（即ち、ＭＰＧを印加しない）ことにより、ｂ＝０に設定することができる。
【００５２】
第１の形態では、シーケンスＶＰ１〜ＶＰｍおよびシーケンスＶＮ１〜ＶＮｎを実行し、シーケンスごとにスライスＳＬの画像を生成する。そして、後述する歪み量マップを生成し、歪み量マップに基づいて、各画像の歪みを補正し、歪みが軽減された画像を生成する。以下に、シーケンスＶＰ１〜ＶＰｍおよびシーケンスＶＮ１〜ＶＮｎを実行し、歪みが軽減された画像を取得するためのＭＲＩ装置の動作フローの一例について説明する。
【００５３】
尚、以下では、第１の形態を理解しやすくするため、ｍ＝ｎ＝２、即ち、２個のシーケンスＶＰ１およびＶＰ２並びに２個のシーケンスＶＮ１およびＶＮ２を実行する例（図５参照）を取り上げて、第１の形態について説明する。
【００５４】
図６は、ＭＲＩ装置の動作フローの一例を示す図である。
【００５５】
ステップＳＴ１では、スキャンＳＣ（図５参照）が実行される。
【００５６】
スキャンＳＣを実行する場合、制御部５（図１参照）は、スキャンＳＣで使用される各シーケンスのＲＦパルスのデータを送信器６に送り、ローカライザスキャンＬＳで使用されるシーケンスの勾配パルスのデータを勾配磁場電源７に送る。送信器６は、制御部５から受け取ったデータに基づいてＲＦコイル２４に電流を供給し、勾配磁場電源７は、制御部５から受け取ったデータに基づいて勾配コイル２３に電流を供給する。したがって、ＲＦコイル２４はＲＦパルスを印加し、勾配コイル２３は勾配パルスを印加する。スキャンＳＣが実行されることにより、撮影部位からＭＲ信号が発生する。ＭＲ信号は受信コイル４（図１参照）で受信される。受信コイル４は、ＭＲ信号を受信し、ＭＲ信号の情報を含むアナログ信号を出力する。受信器８は、受信コイル４から受け取った信号に対して、検波などの信号処理を行い、信号処理により得られたデータを制御部５に出力する。
【００５７】
画像生成手段１０１（図２参照）は、スキャンＳＣにおいて各シーケンスを実行することにより収集されたデータに基づいて、スライスＳＬの画像を生成する。図７に、シーケンスの各々を実行することにより得られたスライスＳＬの画像Ｐ１およびＰ２並びに画像Ｑ1およびＱ２を概略的に示す。
【００５８】
スキャンＳＣを実行した後、ステップＳＴ２に進む。
【００５９】
ステップＳＴ２では、画像選択手段１０２（図２参照）が、画像の歪み補正に使用される歪み量マップを生成するために使用される画像を選択する。ここでは、ｂ＝ｂ１に対応する画像Ｐ１およびＱ１を用いて歪み量マップを生成する例について考える。したがって、画像選択手段１０２は、ｂ＝ｂ１に対応する画像Ｐ１および画像Ｑ１を選択する。画像を選択した後、ステップＳＴ３に進む。
【００６０】
ステップＳＴ３では、位置合わせ手段１０３（図２参照）が画像Ｐ１およびＱ１の位置合わせを行う。ＥＰＩ法を用いたシーケンスでは、静磁場不均一、スキャン中の被検体の体動、シーケンス実行時の中心周波数のずれなどにより、画像Ｐ１と画像Ｑ１との間の位置ずれ量が大きくなることがある。画像Ｐ１およびＱ１の位置ずれ量が大きいままの状態で後述する歪み量マップを生成すると（ステップＳＴ４参照）、画像の歪み量が過大評価されてしまい、適正な歪み補正を行うことができない。そこで、第１の形態では、ステップＳＴ３において、画像Ｐ１およびＱ１の位置合わせを行う。画像の位置合わせの方法としては、例えば、画像Ｐ１の重心と、画像Ｑ１の重心とを求め、画像Ｐ１およびＱ１の重心が一致するように、画像Ｐ１およびＱ１を、相対的に平行移動および回転させる方法を用いることができる。画像Ｐ１およびＱ１の位置合わせを行う場合、画像Ｐ１およびＱ１のうちの一方の画像を位置合わせの基準画像として選択し、基準画像として選択された一方の画像に対して、他方の画像を位置合わせすることができる。図８に、位置合わせ処理が行われた後の画像Ｐ１およびＱ１を概略的に示す。
【００６１】
位置合わせを行った後、ステップＳＴ４に進む。
【００６２】
ステップＳＴ４では、歪み量マップ生成手段１０４（図２参照）が、位置合わせされた後の画像Ｐ１およびＱ１に基づいて、画像Ｐ１およびＱ１の歪み量を表す歪み量マップを生成する。歪み量マップは、非特許文献１に記載された手法を用いることができる。図８に、歪み量マップ生成処理により生成された歪み量マップＤ１を概略的に示す。歪み量マップＤ１を生成した後、ステップＳＴ５に進む。
【００６３】
ステップＳＴ５では、補正手段１０５（図２参照）が、歪み量マップＤ１を用いて、位置合わせ後の画像Ｐ１およびＱ１の各々に対して、画像の歪みを補正する。図８に、補正処理により歪みが補正された後の画像を示す。図８では、歪みが補正された後の画像Ｐ１は符号「Ｐ１’」で示されており、歪みが補正された後の画像Ｑ１は符号「Ｑ１’」で示されている。尚、参考のため、歪みが補正される前の画像Ｐ１およびＱ１が破線で表されている。
【００６４】
補正を行ったら、ステップＳＴ６に進む。
【００６５】
ステップＳＴ６では、全ての画像の歪み補正が終了したか否かを判定する。ここでは、画像Ｐ１およびＱ１の歪み補正は終了したが、画像Ｐ２およびＱ２の歪み補正はまだ終了していない。したがって、ステップＳＴ２に戻る。
【００６６】
ステップＳＴ２では、画像選択手段１０２が画像Ｐ２およびＱ２を選択する。そして、ステップＳＴ３に進む。ステップＳＴ３は、位置合わせ手段１０３が画像Ｐ２およびＱ２の位置合わせを行う。尚、画像Ｐ２およびＱ２の位置合わせは、画像Ｐ１およびＱ１の位置合わせと同じ方法で行うことができる。図９に、位置合わせ処理後の画像Ｐ２およびＱ２を概略的に示す。
【００６７】
位置合わせを行った後、ステップＳＴ４に進む。
【００６８】
ステップＳＴ４では、歪み量マップ生成手段１０４が、位置合わせされた後の画像Ｐ２およびＱ２に基づいて、画像Ｐ２およびＱ２の歪み量を表す歪み量マップを生成する。図９に、歪み量マップ生成処理により生成された歪み量マップＤ２を概略的に示す。歪み量マップＤ２を生成した後、ステップＳＴ５に進む。
【００６９】
ステップＳＴ５では、補正手段１０５が、歪み量マップＤ２を用いて、位置合わせ後の画像Ｐ２およびＱ２の各々に対して、歪み量を補正する。図９に、補正処理により歪みが補正された後の画像を示す。図９では、歪みが補正された後の画像Ｐ２は符号「Ｐ２’」で示されており、歪み量が補正された後の画像Ｑ２は符号「Ｑ２’」で示されている。また、参考のため、歪みが補正される前の画像Ｐ２およびＱ２が破線で表されている。
【００７０】
補正を行ったら、ステップＳＴ６に進む。
【００７１】
ステップＳＴ６では、全ての画像の歪み補正が終了したか否かを判定する。ここでは、画像Ｐ１およびＱ１の歪み補正と、画像Ｐ２およびＱ２の歪み補正とが終了している。したがって、図６に示すフローが終了する。
【００７２】
第１の形態では、歪み量マップＤ１を生成する前に、画像Ｐ１およびＱ１の位置合わせを行う。したがって、静磁場不均一、スキャン中の被検体の体動、シーケンス実行時の中心周波数のずれなどにより、画像Ｐ１と画像Ｑ１との間の位置ずれ量が大きくなっても、位置合わせにより、画像Ｐ１と画像Ｑ１との間の位置ずれ量をゼロに近づけることができる。そして、位置ずれが十分に小さい画像Ｐ１およびＱ１に基づいて歪み量マップＤ１が生成される。したがって、第１の形態では、画像間の位置ずれが生じたままで歪み量マップＤ１を生成した場合に起こり得る歪み量の過大評価を回避することができるので、画像Ｐ１およびＱ１の歪み量を適切に反映した歪み量マップＤ１を得ることができる。また、この歪み量マップＤ１を用いて画像Ｐ１およびＱ１の歪みが補正されるので、歪みが低減された画像Ｐ１’およびＱ１’を得ることができる。
【００７３】
同様に、画像Ｐ２およびＱ２についても、歪み量が適切に反映された歪み量マップＤ２を用いて歪み補正が行われるので、歪みが低減された画像Ｐ２’およびＱ２’を得ることができる。
【００７４】
尚、第１の形態では、画像Ｐ１、Ｐ２、Ｑ１、およびＱ２の全ての画像の歪みを補正している。しかし、全ての画像Ｐ１、Ｐ２、Ｑ１、およびＱ２の歪みを補正する必要はなく、ｂ＝ｂ１に対応する画像Ｐ１およびＱ１のうちの一方の画像のみを歪み補正し、ｂ＝ｂ２に対応する画像Ｐ２およびＱ２のうちの一方の画像のみを歪み補正してもよい。
（２）第２の形態
第１の形態では、正極性の位相エンコード勾配パルスを有するシーケンスの個数を表すｍと、負極性の位相エンコード勾配パルスを有するシーケンスの個数を表すｎが、ｍ＝ｎ（＝２）の例について説明した。第２の形態では、ｍ≠ｎの例について説明する。
【００７５】
図１０は、ｍ≠ｎの場合に実行されるスキャンＳＣの説明図である。第２の形態では、ｍ≠ｎの例として、ｍ＝２およびｎ＝１、即ち、シーケンスＶＰ１およびＶＰ２並びにシーケンスＶＮ１が実行される例について説明する。
【００７６】
図１１は、第２の形態におけるＭＲＩ装置の動作フローの説明図である。
【００７７】
ステップＳＴ１では、スキャンＳＣ（図１０参照）が実行される。
【００７８】
画像生成手段１０１（図２参照）は、スキャンＳＣにおいて各シーケンスを実行することにより収集されたデータに基づいて、スライスＳＬの画像を生成する。図１２に、シーケンスの各々を実行することにより得られたスライスＳＬの画像Ｐ１、Ｐ２、およびＱ１を概略的に示す。
【００７９】
スキャンＳＣを実行した後、ステップＳＴ２に進む。
【００８０】
ステップＳＴ２では、画像選択手段１０２が、歪み量マップを生成するために使用される画像を選択する。ここでは、画像Ｐ１およびＱ１が選択されたとする。画像Ｐ１およびＱ１を選択した後、ステップＳＴ３に進む。
【００８１】
ステップＳＴ３では、位置合わせ手段１０３が、画像Ｐ１およびＱ１の位置合わせを行う。第２の形態では、位置合わせ手段１０３は、画像Ｐ１およびＱ１のうちの一方の画像を位置合わせの基準画像として選択し、基準画像として選択された一方の画像に対して、他方の画像の位置合わせを行う。ここでは、画像Ｐ１が位置合わせの基準画像として選択されたとする。したがって、位置合わせ手段１０３は、画像Ｐ１に対して画像Ｑ１の位置合わせを行う。図１３に、位置合わせ処理が実行された後の画像Ｐ１およびＱ１を概略的に示す。
【００８２】
位置合わせを行った後、ステップＳＴ４に進む。
【００８３】
ステップＳＴ４およびＳＴ５は第１の形態と同じであり、図１３に示すように、ステップＳＴ４において歪み量マップＤ１が生成され、ステップＳＴ５において画像の歪み補正が行われる。
【００８４】
補正を行ったら、ステップＳＴ６に進む。
【００８５】
ステップＳＴ６では、全ての画像の歪み補正が終了したか否かを判定する。ここでは、画像Ｐ１およびＱ１の歪み補正は終了したが、画像Ｐ２の歪み補正はまだ終了していない。したがって、ステップＳＴ７０に進む。
【００８６】
ステップＳＴ７０では、ｂ＝ｂ２に対応する画像Ｑ２の歪みが低減された画像を生成するための処理を実行する。この処理は、以下のように実行される。
【００８７】
先ず、ステップＳＴ７において、画像選択手段１０２が、歪み量マップＤ１を生成するために使用された画像Ｐ１およびＱ１のうち、位置合わせの基準画像として使用された画像を選択する。画像Ｐ１が位置合わせの基準画像として使用された場合、画像選択手段１０２は画像Ｐ１を選択し、画像Ｑ１が位置合わせの基準画像として使用された場合、画像選択手段１０２は画像Ｑ１を選択する。ここでは、画像Ｐ１が位置合わせの基準画像として使用されているので、画像選択手段１０２は、画像Ｐ１を選択する。
【００８８】
ステップＳＴ８では、位置合わせ手段１０３ｈｓ、画像Ｐ１を基準画像にして、画像Ｐ１と画像Ｐ２との位置合わせを行う。尚、画像Ｐ１と画像Ｐ２はｂ値が異なるシーケンスから生成されているので、画像Ｐ１と画像Ｐ２との間には、ＭＰＧのｂ値の違いに起因した歪み量の違いが生じていると考えられる。そこで、画像Ｐ１と画像Ｐ２とを位置合わせする場合には、歪み量の違いを低減できるような位置合わせを行うことが望ましい。このような位置合わせの手法としては、例えば、Affine変換の技術を用いることができる。図１４に、位置合わせ後の画像Ｐ１およびＰ２を概略的に示す。
【００８９】
位置合わせを行った後、ステップＳＴ９に進む。
【００９０】
ステップＳＴ９では、補正手段１０５が、歪み量マップＤ１を用いて、画像Ｐ１に対して位置合わせされた後の画像Ｐ２の歪みを補正する。図１４に、歪みが補正された後の画像を示す。図１４では、歪みが補正された後の画像Ｐ２は符号「Ｐ２’」で示されている。
【００９１】
補正を行ったら、ステップＳＴ６に戻る。
【００９２】
ステップＳＴ６では、全ての画像の歪み補正が終了したか否かを判定する。ここでは、画像Ｐ１およびＱ１の歪み補正と、画像Ｐ２の歪み補正が終了している。したがって、図１１に示すフローが終了する。
【００９３】
第２の形態でも、第１の形態と同様に、画像Ｐ１およびＱ１の位置合わせを行った後で歪み量マップを生成する。したがって、静磁場不均一、スキャン中の被検体の体動、シーケンス実行時の中心周波数のずれなどにより、画像Ｐ１と画像Ｑ１との間の位置ずれ量が大きくても、歪み量の過大評価が回避された歪み量マップを生成することができる。また、画像Ｐ１に対して画像Ｐ２を位置合わせすることにより、画像Ｐ１およびＱ１から求められた歪み量マップＤ１を用いて画像Ｐ２の歪み補正を行うことができる。したがって、ｂ＝ｂ１に対応した歪み量マップＤ１を用いて、別のｂ値（ｂ＝ｂ２）の画像Ｐ２の歪み補正をすることもできる。
【００９４】
尚、第２の形態では、ｍ＝２の場合について説明したが、本発明は、ｍ＝２に限定されることはなく、ｍ≧３の場合にも適用することができる。例えば、ｍ＝ｍ１（≧３）、即ち、ｍ１個のシーケンスＶＰ１〜ＶＰｍ１を実行することができる。この場合、シーケンスＶＰ１〜ＶＰｍ１から得られた画像Ｐ１〜Ｐｍ１のうち、画像Ｐ１を基準画像とし、画像Ｐ１に対して画像Ｐ２〜Ｐｍ１の各々の位置合わせを行い、位置合わせ後の画像Ｐ２〜Ｐｍ１の各々に対して歪み量マップＤ１を用いて歪み補正を行うことができる。したがって、歪み量マップＤ１を用いて、異なるｂ値に対応したｍ１個の画像Ｐ１〜Ｐｍ１の歪み補正を行うことができる。
【００９５】
また、第２の形態では、ｎ＝１の場合について説明したが、本発明は、ｎ＝１に限定されることはなく、ｎ≧２の場合にも適用することができる。例えば、ｎ＝ｎ１（≧２）、即ち、ｎ１個のシーケンスＶＮ１〜ＶＮｎ１を実行することができる。この場合、シーケンスＶＮ１〜ＶＮｎ１から画像Ｑ１〜Ｑｎ１を生成し、画像Ｐ１（基準画像）に対して画像Ｑ１〜Ｑｎ１の各々の位置合わせを行い、位置合わせ後の画像Ｑ１〜Ｑｎ１の各々に対して歪み量マップＤ１を用いて歪み補正を行うことができる。したがって、歪み量マップＤ１を用いて、異なるｂ値に対応したｎ１個の画像Ｑ１〜Ｑｎ１の歪み補正を行うこともできる。尚、ｍ＝ｍ１（≧３）およびｎ＝ｎ１（≧２）の場合、ｂ＝ｂ１に対応した画像Ｐ１およびＱ１を用いて歪み量マップを生成してもよいし、別のｂ値（ｂ＝ｂｉ）に対応した画像ＰｉおよびＱｉを用いて歪み量マップを生成してもよい。更に、異なるｂ値に対応した画像Ｐｉ（ｂ＝ｂｉ）および画像Ｑｊ（ｂ＝ｂｊ）であっても、ｂｉおよびｂｊが近い値であれば、画像Ｐｉおよび画像Ｑｊを用いて歪み量マップを生成することも可能である。
【００９６】
尚、第１および第２の形態では、スライス選択方向ＳＳにＭＰＧが印加される例が示されているが、スライス選択方向ＳＳとは別の方向にＭＰＧを印加してもよい。また、スキャンＳＣにおいて、スライス選択方向ＳＳにＭＰＧが印加されるシーケンス１、位相エンコード方向ＰＥにＭＰＧが印加されるシーケンス２、および読出し方ＲＥにＭＰＧが印加されるシーケンス３を実行し、各シーケンスに対応した画像を生成するようにしてもよい。この場合、シーケンス１、２、および３のうちのいずれか一つのシーケンスから得られた画像Ａと、他のシーケンスから得られた画像Ｂとの位置合わせを行い、位置合わせされた後の画像ＡおよびＢに対して歪み量マップを用いて歪み補正を行うことができる。尚、画像ＡおよびＢの位置合わせを行う場合、Affine変換などの手法を用いて位置合わせすることが望ましい。
（３）第３の形態
第１および第２の形態では、シングルショットで画像を取得する例について説明した。第３の形態では、マルチショットで画像を取得する例について説明する。
【００９７】
図１５は、マルチショットで画像を取得するためのスキャンＳＣの一例を示す図である。
【００９８】
スキャンＳＣでは、ｍ個のシーケンス群ＰＧ１〜ＰＧｍと、ｎ個のシーケンス群ＮＧ１〜ＮＧｎとが実行される。
【００９９】
先ず、ｍ個のシーケンス群ＰＧ１〜ＰＧｍについて説明する。
【０１００】
図１５の左側には、シーケンス群ＰＧ１の具体的な一例が示されている。シーケンス群ＰＧ１は、ＥＰＩ法を用いてマルチショットにより画像を取得するためのα個のシーケンスＰＳ１〜ＰＳαを含むシーケンス群である。例えば、２ショットにより画像を取得する場合、シーケンス群ＰＧ１は、２個のシーケンスＰＳ１およびＰＳ２を含み、４ショットにより画像を取得する場合、シーケンス群ＰＧ１は、４個のシーケンスＰＳ１、ＰＳ２、ＰＳ３、およびＰＳ４を含む。シーケンスＰＳ１〜ＰＳαの各々は、補正勾配パルスＧｐｏ、位相エンコード勾配パルスＧｐ１、リードアウト勾配パルスＧｒが、マルチショット用に設計されている点を除いては、図４に示す左側のシーケンスと同じにすることができる。図１５では、シーケンスＰＳ１〜ＰＳαのうち、代表して、シーケンスＰＳ１の波形が示されている。
【０１０１】
図１５の右側には、シーケンス群ＮＧ１の具体的な一例が示されている。シーケンス群ＮＧ１は、ＥＰＩ法を用いてマルチショットにより画像を取得するためのβ個のシーケンスＮＳ１〜ＮＳβを含むシーケンス群である。例えば、２ショットにより画像を取得する場合、シーケンス群ＮＧ１は、２個のシーケンスＮＳ１およびＮＳ２を含み、４ショットにより画像を取得する場合、シーケンス群ＮＧ１は、４個のシーケンスＮＳ１、ＮＳ２、ＮＳ３、およびＮＳ４を含む。シーケンスＮＳ１〜ＮＳβの各々は、補正勾配パルスＧｐｏ、位相エンコード勾配パルスＧｐ１、リードアウト勾配パルスＧｒが、マルチショット用に設計されている点を除いては、図４に示す右側のシーケンスと同じにすることができる。図１５では、シーケンスＮＳ１〜ＮＳβのうち、代表して、シーケンスＮＳ１の波形が示されている。
【０１０２】
図１６は、第３の形態における歪み補正の説明図である。
【０１０３】
スキャンＳＣを実行し、画像生成手段１０１は、スキャンＳＣの各シーケンス群により得られたデータに基づいて、画像Ｐ１〜Ｐｍおよび画像Ｑ１〜Ｑｎを生成する。
【０１０４】
そして、画像選択手段１０２が画像Ｐ１およびＱ１を選択し、位置合わせ手段１０３が、画像Ｐ１およびＱ１の位置合わせを行う。次に、歪み量マップ生成手段１０４が、位置合わせされた後の画像Ｐ１およびＱ１に基づいて歪み量マップＤ１を生成し、補正手段１０５が、歪み量マップＤ１に基づいて画像Ｐ１およびＱ１を補正する。この歪み補正により、歪み補正後の画像Ｐ１’およびＱ１’が得られる。
【０１０５】
尚、図１６には示されていないが、画像Ｐ２〜Ｐｍ、画像Ｑ２〜Ｐｎも補正される。画像の補正方法としては、第１の形態で説明した方法又は第２の形態で説明した方法を用いることができる。
【０１０６】
第１〜第３の形態では、９０°−１８０°パルスを用いたスピンエコー型ＥＰＩのシーケンスを用いた例について説明されているが、本発明は、当該シーケンスに限定されることはなく、他のシーケンスを実行してもよい。例えば、グラジエントエコー型ＥＰＩのシーケンスを用いてもよいし、反転回復型ＥＰＩのシーケンスを用いてもよい。
【０１０７】
上記のように、第１〜第３の形態では、画像の位置合わせを行った後で歪み量マップを生成しているので、歪みが低減された画像を得ることができる。このことを検証するため、位相エンコード勾配パルスの極性を反転させてファントム画像を取得し、位置合わせが行われないファントム画像に基づいて生成された歪み量マップと、位置合わせが行われたファントム画像に基づいて生成された歪み量マップとを用いて、ファントム画像の歪み補正を行った。以下に、歪み補正の結果を示す。
【０１０８】
図１７および図１８は、歪み補正の結果の説明図である。
【０１０９】
図１７は、位置合わせが行われないファントム画像に基づいて歪み量マップを生成し、当該歪み量マップを用いて歪み補正されたファントム画像を示す図である。
【０１１０】
図１７の上段には、ファントム画像ＶおよびＷ並びに歪み量マップＭ１が示されている。
【０１１１】
画像Ｖは、正極性の位相エンコード勾配パルスを有するシーケンスを用いて生成された画像であり、画像Ｗは、負極性の位相エンコード勾配パルスを有するシーケンスを用いて生成された画像である。
【０１１２】
歪み量マップＭ１は、画像ＶおよびＷの位置合わせをせずに生成されたマップである。
【０１１３】
図１７の下段には、歪み量マップＭ１を用いて画像ＶおよびＷの歪み補正を行うことにより得られた画像Ｖ１およびＷ１が示されている。
【０１１４】
歪み量マップＭ１を参照すると、ファントムの領域が白く描出されており、ファントム画像ＶおよびＷの歪み量が過大評価されていることがわかる。したがって、歪み補正後の画像Ｖ１およびＷ１には、ファントムの周囲にアーチファクトが現れており、歪み量マップＭ１では十分な歪み補正が行われていないことがわかる。
【０１１５】
次に、図１８について説明する。
【０１１６】
図１８は、位置合わせが行われたファントム画像に基づいて歪み量マップを生成し、当該歪み量マップを用いて歪み補正されたファントム画像を示す図である。
【０１１７】
図１８の上段には、画像ＶおよびＷ並びに歪み量マップＭ２が示されている。
【０１１８】
歪み量マップＭ２は、画像ＶおよびＷの位置合わせを行った後に生成されたマップである。
【０１１９】
図１８の下段には、歪み量マップＭ２を用いて画像ＶおよびＷの歪み補正を行うことにより得られた画像Ｖ２およびＷ２が示されている。
【０１２０】
図１８の歪み量マップＭ２では、ファントムの領域の色が背景色に類似している。したがって、歪み量マップＭ２は、歪み量マップＭ１と比較すると、歪み量の過大評価が回避されていることがわかる。このため、歪み補正後の画像Ｖ２およびＷ２は、図１７に示す歪み補正後の画像Ｖ１およびＷ１と比較すると、ファントムの周囲のアーチファクトが低減されており、歪み量マップＭ２を用いることにより歪み補正が改善されていることがわかる。
【符号の説明】
【０１２１】
１  ＭＲＩ装置
２  マグネット
３  テーブル
３ａ  クレードル
４  受信コイル
５  制御部
６  送信器
７  勾配磁場電源
８  受信器
９  記憶部
１０  操作部
１１  表示部
１２  被検体
２１  ボア
２２  超伝導コイル
２３  勾配コイル
２４  ＲＦコイル
１０１  画像生成手段
１０２  画像選択手段
１０３  位置合わせ手段
１０４  歪み量マップ生成手段
１０５  補正手段
【図１】
2019050937.tif 000003
【図２】
2019050937.tif 000004
【図３】
2019050937.tif 000005
【図４】
2019050937.tif 000006
【図５】
2019050937.tif 000007
【図６】
2019050937.tif 000008
【図７】
2019050937.tif 000009
【図８】
2019050937.tif 000010
【図９】
2019050937.tif 000011
【図１０】
2019050937.tif 000012
【図１１】
2019050937.tif 000013
【図１２】
2019050937.tif 000014
【図１３】
2019050937.tif 000015
【図１４】
2019050937.tif 000016
【図１５】
2019050937.tif 000017
【図１６】
2019050937.tif 000018
【図１７】
2019050937.tif 000019
【図１８】
2019050937.tif 000020