(57)【要約】      （修正有）【課題】大動脈内の血液の信号が十分に低減できない場合でも、大動脈を検出することができる技術を提供する。【解決手段】大動脈を検出するための教師データＶを学習させ、識別器Ｃを作成する。大動脈を検出するための教師データＶには、大動脈内の血液の信号が十分に低減されている教師データｖ１１〜ｖ１ｐと、大動脈内の血液の信号があまり低減されていない教師データｖ２１〜ｖ２ｑとが含まれている。したがって、識別器Ｃを用いることにより、血液の信号が十分に低減されている血管だけでなく、血液の信号があまり低減されていない血管も検出することができる。【選択図】図６
2015023944.tif 000002
【特許請求の範囲】
【請求項１】
血管を横切る複数の面の各々の画像であって、前記血管内の血液の信号が低減されている第１の画像と、前記血管内の血液の信号が前記第１の画像よりも低減されていない第２の画像とを含んでいる画像を作成する画像作成手段と、
前記第１の画像から前記血管を検出するための第１の教師データと、前記第２の画像から前記血管を検出するための第２の教師データとに基づいて作成された識別器を用いて、前記画像から前記血管を検出する検出手段と、
を有する血管検出装置。
【請求項２】
前記検出手段は、
前記画像ごとに、前記血管の検索範囲を規定する検索領域を設定し、前記検索領域の中から前記血管を検出する、請求項１に記載の血管検出装置。
【請求項３】
前記検出手段は、
前記検索領域に、前記血管を検出するためのウィンドウを設定し、前記識別器を用いて、前記ウィンドウ内に前記血管が存在している可能性が高いか低いかを判断する、請求項２に記載の血管検出装置。
【請求項４】
前記検出手段は、
前記検索領域内で前記ウィンドウのパラメータを変化させ、前記ウィンドウ内に前記血管が存在している可能性が最も高いときの前記ウィンドウのパラメータの値を求める、請求項３に記載の血管検出装置。
【請求項５】
前記パラメータは、前記ウィンドウの回転角および前記ウィンドウの大きさである、請求項４に記載の血管検出装置。
【請求項６】
前記検出手段は、
前記画像の血管を検出するときに求められた前記ウィンドウの回転角に基づいて、前記画像の隣りに位置する画像に前記検索領域を設定する、請求項５に記載の血管検出装置。
【請求項７】
前記画像ごとに検出された前記血管の中から、血液の信号が低減されている血管を選択する選択手段を有する、請求項１〜６のうちのいずれか一項に記載の血管検出装置。
【請求項８】
前記選択手段により選択された血管の中心位置を特定する特定手段を有する、請求項７に記載の血管検出装置。
【請求項９】
前記特定手段は、
前記選択手段により選択された血管の中心位置に基づいて、前記選択手段により選択されなかった血管の中心位置を特定する、請求項８に記載の血管検出装置。
【請求項１０】
前記特定手段により特定された血管の中心位置に基づいて、造影剤を検出するためのトラッカー領域を設定する設定手段を有する、請求項８又は９に記載の血管検出装置。
【請求項１１】
血管を横切る複数の面の各々の画像であって、前記血管内の血液の信号が低減されている第１の画像と、前記血管内の血液の信号が前記第１の画像よりも低減されていない第２の画像とを含んでいる画像を作成する画像作成手段と、
前記第１の画像から前記血管を検出するための第１の教師データと、前記第２の画像から前記血管を検出するための第２の教師データとに基づいて作成された識別器を用いて、前記画像から前記血管を検出する検出手段と、
を有する磁気共鳴装置。
【請求項１２】
血管を横切る複数の面の各々の画像であって、前記血管内の血液の信号が低減されている第１の画像と、前記血管内の血液の信号が前記第１の画像よりも低減されていない第２の画像とを含んでいる画像を作成する画像作成処理と、
前記第１の画像から前記血管を検出するための第１の教師データと、前記第２の画像から前記血管を検出するための第２の教師データとに基づいて作成された識別器を用いて、前記画像から前記血管を検出する検出処理と、
を計算機に実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、血管を検出する血管検出装置、被検体の画像を取得し、取得した画像から血管を検出する磁気共鳴装置、および血管を検出するためのプログラムに関する。
【背景技術】
【０００２】
従来より、造影剤を用いて被検体を撮影する磁気共鳴イメージング装置が知られている（特許文献１）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００９−２６１９０４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
造影剤を用いて被検体を撮影する方法の一例として、オペレータが、画像データの中から大動脈の位置を見つけて、大動脈の位置に、造影剤を検出するためのトラッカー領域を設定する方法がある。この方法では、トラッカー領域に所定量の造影剤が到達したときに、撮影が実行される。しかし、経験の浅いオペレータは、大動脈の位置を見つけるのに長時間を要したり、大動脈とは別の場所にトラッカー領域を設定しまうことがある。したがって、大動脈を短時間で自動的に検出する手法の開発が望まれている。大動脈を短時間で自動的に検出する手法の一例が、論文“Med Imag Tech Vol.31 No.2 March 2013”で発表されている。
【０００５】
上記の論文では、スライス勾配磁場のＺdephaserの面積を大きくすることより、アキシャル画像に描出される大動脈内の血液の信号をできるだけ低減させて、大動脈を検出している。しかし、心臓の収縮期では、大動脈の血流速度が低下するため、大動脈内の血液が高信号になることがある。したがって、実際には、高信号の血液の影響を受けて、大動脈内の血液の信号を十分に低減することができない場合がある。この場合、大動脈とは別の組織を大動脈として誤検出したり、あるいは大動脈を検出することができないという問題がある。そこで、大動脈内の血液の信号が十分に低減できない場合でも、大動脈を検出することができる技術が望まれている。
【課題を解決するための手段】
【０００６】
本発明の第１の観点は、血管を横切る複数の面の各々の画像であって、前記血管内の血液の信号が低減されている第１の画像と、前記血管内の血液の信号が前記第１の画像よりも低減されていない第２の画像とを含んでいる画像を作成する画像作成手段と、
前記第１の画像から前記血管を検出するための第１の教師データと、前記第２の画像から前記血管を検出するための第２の教師データとに基づいて作成された識別器を用いて、前記画像から前記血管を検出する検出手段とを有する血管検出装置である。
【０００７】
本発明の第２の観点は、血管を横切る複数の面の各々の画像であって、前記血管内の血液の信号が低減されている第１の画像と、前記血管内の血液の信号が前記第１の画像よりも低減されていない第２の画像とを含んでいる画像を作成する画像作成手段と、
前記第１の画像から前記血管を検出するための第１の教師データと、前記第２の画像から前記血管を検出するための第２の教師データとに基づいて作成された識別器を用いて、前記画像から前記血管を検出する検出手段とを有する磁気共鳴装置である。
【０００８】
本発明の第３の観点は、血管を横切る複数の面の各々の画像であって、前記血管内の血液の信号が低減されている第１の画像と、前記血管内の血液の信号が前記第１の画像よりも低減されていない第２の画像とを含んでいる画像を作成する画像作成処理と、
前記第１の画像から前記血管を検出するための第１の教師データと、前記第２の画像から前記血管を検出するための第２の教師データとに基づいて作成された識別器を用いて、前記画像から前記血管を検出する検出処理とを計算機に実行させるためのプログラムである。
【発明の効果】
【０００９】
第１の教師データと第２の教師データとに基づいて作成された識別器を用いることにより、血管内の血液の信号が十分に低減できない場合でも、血管を検出することができる。
【図面の簡単な説明】
【００１０】
【図１】本発明の一形態の磁気共鳴装置の概略図である。
【図２】本形態で実行されるスキャンを示す図である。
【図３】撮影部位を概略的に示す図である。
【図４】本形態において被検体１２を撮影するときのフローを示す図である。
【図５】ローカライザスキャンＬＳの説明図である。
【図６】識別器Ｃを作成するために使用された教師データを概略的に示す図である。
【図７】アキシャル画像ＤＡ１に設定された検索領域Ｒ１を示す図である。
【図８】検索領域Ｒ１に設定されたウィンドウＷを示す図である。
【図９】識別器Ｃを用いて、大動脈ＡがウィンドウＷの先端部に内接するときのウィンドウＷの回転角θ、ピクセル数ｘ、およびピクセル数ｙを求めるときの説明図である。
【図１０】アキシャル画像ＤＡ１の大動脈ＡがウィンドウＷの先端部に内接しているときのウィンドウＷのパラメータ（θ，ｘ，ｙ）＝（θ１，ｘ１，ｙ１）を示す図である。
【図１１】アキシャル画像ＤＡ２に設定された検索領域Ｒ２を示す図である。
【図１２】検索領域Ｒ２に設定されたウィンドウＷを示す図である。
【図１３】識別器Ｃを用いて、アキシャル画像ＤＡ２の大動脈ＡがウィンドウＷの先端部に内接するときのウィンドウＷの回転角θ、ピクセル数ｘ、およびピクセル数ｙを求めるときの説明図である。
【図１４】アキシャル画像ＤＡ２の大動脈ＡがウィンドウＷの先端部に内接しているときのウィンドウＷのパラメータ（θ，ｘ，ｙ）＝（θ２，ｘ２，ｙ２）を示す図である。
【図１５】アキシャル画像ＤＡ３に設定された検索領域Ｒ３を示す図である。
【図１６】検索領域Ｒ３に設定されたウィンドウＷを示す図である。
【図１７】アキシャル画像ＤＡ３の大動脈ＡがウィンドウＷの先端部に内接しているときのウィンドウＷのパラメータ（θ，ｘ，ｙ）＝（θ３，ｘ３，ｙ３）を示す図である。
【図１８】アキシャル画像ＤＡ１〜ＤＡｍの各々に対して決定されたウィンドウＷの回転角θ、ピクセル数ｘ、およびピクセル数ｙを示す図である。
【図１９】識別器Ｄを作成するために使用された教師データを概略的に示す図である。
【図２０】アキシャル画像ＤＡ１における大動脈Ａ内の血液の信号が十分に低減されているか否かを判断する方法の説明図である。
【図２１】アキシャル画像ＤＡ２における大動脈Ａ内の血液の信号が十分に低減されているか否かを判断する方法の説明図である。
【図２２】ウィンドウＷに設定された円Ｃを概略的に示す図である。
【図２３】大動脈Ａの中心位置ｐ１を示す図である。
【図２４】アキシャル画像ＤＡ１、ＤＡ３〜ＤＡｍについて求められた大動脈Ａの中心位置ｐ１、ｐ３〜ｐｍを概略的に示す。
【図２５】仮の中心位置ｋの求め方の一例を示す図である。
【図２６】ウィンドウＷの先端部に内接する円を示す図である。
【図２７】仮の中心位置ｋが円Ｃの外側に位置する場合の一例を示す図である。
【図２８】アキシャル画像ＤＡ２における大動脈Ａの中心位置ｐ２を概略的に示す図である。
【図２９】トラッカー領域Ｒｔを概略的に示す図である。
【発明を実施するための形態】
【００１１】
以下、発明を実施するための形態について説明するが、本発明は、以下の形態に限定されることはない。
【００１２】
図１は、本発明の一形態の磁気共鳴装置の概略図である。
磁気共鳴装置（以下、「ＭＲ装置」と呼ぶ。ＭＲ：Magnetic Resonance）１００は、マグネット２、テーブル３、受信コイル４などを有している。
【００１３】
マグネット２は、被検体１２が収容されるボア２１を有している。また、マグネット２は、超伝導コイル、勾配コイル、およびＲＦコイルなどが内蔵されている。
【００１４】
テーブル３は、クレードル３ａを有している。クレードル３ａは、ボア２１内に移動できるように構成されている。クレードル３ａによって、被検体１２はボア２１に搬送される。
【００１５】
受信コイル４は、被検体１２に取り付けられている。受信コイル４は、被検体１２からの磁気共鳴信号を受信する。
【００１６】
ＭＲ装置１００は、更に、造影剤注入装置５、送信器６、勾配磁場電源７、受信器８、制御部９、操作部１０、および表示部１１などを有している。
【００１７】
造影剤注入装置５は、被検体１２に造影剤を注入する。
送信器６はＲＦコイルに電流を供給し、勾配磁場電源７は勾配コイルに電流を供給する。
受信器８は、受信コイル４から受け取った信号に対して、検波などの信号処理を実行する。
【００１８】
制御部９は、表示部１１に必要な情報を伝送したり、受信器８から受け取ったデータに基づいて画像を再構成するなど、ＭＲ装置１００の各種の動作を実現するように、ＭＲ装置１００の各部の動作を制御する。制御部９は、画像作成手段９１〜トラッカー領域設定手段９５などを有している。
【００１９】
画像作成手段９１は、被検体１２のアキシャル画像ＤＡ１〜ＤＡｍ（図５参照）を作成する。
検出手段９２は、アキシャル画像ＤＡ１〜ＤＡｍの各々から大動脈を検出する。
選択手段９３は、アキシャル画像ごとに検出された大動脈の中から、血液の信号が十分に低減されている大動脈を選択する。
特定手段９４は、アキシャル画像ごとに大動脈の中心位置を特定する。
トラッカー領域設定手段９５は、造影剤を検出するためのトラッカー領域を設定する。
【００２０】
制御部９は、画像作成手段９１〜トラッカー領域設定手段９５を構成する一例であり、所定のプログラムを実行することにより、これらの手段として機能する。制御部９は血管検出装置に相当する。
【００２１】
操作部１０は、オペレータにより操作され、種々の情報を制御部９に入力する。表示部１１は種々の情報を表示する。
ＭＲ装置１００は、上記のように構成されている。
【００２２】
図２は本形態で実行されるスキャンを示す図、図３は撮影部位を概略的に示す図である。
本形態では、ローカライザスキャンＬＳと本スキャンＭＳなどが実行される。
【００２３】
ローカライザスキャンＬＳは、スライス位置やトラッカー領域Ｒｔを設定するときに使用される画像を取得するためのスキャンである。トラッカー領域Ｒｔは、大動脈Ａに流れる造影剤を検出するために設定される領域である。
ローカライザスキャンＬＳの後に本スキャンＭＳが実行される。
【００２４】
本スキャンＭＳでは、被検体に造影剤が注入され、トラッカー領域Ｒｔから造影剤を検出するためのシーケンスが繰り返し実行される。そして、トラッカー領域Ｒｔに所定量の造影剤が注入したときに、肝臓の画像を取得するためのイメージングシーケンスが実行される。以下に、ローカライザスキャンＬＳおよび本スキャンＭＳを実行するときのフローについて説明する。
【００２５】
図４は、本形態において被検体１２を撮影するときのフローを示す図である。
ステップＳＴ１では、ローカライザスキャンＬＳ（図２参照）を実行する。
【００２６】
図５は、ローカライザスキャンＬＳの説明図である。
本形態では、ローカライザスキャンＬＳは、磁場不均一に強いＳＳＦＳＥ（Single
Shot Fast Spin Echo）法によりデータを収集するスキャンである。
【００２７】
ローカライザスキャンＬＳでは、腹部を横切る複数のアキシャル面ＡＸ１〜ＡＸｍをスキャンする。画像作成手段９１（図１参照）は、ローカライザスキャンＬＳにより収集されたデータに基づいて、アキシャル面ＡＸ１〜ＡＸｍの画像ＤＡ１〜ＤＡｍを作成する。以下では、アキシャル面の画像を「アキシャル画像」と呼ぶ。アキシャル面ＡＸ１〜ＡＸｍは大動脈Ａを横切っているので、アキシャル画像ＤＡ１〜ＤＡｍには、大動脈Ａの断面が描出される。アキシャル画像ＤＡ１〜ＤＡｍの大動脈Ａの断面の位置を比較すると、大動脈Ａの断面は、Ｓ側からＩ側に近づくに従って、背骨を中心にして反時計方向ＡＮに回転していることが分かる。
【００２８】
ローカライザスキャンＬＳで使用されるシーケンスは、流速の速い血液のエコー信号ができるだけ収束しないように設計されている。したがって、流速の速い血液の信号は十分に低減することができる。例えば、大動脈Ａを流れる血液の流速について考えると、大動脈Ａを流れる血液の流速は速いので、大動脈Ａの血液の信号は十分に低減することができる。したがって、アキシャル画像ＤＡ１〜ＤＡｍの各々の大動脈Ａ内の血液の信号は、理想的には、十分に低減されているはずである。
【００２９】
しかし、心臓の収縮期では、大動脈Ａの血流速度が低下するため、大動脈Ａ内の血液が高信号になることがある。したがって、実際には、アキシャル画像ＤＡ１〜ＤＡｍの中には、高信号の血液の影響を受けて、大動脈Ａ内の血液の信号があまり低減されていないアキシャル画像も存在する。
【００３０】
図５では、大動脈Ａ内の血液の信号が十分に低減されている場合、大動脈Ａを黒で表しており、大動脈Ａ内の血液の信号があまり低減されていない場合、大動脈Ａを白黒で表している。例えば、アキシャル画像ＤＡ１、ＤＡ３、ＤＡｍの大動脈Ａは黒で示されているので、血液の信号は十分に低減されている。一方、アキシャル画像ＤＡ２の大動脈Ａは白黒で示されているので、血液の信号はあまり低減されていない。
アキシャル画像データＤＡ１〜ＤＡｍを作成した後、ステップＳＴ２に進む。
【００３１】
ステップＳＴ２では、アキシャル画像データＤＡ１〜ＤＡｍの各々から、脳脊髄液ＣＳＦを検出する。脳脊髄液ＣＳＦの検出方法は、例えば、“Med Imag Tech Vol.31 No.2 March 2013”に記載されている方法を用いることができる。アキシャル画像データＤＡ１〜ＤＡｍごとに脳脊髄液ＣＳＦを検出した後、ステップＳＴ３に進む。
【００３２】
ステップＳＴ３では、検出手段９２（図１参照）が、アキシャル画像データＤＡ１〜ＤＡｍの各々から大動脈Ａを検出する。本形態では、大動脈Ａを識別するための識別器を用いて大動脈Ａを検出する。以下に、識別器の作成方法について説明する。
【００３３】
識別器は、被検体を撮影する前に予め準備されている。本形態では、識別器は、機械学習により作成されている。具体的には、教師データを用意し、機械学習で教師データを学習させることにより、大動脈Ａを検出するのに適した識別器Ｃを作成している（図６参照）。
【００３４】
図６は、識別器Ｃを作成するために使用された教師データを概略的に示す図である。
教師データは、実際の複数の人間を撮影することにより得られた腹部のアキシャル画像や胸部のアキシャル画像などを用いて作成することができる。本形態では、教師データとして、２種類の教師データＶおよびＷが用いられている。教師データＶは、アキシャル画像ＤＡ１〜ＤＡｍの中から大動脈Ａを検出するためのデータであり、教師データＷは、大動脈Ａ以外の組織を検出の対象から除外するためのデータである。教師データＶおよびＷは、いずれも、ｘ０×ｙ０ピクセルの長方形領域のデータとして規定されている。以下に、これらの教師データについて順に説明する。
【００３５】
（１）教師データＶについて
教師データＶは、２種類の教師データ、即ち、教師データｖ１１〜ｖ１ｐと、教師データｖ２１〜ｖ２ｑに分けられる。
（１ａ）教師データｖ１１〜ｖ１ｐについて
教師データｖ１１〜ｖ１ｐは、大動脈の断面が長方形領域の右側に内接するように規定されている。教師データｖ１１〜ｖ１ｐは、大動脈内の血液の信号が十分に低減されているデータである。教師データｖ１１〜ｖ１ｐの大動脈は黒色で示されている。
（１ｂ）教師データｖ２１〜ｖ２ｑについて
教師データｖ２１〜ｖ２ｑは、大動脈の断面が長方形領域の右側に内接するように規定されている。教師データｖ２１〜ｖ２ｑは、大動脈内の血液の信号があまり低減されていないデータである。教師データｖ２１〜ｖ２ｑの大動脈は白黒模様で示されている。
【００３６】
（２）教師データＷについて
教師データＷには、大動脈以外の組織の信号を表す教師データｗ１〜ｗｓが含まれている。大動脈以外の組織は、例えば、肝臓や腎臓等である。
【００３７】
このような教師データＶおよびＷを用意し、機械学習により教師データＶおよびＷを学習させることにより、アキシャル画像の中から大動脈を検出するのに適した識別器Ｃを作成する。
【００３８】
上記のようにして作成された識別器Ｃを用いて、アキシャル画像ごとに大動脈を検出する。以下に、大動脈の検出方法について、図７〜図１８を参照しながら説明する。
【００３９】
先ず、検出手段９２は、アキシャル画像ＤＡ１〜ＤＡｍの中からアキシャル画像ＤＡ１を取り出し、アキシャル画像ＤＡ１に、大動脈Ａの検索範囲を規定する検索領域を設定する（図７参照）。
【００４０】
図７は、アキシャル画像ＤＡ１に設定された検索領域Ｒ１を示す図である。
検出手段９２は、脳脊髄液ＣＳＦから水平に延びるラインＬ０と、脳脊髄液ＣＳＦを中心としてラインＬ０に対してθ＝θｅ１だけ傾いたラインＬ１とを規定する。ラインＬ０とラインＬ１とに挟まれる領域が、大動脈Ａの検索範囲を規定する検索領域Ｒ１として設定される。検索領域Ｒ１は、大動脈Ａが含まれるように設定される。θｅ１は、例えば、θｅ１＝６０°である。
【００４１】
次に、検出手段９２は、検索領域Ｒ１に、大動脈Ａを検出するときに使用するウィンドウを設定する（図８参照）。
【００４２】
図８は、検索領域Ｒ１に設定されたウィンドウＷを示す図である。
ウィンドウＷは、ｘ×ｙのピクセルサイズを有している。検出手段９２は、脳脊髄液ＣＳＦを基準にして、ラインＬ０上にウィンドウＷを設定する。ウィンドウＷを設定した後、検出手段９２は、検索領域Ｒ１内で、脳脊髄液ＣＳＦを中心としてウィンドウＷを回転させることによりウィンドウＷの回転角θを変更し、更にウィンドウＷのサイズ（ピクセル数ｘ、ピクセル数ｙ）も変更する。そして、識別器Ｃ（図６参照）を用いて、大動脈ＡがウィンドウＷの先端部に内接するときのウィンドウＷの回転角θ、ピクセル数ｘ、およびピクセル数ｙを求める（図９参照）。
【００４３】
図９は、識別器Ｃを用いて、大動脈ＡがウィンドウＷの先端部に内接するときのウィンドウＷの回転角θ、ピクセル数ｘ、およびピクセル数ｙを求めるときの説明図である。
【００４４】
検出手段９２は、アキシャル画像ＤＡ１に設定されたウィンドウＷの３つのパラメータ（θ，ｘ，ｙ）のうちのいずれかのパラメータを変更するたびに、ウィンドウＷ内のピクセルデータを抽出する。識別器Ｃは、抽出されたデータに基づいて、アキシャル画像ＤＡ１の大動脈ＡがウィンドウＷの先端部に内接している可能性が高いか低いかを判断するための値を出力する。検出手段９２は、識別器Ｃの出力値Ｏｕｔに基づいて、アキシャル画像ＤＡ１の大動脈ＡがウィンドウＷの先端部に内接している可能性が高いか低いかを判断する。本形態では、ウィンドウＷの先端部に大動脈Ａが内接している可能性が高いほど出力値Ｏｕｔが大きくなるように、識別器Ｃが構成されている。したがって、識別器Ｃの出力値Ｏｕｔが最大のときに、アキシャル画像ＤＡ１の大動脈ＡがウィンドウＷの先端部に内接している可能性が最も高いと判断することができる。ここでは、回転角θ＝θ１、ピクセル数ｘ＝ｘ１、およびピクセル数ｙ＝ｙ１のときに、識別器Ｃの出力値Ｏｕｔが最大になったとする。したがって、ウィンドウＷのパラメータ（θ，ｘ，ｙ）＝（θ１，ｘ１，ｙ１）のときに、アキシャル画像ＤＡ１の大動脈ＡがウィンドウＷの先端部に内接していると判断される。図１０に、アキシャル画像ＤＡ１の大動脈ＡがウィンドウＷの先端部に内接しているときのウィンドウＷのパラメータ（θ，ｘ，ｙ）＝（θ１，ｘ１，ｙ１）を示す。
【００４５】
ウィンドウＷの回転角θ１およびピクセル数ｙ１は、アキシャル画像ＤＡ１における大動脈Ａの位置を表しており、ウィンドウＷのピクセル数ｘ１はアキシャル画像ＤＡ１における大動脈Ａの血管径を表している。したがって、ウィンドウＷのパラメータ（θ，ｘ，ｙ）＝（θ１，ｘ１，ｙ１）を求めることによって、アキシャル画像ＤＡ１の中から大動脈Ａを検出することができる。
【００４６】
本形態では、大動脈を検出するための教師データＶに、大動脈内の血液の信号が十分に低減されている教師データｖ１１〜ｖ１ｐが含まれている（図６参照）。したがって、アキシャル画像ＤＡ１のように、大動脈Ａ内の血液の信号が十分に低減されている場合でも、大動脈Ａを検出することができる。
【００４７】
ウィンドウＷのパラメータ（θ，ｘ，ｙ）＝（θ１，ｘ１、ｙ１）を求めた後、検出手段９２は、アキシャル画像ＤＡ２に、大動脈Ａの検索範囲を規定する検索領域を設定する（図１１参照）。
【００４８】
図１１は、アキシャル画像ＤＡ２に設定された検索領域Ｒ２を示す図である。
アキシャル画像ＤＡ２は、アキシャル画像ＤＡ１の隣りに位置している。したがって、アキシャル画像ＤＡ２における大動脈Ａの位置は、アキシャル画像ＤＡ１における大動脈Ａの位置に近いと考えられる。そこで、検出手段９２は、アキシャル画像ＤＡ１における大動脈Ａの位置を表す回転角θ１を基準にして、回転角θ１を含む角度範囲θｓ２≦θ≦θｅ２を規定し、この角度範囲θｓ２≦θ≦θｅ２を検索領域Ｒ２として設定する。
【００４９】
尚、大動脈Ａの断面は、Ｓ側からＩ側に近づくに従って、背骨を中心にして反時計方向ＡＮに回転している。したがって、アキシャル画像ＤＡ２における大動脈Ａは、θ＝θ１に対して反時計方向ＡＮに位置している可能性が高いので、検索領域Ｒ２のうちのθｓ２≦θ≦θ１の角度範囲は狭くし、θ１≦θ≦θｅ２の角度範囲を広く設定してもよい。このように、大動脈Ａの走行方向に基づいて検索領域Ｒ２を設定することにより、大動脈Ａが存在している可能性が低い領域を、検索領域Ｒ２から除外することができる。
【００５０】
次に、検出手段９２は、検索領域Ｒ２に、大動脈Ａを検出するときに使用するウィンドウを設定する（図１２参照）。
【００５１】
図１２は、検索領域Ｒ２に設定されたウィンドウＷを示す図である。
検出手段９２は、脳脊髄液ＣＳＦを基準にして、アキシャル画像ＤＡ２のθ＝θｓ２上にウィンドウＷを設定する。ウィンドウＷを設定した後、検出手段９２は、検索領域Ｒ２内で、脳脊髄液ＣＳＦを中心としてウィンドウＷの回転角θを変更し、更にウィンドウＷのサイズ（ピクセル数ｘ、ピクセル数ｙ）を変更する。そして、識別器Ｃを用いて、アキシャル画像ＤＡ２の大動脈ＡがウィンドウＷの先端部に内接するときのウィンドウＷの回転角θ、ピクセル数ｘ、およびピクセル数ｙを求める（図１３参照）。
【００５２】
図１３は、識別器Ｃを用いて、アキシャル画像ＤＡ２の大動脈ＡがウィンドウＷの先端部に内接するときのウィンドウＷの回転角θ、ピクセル数ｘ、およびピクセル数ｙを求めるときの説明図である。
【００５３】
検出手段９２は、アキシャル画像ＤＡ２に設定されたウィンドウＷの３つのパラメータ（θ，ｘ，ｙ）のうちのいずれかのパラメータを変更するたびに、ウィンドウＷ内のピクセルデータを抽出する。識別器Ｃは、抽出されたデータに基づいて、アキシャル画像ＤＡ２の大動脈ＡがウィンドウＷの先端部に内接している可能性が高いか低いかを判断するための値を出力する。検出手段９２は、識別器Ｃの出力値Ｏｕｔに基づいて、アキシャル画像ＤＡ２の大動脈ＡがウィンドウＷの先端部に内接している可能性が高いか低いかを判断する。本形態では、ウィンドウＷの先端部に大動脈Ａが内接している可能性が高いほど識別器Ｃの出力値Ｏｕｔが大きくなるように、識別器Ｃが構成されている。したがって、識別器Ｃの出力値Ｏｕｔが最大のときに、アキシャル画像ＤＡ２の大動脈ＡがウィンドウＷの先端部に内接している可能性が最も高いと判断することができる。ここでは、回転角θ＝θ２、ピクセル数ｘ＝ｘ２、およびピクセル数ｙ＝ｙ２のときに、識別器Ｃの出力値Ｏｕｔが最大になったとする。したがって、ウィンドウＷのパラメータ（θ，ｘ，ｙ）＝（θ２，ｘ２，ｙ２）のときに、アキシャル画像ＤＡ２の大動脈ＡがウィンドウＷの先端部に内接していると判断される。図１４に、アキシャル画像ＤＡ２の大動脈ＡがウィンドウＷの先端部に内接しているときのウィンドウＷのパラメータ（θ，ｘ，ｙ）＝（θ２，ｘ２，ｙ２）を示す。
【００５４】
ウィンドウＷの回転角θ２およびピクセル数ｙ２は、アキシャル画像ＤＡ２における大動脈Ａの位置を表しており、ウィンドウＷのピクセル数ｘ２はアキシャル画像ＤＡ２における大動脈Ａの血管径を表している。したがって、ウィンドウＷのパラメータ（θ，ｘ，ｙ）＝（θ２，ｘ２，ｙ２）を求めることによって、アキシャル画像ＤＡ２の中から大動脈Ａを検出することができる。
【００５５】
本形態では、大動脈を検出するための教師データＶに、大動脈内の血液の信号があまり低減されていない教師データｖ２１〜ｖ２ｑが含まれている（図６参照）。したがって、アキシャル画像ＤＡ２のように、大動脈Ａ内の血液の信号があまり低減されていない場合でも、大動脈Ａを検出することができる。
【００５６】
ウィンドウＷのパラメータ（θ，ｘ，ｙ）＝（θ２，ｘ２、ｙ２）を求めた後、検出手段９２は、アキシャル画像ＤＡ３に、大動脈Ａの検索範囲を規定する検索領域を設定する（図１５参照）。
【００５７】
図１５は、アキシャル画像ＤＡ３に設定された検索領域Ｒ３を示す図である。
検出手段９２は、アキシャル画像ＤＡ２における大動脈Ａの位置を表す回転角θ２を基準にして、回転角θ２を含む角度範囲θｓ３≦θ≦θｅ３を規定し、この角度範囲θｓ３≦θ≦θｅ３を、検索領域Ｒ３として設定する。
【００５８】
次に、検出手段９２は、検索領域Ｒ３に、大動脈Ａを検出するときに使用するウィンドウを設定する（図１６参照）。
【００５９】
図１６は、検索領域Ｒ３に設定されたウィンドウＷを示す図である。
検出手段９２は、脳脊髄液ＣＳＦを基準にして、アキシャル画像ＤＡ３のθ＝θｓ３上にウィンドウＷを設定する。ウィンドウＷを設定した後、検出手段９２は、検索領域Ｒ３内で、脳脊髄液ＣＳＦを中心としてウィンドウＷの回転角θを変更し、更にウィンドウＷのサイズ（ピクセル数ｘ、ピクセル数ｙ）を変更する。そして、検出手段９２は、図９および図１３を参照しながら説明したように、ウィンドウＷの（θ，ｘ，ｙ）のうちのいずれかのパラメータを変更するたびに、ウィンドウＷ内のピクセルデータを抽出し、識別器Ｃの出力値Ｏｕｔが最大になるときのウィンドウＷのパラメータ（θ，ｘ，ｙ）を求める。ここでは、（θ，ｘ，ｙ）＝（θ３，ｘ３，ｙ３）のときに、識別器Ｃの出力値Ｏｕｔが最大になったとする。したがって、ウィンドウＷのパラメータ（θ，ｘ，ｙ）＝（θ３，ｘ３，ｙ３）のときに、アキシャル画像ＤＡ３の大動脈ＡがウィンドウＷの先端部に内接していると判断される。図１７に、アキシャル画像ＤＡ３の大動脈ＡがウィンドウＷの先端部に内接しているときのウィンドウＷのパラメータ（θ，ｘ，ｙ）＝（θ３，ｘ３，ｙ３）を示す。
【００６０】
ウィンドウＷの回転角θ３およびピクセル数ｙ３は、アキシャル画像ＤＡ３における大動脈Ａの位置を表しており、ウィンドウＷのピクセル数ｘ３はアキシャル画像ＤＡ３における大動脈Ａの血管径を表している。したがって、ウィンドウＷのパラメータ（θ，ｘ，ｙ）＝（θ３，ｘ３，ｙ３）を求めることによって、アキシャル画像ＤＡ３の中から大動脈Ａを検出することができる。
【００６１】
ウィンドウＷのパラメータ（θ，ｘ，ｙ）＝（θ３，ｘ３，ｙ３）を求めた後、以下同様に、残りのアキシャル画像にも検索領域を設定し、大動脈ＡがウィンドウＷの先端部に内接するときのウィンドウＷの回転角θ、ピクセル数ｘ、およびピクセル数ｙを決定する。したがって、全てのアキシャル画像ＤＡ１〜ＤＡｍに対してウィンドウＷの回転角θ、ピクセル数ｘ、およびピクセル数ｙを決定することができる。図１８に、アキシャル画像ＤＡ１〜ＤＡｍの各々に対して決定されたウィンドウＷの回転角θ、ピクセル数ｘ、およびピクセル数ｙを示す。ウィンドウＷの回転角θ、ピクセル数ｘ、およびピクセル数ｙが決定されるので、アキシャル画像ＤＡ１〜ＤＡｍの各々から大動脈Ａを検出することができる。
アキシャル画像ＤＡ１〜ＤＡｍから大動脈Ａを検出した後、ステップＳＴ４０に進む。
【００６２】
ステップＳＴ４０では、アキシャル画像ＤＡ１〜ＤＡｍの大動脈Ａの中心位置を特定する。以下に、ステップＳＴ４０において、大動脈Ａの中心位置を特定する方法について説明する。
【００６３】
本形態では、ステップＳＴ３において、アキシャル画像ＤＡ１〜ＤＡｍごとに、大動脈ＡがウィンドウＷの先端部に内接するときのパラメータ（θ，ｘ，ｙ）が求められている。したがって、パラメータ（θ，ｘ，ｙ）の値から、アキシャル画像ＤＡ１〜ＤＡｍの各々の大動脈Ａの中心位置を求めることが可能である。大動脈Ａの中心位置を求める方法としては、Mean Shift法を用いることが考えられる。ただし、アキシャル画像ＤＡ２のように大動脈Ａ内の血液の信号があまり低減できていない場合、大動脈内の信号分布のばらつきが大きくなり、大動脈内に信号ムラが発生する。大動脈内に信号ムラが発生した場合、Mean Shift法では、大動脈の中心位置を十分な精度で検出することができないことがある。そこで、本形態では、大動脈内に信号ムラが現れる場合であっても、大動脈の中心位置の検出精度を高めることができるように、大動脈の中心位置を求める方法を工夫している。以下に、ステップＳＴ４０において大動脈の中心位置を求める方法について説明する。尚、ステップＳＴ４０は、ステップＳＴ４〜ＳＴ６を有しているので、ステップＳＴ４〜ＳＴ６について順に説明する。
【００６４】
ステップＳＴ４では、選択手段９３（図１参照）が、アキシャル画像ＤＡ１〜ＤＡｍの各々から検出された大動脈Ａの中から、血液の信号が十分に低減されている大動脈Ａを選択する。本形態では、識別器を用いて、血液の信号が十分に低減されている大動脈Ａを選択する。以下に、血液の信号が十分に低減されている大動脈Ａを選択するための識別器の作成方法について説明する。
【００６５】
識別器は、被検体を撮影する前に予め準備されている。本形態では、識別器は、機械学習により作成されている。具体的には、教師データを用意し、機械学習で教師データを学習させることにより、血液の信号が十分に低減されている大動脈Ａを検出するための識別器Ｄを作成している（図１９参照）。
【００６６】
図１９は、識別器Ｄを作成するために使用された教師データを概略的に示す図である。
教師データは、実際の複数の人間を撮影することにより得られた腹部のアキシャル画像や胸部のアキシャル画像などを用いて作成することができる。本形態では、教師データとして、２種類の教師データＦおよびＧが用いられている。
【００６７】
教師データＦは、血液の信号が十分に低減されている大動脈を検出するための教師データｖ１１〜ｖ１ｐを有している。図１９に示す教師データｖ１１〜ｖ１ｐは、図６に示す教師データｖ１１〜ｖ１ｐと同じである。
【００６８】
一方、教師データＧは、血液の信号があまり低減されていない大動脈を検出の対象から除外するための教師データｖ２１〜ｖ２ｑを有している。図１９に示す教師データｖ２１〜ｖ２ｑは、図６に示す教師データｖ２１〜ｖ２ｑと同じである。
【００６９】
このような教師データＦおよびＧを用意し、機械学習により教師データＦおよびＧを学習させることにより、血液の信号が十分に低減されている大動脈を選択するための識別器Ｄを作成する。
【００７０】
上記のようにして作成された識別器Ｄを用いて、アキシャル画像ＤＡ１〜ＤＡｍにおける大動脈Ａごとに、血液の信号が十分に低減されているか否かを判断する。以下に、この判断方法について説明する。
【００７１】
選択手段９３は、先ず、アキシャル画像ＤＡ１における大動脈Ａ内の血液の信号が十分に低減されているか否かを判断する（図２０参照）。
【００７２】
図２０は、アキシャル画像ＤＡ１における大動脈Ａ内の血液の信号が十分に低減されているか否かを判断する方法の説明図である。
選択手段９３は、アキシャル画像ＤＡ１の大動脈ＡがウィンドウＷの先端部に内接するときのウィンドウＷのパラメータ（θ，ｘ，ｙ）＝（θ１，ｘ１，ｙ１）を確認し、パラメータ（θ１，ｘ１，ｙ１）のウィンドウＷ内から、ピクセルデータを抽出する。そして、識別器Ｄは、抽出されたデータに基づいて、アキシャル画像ＤＡ１の大動脈Ａ内の血液の信号が十分に低減されているか否かを判断するための値を出力する。選択手段９３は、識別器Ｄの出力値Ｏｕｔに基づいて、血液の信号が十分に低減されているか否かを判断する。本形態では、血液の信号が十分に低減されている場合、出力値Ｏｕｔが閾値ＴＨ以上になり、血液の信号があまり低減されていない場合、出力値Ｏｕｔが閾値ＴＨより小さくなるように、識別器Ｄが構成されている。したがって、識別器Ｄの出力値Ｏｕｔを閾値ＴＨと比較することにより、大動脈Ａ内の血液の信号が十分に低減されているか否かを判断することができる。ここでは、識別器Ｄの出力値Ｏｕｔが閾値ＴＨ以上であるとする。したがって、選択手段９３は、アキシャル画像ＤＡ１の大動脈Ａ内の血液の信号は十分に低減されていると判断する。
【００７３】
次に、選択手段９３は、アキシャル画像ＤＡ２における大動脈Ａ内の血液の信号が十分に低減されているか否かを判断する。
【００７４】
図２１は、アキシャル画像ＤＡ２における大動脈Ａ内の血液の信号が十分に低減されているか否かを判断する方法の説明図である。
選択手段９３は、アキシャル画像ＤＡ２の大動脈ＡがウィンドウＷの先端部に内接するときのウィンドウＷのパラメータ（θ，ｘ，ｙ）＝（θ２，ｘ２，ｙ２）を確認し、パラメータ（θ２，ｘ２，ｙ２）のウィンドウＷ内から、ピクセルデータを抽出する。そして、識別器Ｄは、抽出されたデータに基づいて、アキシャル画像ＤＡ２の大動脈Ａ内の血液の信号が十分に低減されているか否かを判断するための値を出力する。選択手段９３は、識別器Ｄの出力値Ｏｕｔと閾値ＴＨとを比較する。そして、その比較結果に基づいて、血液の信号が十分に低減されているか否かを判断する。ここでは、識別器Ｄの出力値Ｏｕｔが閾値ＴＨより小さいとする。したがって、選択手段９３は、アキシャル画像ＤＡ２の大動脈Ａ内の血液の信号はあまり低減されていないと判断する。
【００７５】
以下同様に、選択手段９３は、アキシャル画像ＤＡ３〜ＤＡｍごとに、大動脈Ａ内の血液の信号が十分に低減されているか否かを判断する。したがって、アキシャル画像ＤＡ１〜ＤＡｍの各々から検出された大動脈Ａの中から、血液の信号が十分に低減されている大動脈Ａを選択することができる。ここでは、アキシャル画像ＤＡ１、ＤＡ３〜ＤＡｍの大動脈Ａ内の血液の信号が十分に低減されていると判断されたとする。血液の信号が十分に低減されている大動脈Ａを選択した後、ステップＳＴ５に進む。
【００７６】
ステップＳＴ５では、特定手段９４（図１参照）が、アキシャル画像ＤＡ１、ＤＡ３〜ＤＡｍの各々の大動脈Ａの中心位置を特定する。以下では、先ず、アキシャル画像ＤＡ１の大動脈Ａの中心位置を特定する方法について説明する。
【００７７】
図２２および図２３は、アキシャル画像ＤＡ１の大動脈Ａの中心位置を特定する方法の説明図である。
特定手段９４は、先ず、ウィンドウＷの先端部に内接する円Ｃを設定する。図２２に、設定された円Ｃを概略的に示す。一般的に大動脈Ａの輪郭は円形に近いので、大動脈Ａの輪郭は円Ｃに近似することができる。円Ｃを設定した後、特定手段９４は、円Ｃの中心位置Ｃ０を開始点としてMean Shift法により、アキシャル画像ＤＡ１における大動脈Ａの中心位置を求める。図２３に、求められた大動脈Ａの中心位置ｐ１を示す。
【００７８】
以下同様に、血液の信号が十分に低減されている他のアキシャル画像ＤＡ３〜ＤＡｍについても、Mean Shift法により、大動脈Ａの中心位置を求める。図２４に、アキシャル画像ＤＡ１、ＤＡ３〜ＤＡｍについて求められた大動脈Ａの中心位置ｐ１、ｐ３〜ｐｍを概略的に示す。
【００７９】
このように、ステップＳＴ５を実行することにより、アキシャル画像データＤＡ１、ＤＡ３〜ＤＡｍにおける大動脈Ａの中心位置ｐ１、ｐ３〜ｐｍを特定することができる。しかし、アキシャル画像ＤＡ２における大動脈Ａの中心位置はまだ特定されていない。そこで、アキシャル画像ＤＡ２における大動脈Ａの中心位置を特定するために、ステップＳＴ６に進む。
【００８０】
ステップＳＴ６では、アキシャル画像ＤＡ２における大動脈Ａの中心位置を求める。
図２５〜図２７は、アキシャル画像ＤＡ２における大動脈Ａの中心位置を求める方法の説明図である。
特定手段９４は、先ず、アキシャル画像ＤＡ２における大動脈Ａの仮の中心位置ｋを求める。図２５に、仮の中心位置ｋの求め方の一例を示す。仮の中心位置ｋは、例えば、アキシャル画像ＤＡ１の大動脈Ａの中心位置ｐ１と、アキシャル画像ＤＡ３の大動脈Ａの中心位置ｐ３と用いて、線形補間により求めることができる。
【００８１】
仮の中心位置ｋを求めた後、アキシャル画像ＤＡ２のウィンドウＷの先端部に内接する円を考える（図２６参照）。
【００８２】
図２６は、ウィンドウＷの先端部に内接する円を示す図である。
一般的に大動脈Ａの輪郭は円形に近いので、大動脈Ａの輪郭は円Ｃに近似できる。そこで、特定手段９４は、仮の中心位置ｋが円Ｃの内側に位置している場合、仮の中心位置ｋを、アキシャル画像ＤＡ２における大動脈Ａの中心位置と決定する。図２６では、仮の中心位置ｋが円Ｃの内側に位置しているので、仮の中心位置ｋを、アキシャル画像ＤＡ２における大動脈Ａの中心位置と決定する。
【００８３】
一方、仮の中心位置ｋが円Ｃの外側に位置する場合がある（図２７参照）。
【００８４】
図２７は、仮の中心位置ｋが円Ｃの外側に位置する場合の一例を示す図である。
仮の中心位置ｋが円Ｃの外側に位置する場合、円Ｃの円周上の点の中で、仮の中心位置ｋに最も近い位置を、アキシャル画像ＤＡ２における大動脈Ａの中心位置ｐ２と決定する。
【００８５】
したがって、図２６および図２７に示すように、仮の中心位置ｋに基づいて、アキシャル画像ＤＡ２における大動脈Ａの中心位置ｐ２を決定することできる。図２８に、アキシャル画像ＤＡ２における大動脈Ａの中心位置ｐ２を概略的に示す。中心位置ｐ２を決定したら、ステップＳＴ７進む。
【００８６】
ステップＳＴ７では、トラッカー領域設定手段９５（図１参照）が、アキシャル画像ＤＡ１〜ＤＡｍの各々の大動脈Ａの中心位置ｐ１〜ｐｍに基づいて、大動脈Ａに、造影剤を検出するためのトラッカー領域を設定する。図２９に、設定されたトラッカー領域Ｒｔを概略的に示す。トラッカー領域Ｒｔは、大動脈Ａの血管壁に内接するように設定される。トラッカー領域Ｒｔを設定した後、ステップＳＴ８に進む。
【００８７】
ステップＳＴ８０では、本スキャンＭＳ（図２参照）が実行される。本スキャンＭＳでは、被検体に造影剤が注入され、トラッカー領域Ｒｔから造影剤を検出するためのシーケンスが繰り返し実行される。そして、トラッカー領域Ｒｔに所定量の造影剤が注入したときに、肝臓の画像を取得するためのイメージングシーケンスが実行され、フローが終了する。
【００８８】
本形態では、大動脈を検出するための教師データＶに、大動脈内の血液の信号が十分に低減されている教師データｖ１１〜ｖ１ｐと、大動脈内の血液の信号があまり低減されていない教師データｖ２１〜ｖ２ｑとが含まれている（図６参照）。したがって、アキシャル画像ＤＡ１〜ＤＡｍの各々の大動脈Ａ内の血液の信号が十分に低減されているか否かに拘らず、アキシャル画像ＤＡ１〜ＤＡｍの各々の大動脈Ａを検出することができる。
【００８９】
また、本形態では、アキシャル画像ＤＡ１〜ＤＡｍの各々の大動脈Ａの中心位置を特定する場合、血液の信号が十分に低減されている大動脈Ａについては、Mean Shift法により大動脈の中心位置を特定する。しかし、血液の信号があまり低減されていない大動脈Ａについては、線形補間を用いて中心位置を決定している。したがって、信号ムラが大きくなる大動脈であっても、十分な精度で大動脈の中心位置を決定することができる。尚、線形補間により求めた仮の中心位置ｋが、大動脈Ａの実際の中心位置に十分近いと考えられる場合は、ウィンドウＷの先端部に内接する円Ｃを考えることなく、仮の中心位置ｋを大動脈Ａの中心位置と決定してもよい。
【００９０】
本形態では、脳脊髄液ＣＳＦを基準にして大動脈Ａを検出しているが、脳脊髄液ＣＳＦとは別の部位を基準にして、大動脈Ａを検出してもよい。また、本形態では大動脈Ａを検出しているが、本発明は、大動脈Ａとは別の血管を検出する場合にも適用することができる。
【００９１】
本形態では、アキシャル画像に基づいて、大動脈を検出しているが、アキシャル面とは別の面（例えば、アキシャル面に対して斜めに交差するオブリーク面）の画像データに基づいて大動脈を決定してもよい。
【００９２】
本形態では、脳脊髄液ＣＳＦを含むようにウィンドウＷが設定されている。しかし、大動脈Ａの位置を検出することができるのであれば、脳脊髄液ＣＳＦはウィンドウＷから外れていてもよい。また、ウィンドウＷは矩形状であるが、別の形状（例えば、楕円形状）であってもよい。
【符号の説明】
【００９３】
２  マグネット
３  テーブル
３ａ  クレードル
４  受信コイル
５  造影剤注入装置
６  送信器
７  勾配磁場電源
８  受信器
９  制御部
１０  操作部
１１  表示部
１２  被検体
２１  ボア
９１  画像作成手段
９２  検出手段
９３  選択手段
９４  特定手段
９５  トラッカー領域設定手段
【図１】
2015023944.tif 000003
【図２】
2015023944.tif 000004
【図３】
2015023944.tif 000005
【図４】
2015023944.tif 000006
【図５】
2015023944.tif 000007
【図６】
2015023944.tif 000008
【図７】
2015023944.tif 000009
【図８】
2015023944.tif 000010
【図９】
2015023944.tif 000011
【図１０】
2015023944.tif 000012
【図１１】
2015023944.tif 000013
【図１２】
2015023944.tif 000014
【図１３】
2015023944.tif 000015
【図１４】
2015023944.tif 000016
【図１５】
2015023944.tif 000017
【図１６】
2015023944.tif 000018
【図１７】
2015023944.tif 000019
【図１８】
2015023944.tif 000020
【図１９】
2015023944.tif 000021
【図２０】
2015023944.tif 000022
【図２１】
2015023944.tif 000023
【図２２】
2015023944.tif 000024
【図２３】
2015023944.tif 000025
【図２４】
2015023944.tif 000026
【図２５】
2015023944.tif 000027
【図２６】
2015023944.tif 000028
【図２７】
2015023944.tif 000029
【図２８】
2015023944.tif 000030
【図２９】
2015023944.tif 000031