(57)【要約】【課題】光音響装置において、被検体内の画像のリアルタイム表示を良好に行う。【解決手段】第１および第２の波長の光を照射する光照射部と、音響波を受信して信号を出力する受信部と、被検体の構造情報および機能情報を生成する処理部と、表示制御部を有し、光照射部は、第１の波長での第１の照射、２の波長での第２の照射、第１の波長での第３の照射を行い、受信部は、第１から第３の信号を出力し、処理部は、第１および第２の信号に基づいて第１の機能情報を生成し、第２および第３の信号に基づいて第２の機能情報を生成し、表示制御部は、第１および第２の機能情報に基づく第１および第２の表示画像を順次、表示する、被検体情報取得装置を用いる。【選択図】図４
2018126454.tif 000002
【特許請求の範囲】
【請求項１】
第１の波長の光と、前記第１の波長とは異なる第２の波長の光とを被検体に対して照射する光照射部と、
前記被検体から伝搬した音響波を受信して信号を出力する受信部と、
前記音響波に基づいて前記被検体の構造情報および機能情報を生成する処理部と、
前記構造情報および機能情報に基づく表示画像を表示部に表示させる表示制御部と、を有し、
前記光照射部は、前記第１の波長の光を前記被検体に照射する第１の照射を行い、前記第１の照射の後に前記第２の波長の光を前記被検体に照射する第２の照射を行い、前記第２の照射の後に前記第１の波長の光を前記被検体に照射する第３の照射を行い、
前記受信部は、前記第１から第３の照射のそれぞれに由来する前記音響波を受信して、第１から第３の信号を出力し、
前記処理部は、
前記第１から第３の信号のそれぞれに基づく第１から第３の構造情報の、少なくともいずれかを生成し、
前記第１および第２の信号に基づいて第１の機能情報を生成し、前記第２および第３の信号に基づいて第２の機能情報を生成し、
前記第１および第２の機能情報のそれぞれに基づいて第１および第２の表示画像を生成し、前記第１および第２の表示画像の生成には、前記第１から第３の構造情報の少なくともいずれかが用いられており、
前記表示制御部は、前記第１および第２の表示画像を順次、前記表示部に表示させる
ことを特徴とする被検体情報取得装置。
【請求項２】
前記処理部は、
前記第１から第３の構造情報のうち、少なくとも、前記第２および第３の構造情報を生成し、
前記第１の機能情報と前記第２の構造情報に基づいて前記第１の表示画像を生成し、前記第２の機能情報と前記第３の構造情報に基づいて前記第２の表示画像を生成する
ことを特徴とする請求項１に記載の被検体情報取得装置。
【請求項３】
前記機能情報は前記被検体の酸素飽和度分布を示すものであり、前記構造情報は前記被検体における血管の構造を示すものであり、
前記処理部は、前記酸素飽和度分布を前記血管の構造を用いてトリミングすることにより、前記表示画像を生成する
ことを特徴とする請求項１または２に記載の被検体情報取得装置。
【請求項４】
前記第１および第２の波長は、一方がオキシヘモグロビンとデオキシヘモグロビンの吸収係数が略等しい波長であり、他方がオキシヘモグロビンとデオキシヘモグロビンの吸収係数が異なる波長であり、
前記処理部は、
前記第１の信号に基づいて第１の吸収係数分布を生成し、前記第２の信号に基づいて第２の吸収係数分布を生成し、前記第１および第２の吸収係数分布を用いて前記酸素飽和度分布を生成し、かつ、
前記第１および第２の吸収係数分布の少なくともいずれかに対する閾値処理によって前記構造情報を生成するものである
ことを特徴とする請求項３に記載の被検体情報取得装置。
【請求項５】
前記処理部は、前記第１および第２の波長のうち、オキシヘモグロビンとデオキシヘモグロビンの吸収係数が異なる波長に由来する吸収係数分布を用いて前記構造情報を生成す
るときに、動脈である血管の幅が太くなるような補正を行う
ことを特徴とする請求項４に記載の被検体情報取得装置。
【請求項６】
前記表示制御部は、前記構造情報に基づく画像と前記機能情報に基づく画像とを重畳して前記表示部に表示させる
ことを特徴とする請求項１または２に記載の被検体情報取得装置。
【請求項７】
前記表示制御部は、前記構造情報に基づく画像と前記機能情報に基づく画像とを異なる色相で前記表示部に表示させる
ことを特徴とする請求項６に記載の被検体情報取得装置。
【請求項８】
前記受信部が、前記第１から第３の照射のそれぞれに由来する前記音響波を受信する位置は、互いに異なる
ことを特徴とする請求項１ないし７のいずれか１項に記載の被検体情報取得装置。
【請求項９】
前記受信部は、ハンドヘルド型プローブの中に配置されている
ことを特徴とする請求項８に記載の被検体情報取得装置。
【請求項１０】
前記受信部を移動させる走査部をさらに有する
ことを特徴とする請求項８に記載の被検体情報取得装置。
【請求項１１】
前記処理部は、前記受信部が前記第１から第３の照射のそれぞれに由来する前記音響波を受信した位置を用いて、前記第１および第２の機能情報を生成する
ことを特徴とする請求項８ないし１０のいずれか１項に記載の被検体情報取得装置。
【請求項１２】
第１の波長の光と、前記第１の波長とは異なる第２の波長の光とを被検体に対して照射する光照射部と、
前記被検体から伝搬した音響波を受信して信号を出力する受信部と、
前記音響波に基づく再構成を行う処理部と、
前記再構成により生成された画像データに基づく表示画像を表示部に表示させる表示制御部と、を有し、
前記処理部は、
前記第１の波長の光に由来する音響波から再構成される第１の吸収係数分布、および、前記第２の波長の光に由来する音響波から再構成される第２の吸収係数分布から、酸素飽和度分布を生成するものであり、かつ、
前記酸素飽和度分布を生成するために用いられた前記第１の吸収係数分布に関連する前記第１の波長での発光の回数を第１の回数とし、前記酸素飽和度分布を生成するために用いられた前記第２の吸収係数分布に関連する前記第２の波長での発光の回数を第２の回数としたときに、以下の式（１）〜式（３）を満たす生成回数だけ前記酸素飽和度分布を生成するものであり、
前記表示制御部は、前記表示画像として、前記酸素飽和度分布に基づく画像を前記表示部に表示させる
ことを特徴とする被検体情報取得装置。
第１の回数  ＜  生成回数    …（１）
第２の回数  ＜  生成回数    …（２）
生成回数  ≦  第１の回数＋第２の回数    …（３）
【請求項１３】
前記処理部は、前記第１の波長での発光ごと、および、前記第２の波長での発光ごとに、前記酸素飽和度分布を生成し、
前記表示制御部は、前記酸素飽和度分布が生成されるごとに、前記酸素飽和度分布に基
づく画像を順次、前記表示部に表示する
ことを特徴とする請求項１２に記載の被検体情報取得装置。
【請求項１４】
前記処理部は、前記酸素飽和度分布を生成するために用いられた前記第１の吸収係数分布、または、前記酸素飽和度分布を生成するために用いられた前記第２の吸収係数分布を用いて、前記酸素飽和度分布から血管を抽出し、
前記表示制御部は、前記表示画像として、前記血管が抽出された前記酸素飽和度分布を前記表示部に表示する
ことを特徴とする請求項１２または１３に記載の被検体情報取得装置。
【請求項１５】
第１の波長での第１の照射により被検体から発生した音響波に由来する第１の信号を取得するステップと、
前記第１の波長と異なる第２の波長での第２の照射により前記被検体から発生した音響波に由来する第２の信号を取得するステップと、
前記第１の波長での第３の照射により前記被検体から発生した音響波に由来する第３の信号を取得するステップと、
前記第１の信号および前記第２の信号に基づいて第１の機能情報を生成するステップと、
前記第２の信号および前記第３の信号に基づいて第２の機能情報を生成するステップと、
前記第１の機能情報に基づく第１の表示画像と、前記第２の機能情報に基づく第２の表示画像を順次、表示するステップと、
を有することを特徴とする表示方法。
【請求項１６】
前記第１から第３の信号のそれぞれに基づく第１から第３の構造情報のうち、少なくとも第２および第３の構造情報を生成するステップをさらに有し、
前記第１の表示画像の生成には、前記第２の構造情報が用いられており、前記第２の表示画像の生成には、前記第３の構造情報が用いられて
ことを特徴とする請求項１５に記載の表示方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被検体情報取得装置および表示方法に関する。
【背景技術】
【０００２】
光を用いて被検体の内部を画像化する光音響装置が研究・開発されている。光音響装置は、被検体に照射された光のエネルギーを吸収した光吸収体から光音響効果により発生する音響波（光音響波）を用いて再構成を行い、吸収係数分布画像を形成する。さらに、吸収係数分布画像から、被検体内の構造画像や機能画像を生成する。構造画像の例としては、被検体内の血管の位置を示す画像がある。機能画像の例としては、被検体内の光学特性に対応する特性情報分布を示す画像がある。機能画像として特に、複数の波長の光を用いることで取得される酸素飽和度分布画像が注目されている。
【０００３】
また、超音波診断装置に用いられるものと同様のハンドヘルド型プローブを用いて、容易に観察部位にアクセスできる光音響装置が研究・開発されている。このようなハンドヘルド型プローブを備える光音響装置において、リアルタイムに被検体内の構造画像や機能画像を観察可能とするための研究・開発が進められている。
【０００４】
特許文献１には、発光時刻が異なる複数の波長の光パルスを用いる光音響装置において、発光間の動きを補正して機能画像（酸素飽和度分布画像）を得る方法が開示されている。また、特許文献２には、複数の波長のうち、オキシヘモグロビンとデオキシヘモグロビンの吸収係数が等しい第１の波長の吸収係数分布画像を用いて、構造画像（血管位置を特定した画像）を得る方法が開示されている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０１６−０１３４２１号公報
【特許文献２】特開２０１５−１４２７４０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
従来、ステージに設けたプローブを機械的にスキャンし、複数の光パルスに対する吸収係数分布画像を形成する、機械走査方式の光音響装置が知られている。このような機械走査方式の光音響装置に特許文献１，２の技術を適用した場合、良好な画像が得られる。しかし、特許文献１，２の技術においては、複数の波長での発光が終了した後に表示画像が得られる。そのため、リアルタイムに被検体内の構造画像や機能画像を観察する場合、発光回数よりも、表示画像の更新回数が少なくなるという課題がある。すなわち、リアルタイムに構造画像や機能画像を表示する際に、プローブの動きや体動に対する画像表示の時間追従性が悪くなるという課題がある。特に、ハンドヘルド型プローブを備える光音響装置においてはリアルタイム表示への要望が強いため、時間追従性の改善が求められている。また、機械走査式の場合であっても、画像表示の時間追従性が良いほうが望ましい。
【０００７】
本発明は上記課題に鑑みてなされたものである。本発明の目的は、光音響装置において、被検体内の画像のリアルタイム表示を良好に行うことにある。
【課題を解決するための手段】
【０００８】
本発明は、以下の構成を採用する。すなわち、
第１の波長の光と、前記第１の波長とは異なる第２の波長の光とを被検体に対して照射
する光照射部と、
前記被検体から伝搬した音響波を受信して信号を出力する受信部と、
前記音響波に基づいて前記被検体の構造情報および機能情報を生成する処理部と、
前記構造情報および機能情報に基づく表示画像を表示部に表示させる表示制御部と、を有し、
前記光照射部は、前記第１の波長の光を前記被検体に照射する第１の照射を行い、前記第１の照射の後に前記第２の波長の光を前記被検体に照射する第２の照射を行い、前記第２の照射の後に前記第１の波長の光を前記被検体に照射する第３の照射を行い、
前記受信部は、前記第１から第３の照射のそれぞれに由来する前記音響波を受信して、第１から第３の信号を出力し、
前記処理部は、
前記第１から第３の信号のそれぞれに基づく第１から第３の構造情報の、少なくともいずれかを生成し、
前記第１および第２の信号に基づいて第１の機能情報を生成し、前記第２および第３の信号に基づいて第２の機能情報を生成し、
前記第１および第２の機能情報のそれぞれに基づいて第１および第２の表示画像を生成し、前記第１および第２の表示画像の生成には、前記第１から第３の構造情報の少なくともいずれかが用いられており、
前記表示制御部は、前記第１および第２の表示画像を順次、前記表示部に表示させる
ことを特徴とする被検体情報取得装置である。
【０００９】
本発明はまた、以下の構成を採用する。すなわち、
第１の波長の光と、前記第１の波長とは異なる第２の波長の光とを被検体に対して照射する光照射部と、
前記被検体から伝搬した音響波を受信して信号を出力する受信部と、
前記音響波に基づく再構成を行う処理部と、
前記再構成により生成された画像データに基づく表示画像を表示部に表示させる表示制御部と、を有し、
前記処理部は、
前記第１の波長の光に由来する音響波から再構成される第１の吸収係数分布、および、前記第２の波長の光に由来する音響波から再構成される第２の吸収係数分布から、酸素飽和度分布を生成するものであり、かつ、
前記酸素飽和度分布を生成するために用いられた前記第１の吸収係数分布に関連する前記第１の波長での発光の回数を第１の回数とし、前記酸素飽和度分布を生成するために用いられた前記第２の吸収係数分布に関連する前記第２の波長での発光の回数を第２の回数としたときに、以下の式（１）〜式（３）を満たす生成回数だけ前記酸素飽和度分布を生成するものであり、
前記表示制御部は、前記表示画像として、前記酸素飽和度分布に基づく画像を前記表示部に表示させる
ことを特徴とする被検体情報取得装置。
第１の回数  ＜  生成回数    …（１）
第２の回数  ＜  生成回数    …（２）
生成回数  ≦  第１の回数と第２の回数の合計数  …（３）である。
【００１０】
本発明はまた、以下の構成を採用する。すなわち、
第１の波長での第１の照射により被検体から発生した音響波に由来する第１の信号を取得するステップと、
前記第１の波長と異なる第２の波長での第２の照射により前記被検体から発生した音響波に由来する第２の信号を取得するステップと、
前記第１の波長での第３の照射により前記被検体から発生した音響波に由来する第３の
信号を取得するステップと、
前記第１の信号および前記第２の信号に基づいて第１の機能情報を生成するステップと、
前記第２の信号および前記第３の信号に基づいて第２の機能情報を生成するステップと、
前記第１の機能情報に基づく第１の表示画像と、前記第２の機能情報に基づく第２の表示画像を順次、表示するステップと、
を有することを特徴とする表示方法である。
【発明の効果】
【００１１】
本発明によれば、光音響装置において、被検体内の画像のリアルタイム表示を良好に行うことができる。
【図面の簡単な説明】
【００１２】
【図１】第１の実施形態に係る光音響装置のブロック図
【図２】第１の実施形態に係るハンドヘルド型プローブの模式図
【図３】第１の実施形態に係るコンピュータと周辺構成を示すブロック図
【図４】第１の実施形態における動作を説明するためのタイミング図
【図５】本発明の第１の実施形態における動作を説明するためのフローチャート
【図６】第１の波長λ１と第２の波長の吸収係数の一例を示すグラフ
【図７】閾値を下げる補正の説明のための図
【図８】「忠実モード」のタイミング図
【発明を実施するための形態】
【００１３】
以下に図面を参照しつつ、本発明の好適な実施の形態について説明する。ただし、以下に記載されている構成部品の寸法、材質、形状およびそれらの相対配置などは、発明が適用される装置の構成や各種条件により適宜変更されるべきものである。よって、この発明の範囲を以下の記載に限定する趣旨のものではない。
【００１４】
本発明は、被検体から伝搬する音響波を検出し、被検体内部の特性情報を生成し、取得する技術に関する。よって本発明は、被検体情報取得装置またはその制御方法、あるいは被検体情報取得方法や信号処理方法として捉えられる。本発明はまた、被検体内部の特性情報を示す画像を生成し表示する表示方法として捉えられる。本発明はまた、これらの方法をＣＰＵやメモリ等のハードウェア資源を備える情報処理装置に実行させるプログラムや、そのプログラムを格納した、コンピュータにより読み取り可能な非一時的な記憶媒体としても捉えられる。
【００１５】
本発明の被検体情報取得装置には、被検体に光（電磁波）を照射することにより被検体内で発生した音響波を受信して、被検体の特性情報を画像データとして取得する光音響効果を利用した光音響イメージング装置を含む。この場合、特性情報とは、受信された光音響波に由来する信号を用いて生成される、被検体内の複数位置のそれぞれに対応する特性値の情報である。
【００１６】
本発明に係る光音響画像データは、光照射により発生した光音響波に由来するあらゆる画像データを含む概念である。例えば、光音響画像データは、光音響波の発生音圧（初期音圧）、吸収エネルギー密度、及び吸収係数、被検体を構成する物質の濃度（酸素飽和度など）などの少なくとも１つの被検体情報の空間分布を表す画像データである。なお、互いに異なる複数の波長の光照射により発生する光音響波に基づいて、被検体を構成する物質の濃度などの、分光情報を示す光音響画像データが得られる。分光情報を示す光音響画像データは、酸素飽和度、酸素飽和度に吸収係数等の強度を重み付けした値、トータルヘ
モグロビン濃度、オキシヘモグロビン濃度、またはデオキシヘモグロビン濃度であってもよい。また、分光情報を示す光音響画像データは、グルコース濃度、コラーゲン濃度、メラニン濃度、または脂肪や水の体積分率であってもよい。
【００１７】
被検体内の各位置の特性情報に基づいて、二次元または三次元の特性情報分布が得られる。分布データは画像データとして生成され得る。特性情報は、数値データとしてではなく、被検体内の各位置の分布情報として求めてもよい。すなわち、初期音圧分布、エネルギー吸収密度分布、吸収係数分布や酸素飽和度分布などの分布情報である。
【００１８】
本発明でいう音響波とは、典型的には超音波であり、音波、音響波と呼ばれる弾性波を含む。トランスデューサ等により音響波から変換された電気信号を音響信号とも呼ぶ。ただし、本明細書における超音波または音響波という記載は、それらの弾性波の波長を限定する意図ではない。光音響効果により発生した音響波は、光音響波または光超音波と呼ばれる。光音響波に由来する電気信号を光音響信号とも呼ぶ。分布データは、光音響画像データや再構成画像データとも呼ばれる。
【００１９】
以下の実施形態では、被検体情報取得装置として、被検体に互いに異なる複数の波長のパルス光を照射し、被検体からの光音響波を受信し、被検体内の血管画像（構造画像）や酸素飽和度分布画像（機能画像）を生成する光音響装置を取り上げる。以下の実施形態ではまた、ハンドヘルド型プローブを有する光音響装置を取り上げているが、本発明は、ステージにプローブを設けて機械的にスキャンする光音響装置にも適用できる。
【００２０】
＜第１の実施形態＞
（装置構成）
以下、図１の概略ブロック図を用いて、本実施形態に係る光音響装置１の構成を説明する。光音響装置１は、プローブ１８０、信号収集部１４０、コンピュータ１５０、表示部１６０、入力部１７０、及び、光源部２００を有する。プローブ１８０には、光照射部１１３と受信部１２０が含まれる。コンピュータ１５０には、演算部１５１、記憶部１５２、制御部１５３、表示制御部１５４が含まれる。
【００２１】
光源部２００は、光ファイバ（バンドルファイバ）等の光学系１１２を経由して、光照射部１１３に光パルスを供給する。光照射部１１３は、供給された光を被検体１００に照射する。受信部１２０は、被検体１００から発生した光音響波を受信して、アナログ信号としての電気信号（光音響信号）を出力する。信号収集部１４０は、受信部１２０から出力されたアナログ信号をデジタル信号に変換し、コンピュータ１５０に出力する。
【００２２】
コンピュータ１５０は、信号収集部１４０から出力されたデジタル信号を、光音響波に由来する電気信号（光音響信号）としてメモリに記憶する。コンピュータ１５０は、記憶されたデジタル信号に対して画像再構成などの処理を行うことにより光音響画像データを生成する。また、コンピュータ１５０は、得られた光音響画像データに対して表示のための画像処理を施した後に、画像データを表示部１６０に出力する。コンピュータ１５０は、また、光音響装置１の全体の制御を行う。
【００２３】
表示部１６０は、光音響画像データに基づき光音響画像を表示する。ユーザー（医師や技師等）は、表示部１６０に表示された光音響画像を確認することにより、診断を実施できる。表示画像は、ユーザーやコンピュータ１５０からの保存指示に基づいて、コンピュータ１５０内のメモリや、光音響装置とネットワークで接続されたデータ管理システムなどに保存されても良い。入力部１７０は、ユーザーからの指示などを受け付ける。
【００２４】
（各ブロックの詳細構成）
続いて、各ブロックの好ましい構成について詳細に述べる。
【００２５】
（プローブ１８０）
図２は、本実施形態に係るプローブ１８０の模式図である。プローブ１８０は、光照射部１１３、受信部１２０、及び、ハウジング１８１を含む。ハウジング１８１は、光照射部１１３及び受信部１２０を囲う筺体である。ユーザーは、ハウジング１８１を把持することにより、プローブ１８０をハンドヘルド型プローブとして利用できる。光照射部１１３は、光学系１１２により伝搬された光パルスを被検体に照射する。なお、図中のＸＹＺ軸は、プローブを静置した場合の座標軸を示すものであり、プローブ使用時の向きを限定するものではない。
【００２６】
図２に示すプローブ１８０は、ケーブル１８２を介して、信号収集部１４０、及び光源部２００と繋がっている。ケーブル１８２は、光源部２００から光照射部１１３に光パルスを供給する光学系１１２と、受信部１２０から出力されたアナログ信号を信号収集部１４０に出力する配線（不図示）を含む。ケーブル１８２にコネクタを設け、プローブ１８０と光音響装置のその他の構成とを分離できる構成としてもよい。
【００２７】
（光学系１１２）
光学系１１２は、光源部２００が発生させた光を光照射部１１３に伝送する光学部材である。取り回しを良好にする観点からは、バンドルファイバを用いるとよい。他に、プリズム、ミラーなど、光を伝達させるための任意の部材を利用できる。
【００２８】
（光照射部１１３）
光照射部１１３は、光を被検体に照射するための出射端である。光照射部１１３としては、バンドルファイバの終端部を利用できる。また、生体の一部（乳房など）を被検体１００とする場合、パルス光のビーム径を広げて照射するために、光を拡散させる拡散板等を用いてもよい。
【００２９】
（受信部１２０）
受信部１２０は、音響波を受信して電気信号を出力するトランスデューサと、トランスデューサを支持する支持体とを含む。トランスデューサを構成する部材として例えば、圧電材料、静電容量型トランスデューサ（ＣＭＵＴ：Ｃａｐａｃｉｔｉｖｅ  Ｍｉｃｒｏ−ｍａｃｈｉｎｅｄ  Ｕｌｔｒａｓｏｎｉｃ  Ｔｒａｎｓｄｕｃｅｒｓ）、ファブリペロー干渉計を用いたトランスデューサなどを使用できる。圧電材料として例えば、ＰＺＴ（チタン酸ジルコン酸鉛）等の圧電セラミック材料や、ＰＶＤＦ（ポリフッ化ビニリデン）等の高分子圧電膜材料がある。
【００３０】
トランスデューサにより得られる電気信号は時間分解信号である。そのため、電気信号の振幅は、各時刻にトランスデューサで受信される音圧に基づく値（例えば、音圧に比例した値）を表している。なお、トランスデューサとしては、光音響波を構成する周波数成分（典型的には１００ＫＨｚから１００ＭＨｚ）を検出できるものが好ましい。また、支持体に複数のトランスデューサを並べて配置して、１Ｄアレイ、１．５Ｄアレイ、１．７５Ｄアレイ、または２Ｄアレイと呼ばれるような平面や曲面を形成することも好ましい。
【００３１】
受信部１２０が、トランスデューサから出力される時系列のアナログ信号を増幅する増幅器を備えてもよい。また、受信部１２０が、トランスデューサから出力される時系列のアナログ信号を時系列のデジタル信号に変換するＡ／Ｄ変換器を備えてもよい。すなわち、受信部１２０が信号収集部１４０を備えてもよい。
【００３２】
なお、音響波を様々な角度から検出して画像精度を向上させるためには、被検体１００
を全周囲から囲むようなトランスデューサ配置が好ましい。また、全周囲を囲めないほど被検体１００が大きい場合は、半球状の支持体上にトランスデューサを配置してもよい。このような形状の受信部１２０を備えるプローブ１８０は、ハンドヘルド型ではなく、プローブを被検体１００に対して相対移動させる機械走査型の光音響装置に好適である。プローブの移動には、ＸＹステージなどの走査部を用いれば良い。なお、トランスデューサの配置および数、ならびに支持体の形状は、上記に限定されず、被検体１００に応じて最適化すればよい。
【００３３】
受信部１２０と被検体１００との間の空間には、光音響波を伝搬させる媒質を配置すると良い。これにより、被検体１００とトランスデューサの界面における音響インピーダンスが整合する。媒質として例えば、水、油、超音波ジェルなどがある。
光音響装置１は、被検体１００を保持して形状を安定させる保持部材を備えていても良い。保持部材としては光透過性と音響波透過性がともに高いものが好ましい。例えば、ポリメチルペンテンやポリエチレンテレフタレート、アクリルなどを利用できる。
【００３４】
本実施形態に係る装置が、光音響画像に加えて、音響波の送受信により超音波画像も生成する場合、トランスデューサは、音響波を送信する送信手段として機能してもよい。受信手段としてのトランスデューサと送信手段としてのトランスデューサとは、単一（共通）のトランスデューサでもよいし、別々の構成であってもよい。
【００３５】
（光源部２００）
光源部２００は、被検体１００に照射するための光を発生させる装置である。光源部２００としては、大出力のパルス光を発生させ、かつ、酸素飽和度などの物質濃度を取得するために、波長可変な固体レーザー装置が好適である。ただし、半導体レーザー装置や、レーザー以外の光源装置（発光ダイオードやフラッシュランプなど）を利用しても良い。また、波長可変とするために、互いに異なる波長の光を発生させる複数の光源装置を組み合わせて利用しても良い。
【００３６】
光源部２００が発する光のパルス幅は、例えば１ｎｓ以上、１００ｎｓ以下である。また、光の波長としては、４００ｎｍ以上、１６００ｎｍ以下が好適であるが、画像化したい光吸収体の光吸収特性に応じて波長を決定して良い。血管を高解像度でイメージングする場合は、血管での吸収が大きい波長（４００ｎｍ以上、７００ｎｍ以下）を用いてもよい。生体の深部をイメージングする場合には、生体の背景組織（水や脂肪など）において吸収が少ない波長（７００ｎｍ以上、１１００ｎｍ以下）の光を用いてもよい。
【００３７】
本実施形態では、光源部２００として、Ｔｉ：Ｓ（チタンサファイアレーザー）を用いる。Ｔｉ：Ｓの励起にはＮｄ：ＹＡＧレーザー光（波長１０６４ｎｍのナノ秒オーダーのパルス光）を用いる。光源部２００は、被検体の深部まで届く７００ｎｍ以上の２つの波長で発光する。第１の波長λ１は、７９７ｎｍである。第１の波長においては、オキシヘモグロビンとデオキシヘモグロビンの吸収係数が略等しい。なお、波長を選択するときには、構造画像が良好に表示できる波長を用いるとよい。例えば、７７８ｎｍ以上、９５０ｎｍ以下の波長を用いることで、再構成した吸収係数分布画像に対して、画像中の最大吸収係数の３割の値を閾値として２値化処理とトリミングを施した場合に、血管太さの変化が±１０％以内にとどめられる。また、第１の波長と第２の波長という名称は便宜上のものであり、一方の波長と他方の波長のいずれを第１の波長と呼んでも良い
【００３８】
第２の波長λ２は、７５６ｎｍである。この波長には、デオキシヘモグロビンの吸収係数のピークがあるため、オキシヘモグロビンとデオキシヘモグロビンの吸収係数差が大きい。酸素飽和度を計算する際には、この第２の波長７５６ｎｍの光のように、２種類のヘモグロビンの吸収係数の差が大きく、第１の波長７９７ｎｍの光に対する吸収係数と同程
度の吸収係数を持つ波長を選ぶと、精度よく酸素飽和度を取得できる。
【００３９】
（信号収集部１４０）
信号収集部１４０は、受信部１２０から出力されたアナログ信号である電気信号を増幅するアンプと、アンプから出力されたアナログ信号をデジタル信号に変換するＡ／Ｄ変換器とを含む。信号収集部１４０は、ＦＰＧＡ（Ｆｉｅｌｄ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ）チップなどで構成されてもよい。信号収集部１４０から出力されるデジタル信号は、コンピュータ１５０内の記憶部１５２に記憶される。信号収集部１４０は、Ｄａｔａ  Ａｃｑｕｉｓｉｔｉｏｎ  Ｓｙｓｔｅｍ（ＤＡＳ）とも呼ばれる。本明細書において電気信号は、アナログ信号もデジタル信号も含む概念である。なお、信号収集部１４０に、光源部２００からの光を検出する光検出センサを接続することで、光照射と信号収集処理を同期させられる。
【００４０】
上述したように、信号収集部１４０をプローブ１８０のハウジング１８１の内部に配置してもよい。このような構成であれば、プローブ１８０とコンピュータ１５０との間の情報がデジタル信号で伝搬されるため、耐ノイズ性が向上する。また、アナログ信号を伝送する場合に比べ、高速デジタル信号を用いることによって、配線数を少なくすることが可能となり、プローブ１８０の操作性が向上する。
【００４１】
（コンピュータ１５０）
コンピュータ１５０は、演算部１５１、記憶部１５２、制御部１５３、表示制御部１５４を含む。演算部１５１としての演算機能を担うユニットは、ＣＰＵやＧＰＵ（Ｇｒａｐｈｉｃｓ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）等のプロセッサ、ＦＰＧＡ（Ｆｉｅｌｄ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ）チップ等の演算回路で構成できる。これらのユニットは、単一のプロセッサや演算回路から構成されても良いし、複数のプロセッサや演算回路から構成されても良い。演算部１５１は、画像再構成による光音響画像データ（構造画像や機能画像）の生成や、その他各種の演算処理を実行する。演算部１５１は、入力部１７０から、被検体音速や保持部の構成などの各種パラメータ入力を受け付けて、演算に用いてもよい。
【００４２】
演算部１５１が電気信号を３次元のボリュームデータに変換するときの再構成アルゴリズムとしては、タイムドメインでの逆投影法、フーリエドメインでの逆投影法、モデルベース法（繰り返し演算法）など、任意の手法を採用できる。タイムドメインでの逆投影法として、Ｕｎｉｖｅｒｓａｌ  ｂａｃｋ−ｐｒｏｊｅｃｔｉｏｎ（ＵＢＰ）、Ｆｉｌｔｅｒｅｄ  ｂａｃｋ−ｐｒｏｊｅｃｔｉｏｎ（ＦＢＰ）、または整相加算（Ｄｅｌａｙ−ａｎｄ−Ｓｕｍ）などが挙げられる。
【００４３】
２波長を用いる場合、演算部１５１は、画像再構成処理によって、第１の波長の光に由来する光音響信号から第１の初期音圧分布を、第２の波長の光に由来する光音響信号から第２の初期音圧分布を生成する。さらに、第１の初期音圧分布を第１の波長の光の光量分布で補正することによって第１の吸収係数分布を、第２の初期音圧分布を第２の波長の光の光量分布で補正することによって第２の吸収係数分布を取得する。さらに、第１および第２の吸収係数分布から、酸素飽和度分布を取得する。なお、最終的に酸素飽和度分布が得られれば良いので、演算の内容や順序はこれに限られない。
【００４４】
記憶部１５２は、ＲＯＭ（Ｒｅａｄ  ｏｎｌｙ  ｍｅｍｏｒｙ）、磁気ディスクやフラッシュメモリなどの非一時記憶媒体で構成できる。また、記憶部１５２は、ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）などの揮発性の媒体であってもよい。なお、プログラムが格納される記憶媒体は、非一時記憶媒体である。また、記憶部１５２は、１つの記憶媒体から構成されるだけでなく、複数の記憶媒体から構成されていてもよい。
記憶部１５２は、演算部１５１により生成される光音響画像データや、光音響画像データに基づいた表示画像など、各種のデータを保存できる。
【００４５】
制御部１５３は、ＣＰＵなどの演算素子で構成される。制御部１５３は、光音響装置の各構成の動作を制御する。制御部１５３は、入力部１７０からの測定開始などの各種操作による指示信号を受けて、光音響装置の各構成を制御してもよい。また、制御部１５３は、記憶部１５２に格納されたプログラムコードを読み出し、光音響装置の各構成の動作を制御する。
【００４６】
表示制御部１５４は、制御部１５３と共通の、または同様の構成を持つ。表示制御部１５４は表示部１６０に対する画像データ出力や画像の調整などを行う。これにより、プローブの移動と光音響測定に伴い順次、酸素飽和度分布画像が表示される。
【００４７】
コンピュータ１５０は本発明に専用に設計されたワークステーションであってもよい。コンピュータ１５０はまた、汎用的なＰＣやワークステーションを、記憶部１５２に格納されたプログラムの指示に従って動作させたものであっても良い。また、コンピュータ１５０の各構成は異なるハードウェアによって構成されてもよい。また、コンピュータ１５０の少なくとも一部の構成は単一のハードウェアで構成されてもよい。
【００４８】
図３は、本実施形態に係るコンピュータ１５０の具体的な構成例を示す。本実施形態に係るコンピュータ１５０は、ＣＰＵ１５４、ＧＰＵ１５５、ＲＡＭ１５６、ＲＯＭ１５７、外部記憶装置１５８から構成される。また、コンピュータ１５０には、表示部１６０としての液晶ディスプレイ１６１、入力部１７０としてのマウス１７１、キーボード１７２が接続されている。
【００４９】
また、コンピュータ１５０および受信部１２０は、共通の筺体に収められた構成で提供されてもよい。また、筺体に収められたコンピュータで一部の信号処理を行い、残りの信号処理を筺体の外部に設けられたコンピュータで行ってもよい。この場合、筺体の内部および外部に設けられたコンピュータを総称して、本実施形態に係るコンピュータとすることができる。すなわち、コンピュータを構成するハードウェアが一つの筺体に収められていなくてもよい。コンピュータ１５０として、クラウドコンピューティングサービスなどで提供される、遠隔地に設置された情報処理装置を用いても構わない。
コンピュータ１５１は、本発明の処理部に相当する。特に、演算部１５１が中心となって処理部の機能を実現する。
【００５０】
（表示部１６０）
表示部１６０は、液晶ディスプレイや有機ＥＬ（Ｅｌｅｃｔｒｏ  Ｌｕｍｉｎｅｓｃｅｎｃｅ）などのディスプレイである。コンピュータ１５０により得られた被検体情報等に基づく画像や特定位置の数値等を表示する装置である。表示部１６０は、画像や装置を操作するためのＧＵＩを表示してもよい。表示部１６０またはコンピュータ１５０において画像処理（輝度値の調整等）を行ってもよい。
【００５１】
（入力部１７０）
入力部１７０としては、ユーザーが操作可能な、マウスやキーボードなどで構成される操作コンソールを採用できる。また、表示部１６０をタッチパネルで構成し、表示部１６０を入力部１７０として利用してもよい。入力部１７０は、ユーザーからの指示や数値などの入力を受け付け、コンピュータ１５０に伝達する。
【００５２】
なお、光音響装置の各構成はそれぞれ別の装置として構成されてもよいし、一体となった１つの装置として構成されてもよい。また、光音響装置の少なくとも一部の構成が一体
となった１つの装置として構成されてもよい。
【００５３】
また、コンピュータ１５０は、制御部１５３により、光音響装置に含まれる構成の駆動制御も行う。また、表示部１６０は、コンピュータ１５０で生成された画像の他にＧＵＩなどを表示してもよい。入力部１７０は、ユーザーが情報を入力できるように構成されている。ユーザーは、入力部１７０を用いて測定開始や終了、作成画像の保存指示などの操作を行うことができる。
【００５４】
（被検体１００）
被検体１００は光音響装置を構成するものではないが、以下に説明する。本実施形態に係る光音響装置は、人や動物の悪性腫瘍や血管疾患などの診断や化学治療の経過観察などを目的として使用できる。よって、被検体１００としては、生体、具体的には人体や動物の乳房や各臓器、血管網、頭部、頸部、腹部、手指および足指を含む四肢などの診断の対象部位が想定される。例えば、人体が測定対象であれば、オキシヘモグロビンあるいはデオキシヘモグロビンやそれらを含む多く含む血管あるいは腫瘍の近傍に形成される新生血管などを光吸収体の対象としてもよい。また、頸動脈壁のプラークなどを光吸収体の対象としてもよい。また、メチレンブルー（ＭＢ）、インドシニアングリーン（ＩＣＧ）などの色素、金微粒子、またはそれらを集積あるいは化学的に修飾した外部から導入した物質を光吸収体としてもよい。また、穿刺針や穿刺針に付された光吸収体を観察対象としてもよい。被検体は、ファントムや試験対象物などの無生物であっても良い。
【００５５】
（実施形態の動作）
図４は、第１の実施形態における動作を説明するタイミング図である。図４において横軸は時間軸である。下記の様々な制御は、基本的にはコンピュータ１５０が行う。以下の記載における「λ２／ｔ０」などの記載は、「発光される波長／発光タイミング」を示す。
【００５６】
（タイミングＡ〜Ｄ：発光（λ２）・再構成）
光源部２００は時刻ｔ０に第２の波長λ２で発光する（Ａ）。被検体１００からの光音響波（第２の光音響波）は、受信部１２０および信号収集部１４０によりデジタル信号（第２の信号）に変換され、コンピュータ１５０内の記憶部１５２に記憶される（Ｂ）。続いて、コンピュータ１５０は、記憶部１５２に記憶された光音響波に基づくデジタル信号（第２の信号）を用い、３次元のボリュームデータに変換する再構成アルゴリズムにより、初期音圧分布画像（第２の初期音圧分布）を算出する。さらに、初期音圧分布画像に対して、第２の波長λ２の光の光量分布（第２の光量分布）で補正し吸収係数分布画像（第２の吸収係数分布）を得る（Ｃ）。吸収係数分布画像は、次の光源部２００の第２の波長λ２の発光（時刻ｔ２）に伴う吸収係数分布画像取得まで、記憶部１５２に記憶される（Ｄ）。
【００５７】
（タイミングＥ〜Ｈ：発光（λ１）・再構成）
次に、光源部２００は時刻ｔ１に第１の波長λ１で発光する（Ｅ）。被検体１００からの光音響波（第１の光音響波）は、受信部１２０および信号収集部１４０によりデジタル信号（第１の信号）に変換され、コンピュータ１５０内の記憶部１５２に記憶される（Ｆ）。続いて、コンピュータ１５０は、記憶部１５２に記憶された光音響波に基づくデジタル信号（第１の信号）を用い、３次元のボリュームデータに変換する再構成アルゴリズムにより、初期音圧分布画像（第１の初期音圧分布）を算出する。さらに、初期音圧分布画像に対して、第１の波長λ１の光の光量分布（第１の光量分布）で補正し吸収係数分布画像（第１の吸収係数分布）を得る（Ｇ）。吸収係数分布画像は、次の光源部２００の第１の波長λ１の発光（時刻ｔ３）に伴う吸収係数分布画像取得まで、記憶部１５２に記憶される（Ｈ）。
【００５８】
（タイミングＪ〜Ｌ：酸素飽和度分布の生成・表示）
次に、コンピュータ１５０は、第２の波長λ２の吸収係数分布画像（Ｄ）と、第１の波長λ１の吸収係数分布画像（Ｈ）から、機能情報に基づく画像（機能画像）である酸素飽和度分布画像を得る（Ｊ）。このとき、被検体の体動や受信部１２０の移動により、第１の波長の光に由来する光音響波の取得位置と、第２の波長の光に由来する光音響波の取得位置とがずれる場合がある。その位置ずれを補正するために、コンピュータ１５０は、得られた吸収係数分布画像の相関を利用してプローブ移動量を推定し、吸収係数分布画像（Ｈ）と吸収係数分布画像（Ｄ）を位置合わせしてから、精度の良い酸素飽和度分布画像を得ることが好ましい。あるいは、時刻ｔ０から時刻ｔ１におけるプローブの移動量の測定のために、プローブ１８０にジャイロを搭載しても良い。
【００５９】
補正の別手法として、吸収係数分布画像（Ｈ）と吸収係数分布画像（Ｄ）に対して、平滑化（移動平均）フィルタや、ガウシアンフィルタや、メディアンフィルタを用い、ぼかす処理を行っても良い。ぼかし処理を行うことで、吸収係数分布画像（Ｈ）と吸収係数分布画像（Ｄ）の間の位置合わせが上手く行えない場合でも、血管位置の領域内で、機能情報である酸素飽和度を精度良く計算できる。この様な処理を行うため、得られた酸素飽和度分布画像（Ｊ）は、ぼけた画像であることが多い。
【００６０】
一方、コンピュータ１５０は、構造情報に基づく画像（構造画像）である血管画像を、吸収係数分布画像に対する画像処理によって取得する（Ｉ）。画像処理方法として例えば、吸収係数分布画像中の最大吸収係数の３割の値を閾値として２値化を行い、閾値より大きい吸収係数値を示す位置を血管だと判定する方法がある。
【００６１】
コンピュータ１５０は次に、酸素飽和度分布画像（Ｊ）に対して、血管画像（Ｉ）でトリミングを行い、血管画像内の酸素飽和度がわかる表示画像を生成する（Ｋ）。以上の処理を、次の光源部２００の発光（ｔ２）前に完了し、得られた表示画像を表示部１６０に対して出力する（Ｌ）。
【００６２】
（タイミングＭ〜Ｐ：発光（λ２）・再構成）
続いて、時刻ｔ０のときと同様に、光源部２００は時刻ｔ２に第２の波長λ２で発光する（Ｍ）。被検体１００からの光音響波は、受信部１２０および信号収集部１４０によりデジタル信号に変換され、コンピュータ１５０内の記憶部１５２に記憶される（Ｎ）。続いて、コンピュータ１５０は、記憶部１５２に記憶された光音響波に基づくデジタル信号を用い、３次元のボリュームデータに変換する再構成アルゴリズムにより、初期音圧分布画像を算出する。さらに、初期音圧分布画像に対して、第２の波長λ２の光の光量分布で補正し吸収係数分布画像を得る（Ｏ）。吸収係数分布画像は、次の光源部２００の発光（時刻ｔ４）に伴う吸収係数分布画像取得まで、記憶部１５２に記憶される（Ｐ）。
【００６３】
（タイミングＱ〜Ｔ：酸素飽和度分布の生成・表示）
次に、コンピュータ１５０は、第１の波長λ１の吸収係数分布画像（Ｈ）と、第２の波長λ２の吸収係数分布画像（Ｐ）から、機能画像である酸素飽和度分布画像を得る（Ｒ）。前述したように、得られた酸素飽和度分布画像（Ｒ）は、ぼけた画像であることが多い。一方、コンピュータ１５０は、吸収係数分布画像に対する閾値処理等の画像処理ことにより血管画像を得る（Ｑ）。次に、酸素飽和度分布画像（Ｒ）に対して、血管画像（Ｑ）でトリミングを行い、血管画像内の酸素飽和度がわかる表示画像を算出する（Ｓ）。以上の処理を、次の光源部２００の発光（ｔ３）前に完了し、得られた表示画像を表示部１６０に対して出力する（Ｔ）。
【００６４】
上記のように、順次、ある光照射の組に由来する表示画像の生成処理を、光源部２００
の次の発光前に完了し、得られた表示画像を表示部１６０に対して出力することにより、酸素飽和度を示す画像をリアルタイムに表示できる。このとき、第１の波長λ１の吸収係数分布画像（Ｈ）を、表示画像（Ｌ）の生成と表示画像（Ｔ）の生成の両方で用いることにより、２波長の発光が完了するのを待つことなく、発光ごとに表示画像のフレームを生成できる。
【００６５】
上記のタイミングチャートに従った被検体情報取得方法においては、２つの波長（第１および第２の波長）での複数回の発光が行われる。複数回の発光では、どの波長の光が先に照射されても良い。例えば、ｔ０での発光を第１の照射、ｔ１での発光を第２の照射、ｔ２での発光を第３の照射…、と呼んでも良いし、ｔ１での発光を第１の照射、ｔ２での発光を第２の照射、ｔ３での発光を第３の照射…と呼んでも良い。
また、第１の照射により励起された光音響波に由来する信号を第１の信号、第２の照射により励起された光音響波に由来する信号を第２の信号、第３の照射により励起された光音響波に由来する信号を第３の信号、…と呼ぶことができる。
また、第１の信号に基づく吸収係数分布を第１の吸収係数分布、第２の信号に基づく吸収係数分布を第２の吸収係数分布、第３の信号に基づく吸収係数分布を第３の吸収係数分布、…と呼ぶことができる。
また、第１の信号に基づく吸収係数分布から抽出された構造情報を第１の構造情報、第２の信号に基づく吸収係数分布から抽出された構造情報を第２の構造情報、第３の信号に基づく吸収係数分布から抽出された構造情報を第３の構造情報、…と呼ぶことができる。
また、第１の信号に基づく吸収係数分布と第２の信号に基づく吸収係数分布から生成された機能情報を第１の機能情報、第２の信号に基づく吸収係数分布と第３の信号に基づく吸収係数分布から生成された機能情報を第２の機能情報、…と呼ぶことができる。
また、第１の機能情報を第１ないし第３の構造情報の少なくともいずれか（好ましくは、第１または第２の構造情報）に基づいてトリミングした画像を第１の表示画像、第２の機能情報を第１ないし第３の構造情報の少なくともいずれか（好ましくは、第２または第３の構造情報）に基づいてトリミングした画像を第２の表示画像、…と呼ぶことができる。
【００６６】
上記のタイミングチャートに従った被検体情報取得方法においては、従来の、複数の波長の発光が終了した後に表示画像が得られるような手法に比べて、表示画像の生成回数が多くなる。そのため、発光回数に対する表示画像の更新回数が、従来技術の場合と比べて多くなる。そのため、プローブ移動に対する時間追従性が向上する。
例えば、酸素飽和度分布を生成するために用いられた第１の吸収係数分布に関連する第１の波長での発光の回数を「第１の回数」とし、酸素飽和度分布を生成するために用いられた第２の吸収係数分布に関連する第２の波長での発光の回数を「第２の回数」とする。このとき、従来の手法においては、「表示画像の生成回数＝第１の回数＝第２の回数」であった。
一方、上記タイミングチャートの方法においては、あるフレームでの表示画像生成に用いた吸収係数分布を、後のフレームでの表示画像生成にも利用することにより、「表示画像の生成回数＝第１の回数＋第２の回数」となる。この表示画像を順次、表示することにより、スムーズで違和感の少ないリアルタイム表示が実現できる。
【００６７】
なお、上記タイミングチャートの場合は、第１または第２の波長のいずれかで発光がなされるたびに、表示画像が順次、更新される。ただし、発光のたびに表示画像を生成しない場合でも、従来技術の場合よりも多い回数の表示画像生成を行えば、本発明の効果は得られる。この関係は、表示画像の生成回数が以下の式（１）〜式（３）を満たす場合として表現できる。
第１の回数  ＜  生成回数    …（１）
第２の回数  ＜  生成回数    …（２）
生成回数  ≦  第１の回数＋第２の回数    …（３）
【００６８】
（フローチャート）
図５は本発明の第１の実施形態における動作を説明するためのフローチャートである。図４のタイミング図を用いた説明と重複する部分があるが、図５を用いて、コンピュータ１５０が行う動作を説明する。
【００６９】
ステップＳ１００において、光源部２００が波長λ２で発光し、波長λ２の吸収係数分布画像を取得する（図４、Ａ〜Ｄ）。
ステップＳ１０１において、コンピュータ１５０は、変数“ｎ”に“０”を代入する。
ステップＳ１０２において、光源部２００が波長λ１で発光し、波長λ１の吸収係数分布画像を取得する（図４、Ｅ〜Ｈ）。
【００７０】
ステップＳ１０３において、コンピュータ１５０は、波長λ２の吸収係数分布画像をずらす。
ステップＳ１０４において、ずらした波長λ２の吸収係数分布画像とステップＳ１０２で求めた波長λ１の吸収係数分布画像から、酸素飽和度分布画像を算出する（図４、Ｊ）。なお、発光タイミング間の光照射位置のずれを補正する方法は任意である。例えば画像生成位置を波長λ２の照射位置に合わせても良いし、波長λ１とλ２それぞれの照射位置の中間に合わせても良い。いかなる方法であれ、２つの画像における、被検体内部の同じ位置を示す点同士を対応付けられれば良い。
ステップＳ１０５において、コンピュータ１５０は、吸収係数分布画像から血管画像（構造画像）を算出する（図４、Ｉ）。
ステップＳ１０６において、コンピュータ１５０は、血管画像（構造画像）で酸素飽和分布をトリミングし表示画像データを算出する（図４、Ｋ）。
ステップＳ１０７において、表示画像データを表示する（図４、Ｌ）。
以上の処理により、リアルタイム画像表示の１単位分の画像処理が完了する。表示単位は発光間隔によって規定されるものであり、本実施例では、波長を問わず、１回の発光が１単位に相当する。
【００７１】
ステップＳ１０８において、光源部２００が波長λ２で発光し、波長λ２の吸収係数分布画像を取得する（図４、Ｍ〜Ｐ）。
ステップＳ１０９において、コンピュータ１５０は、波長λ１の吸収係数分布画像をずらす。
ステップＳ１１０において、コンピュータ１５０は、ずらした波長λ１の吸収係数分布画像とステップＳ１０８で求めた波長λ２の吸収係数分布画像から酸素飽和度分布画像を算出する（図４、Ｒ）。
ステップＳ１１１において、コンピュータ１５０は、吸収係数分布画像から血管画像（構造画像）を算出する（図４、Ｑ）。
ステップＳ１１２において、コンピュータ１５０は、血管画像（構造画像）で酸素飽和分布をトリミングし表示画像データを算出する（図４、Ｓ）。
ステップＳ１１３において、表示画像データを表示する（図４、Ｔ）。
以上の処理により、リアルタイム画像表示の次の１単位分の画像処理が完了する。
【００７２】
ステップＳ１１４において、コンピュータ１５０は、終了条件が満たされているかどうかを判定する。終了条件として例えば、ユーザーからの入力部を用いた停止指示の検知、接触センサ等によるユーザーがプローブ１８０を手から離したことの検知、測定開始から所定の時間の経過、などが考えられる。終了条件を満たしている場合、光音響測定およびリアルタイム表示は終了する。
一方、終了条件が満たされていない場合、ステップＳ１１５に進む。コンピュータ１５
０はｎのインクリメントを行い、ステップＳ１０２に戻って上述の工程を繰り返す。
上記のように、本フローチャートに従って、光源部２００の発光に合わせた、表示画像データの更新および表示が行われることで、図４で説明したタイミング図の動作が実行される。
【００７３】
本実施形態の光音響装置を用いた表示方法によれば、互いに異なる複数の波長の光を用いて酸素飽和度を取得するときに、光パルスの照射間隔内で構造画像（血管画像）と機能画像を精度良く求めて表示画像を形成し、リアルタイムに表示することが可能である。
【００７４】
＜変形例１＞
図５のフローチャートのステップＳ１０５、ステップＳ１１１の工程では、それぞれ第１，第２の波長に対応する吸収係数分布画像に対して、閾値を設けて２値化処理を行うことで、血管画像（構造情報）を求める。本変形例では、血管幅の誤差を補正して良好な血管画像（構造情報）を求める手法を説明する。
【００７５】
（血管画像の誤差について）
図６を見る前提として、光音響測定対象となる被検体内部には動脈と静脈が走行しており、両者の血管幅は同程度であるとする。図６は、上記被検体について、第１の波長λ１（７９７ｎｍ）に由来する音響波から生成した吸収係数分布を破線で、第２の波長λ２（７５６ｎｍ）に由来する吸収係数分布を実線で、それぞれ示したものである。図中左側は、比較的酸素飽和度が高い動脈に対応し、図中右側は、酸素飽和度が比較的低い（例えば６０％程度）の静脈に対応する。なお、静脈に対応する図中右側では、第１の波長のグラフと第２の波長のグラフが重なっている。
【００７６】
ここで、２種類のヘモグロビンの吸収係数差に起因して、抽出された血管画像中における動脈の血管幅が光の波長によって変化する現象が知られている。
上述したように、第１の波長λ１ではオキシヘモグロビンとデオキシヘモグロビンの吸収係数が略等しい。そのため図６に示されるように、第１の波長に由来する吸収係数分布画像においては、血管が動脈であるか静脈であるかに関わらず、同じ太さの血管は同じ太さで抽出される。よって、第１の波長に対応する吸収係数分布画像から血管を抽出する場合、その血管が動脈であっても静脈であっても、閾値処理（例えば最大吸収係数の３０％を閾値とした２値化処理）によって良好に血管画像（構造画像）を取得できる。
【００７７】
一方、第２の波長λ２ではデオキシヘモグロビンの吸収係数のピークがあるので、オキシヘモグロビンとデオキシヘモグロビンの吸収係数差が大きい。そのため、第２の波長に由来する吸収係数分布画像においては、実際には血管幅が同じであっても、静脈のほうがより太く、動脈のほうがより細く抽出される。
【００７８】
（補正方法）
動画表示であるリアルタイム表示においては、血管幅の誤差がユーザーに与える妨害感は少ないため、血管幅の補正は必要ない、という考え方もある。しかし、より高精度な画像を表示するために補正を行う場合は、次に示す方法が有効である。まず、吸収係数分布画像の注目する領域の近傍において、吸収係数のピークの値を取得する。近傍は、ピクセルやボクセルなどの画素の数に基づき決定しても良いし、距離に基づき決定しても良い。続いて、注目領域の近傍におけるピーク値（近傍ピーク値）を取得し、近傍ピーク値の３０％を第２の閾値として設定し、第２の閾値よりも高い吸収係数値を持つ部分を血管として抽出する。一般的に、近傍ピーク値は画像全体のピーク値と比べて低いため、血管として抽出される画素数が増加し、血管幅が太い動脈が抽出される。
【００７９】
＜変形例２＞
また、波長λ２における構造画像（血管画像）として、波長λ１で得られた構造画像（血管画像）を補正して利用しても良い。補正の際には、プローブの動きや体動を基にした演算を行う。ただし、体動の様な不規則な動きがある場合、波長λ１の構造画像に対して、複雑な処理（例えば、画像を分割して、各々の分割領域の動きを求めて補正する処理）が必要であり、計算量が膨大となる。また動き推定に誤りがある場合、動画が不自然になる。そのため、限られた時間内に表示画像を作成しなくてはならないリアルタイム表示には、波長の異なる光源部２００の光パルスの照射毎に得られた吸収係数分布画像から構造画像（血管画像）を得る方法が、簡便な構成で自然な表示が可能なため好適である。
【００８０】
＜第２の実施形態＞
２種類のヘモグロビンの吸収係数差が大きい波長の吸収係数分布画像に対する閾値処理によって構造画像（血管画像）を求める場合、図６に示したように、酸素飽和度に応じて抽出された血管幅に違いが生じる可能性がある。第２の実施形態では、このような吸収係数分布画像からでも、構造画像（血管画像）を高精度に求める方法を説明する。なお、第１の実施形態と同じ構成や処理については説明を省略する。
【００８１】
図５のステップＳ１０５において、吸収係数分布画像からの血管画像（構造画像）の算出に用いられる波長λ１（７９７ｎｍ）は、２種のヘモグロビンの吸収係数が略等しい波長である。そのため、吸収係数分布画像中のピーク値をもとに閾値を決定して２値化処理を行うことで、構造情報（血管幅）を良好に決定できる。
【００８２】
一方、ステップＳ１１１で血管抽出に用いられる第２の波長λ２（７５６ｎｍ）では、２種のヘモグロビンの間で吸収係数が異なる。そのため、Ｓ１０５と同じ方法で構造情報（血管幅）を取得しようとすると、実際は同じ血管幅の動脈と静脈であっても、抽出された血管幅が異なってしまう。この問題に対し、第１の実施形態の変形例１では、ステップＳ１１１の工程で、注目する領域の近傍におけるピークの値をもとに閾値（第２の閾値）を決定し２値化処理を行い、構造情報（血管幅）を取得した。
【００８３】
第２の実施形態では、構造情報（血管幅）を精度良く算出するために、ステップＳ１１０において求められた機能画像である酸素飽和度分布画像も用いて、ステップＳ１１１で用いる閾値を補正する。具体的には、コンピュータ１５０はまず、ステップＳ１１０で求めた酸素飽和度分布画像から、動脈のような酸素飽和度が高い領域を抽出する。続いて、抽出した領域について、閾値を下げる補正を行う。一方、静脈のような酸素飽和度の低い領域では、閾値の補正はしない。この様子を図７に示す。第２の波長における動脈画像において、閾値が３０％から低減されて２０％程度になっていることが分かる。このように補正された閾値（第３の閾値）を用いることにより、波長によらず正しい構造画像（血管画像）を取得できる。
【００８４】
閾値の補正方法としては、酸素飽和度に比例する定数を元の閾値に乗じて新たな閾値を求める方法や、酸素飽和度に対する補正値をテーブルに記憶し、酸素飽和度で参照した補正値を閾値に乗じて新たな閾値を求める方法などがある。また、酸素飽和度に対する閾値（ピークに対して何％を閾値にする等）自体をテーブルに直接記憶し、酸素飽和度で参照し閾値を決定しても良い。このような方法で、酸素飽和度分布画像を基に適切な閾値を簡便に取得できる。
【００８５】
また、第１の実施形態では、波長λ１として２種類のヘモグロビンの吸収係数差が小さい波長の光源を用いた。しかし、上述した閾値の補正を行うのであれば、波長λ１として必ずしも吸収係数差が小さい波長を使用する必要はない。
【００８６】
第２の実施形態によれば、波長の異なる光源部２００の光パルスの照射の間隔内で、構
造画像と機能画像を精度良く求めて表示画像を形成し、リアルタイムに表示できる。
＜第３の実施形態＞
【００８７】
次に、第３の実施形態を説明する。第１の実施形態では、波長の異なる光パルスの照射ごとに表示画像が更新される。便宜上、このモードを「リアルタイムモード」と呼ぶことにする。第３の実施形態に係る光音響装置は、「リアルタイムモード」の表示が好ましくない場合に、より忠実な表示が行われる「忠実モード」を有し、各モードの切り換えが可能となっている。なお、上記各実施形態と同じ構成や処理については説明を省略する。
【００８８】
図８は、本発明の第３の実施形態における構造画像を高精度に表示できるモード（「忠実モード」）のタイミング図である。図８においても横軸は時間軸である。これらの制御は、コンピュータ１５０が行う。
【００８９】
第１の実施形態同様、初めに、光源部２００が時刻ｔ０に発光した第２の波長λ２の光に基づき、吸収係数分布画像が取得される（Ａ〜Ｃ）。吸収係数分布画像は、次の光源部２００の第２の波長λ２の発光（時刻ｔ２）に伴う吸収係数分布画像取得まで、記憶部１５２に記憶される（Ｄ）。次に、光源部２００が時刻ｔ１に発光した第１の波長λ１の光に基づき、吸収係数分布画像が取得される（Ｅ〜Ｇ）。吸収係数分布画像は、次の光源部２００の発光（時刻ｔ３）に伴う吸収係数分布画像取得まで、記憶部１５２に記憶される（Ｈ）。
【００９０】
続いて、第１の実施形態同様に、第２の波長λ２の吸収係数分布画像（Ｄ）と、第１の波長λ１の吸収係数分布画像（Ｈ）から、機能画像である酸素飽和度分布画像が得られる（Ｊ）。一方、構造画像（血管画像）は、吸収係数分布画像に対する画像処理（例えば、最大吸収係数の３割を閾値とする２値化処理）により取得される（Ｉ）。次に、酸素飽和度分布画像（Ｊ）に対して血管画像（Ｉ）でトリミングを行うことで、血管画像内の酸素飽和度が分かる表示画像が算出される（Ｋ）。そして、得られた表示画像を表示部１６０に対して出力する（Ｌ）。
【００９１】
次に、光源部２００が時刻ｔ２に発光した第２の波長λ２の光に基づき、吸収係数分布画像が取得される（ａ〜ｃ）。吸収係数分布画像は、次の光源部２００の第２の波長λ２の発光（時刻ｔ４）に伴う吸収係数分布画像取得まで、記憶部１５２に記憶される（ｄ）。次に、光源部２００が時刻ｔ３に発光した第１の波長λ１の光に基づき、吸収係数分布画像が取得される（ｅ〜ｇ）。吸収係数分布画像は、次の光源部２００の発光（時刻ｔ５）に伴う吸収係数分布画像取得まで、記憶部１５２に記憶される（ｈ）。
すなわち、本実施形態においては、図４の（Ｓ）に相当する表示画像の生成が行われない。言い換えると、第２の波長λ２に対応する吸収係数分布画像を用いた閾値処理やトリミングが行われない。
【００９２】
続いて、コンピュータ１５０は、第２の波長λ２の吸収係数分布画像（ｄ）と、第１の波長λ１の吸収係数分布画像（ｈ）から、機能画像である酸素飽和度分布画像を得る（ｊ）。一方、コンピュータ１５０は、吸収係数分布画像に対する閾値処理（例えば、最大吸収係数の３割を閾値としたトリミング）により、構造画像（血管画像）を得る（ｉ）。次に、酸素飽和度分布画像（ｊ）に対して、血管画像（ｉ）でトリミングを行い、血管画像内の酸素飽和度がわかる表示画像を算出する（ｋ）。そして、得られた表示画像を表示部１６０に対して出力する（ｌ）。
【００９３】
以上の処理を、順次、光源部２００の波長λ２、λ１の発光の組ごとに行い、得られた表示画像を表示部１６０に対して出力することにより、表示画像をリアルタイムに表示できる。この場合、表示部１６０の表示画像の更新は、波長λ２、λ１の発光の組ごとにな
る。
【００９４】
以上説明した様に、「忠実モード」では、波長間での吸収係数差が無い第１の波長によるトリミングのみが行われるので、構造画像の精度が向上する。したがって、表示画像において血管幅の変動がなくなる。一方、「忠実モード」は「リアルタイムモード」と比べて表示部１６０の表示画像の更新頻度が１／２になり、リフレッシュ周波数が低い表示となる。その結果、プローブ１８０の動きに合わせた表示に、遅延やカクカク感などの妨害感が発生する。
【００９５】
本実施形態の光音響装置は、プローブの移動に対する追従性が高い反面、構造画像の精度が低い「リアルタイムモード」と、構造画像の精度が高いものの追従性が低い「忠実モード」という２つのモードを有する。このような光音響装置におけるモードの指定や切り替えには、様々な方法が考えられる。例えば、ユーザーの好みや被検体の状況に応じて、ユーザーが入力部１７０を介してモードを指定しても良い。また、プローブ１８０の動きに応じて、コンピュータ１５０が自動でモードを切り換えるようにしても良い。具体的には、速度が遅い時（又は止まった時）は「忠実モード」、速度が速い時は「リアルタイムモード」に自動で切り替えるように制御しても良い。速さの検出には、プローブ１８０にジャイロを搭載する方法や、光源部２００の発光間隔と吸収係数分布画像の相関を利用して移動量を算出し、速さを求める方法がある。また、自動でモードが切り換わる場合は特に、ユーザーに分かるように、表示部１６０に現在のモードを文字やアイコン等で表示すると良い。
【００９６】
本発明の第３の実施形態によれば、２つのモードをユーザーの意思により、あるいは自動的に切り替えて利用できる。その結果、ゆっくりしたプローブ１８０の動きでは構造画像がより精度良く表示され、速いプローブ１８０の動きではリフレッシュ周波数が高い動きに追従した表示を行うことが可能となる。
【００９７】
＜その他の実施例＞
（プローブと光源部に関する変形例）
本発明は、プローブ１８０の中に光源部２００を有する構成であっても実現できる。本変形例のプローブ１８０は、内部にレーザーダイオード（ＬＤ）及び発光ダイオード（ＬＥＤ）の光源とドライバ回路からなる光源部２００を有する。酸素飽和度を測定するために、発光波長の異なる複数のＬＤやＬＥＤを組み合わせて光源部２００を構成すると良い。本変形例の光学系１１２は、プローブ１８０内で、ＬＤ及びＬＥＤの光パルスを光照射部１１３に導く、光ファイバやレンズやミラー光学系などの光学部材である。なお、これらの光学部材を介さずに、ＬＤ及びＬＥＤから光パルスを被検体１００に直接照射する構造でもかまわない。
【００９８】
本変形例のような、プローブ１８０に光源部２００を一体化した構成の場合、発熱の問題から、光源の発光間隔を短くすることが難しい。しかし、本発明によれば、光源の発光周波数と表示画像のリフレッシュ周波数を同じにできるため、良好なリアルタイム表示を実現できる。
【００９９】
（照射タイミングに関する変形例）
上記各実施形態では、２波長を交互に発光する形態について説明した。しかし本発明は、交互発光以外の発光方式にも適用できる。例えば、第１の波長λ１での発光を２回繰り返した後、第２の波長λ２での発光を２回繰り返すような発光方式が考えられる。この場合、コンピュータ１５０は、まず、ある発光タイミング（タイミング１）で得られた光音響波に基づいて吸収係数分布画像を生成し、閾値処理により構造情報を求める。続いて、タイミング１の発光に由来する吸収係数分布画像と、タイミング１の直近で、異なる波長
の光が発光されたタイミングを検出する（タイミング２）そして、タイミング１の発光に由来する吸収係数分布画像と、タイミング２の発光に由来する吸収係数分布画像を用いて、機能画像（酸素飽和度分布画像）を求める。そして、先ほど求めた構造情報を用いて表示画像を生成する。
【０１００】
本変形例によれば、２波長の交互照射以外の場合でも、良好なリアルタイム表示を実現できる。
【０１０１】
（波長の数に関する変形例）
上記各実施形態では２種類の波長を用いたが、本発明は３波長以上の光を順次発光する場合にも適用できる。この場合、波長によらず各時刻の発光に伴う吸収係数分布画像から形態情報を求め、各時刻の発光に伴う吸収係数分布画像と、その直近の異なる波長の発光に伴う吸収係数分布画像を用いて、機能画像を求めると良い。必要に応じて、波長によらず各時刻の発光に伴う吸収係数分布画像と、その直近の異なる複数の波長の発光に伴う吸収係数分布画像を用いて機能画像を求めても良い。
【０１０２】
（表示画像に関する変形例）
上記各実施形態では、機能情報を構造情報に基づいて処理して表示画像を生成する手法として、被検体内部の血管の位置を抽出した後、酸素飽和度分布を血管の位置に基づいてトリミングする方法について説明した。しかし、画像処理の手法はこれに限られない。例えば、構造情報に基づく画像と、機能情報に基づく顔図を重畳して表示しても良い。また重畳して表示する際には、構造情報に基づく画像と機能情報に基づく画像を区別して視認性を高めるために、両者を異なる色相で表示しても良い。
【０１０３】
＜その他の実施形態＞
また、本発明は、以下の処理を実行することによっても実現される。即ち、上述した実施形態の機能を実現するソフトウェア（プログラム）を、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータ（またはＣＰＵやＭＰＵ等）がプログラムを読み出して実行する処理である。
【符号の説明】
【０１０４】
１１３：光照射部、１２０：受信部、１５０：コンピュータ、２００：光源部
【図１】
2018126454.tif 000003
【図２】
2018126454.tif 000004
【図３】
2018126454.tif 000005
【図４】
2018126454.tif 000006
【図５】
2018126454.tif 000007
【図６】
2018126454.tif 000008
【図７】
2018126454.tif 000009
【図８】
2018126454.tif 000010