(57)【要約】【課題】光音響測定装置において、計測光の発光時刻と計測データのサンプリング開始時刻との差がばらつくことによる、測定の精度低下を防ぐ。【解決手段】光照射手段と、パルス光の照射によって被検体内で発生した音響波を受信し、電気信号に変換する音響波探触子と、電気信号をデジタルデータに変換する変換手段と、変換手段を駆動するためのサンプリングクロックを生成するクロック生成手段と、パルス光の発光タイミングを表す発光タイミング信号を生成する通知手段と、サンプリングクロックと発光タイミング信号とのずれを表す位相データを生成する位相検出手段と、デジタルフィルタを用いてデジタルデータを第二のデジタルデータに変換するフィルタ手段と、第二のデジタルデータに基づいて、被検体内の情報を表す画像を生成する画像生成手段と、を有し、フィルタ手段は、位相データに基づいてデジタルフィルタが用いるフィルタ係数を決定する。【選択図】図５
2015123322.tif 000002
【特許請求の範囲】
【請求項１】
被検体にパルス光を照射する光照射手段と、
前記パルス光の照射によって被検体内で発生した音響波を受信し、電気信号に変換する音響波探触子と、
前記電気信号をデジタルデータに変換する変換手段と、
前記変換手段を駆動するためのサンプリングクロックを生成するクロック生成手段と、
前記パルス光の発光タイミングを表す発光タイミング信号を生成する通知手段と、
前記サンプリングクロックと前記発光タイミング信号とのずれを検出し、当該ずれを表す位相データを生成する位相検出手段と、
デジタルフィルタを用いて前記デジタルデータを第二のデジタルデータに変換するフィルタ手段と、
前記第二のデジタルデータに基づいて、前記被検体内の情報を表す画像を生成する画像生成手段と、
を有し、
前記フィルタ手段は、前記位相データに基づいて、デジタルフィルタが用いるフィルタ係数を決定する
ことを特徴とする、被検体情報取得装置。
【請求項２】
前記フィルタ手段は、前記デジタルフィルタを用いて、前記デジタルデータがサンプリングされた時刻以外におけるデータを推定し、当該推定したデータに基づいて、前記発光タイミング信号を基準とした第二のデジタルデータを生成する
ことを特徴とする、請求項１に記載の被検体情報取得装置。
【請求項３】
前記光照射手段は、被検体に対するパルス光の照射を複数回行い、
複数回のパルス光の照射によって、被検体内の同一の位置から発生した音響波に対応するデジタルデータを収集し、当該デジタルデータを平均化する平均化手段をさらに有する
ことを特徴とする、請求項１または２に記載の被検体情報取得装置。
【請求項４】
前記位相検出手段は、前記サンプリングクロックに同期した、より高い周波数のクロックをカウントすることで位相データを生成し、
前記フィルタ手段は、前記位相データに基づいて、予め記憶された複数のフィルタ係数の組から使用するフィルタ係数を選択する
ことを特徴とする、請求項１から３のいずれか１項に記載の被検体情報取得装置。
【請求項５】
前記デジタルフィルタはキュービック法による補間を行うフィルタである
ことを特徴とする、請求項１から４のいずれか１項に記載の被検体情報取得装置。
【請求項６】
前記通知手段は、被検体に照射された前記パルス光を検出し、当該パルス光を検出した場合に、発光タイミング信号を生成する
ことを特徴とする、請求項１から５のいずれか１項に記載の被検体情報取得装置。
【請求項７】
前記光照射手段を発光させるためのトリガを生成する発振制御手段をさらに有し、
前記通知手段は、前記トリガに基づいて発光タイミング信号を生成する
ことを特徴とする、請求項１から５のいずれか１項に記載の被検体情報取得装置。
【請求項８】
前記光照射手段は、被検体に対するパルス光の照射を複数回行い、
前記画像生成手段は、複数回のパルス光の照射によって、被検体内の異なる位置から発生した音響波に対応するデジタルデータを合成して画像を生成する
ことを特徴とする、請求項１から７のいずれか１項に記載の被検体情報取得装置。
【請求項９】
被検体にパルス光を照射する光照射手段と、
前記パルス光の照射によって被検体内で発生した音響波を受信し、電気信号に変換する音響波探触子と、
を有する被検体情報取得装置の制御方法であって、
前記光照射手段からパルス光を発生させる光照射ステップと、
サンプリングクロックを生成するクロック生成ステップと、
前記音響波探触子によって変換された電気信号を、前記サンプリングクロックを用いてデジタルデータに変換する変換ステップと、
前記パルス光の発光タイミングを表す発光タイミング信号を生成する通知ステップと、
前記サンプリングクロックと前記発光タイミング信号とのずれを検出し、当該ずれを表す位相データを生成する位相検出ステップと、
デジタルフィルタを用いて前記デジタルデータを第二のデジタルデータに変換するフィルタステップと、
前記第二のデジタルデータに基づいて、前記被検体内の情報を表す画像を生成する画像生成ステップと、
を含み、
前記フィルタステップでは、前記位相データに基づいて、デジタルフィルタが用いるフィルタ係数を決定する
ことを特徴とする、被検体情報取得装置の制御方法。
【請求項１０】
前記フィルタステップでは、前記デジタルフィルタを用いて、前記デジタルデータがサンプリングされた時刻以外におけるデータを推定し、当該推定したデータに基づいて、前記発光タイミング信号を基準とした第二のデジタルデータを生成する
ことを特徴とする、請求項９に記載の被検体情報取得装置の制御方法。
【請求項１１】
前記光照射ステップでは、パルス光の照射を複数回行い、
複数回のパルス光の照射によって、被検体内の同一の位置から発生した音響波に対応するデジタルデータを収集し、当該デジタルデータを平均化する平均化ステップをさらに含む
ことを特徴とする、請求項９または１０に記載の被検体情報取得装置の制御方法。
【請求項１２】
前記位相検出ステップでは、前記サンプリングクロックに同期した、より高い周波数のクロックをカウントすることで位相データを生成し、
前記フィルタステップでは、前記位相データに基づいて、予め記憶された複数のフィルタ係数の組から使用するフィルタ係数を選択する
ことを特徴とする、請求項９から１１のいずれか１項に記載の被検体情報取得装置の制御方法。
【請求項１３】
前記デジタルフィルタはキュービック法による補間を行うフィルタである
ことを特徴とする、請求項９から１２のいずれか１項に記載の被検体情報取得装置の制御方法。
【請求項１４】
前記通知ステップでは、被検体に照射された前記パルス光を検出し、パルス光を検出した場合に、発光タイミング信号を生成する
ことを特徴とする、請求項９から１３のいずれか１項に記載の被検体情報取得装置の制御方法。
【請求項１５】
前記光照射ステップでは、前記光照射手段を発光させるためのトリガを生成し、
前記通知ステップでは、前記トリガに基づいて発光タイミング信号を生成する
ことを特徴とする、請求項９から１３のいずれか１項に記載の被検体情報取得装置の制御方法。
【請求項１６】
前記光照射ステップでは、パルス光の照射を複数回行い、
前記画像生成ステップでは、複数回のパルス光の照射によって、被検体内の異なる位置から発生した音響波に対応するデジタルデータを合成して画像を生成する
ことを特徴とする、請求項９から１５のいずれか１項に記載の被検体情報取得装置の制御方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被検体内の情報を取得する被検体情報取得装置に関する。
【背景技術】
【０００２】
近年、光音響効果を利用して被検体内の情報を画像化する光イメージング装置の研究が、医療分野で積極的に進められている。
パルスレーザ光などの計測光を被検体である生体に照射すると、計測光が被検体内の生体組織で吸収される際に音響波が発生する。この現象を光音響効果と呼び、光音響効果により発生した音響波を光音響波と呼ぶ。被検体を構成する組織は、光エネルギーの吸収率がそれぞれ異なるため、発生する光音響波の音圧も異なったものとなる。この光音響波を受信し、受信信号を解析することにより、被検体内で発生した光音響波の音圧分布を取得することができ、当該音圧分布に基づいて、被検体内の情報を表す画像を生成することができる。
【０００３】
具体的には、生体組織から発生した光音響波を、音響波探触子（トランスデューサ）を用いてアナログの電気信号に変換し、Ａ／Ｄ変換器によってデジタルデータに変換する。そして、変換されたデジタルデータを処理することによって、画像を再構成する。
光音響効果を利用した測定装置として、例えば特許文献１には、測定対象領域を複数の領域に分割して複数回の測定を行い、それぞれ取得したデータを結合して一枚の画像を再構成する装置が記載されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１２−００５６２２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
光音響効果を利用して被検体の測定を行う場合、パルス光の発光タイミングにばらつきがあると、データの発生タイミングにずれが生じ、再構成された画像がずれるおそれがある。パルス光の発光タイミングのばらつきは、例えばレーザ装置のジッタによって発生する。特許文献１に記載の発明では、この問題を解決するために、パルス光の発光タイミングを検出し、検出したタイミングに基づいて取得したデータを補正することで、ずれの無い画像生成を可能にしている。
【０００６】
しかし、既存の光音響測定装置には、以下のような別の課題がある。
すなわち、Ａ／Ｄ変換回路のサンプリング開始時刻とパルス光の発光時刻との差が一定ではないため、再構成された画像がずれてしまうという課題である。当該ばらつきは、特許文献１に記載されているような、パルス光のジッタに起因するばらつきよりも小さい。具体的には、装置のサンプリングクロック周期以下のずれである。従って、画像上に生じるずれも微小なものとなる。しかし、このような微小なずれに起因するアーチファクトは、アーチファクトであることを認識しづらいため、医用の光音響イメージング装置においては誤診につながるおそれがある。
すなわち、当該課題により発生するずれの量は、特許文献１に記載の装置が対応可能なずれよりも小さいが、その影響はむしろ大きいと言える。
【０００７】
また、Ａ／Ｄ変換回路のサンプリング開始時刻とパルス光の発光時刻との差にばらつきがあると、取得したデータを平均化した場合に、高周波成分が減衰してしまう。光音響測
定装置においては、信号のＳＮ比を向上させるため、複数回のパルス光を照射し、得られたデータを平均化する手法が多くとられる。しかし、サンプリング開始時刻とパルス光の発光時刻との差にばらつきがある、すなわち、データの先頭がずれている複数のデータに対して平均化処理を行うと、高周波成分が減衰してしまい、測定の精度が低下してしまう。
【０００８】
本発明は、このような従来技術の課題に鑑みてなされたものであり、光音響測定装置において、計測光の発光時刻と計測データのサンプリング開始時刻との差がばらつくことによる、測定の精度低下を防ぐ技術を提供することを目的とする。
【課題を解決するための手段】
【０００９】
上記課題を解決するために、本発明に係る被検体情報取得装置は、
被検体にパルス光を照射する光照射手段と、前記パルス光の照射によって被検体内で発生した音響波を受信し、電気信号に変換する音響波探触子と、前記電気信号をデジタルデータに変換する変換手段と、前記変換手段を駆動するためのサンプリングクロックを生成するクロック生成手段と、前記パルス光の発光タイミングを表す発光タイミング信号を生成する通知手段と、前記サンプリングクロックと前記発光タイミング信号とのずれを検出し、当該ずれを表す位相データを生成する位相検出手段と、デジタルフィルタを用いて前記デジタルデータを第二のデジタルデータに変換するフィルタ手段と、前記第二のデジタルデータに基づいて、前記被検体内の情報を表す画像を生成する画像生成手段と、を有し、前記フィルタ手段は、前記位相データに基づいて、デジタルフィルタが用いるフィルタ係数を決定することを特徴とする。
【００１０】
また、本発明に係る被検体情報取得装置の制御方法は、
被検体にパルス光を照射する光照射手段と、前記パルス光の照射によって被検体内で発生した音響波を受信し、電気信号に変換する音響波探触子と、を有する被検体情報取得装置の制御方法であって、前記光照射手段からパルス光を発生させる光照射ステップと、サンプリングクロックを生成するクロック生成ステップと、前記音響波探触子によって変換された電気信号を、前記サンプリングクロックを用いてデジタルデータに変換する変換ステップと、前記パルス光の発光タイミングを表す発光タイミング信号を生成する通知ステップと、前記サンプリングクロックと前記発光タイミング信号とのずれを検出し、当該ずれを表す位相データを生成する位相検出ステップと、デジタルフィルタを用いて前記デジタルデータを第二のデジタルデータに変換するフィルタステップと、前記第二のデジタルデータに基づいて、前記被検体内の情報を表す画像を生成する画像生成ステップと、を含み、前記フィルタステップでは、前記位相データに基づいて、デジタルフィルタが用いるフィルタ係数を決定することを特徴とする。
【発明の効果】
【００１１】
本発明によれば、光音響測定装置において、計測光の発光時刻と計測データのサンプリング開始時刻との差がばらつくことによる、測定の精度低下を防ぐことができる。
【図面の簡単な説明】
【００１２】
【図１】従来技術に係る光音響測定装置の構成を示す図。
【図２】従来技術における測定タイムチャート。
【図３】従来技術に係る光音響測定装置の構成を示す第二の図。
【図４】従来技術において、高周波領域に発生する減衰を説明する図。
【図５】第一の実施形態に係る光音響測定装置の構成を示す図。
【図６】第一の実施形態において記憶されるデジタルデータの例。
【図７】第一の実施形態に係る適応フィルタ回路の構成例。
【図８】適応フィルタ回路が利用する係数の一例を表す図。
【図９】第一の実施形態における測定タイムチャート。
【図１０】第二の実施形態に係る光音響測定装置の構成を示す図。
【図１１】第二の実施形態において、高周波領域に発生する減衰を説明する図。
【図１２】第三の実施形態に係る光音響測定装置の構成を示す図。
【図１３】第四の実施形態に係る光音響測定装置の構成を示す図。
【図１４】第四の実施形態において記憶されるデジタルデータの例。
【発明を実施するための形態】
【００１３】
以下、図面を参照しながら、本発明の実施形態を詳細に説明する。なお、同一の構成要素には原則として同一の参照番号を付して、説明を省略する。また、実施形態の説明で用いる数値や材質等は、発明の範囲を限定するものではない。
【００１４】
（課題の説明）
最初に、前述した課題について、従来技術に係る光音響測定装置のシステム構成図である図１を参照しながら説明する。
光音響測定装置１００は、従来の光音響測定装置であり、レーザパルス送信部１、光音響受信部３、システム制御部４からなる。なお、符号２は被検体であり、符号５は画像データの出力端子である。
レーザパルス送信部１は、被検体に対してレーザパルス光を照射する光照射手段であり、送信基準クロック回路１１、レーザ発光制御回路１２、Ｑスイッチ１３、レーザ装置１４からなる。
送信基準クロック回路１１は、レーザ発光制御回路１２に周期的なクロック信号を供給する手段であり、レーザ発光制御回路１２は、Ｑスイッチ１３に対して、パルス光を発光させる発光トリガ（以下、発振開始信号）を送信する回路である。
符号２１は、被検体表面にレーザ光を導く光ファイバであり、符号２２は、光音響受信部３にレーザの発光時刻を通知するための光ファイバである。
【００１５】
光音響受信部３は、光音響波をアナログの電気信号（以下、光音響信号）に変換する超音波変換器３１、ＡＤ変換器３２、受信基準クロック回路３３、光検出器３４、信号処理部３５からなる。また、ＡＤ変換器３２は、増幅器３２１とＡＤ変換回路３２２からなる。
超音波変換器３１は、受信した音響波をアナログ電気信号（光音響信号）に変換する手段である。超音波変換器は、超音波探触子（音響波探触子）とも呼ばれ、例えばＣＭＵＴと呼ばれる静電容量型のセンサ等であるが、音響波を電気信号に変換することができれば、どのようなものを用いてもよい。例えば、圧電セラミックス（ＰＺＴ）を利用した変換素子や、磁性膜を用いるＭＭＵＴ（Magnetic MUT）、圧電薄膜を用いるＰＭＵＴ（Piezoelectric MUT）などを用いることができる。
受信基準クロック回路３３は、ＡＤ変換回路３２２を駆動するためのサンプリングクロックである受信基準クロックＳ３を生成するクロック生成手段である。また、光検出器３４は、レーザパルスの発光タイミングを検出するための手段である。また、信号処理部３５は、ＡＤ変換器３２で変換されたデジタルデータ（以下、光音響データ）を処理し、画像に変換する画像生成手段である。信号処理部３５は、書き込み制御回路３５１、メモリ３５２、信号処理回路３５３からなる。
【００１６】
当該装置において、システム制御部４から、レーザ発光指示がレーザパルス送信部１に送信されると、送信基準クロック回路１１が生成したクロックと同期したタイミングで、レーザ発光制御回路１２が、発振開始信号Ｓ１をＱスイッチ１３に出力する。そして、レーザ装置１４が、当該信号に応答してレーザパルス光を出射させる。
なお、本例ではレーザの発振にＱスイッチを用いる例を示すが、他の発振制御手段を用いてもよい。例えば、半導体レーザを用いる場合は、直接変調で十分に高速な応答が可能
なため、Ｑスイッチではなく変調ドライバを使用することもできる。
【００１７】
出射したパルス光は、光ファイバ２１によって被検体の表面に導かれ、被検体２に照射される。被検体にパルス光が照射されると、光音響効果によって被検体内の組織から音響波が発生する。この音響波は、超音波変換器３１によって受信され、光音響信号に変換される。変換された光音響信号は、所望の振幅になるよう増幅器３２１で増幅され、ＡＤ変換回路３２２で光音響データＳ４に変換される。
【００１８】
ＡＤ変換回路３２２に入力されるサンプリングクロックである受信基準クロックＳ３は、受信基準クロック回路３３で作られた、ジッタの少ない安定した基準クロックである。
ＡＤ変換回路３２２からは、光音響データが順次出力されるため、光検出器３４から出力される受光トリガ信号Ｓ２（発光タイミング信号）を用いてタイミングをとり、所望の数のデータをメモリ３５２に格納する。すなわち、書き込み制御回路３５１は、受光トリガ信号を検知した時刻を起点として、受信基準クロックＳ３の立ち上がりタイミングごとに、ＡＤ変換回路３２２が変換した光音響データを、所定の数だけ連続してメモリ３５２に格納する。
このようにして、メモリ３５２には、レーザパルス光の発光に対応する光音響データが格納される。
そして、信号処理回路３５３が、メモリ３５２に記憶された光音響データを読み出し、信号処理（画像再構成）を行い、画像データを生成する。生成された画像データは、出力端子５から出力される。
なお、図１では、わかりやすく説明するために、超音波変換器３１が一つである場合の構成を図示しているが、特許文献１にあるように、複数の超音波変換器を用いて同時に複数の光音響データを取得するようにしてもよい。この場合、超音波変換器３１からメモリ３５２までの構成が、複数個並列に設けられる。
【００１９】
図２は、前述した光音響測定装置１００において、信号またはデータが発生するタイミングを説明する図である。図２中、Ｓ１は発振開始信号、Ｓ２は受光トリガ信号、Ｓ３はサンプリングクロック、Ｓ４は光音響データである。
Ｓ１は、Ｑスイッチ１３に対して入力される発振開始信号であり、当該信号に応答してパルス光が照射される。Ｓ２は、光検出器３４が出力する受光トリガ信号であり、パルス光の発光を検出すると立ち上がる。受光トリガ信号Ｓ２を検出すると、書き込み制御回路３５１が、ＡＤ変換回路３２２で変換された光音響データを所定の数だけメモリ３５２に書き込む。以降の説明において、光音響データとは、サンプリングされたデジタルデータもしくは当該データの集合（ビット列）のことを指す。
【００２０】
図２（Ａ）の例では、発振開始信号Ｓ１の生成から時間Ｔ１が経過した後に受光トリガ信号Ｓ２が発生する。すなわち、時間Ｔ１は、レーザ装置１４が発光開始信号Ｓ１を取得してからレーザパルスを発光させるまでの時間と、光ファイバ２２をパルス光が伝搬する時間と、光検出器３４のディレイの合計時間である。
また、時間Ｔ２は、受光トリガ信号Ｓ２が発生してから、直後の受信基準クロックＳ３の立ち上がり（すなわちサンプリング開始時刻）までの時間である。
メモリ３５２に書き込むデータの数は、測定すべき被検体の奥行き方向（Ｚ方向）の距離の仕様から決定する。奥行き方向の距離と被検体内の音速から、被検体内で発生した音響波が超音波変換器３１に到達するまでの時間がわかる。また、サンプリングクロックである受信基準クロックＳ３の周期と前記時間から、書き込むデータ数が決定できる。
【００２１】
以上のような構成を持つ光音響測定装置の課題について、図２（Ｂ）を参照しながら説明する。図２（Ａ）との違いは、受光トリガ信号Ｓ２が発生してから、サンプリングクロックである受信基準クロックＳ３が立ち上がるまでの時間が異なっている点である。
すなわち、図２（Ａ）では、受光トリガ信号Ｓ２から時間Ｔ２後に光音響データのサンプリングが始まるが、図２（Ｂ）では、時間Ｔ３後までサンプリングが始まらない。この時間の違い（ばらつき）は、サンプリングクロックである受信基準クロックＳ３の周期より短い時間である。このばらつきは、発振開始信号Ｓ１を生成するためのクロックと、サンプリングクロックである受信基準クロックＳ３とが互いに同期していないために発生し、その時間は、パルス光の発光ごとに異なる値となる。
【００２２】
このため、複数回のパルス光を照射すると、パルス光が照射されてから光音響データのサンプリングが開始されるまでの時間が毎回異なるため、生成した光音響データに時間的なずれが生じる。すなわち、複数回のパルス光の照射に対応する複数の光音響データに基づいて画像データを生成すると、画像にずれが生じてしまう。
【００２３】
次に、取得した信号の高周波成分の減衰について説明する。
一般に、人体に対して障害を起こさないようにするため、人体に入射を許容できる単位面積当たりの光エネルギーの上限は決まっている。この値はＭＰＥ（Maximum Permissible Exposure）と呼ばれている。すなわち、被検体が人体である場合、発生させる光音響波のレベルを大きくすることには上限がある。
【００２４】
このような状況の中で光音響信号のＳＮ比を向上させるため、被検体に対して複数回のパルス光を照射し、取得した光音響データを平均化する手法が知られている。具体的には、得られる光音響データが同じ波形となるように、超音波変換器３１の位置を制御しながら、複数回のパルス光の発光を行い、複数の光音響データを取得する。
平均化は、例えば取得した光音響データを加算し、発光回数で除算してもよいし、単純に加算のみを行ってもよい。得られる光音響信号は相対的なものであるため、加算のみを行っても信号処理上の問題は生じない。本明細書では、以降、平均または平均化という言葉を使用して実施形態を説明するが、平均化は加算処理のみによって行ってもよい。
【００２５】
図３に、従来技術に係る光音響測定装置であって、取得した光音響データの平均化を行う機能を有する光音響測定装置３００の構成図を示す。当該装置は、図１に示した光音響測定装置１００に、光音響データの平均化を行う回路である平均回路３５４を追加した形態である。他の構成要素は同一である。
平均回路３５４は、メモリ３５２に記憶された光音響データ（すなわち、パルス光の発光毎に対応する光音響データ）を平均化する手段である。ここで行う平均化とは、発光毎にそれぞれ得られた光音響データのデータ列同士の平均をとることを意味している。すなわち、パルス光の発光毎にＡＤ変換回路３２２で変換された光音響データの列を、先頭から順に所定の数だけ平均化する。
そして、当該平均化したデータを用いて、信号処理回路３５３が画像データを生成する。
【００２６】
ここで、光音響データを平均化することによって発生する問題について説明する。すなわち、レーザパルスの発光毎に得られた光音響データＳ４を平均化した結果、光音響データの高周波成分が減衰するという問題である。光音響データの高周波成分とは、得られた光音響データを得られた順に時間軸上に並べた場合、周波数の高い変動成分を指す。
【００２７】
レーザパルスを複数回発光させる場合、前述したように、光音響データが同じ波形となるように、超音波変換器３１の位置が制御される。よって、複数回のレーザパルスに対応する光音響信号は、レーザパルス発光時刻を基準とすれば全て同じ波形となる。もちろん、この波形そのものを平均化した場合は、ノイズ以外に起因する波形の変化は起きないため、信号成分の減衰は基本的に発生しない。
【００２８】
しかし、前述したように、複数回パルス光を発光し、その都度光音響データを取得すると、レーザパルスの発光時刻とサンプリング開始時刻にばらつきが発生する。すなわち、取得した光音響データの先頭がずれてしまうため、当該光音響データ同士を平均すると、当該ばらつきに起因したローパスフィルタが形成され、高周波成分が減衰してしまう。
【００２９】
図４を参照しながら具体的に説明する。ここでは、レーザパルスの発光時刻とサンプリング開始時刻とのばらつきがランダムである場合を例に説明を行う。図４（Ａ）〜（Ｈ）は、信号の時間波形と周波数特性を模式的に示した図である。
【００３０】
レーザパルスの発光回数が多い場合、サンプリング時刻のばらつきがランダムだと仮定すれば、平均化した光音響データは、レーザパルスの発光毎に−Ｔ／２から＋Ｔ／２だけ時刻がずれた光音響データを加算し、発光回数で除算したものとなる。（Ｔはサンプリングクロックの周期である）
すなわち、図４（Ａ）に示した入力波形を、図４（Ｂ）で示した波形で畳み込み積分した結果得られる、図４（Ｃ）に示した波形を時間Ｔおきにサンプリングして得られたデータ（図４（Ｄ））と等しい。
もちろん、発光回数は有限であるので、図４（Ｂ）の波形は連続関数では無い。しかし、発光回数が多い場合、図４（Ｂ）に近い形となり、畳み込み積分した結果も、図４（Ｃ）に示した波形に近い形状となる。
図４（Ｅ）、図４（Ｆ）、図４（Ｇ）、図４（Ｈ）は、それぞれ、図４（Ａ）、図４（Ｂ）、図４（Ｃ）、図４（Ｄ）をフーリエ変換した結果である。図４（Ａ）と図４（Ｂ）の波形を畳み込み積分して、図４（Ｃ）の波形を計算するため、周波数軸上では、図４（Ｅ）の周波数特性と、図４（Ｆ）の周波数特性を乗算することによって、図４（Ｇ）の周波数特性が計算できる。このようにすると、図４（Ｇ）において、（１）で示した灰色部分が減衰していることがわかる。
そして、図４（Ｇ）の周波数特性の波形を時間Ｔの間隔でサンプリングすると、図４（Ｈ）のような周波数特性が得られる。このように、レーザパルス発光回数が多い場合、サンプリング時刻のばらつきがランダムだと仮定すれば、平均化された光音響データの周波数特性は、図４（Ｈ）に示したように、高い周波数領域が減衰することがわかる。
【００３１】
実際には、レーザパルス発光回数は多くない場合もあり、サンプリング時刻のばらつきがランダムに生じない場合もある。しかし、実施する平均処理が、レーザパルス発光毎にサンプリング時刻のばらついたデータを加算するフィルタと等価であることを考えれば、周波数特性のカーブの形状はさまざまではあるが、少なくとも高周波成分が減衰することは明らかである。もちろん、サンプリング時刻のばらつきが無いデータを加算する場合は、フィルタ効果が生じないので高周波成分は減衰することは無い。
【００３２】
（第一の実施形態）
第一の実施形態に係る光音響測定装置は、レーザパルス光を被検体に照射し、当該パルス光に起因して被検体内で発生した光音響波を受信および解析することで、被検体内の光学特性に関連した機能情報を可視化、すなわち画像化する装置である。
また、前述した課題を解決するために、サンプリングクロックとパルス光の発光タイミングとのずれを検出し、当該ずれに対応するデジタルフィルタを用いてデータを変換することで、ずれた位相を補正する機能を有する。これにより、レーザパルスの発光時刻とサンプリング開始時刻のばらつきに起因する測定精度の低下を抑制することができる。
【００３３】
<システム構成>
まず、図５を参照しながら、第一の実施形態に係る光音響測定装置５００の構成を説明する。本実施形態に係る光音響測定装置５００は、位相検出回路３６を有しているという点と、適応フィルタ回路３５６を有しているという点において、従来技術に係る光音響測
定装置１００と相違する。
【００３４】
本実施形態では、受信基準クロックＳ３は、受信基準クロック回路３３によって生成された、所望のサンプリング周波数のクロックである。
位相検出回路３６は、受信基準クロックＳ３をｍ倍に逓倍し、さらに、ｍ倍に逓倍したクロックをカウントする手段である。また、受光トリガ信号Ｓ２の立ち上がりパルスを検出した場合に、検出時刻における当該カウント値を、書き込み制御回路３５５に出力する。当該カウント値は、受光トリガ信号Ｓ２とサンプリングクロック（受信基準クロックＳ３）とのずれ（位相）を表すデータである。以降、当該情報を位相データＣ１と称する。具体的なカウント方法と位相の検出方法については後述する。
【００３５】
書き込み制御回路３５５は、光音響測定装置１００が有する書き込み制御回路３５１に対応する回路であるが、ＡＤ変換された光音響データをメモリ３５２に書き込む前に、取得した位相データＣ１を書き込むという点において相違する。その後、所望の数の光音響データを書き込むという点については同様である。
図６は、メモリ３５２に記憶されるデータのアドレスマップ（ＦＩＦＯを用いた場合の例）である。図６に示したように、データの先頭に位相データＣ１が記憶され、続いて、ＡＤ変換回路３２２が出力する光音響データが順次記憶される。
【００３６】
適応フィルタ回路３５６は、メモリ３５２に記憶された光音響データに対してフィルタ演算を行う手段である。
【００３７】
ここで、図７を参照しながら、適応フィルタ回路３５６について説明する。図７は、適応フィルタ回路３５６の構成の一例を表すブロック図である。
図７において、符号３５６１は適応フィルタ回路３５６の入力端子であり、符号３５６７は出力端子である。また、符号３５６２は位相データＣ１をラッチするラッチ回路であり、符号３５６３ａ〜３５６３ｃは、デジタルフィルタを構成する遅延素子である。また、符号３５６４ａ〜３５６４ｄは、デジタルフィルタが用いるフィルタ係数を記憶している係数レジスタであり、符号３５６５ａ〜３５６５ｄは、デジタルフィルタの係数を乗算する乗算器である。また、符号３５６６は、デジタルフィルタを構成する加算器である。
【００３８】
図７を参照しながら、適応フィルタ回路３５６の動作について説明する。
適応フィルタ回路３５６は、まず、メモリ３５２に記憶された光音響データを順次読み出す。前述した通り、記憶された光音響データの先頭は位相データＣ１であり、これをラッチ回路３５６２が記憶する。続く光音響データを読み出している間は、ラッチ回路３５６２に記憶された位相データＣ１は保持され、係数レジスタ３５６４ａ〜ｄにそれぞれ供給される。
係数レジスタ３５６４ａ〜ｄは、入力された位相データＣ１に対応するフィルタ係数を出力する。図８は、各係数レジスタが出力する係数を表した表である。例えば、位相データＣ１が０である場合は、Ｋ０＝０、Ｋ１＝０、Ｋ２＝１、Ｋ３＝０という係数がそれぞれ出力される。このように、各係数レジスタは、複数の組の係数の中から、位相データＣ１に応じた係数を選択する。当該係数は、キュービック補間（キュービック法）と呼ばれる３次関数のデジタルフィルタの係数である。
【００３９】
メモリ３５２に記憶されたデータは、受光トリガ信号Ｓ２が立ち上がった後、受信基準クロックＳ３が立ち上がるごとに取得したデータである。一方、本来取得すべきデータは、受光トリガ信号Ｓ２が立ち上がった時刻を開始時刻として、サンプリング周期ごとに取得したデータである。そこで、本実施形態では、実際にサンプリングされたデータにキュービック補間を適用することで、本来取得すべきデータを求める。
キュービック補間で用いる重みは、式（１）によって表される。図８に示した係数は、
式（１）を用いて算出した値である。
【数１】
2015123322.tif 000003
式（１）中のｔは、求めるべきデータからサンプルされたデータまでの時間差であり、サンプリングクロック周期で正規化された値である。すなわち、位相データＣ１に対応したｔによって重みが定まるため、サンプリング時刻以外の任意の時刻における光音響データを補間によって求めることができる。
なお、変数ａは、−０．５〜−２．０の範囲で選ぶと好適であるが、図８の表では、ａ＝−１．０として計算した値を用いている。
【００４０】
順次読み出された光音響データは、遅延素子３５６３ａ〜３５６３ｃによって遅延され、選択された係数Ｋ０〜Ｋ３と乗算器３５６５ａ〜３５６５ｄによって乗算され、加算器３５６６によって加算された後に、出力端子３５６７から出力される。
本実施形態では、以上のようにして形成したデジタルフィルタによって処理された光音響データを出力する。見方を変えれば、出力される光音響データは、異なる時刻の光音響データ（実際にＡＤ変換された光音響データ）にデジタルフィルタを適用することで生成された、受光トリガ信号Ｓ２の立ち上がりパルスが入力された時刻から始まる光音響データである。
そして、信号処理回路３５３が、適応フィルタ回路３５６が出力した光音響データを取得し、断層画像を再構成して、出力端子５より出力する。
【００４１】
図９に、第一の実施形態におけるタイミングチャートを示す。図９において、Ｓ１は発振開始信号、Ｓ２は受光トリガ信号である。また、Ｓ３は受信基準クロック（ＡＤ変換回路３２２のサンプリングクロック）である。すなわち、サンプリングクロックの立ち上がり時刻に、光音響信号から光音響データへの変換が開始される。Ｓ４は、変換後の光音響データである。
なお、図９（Ａ）は図２（Ａ）に対応し、図９（Ｂ）は図２（Ｂ）に対応する。
【００４２】
図９を参照しながら、位相検出回路３６の動作をよりわかりやすく説明する。図９において、Ｓ６は受信基準クロックＳ３をｍ倍に逓倍したクロック（以下、ｍクロック）である。カウント値Ｃｎは、受信基準クロックＳ３の立ち上がり時刻を基準としてｍクロックＳ６をダウンカウントするカウンタの値である。ここでは、ｍ＝８である場合を例に説明する。すなわち、ｍクロックＳ６は、受信基準クロックＳ３の８倍の周波数を持つクロックである。ｍクロックは、例えばＰＬＬ回路などによって生成することができる。
位相検出回路３６は、ｍクロックＳ６が立ち上がることにカウンタを減算し、受信基準クロックＳ３が立ち上がるごとにカウンタを７にリセットする。すなわち、受信基準クロックＳ３を８分割し、７から０に向けて繰り返しカウントする。
【００４３】
位相検出回路３６は、受光トリガ信号Ｓ２が立ち上がった時刻におけるカウンタの値Ｃｎをラッチし、位相データＣ１として出力する。例えば、図９（Ａ）の例では、受光トリガ信号Ｓ２が立ち上がった時刻のカウンタ値は１であり、図９（Ｂ）の例では、受光トリガ信号Ｓ２が立ち上がった時刻のカウンタ値は５である。すなわち、位相データＣ１として１および５をそれぞれ出力する。位相データが１および５であるということは、ＡＤ変
換時刻が、受光トリガ信号Ｓ２の立ち上がり時刻から、それぞれサンプリングクロック周期の１／８および５／８だけ遅れていることを意味している。
【００４４】
以上説明したように、第一の実施形態に係る光音響測定装置は、位相データＣ１をＡＤ変換したデータに付加し、当該位相データＣ１をもとにデジタルフィルタの係数を選択する。また、フィルタリング演算により、時刻をサンプリングクロック周期の１／ｍ単位でずらした光音響データを求める。すなわち、ずれ量に応じた重みをつけて補間処理を行うことで、本来取得すべき、受光トリガ信号Ｓ２を基準とした、ずれのない光音響データを推定する。この結果、レーザパルスの発光時刻とサンプリング開始時刻の差のばらつきを１／ｍに抑えることができる。
【００４５】
さらに、光音響データを複数回に分けて取得する場合であっても、レーザパルスの発光時刻とサンプリング開始時刻との差が小さくなるため、精度の高い断層画像を取得することができる。また、測定対象領域を複数領域に分割して測定を行い、得られた複数の画像を合成する場合、接続部分に生じるずれを抑制することができる。
【００４６】
本来、高周波成分の減少を防ぐためには、ｓｉｎｃ関数等を利用したタップ数の多いデジタルフィルタを用いることが好ましい。
しかし、タップ数を多くすると、遅延素子や係数レジスタ、乗算器で構成される適応型デジタルフィルタの段数が多く必要となる。すなわち、ハードウェア構成や、ソフトウェアによる処理量が増加し、それに伴う速度低下等の問題が発生する。また、タップ数の多いデジタルフィルタほど、遅延素子の段数が多いため、正しいデータが計算されるまでの遅延が大きくなる。換言すれば、データ列の両端に無駄なデータが多く生じてしまう。また、必要なデータ数以上のデータを取得する必要がある。
このように、タップ数の多いデジタルフィルタを用いると、精度は向上するが、コスト的な問題が生じることがある。一方、本実施形態のようにキュービック補間を用いると、比較的少ないタップ数（処理量）で高周波成分の減少を防ぐことができる。すなわち、コストと効果のバランスのとれたデジタルフィルタを提供することができる。
【００４７】
なお、本発明は、測定対象領域を分割して複数回の測定を行い、複数の画像データを合成する場合に特に効果が得られる。しかし、領域分割をしない場合であっても、再構成された画像データのＺ方向のずれが小さくなるため、別々の測定で得られた画像データが比較しやすくなるという効果を得ることができる。
【００４８】
なお、実施形態の説明では、ｍクロックＳ６をダウンカウントする構成を述べたが、カウント方法はこれに限定されない。例えば、カウンタをアップカウンタとしても良い。この場合であっても、図８と等価なデータを用いることによって、同一の適応型デジタルフィルタを実現できる。このように、レーザパルスの発光時刻とサンプリング開始時刻のずれを検出できれば、カウンタはどのように構成してもよい。
【００４９】
また、ｍは２以上の値であるが、ｍの上限は６４程度であることが好ましい。ｍの値を大きくすると、レーザパルスの発光時刻とサンプリング開始時刻のばらつきの影響をより小さくできるが、ｍクロックＳ６の周波数が高くなり、回路を安定して動作させることが難しくなるためである。ＦＰＧＡやＡＳＩＣを用いて当該回路を実装する場合、例えば、サンプリング周波数が１０ＭＨｚであれば、ｍクロックＳ６の最大周波数は６４０ＭＨｚとなり、十分実現可能である。
【００５０】
（第二の実施形態）
第二の実施形態は、被検体内の同一の位置に対して複数回の測定を行い、得られた光音響データを平均化することでＳＮ比を向上させる実施形態である。
【００５１】
図１０に、第二の実施形態に係る光音響測定装置１０００の構成図を示す。当該装置は、第一の実施形態に係る光音響測定装置５００に、光音響データの平均化を行う回路である平均回路３５４を追加した形態である。他の構成要素は第一の実施形態と同一である。
【００５２】
第二の実施形態でも、第一の実施形態と同様に、レーザパルスの発光時刻とサンプリング開始時刻の差のばらつきを１／ｍにすることができる。ここで、レーザパルスの発光時刻とサンプリング開始時刻の差のばらつきが１／ｍになった場合の、平均化した光音響データの周波数特性を、図１１を参照しながら説明する。
図１１は、レーザパルスの発光時刻とｍクロックＳ６の立ち上がり時刻との差がランダムに生じる場合の、信号の時間波形と周波数特性を模式的に示した図である。
【００５３】
図１１は、図４（Ａ）と同じ入力波形のデータを平均化する過程を示したものであるが、光音響データのずれが、−Ｔ／２ｍから＋Ｔ／２ｍの範囲に分布しているという点において相違する。
図４の説明と同様に、図１１（Ａ）で示した入力波形を図１１（Ｂ）で示した波形で畳み込み積分した結果得られる、図１１（Ｃ）に示した波形を時間Ｔおきにサンプリングして得られたデータ（図１１（Ｄ））が、平均化された光音響データである。
図１１（Ｅ）、図１１（Ｆ）、図１１（Ｇ）、図１１（Ｈ）は、それぞれ、図１１（Ａ）、図１１（Ｂ）、図１１（Ｃ）、図１１（Ｄ）をフーリエ変換した結果である。図１１（Ａ）と図１１（Ｂ）の波形を畳み込み積分して、図１１（Ｃ）の波形を計算するため、周波数軸上では、図１１（Ｅ）の周波数特性と、図１１（Ｆ）の周波数特性を乗算することによって、図１１（Ｇ）の周波数特性が計算できる。このようにすると、図１１（Ｇ）において、（２）で示した灰色部分が減衰した部分となるが、減衰量が、図４（Ｇ）と比較して少なくなっていることがわかる。
【００５４】
第二の実施形態も第一の実施形態と同様に、フィルタリング演算で求めた光音響データを加算あるいは平均することにより、光音響データの高周波成分の減衰を低減することができる。
また、その結果、高周波成分の減衰が少ないデータ（すなわち劣化が少ないデータ）をもとに、精度の良い解析を行うことができる。また、画像を再構成する場合、画像の劣化を抑制することができる。
このように、本発明は、レーザパルスの発光を複数回行い、発光毎に得られた光音響データを収集し平均化処理することでＳＮ比を向上する光音響イメージング装置に好適に適用することができる。
【００５５】
（第三の実施形態）
第一ないし第二の実施形態では、光検出器３４の出力を受光トリガ信号Ｓ２とした。これに対し、第三の実施形態は、Ｑスイッチ（発振制御手段）１３に入力される発振開始信号Ｓ１を受光トリガ信号Ｓ２として利用する形態である。
図１２は、第三の実施形態に係る光音響測定装置１２００の構成図である。当該装置は、光検出部３４を有しておらず、書き込み制御回路３５５が発振開始信号Ｓ１を取得してパルス光の発光タイミングを判断するという点において第一の実施形態と相違する。
【００５６】
第三の実施形態に係る光音響測定装置によれば、レーザパルスの発光タイミングを通知する信号を直接用いるため、高価な光ファイバや光検出器を用いることなく、第一の実施形態と同様の効果を得ることができる。
なお、第三の実施形態では、発振開始信号Ｓ１を利用する形態について述べたが、レーザパルスの発光タイミングがわかる信号であれば、他の信号を利用してもよい。また、本実施形態の構成は、第二の実施形態に適用してもよい。
【００５７】
（第四の実施形態）
第一ないし第三の実施形態では、受光トリガ信号とサンプリングクロックとのずれ量を検出し、本来取得すべき光音響データを演算によって求めた。一方で、このような場合は、ずれが生じたデータを、正規のサンプリング時刻とは異なる時刻にサンプリングされたデータとみなすようにしてもよい。
第四の実施形態は、ＡＤ変換された光音響データを、発生したずれ量ごとに分類して記憶することで、別時刻にサンプリングされたデータとみなす（すなわち、オーバーサンプリングと同様の効果を得る）実施形態である。
【００５８】
図１３に、第四の実施形態に係る光音響測定装置１３００の構成図を示す。当該装置は、第一の実施形態に係る光音響測定装置５００が有する適応フィルタ回路３５６を、レート変換回路３５９に置き換えた形態である。他の構成要素は第一の実施形態と同一である
【００５９】
第四の実施形態においても、第一の実施形態同様に、レーザパルスの発光毎に、ＡＤ変換された光音響データが、位相データＣ１と共にメモリ３５２に記憶される。例えば、メモリ３５２には、図６に示したようなアドレスマップでデータが記憶される。
【００６０】
第四の実施形態において、レート変換回路３５９が実行する処理について説明する。
ここでは、位相検出回路３６が生成するｍクロックがｍ＝２である場合を例に説明を行う。この場合、位相検出回路３６の出力である位相データＣ１は、０または１の値をとる。すなわち、本来のサンプリング時刻でサンプリングが行われた場合は０を出力し、本来の時刻とずれてサンプリングが行われた場合は１を出力する。レート変換回路３５９は、当該位相データＣ１を用いて、以下の様な処理を行い、不図示の第二のメモリに光音響データを記憶する。
【００６１】
例えば、レーザパルスの発光回数が二回である場合を例に説明をする。ここでは、１回目のレーザパルスの発光により得られた光音響データの位相データＣ１が０であったとする。位相データが０である場合、図１４に示すように、第２のメモリの奇数アドレスに順次記憶する（D0-1、D1-1、D2-1・・・）。次に、２回目のレーザパルスの発光により得られた光音響データの位相データが１であったとする。位相データが１である場合、図１４に示すように、第２のメモリの偶数アドレスに順次記憶する（D0-2、D1-2、D2-2・・・）。
【００６２】
このようにすると、取得した光音響データを、発生したずれ量ごとにシフトして記憶することができる。すなわち、本来のタイミングでサンプリングを行った光音響データと、本来のタイミングとずれたタイミングでサンプリングを行った二種類の光音響データを発生時間順に取得できる。
発光タイミング信号に対するサンプリング時刻のばらつきはランダムである。従って、レーザパルスを複数回発光させることによって、所望のサンプリング周波数よりも高いサンプリングクロックでサンプルした場合と等価なデータを得ることができる。
【００６３】
なお、発光タイミング信号に対するサンプリング時刻のばらつきはランダムであるため、格納先のアドレスが衝突する場合がある。
例えば、レーザパルスの２回目の発光で、得られた光音響データの位相データが０であったとする。このような場合は、１回目で記憶されている光音響データと平均を計算し平均値を同一アドレスに書き込む。この場合、奇数アドレスのデータは、２回のレーザパルスの発光により得られた光音響データの平均値となり、偶数アドレスのデータは無しとなる。データが無い偶数アドレスのデータは、奇数アドレスのデータから補間して求める。この様に、データが無いアドレスに対しては、全数のレーザパルス発光終了後に、データ
があるアドレスのデータから補間により算出する。
【００６４】
ここで、ｍクロックのｍが２以外の値である場合を説明する。位相検出回路３６が出力する位相データＣ１は、０からｍ−１の値となる。また、レート変換回路３５９内の不図示の第２のメモリには、レーザパルス１回あたりに得られる所望の光音響データの数のｍ倍のアドレス空間（容量）が必要となる。
そして、前述したように、前記位相データＣ１の値に応じて、ｋ番目の所望の光音響データを、第２のメモリのアドレス(ｋ×ｍ＋Ｃ１)に書きこむ。また、前述した様に、もし、すでにデータが書かれている場合は、書かれているデータと書き込むデータの平均を計算して同じアドレスに書き込む。例えば、書き込まれているデータが３回のレーザパルスの発光に対する平均値である場合は、書き込まれているデータの重みを３、書き込むデータの重みを１として平均をとる。
【００６５】
なお、レーザパルスの発光回数をｎとすると、ｍクロックのｍがｎより大きい場合、必ず書き込みを行わないアドレスが生じるため、補間計算が必ず必要となる。この場合、より高いサンプリングクロック相当のデータ列を得ることができるが、補間して得られたデータが多くなるため、データの精度が高くならない。
そのため、ｍの値はレーザパルスの発光回数ｎ以下の値とすることが好ましい。これにより、書き込みを行わないアドレスを少なくすることができる。この場合、補間処理が少なくなることが期待できるので、光音響データのサンプリングクロックをｍ倍に高くしたデータに近づけることが可能となる。
【００６６】
すなわち、ｍクロックのｍをｍ≦ｎである整数とし、前記位相検出部の出力（位相データＣ１）に基づき、サンプリングクロックの周期をｍ分割した時刻に対応する光音響データとして取得することにより、補間処理を行うデータを少なくすることができる。換言すると、光音響データのサンプリングクロックをｍ倍に高くしたデータに近づけることが可能となる。
【００６７】
以上、説明したように、第四の実施形態では、受信基準クロックＳ３の逓倍数ｍを、レーザパルスの発光回数ｎ以下とし、検出した位相データに基づいて、光音響データの順序を調整する。これにより、実質的に、より高周波のサンプリングクロックでＡＤ変換したデータと等価なデータを得ることが可能となる。
【００６８】
（変形例）
なお、各実施形態の説明は本発明を説明する上での例示であり、本発明は、発明の趣旨を逸脱しない範囲で適宜変更または組み合わせて実施することができる。
例えば、本発明は、上記処理の少なくとも一部を含む被検体情報取得装置として実施することもできる。また、上記処理の少なくとも一部を含む被検体情報取得装置の制御方法として実施することもできる。上記処理や手段は、技術的な矛盾が生じない限りにおいて、自由に組み合わせて実施することができる。
【００６９】
例えば、実施形態の説明では、再構成された画像データを出力する出力端子５を例示したが、出力手段は、ネットワーク入出力端子や、単に画像を記憶するための不揮発性メモリ等であってもよい。また、分割した領域に対して複数回の測定を行い、領域毎に生成した画像データを結合する場合は、出力端子５の外部に、生成された画像データを結合する手段を設けてもよい。
また、実施形態の説明では、信号処理回路としてハードウェアで処理を行う構成を例示したが、信号処理はソフトウェアを用いて行ってもよい。また、再構成された画像データを合成する手段を用いる場合も、同様にソフトウェアで処理するようにしてもよい。特に、近年マルチコア化されたＣＰＵを用いることによって、比較的短時間で画像再構成を行
うことができる。
【符号の説明】
【００７０】
１・・・レーザパルス送信部、３・・・光音響受信部、４・・・システム制御部
【図１】
2015123322.tif 000004
【図２】
2015123322.tif 000005
【図３】
2015123322.tif 000006
【図４】
2015123322.tif 000007
【図５】
2015123322.tif 000008
【図６】
2015123322.tif 000009
【図７】
2015123322.tif 000010
【図８】
2015123322.tif 000011
【図９】
2015123322.tif 000012
【図１０】
2015123322.tif 000013
【図１１】
2015123322.tif 000014
【図１２】
2015123322.tif 000015
【図１３】
2015123322.tif 000016
【図１４】
2015123322.tif 000017