(57)【要約】【課題】精度の高い被検体情報を取得できる被検体情報取得装置を提供する。【解決手段】光を被検体に照射することにより前記被検体で発生する音響波を音響素子によって変換して得られる信号を取得する信号取得手段と、前記取得した信号に対して信号処理を行い、処理後の信号を用いて前記被検体内部の特性情報を取得する情報取得手段と、前記被検体の形態情報を取得する形態取得手段と、を有し、前記情報取得手段は、前記特性情報を取得する対象領域から到来する音響波に対応する信号に対して行う前記信号処理の内容を、前記被検体の形態情報に基づいて決定する。【選択図】図１
2018110734.tif 000002
【特許請求の範囲】
【請求項１】
光を被検体に照射することにより前記被検体で発生する音響波を音響素子によって変換して得られる信号を取得する信号取得手段と、
前記取得した信号に対して信号処理を行い、処理後の信号を用いて前記被検体内部の特性情報を取得する情報取得手段と、
前記被検体の形態情報を取得する形態取得手段と、
を有し、
前記情報取得手段は、前記特性情報を取得する対象領域から到来する音響波に対応する信号に対して行う前記信号処理の内容を、前記被検体の形態情報に基づいて決定する
ことを特徴とする、被検体情報取得装置。
【請求項２】
前記信号処理は、前記信号処理手段が取得した信号を増幅する処理を含み、
前記情報取得手段は、前記被検体の形態情報に基づいて、前記対象領域から到来する音響波に対応する信号から情報が失われない範囲で前記信号の増幅率を決定する
ことを特徴とする、請求項１に記載の被検体情報取得装置。
【請求項３】
前記対象領域は、前記被検体内の表面血管を含んだ領域であり、
前記形態情報は、前記表面血管の位置を表す情報である
ことを特徴とする、請求項１または２に記載の被検体情報取得装置。
【請求項４】
前記形態情報は、前記被検体の外形を表す情報である
ことを特徴とする、請求項１または２に記載の被検体情報取得装置。
【請求項５】
前記対象領域は、前記被検体内の表面血管を含んだ領域であり、
前記情報取得手段は、前記被検体の外形を表す情報に基づいて、前記表面血管の位置を推定する
ことを特徴とする、請求項４に記載の被検体情報取得装置。
【請求項６】
前記形態取得手段は、前記被検体に超音波を送信し、前記被検体にて反射した前記超音波を解析することで前記形態情報を取得する、
ことを特徴とする、請求項１から５のいずれか１項に記載の被検体情報取得装置。
【請求項７】
前記情報取得手段は、前記形態情報に基づいて、前記対象領域から到来する音響波のシミュレーションを行い、前記シミュレーションの結果に基づいて、前記対象領域から到来する音響波に対応する信号を抽出する
ことを特徴とする、請求項１から６のいずれか１項に記載の被検体情報取得装置。
【請求項８】
前記信号処理は、前記音響素子によって変換されたアナログ信号に対して行う、入力インピーダンスマッチング処理、フィルタ処理、増幅処理、減算処理のうちの少なくともいずれかである
ことを特徴とする、請求項１から７のいずれか１項に記載の被検体情報取得装置。
【請求項９】
前記信号処理は、前記音響素子によって変換されたアナログ信号をさらに変換して得られたデジタル信号に対して行う、フィルタ処理、増幅処理、信号抽出処理の少なくともいずれかである
ことを特徴とする、請求項１から８のいずれか１項に記載の被検体情報取得装置。
【請求項１０】
前記音響素子は、グループ分けされた複数の素子からなり、
前記信号取得手段は、前記複数の素子から出力された信号に対して、前記グループごと
に異なる内容の信号処理を行う
ことを特徴とする、請求項１から９のいずれか１項に記載の被検体情報取得装置。
【請求項１１】
装置のシステムノイズを計測するノイズ計測手段をさらに有する
ことを特徴とする、請求項１から１０のいずれか１項に記載の被検体情報取得装置。
【請求項１２】
前記被検体に照射される光を発生させる光源をさらに有する
ことを特徴とする、請求項１から１１のいずれか１項に記載の被検体情報取得装置。
【請求項１３】
一つ以上の前記音響素子を備えた探触子をさらに有する
ことを特徴とする、請求項１から１２のいずれか１項に記載の被検体情報取得装置。
【請求項１４】
前記探触子は、ハンドヘルド型である
ことを特徴とする、請求項１３に記載の被検体情報取得装置。
【請求項１５】
前記探触子は、走査機構によって移動される
ことを特徴とする、請求項１３に記載の被検体情報取得装置。
【請求項１６】
光を被検体に照射することにより前記被検体で発生する音響波を音響素子によって変換して得られる信号を取得する信号取得ステップと、
前記取得した信号に対して信号処理を行い、処理後の信号を用いて前記被検体内部の特性情報を取得する情報取得ステップと、
前記被検体の形態情報を取得する形態取得ステップと、
を含み、
前記情報取得ステップでは、被検体情報を取得する対象領域から到来する音響波に対応する信号に対して行う前記信号処理の内容を、前記被検体の形態情報に基づいて決定する
ことを特徴とする、被検体情報取得方法。
【請求項１７】
請求項１６に記載の被検体情報取得方法の各ステップをコンピュータに実行させるプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被検体内の情報を取得する装置に関する。
【背景技術】
【０００２】
近年、医療分野において、被検体内の構造情報や、生理的情報、すなわち機能情報をイメージングするための研究が進められている。このような技術の一つとして、近年、光音響トモグラフィ（ＰＡＴ：PhotoAcoustic Tomography）が提案されている。
レーザ光などの光を被検体である生体に照射すると、光が被検体内の生体組織で吸収される際に音響波（典型的には超音波）が発生する。この現象を光音響効果と呼び、光音響効果により発生した音響波を光音響波と呼ぶ。被検体を構成する組織は、光エネルギーの吸収率がそれぞれ異なるため、発生する光音響波の音圧も異なったものとなる。ＰＡＴでは、発生した光音響波を探触子で受信し、受信信号を数学的に解析することにより、被検体内の特性情報を取得することができる。
【０００３】
音響波を利用して被検体情報を取得する装置においては、情報の精度を確保するため、適切な信号処理を施すことが好ましい。例えば、特許文献１には、増幅器のゲインを被検体の状況に応じて適切に調整し、Ｓ／Ｎ比を向上させる技術が開示されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１１−２２９８１５号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
特許文献１に記載された被検体情報取得装置では、被検体から発生した光音響波の最大音圧を検出し、音圧が飽和しないように増幅率（ゲイン）を設定している。
【０００６】
ところで、光音響効果を利用した被検体情報取得装置では、血管を対象として測定を行い、得られた結果に基づいて血液の酸素飽和度等を算出することができる。しかし、実際は、生体内部の血管において発生する音響波よりも、皮膚において発生する音響波のほうが強くなるために信号が飽和し、血管に由来する信号が正確に取得できなくなってしまうことがある。この問題を解決するためには、所望する領域（例えば、血管が存在する領域）から発生する信号を飽和させないような処理を行うことが考えられる。しかし、皮膚と表面血管は近接して存在しているため、光音響信号が、どの部位に由来するものであるかを判定することが難しい。
【０００７】
このように、従来の被検体情報取得装置では、皮膚に由来する音響波と、血管に由来する音響波を区別することができないため、取得した光音響信号が飽和しないように、ゲイン設定を低くせざるを得ないという問題があった。また、これにより、血管に由来する信号が十分な精度で得られないという問題があった。
【０００８】
本発明はこのような従来技術の課題に鑑みてなされたものであり、精度の高い被検体情報を取得できる被検体情報取得装置を提供することを目的とする。
【課題を解決するための手段】
【０００９】
本発明に係る被検体情報取得装置は、
光を被検体に照射することにより前記被検体で発生する音響波を音響素子によって変換
して得られる信号を取得する信号取得手段と、前記取得した信号に対して信号処理を行い、処理後の信号を用いて前記被検体内部の特性情報を取得する情報取得手段と、前記被検体の形態情報を取得する形態取得手段と、を有し、前記情報取得手段は、前記特性情報を取得する対象領域から到来する音響波に対応する信号に対して行う前記信号処理の内容を、前記被検体の形態情報に基づいて決定することを特徴とする。
【００１０】
また、本発明に係る被検体情報取得方法は、
光を被検体に照射することにより前記被検体で発生する音響波を音響素子によって変換して得られる信号を取得する信号取得ステップと、前記取得した信号に対して信号処理を行い、処理後の信号を用いて前記被検体内部の特性情報を取得する情報取得ステップと、前記被検体の形態情報を取得する形態取得ステップと、を含み、前記情報取得ステップでは、前記特性情報を取得する対象領域から到来する音響波に対応する信号に対して行う前記信号処理の内容を、前記被検体の形態情報に基づいて決定することを特徴とする。
【発明の効果】
【００１１】
本発明によれば、精度の高い被検体情報を取得できる被検体情報取得装置を提供することができる。
【図面の簡単な説明】
【００１２】
【図１】第一の実施形態に係る被検体情報取得装置の構成を示すブロック図。
【図２】コンピュータ７００の構成を詳細に示した図。
【図３】第一の実施形態に係る被検体情報取得装置の動作フローチャート図。
【図４Ａ】音響波受信素子と被検体との位置関係を説明する図。
【図４Ｂ】音響波受信素子と被検体との位置関係を説明する図。
【図４Ｃ】音響波受信素子と被検体との位置関係を説明する図。
【図５】シミュレーションによって得られた光音響信号の例。
【図６】アナログ回路部７３５の構成を詳細に示した図。
【図７】ＦＩＦＯ７１６に対する書き込みイネーブル信号を説明する図。
【図８Ａ】システムノイズのチェック方法を説明する図。
【図８Ｂ】システムノイズのチェック方法を説明する図。
【図９】第三の実施形態に係る被検体情報取得装置の構成を示すブロック図。
【図１０】コンピュータ１７００の構成を詳細に示した図。
【図１１】第三の実施形態における音響波受信素子の配置図。
【図１２】第三の実施形態に係る被検体情報取得装置の動作フローチャート図。
【図１３】超音波測定と光音響測定の実施タイミングを説明する図。
【発明を実施するための形態】
【００１３】
以下に図面を参照しつつ、本発明の好適な実施の形態について説明する。ただし、以下に記載されている構成部品の寸法、材質、形状およびそれらの相対配置などは、発明が適用される装置の構成や各種条件により適宜変更されるべきものである。よって、この発明の範囲を以下の記載に限定する趣旨のものではない。
【００１４】
本発明は、被検体から伝播する音響波を検出し、被検体内部の特性情報を生成し、取得する技術に関する。よって本発明は、被検体情報取得装置またはその制御方法、あるいは被検体情報取得方法や信号処理方法として捉えられる。本発明はまた、これらの方法をＣＰＵやメモリ等のハードウェア資源を備える情報処理装置に実行させるプログラムや、そのプログラムを格納した記憶媒体や、情報処理装置としても捉えられる。
【００１５】
本発明の被検体情報取得装置には、被検体に光（電磁波）を照射することにより被検体内で発生した音響波を受信して、被検体の特性情報を画像データとして取得する光音響効
果を利用した装置を含む。この場合、特性情報とは、光音響波を受信することにより得られる受信信号を用いて生成される、被検体内の複数位置のそれぞれに対応する特性値の情報である。
【００１６】
光音響測定により取得される特性情報は、光エネルギーの吸収率を反映した値である。例えば、光照射によって生じた音響波の発生源、被検体内の初期音圧、あるいは初期音圧から導かれる光エネルギー吸収密度や吸収係数、組織を構成する物質の濃度を含む。また、物質濃度として酸化ヘモグロビン濃度と還元ヘモグロビン濃度を求めることにより、酸素飽和度分布を算出できる。また、グルコース濃度、コラーゲン濃度、メラニン濃度、脂肪や水の体積分率なども求められる。
【００１７】
本発明に係る被検体情報取得装置には、被検体に超音波を送信し、被検体内部で反射した反射波（エコー波）を受信して、被検体情報を画像データとして取得する超音波エコー技術を利用した装置を含む。この場合、取得される被検体情報とは、被検体内部の組織の音響インピーダンスの違いを反映した情報である。
【００１８】
被検体内の各位置の特性情報に基づいて、二次元または三次元の特性情報分布が得られる。分布データは画像データとして生成され得る。特性情報は、数値データとしてではなく、被検体内の各位置の分布情報として求めてもよい。すなわち、初期音圧分布、エネルギー吸収密度分布、吸収係数分布や酸素飽和度分布などの分布情報である。
【００１９】
本明細書における音響波とは、典型的には超音波であり、音波、音響波と呼ばれる弾性波を含む。探触子等により音響波から変換された電気信号を音響信号とも呼ぶ。ただし、本明細書における超音波または音響波という記載には、それらの弾性波の波長を限定する意図はない。光音響効果により発生した音響波は、光音響波または光超音波と呼ばれる。光音響波に由来する電気信号を光音響信号とも呼ぶ。超音波エコーに由来する電気信号を超音波信号とも呼ぶ。
【００２０】
（第一の実施形態）
<システム構成>
図１は、第一の実施形態に係る被検体情報取得装置の構成を示す機能ブロック図である。
第一の実施形態に係る被検体情報取得装置は光音響装置であり、光源１００、光検出器１２０、光学系２００、音響波受信素子３００、支持体４００、走査部５００、走査位置センサ５１０を有している。
さらに、本実施形態における被検体情報取得装置は、形態取得部６００、コンピュータ７００、給排液ユニット８１０、表示装置９００、入力部１０００、形態保持部１１００、取り付け部１２００を有している。
【００２１】
<<被検体Ｅ>>
被検体Ｅは、本発明における被検体情報取得装置の一部を構成するものではないが、以下に説明する。本発明に係る被検体情報取得装置は、人や動物の悪性腫瘍や血管疾患、血糖値などの診断や化学治療の経過観察を行う装置である。よって、被検体としては、生体、具体的には人体や動物の乳房や指、手足などが想定される。
被検体の内部には、生体であれば水や脂肪、タンパク質、酸化ヘモグロビン、還元ヘモグロビンなど、光吸収係数の大きい光吸収体が存在するため、光照射に起因して光音響波が発生する。被検体としてファントムを用いる場合、光学特性を模擬した物質を光吸収体として内部に封入することで、光音響波の発生および測定が可能となる。なお、便宜上、図１において被検体Ｅは点線で示している。
【００２２】
<<光源１００>>
光源１００は、被検体に照射するパルス光を発生させる装置である。光源は、大出力を得るためレーザ光源であることが望ましいが、レーザの代わりに発光ダイオードやフラッシュランプ等を用いることもできる。光源としてレーザを用いる場合、固体レーザ、ガスレーザ、色素レーザ、半導体レーザなど様々なものが使用できる。
また、パルス光の波長は、被検体を構成する成分のうち特定の成分に吸収される特定の波長であって、被検体内部まで光が伝搬する波長であることが望ましい。具体的には、被検体が生体である場合、７００ｎｍ以上１２００ｎｍ以下であることが望ましい。この領域の光は比較的生体深部まで到達することができ、深部の情報を得ることができる。生体表面部の測定に限定すれば、５００〜７００ｎｍ程度の可視光から近赤外領域を使用してもよい。
また、光音響波を効果的に発生させるためには、被検体の熱特性に応じて十分短い時間に光を照射させなければならない。被検体が生体である場合、光源から発生するパルス光のパルス幅は数１０ナノ秒以下が好適である。
なお、光照射のタイミング、波形、強度等は、後述するコンピュータ７００によって制御される。
【００２３】
<<光学系２００>>
光学系２００は、光源１００から発せられたパルス光を伝送する部材である。光源から出射された光は、レンズやミラーなどの光学部品により、所定の光分布形状に加工されながら被検体に導かれ、照射される。なお、光ファイバなどの光導波路などを用いて光を伝搬させることも可能である。
光学系２００は、例えば、レンズ、ミラー、プリズム、光ファイバ、拡散板、シャッターなどの光学機器などを含んでいてもよい。このような光学部品は、光源から発せられた光が被検体に所望の形状で照射されれば、どのようなものを用いてもかまわない。なお、光はレンズで集光させるより、ある程度の面積に広げる方が、生体への安全性ならびに診断領域を広げられるという観点で好ましい。
本実施形態では、光学系２００は、半球状の支持体４００（後述）に配置された音響波探触子の曲率中心となる領域を照射するように構成されている。
なお、光学系２００におけるシャッター等の開閉は、コンピュータ７００（後述）によって行われる。また、光源１００および光学系２００は、後述する支持体４００の底部に設けられており、支持体４００の動きにあわせて移動する。
【００２４】
<<光検出器１２０>>
光検出器１２０は、光源１００から光が照射されたことを検出する装置である。本実施形態では、光学系２００が有する光ファイバ（不図示）を一部分岐してフォトダイオードに入力し、検出信号をコンピュータ７００に伝達する。なお、光源１００が発光したタイミングを検出できれば、光検出器１２０にはどのようなものを用いてもよい。また、検出手段もフォトダイオードに限定されない。
【００２５】
<<音響波受信素子３００>>
音響波検出素子３００は、被検部Ｅの内部から到来する音響波を受信して、電気信号に変換する手段である。音響波検出素子は、探触子、音響波探触子、音響波検出器、音響波受信器、トランスデューサとも呼ばれる。
生体から発生する音響波は、１００ＫＨｚから１００ＭＨｚの超音波であるため、音響波検出素子には、上記の周波数帯を受信できる素子を用いる。具体的には、圧電現象を用いたトランスデューサ、光の共振を用いたトランスデューサ、容量の変化を用いたトランスデューサなどを用いることができる。
【００２６】
また、音響素子は、感度が高く、周波数帯域が広いものを用いることが望ましい。具体
的にはＰＺＴ（チタン酸ジルコン酸鉛）などを用いた圧電素子、ＰＶＤＦ（ポリフッ化ビニリデン）などの高分子圧電膜材料、ＣＭＵＴ（容量性マイクロマシン超音波トランスデューサ）、ファブリペロー干渉計を用いたものなどが挙げられる。ただし、ここに挙げたものだけに限定されず、探触子としての機能を満たすものであれば、どのようなものであってもよい。
【００２７】
音響波受信素子３００は、受信面の法線方向から光音響波が入射する場合において、受信感度が最も高くなり、受信面の法線方向と光音響波の入射角とがなす角度が大きくなるほど受信感度が低くなる。なお、本実施形態に係る音響波受信素子３００は、円形の平面形状の受信面を有しているものとする。
【００２８】
<<支持体４００>>
支持体４００は、光学系２００および複数の音響波受信素子３００を支持する、略半球形状の容器である。本実施形態では、半球の内側の面に、複数の音響波受信素子３００が設置され、半球の底部（極）に光学系２００が設置されている。また、半球の内側には、後述する音響マッチング材８００が貯留される。支持体４００は、これらの部材を支持するため、機械的強度が高い金属材料などを用いて構成することが好ましい。
【００２９】
複数の音響波受信素子３００は、素子の受信方向が半球の曲率中心に向かうように、半球面上にアレイ状に配置されている。なお、図１は、半球状の支持体４００を中心軸で切断した断面を表す。また、被検体Ｅ内の一部の領域に集束する一点鎖線が、音響波受信素子３００の受信方向を示している。
このように、複数の音響波受信素子３００のそれぞれは、特定の領域で発生する光音響波を高感度に受信することができるように支持体４００上に配置されている。本実施形態においては、この特定の領域（図１における領域Ｇ）を「高感度領域」と呼ぶ。
【００３０】
複数の音響波受信素子３００をこのように配置して被検体情報を取得した場合、半球の曲率中心の分解能が最も高くなり、中心から離れるにつれて分解能が低くなる。本実施形態では、最も高い分解能が得られる点から、当該分解能の半分の分解能が得られる点までの領域を高感度領域とする。
なお、所定の高感度領域を形成できれば、各音響波受信素子の最も感度の高い方向は必ずしも交わらなくてもよい。すなわち、複数の音響波受信素子３００の少なくとも一部の素子が、所定の領域で発生する光音響波を感度よく受信することができるように配置されていればよい。
【００３１】
<<走査部５００>>
走査部５００は、支持体４００の位置を三次元方向に移動することにより、被検体Ｅに対する支持体４００の相対位置を変更する装置（走査機構）である。走査部５００は、ガイド機構と、駆動機構と、走査位置センサ５１０と、をＸ，Ｙ，Ｚ軸の三方向それぞれに有している。
図１に示すように、走査部５００の上に支持体４００が積載されるため、ガイド機構は大きな荷重に耐えることが可能なリニアガイドなどであることが好ましい。また、駆動機構として、リードスクリュー機構、リンク機構、ギア機構、油圧機構などを用いることができる。また、駆動力はモーターなどによって得ることができる。
【００３２】
なお、支持体４００を固定し、被検体Ｅを移動させてもよい。被検体Ｅを移動させる場合は、被検体Ｅを支持する支持部（不図示）、もしくは取り付け部１２００を動かすことで被検体Ｅを移動させる構成が考えられる。さらに、被検体Ｅと支持体４００の両方を移動可能としてもよい。
【００３３】
また、移動は連続的に行うことが望ましいが、一定のステップで繰り返してもよい。また、走査部５００は電動ステージであることが望ましいが、手動ステージであってもよい。ただし、ここに挙げたものだけに限定されず、被検体Ｅと支持体４００のうち少なくとも一方を移動させることができれば、どのようなものであってもよい。
【００３４】
<<走査位置センサ５１０>>
走査位置センサ５１０は、支持体４００の位置（座標）を取得する手段である。走査位置センサ５１０は、装置の構成に応じて、１次元、２次元、もしくは３次元の位置座標情報を取得する。走査位置センサ５１０は、支持体４００の位置を取得することができれば、どのようなセンサであってもよい。例えば、リニアスケール、磁気センサ、赤外線センサ、超音波センサ、エンコーダー、可変抵抗器などを用いたポテンショメータなどを用いることができる。
【００３５】
<<形態取得部６００>>
形態取得部６００は、被検体Ｅの形態情報を取得する手段である。被検体の形態情報とは、被検体を構成する構造物の位置を表す情報である。被検体の形態情報は、被検体表面の形状を表す情報（以下、形状情報）であってもよいし、被検体内部の構造を表す情報（以下、構造情報）であってもよい。
形態取得部６００が取得した形態情報は、後述するコンピュータ７００（演算部７１０）に送信される。
【００３６】
ここで、形態取得部６００が形態情報を取得する方法として、以下の三種類を例示する。
第一の方法は、被検体表面の形状情報を予め記憶しておき、読み出して利用する方法である。例えば、複数の被検体の形状についての情報を、形態取得部６００（あるいはコンピュータ７００）に格納しておき、適切なものを自動または手動で読み出してもよい。また、形態保持部１１００（後述）の形状を表す情報を、被検体表面の形状を表す情報として扱ってもよい。
【００３７】
第二の方法は、被検体Ｅの表面形状を測定して形状情報を取得する方法である。
例えば、形態取得部６００が被検体Ｅに超音波を送信し、反射波を解析することで、被検体Ｅの形状情報を取得することができる。この場合、形態取得部６００が、音響波受信素子３００とは別に、被検体に対して超音波を送受信するトランスデューサを有していてもよい。また、音響波受信素子３００のうちの一部を、超音波を送信可能に構成したうえで、光音響波の受信と、形状を測定するための超音波の送受信とを行うようにしてもよい。なお、反射波を用いた被検体形状の測定は、形態取得部６００が行ってもよいが、コンピュータ７００（演算部７１０）が行うようにしてもよい。
【００３８】
第三の方法は、被検体Ｅの内部に存在する血管の位置（血管形状）を測定によって取得する方法である。
例えば、前述した構成を用いて、超音波測定（ドップラーエコー法）によって、被検体の内部において血流が存在する部分を推定することができる。なお、血流が存在する部分は、２次元情報として取得してもよいし、３次元情報として取得してもよい。さらには、血流の変化の状態をリアルタイムに取得するようにしてもよい。
なお、ドップラーエコー法を用いる場合、別系統の処理回路を設けてもよい。例えば、ダイナミックレンジが広く、分解能が高いＡ／Ｄコンバータを専用に設けてもよい。
【００３９】
なお、第二および第三の方法では、超音波測定を利用して形態情報を取得する例を挙げたが、これ以外の手段（例えばＭＲＩ等）を用いてもよい。また、造影剤と、ＭＲＩやＸ線ＣＴ、その他、血流が存在する部分を確認できるモダリティを併用し、被検体内部の構
造情報を取得してもよい。
また、カメラなどを利用して形態情報を取得してもよい。例えば、複数の方向から撮像した画像に基づいて、ステレオ法などの三次元計測技術を用いて被検体の形態情報を取得することもできる。
【００４０】
また、形態取得部６００は、取得した形態情報を装置の利用者に提示してもよい。例えば、取得した形態情報を、カメラで撮影した画像や、装置が生成した被検体情報（機能情報）などと併せて表示（並列表示、重畳表示、交互表示など）してもよい。
また、形態取得部６００は、複数種類の情報を取得してもよい。例えば、形態保持部１１００の形状情報、被検体表面の形状情報、被検体内部の構造情報（血管形状に関する情報など）などを取得してもよい。また、取得した複数の情報から、利用する情報を選択してもよい。
【００４１】
<<コンピュータ７００>>
コンピュータ７００は、本実施形態に係る被検体情報取得装置の制御を司る手段であり信号取得手段としての機能も備える。
図２は、コンピュータ７００の構成をより詳細に示した図である。図２に示したように、コンピュータ７００は、演算部７１０および記憶部７２０を有しており、演算部７１０は、データ取得部７１８と演算処理部７２８から構成される。さらに、データ取得部７１８は、アナログ回路部７３５とデジタル回路部７３６から構成される。
【００４２】
演算部７１０は、複数の音響波受信素子３００から取得した信号を解析することで、被検体内部の情報を取得する手段である。また、演算部７１０は、形態取得部６００が取得した形態情報に基づいて、複数の音響波受信素子３００に到達する光音響波の波形をシミュレーションし、当該シミュレーション結果を用いて、アナログおよびデジタル信号処理の内容を決定する機能を有している。これにより、被検体内の所望の領域から到来する光音響波を最適な方法で処理することができる。具体的な方法については後述する。
【００４３】
次に、演算部７１０を構成する各手段について説明する。
アナログ回路部７３５は、複数の音響波受信素子３００から出力された電気信号にアナログ信号処理を施し、処理後のアナログ信号をデジタル回路部７３６へ出力する回路である。信号処理の内容として、例えば、入力インピーダンス回路による入力インピーダンスマッチング処理、フィルタによるフィルタ処理、アンプによる増幅処理などのアナログ回路による信号処理などがある。ただし、アナログ信号処理の内容は、例示したものに限定されない。
デジタル回路部７３６は、アナログ回路部７３５から入力されたアナログ信号をサンプリングし、デジタル信号の生成を行う回路である。また、デジタル信号に対し、デジタルゲインの適用、さらにデジタルフィルタの適用といった、デジタル信号処理を行うことができる。ただし、デジタル信号処理の内容は、例示したものに限定されない。
【００４４】
演算処理部７２８は、上位演算器７１５とシミュレーション部７５０から構成される。
上位演算器７１５は、データ取得部７１８において生成されたデジタル信号に対して、画像再構成アルゴリズムに基づく処理を施し、被検体情報を取得する手段である。上位演算器７１５によって、被検体の表面や内部の機能情報、被検体の内部血管における酸素飽和度、その他血中成分の濃度といった情報が取得される。
上位演算器７１５は、典型的にはＣＰＵ、ＧＰＵ、ＦＰＧＡ、ＡＳＩＣなどの回路から構成される。なお、上位演算器７１５は、複数の回路から構成されていてもよい。また、コンピュータ７００が行う各処理をいずれの回路が実行してもよい。また、デジタル回路部７３６は行うデジタル信号処理の一部または全部を、演算処理部７２８が行うようにしてもよい。
【００４５】
シミュレーション部７５０は、形態取得部６００が取得した被検体の形態情報を用いて、光が照射された被検体から発生し、複数の音響波受信素子３００に到達する光音響波の波形をシミュレーションする手段である。また、当該音響波の波形をフーリエ解析し、周波数成分の特定を行う機能を有している。
シミュレーション部７５０が生成したシミュレーション結果は、アナログ回路部７３５およびデジタル回路部７３６へフィードバックされ、各回路において施される信号処理の用に供される。これにより、シミュレーション結果に応じて適切な内容の信号処理を行うことができる。
なお、本実施形態では、シミュレーション部７５０と上位演算器７１５が分かれているが、一つのモジュールとしてもよい。
【００４６】
記憶部７２０は、典型的にはＲＯＭ、ＲＡＭ、およびハードディスクなどの記憶媒体から構成される記憶手段である。なお、記憶部は、複数の記憶媒体から構成されていてもよい。
【００４７】
以上説明したように、演算部７１０は、複数の音響波受信素子３００から出力された電気信号に対して、アナログ信号処理、デジタル信号処理の双方を施したうえで、被検体情報を生成する。
また、演算部７１０は、バス２０００を介して、被検体情報取得装置を構成する各要素の動作を制御する役割も果たす。例えば、演算部７１０は、光源１００がパルス光を発生させるタイミングを制御する。
【００４８】
なお、コンピュータ７００は、同時に複数の信号をパイプライン処理できるように構成されていることが好ましい。これにより、被検体情報を取得するまでの時間を短縮することができる。また、コンピュータ７００が行うそれぞれの処理を、演算部７１０に実行させるプログラムとして記憶部７２０に保存しておくこともできる。この場合、プログラムが保存される記憶部７２０は、非一時的な記録媒体である。
【００４９】
<<音響マッチング材８００>>
音響マッチング材８００は、被検体Ｅと音響波受信素子３００との間の空間を満たし、被検体Ｅと音響波受信素子３００とを音響的に結合させるための部材である。音響マッチング材８００は、音響波受信素子３００と形態保持部１１００との間や、形態保持部１１００と被検体との間に配置してもよい。なお、音響マッチング材８００が複数の箇所に分散する場合、それぞれに異なる部材を用いてもよい。
音響マッチング材８００は、被検体や音響波受信素子３００と音響インピーダンスが近い（例えば、両者の中間）材料であることが好ましい。また、音響マッチング材８００は、光源１００で発生するパルス光を透過する材料であることが好ましい。また、音響マッチング材８００は液体であることが好ましい。典型的には、水、ひまし油、ジェルなどを用いることができる。
【００５０】
<<給排液ユニット８１０>>
給排液ユニット８１０は、音響マッチング材８００を供給および排出するユニットである。給排液ユニット８１０は、測定を実施する前の適切なタイミングで、支持体４００の内部に音響マッチング材８００を供給する。そして、不要となったタイミングで、支持体４００から音響マッチング材８００を回収する。なお、支持体４００に音響マッチング材８００を貯めたままにした場合、室内の湿度が上昇し装置が不具合を起こす可能性があるため、装置不使用時は回収することが望ましい。
【００５１】
<<表示装置９００>>
表示装置９００は、コンピュータ７００から出力される被検体情報を、画像や数値などを用いて表示する装置である。表示装置９００は、典型的には液晶ディスプレイなどであるが、プラズマディスプレイ、有機ＥＬディスプレイ、ＦＥＤなど、他の方式のディスプレイであってもよい。また、表示装置９００は、本実施形態に係る被検体情報取得装置とは別に提供されてもよい。
【００５２】
<<入力部１０００>>
入力部１０００は、装置の利用者がコンピュータ７００に対して情報を入力するための手段である。入力部１０００として、例えば、キーボード、マウス、タッチパネル、ダイヤル、およびボタンなどを用いることができる。また、入力部１０００としてタッチパネルを採用し、表示装置９００と入力部１０００を兼ねることもできる。
【００５３】
<<形態保持部１１００>>
形態保持部１１００は、被検体の形状を一定に保つために被検体と当接される部材である。
形態保持部１１００は、取り付け部１２００に取り付けられている。なお、被検体の形状や保持状態にあわせて、複数の形態保持部を交換可能にしてもよい。この場合、取り付け部１２００は、複数の形態保持部を取り付け可能に構成されていることが好ましい。
本実施形態では、形態保持部１１００を介して被検体Ｅに光を照射するため、形態保持部１１００は、光を透過する部材であることが好ましい。例えば、ポリメチルペンテンやポリエチレンテレフタラートなどを好適に用いることができる。
なお、被検体が乳房である場合、乳房形状の変形を少なくして形状を一定にした状態で保持することが好ましい。このため、形態保持部１１００の形状を、球をある断面で切った形状にするとよい。形態保持部１１００の形状は、被検体の体積や保持後の所望の形状に応じて適宜設計することができる。
【００５４】
<処理フローチャート>
次に、図３に示すフローチャートを参照しながら、本実施形態に係る被検体情報取得装置の動作について説明する。
【００５５】
ステップＳ１００にて、測定を開始する。ここでは、形態保持部１１００に被検体を挿入し、支持体４００と形態保持部１１００との間、および形態保持部１１００と被検体との間に音響マッチング材８００を充填する。
【００５６】
次に、ステップＳ２００にて、形態取得部６００が、被検体の形態情報を取得するための超音波測定（超音波エコー測定）を実施する。本ステップでは、１次元のＡモード撮像、または２次元もしくは３次元のＢモード撮像、またはドップラー撮像を行って、被検体についての形状情報（例えば、被検体の外形に関する情報）や、構造情報（例えば、被検体内の血管形状に関する情報）を得る。本ステップを行うことで、被検体における皮膚と表面血管との位置関係を把握することができる。
なお、被検体の形状情報を取得したうえで、任意の手法によって血管の位置を推定する処理を行ってもよい。
【００５７】
次に、ステップＳ３００にて、演算部７１０（シミュレーション部７５０）が、形態取得部６００が取得した情報に基づいて、光が照射された被検体から発生し、複数の音響波受信素子３００に到達する音響波の波形をシミュレーションする。
【００５８】
ここで、シミュレーションを行う具体的な方法について説明する。
まず、シミュレーションの最初のステップとして、被検体の表面に照射される光の分布（照射光分布）、被検体形状、被検体の光学係数分布から、被検体内部に到達する光量の
分布を計算する。被検体内部の光量分布は、基本的には、式（１）のような拡散方程式を用いて計算することができる。なお、用いる拡散方程式は必ずしも例示したものでなくてもよい。例えば、三次元の光の分布を求めるためや、演算量を適切な量にするため、適切な近似等を用いて、式（１）を適宜変形してもよい。
【数１】
2018110734.tif 000003
【００５９】
次に、前ステップで算出した被検体内部の光量分布に、被検体の吸収係数分布とグリ
ュナイゼン係数を乗算し、被検体内部の初期音圧分布を算出する。
ここで、光音響効果により生じる音響波の音圧と、被検体内の光強度の関係について説明する。光音響効果により生じる音響波の音圧ｐ０［Ｐａ］は、式（２）によって表される。
【数２】
2018110734.tif 000004
【００６０】
式（２）において、μaは光吸収体（腫瘍等）の吸収係数［／mm］である。Γはグリ
ュナイゼン係数である。Φは光吸収体の位置の光強度［J／mm２］である。グリュナイゼ
ン係数Γは、体積膨張係数に音速の２乗をかけた値を定圧比熱で除した値であって、生体では略一定値をとる。式（２）からわかるように、発生する音圧は光の強度に比例する。ここで、式（２）で求められたｐ０の被検体内部における分布が、被検体内部の初期音圧分布ｐ0（ｒ）となる。
【００６１】
次に、被検体内部で発生した光音響波が、どのように生体内部を伝搬するか計算する。一般的に、生体などの非粘性媒体を伝搬する光音響波の伝搬は、以下の光音響波動方程式によって表される。
【数３】
2018110734.tif 000005
【００６２】
式（３）において、ｐ（ｒ，ｔ）は音圧、ρ（ｒ）は密度、c（ｒ）は音速 、Γはグリューナイセン係数である。また、Ｈ（ｒ，ｔ）は単位時間および単位体積当たりの熱量、ｔは時間、ｒは位置である。単位時間および単位体積当たりの熱量Ｈ（ｒ，ｔ）は、熱閉じ込め条件のもとでは、単位体積当たりの熱量Ｈ（ｒ）と、光パルス関数Ｉ（ｔ）に分離することができる。
【数４】
2018110734.tif 000006
【００６３】
また、単位体積当たりの熱量Ｈ（ｒ）と、発生する初期音圧分布p0（ｒ）は、以下の関係にある。
【数５】
2018110734.tif 000007
【００６４】
ここで、μa（ｒ）は光吸収体（腫瘍等）の吸収係数［／mm］であり、Φ（ｒ）は光吸
収体における光強度［J／mm２］である。
【００６５】
以上の結果をまとめると、熱閉じ込め条件が成り立つ条件での波動方程式は以下のようになる。
【数６】
2018110734.tif 000008
【００６６】
この式を解くことで、光音響波の伝搬過程を再現できる。つまり、音響波受信素子３００の存在する位置Ｒにおける、時間Ｔでの光音響波の音圧p（Ｒ，Ｔ）をシミュレートす
ることができる。
【００６７】
なお、シミュレーションに用いる波動方程式は、必ずしも式（６）に限定されない。演算量を適切な量にするために条件を限定する、もしくは適切な近似等を用いることで、式（３）を適宜変形、簡略化させたものを用いて光音響波の伝搬過程を計算してもよい。
【００６８】
また、シミュレーションを実施する対象領域は、被検体の測定領域全てであってもよいし、関心領域が存在する領域に絞ってもよい。本発明の主目的は、皮膚に由来する光音響信号と、表面血管などの関心領域に由来する光音響信号とを区別することであるため、関心領域に絞ってシミュレーションを行うことで、処理時間を短縮できる。
【００６９】
なお、シミュレーションにおいて必要な処理を、全てシミュレーション部７５０が行う必要はない。例えば、別に設けられた外部の計算機（不図示）でシミュレーション（一部または全部）を実施し、結果を取り込むようにしてもよい。このように、処理を分散することで、シミュレーション部７５０にかかる計算負荷を削減することができる。
【００７０】
また、本例では演算によって音響波のシミュレーションを行ったが、サンプルデータに基づいて近似を行ってもよい。例えば、生体内部の代表的な形状から発生する信号波形のサンプルデータを用いて、シミュレーション結果を得るようにしてもよい。かかる手法によると、演算量を削減し、処理の高速化を図ることができる。なお、サンプルデータは、装置の内部に保存しておいてもよいし、装置の外部から取得するようにしてもよい。
【００７１】
次に、図４Ａを参照しながら、シミュレーションを行う対象領域について説明する。
前述したように、シミュレーションは、測定領域全体に対して行ってもよいし、表面血管が存在すると認められる領域（例えば、領域１５０１）に対してのみ行ってもよい。
シミュレーションは、複数の音響波受信素子３００の帯域や指向性を考慮し、表面血管から発生して複数の音響波受信素子３００に到達する光音響波の波形を演算することで行う。なお、図４Ａにおいては、音響波受信素子３００−６と被検体、表面血管が存在する領域１５０１との位置関係を示しているが、このように、素子と関心領域との位置関係を考慮して、シミュレーションを行う領域を決定してもよい。
また、同時に、被検体の表面（すなわち皮膚）において発生し、複数の音響波受信素子３００に到達する光音響信号の波形をシミュレーションしてもよい。このとき、被検体の肌の色情報をパラメータとして用いてもよい。
【００７２】
このようにすることで、表面血管や皮膚から発生し、複数の音響波受信素子３００に到達する音響波の波形や強度、周波数特性を予測することができる。すなわち、表面血管に由来する音響波が複数の音響波受信素子３００に到達するタイミングと、当該音響波の波形、強度が把握できるため、表面血管に由来する光音響信号を飽和させずに精度良く抽出するための条件を正確に把握できる。
シミュレーション結果の利用方法については、ステップＳ７００において説明する。
【００７３】
次に、ステップＳ４００にて利用者が測定条件を入力する。具体的には、例えば、被検体情報を取得する対象領域、光の照射間隔、光の波長、光の照射回数、支持体４００の移動速度、吸収係数、散乱係数、被検体内で発生した光音響波を取得する期間（以下、受信期間）などのパラメータが入力される。
なお、本ステップでは、シミュレーションの結果を参照しつつ、フィルタの帯域特性、入力インピーダンス値、ゲイン値、サンプリングクロック周波数などを入力するようにしてもよい。また、測定条件を表すパラメータは手動で入力してもよいし、シミュレーションの結果に基づいて自動的に入力されるようにしてもよい。
【００７４】
次に、ステップＳ５００にて、走査を開始する。
走査が開始されると、光源１００から所定のタイミングで光が出射するように、コンピュータ７００が制御信号を出力する。出射した光は光学系２００によって導かれ、音響マッチング材８００を介して被検体Ｅに照射される（ステップＳ６００）。そして、被検体Ｅに照射された光が被検体Ｅ内で吸収され、光音響波が発生する。
【００７５】
ステップＳ７００では、コンピュータ７００内部に存在するデータ取得部７１８が、利用者が指定した受信期間において光音響信号を取得する。また、光の照射回数が複数回指定されている場合、当該回数分の光音響信号を取得する。
【００７６】
<<シミュレーション結果に基づく信号処理>>
ここで、シミュレーション結果に基づいた信号処理の内容について説明する。
音響波受信素子３００が受信した音響波を変換して得られるアナログ信号は、データ取得部７１８において処理され、デジタル信号に変換される。
アナログ信号に対して処理を行うアナログ回路部７３５は、例えば、以下のような回路を有している。
・音響波受信素子３００とのインピーダンス整合を取り、アナログ信号の反射に起因する波形劣化を防ぐ入力インピーダンス回路
・アナログ信号の高周波成分やＤＣ成分を除去するフィルタ回路
・アナログ信号を増幅するアンプ回路
・アナログ信号をデジタル信号に変換するＡ／Ｄコンバータ回路
音響波受信素子３００から出力されたアナログ信号の電圧値は、受信された光音響信号の音圧の時系列変化に応じ、同じく時系列で変化する。そして、入力インピーダンス回路部で分圧され、フィルタ回路で不要成分が除去され、アンプ回路で増幅され、Ａ／Ｄコンバータ回路にて、ある所定のクロック周波数でサンプリングされ、デジタル信号化される。
【００７７】
ここで、光音響信号を飽和させずに、表面血管に由来する信号を精度よく取得するための条件について説明する。
本例では、光音響信号を受信する受信回路において、Ａ／Ｄコンバータの入力許容電圧範囲が２Ｖｐ−ｐである場合を想定する。信号を飽和させないためには、光音響信号の振幅（ピーク−ピーク値）が２Ｖを超えない（アナログ信号を飽和させない）ような入力インピーダンス値（分圧によってアナログ信号の振幅を変更する）や、プリアンプの増幅率を決定する必要がある。加えて、受信回路においてフィルタ処理を行う場合、光音響信号の持つ周波数特性のうち、重要な特性を劣化させないようにフィルタの通過帯域特性を決定する必要がある。
【００７８】
シミュレーションによって、表面血管に由来する音響波が複数の音響波受信素子３００に到達するタイミングと、当該光音響信号の波形および強度が予測できる。また、表面血管から発生し、複数の音響波受信素子３００に到達する光音響波の周波数帯域特性も予測
できる。よって、所望の（表面血管に由来する）光音響波を飽和させずに取得するために、受信回路において行うべき処理の内容（最適な入力インピーダンス、フィルタ帯域、ゲインなどの設定値）を決定することが可能となる。
これらの条件を決定するために参照する情報は、光音響信号の波形情報（波形のピーク値、周波数特性）であってもよいし、シミュレーションの結果得られる他の情報であってもよい。結果として、受信した音響波に対応する光音響信号を飽和させず（すなわち、情報を劣化または喪失させず）にサンプリングすることが可能となる。
【００７９】
なお、シミュレーションから得た情報は画像再構成の際にも活用できる。これについては、ステップＳ９００において説明する。
【００８０】
次に、シミュレーション結果に基づいて、取得した光音響信号から、関心領域に由来する信号を抽出する処理について説明する。
図４Ｂに、光音響信号の受信状況の一例を示す。一点鎖線で示した領域Ａは、音響波受信素子３００−６の指向性の範囲を概念的に示したものである。図４Ｂの例の場合、音響波受信素子３００−６に対し、被検体表面（つまり皮膚）にある点Ｐ１に由来する音響波が、表面血管上にある点Ｐ２に由来する音響波よりも先に音響波受信素子３００−６に到達する。言い換えると、点Ｐ１に由来する音響波と、点Ｐ２に由来する音響波が異なるタイミングで音響波受信素子３００−６に到達する。
よって、シミュレーションから得られた、表面血管に由来する音響波の特性を考慮し、表面血管に由来する光音響信号を飽和させず、精度よく取得できるような適切な条件を設定したうえで、光音響波の受信および処理を行えばよい。
【００８１】
しかし、音響波受信素子３００−６の指向性の範囲が変わると、状況が変化する。
図４Ｃは、図４Ｂの例に比べて、音響波受信素子３００−６の指向性の範囲が広い場合を示した図である。一点鎖線で示した領域Ｂは、音響波受信素子３００−６の指向性の範囲を概念的に示したものである。図４Ｂと比較すると、音響波受信素子３００−６の指向性の範囲が広がっている。この例では、被検体表面（つまり皮膚）にある点Ｐ３，Ｐ４に由来する音響波が、表面血管上にある点Ｐ５に由来する音響波と同時に音響波受信素子３００−６に到達する。つまり、各点から生じた音響波が加算された状態で音響波受信素子３００−６に到達することになる。
【００８２】
このような場合、Ｐ３、Ｐ４、Ｐ５に由来する音響波（加算された音響波）の特性を損なわないような条件で光音響信号を一旦取得したうえで、点Ｐ５に由来する光音響波信号を抽出する処理を行えばよい。
信号抽出処理は、データ取得部７１８が行ってもよいし、後段（ステップＳ９００）で演算処理部７２８が行ってもよいが、ここで説明する。
【００８３】
<<光音響信号の抽出をデジタル信号から行う例>>
図５は、図４Ｃの例において、シミュレーションによって得られた光音響信号の例である。ここでは、点Ｐ３，Ｐ４，Ｐ５にそれぞれ由来する光音響信号（音響波受信素子３００が出力したアナログ信号波形）の波形情報Ｐ３（ｔ），Ｐ４（ｔ），Ｐ５（ｔ）が得られたものとする。そしてＰ３，Ｐ４，Ｐ５に由来する光音響波が加算された光音響波の波形情報Ｐａ（ｔ）が得られたものとする。（ｔは時間）
このとき、Ｐａ（ｔ）の最大値ｍａｘＡから、Ｐａ（ｔ）波形が飽和しないゲイン値、入力インピーダンス値、かつ、Ｐａ（ｔ）の持つ周波数帯域を損しないフィルタ帯域を設定し、Ｐａ（ｔ）の波形をサンプリングして、デジタル信号Ｐｄ（ｔ）を得る。
【００８４】
そして、点Ｐ３，Ｐ４，Ｐ５に由来する光音響信号の波形情報Ｐ３（ｔ），Ｐ４（ｔ），Ｐ５（ｔ）と、Ｐａ（ｔ）の情報に基づいて、関心領域（本例では点Ｐ５）に由来する
波形を、デジタル信号Ｐｄ（ｔ）から抽出する。
なお、デジタル信号Ｐｄ（ｔ）に含まれるＰ３（ｔ），Ｐ４（ｔ），Ｐ５（ｔ）に由来するデジタル信号成分Ｐ３ｄ（ｔ），Ｐ４ｄ（ｔ），Ｐ５ｄ（ｔ）、は、Ｐａ（ｔ）の信号振幅に占めるＰ３（ｔ），Ｐ４（ｔ），Ｐ５（ｔ）の信号振幅の割合により算出できる。
【００８５】
例えば、時刻ｔ１における、デジタル信号Ｐｄ（ｔ１）に含まれるＰ５由来のデジタル信号成分Ｐ５ｄ（ｔ１）は、以下のように算出できる。
Ｐ５ｄ（ｔ１）＝ Ｐｄ（ｔ１）×｛ｍａｘ５／（ｍａｘ３ + ｍａｘ４ + ｍａｘ５）
｝ ＝ Ｐｄ（ｔ１）×｛ｍａｘ５／（ｍａｘＡ）｝
時刻ｔ１におけるＰａ（ｔ）の振幅値Ｐａ（ｔ１）＝ｍａｘＡは、Ｐ３（ｔ），Ｐ４（ｔ），Ｐ５（ｔ）の時刻ｔ１における振幅値ｍａｘ３、ｍａｘ４、ｍａｘ５の和として表されるためである。
【００８６】
また、時刻ｔ２における、デジタル信号Ｐｄ（ｔ２）に含まれるＰ５由来のデジタル信号成分Ｐ５ｄ（ｔ２）は、以下のように算出できる。
Ｐ５ｄ（ｔ２）＝ Ｐｄ（ｔ２）×｛ｍｉｎ５／（ｍｉｎ３ + ｍｉｎ４ + ｍｉｎ５）
｝ ＝ Ｐｄ（ｔ２）×｛ｍｉｎ５／（ｍｉｎＡ）｝
時刻ｔ２におけるＰａ（ｔ）の振幅値Ｐａ（ｔ２）＝ｍｉｎＡは、Ｐ３（ｔ），Ｐ４（ｔ），Ｐ５（ｔ）の時刻ｔ２における振幅値ｍｉｎ３、ｍｉｎ４、ｍｉｎ５の和として表されるためである。
【００８７】
このように、複数の点に由来する信号が同時に音響波受信素子に到達する場合であっても、シミュレーションから得た事前情報をもとに、デジタル信号から定量性良く関心領域に由来する信号波形を抽出することが可能である。
【００８８】
<<光音響信号の抽出をアナログ信号から行う例>>
シミュレーションから得た事前情報をもとに、関心領域に由来する信号波形をアナログ信号から抽出することも可能である。
図６に、関心領域に由来する信号波形をアナログ信号から抽出するアナログ回路部７３５を示す。図６に示したアナログ回路部７３５は、入力インピーダンス回路７３３、フィルタ７３０、７３２、ローノイズアンプ（Low Noise Amplifier、以下ＬＮＡ）７３１、
ゲインアンプ（Variable Gain Amplifier、以下ＶＧＡ）７１９から成る。加えて、ＤＡ
コンバータ７３９、減算器７３８から構成される。
【００８９】
例えば、音響波受信素子３００が出力するアナログ電気信号波形Ｐａ（ｔ）がフィルタで処理され、α倍に増幅される場合に、Ｐ３（ｔ）とＰ４（ｔ）の和を同じくα倍増幅した値｛Ｐ３（ｔ） + Ｐ４（ｔ）｝×α をＰａ（ｔ）から減算する。
すると、以下のようになり、アナログ信号の減算処理によってＰ５（ｔ）を抽出することができる。
Ｐａ（ｔ）×α − ｛Ｐ３（ｔ）+Ｐ４（ｔ）｝×α 
＝ ｛Ｐ３（ｔ）+Ｐ４（ｔ）+Ｐ５（ｔ）｝×α− ｛Ｐ３（ｔ）+Ｐ４（ｔ）｝×α
＝ Ｐ５（ｔ）×α
【００９０】
例示したコンピュータ７００の構成に追加して、図６に示したように減算器７３８およびＤ／Ａコンバータ７３９を追加することで、関心領域に由来する信号波形をアナログ信号から抽出することができるようになる。
例えば、制御部７１１が出力する、アナログ信号Ｐ３＿４ａｇ（ｔ）（＝｛Ｐ３（ｔ）+Ｐ４（ｔ）｝×α）に相当するデジタル信号データＰ３＿４ｄｇ（ｔ）を、ＤＡコンバ
ータ７３９でアナログ信号化して、Ｐ３＿４ａｇ（ｔ）を再現する。また、減算器７３８
を用いてＰａ（ｔ）x αからＰ３＿４ａｇ（ｔ）を減算して、Ｐ５（ｔ）×αを得る。
【００９１】
このような処理を行うことで、シミュレーションから得た事前情報をもとに、アナログ信号から関心領域に由来する信号波形を抽出することも可能である。なお、アナログ信号減算処理は、必ずしも増幅されたアナログ信号に対して行う必要はない。増幅前の信号に対して信号減算処理を行ってもよいし、アナログ信号が複数回増幅される場合には、複数回の増幅段階のうち任意の段階で信号減算処理を行ってもよい。目的の信号を適切に抽出できる限り、アナログ信号の減算処理はどのように行ってもよい。
【００９２】
ステップＳ７００の説明を続ける。
<<光音響信号の取得処理の詳細>>
ここで、ステップＳ７００において、光音響信号を取得する処理の詳細について、図２を参照しながら説明する。なお、本例では、音響波受信素子が８つである場合の例を示すが、音響波受信素子の数はこれに限定されない。前述したように、データ取得部７１８は、アナログ回路部７３５およびデジタル回路部７３６から構成される。
【００９３】
まず、アナログ回路部７３５にて行われる処理について説明する。
アナログ回路部７３５は、は、入力インピーダンス回路７３３−１〜７３３−８、フィルタ７３０−１〜７３０−８、ローノイズアンプ７３１−１〜７３１−８、ゲインアンプ７１９−１〜７１９−８、フィルタ７３２−１〜７３２−８から構成される。
また、デジタル回路部７３６は、Ａ／Ｄ変換器（以下、ＡＤＣ）７１７−１〜７１７−８、ＦＩＦＯ７１６−１〜７１６−８、クロック供給部７１３、ＦＩＦＯ制御部７１２、制御部７１１、そしてセレクトスイッチ７１４から構成される。
なお、以降の説明において、特に区別する必要がない場合は、参照符号のうちハイフン以降を省略する。
【００９４】
次に、データ取得部７１８が光音響信号を取得する動作について説明する。
まず、複数の音響波受信素子３００が、光音響波を受信して電気信号に変換し、フィルタ７３０へ入力する。
入力インピーダンス回路７３３は、複数の音響波受信素子３００の出力インピーダンス値との関係によって、各出力信号を分圧してフィルタ７３０へ入力する。
【００９５】
また、フィルタ７３０は、電気信号をフィルタ処理して、ＬＮＡ７３１に出力する。ＬＮＡ７３１で増幅された電気信号は、ＶＧＡ７１９に出力される。また、ＶＧＡ７１９で増幅された電気信号は、フィルタ７３２に出力され、フィルタ処理された電気信号がＡＤＣ７１７へ出力され、デジタル信号へと変換される。
【００９６】
なお、入力インピーダンス回路７３３は、ある固定のインピーダンス値を持つものであってもよいが、制御部７１１によって、任意のインピーダンスに設定できるようにしてもよい。もしくは、複数の入力インピーダンス値を選択できるようにしておき、制御部７１１からの指示により、特定の入力インピーダンス値を選択および設定するようにしてもよい。
【００９７】
フィルタ７３０および７３２は、複数の音響波受信素子３００から入力された信号のＤＣ成分や、各音響波受信素子の帯域より低い周波数成分を除去するためのハイパスフィルタである。さらには、サンプリングの折り返し雑音（エイリアシング）や、各音響波受信素子の帯域よりも高い周波数成分を除去するためのローパスフィルタである。
例えば、フィルタ７３０をローパスフィルタ、フィルタ７３２をハイパスフィルタとしてもよいが、必ずしもフィルタの構成は例示したものに限定されない。例えば、ハイパスフィルタとローパスフィルタを個別に設けるのではなく、バンドパスフィルタを用いて、
低周波成分と高周波成分の双方を除去してもよい。
【００９８】
例えば、サンプリング周波数が２０ＭＨｚである場合、ローパスフィルタのカットオフ周波数は、サンプリング定理により１０ＭＨｚ以下である必要がある。ハイパスフィルタのカットオフ周波数は、複数の音響波受信素子３００の帯域以内の信号を除去しないよう留意して設定する必要がある。
なお、フィルタ７３０および７３２のカットオフ周波数は固定であってもよいが、制御部７１１によって、任意のカットオフ周波数に設定できるようにしてもよい。もしくは、複数のカットオフ周波数を選択できるようにしておき、制御部７１１からの指示により、カットオフ周波数を選択して設定するようにしてもよい。
【００９９】
ＬＮＡ７３１、ＶＧＡ７１９は、複数の音響波受信素子３００が出力する電気信号を増幅する。ＬＮＡ７３１、ＶＧＡ７１９における増幅率は、制御部７１１によって設定される。以上が、アナログ回路部７３５における信号処理の説明である。
【０１００】
次に、デジタル回路部７３６にて行われる処理について説明する。
ＡＤＣ７１７は、クロック供給部７１３が出力するクロックに従って、ある周波数で電気信号をサンプリングし、電気信号をデジタル信号に変換して、ＦＩＦＯメモリ（先入れ先出しメモリ。以下、ＦＩＦＯ）７１６へ出力する。ＦＩＦＯ７１６は、クロック供給部７１３が出力するクロックと、ＦＩＦＯ制御部７１２が出力する書き込みイネーブル信号に従って、ＡＤＣ７１７が出力したデジタル信号を記憶する。
【０１０１】
ＦＩＦＯ制御部７１２は、ステップＳ２００で取得した受信期間パラメータに基づいて、例えば図７に示されるように、ＦＩＦＯ７１６に対する書き込みイネーブル信号[7:0]
の出力タイミングを調整する。図７において、書き込みイネーブル信号のレベルがＬｏｗである時は、ＦＩＦＯ７１６に対する書き込みは行われず、書き込みイネーブル信号のレベルがＨｉｇｈである時は、ＦＩＦＯ７１６に対する書き込みが行われる。図７の例では、光照射が行われたタイミングＴ０の後に存在する、Ｔ１からＴ２までの期間が、光音響信号をＦＩＦＯ７１６に書き込んでいる期間である。
なお、クロック供給部７１３が出力するクロック（サンプリングクロック）の周波数は、利用者が予め設定することもできる。
【０１０２】
ＦＩＦＯ７１６は、クロック供給部７１３が出力するクロックと、ＦＩＦＯ制御部７１２が出力する読み出しイネーブル信号に従って、記憶したデジタル信号を上位演算器７１５へ転送する。
セレクトスイッチ７１４は、ＦＩＦＯ７１６−１〜８のうち、１つを選択して上位演算器７１５へ接続し、デジタル信号を上位演算器７１５内へ転送する。また、上位演算器７１５へ接続するＦＩＦＯの選択と、上位演算器７１５へのデジタル信号の転送を繰り返し、ＦＩＦＯ７１６−１〜８に記憶された全てのデジタル信号を転送する。なお、デジタル信号の転送先は、必ずしも上位演算器７１５である必要はなく、記憶部７２０であってもよい。上位演算器７１５が適切なタイミングでデジタル信号にアクセスできれば、デジタル信号の転送先はどこにあってもよい。
【０１０３】
また、図２では、セレクトスイッチ７１４が、ＦＩＦＯ７１６−１〜８のうち、１つを選択してデジタル信号を転送する構成を示したが、必ずしもこれに限定されない。例えば、ＦＩＦＯ７１６−１〜８のうち複数、もしくは全てに記憶されているデジタル信号を同時に上位演算器７１５へ転送してもよい。ＦＩＦＯ７１６−１〜８に記憶されたデジタル信号が適切に上位演算器７１５に転送されれば、特定の構成に限定されない。
【０１０４】
以上に説明したように、ステップＳ７００では、複数の音響波受信素子３００から出力
された電気信号が、アナログ回路部７３５およびデジタル回路部７３６によって処理され、デジタル信号として保存される。
なお、本実施形態では、書き込みイネーブル信号によってＦＩＦＯ７１６に光音響信号を書き込んだが、必ずしもこの手法に限定されない。例えば、ＦＩＦＯ７１６にクロックを常時供給する構成としたうえで、光音響信号の取得時にのみＡＤＣ７１７にクロックを供給するようにしてもよい。また、ＡＤＣ７１７にクロックを常時供給する構成としたうえで、光音響信号の取得時にのみＦＩＦＯ７１６にクロックを供給するようにしてもよい。また、光音響信号の取得時にのみ、ＦＩＦＯ７１６とＡＤＣ７１７の双方にクロックを供給するようにしてもよい。
また、本実施形態では、ＡＤＣ７１７が出力したデジタル信号を記憶する手段としてＦＩＦＯを例示したが、必ずしもこれに限定されない。例えば、ランダムアクセスメモリ等を用いてもよい。
【０１０５】
なお、不図示の接続スイッチを、複数の音響波受信素子３００とアナログ回路部７３５との間に設け、複数の音響波受信素子３００とアナログ回路部７３５との接続を切り替えながら信号を取得するようにしてもよい。この場合、アナログ回路部７３５のチャンネル数を複数の音響波受信素子３００より少なくすることができ、装置の小規模化が可能になる。
【０１０６】
また、不図示の接続スイッチをアナログ回路部７３５とデジタル回路部７３６の間に設け、アナログ回路部７３５の出力と、デジタル回路部７３６のＡＤＣ７１７との接続を切り替えるようにしてもよい。例えば、アナログ回路部７３５の複数の出力を、不図示の接続スイッチを介してデジタル回路部７３６の１つのＡＤＣ７１７に接続する。
そして、サンプリング波数の数倍以上の頻度で、アナログ回路部７３５の出力とＡＤＣ７１７との接続先を切り替えて、１回のサンプリング周期内で、タイミングをずらしながらサンプリングを実施する。このようにすることで、アナログ回路部７３５に用意されたチャンネル数より少ない数のＡＤＣ７１７で受信を行うことが可能になる。
以上が、ステップＳ７００における処理の説明である。
【０１０７】
予め設定された測定条件におけるデータ取得が完了すると、走査は完了となる（ステップＳ８００）。
【０１０８】
次に、ステップＳ９００にて、演算処理部７２８が、ステップＳ７００で取得したデジタル信号に対して画像再構成アルゴリズムを適用することで、被検体の光学特性に関する情報（吸収係数分布、酸素飽和度等。被検体情報と称する）を取得する。
【０１０９】
被検体情報を取得するための画像再構成アルゴリズムとして、トモグラフィ技術で通常に用いられるタイムドメインあるいはフーリエドメインでの逆投影などが知られている。なお、画像再構成の時間を多く取れる場合は、繰り返し処理による逆問題解析法などの画像再構成手法を用いることもできる。所望の画像再構成が実施されうる限り、画像再構成アルゴリズムは特定の範囲に限定されない。
【０１１０】
本実施形態に係る被検体情報取得装置においては、生体組織内の血液中に含まれるＨｂとＨｂＯ２の濃度比を、複数の波長を用いて測定し、酸素飽和度を算出する。
酸素飽和度を取得する場合、被検体の表面血管に該当する部分に対応する信号を精度よく抽出することが望ましい。前述したように、シミュレーションによって、表面血管に由来する信号の周波数帯域を算出しておけば、当該周波数帯域にマッチした通過帯域特性を持つフィルタを、表面血管存在部分に該当する信号に適用することができる。この他にも、表面血管が存在する部分に画像処理フィルタをかけてもよい。結果的に、表面血管に由来する信号を精度よく抽出し、視認性の優れた画像を生成することが可能となる。
【０１１１】
なお、被検体情報を取得する際には、ステップＳ３００で行ったシミュレーションの結果に基づいて、受信した光音響信号の補正を行うことが好ましい。
例えば、複数の音響波受信素子３００のそれぞれに対して異なるサンプリング周波数を設定した場合、サンプリング周波数の違いを考慮して画像再構成に用いるデジタル信号を適切に選択することが必要である。
また、複数の音響波受信素子３００のそれぞれに対して異なる入力インピーダンスを適用した場合、入力インピーダンスの差異に起因して、ＡＤＣ７１７における信号レベルがそれぞれ変わってしまう。そこで、ＡＤＣ７１７におけるＡ／Ｄ変換後に、複数の音響波受信素子３００のそれぞれに対応するデジタル信号を適宜増幅するようにしてもよい。
【０１１２】
さらに、複数の音響波受信素子３００のそれぞれに対してＬＮＡ７３１、ＶＧＡ７１９で異なるゲインを適用した場合、ＡＤＣ７１７においてＡ／Ｄ変換後に、当該ゲインの違いを補正する処理を加えてもよい。こうすることで、複数の音響波受信素子３００の全てについて実質的にゲインを揃えることができる。複数の音響波受信素子３００ごとにゲイン値がばらついたままでは、取得した形態情報や機能情報の定量性に劣化が生じる。複数の音響波受信素子３００のそれぞれに対応するゲイン値を正しく把握し、かつゲイン値のばらつきを補正する処理を行うことで、より正しい濃度比測定が行える。
【０１１３】
また、複数の音響波受信素子３００とアナログ回路部７３５（入力インピーダンス回路７３３、フィルタ７３０、フィルタ７３２、ＬＮＡ７３１、ＶＧＡ７１９、ＡＤＣ７１７）を連結した"系"の伝達関数の逆関数を利用した処理も可能である。例えば、伝達関数の逆関数を用いて、取得したデジタル信号から複数の音響波受信素子３００が受信したアナログの光音響信号波形を復元するデコンボリューションと呼ばれる処理がある。
このような場合、アナログ回路部７３５の受信条件設定に応じて伝達関数が変わると、対応する逆関数にも変化が生じるため、信号波形の復元処理結果にも変化が生じる。伝達関数が変わった場合、対応する逆関数の変化を考慮せずに逆関数を固定で信号波形の復元処理を実施すると、正しい復元処理ができない。そこで、アナログ回路部７３５の設定条件が更新されるごとに、アナログ回路部７３５の伝達関数および逆関数を導出し、信号波形の復元処理を正確に実施することが好ましい。
【０１１４】
次に、ステップＳ９００にて、表示装置９００が、取得された被検体情報を表示し、処理は終了する（ステップＳ１１００）。
【０１１５】
以上説明したように、第一の実施形態に係る被検体情報取得装置は、皮膚に由来する光音響波と、表面血管に由来する光音響波とを区別することで、第一に、再構成が行われる前の光音響信号に対して適切な処理を行うことができるようになる。また、第二に、表面血管に由来する信号を抽出できるようになる。これにより、被検体情報の算出を精度よく行うことができる。
【０１１６】
なお、本実施形態では、光音響画像のみを取得したが、光音響画像の取得と超音波画像の取得を同時に行い、被検体情報を表示する際に、双方を重畳表示ないし並列表示するようにしてもよい。また、光音響画像と超音波画像を個別に表示するようにしてもよい。
【０１１７】
また、光音響信号を処理するための条件は、複数の音響波受信素子３００全てに対して同一である必要はなく、複数の音響波受信素子の１つ１つに対し、個別に設定されてもよい。もしくは、同一の条件をある数の素子にまとめて割り当ててもよい。つまり、複数の音響波受信素子３００を複数のグループに分け、各グループに個別に条件を割り当ててもよい。
【０１１８】
また、複数の音響波受信素子３００のそれぞれに対して、異なる帯域特性のフィルタ（フィルタ７３０，７３２）を適用した場合、同じ帯域特性のフィルタを用いた音響波受信素子同士をグループ分けして、個別に画像再構成するようにしてもよい。
【０１１９】
また、本実施形態では、シミュレーションから得た事前情報に基づいて、関心領域に由来する信号波形を抽出するものとしたが、デジタル信号波形そのものから、関心領域に由来する信号波形を分離抽出する手法を用いてもよい。例えば、特異スペクトル解析、独立成分分析、主成分分析といった重畳信号の分離手法や、統計的信号処理手法に基づく関心領域に由来する信号波形推定法、もしくは多変量解析手法に基づく信号抽出法、またはウェーブレット変換などを適用してもよい。また、それらの手法を組み合わせて適用してもよい。
また、関心領域の内外にそれぞれ由来する光音響波を分離抽出する手法は、被検体内部の表面血管の点に由来する光音響波を適切に分離抽出できればよく、特定の手法に限定されない。
【０１２０】
また、本実施形態では、表面（つまり皮膚）に由来する光音響波と、被検体内部の表面血管の点に由来する光音響波を区別する例を挙げたが、これに限られない。例えば、被検体内部の表面血管上の異なる点にそれぞれ由来する光音響波が、同時に音響波受信素子３００に到達するような場合であっても、同様の手法によって光音響波の分離抽出を行うことができる。このようにすることで、被検体内部の表面血管における、局在的な酸素飽和度等の血中濃度成分測定の定量性が向上する。
いずれにせよ、指向性が異常に広い場合、それだけでなく指向性にサイドローブが生じる場合などでは、光音響波を受信する条件設定や、受信後の処理方法を考慮することで、血中成分濃度測定の定量性が向上する。
【０１２１】
（第二の実施形態）
光音響装置においては、装置のシステムノイズが適切に保たれているかを監視することが好ましい。システムノイズが何らかの不具合によって大きくなってしまうと、取得した光音響信号のＳ／Ｎ比が悪化し、機能情報の定量性が劣化するなどの不具合が生じてしまうためである。第二の実施形態は、これに対応するため、装置のシステムノイズを計測する実施形態である。
なお、システムノイズとは、被検体の測定を行っていないタイミングにおいて、データ取得部７１８によって取得されるデータを指す。第二の実施形態に係る被検体情報取得装置は、コンピュータ７００が、システムノイズを計測する機能をさらに有している。
【０１２２】
図８Ａ〜図８Ｂに示したタイミングチャートを参照しながら、本実施形態に係る被検体情報取得装置における、システムノイズのチェック方法を説明する。なお、図８Ａ〜図８Ｂにおいて、Ｈは、被検体情報取得装置を構成する各モジュールの電源が投入されている状態を表し、Ｌは、電源が切断されている状態を表す。
【０１２３】
システムノイズの計測は、図８Ａに示したように、装置の構成モジュールの電源が全て投入された後で、光音響波を受信していないタイミングＴｎで行うことができる。
もしくは、各構成モジュールの電源を順次投入しながら、随時行うようにしてもよい。例えば、Ｔｎ０，Ｔｎ１，Ｔｎ２，Ｔｎ３といったタイミングで行ってもよい。構成モジュールのいずれかが不具合により大きなノイズを発生し、システムノイズを悪化させていた場合、このように測定することで、どの構成モジュールに不具合があるかを特定しやすくなる。
【０１２４】
または、図８Ｂに示したように、特定の構成モジュールの電源を切断した状態で、他の構成モジュールの電源を順次投入して測定を行うようにしてもよい。図８Ｂは、給排液ユ
ニットの電源を切断した状態（Phase A）と、給排液ユニットの電源を投入した状態（Phase B）に分割して測定を行う場合の例である。このようにすることで、支持体４００に音響マッチング材８００を入れたことによるシステムノイズの変化を確認することができる。
なお、図８Ａおよび図８Ｂの例では、構成モジュールの電源を順次投入したが、各構成モジュールのうち、それぞれ一つだけ電源を投入した状態で測定を行うようにしてもよい。このようにすることで、各構成モジュール自体がシステムノイズに及ぼす影響を確認することができる。
【０１２５】
また、複数の音響波受信素子３００として、ＣＭＵＴなどの静電容量型の素子を用いる場合、素子に対して複数の電圧が印加される。このように、音響波受信素子のそれぞれに対して電源供給が必要なシステムである場合、素子ごとに電源供給を制御しながらシステムノイズ計測を実施するようにしてもよい。このようにすることで、複数の音響波受信素子のそれぞれがシステムノイズに及ぼす影響を確認することができる。
なお、システムノイズの計測方法は、前述したものに限られない。また、構成モジュールの電源を投入する順番も、特定の順番に限られない。
【０１２６】
（第三の実施形態）
第三の実施形態に係る被検体情報取得装置は、音響波受信素子が半球状の支持体に取り付けられているのではなく、ハンドヘルド型のプローブに取り付けられているという点において、第一の実施形態に係る被検体情報取得装置と相違する。
【０１２７】
第三の実施形態に係る被検体情報取得装置を、図９を参照しながら説明する。
第三の実施形態に係る被検体情報取得装置は、光源１５４０、光学系１５６０、光音響プローブ１５００、バンドルファイバ１５５０、光検出器１５９０を有して構成される。さらに、第三の実施形態に係る被検体情報取得装置は、形態取得部６００、コンピュータ１７００、表示装置９００、入力部１０００を有している。
光音響プローブ１５００は、出射端１５１０と、光音響波を受信する受信部１５２０と、走査位置センサ１５３０からなる。受信部１５２０には、複数の音響波受信素子１６５０が配置される。
また、コンピュータ１７００は、演算部１７１０および記憶部１７２０を有している。
【０１２８】
図１０は、第三の実施形態におけるコンピュータ１７００の構成をより詳細に示した図である。図示したように、演算部１７１０は、データ取得部１７１８と、演算処理部７２８から構成される。
【０１２９】
光源１５４０は、コンピュータ１７００の制御に従って照射光を発し、光学系１５６０を介してバンドルファイバ１５５０へ入射させる。入射した光は光音響プローブ１５００まで伝送され、出射端１５１０から出射される。
【０１３０】
出射端１５１０から被検体へ照射された光は被検体の内部で拡散し、光吸収体から光音響波が発せられる。
受信部１５２０には、ピエゾ素子やＣＭＵＴなどを用いた複数の音響波受信素子１６５０がアレイ状に配置されている。被検体から発せられた光音響波は、複数の音響波受信素子１６５０によって受信され、電気信号に変換されたのち、コンピュータ１７００へ送られる。
【０１３１】
なお、図９の例では、バンドルファイバ１５５０を途中で分岐させ、受信部１５２０を挟んで取り付けられた二つの出射端１５１０にそれぞれ接続したが、分岐の個数はこれに限定されない。例えば、出射端１５１０を一つとしてもよい。また、本例では、光の伝送
にバンドルファイバ１５５０を用いたが、ミラーやプリズムなどの光学素子を介して光を伝送してもよい。
光検出器１５９０は、第一の実施形態における光検出器１２０と同様であるため、説明は省略する。
【０１３２】
第三の実施形態における形態取得部６００は、光音響プローブ１５００に配置された、複数の音響波受信素子１６５０から超音波を送受信可能な回路構成となっており、被検体表面（もしくは内部）の形態情報を取得することができる。具体的には、光音響プローブ１５００に配置される複数の音響波受信素子１６５０の一部、もしくは全てを用いて音響波を送信し、反射波を解析することで、被検体の形態情報を取得する。
なお、反射波の解析処理は、形態取得部６００が行ってもよいし、コンピュータ１７００が有する演算処理部１７２８が行ってもよい。
取得された形態情報は、第一の実施形態と同様に一時的に記憶される。
【０１３３】
また、形態取得部６００は、光音響プローブ１５００に配置される複数の音響波受信素子１６５０から音響波を送受信し、第一の実施形態と同様にドップラーエコー法を用いて、被検体内部において血流が存在する部分を確認する機能を有する。血流が存在する部分についての情報は、２次元的、３次元的に取得することが可能であり、表示させることができる。さらには、血流の変化の状態を、従来のドップラーエコー法と同様、リアルタイムに取得することができる。ドップラーエコー法の演算は、形態取得部６００が行ってもよいし、コンピュータ１７００（演算処理部１７２８）が行ってもよい。
このように、第三の実施形態における光音響プローブ１５００は、光音響信号と、超音波エコー信号（超音波信号）の両方を取得可能な構成となっている。
【０１３４】
光音響プローブ１５００に配置される音響波受信素子のアレイ形状は特定のものに限定されない。例えば、１Ｄのリニアアレイ、セクタアレイ、コンベックスアレイ、２Ｄアレイ、スパースアレイなどを用いてもよい。
【０１３５】
第三の実施形態では、複数の音響波受信素子１６５０を用いて、光音響波と反射波（超音波エコー）の双方を受信する。これを実現する方法として、以下のようなものがある。
第一の方法は、複数の音響波受信素子１６５０の全てを用いて、光音響波と、超音波エコーを別のタイミングで受信する方法である。
また、第二の方法は、複数の音響波受信素子１６５０の少なくとも一部を光音響波の受信用に割り当て、残りを超音波エコーの受信用に割り当てる方法である。
【０１３６】
例えば、図１１は、複数の音響波受信素子１６５０が１０素子×１０素子の２Ｄアレイ形状である場合の例である。
例えば、１００個の素子全てを用いて、光音響波の受信および超音波エコーの受信を異なるタイミングでそれぞれ行うようにしてもよい。また、領域Ａに存在する素子が光音響波を受信し、領域Ｂに存在する素子が超音波エコーを受信するようにしてもよい。
このように、素子の分割は時間的に行ってもよいし、空間的に行ってもよい。
【０１３７】
また、図９では、光音響プローブ１５００を直接被検体に接触させて測定する様態を示したが、測定方法はこれに限定されない。例えば、被検体を保持部材で挟み込むなどして固定し、保持部材越しに測定を行うようにしてもよい。もしくは、保持部材を光音響プローブ１５００と被検体との間に設置し、保持部材越しに測定を行うようにしてもよい。保持部材は、樹脂やゲル状部材で構成することも可能で、光と超音波の両方を良好に透過させる材質のものが望ましい。このようにすることで、被検体を安定した状態で測定することが可能となる。
【０１３８】
次に、図１２に示すフローチャートを参照しながら、第三の実施形態に係る被検体情報取得装置の動作について説明する。
【０１３９】
まず、ステップＳ２１００にて、被検体の測定を開始する。なおこの際、光音響プローブ１５００と被検体との間に、音響マッチングを改善するためのジェルなどを塗布することで、取得できる光音響画像（超音波エコー画像）の画質を向上させることができる。
【０１４０】
ステップＳ２２００の処理は、第一の実施形態におけるステップＳ２００と同様である。
なお、被検体内部における血管形状に関する情報を取得できれば、ドップラー撮像以外の手法を利用してもよい。例えば、他のモダリティ（Ｘ線ＣＴ、ＭＲＩなど）で事前に取得した被検体内外の形態情報を取得するようにしてもよい。
また、超音波エコー信号を取得する回数は、特定の回数に制限されない。光音響測定に影響のない範囲で、いかなる設定も可能であるものとする。
【０１４１】
ステップＳ２３００，Ｓ２４００の処理は、第一の実施形態におけるステップＳ３００，Ｓ４００とそれぞれ同様である。
【０１４２】
第一の実施形態では、被検体の表面（皮膚）に由来する光音響波と、被検体内部の表面血管に由来する光音響波が同時に音響波受信素子に到達する場合を想定して処理を行った。
しかし、第三の実施形態では、被検体の表面に音響波受信素子を直接接触させるため、このようなケースが生じにくい。
ただし、第三の実施形態においても、音響波受信素子１６５０と被検体との間に一定以上の厚みをもつ部材をはさんで測定を行う場合は状況が異なる。このよう場合、被検体の表面（皮膚）に由来する光音響波と、被検体内部の表面血管に由来する光音響波が同時に音響波受信素子１６５０に到達する状況が生じ得るため、光音響波を受信する条件設定や、受信後の処理方法を考慮することが好ましい。
【０１４３】
ステップＳ２５００，Ｓ２６００の処理は、第一の実施形態におけるステップＳ５００，Ｓ６００とそれぞれ同様である。
ステップＳ２７００の処理は、第一の実施形態におけるステップＳ７００と同様であるが、本実施形態は、一回の光照射を行うごとに、演算部１７１０が被検体情報を取得して表示するという点において、第一の実施形態と相違する。
【０１４４】
予め設定された測定条件におけるデータ取得が完了すると、走査は完了となり（ステップＳ２８００）、処理は終了する（ステップＳ２９００）
【０１４５】
第三の実施形態では、超音波エコーによる測定と、光音響測定とを時分割にて行う。これについて、図１３を参照しながら説明する。図１３は、超音波エコー信号の取得、光照射、光音響信号の取得の順番を表した図である。
【０１４６】
測定が開始されると、被検体内外の形態情報を取得するための超音波エコー信号の取得が行われる。超音波測定はＮ（Ｎ≧０）回実施され、１回の超音波測定で、超音波の送信と受信が実行される。この超音波測定によって得られた超音波エコー信号から、利用者が設定した幅、深さ、フレーム数に応じた超音波エコー画像が取得される。そして、超音波エコー画像を解析した結果に基づいて、光音響測定（光照射および光音響信号の取得）が行われる。この一連の動作を１つの測定サイクルと定義する。
【０１４７】
なお、本例では１つの測定サイクル中に、超音波の送信と受信がＮ（Ｎ≧０）回ずつ実
施されるが、Ｎの値は、光音響信号の取得に影響のない範囲であれば特定の回数に制限されない。また、各測定サイクルで、Ｎは同一である必要はない。
【０１４８】
第三の実施形態においては、前述したような測定サイクルを繰り返すことにより、精度よく光音響測定を行うために必要な情報をリアルタイムで取得することができる。
【０１４９】
なお、光音響トモグラフィでは、被検体の過度の温度上昇を避けるため、光の照射間隔をある一定時間（数十ミリ秒）以上に設定する必要がある。よって、光の照射間隔以内に超音波測定を完了することができれば、光音響画像のフレームレートには影響しない。
【０１５０】
なお、図１３では、超音波測定を行ったのちに光音響測定を行っているが、必ずしもこの順でなくてもよい。例えば、各測定サイクルにおいて、超音波測定を０回に設定することも可能である。例えば、ある測定サイクルにおいて、超音波測定を省略して光音響測定を行ってもよい。
また、本例では、サイクル毎に超音波測定を行っているが、光音響信号を処理するための条件は毎回更新しなくてもよい。例えば、利用者が指定したタイミングで条件を更新するようにしてもよい。このため、例えば、入力部１０００もしくは光音響プローブ１５００にボタンなどを設けておき、入力に応じて条件を更新するようにしてもよい。
【０１５１】
また、本実施形態では、１回の光照射ごとに被検体情報を取得し、表示する例を示したが、必ずしもこれに限定されない。例えば、複数回の光照射で取得したデジタル信号を蓄積し、被検体情報をまとめて取得するようにしてもよい。例えば、何回光照射を行うごとに被検体情報を取得するかを設定できるようにしてもよい。
また、取得した光音響信号を一時的に記憶し、被検体情報の取得および表示を、走査が完了した後で実施するようにしてもよい。このようにすることで、演算に時間がかかり、リアルタイムでは表示できないような緻密な再構成画像を得られるようになる。
【０１５２】
ここで、図１０を参照しながら、第三の実施形態におけるコンピュータ１７００の構成について説明する。なお、図１０の例では受信回路が８系統あるが、受信回路数はこれに限定されない。
デジタル回路部７３６、演算処理部７２８については、第一の実施形態と同様であるため説明は省略し、ここでは、第三の実施形態におけるデータ取得部１７１８と、第一の実施形態におけるデータ取得部７１８との相違点について説明する。
【０１５３】
データ取得部１７１８は、接続スイッチ１７１１を有しているという点において、第一の実施形態におけるデータ取得部７１８と相違する。それ以外の構成要素については、第一の実施形態で説明したものと同様である。
接続スイッチ１７１１は、光音響プローブ１５００とアナログ回路部１７３５との接続を切り替える手段である。具体的には、光音響プローブ１５００に配置された複数の音響波受信素子１６５０のうち、所望の音響波受信素子を選択してアナログ回路部１７３５と接続することができる。
【０１５４】
なお、図１３では、データ取得部１７１８が、同一の回路で超音波信号と光音響信号を取得する構成となっているが、必ずしもこれに限定されない。システム仕様の要求に応じて、超音波信号と光音響信号の取得を別の回路で行うようにしてもよい。
【０１５５】
第三の実施形態では、形態取得部６００が、配線６１０、バス２１００、接続スイッチ１７１１を介して、超音波を送信するための送信波形電圧を光音響プローブ１５００に伝送する。このようにすることで、複数の音響波受信素子１６５０のうち、特定の音響波受信素子を選択して超音波を送信することができる。
【０１５６】
なお、本実施形態の説明では、ハンドヘルド型のプローブを手動で走査する例を挙げたが、たとえばガイドに沿ってプローブを機械的に走査させるようにしてもよい。
【０１５７】
（変形例）
なお、各実施形態の説明は本発明を説明する上での例示であり、本発明は、発明の趣旨を逸脱しない範囲で適宜変更または組み合わせて実施することができる。
例えば、本発明は、上記処理の少なくとも一部を含む被検体情報取得装置や情報処理装置として実施することもできる。また、本発明は、上記処理の少なくとも一部を含む被検体情報取得方法や情報処理方法として実施することもできる。上記処理や手段は、技術的な矛盾が生じない限りにおいて、自由に組み合わせて実施することができる。
【０１５８】
また、第二の実施形態で述べたように、システムノイズの計測は手動で行ってもよいが、所定のタイミング（例えば、装置の起動時や設定されたタイミング）で、自動で行ってもよい。また、外部からの指示に基づいてシステムノイズの計測を実行するようにしてもよい。例えば、バス２０００やバス２１００に、リモートアクセスのためのモジュールを接続し、当該モジュール経由で遠隔で実行してもよい。このようにすることで、メンテナンス性が向上する。また、計測されたシステムノイズ情報は、被検体情報取得装置に累積保存され、随時参照可能にしてもよい。
【０１５９】
また、システムノイズに関する情報以外を、リモート伝送してもよい。例えば、被検体情報取得装置が取得した情報（光音響信号、超音波信号、演算されたパラメータ値、光音響画像、超音波画像、装置の動作状態、メンテナンス情報など）を記憶し、無線やネットワーク経由で伝送できるようにしてもよい。また、無線やインターネット経由で装置のオペレーションを可能にしてもよい。
【０１６０】
また、実施形態の説明では、測定対象の領域（関心領域）として表面血管を挙げたが、測定の対象は必ずしも表面血管でなくてもよい。
また、実施形態の説明では、機能情報として酸素飽和度を計測する例を挙げたが、機能情報は必ずしも酸素飽和度でなくてもよい。例えば、糖化タンパク質、グルコース、脂肪、コレステロール、コラーゲンなど、電磁波の波長に対して異なる吸収係数を持ち、光音響技術にて判別可能ないかなる物質に対しても本発明は適用可能である。
【０１６１】
本発明は、以下の処理を実行することによっても実現される。即ち、上述した各実施形態の１以上の機能を実現するプログラムを、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータにおける１つ以上のプロセッサがプログラムを読み出して実行する処理でも実現可能である。また、１以上の機能を実現する回路（例えば、ＦＰＧＡやＡＳＩＣ）によっても実現可能である。
【符号の説明】
【０１６２】
１００：光源、３００：音響波受信素子、４００：支持体、６００：形態取得部、７００：コンピュータ
【図１】
2018110734.tif 000009
【図２】
2018110734.tif 000010
【図３】
2018110734.tif 000011
【図４Ａ】
2018110734.tif 000012
【図４Ｂ】
2018110734.tif 000013
【図４Ｃ】
2018110734.tif 000014
【図５】
2018110734.tif 000015
【図６】
2018110734.tif 000016
【図７】
2018110734.tif 000017
【図８Ａ】
2018110734.tif 000018
【図８Ｂ】
2018110734.tif 000019
【図９】
2018110734.tif 000020
【図１０】
2018110734.tif 000021
【図１１】
2018110734.tif 000022
【図１２】
2018110734.tif 000023
【図１３】
2018110734.tif 000024