(57)【要約】【課題】光音響測定において、簡易な手法で安定した濃度情報を取得する。【解決手段】光出射部から出射された第１の波長を持つ第１の光と第１の波長と異なる第２の波長を持つ第２の光が照射された被検体内部の物質から発生する光音響波を探触子が受信して出力される、第１の光に対応する第１の信号と、第２の光に対応する第２の信号と、を取得する信号取得部と、光出射部から出射された光が照射される照射面における、第１の光と第２の光の光量比を示す光量比分布を取得する光量比分布取得部と、光量比分布と、第１の信号および第２の信号を用いて、被検体内部における物質の三次元の濃度情報を取得する濃度情報取得部を備え、濃度情報取得部は、被検体内部の第１の位置における光量比として、照射面内において第１の位置に対応する位置である第２の位置における光量比に基づく値を割り当てる被検体情報取得装置。【選択図】図１
2019193732.tif 000002
【特許請求の範囲】
【請求項１】
光出射部から出射された、第１の波長を持つ第１の光と、前記第１の波長と異なる第２の波長を持つ第２の光と、が照射された被検体内部の物質から発生する光音響波を探触子が受信して出力される、前記第１の光に対応する第１の信号と、前記第２の光に対応する第２の信号と、を取得する信号取得部と、
前記光出射部から出射された光が照射される照射面における、前記第１の光と前記第２の光の光量比を示す光量比分布を取得する光量比分布取得部と、
前記光量比分布と、前記第１の信号および前記第２の信号を用いて、前記被検体内部における前記物質の三次元の濃度情報を取得する濃度情報取得部と、
を備え、
前記濃度情報取得部は、前記被検体内部の第１の位置における前記光量比として、前記照射面内において前記第１の位置に対応する位置である第２の位置における前記光量比に基づく値を割り当てる
ことを特徴とする被検体情報取得装置。
【請求項２】
前記被検体内部の前記第１の位置は、前記照射面の前記第２の位置を入射位置とする光の進行方向の位置である
ことを特徴とする請求項１に記載の被検体情報取得装置。
【請求項３】
前記濃度情報取得部は、前記第２の位置を入射位置とする光の進行方向において、前記入射位置からの距離によらず略一定の前記光量比を用いる
ことを特徴とする請求項２に記載の被検体情報取得装置。
【請求項４】
前記被検体内部の前記第１の位置は、前記照射面の前記第２の位置が光の入射位置であるときに、前記照射面に略垂直な方向の位置である
ことを特徴とする請求項１に記載の被検体情報取得装置。
【請求項５】
前記濃度情報取得部は、前記第２の位置を起点として、前記照射面に略垂直な方向にある位置において、前記入射位置からの距離によらず略一定の前記光量比を用いる
ことを特徴とする請求項４に記載の被検体情報取得装置。
【請求項６】
前記濃度情報取得部は、前記入射位置からの光の伝搬距離が所定の閾値より小さい前記被検体内部の領域においては、前記光量比分布を用いて前記濃度情報を取得し、前記伝搬距離が前記所定の閾値以上の前記被検体内部の領域においては、前記光量比分布を用いずに前記濃度情報を取得する
ことを特徴とする請求項２から５のいずれか１項に記載の被検体情報取得装置。
【請求項７】
前記濃度情報取得部は、
ユーザーにより指定された前記濃度情報に関して許容できる誤差の範囲を取得し、
前記所定の閾値を、前記光量比分布を用いる手法で取得した前記濃度情報の誤差が前記許容できる誤差の範囲を超えるかどうかに基づいて決定する
ことを特徴とする請求項６に記載の被検体情報取得装置。
【請求項８】
前記濃度情報取得部は、ユーザーにより指定された前記被検体内部における関心領域を示す情報を取得し、前記関心領域の範囲と、前記伝搬距離に関する前記所定の閾値とに基づいて、前記濃度情報を取得する
ことを特徴とする請求項６または７に記載の被検体情報取得装置。
【請求項９】
前記濃度情報取得部は、前記濃度情報として前記被検体内部のヘモグロビンの酸素飽和
度を取得するものであり、前記入射位置からの光の伝搬距離が１８ｍｍ以下の場合に、前記光量比分布を用いて前記濃度情報を取得する
ことを特徴とする請求項６から８のいずれか１項に記載の被検体情報取得装置。
【請求項１０】
前記濃度情報取得部は、前記濃度情報として前記被検体内部のヘモグロビンの酸素飽和度を取得する
ことを特徴とする請求項１から９のいずれか１項に記載の被検体情報取得装置。
【請求項１１】
前記濃度情報取得部は、
前記第１の信号から前記第１の光に対応する第１の初期音圧分布を取得し、前記第２の信号から前記第２の光に対応する第２の初期音圧分布を取得し、
前記第１の初期音圧分布および前記第２の初期音圧分布ならびに前記光量比を用いて前記被検体内部における前記第１の光と前記第２の光のそれぞれに対応する吸収係数の比である吸収係数比を取得し、
前記吸収係数比と、前記ヘモグロビンの光学特性に基づいて、前記酸素飽和度を取得する
ことを特徴とする請求項１０に記載の被検体情報取得装置。
【請求項１２】
前記光量比分布取得部は、前記光出射部から出射された光の波長ごとの光量プロファイルに基づいて、前記光量比分布を取得する
ことを特徴とする請求項１から１１のいずれか１項に記載の被検体情報取得装置。
【請求項１３】
前記被検体を平面形状に保持する保持部をさらに備え、
前記照射面は、前記保持部が前記被検体を保持する面である
ことを特徴とする請求項１から１２のいずれか１項に記載の被検体情報取得装置。
【請求項１４】
前記被検体は手または足であり、
前記保持部は、前記手または足の指の間を突き当てるための突き当て部を有する
ことを特徴とする請求項１３に記載の被検体情報取得装置。
【請求項１５】
前記第１の信号および前記第２の信号は、前記光出射部および前記探触子が移動機構の走査によって前記被検体に対して移動しながらなされた光音響測定によって取得されたものであり、
前記光量比分布取得部は、前記光出射部および前記探触子が移動した各位置における照射面を合わせた仮想的な照射面における前記光量比分布を取得し、
前記濃度情報取得部は、前記仮想的な照射面における前記光量比分布を用いて、前記濃度情報を取得する
ことを特徴とする請求項１から１４のいずれか１項に記載の被検体情報取得装置。
【請求項１６】
光出射部から出射された、第１の波長を持つ第１の光と、前記第１の波長と異なる第２の波長を持つ第２の光と、が照射された被検体内部の物質から発生する光音響波を探触子が受信して出力される、前記第１の光に対応する第１の信号と、前記第２の光に対応する第２の信号と、を取得する信号取得部と、
前記光出射部から出射された光が照射される照射面における、前記第１の光と前記第２の光の光量比を示す光量比分布を取得する光量比分布取得部と、
前記光量比分布と、前記第１の信号および前記第２の信号を用いて、前記被検体内部における前記物質の三次元の濃度情報を取得する濃度情報取得部と、
を備え、
前記濃度情報取得部は、
前記第１の信号から前記第１の光に対応する第１の初期音圧分布を取得し、前記第２
の信号から前記第２の光に対応する第２の初期音圧分布を取得し、
前記被検体内部の第１の位置における光量として、前記照射面内において前記第１の位置に対応する位置である第２の位置における光量に基づく値を割り当て、
割り当てられた前記光量によって前記第１の初期音圧分布および前記第２の初期音圧分布を補正し、補正された前記第１の初期音圧分布および前記第２の初期音圧分布を用いて前記濃度情報を取得する
ことを特徴とする被検体情報取得装置。
【請求項１７】
光出射部から出射された、第１の波長を持つ第１の光と、前記第１の波長と異なる第２の波長を持つ第２の光と、が照射された被検体内部の物質から発生する光音響波を探触子が受信して出力される、前記第１の光に対応する第１の信号と、前記第２の光に対応する第２の信号と、を取得する信号取得部と、
前記光出射部から出射された光が照射される照射面における、前記第１の光と前記第２の光の光量比を示す光量比分布を取得する光量比分布取得部と、
前記光量比分布と、前記第１の信号および前記第２の信号を用いて、前記被検体内部における前記物質の三次元の濃度情報を取得する濃度情報取得部と、
を備えることを特徴とする被検体情報取得装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被検体情報取得装置に関する。
【背景技術】
【０００２】
光を用いたイメージング技術の一つとして、光音響イメージング技術がある。光音響イメージングでは、まず、光源から発生したパルス光が被検体に照射される。被検体内の光吸収体が、被検体内で伝搬・拡散した照射光のエネルギーを吸収すると、音響波（以降、光音響波と呼ぶ）が発生する。トランスデューサは、この光音響波を受信して受信信号に変換する。そして、信号処理装置が受信信号を解析処理することで、被検体内部の光学特性値に関する情報を取得する。これにより、被検体内の光吸収に関連する特性値分布が可視化される。
【０００３】
また、互いに異なる波長を有する光を被検体に照射することにより、被検体内に存在する物質（光吸収体）の濃度に関する分布を取得できる。特に、互いに異なる波長を有する光を被検体に照射し、各波長において血管内の血液の光吸収に関する情報分布を取得することにより、オキシヘモグロビンＨｂＯとデオキシヘモグロビンＨｂとの濃度を取得し、さらに血液の酸素飽和度を取得できる。例えば、２波長の光を用いた場合、酸素飽和度分布ＳＯ２（ｒ）は次式（１）で求められる。
【数１】
2019193732.tif 000003
【０００４】
ここで、μａ(λ１)（ｒ）は、ある位置ｒでの波長λ１における吸収係数、μａ(λ２)（ｒ）は、ある位置ｒでの波長λ２における吸収係数を示す。ここで、ｒは位置座標を示す。また、εＨｂ０(λ１)は、波長λ１におけるオキシヘモグロビンのモル吸光係数、εＨｂ(λ１)は、波長λ１におけるデオキシヘモグロビンのモル吸光係数を示す。εＨｂ０(λ２)は、波長λ２におけるオキシヘモグロビンのモル吸光係数、εＨｂ(λ２)は、波長λ２におけるデオキシヘモグロビンのモル吸光係数を示す。
【０００５】
ここで、ヘモグロビンの光学特性に関する値、すなわち、εＨｂ０(λ１)、εＨｂ(λ
１)、εＨｂ０(λ２)、εＨｂ(λ２)は既知の値である。したがって、式（１）の形で表
現された血液の酸素飽和度を求めるためには、２波長における吸収係数の比が分かれば良い。
【０００６】
一方、光吸収により被検体内の吸収体から発生する光音響波の初期音圧分布（Ｐ０（ｒ））は次式（２）で表される。
【数２】
2019193732.tif 000004
ここで、Γ（ｒ）は、ある位置ｒでのグリューナイセン（Ｇｒｕｎｅｉｓｅｎ）係数であり、体積膨張係数（β）と音速（ｃ）の二乗との積を定圧比熱（Ｃｐ）で除したものであり、一般的に、位置に依存するが光の波長には依存しない。μａ（ｒ）は、ある位置ｒ
での吸収係数である。Φ（ｒ）は、ある位置ｒでの光量である。ここで言う光量は、吸収体に照射された後の被検体内での光量であり、光フルエンスとも呼ばれる。
式（２）より、ある位置ｒでの初期音圧（Ｐ０（ｒ））は、光音響波を受信した探触子から出力される受信信号（ＰＡ信号）を用いて算出できることが分かる。
【０００７】
式（２）より、２波長における吸収係数比の値は以下のように表現される。
【数３】
2019193732.tif 000005
式（３）より、ある位置ｒでの吸収係数の比の値を得るためには、波長λ１における光量Φλ１と波長λ２におけるΦλ２との比である係数αを求める必要がある。
【０００８】
特許文献１には、波長ごとの光量を含む関係式で表現可能な値を係数αと規定し、係数αと受信信号から初期音圧であるＰ０を求めることで、酸素飽和度が算出できることが開示されている。
【先行技術文献】
【特許文献】
【０００９】
【特許文献１】特開２０１５−２０５１３６号公報
【発明の概要】
【発明が解決しようとする課題】
【００１０】
特許文献１のような方法の場合、光量比がほぼ同一である領域においては、一つの係数αにより酸素飽和度を算出できる。特許文献１においては、例えば測定対象が動脈であれば酸素飽和度の値を９８％などと予測した上で、ユーザーが係数αを任意に設定する。そして、係数αの値を振りながら測定対象の値が予測値になるように調整する。したがってこの方法では、ユーザーの技量や熟練度によって酸素飽和度の値が変わるおそれがある。
【００１１】
また、係数を用いずに精度良く酸素飽和度を算出するためには、被検体内の光量分布を波長ごとに算出する方法もある。しかしこの方法では、被検体内部での散乱や吸収を考慮した複雑な演算が必要である。
【００１２】
本発明は上記課題に鑑みてなされたものである。本発明の目的は、光音響測定において、簡易な手法で安定した濃度情報を取得することにある。
【課題を解決するための手段】
【００１３】
本発明は、以下の構成を採用する。すなわち、
光出射部から出射された、第１の波長を持つ第１の光と、前記第１の波長と異なる第２の波長を持つ第２の光と、が照射された被検体内部の物質から発生する光音響波を探触子が受信して出力される、前記第１の光に対応する第１の信号と、前記第２の光に対応する第２の信号と、を取得する信号取得部と、
前記光出射部から出射された光が照射される照射面における、前記第１の光と前記第２の光の光量比を示す光量比分布を取得する光量比分布取得部と、
前記光量比分布と、前記第１の信号および前記第２の信号を用いて、前記被検体内部における前記物質の三次元の濃度情報を取得する濃度情報取得部と、
を備え、
前記濃度情報取得部は、前記被検体内部の第１の位置における前記光量比として、前記
照射面内において前記第１の位置に対応する位置である第２の位置における前記光量比に基づく値を割り当てる
ことを特徴とする被検体情報取得装置である。
本発明はまた、以下の構成を採用する。すなわち、
光出射部から出射された、第１の波長を持つ第１の光と、前記第１の波長と異なる第２の波長を持つ第２の光と、が照射された被検体内部の物質から発生する光音響波を探触子が受信して出力される、前記第１の光に対応する第１の信号と、前記第２の光に対応する第２の信号と、を取得する信号取得部と、
前記光出射部から出射された光が照射される照射面における、前記第１の光と前記第２の光の光量比を示す光量比分布を取得する光量比分布取得部と、
前記光量比分布と、前記第１の信号および前記第２の信号を用いて、前記被検体内部における前記物質の三次元の濃度情報を取得する濃度情報取得部と、
を備え、
前記濃度情報取得部は、
前記第１の信号から前記第１の光に対応する第１の初期音圧分布を取得し、前記第２の信号から前記第２の光に対応する第２の初期音圧分布を取得し、
前記被検体内部の第１の位置における光量として、前記照射面内において前記第１の位置に対応する位置である第２の位置における光量に基づく値を割り当て、
割り当てられた前記光量によって前記第１の初期音圧分布および前記第２の初期音圧分布を補正し、補正された前記第１の初期音圧分布および前記第２の初期音圧分布を用いて前記濃度情報を取得する
ことを特徴とする被検体情報取得装置である。
本発明はまた、以下の構成を採用する。すなわち、
光出射部から出射された、第１の波長を持つ第１の光と、前記第１の波長と異なる第２の波長を持つ第２の光と、が照射された被検体内部の物質から発生する光音響波を探触子が受信して出力される、前記第１の光に対応する第１の信号と、前記第２の光に対応する第２の信号と、を取得する信号取得部と、
前記光出射部から出射された光が照射される照射面における、前記第１の光と前記第２の光の光量比を示す光量比分布を取得する光量比分布取得部と、
前記光量比分布と、前記第１の信号および前記第２の信号を用いて、前記被検体内部における前記物質の三次元の濃度情報を取得する濃度情報取得部と、
を備えることを特徴とする被検体情報取得装置である。
【発明の効果】
【００１４】
本発明によれば、光音響測定において、簡易な手法で安定した濃度情報を取得することができる。
【図面の簡単な説明】
【００１５】
【図１】実施形態１を適用できる光音響装置の全体構成を示す模式図
【図２】実施形態１を適用できる信号処理部を示す模式図
【図３】実施形態１を適用できる光量比分布取得部を示す模式図
【図４】実施形態１の処理フローの一例を示すフローチャート
【図５】実施形態２を適用できる光音響装置の全体構成を示す模式図
【図６】実施形態２を適用できる信号処理部を示す模式図
【図７】実施形態２を適用できる照射光量取得部を示す模式図
【図８】実施形態２の処理フローの一例を示すフローチャート
【図９】実施形態１に記載の保持部の一例の構成を示す図
【図１０】実施形態１に記載の撮像領域内の各位置に対する補正値を示した図
【図１１】ある乳房の深さに応じた酸素飽和度算出値の変化を示す図
【図１２】実施形態２に記載の撮像領域内の各位置に対する補正値を示した図
【発明を実施するための形態】
【００１６】
以下に図面を参照しつつ、本発明の好適な実施の形態について説明する。ただし、以下に記載されている構成部品の寸法、材質、形状およびそれらの相対配置などは、発明が適用される装置の構成や各種条件により適宜変更されるべきものである。よって、この発明の範囲を以下の記載に限定する趣旨のものではない。
【００１７】
本発明は、被検体から伝搬する音響波を検出し、被検体内部の特性情報を生成し、取得する技術に関する。よって本発明は、被検体情報取得装置またはその制御方法、あるいは被検体情報取得方法や信号処理方法として捉えられる。本発明はまた、これらの方法をＣＰＵやメモリ等のハードウェア資源を備える情報処理装置に実行させるプログラムや、そのプログラムを格納した記憶媒体としても捉えられる。記憶媒体は、コンピュータにより読み取り可能な非一時的な記憶媒体であって良い。
【００１８】
本発明の被検体情報取得装置には、被検体に光（電磁波）を照射することにより被検体内で発生した音響波を受信して、被検体の特性情報を画像データとして取得する光音響効果を利用した光音響イメージング装置を含む。この場合、特性情報とは、光音響波を受信することにより得られる受信信号を用いて生成される、被検体内の複数位置のそれぞれに対応する特性値の情報である。
【００１９】
光音響測定により取得される特性情報は、光エネルギーの吸収率を反映した値である。例えば、単一の波長の光照射によって生じた音響波の発生源、被検体内の初期音圧、あるいは初期音圧から導かれる光エネルギー吸収密度や吸収係数を含む。また、互いに異なる複数の波長により得られる特性情報から、組織を構成する物質の濃度を取得できる。物質濃度としてオキシヘモグロビン濃度とデオキシヘモグロビン濃度を求めることにより、酸素飽和度分布情報を算出できる。また、物質濃度としては、トータルヘモグロビン濃度、グルコース濃度、コラーゲン濃度、メラニン濃度、脂肪や水の体積分率なども求められる。
【００２０】
被検体内の複数の位置の特性情報に基づいて、二次元または三次元の特性情報分布が得られる。分布データは画像データとして生成され得る。特性情報は、数値データとしてではなく、被検体内の各位置の分布情報として求めてもよい。すなわち、初期音圧分布、エネルギー吸収密度分布、吸収係数分布や酸素飽和度分布などの分布情報である。
【００２１】
本発明でいう音響波とは、典型的には超音波であり、音波、音響波と呼ばれる弾性波を含む。トランスデューサ等により音響波から変換された電気信号を音響信号とも呼ぶ。ただし、本明細書における超音波または音響波という記載は、それらの弾性波の波長を限定する意図ではない。光音響効果により発生した音響波は、光音響波または光超音波と呼ばれる。光音響波に由来する電気信号を光音響信号とも呼ぶ。
【００２２】
なお、以下の実施形態における被検体情報取得装置は、人や動物の悪性腫瘍や血管疾患などの診断や化学治療の経過観察などを主な目的とする。よって、被検体としては生体の一部、具体的には人や動物の一部位（乳房、臓器、循環器、消化器、骨、筋肉、脂肪等）の検査対象が想定される。また、検査対象の物質としては、ヘモグロビン、グルコース、また、体内に存在する水、メラニン、コラーゲン、脂質などを含む。さらには、体内に投与されたＩＣＧ（インドシアニン・グリーン）等の造影剤等、光の吸収スペクトルが特徴的な物質であればよい。また、ファントムなどの無生物を被検体としても良い。
【００２３】
以下の実施形態では、被検体情報取得装置として、被検体にパルス光を照射し、光音響効果により被検体から音響波を受信し解析する事により被検体内の光吸収体の分布情報を
取得する光音響装置を取り上げる。被検体として被検者の乳房を想定する。ただし被検体は乳房に限られず、被検者の四肢など他の部位、ヒト以外の動物、無生物、ファントムなども考えられる。以下の実施形態における被検体情報取得装置は、人や動物の悪性腫瘍や血管疾患などの診断や化学治療の経過観察などに好適である。
【００２４】
［実施形態１］
以下、第１の実施形態の被検体情報取得装置の構成及び処理について説明する。図中、同一の構成要素には原則として同一の符号を付して、説明を省略する。
【００２５】
（全体的な装置構成）
図１は本実施形態の光音響装置の構成を示す模式図である。光音響装置は、少なくとも、光源００１、光音響波を受信する変換素子００６を備える探触子０１１、変換素子００６から出力される受信信号を用いて被検体００４内の特性値情報を取得する信号処理部０１０、照射される異なる波長の光量比の分布を算出する光量比分布取得部００７を備える。
【００２６】
光源００１からの光は光導波部００２により光出射部００３に導かれ、光出射部００３から出射される。酸素飽和度を測定する場合、光源００１は、互いに異なる波長を有する複数のパルス光を別々のタイミングで出力する。光出射部００３から出射された照射光は被検体００４に照射される。照射光は、被検体００４内を伝搬し、測定対象部位としての光吸収体００５に到達する。ここでは光吸収体００５は、オキシヘモグロビン、デオキシヘモグロビン、血液、血管、血管を多く含む腫瘍、などである。光吸収体００５が互いに異なる波長それぞれの光のエネルギーを吸収するごとに、光音響波が発生する。発生した光音響波は、被検体００４内を伝搬し変換素子００６に到達し、受信信号に変換される。
【００２７】
光量比分布取得部００７は、被検体００４に照射された複数の互いに異なる波長の光の、光量比の分布を算出する。被検体００４の光音響波を効率よく取得するため、被検体００４を保持する保持部００８が変換素子００６と被検体００４の間に設置され、保持部００８と変換素子００６の間には音響整合液００９が満たされている。複数の変換素子００６の各々は、光音響波を受信して時系列のアナログ受信信号を出力する。
【００２８】
信号取得部０２２は、出力されたアナログ受信信号に対して、増幅器によるアナログ信号の増幅やＡＤ変換器によるデジタル変換を行う。信号処理部０１０には、照射されたパルス光分のデジタル受信信号（以下、受信信号）が順次入力される。受信信号は、いったんメモリに保存し、必要に応じて信号処理部０１０が読み出すようにしても良い。信号処理部０１０は、入力された受信信号を用いて、被検体００４内の特性値情報を生成する。なお、探触子０１１内にアナログ信号の増幅器を設けてもよい。これにより、少なくとも、第１の光に対応する第１の信号と、第２の光に対応する第２の信号と、が取得できる
なお、受信信号の符号化を行う場合、光源００１が、互いに異なる複数波長を有する複数のパルス光を同時に出力するタイミングが存在してもよい。光源００１は、互いに異なる波長のそれぞれを異なる符号系列に対応する強度変調光を出力することにより、受信信号の符号化を行うことができる。この場合、信号処理部０１０が、符号化された受信信号に対して、上記符号系列に関する情報に基づいた復号化処理を行うことにより、複数波長のそれぞれに対応する受信信号を復号することができる。
【００２９】
なお、光音響装置が光音響顕微鏡等である場合は、探触子０１１が備える変換素子００６は１つでもよい。ただし、光音響装置が乳房等の被検体００４を検査する被検体情報取得装置である場合は、探触子０１１が複数の変換素子００６を備えることが好ましい。特に、被検体００４内の特性値情報の算出精度を高めるには、複数の変換素子００６を、球状、半球状あるいは円柱状に、三次元的に密に配置することが好ましい。ただし本発明は
、ハンドヘルド型などの光音響装置にも適用可能である。
【００３０】
（光源００１）
光源００１としては、ナノ秒からマイクロ秒オーダーのパルス光を発生可能なパルス光源が好ましい。具体的なパルス幅としては、１〜１００ナノ秒程度が好ましい。また、波長としては４００ｎｍから１６００ｎｍ程度の範囲が好適に使われる。特に、生体の深部をイメージングする際には、生体を構成する成分のうち検査対象とする特定の物質（例えばヘモグロビン）に吸収され、かつ、他の物質において吸収が少ない波長の光を用いる。具体的には、７００ｎｍから１１００ｎｍの領域が好ましい。一方、生体表面近傍の血管を高解像度でイメージングする際は可視光領域を用いることが好ましい。ただし、テラヘルツ波、マイクロ波、ラジオ波領域の使用も可能である。
【００３１】
具体的な光源００１としては、高出力が可能となるレーザーが好ましい。本実施形態では複数波長の光を用いるため、発振する波長の変換が可能なレーザーがより好ましい。ただし、互いに異なる波長の光を発振する複数台のレーザーを、発振切り替えを行いながら用いることも可能である。複数台のレーザーを用いた場合もそれらをまとめて光源００１として表現する。光源は、少なくとも、第１の波長を持つ第１の光と、前記第１の波長と異なる第２の波長を持つ第２の光とを出力できる。また、３つ以上の波長の光を出射可能であってもよい。
【００３２】
レーザーとしては、固体レーザー、ガスレーザー、色素レーザー、半導体レーザーなど様々なレーザーを使用できる。特に、Ｎｄ：ＹＡＧレーザーやアレクサンドライトレーザーなどのパルスレーザーが好ましい。また、Ｎｄ：ＹＡＧレーザー光を励起光とするＴｉ：ｓａレーザーやＯＰＯ（Optical Parametric Oscillators）レーザーを用いてもよい。また、レーザーの代わりに発光ダイオードやフラッシュランプなども利用できる。
【００３３】
（光導波部００２、光出射部００３）
光導波部００２と光出射部００３は、光源００１から被検体００４に光を伝達する。光導波部００２及び光出射部００３としては、レンズ、ミラー、光ファイバ等の光学素子を使用できる。ただし、光源００１から直接被検体００４に光照射することも可能である。また、乳房等を検査対象とする生体情報取得装置においては、光出射部００３はレンズ等によりビーム光の径を広げて照射することが好ましい。一方、光音響顕微鏡においては、解像度を上げるためにレンズ等でビーム光の径をフォーカスして照射することが好ましい。また、光出射部００３を被検体００４に対して移動してもよく、これにより被検体００４の広い範囲のイメージングが可能になる。
【００３４】
（探触子０１１）
探触子０１１は１つ以上の変換素子００６を備える。変換素子００６は、チタン酸ジルコン酸鉛（ＰＺＴ）などの圧電現象を用いた圧電素子、光の共振を用いた変換素子００６、ＣＭＵＴ等の静電容量型の変換素子００６など、音響波を受信して電気信号に変換できるものであればどのようなタイプでもよい。探触子０１１が複数の変換素子００６を備える場合は、１Ｄアレイ、１．５Ｄアレイ、１．７５Ｄアレイ、２Ｄアレイ、アーク型アレイ、半球型アレイと呼ばれる、平面又は曲面内に複数の素子を並べる配置が好ましい。
【００３５】
また、乳房等を検査する生体情報取得装置の場合は、広い範囲のイメージングを行うため、探触子０１１が、被検体００４に対して機械的に移動可能であることが好ましい。ハンドヘルド型の探触子０１１の場合は、ユーザーが探触子０１１を把持して動かしてもよい。光音響顕微鏡の場合は、探触子０１１はフォーカス型の探触子０１１とすることが好ましく、探触子０１１は被検体００４の表面に沿って機械的に移動することが好ましい。
【００３６】
また、照射光の照射位置と探触子０１１とは同期して移動することが好ましい。本実施形態では、光出射部００３は半球形の探触子０１１の中央部に設置されている。探触子０１１は、移動制御部０１３の制御に従い、移動機構０１２によって二次元平面上を移動する。移動機構０１２は、モータやソレノイドなどの動力源と、ＸＹステージや直動ボールねじなどの動作部を備える。これにより、被検体００４の表面と、それぞれの変換素子００６との相対位置を、位置に制御できる。
【００３７】
（保持部００８）
検査対象の情報を正確に取得するために、保持部００８を設置し撮像中の被検体００４を固定するのが好ましい。保持部００８の材料としては、照射光や光音響波を伝搬するものが好ましい。例えば、アクリル、ポリカーボネート、ＰＥＴ等の透明の樹脂や、シリコーンゴム、天然ゴム、ウレタンゴム等のゴムシート、また、広範囲に均一な開口が設けられているメッシュ素材などの材料が使用できる。薄いシート状の部材と、金属等で形成された強度の高い骨格部材とを組み合わせてもよい。
【００３８】
また、光の照射量の見積もりや、光吸収体００５からの音響波の変化の推定を容易にするためには、被検体００４を任意の形状に固定できることが好ましい。これにより検査対象の情報をより正確に取得できる。たとえば、乳房等の凸形状の被検体００４に対しては、探触子０１１側に突き出たカップ形状の保持部００８が好ましい。また、手のひらや足裏等の平面形状の被検体００４に対しては、平板形状の保持部００８が好ましい。本実施形態では、移動機構０１２の走査面と並行な平面ＰＥＴ板を保持部００８として使用しているため、変換素子００６と光出射部００３は被検体００４と常に同じ距離に保たれる。
【００３９】
なお、保持部００８は必ず必要というものではない。また、被検体００４の形状情報は、保持部００８を用いる以外の手法でも取得可能である。別の手法として例えば、カメラ撮影による画像処理や、被検体００４の皮膚や皮膚表面付近の血管からの音響情報から被検体００４形状を推定する手法や、超音波を送受信する方法がある。
【００４０】
（音響整合液００９）
音響整合液００９は、被検体００４（または保持部００８）と探触子０１１を音響的に整合させる。したがって音響整合液００９としては、音響波を伝搬し、かつ探触子０１１の走査を妨げないものが好ましい。例えば、水、ＤＩＤＳ（セバシン酸ジイソデシルエステル）、ＰＥＧ（ポリエチレングリコール）、シリコーンオイル、ひまし油などの液体が挙げられる。本実施形態では水を使用している。また、保持部００８と被検体００４の間にも、水やジェルなどの音響整合材を配置してよい。
【００４１】
（表示部０１４）
表示部０１４である表示部０１４には、ＬＣＤ（Liquid Crystal Display）やＣＲＴ（Cathode Ray Tube）、有機ＥＬディスプレイ等のディスプレイを利用できる。なお、表示部０１４は、本実施形態の光音響装置が備える構成とはせずに、別に用意して光音響装置に接続しても良い。
【００４２】
（情報取得部０１５の内部構成）
次に、本実施形態の信号処理部０１０内の構成を、図２を用いて説明する。本実施形態の信号処理部０１０は、情報取得部０１５、表示制御部０１６、係数決定部０１７を備える。
【００４３】
情報取得部０１５は、信号取得部０２２から出力される受信信号を用いて、被検体００４内の特性値情報を位置毎に取得する。具体的には、変換素子００６毎の時系列の受信信号を用いて画像再構成を行うことにより、二次元又は三次元の空間座標上の位置に対応す
る特性値のデータ（分布情報）を求める。再構成の単位領域は、ピクセルまたはボクセルと呼ばれる。画像再構成手法としては、ＦＢＰ（Filtered Back Projection）やタイムリバーサル法、モデルベース法、フーリエ変換法等の公知の再構成手法を利用できる。また、超音波イメージングで使われている整相加算（Delay and Sum）処理を用いてもよい。
【００４４】
また、光フォーカス型の光音響顕微鏡や、フォーカス型の探触子０１１を用いた光音響顕微鏡の場合は、画像再構成処理を行わずに分布情報を生成することもできる。具体的には、走査機構により、探触子０１１と光照射スポットとを被検体００４に対して相対移動させる。探触子０１１は複数の走査位置で光音響波を受信する。そして、情報取得部０１５は、得られた受信信号を時間変化に対して包絡線検波した後、受信信号における時間軸方向を奥行き方向に変換して、空間座標上にプロットする。これを走査位置毎に行うことにより、分布情報を構成できる。
【００４５】
表示制御部０１６は、情報取得部０１５により生成された特性情報や分布情報を基に、表示部０１４に表示するための画像データを生成する。具体的には、分布情報を基に、輝度変換、歪補正、注目領域の切り出し、血管抽出処理、動脈・静脈分離処理、対数圧縮処理などの画像処理を行う。
【００４６】
係数決定部０１７は、情報取得部０１５が受信信号を基に特性値情報を求める際に用いる係数βを決定する。係数βは光量比分布取得部００７から出力された、異なる波長の光量比情報を用いて算出される。係数βについては後に詳述する。
【００４７】
（光量比分布取得部００７の内部構成）
次に、本実施形態の光量比分布取得部００７内の構成を、図３を用いて説明する。本実施形態の光量比分布取得部００７は、照射面決定部０１８、光量プロファイル取得部０１９、光量分布取得部０２０を備える。
【００４８】
照射面決定部０１８は、撮像領域における波長間の光量比を算出する光量比取得面を決定する。この光量比取得面として本実施形態では、被検体表面における光の照射面を用いる。以降、照射面が光量比取得面となる場合について検討する。保持部００８によって被検体００４が保持される場合は、光量比取得面として保持部００８と被検体００４の接触面を利用できる。特に本実施形態では、探触子０１１が移動する走査平面と平行な二次元平面において被検体００４が保持されるため、この二次元平面が利用される。ただし、ビームプロファイル中の任意の距離における切断面を照射面とする場合でも、当該切断面に光が照射されていることには変わりないため、当該切断面を照射面と呼んで構わない。
【００４９】
なお、保持部００８が交換可能な場合、設置した保持部００８に応じて照射面を自動、もしくは手動で変更できることが好ましい。自動で変更する場合は、保持部００８の側にＩＤタグ等の認識システムを搭載し、装置本体の側にタグ情報の読み取り機構を設けることが好ましい。この場合、照射面決定部０１８は、読み取り機構からの情報に基づいて光量比取得面を決定する。また、図１における乳房の側部のように、被検体００４の部位によっては、保持部００８に接触しない場合がある。その場合は被検体００４表面位置を計測するか、ユーザーが任意断面を入力することで、照射面を決定できる。
【００５０】
光量プロファイル取得部０１９は、光出射部００３から、照射面決定部０１８により規定された照射面に照射される光の、光量プロファイルを算出する。光量プロファイルは、任意の内における光の強度の分布を表し、波長ごとに求められる。算出には、あらかじめ取得されている光量プロファイル情報と、光出射部００３と任意の面の距離関係と、が利用される。任意の面は、被検体表面の照射面でも良いし、ビームプロファイルを出射部から任意の距離で切断した断面であっても良い。
【００５１】
光量プロファイル算出方法の一例を述べる。まず、あらかじめ、複数の異なる距離での光量プロファイル情報を取得して、メモリ等の記憶手段（不図示）に保存しておく。そして、異なる距離での光量プロファイル情報に基づいて三次元的な変化状況を見積り、補間を行う。そして、光出射部００３から照射面までの距離関係を用いて三次元的な光量プロファイルを参照することで、照射面における光量プロファイルを取得する。あるいは、光量プロファイル情報として、光出射部００３から出射された光の各距離における面内強度分布を示す三次元的なビームプロファイルを、メモリに保存しておいてもよい。
【００５２】
本実施形態では、光出射部００３と光量比取得面（この場合は照射面）との間には音響整合液００９のみが存在すると仮定している。光量プロファイル取得部０１９は、波長ごとに、音響整合液００９による光量の減衰を考慮して光量プロファイルを算出する。なお、音響整合液００９が存在することを前提とした光量プロファイル情報を、メモリ内に保存しておいても良い。照射面が平面でない場合でも、同様の手法で光量比を取得できる。また、特に本実施形態のように照射面と光出射部００３の位置が一定値で規定される場合は、あらかじめ見積もられた光量プロファイルをメモリに保管しておくのも効果的である。
【００５３】
光量分布取得部０２０は、光量プロファイル取得部０１９で算出された光量プロファイルと、光源００１から照射された光量値と、光が照射された時の移動機構０１２の座標情報より、被検体００４の情報取得時に照射面に照射された光量の分布を算出する。なお、異なる波長で算出された各々の光量分布の比を取得することで、光量比分布が算出される。
【００５４】
（信号処理部０１０における処理）
本実施形態では情報取得部０１５が特性値情報として、光音響波の音圧に関する情報と、酸素飽和度に関する情報と、を少なくとも求める。ここで、「酸素飽和度」とは、本明細書における「濃度に関する情報」のうちの１つであり、赤血球中のヘモグロビンのうち、酸素と結合しているヘモグロビンの割合を示す。したがって本実施形態における情報取得部０１５は、濃度情報取得部としての機能を有する。
【００５５】
式（１）によって酸素飽和度を求めるためには、式（３）に示すような複数の波長（少なくとも、互いに異なる２波長）における吸収係数の比が必要となる。式（３）における初期音圧（Ｐ０）は、実際に被検体００４内で発生した光音響波の発生圧力の相対値を示している。一般に、光量分布情報Φ（ｒ）は、無限媒体の拡散方程式の解析解から、簡易的に以下の式（４）で表される。
【数４】
2019193732.tif 000006
【００５６】
ここでΦ０は単位面積当たりの光照射エネルギーである。μｅｆｆは有効減衰係数（effective attenuation coefficient）である。被検体００４背景の等価散乱係数をμｓ’
（reduced scattering coefficient）、吸収係数をμａとすると、以下のように表現される。
μｅｆｆ＝｛３μａ（μａ＋μｓ’）｝１／２
ｄ（ｒ）は、被検体００４内のある位置ｒと、被検体００４表面における光照射領域との距離である。なお、ここでは、被検体００４の有効減衰係数の位置による変化は無視できるほど小さく、被検体００４の有効減衰係数は被検体００４内の位置に依存しない係数であると仮定している。
【００５７】
これより、式（３）を変形すると以下の式（５）のようになる。
【数５】
2019193732.tif 000007
この式（５）が表す吸収係数比を式（１）に代入することで酸素飽和度を算出できる。すなわち、酸素飽和度は異なる２波長の初期音圧Ｐ０と、光照射エネルギーΦの比から算出される。光照射エネルギーΦの比は、光量比とも呼ばれる、係数βとして係数決定部で決定される値であり、下記の（Ａ）のように表現される。ここで、第１の光に対応する初期音圧分布を第１の初期音圧分布、第２の光に対応する初期音圧分布を第２の初期音圧分布と呼んでもよい。
Φ（λ２）／Φ（λ１）  …（Ａ）
この値は、初期の単位面積当たりの光照射エネルギー比（初期光照射エネルギー比）と、各波長での減衰から求められる。具体的には、下記の値（Ｂ）と値（Ｃ）の積である。
Φ０（λ２）／Φ０（λ１）  …（Ｂ）
ｅｘｐ（−μｅｆｆ（λ２）・ｄ（ｒ））／ｅｘｐ（−μｅｆｆ（λ１）・ｄ（ｒ））
…（Ｃ）
【００５８】
前者の値（Ｂ）「Φ０（λ２）／Φ０（λ１）」は、光量比分布取得部００７から精度よく算出できる。しかし、後者の値（Ｃ）「ｅｘｐ（−μｅｆｆ（λ２）・ｄ（ｒ））／ｅｘｐ（−μｅｆｆ（λ１）・ｄ（ｒ））」は、生体内の光学特性が関係するため、算出が難しい。特に脂肪層が多い乳房のような被検体００４に比べ、手脚のような複数の組織が入り組んでいる被検体００４における光学特性を算出するのは困難である。
【００５９】
本実施形態では係数βを求める際に、上記の値（Ｃ）「ｅｘｐ（−μｅｆｆ（λ２）・ｄ（ｒ））／ｅｘｐ（−μｅｆｆ（λ１）・ｄ（ｒ））」として簡易的な値を設定する。ここでは、値（Ｃ）として、被検体００４の深さ方向ｄ（ｒ）に関して、深さによらず一定の値を使用する。特にここでは、値（Ｃ）を「１」とする。このとき係数βは、初期の光量比から変化せず、「Φ０（λ２）／Φ０（λ１）」となる。このように、入射位置からの距離にかかわらず略一定の減衰率比および光量比を用いることにより、簡易的に、かつユーザーの技量にかかわらず酸素飽和度を算出できる。
【００６０】
なお、深さ方向とは、光の照射位置を起点として光の照射方向と略同一の方向を指してもよいし、光の照射位置を起点として被検体に略垂直な方向を指しても良い。本実施形態では、照射面における光量比分布を、そのまま被検体内部の深さ方向において対応する位置に適用する。言い換えると、本実施形態では、照射面内の任意の入射位置における光量比を、当該入射位置を起点とした光の進行方向に位置する被検体の領域、または、当該入射位置を起点とした被検体表面に略垂直な方向に位置する被検体の領域に割り当てる。照射面（光量比算出面）において照射位置を第２の位置とし、第２の位置に対応する位置を第１の位置とする。本実施形態では、第１の位置は第２の位置から深さ方向に存在する。
【００６１】
図１０を用いて説明する。ここで、被検体に円筒状の撮像領域を設定する。また、光出射部００３からの光が照射される照射面Ｓを光量比取得面とし、照射面Ｓにおける便宜上の小領域として、第１部分Ｓ１、第２部分Ｓ２、第３部分Ｓ３を仮定する。このとき、照射面Ｓから、第１〜第３部分の各々を起点として照射方向に伸びる（この例では照射面に略垂直な）撮像領域を、領域Ｒ１〜Ｒ３とする。また、光量比分布取得部００７にて算出された、第１〜第３の各部分での初期光照射エネルギー比（Ｂ）「Φ０（λ２）／Φ０（
λ１）」の値は、Ａ１〜Ａ３であるものとする。このとき、領域Ｒ１〜Ｒ３内の全ての深さにおいて、係数βの値もそれぞれ、Ａ１〜Ａ３と設定する。
【００６２】
このように、係数βとして深さ方向の位置に関わらず一定の値を適用すると、特性情報の取得対象となる被検体内部の位置が光照射面から深く（遠く）なればなるほど、酸素飽和度の誤差も大きくなることが懸念される。しかし、被検体における浅い箇所の酸素飽和度は比較的安定して計測できる。
【００６３】
図１１を用いて、深さに応じた誤差について説明する。ここで、被検体００４は乳房である。また、照射される２波長はそれぞれ、λ１＝７５６ｎｍ、λ２＝７９７ｎｍとする。また、各波長の光に対する被検体の代表的な有効減衰係数μｅｆｆ（λ１）、μｅｆｆ（λ２）を適宜設定する。図１１は、動脈を想定した吸収係数９８％の吸収体の酸素飽和度ＳＯ２の算出値を、深さｄ（ｒ）ごとに示す。通常、一本の動脈における酸素飽和度値は略同じ値になる。
【００６４】
この図より、皮膚表面から深さ数ｍｍの位置での酸素飽和度の誤差は１％程度であることが分かった。また、深さが５ｍｍでは誤差２〜３％程度、深さ１０ｍｍでは誤差５〜６％程度、深さ２０ｍｍでは誤差１０〜１１％である。このことから、被検体表面に照射される光量比を、そのまま照射面から垂直方向における係数βの取得に利用しても、深さが５ｍｍ以下の浅い箇所の酸素飽和度をほとんど誤差なく求められることが分かった。また、深さが１０ｍｍ以下程度までなら、動脈・静脈の識別に問題なく利用できる程度の算出精度を期待できることが分かった。また、酸素飽和度９８％の動脈について、深さ１８ｍｍ以下程度であれば、酸素飽和度９０％台の値を導出できることが分かった。
【００６５】
（処理フロー）
図４を参照して、信号処理部０１０が酸素飽和度分布を求める処理フローを説明する。図４のフローは、信号処理部０１０がデジタル信号を取得した状態からスタートする。すなわちフロー開始前に、光出射部００３が、波長ごとに、光源００１からの光を出射し、変換素子００６が、被検体００４から発生する光音響波を受信してアナログの受信信号に変換し、信号取得部０２２が、照射された光の波長ごとに順次、受信信号に対するＡＤ変換や増幅などを行う。変換後のデジタル信号は、メモリに保存されたのち信号処理部０１０に入力されるか、信号取得部０２２から信号処理部０１０に直接入力される。
【００６６】
ステップＳ１０１では、情報取得部０１５が、入力された受信信号を用いて、第１波長λ１における音圧（Ｐ０λ１（ｒ））分布情報と、第２波長λ２における音圧（Ｐ０λ２（ｒ））分布情報と、を取得する。
【００６７】
ステップＳ１０２では、表示制御部０１６が、情報取得部０１５により生成された複数波長分の音圧分布情報のうち少なくとも１波長分の音圧分布情報を基に画像処理を行い、音圧分布を示す画像を表示部０１４に表示する。このとき表示制御部０１６は、表示部０１４に、音圧分布を示す画像に基づいて生成された画像（例えば、動脈あるいは静脈など特定の血管のみを表示する画像、２波長間の画像の差や比の画像、擬似的な酸素飽和度分布画像など）を表示しても良い。
【００６８】
ステップＳ１０３では、光量比分布取得部００７が、任意の照射面における光照射領域の２波長の照射光量比を算出する。任意の照射面は、典型的には被検体表面であり、保持部００８を用いる場合は当該保持部００８の形状から求められる。その他、上述した各種の方法で表面形状情報を取得できる。また、被検体００４の形状にかかわらず、ユーザーが任意の位置に光量比取得面を設定しても構わない。
【００６９】
ステップＳ１０４では、情報取得部０１５が、光量比分布取得部００７において算出された任意の照射面における照射光量比と、波長λ１における音圧（Ｐ０λ１（ｒ））分布情報と、波長λ２における音圧（Ｐ０λ２（ｒ））分布情報と、既知のオキシヘモグロビン・デオキシヘモグロビンのスペクトル情報を用いて、酸素飽和度の分布情報を生成する。基地の情報は、メモリに保存しておき必要に応じて読み出すことが好ましい。
【００７０】
ステップＳ１０５では、表示制御部０１６が、情報取得部０１５により生成された酸素飽和度の分布情報を基に画像データを生成し、表示部０１４に表示する。このように、ユーザーは、表示された音圧分布を示す画像あるいは音圧分布を示す画像に基づいて生成された画像から、酸素飽和度分布の画像を取得できる。
【００７１】
以上説明したように、本実施形態では、酸素飽和度等の濃度に関する情報を求める際に、被検体００４表面に代表される規定面に照射される異なる波長の光量比の分布を算出し用いる。すなわち、特許文献１に示されたようなユーザーの技量によって変化する値を用いる必要がない。また、光量分布を算出するための散乱や吸収を考慮した複雑な演算を行う必要がない。その結果、簡易的な計算で、安定して、酸素飽和度などの濃度情報の分布を取得できる。
【００７２】
［変形例１］
本変形例では、上で図１１を用いて説明したような、被検体内部の比較的深い部分、すなわち照射位置からの光の伝搬距離が所定の閾値以上の長さの部分における酸素飽和度の算出誤差について検討する。深さに関する所定の閾値は、誤差がユーザーの許容範囲を超えるかどうかに応じて決定される。所定の閾値以上の範囲での酸素飽和度を取得する場合、値（Ｃ）「ｅｘｐ（−μｅｆｆ（λ２）・ｄ（ｒ））／ｅｘｐ（−μｅｆｆ（λ１）・ｄ（ｒ））」として照射面における値を割り当てることはできない。この場合、被検体に合わせた光学係数を元に値（Ｃ）を算出する必要がある。よって、信号処理部０１０としては、かかる減衰補正処理を実行できる装置を用いる。
【００７３】
また、ユーザーが、簡易的な算出方法と一般的な算出方法を選択可能にすることも好ましい。ユーザーによる選択には、例えば情報処理装置のマウスやキーボード等のユーザインタフェースを利用できる。選択は予め行われてもよいし、表示部０１４に表示された画像を参照してユーザーが動的に指定してもよい。
【００７４】
［変形例２−１］
また、あらかじめ設定された、許容できる算出誤差を示す所定の値に基づいて、信号処理部０１０が、一般的な算出方法を使うか簡易的な算出方法を使うかを自動的に決定してもよい。一般的な算出方法の場合、信号処理部が値（Ｃ）「ｅｘｐ（−μｅｆｆ（λ２）・ｄ（ｒ））／ｅｘｐ（−μｅｆｆ（λ１）・ｄ（ｒ））」を計算し係数βに反映させる。また、簡易的な算出方法の場合、Φ０（λ２）／Φ０（λ１）を係数βとする。例えば、ユーザーが、動脈・静脈を識別するために、許容算出誤差を５％と設定した場合、信号処理部０１０は、許容算出誤差の値と、図１１のグラフが示す情報を比較して、深さ１０ｍｍの位置を境界として設定する。すなわち、深さ１０ｍｍより深い領域の酸素飽和度を算出する場合には値（Ｃ）「ｅｘｐ（−μｅｆｆ（λ２）・ｄ（ｒ））／ｅｘｐ（−μｅｆｆ（λ１）・ｄ（ｒ））」を計算して係数βに反映させる。
【００７５】
さらに、表示部０１４に表示される画像中に、簡易的な手法で算出された部分と一般的な手法で算出された部分が混在する場合、どちらの手法で求められたかをユーザーが判別可能とすることが好ましい。例えば、色度や明度、彩度、解像度を変えたり、文字情報を提示したりしても良い。
【００７６】
［変形例２−２］
上記のような、照射位置からの深さに応じた算出方法の使い分けを、光音響測定と、光音響測定以外のモダリティとを組み合わせた測定方法に適用してもよい。例えば、超音波測定と光音響測定を組み合わせて乳がん検診を行う場合について述べる。ここで、本発明の被検体情報取得装置が、光音響測定と超音波エコー測定の両方を実施可能だとする。かかる兼用の構成は、探触子０１１の変換素子００６が被検体００４に超音波を送信し、被検体表面や内部で反射したエコー波を受信することにより実現できる。
【００７７】
被検体情報取得装置は、まず超音波エコー測定を行い、被検体内部の構造情報を示す機能画像を作成し、表示部０１４に表示する。超音波画像を見たユーザーは、マウス等のユーザインタフェースを用いて、特に詳しく観察したい関心領域（ＲＯＩ：Region of Interest）を指定して酸素飽和度画像を要求する。続いて信号処理部０１０は、当該関心領域が、簡易的な手法で酸素飽和度を取得できる範囲か、一般的な手法を使う必要がある範囲かを判断する。そして、決定された手法を用いて酸素飽和度を算出する。なお、関心領域が両方の範囲にまたがっている場合は、関心領域を分割しても良い。
【００７８】
［変形例３］
また、被検体００４を複数回計測する場合は、各回での計測条件を同じにすることで、深い箇所の濃度情報の再現性が高くなる。計測条件としては、光出射部００３と被検体００４との相対位置関係や傾きなどがある。これらの位置関係を同じにすることで、光照射面から撮像位置までの波長ごとの光学特性や経路が再現され、深い位置における酸素飽和度の誤差も再現される。図１１を用いて説明すると、乳房を複数回計測したとき、乳房の設置位置を同じにできれば、深さ１０ｍｍに設置されていた光吸収体００５は再び深さ１０ｍｍに設置されるため、酸素飽和後のずれ量も５％程度の同じ値に再現される。そのため、酸素飽和度としては誤差が生じても、再計測時の誤差は同じになり、濃度に起因する情報の再現性は確保される。この方法は、各回の計測の間で濃度情報の変化を確認したい場合に特に有効である。
【００７９】
上述した、複数回計測時の濃度情報の変化の再現性を高めるために、以前の計測条件を再現させることが望ましい。たとえば、保持部００８を用いて被検体００４の形状を規定することで、光出射部００３と被検体００４との距離が再現される。また、探触子０１１の走査方向の位置再現性を確保するために、光音響装置に撮像時の被検体００４の設置状況を撮影するカメラと設置状況を表示するディスプレイを設けるのが好ましい。また、設置位置の再現性を高めるため、以前の設置状況を撮像した画像を表示するのも好ましい。また、グリッドや、以前の被検体の輪郭や特徴点を抽出した画像を表字するのが好ましい。これらの施策により、ユーザーは前回と同じ位置に被検体００４を設置しやすくなる。
【００８０】
［変形例４］
光音響装置に、被検体００４の走査方向の位置を規定される部材を設けることも好ましい。たとえば図９は、被検体００４が手や足である場合に好適な例を示す。図９（ａ）の保持部００８には、被検体００４側に突き出た突き当て部０２３が設けられている。被検者は、図９（ｂ）のように、指の間に突き当て部０２３を押し付けることで被検体００４を固定する。これにより、被検体００４の設置位置を三次元的に規定できる。
【００８１】
［変形例５］
ここでは、移動機構０１２の走査によって探触子０１１が移動しながら光音響測定を行う場合について説明する。図１に示すように、光出射部００３も探触子０１１と同時に移動する。本変形例では、光出射部０１１が波長λ１の光パルスと波長λ２の光パルスを交互に出射するものとする。信号取得部は、光の出射ごとに、探触子が受信した受信信号を波長および照射位置と関連付けてメモリに保存する。照射位置は、装置に固有のパラメー
タで表現されてもよいし、被検体との相対的な位置関係で表現されても良い。
【００８２】
本変形例の光量比分布取得部００７は、照射面（光量比算出面）として、探触子０１１が移動する各位置における照射面を足し合わせて仮想的な照射面を作成する。そして光出射部００３の波長ごとの光量プロファイルと、探触子０１１の軌跡に基づいて、仮想的な光量比分布を算出する。本変形例のように波長ごとに交互に光が出射されている場合、光量分布取得部０２０は、１パルスおきに光量を足し合わせて仮想的な光量分布を算出する。そして、仮想的な照射面における光量分布を深さ方向に割り当てることにより、簡易的に酸素飽和度を算出する。
【００８３】
［実施形態２］
実施形態２について説明する。本実施形態の光音響装置は、実施形態１の光音響装置とほぼ同様の構成であるため、相違点を中心として説明する。構成上の相違点としては、図５に記載のように、光量比分布取得部００７の位置に照射光量取得部０２１を設けたことと、図６の記載のように、信号処理部０１０が初期音圧取得部０２４、初期音圧補正部０２５、濃度情報取得部０２６、表示制御部０１６から構成されている点である。
【００８４】
（照射光量取得部０２１の内部構成）
本実施形態の照射光量取得部０２１内の構成を、図７を用いて説明する。本実施形態の照射光量取得部０２１は、照射面決定部０１８、光量プロファイル取得部０１９、光量分布取得部０２０を備える。
【００８５】
照射面決定部０１８は、撮像領域における各波長の光の光量分布を算出する、光量分布取得面を決定する。光量分布取得面としては照射面、すなわち光が照射される被検体表面が好ましいが、それ以外の任意の面でも構わない。保持部００８によって被検体００４が保持される場合は、保持部００８と被検体００４の接触面が好適である。本実施形態では、探触子０１１の走査平面と平行な平面形状に被検体００４が保持されるため、この二次元平面を照射面（光量分布取得面）とする。なお、保持部００８が交換可能な場合は、設置した保持部００８によって照射面を自動、もしくは手動で変更できることが好ましい。自動でできる場合は保持部００８にＩＤタグ等の認識システムを搭載し、この情報を読み取ることで照射面決定部０１８が光量分布算出面を決定することも好ましい。
【００８６】
光量プロファイル取得部０１９は、照射面決定部０１８により規定された照射面に光出射部００３から照射される光の、光量プロファイルを算出する。光量プロファイル取得部０１９は、あらかじめ取得されている光量プロファイル情報をもとに、光出射部００３と照射面の距離関係から、光量プロファイルを算出する。たとえば異なる距離であらかじめ取得された光量プロファイルから三次元的にプロファイルの変化状況を見積ることで、任意の距離関係における照射面での光量プロファイルを算出できる。また、図示したように照射面と光出射部００３の位置が一定の場合は、あらかじめ見積もられてメモリに保存された光量プロファイルを読み出しても良い。
【００８７】
光量分布取得部０２０は、光量プロファイル取得部０１９が取得した光量プロファイルと、光源００１から照射された光量値と、光が照射された時の移動機構０１２の座標情報より、被検体情報を取得するために照射面に照射された光量の分布を算出する。なお、本実施形態では、照射光は、光出射部００３と照射面の間の音響整合液００９を通過する。そこで、音響整合液００９による光の減衰を考慮に入れて光量分布を算出する。
【００８８】
（信号処理部０１０における処理）
本実施形態では、信号処理部０１０が特性値情報として、光音響波の音圧に関する情報と、算出した照射光量から補正された補正音圧に関する情報と、酸素飽和度に関する情報
と、を少なくとも求める。
【００８９】
初期音圧取得部０２４は、実施形態１で説明したような初期音圧（Ｐ０（λ１））を算出する。初期音圧補正部０２５は、式（６）で記載したように、初期音圧（Ｐ０（λ１））を照射光量取得部０２１で算出された照射光量（Φ（λ１））で補正することで、補正初期音圧（Ｐ１（λ１））が算出する。同様の処理を波長λ２についても行う。
【数６】
2019193732.tif 000008
【００９０】
本実施形態では、異なる２波長の初期音圧（Ｐ０（λ１），Ｐ０（λ２））を、各波長の光量分布で補正し、規格化する。これにより、初期音圧の照射光量分布によるムラが補正される。なお、式（６）において、ｅｘｐ（−μｅｆｆ（λ１）・ｄ（ｒ））として、被検体００４の深さ方向ｄ（ｒ）に関わらず一定の値を使用する。また、第２波長についても同様とする。結果、本実施形態では、ｅｘｐ（−μｅｆｆ（λ１）・ｄ（ｒ））、および、ｅｘｐ（−μｅｆｆ（λ２）・ｄ（ｒ））を、「１」と規定する。
【００９１】
図１２を用いて説明する。被検体に円筒状の撮像領域を設定し、光出射部００３からの光が照射される照射面Ｓを光量比取得面とする。照射面Ｓにおける便宜上の小領域として、第１部分Ｓ１、第２部分Ｓ２、第３部分Ｓ３を仮定し、照射面Ｓから、第１〜第３部分の各々を起点とした照射面に垂直な撮像領域を、領域Ｒ１〜Ｒ３とする。そして、照射光量取得部０２１にて算出された、Ｓ１〜Ｓ３でのΦ０（λ１）を、Ｂ１〜Ｂ３とする。このとき、本実施形態では、領域Ｒ１〜Ｒ３のそれぞれにおいて、値Ｂ１〜Ｂ３を用いて初期音圧Ｐ０（λ１）を補正する。
【００９２】
式（５）に式（６）を代入することで、異なる波長の吸収係数の比が、式（７）に示すに異なる波長の補正初期音圧（Ｐ１）の比として表される。
【数７】
2019193732.tif 000009
この式（７）を係数として式（１）に代入することで、酸素飽和度が算出される。
【００９３】
なお、上記のように、ｅｘｐ（−μｅｆｆ（λ１）・ｄ（ｒ））として深さにかかわらず一定の値を用いると、特性情報の取得対象となる被検体内部の位置が光照射面から深くなるほど（遠くなるほど）、光量算出の誤差が大きくなることが懸念される。このような誤差の拡大は、初期音圧の光量補正時の誤差拡大につながり、ひいては酸素飽和度の誤差が大きくなるおそれがある。しかし、被検体の浅い箇所（照射面から近い箇所）の酸素飽和度は、比較的安定して計測できる。かかる酸素飽和度の誤差に関しては、実施形態１と同様に、深さ（照射面からの距離）が１８ｍｍであれば、十分に実用的な精度で酸素飽和度を導出できる。さらに、実施形態１と同様に、許容できる算出誤差の値に基づいて、簡易的に光量補正を行う方法と、一般的だが計算量の多い補正する方法とを使い分けても良い。
【００９４】
さらに、実施形態１で図１１を用いて説明した場合と同様に、被検体００４複数回を計測する場合において、各回での計測条件を同じにすることで、深い箇所の濃度情報の再現
性が高くなる。この方法は、各回の間での変化を確認したい場合に特に好適である。
【００９５】
（処理フロー）
図８を参照して、信号処理部０１０が酸素飽和度分布を求める処理フローを説明する。図８のフローは、信号処理部０１０がデジタル信号を取得した状態からスタートする。すなわちフロー開始前に、光出射部００３が、波長ごとに、光源００１からの光を出射し、変換素子００６が、被検体００４から発生する光音響波を受信してアナログの受信信号に変換し、信号取得部０２２が、照射された光の波長ごとに順次、受信信号に対するＡＤ変換や増幅などを行う。
【００９６】
ステップＳ２０１では、初期音圧取得部０２４が、入力された受信信号を用いて、第１波長λ１における音圧（Ｐ０λ１（ｒ））分布情報と、第２波長λ２における音圧（Ｐ０λ２（ｒ））分布情報と、を取得する。
【００９７】
ステップＳ２０２では、光量分布取得部０２０が、任意の照射面における光照射領域での、２波長の照射光量分布を算出する。任意の照射面は、典型的には被検体表面であり、保持部００８を用いる場合は当該保持部００８の形状から求められる。その他、上述した各種の方法で表面形状情報を取得できる。また、被検体００４の形状にかかわらず、ユーザーが任意の位置に光量比取得面を設定しても構わない。
【００９８】
ステップＳ２０３では、初期音圧補正部０２５が、各波長において、Ｓ２０１で算出された音圧分布をＳ２０２で算出された光量分布で補正して、補正音圧分布情報を算出する。これにより、第１波長での補正音圧分布情報（Ｐ１λ１（ｒ））と、第２波長での補正音圧分布情報（Ｐ１λ２（ｒ））が得られる。
【００９９】
ステップＳ２０４では、表示制御部０１６が、複数波長分の補正音圧分布情報のうち少なくとも１波長分の補正音圧分布情報を基に画像処理を行い、補正音圧分布を示す画像を生成して表示部０１４に表示する。あるいは、補正音圧分布を示す画像に基づいて生成された画像を表示部０１４に表示する。なお、補正音圧分布（Ｐ１）ではなく補正前の音圧分布（Ｐ０）を用いた画像を表示しても構わない。
【０１００】
ステップＳ２０５では、濃度情報取得部０２６が、波長λ１における補正音圧（Ｐ１λ１（ｒ））分布情報と、波長λ２における補正音圧（Ｐ１λ２（ｒ））分布情報と、既知のオキシヘモグロビン・デオキシヘモグロビンのスペクトル情報を用いて、酸素飽和度の分布情報を生成する。
【０１０１】
ステップＳ２０６では、表示制御部０１６が、濃度情報取得部０２６により生成された酸素飽和度の分布情報を基に画像データを生成し、表示部０１４に表示する。これによりユーザーは、表示された音圧分布を示す画像あるいは音圧分布を示す画像に基づいて生成された画像から、酸素飽和度分布の画像を取得できる。
【０１０２】
以上説明したように、本実施形態では、酸素飽和度等の濃度に関する情報を求める際に、被検体表面に照射される異なる波長の光量を用いて補正された補正音圧分布が利用される。したがって、特許文献１の場合のようにユーザーの技量によって値が変化することがない。また、光量分布を算出するための散乱や吸収を考慮した複雑な演算を行う必要がない。その結果、簡易的な計算で安定した酸素飽和度の分布を取得できる。
【０１０３】
（その他の実施例）
本発明は、上述の実施例の１以上の機能を実現するプログラムを、ネットワーク又は記憶媒体を介してシステム又は装置に供給し、そのシステム又は装置のコンピュータにおけ
る１つ以上のプロセッサがプログラムを読出し実行する処理でも実現可能である。また、１以上の機能を実現する回路（例えば、ＡＳＩＣ）によっても実現可能である。
【符号の説明】
【０１０４】
００３：光出射部、００６：変換素子、００７：光量比分布取得部、０１０：信号処理部、０１１：探触子
【図１】
2019193732.tif 000010
【図２】
2019193732.tif 000011
【図３】
2019193732.tif 000012
【図４】
2019193732.tif 000013
【図５】
2019193732.tif 000014
【図６】
2019193732.tif 000015
【図７】
2019193732.tif 000016
【図８】
2019193732.tif 000017
【図９】
2019193732.tif 000018
【図１０】
2019193732.tif 000019
【図１１】
2019193732.tif 000020
【図１２】
2019193732.tif 000021