(57)【特許請求の範囲】
【請求項１】
光源と、
前記光源からの光が被検体に照射されることにより発生した光音響波を検出して複数の検出信号を出力する音響波検出器と、
前記複数の検出信号に基づいて関心領域における光学特性値を取得する信号処理手段と、を有し、
前記信号処理手段は、
前記複数の検出信号の一部の検出信号を選択し、
前記一部の検出信号と、前記一部の検出信号に対応する前記関心領域における光量値と、を用いて前記関心領域における前記光学特性値を取得する
ことを特徴とする被検体情報取得装置。
【請求項２】
前記音響波検出器は、光音響波を互いに異なる複数の位置で検出して前記複数の検出信号を出力する
ことを特徴とする請求項１に記載の被検体情報取得装置。
【請求項３】
前記光源は、光を異なる複数のタイミングで出射し、
前記音響波検出器は、前記複数のタイミングの光照射によりそれぞれ発生した光音響波を検出して前記複数の検出信号を出力する
ことを特徴とする請求項１または２に記載の被検体情報取得装置。
【請求項４】
前記音響波検出器を前記被検体に対して相対的に移動させる検出器走査機構を更に有し、
前記検出器走査機構は、前記複数のタイミングで前記音響波検出器の位置がそれぞれ異なるように前記音響波検出器を移動させる
ことを特徴とする請求項３に記載の被検体情報取得装置。
【請求項５】
前記光源からの光を前記被検体に導く光学系と、
前記光学系を前記被検体に対して相対的に移動させる光学系走査機構と、を更に有し、
前記光学系走査機構は、前記複数のタイミングで前記光学系の位置がそれぞれ異なるように前記光学系を移動させることを特徴とする請求項３または４に記載の被検体情報取得装置。
【請求項６】
前記音響波検出器は複数の音響波検出素子を備える
ことを特徴とする請求項１から５のいずれか１項に記載の被検体情報取得装置。
【請求項７】
前記信号処理手段は、
前記一部の検出信号に対応する前記関心領域における光量値を積算して前記関心領域における積算光量値を取得し、
前記一部の検出信号と、前記関心領域における積算光量値とを用いて、前記関心領域における前記光学特性値を取得する
ことを特徴とする請求項１から６のいずれか１項に記載の被検体情報取得装置。
【請求項８】
前記信号処理手段は、
前記一部の検出信号を用いて前記関心領域における初期音圧を取得し、
前記関心領域における前記初期音圧および前記積算光量値を用いて前記関心領域における前記光学特性値を取得する
ことを特徴とする請求項７に記載の被検体情報取得装置。
【請求項９】
前記信号取得手段は、前記音響波検出器の感度分布に基づいて前記複数の検出信号から前記一部の検出信号を選択する
ことを特徴とする請求項１から８のいずれか１項に記載の被検体情報取得装置。
【請求項１０】
前記信号取得手段は、
前記複数のタイミングのそれぞれにおいて、前記関心領域における前記音響波検出器の感度が所定の値以上のときに前記音響波検出器が検出した検出信号を選択し、選択された検出信号を前記一部の検出信号とすることを特徴とする請求項３から５のいずれか１項に記載の被検体情報取得装置。
【請求項１１】
前記信号取得手段は、
前記複数の音響波検出素子のそれぞれにおいて、前記関心領域における前記音響波検出素子の感度が所定の値以上のときに前記音響波検出素子が検出した検出信号を選択し、選択された検出信号を前記一部の検出信号とする
ことを特徴とする請求項６に記載の被検体情報取得装置。
【請求項１２】
前記音響波検出器の感度は、前記音響波検出器の光音響波から検出信号への変換効率、及び前記関心領域から前記音響波検出器までの光音響波の減衰のうち少なくとも１つに基づいて決定される
ことを特徴とする請求項９から１１のいずれか１項に記載の被検体情報取得装置。
【請求項１３】
前記信号処理手段は、前記関心領域における前記光学特性値として前記関心領域における吸収係数を取得する
ことを特徴とする請求項１から１２のいずれか１項に記載の被検体情報取得装置。
【請求項１４】
前記信号処理手段は、
前記一部の検出信号を重み係数に基づいて補正して補正信号を取得し、
前記一部の検出信号に対応する前記関心領域における光量値を前記重み係数に基づいて補正して、前記一部の検出信号に対応する前記関心領域における補正光量値を取得し、
前記補正信号と前記関心領域における前記補正光量値とを用いて前記関心領域における前記光学特性値を取得する
ことを特徴とする請求項１から１３のいずれか１項に記載の被検体情報取得装置。
【請求項１５】
光を異なる複数のタイミングで出射する光源と、
前記複数のタイミングで光が被検体に照射されることにより発生した光音響波を検出して複数の検出信号を出力する音響波検出器と、
前記複数の検出信号に基づいて関心領域における光学特性値を取得する信号処理手段と、を有し、
前記信号処理手段は、
前記複数の検出信号の少なくとも一部の検出信号を選択し、
前記少なくとも一部の検出信号と、前記少なくとも一部の検出信号に対応する前記関心領域における光量値と、を用いて前記関心領域における前記光学特性値を取得する
ことを特徴とする被検体情報取得装置。
【請求項１６】
前記音響波検出器を前記被検体に対して相対的に移動させる検出器走査機構を更に有し、
前記検出器走査機構は、前記複数のタイミングで前記音響波検出器の位置がそれぞれ異なるように前記音響波検出器を移動させる
ことを特徴とする請求項１５に記載の被検体情報取得装置。
【請求項１７】
前記光源からの光を前記被検体に導く光学系と、
前記光学系を前記被検体に対して相対的に移動させる光学系走査機構と、を更に有し、
前記光学系走査機構は、前記複数のタイミングで前記光学系の位置がそれぞれ異なるように前記光学系を移動させることを特徴とする請求項１５または１６に記載の被検体情報取得装置。
【請求項１８】
光が被検体に照射されることにより発生した光音響波を検出して得られた複数の検出信号に基づいて関心領域における光学特性値を取得する信号処理方法であって、
前記複数の検出信号の一部の検出信号を選択する工程と、
前記一部の検出信号と、前記一部の検出信号に対応する前記関心領域における光量値と、を用いて前記関心領域における前記光学特性値を取得する工程と、
を有する
ことを特徴とする信号処理方法。
【請求項１９】
請求項１８に記載の信号処理方法をコンピュータに実行させるプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被検体に光を照射することにより発生した光音響波を検出することにより、被検体情報を取得する被検体情報取得装置及び被検体情報取得方法に関する。
【背景技術】
【０００２】
レーザーなどの光源から被検体に照射した光を被検体内に伝播させ、被検体内の情報を得る光イメージング装置の研究が医療分野を中心に積極的に進められている。このような光イメージング技術の一つとして、光音響イメージング（ＰＡＩ：Ｐｈｏｔｏａｃｏｕｓｔｉｃ  Ｉｍａｇｉｎｇ）がある。光音響イメージングとは、光源から発生したパルス光を被検体（生体）に照射し、被検体内を伝播・拡散した光が被検体内で吸収されて発生する光音響波を検出し、検出した音響波を解析処理することで、被検体内の光学特性に関連した情報を可視化する技術である。これにより、被検体内の光学特性値分布、特に、吸収係数分布、酸素飽和度分布などを得ることができる。【０００３】
光音響イメージングでは、被検体内における関心領域から発生する光音響波の初期音圧Ｐ０は、次式で表すことができる。【０００４】
【数１】
0005871958.tif 000002
【０００５】
ここで、Γはグルナイゼン係数であり、体積膨張係数βと音速ｃの二乗の積を定圧比熱ＣＰで割ったものである。Γは被検体が決まれば、ほぼ一定の値をとることが知られている。そして、μａは関心領域の吸収係数、Φは関心領域における積算光量値である。【０００６】
特許文献１には、被検体中を伝搬してきた光音響波の音圧Ｐの時間変化を音響波検出器で検出し、その検出結果から被検体内の初期音圧分布を算出する技術が記載されている。特許文献１によると、算出された初期音圧をグルナイゼン係数Γで除することにより、μａとΦの積、つまり光エネルギー吸収密度を得ることができる。そして、式（１）で示されるように、初期音圧Ｐ０から吸収係数μａを得るためには、光エネルギー吸収密度を光量Φで除することが必要である。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２０１０−８８６２７号公報
【特許文献２】特開２００６−５１３５５号公報
【発明の概要】
【発明が解決しようとする課題】
【０００８】
しかしながら、特許文献１に記載の光音響イメージングにおいては、さらに光学特性値を精度良く取得することが望まれていた。【０００９】
そこで、本発明は、光音響イメージングにおいて、光学特性値をより精度良く取得することのできる被検体情報取得装置及び被検体情報取得方法を提供することを目的とする。
【課題を解決するための手段】
【００１０】
本発明に係る被検体情報取得装置は、光源と、光源からの光が被検体に照射されることにより発生した光音響波を検出して複数の検出信号を出力する音響波検出器と、複数の検出信号に基づいて関心領域における光学特性値を取得する信号処理手段と、を有し、信号処理手段は、複数の検出信号の一部の検出信号を選択し、一部の検出信号と、一部の検出信号に対応する関心領域における光量値と、を用いて関心領域における前記光学特性値を取得する。
【発明の効果】
【００１１】
本発明によれば、光学特性値をより精度良く取得することのできる被検体情報取得装置及び被検体情報取得方法を提供することができる。
【図面の簡単な説明】
【００１２】
【図１】第１の実施形態に係る被検体情報取得装置の模式図である。
【図２】第１の実施形態に係る被検体情報取得方法のフローチャート図である。
【図３Ａ】第２の実施形態に係る被検体情報取得装置の模式図である。
【図３Ｂ】第２の実施形態に係る被検体情報取得装置の模式図である。
【図３Ｃ】第２の実施形態に係る被検体情報取得装置の模式図である。
【図４Ａ】第４の実施形態に係る被検体情報取得装置の模式図である。
【図４Ｂ】第４の実施形態に係る被検体情報取得装置の模式図である。
【図４Ｃ】第４の実施形態に係る被検体情報取得装置の模式図である。
【図５Ａ】第４の実施形態に係る別の被検体情報取得装置の模式図である。
【図５Ｂ】第４の実施形態に係る別の被検体情報取得装置の模式図である。
【図５Ｃ】第４の実施形態に係る別の被検体情報取得装置の模式図である。
【発明を実施するための形態】
【００１３】
光音響イメージングにおいて、光音響波を検出して取得した検出信号は、バックグランドノイズを含んでいる。そのため、光音響イメージングにおいては、バックグランドノイズを含むＳ／Ｎが低い検出信号を用いずに、関心領域における初期音圧を取得することが望ましい。例えば、特許文献２には、超音波イメージングではあるが、関心領域と音響波検出素子とのなす角度が所定の値以下である場合（音響波検出素子に対応する所定の感度領域に関心領域が含まれていない場合）、音響波検出素子が関心領域からの音響波の受信を阻止することが記載されている。このような方法を行うことにより、Ｓ／Ｎが低い検出信号を用いずに超音波画像を取得している。
【００１４】
そこで、本発明者は、特許文献２に記載の技術を光音響イメージングに適用した。具体的には、シミュレーションにより、所定の感度領域に関心領域が含まれていない音響波検出素子が取得した検出信号を用いずに再構成を行い、関心領域における初期音圧を取得した。このようにして得られた初期音圧は、Ｓ／Ｎが低い検出信号を用いずに再構成された初期音圧であるため、ノイズによる誤差が小さい。そして、本発明者は、この初期音圧を用いて、特許文献１に記載の方法で、関心領域における吸収係数を求めた。しかしながら、上記方法で算出された吸収係数の値は、シミュレーションで設定した吸収係数の値とは異なっていた。
【００１５】
そこで、上記課題に鑑み、本発明者が鋭意検討した結果、初期音圧を取得する際には、音響波検出素子の感度に基づいて用いる検出信号を選択していたのに対し、積算光量値を取得する際には、音響波検出素子の感度を考慮していなかったことが原因であることがわかった。
【００１６】
そこで、本発明者は、吸収係数の取得の際に、音響波検出素子の感度に基づいて用いる検出信号を選択することに加え、音響波検出素子の感度に基づいて積算光量値を取得することにより、精度良く光学特性値としての吸収係数を取得できることを見出した。
【００１７】
以下に、シミュレーションによる本発明に係る実施形態を説明する。
【００１８】
（第１の実施形態）
図１は、本実施形態に係る被検体情報取得装置の模式図である。光源１０から発したパルス光は、光学系１１に導かれ、照射光１２として被検体３０に照射される。被検体３０内の光吸収体３１から発生した光音響波３２は、音響波検出素子ｅ１，ｅ２，ｅ３を備えた音響波検出器２０で検出される。そして、音響波検出器２０が取得した複数の検出信号は、信号収集器４７で増幅、デジタル変換され、信号処理装置４０のメモリに格納される。そして、信号処理手段としての信号処理装置４０に備えられた初期音圧取得部としての初期音圧取得モジュール４２は、複数の検出信号を用いて画像再構成することにより、被検体３０内の関心領域３３における初期音圧を取得する。また、信号処理装置４０に備えられた光量値取得部としての光量値取得モジュール４３は、関心領域３３における積算光量値を取得する。そして、信号処理装置４０に備えられた光学特性値取得部としての光学特性値取得モジュール４４は、関心領域３３における初期音圧及び光量値を用いて、関心領域３３における光学特性値を取得する。そして、取得した光学特性値を表示手段としての表示装置５０に表示させる。
【００１９】
ここで、関心領域とは、初期音圧取得モジュール４２により再構成される領域の最小単位であるボクセルのことを指す。なお、初期音圧取得モジュール４２は、被検体３０の全領域にわたって関心領域を設定することにより、被検体全体の初期音圧分布を取得することができる。また、同様に、光量値取得モジュール４３および光学特性値取得モジュール４４は、被検体の全領域にわたって関心領域を設定することにより、被検体全体の積算光量値分布および吸収係数分布を取得することができる。
【００２０】
ここで、図１に示す音響波検出素子ｅ１，ｅ２，ｅ３のそれぞれが取得した関心領域３３に対応する検出信号をＰｄ１（ｒＴ）、Ｐｄ２（ｒＴ）、Ｐｄ３（ｒＴ）とする。また、音響波検出素子の正面から入ってきた光音響波に対して、音響波検出素子の正面からθだけ角度を持って入ってきた光音響波の検出信号への変換効率をＡ（θ）とする。そして、関心領域３３に対してそれぞれの音響波検出素子のなす角度をθ１、θ２、θ３とすると、それぞれの音響波検出素子の指向性による変換効率は、Ａ（θ１），Ａ（θ２），Ａ（θ３）と表現できる。また、検出信号Ｐｄ１（ｒＴ）、Ｐｄ２（ｒＴ）、Ｐｄ３（ｒＴ）のそれぞれに対応する関心領域３３における光量値Φ１（ｒＴ）、Φ２（ｒＴ）、Φ３（ｒＴ）とする。ここで、本実施形態では、関心領域３３を光吸収体３１の位置ｒＴに設定している。
【００２１】
ここで、音響波検出素子から関心領域３３までの距離をｒとし、被検体内の光音響波の伝搬速度をｃとし、被検体３０に照射光１２を照射した時間をｔ＝０とする。この場合、関心領域に対応する検出信号とは、ｔ＝ｒ／ｃの時間に音響波検出素子が取得した検出信号のことを指す。また、関心領域に対応する検出信号に対応する、関心領域３３における光量値とは、ｔ＝０の時間に照射された照射光１２の関心領域３３における光量値のことを指す。
【００２２】
（一部の検出信号を用いないシミュレーション例）
以下に、図１を用いて、音響波検出素子の感度に基づき、一部の検出信号を用いずに取得した初期音圧より、吸収係数を取得するシミュレーションの例を説明する。本シミュレーションでは、光吸収体３１の吸収係数を、μａ＝０．０８８／ｍｍと設定した。
【００２３】
まず、初期音圧取得モジュール４２は、式（２）に示すように検出信号Ｐｄ１（ｒＴ）、Ｐｄ２（ｒＴ）、Ｐｄ３（ｒＴ）と、変換効率Ａ（θ１），Ａ（θ２），Ａ（θ３）を用いて、関心領域３３における初期音圧Ｐ０（ｒＴ）を取得する。
【００２４】
【数２】
0005871958.tif 000003
【００２５】
ここで、シミュレーションにより得られた検出信号、および、シミュレーションで設定した変換効率は、以下のとおりである。
Ｐｄ１（ｒＴ）＝１３２Ｐａ
Ｐｄ２（ｒＴ）＝２３１Ｐａ
Ｐｄ３（ｒＴ）＝１９８Ｐａ
Ａ（θ１）＝０．４
Ａ（θ２）＝０．７
Ａ（θ３）＝０．６
【００２６】
そして、これらのパラメータを用いて、式（２）より算出した初期音圧は、Ｐ０（ｒＴ）＝９９０となる。
【００２７】
また、図１においては、音響波検出素子の変換効率が所定の値より大きい領域（所定の感度領域）を、点線の三角形の領域で示す。ここでは、変換効率Ａ（θ）＝０．５を所定の値として設定している。
【００２８】
ここで、本実施形態においては、音響波検出素子ｅ１に対応する三角形の領域（所定の感度領域）に関心領域３３が含まれていない。そのため、初期音圧取得モジュール４２は、音響波検出素子ｅ１が取得した関心領域３３に対応する検出信号Ｐｄ１（ｒＴ）を用いずに、式（３）で示される関心領域３３における初期音圧Ｐ０’（ｒＴ）を取得する。
【００２９】
【数３】
0005871958.tif 000004
【００３０】
そして、前述したパラメータを用いて、式（３）から算出される関心領域３３における初期音圧は、Ｐ０’（ｒＴ）＝６６０となった。
【００３１】
次に、光量値取得モジュール４３は、被検体の背景光学係数などから光伝搬モンテカルロ法・輸送方程式・光拡散方程式等を用いて、被検体内の積算光量値を取得する。
【００３２】
例えば、光量値取得モジュール４３は、検出信号Ｐｄ１（ｒＴ）、Ｐｄ２（ｒＴ）、Ｐｄ３（ｒＴ）のそれぞれに対応する関心領域３３における光量値Φ１（ｒＴ）、Φ２（ｒＴ）、Φ３（ｒＴ）を算出する。
【００３３】
そして、光量値取得モジュール４３は、これらを用いて式（４）で示される関心領域３３における積算光量値Φ（ｒＴ）を取得する。
【００３４】
【数４】
0005871958.tif 000005
【００３５】
ここで、シミュレーションで得られた関心領域における光量値は、以下のとおりであった。
Φ１（ｒＴ）＝３７５０ｍＪ／ｍ２
Φ２（ｒＴ）＝３７５０ｍＪ／ｍ２
Φ３（ｒＴ）＝３７５０ｍＪ／ｍ２
【００３６】
そして、これらのパラメータを用いて、式（４）より関心領域における積算光量値を算出すると、Φ（ｒＴ）＝１１２５０ｍＪ／ｍ２となる。
【００３７】
次に、光学特性値取得モジュール４４は、式（３）に示す関心領域３３における初期音圧Ｐ０’（ｒＴ）と、式（４）に示す関心領域３３における積算光量値Φ（ｒＴ）とを用いて、式（５）に示す関心領域３３における吸収係数μａ（ｒＴ）を取得する。
【００３８】
ここで、グリュナイゼン係数Γ＝１としている。
【００３９】
【数５】
0005871958.tif 000006
【００４０】
ここで、前述したパラメータを用いて式（５）より算出される、光吸収体の位置ｒＴに設定された関心領域３３における吸収係数は、μａ＝０．０５９／ｍｍとなる。一方、シミュレーションで設定された光吸収体３１の吸収係数は、μａ＝０．０８８／ｍｍである。これより、式（５）より求めた吸収係数は設定値より小さくなっていることがわかる。すなわち、上記方法により取得した初期音圧を用いて吸収係数を取得する場合には、積算光量値の取得にさらなる工夫が必要である。
【００４１】
（一部の検出信号および一部の光量値を用いないシミュレーション例）
そこで、以下に、本発明者が見出した本実施形態に係る被検体情報取得方法を、図２のフローチャートを用いて説明する。以下の番号は、図２に示す処理番号と一致する。
【００４２】
（Ｓ１００：音響波検出素子の感度分布に基づいて、所定の感度領域を設定する工程）
この工程では、複数の音響波検出素子の感度分布に基づいて、複数の音響波検出素子のそれぞれに対応する所定の感度領域を設定する。それぞれの音響波検出素子に対応する所定の感度領域のテーブルを、信号処理装置４０のメモリに格納する。
【００４３】
ここでは、信号処理装置４０に備えられた設定部としての設定モジュール４１が、音響波検出素子の感度が所定の値より大きい領域を、所定の感度領域として設定してもよい。なお、所定の値は、システムノイズに基づき、設定モジュール４１により自動で設定されてもよい。また、所定の値は、音響波検出素子の感度をヒストグラムとして表示装置５０に表示させ、作業者がヒストグラムに基づいて所定の値を選択してもよい。このとき、所定の値は、システムノイズを考慮して選択することが好ましい。
【００４４】
ここで、音響波検出素子の感度とは、例えば、音響波検出素子の変換効率や、関心領域から音響波検出素子までの光音響波の拡散や散乱による減衰を示す減衰率などによって決定される。なお、変換効率は、光音響波が音響波検出素子に入射する角度などにより決定される。また、減衰率は、関心領域と音響波検出素子との距離などにより決定される。
【００４５】
例えば、上記で示したシミュレーションの例の場合、設定モジュール４１は、変換効率Ａ（θ）＝０．５を所定の値として設定した。そして、音響波検出素子ｅ１，ｅ２，ｅ３のそれぞれについて、変換効率Ａ（θ）が０．５より大きい領域を三角形の領域で示した。その結果、音響波検出素子ｅ１の変換効率Ａ（θ）が０．５より大きい領域（所定の感度領域）に関心領域３３が含まれていなかった。
【００４６】
また、音響波検出素子の感度分布の画像から選択された任意の領域に基づき、所定の感度領域を設定することもできる。
【００４７】
例えば、まず、表示装置５０に、信号処理装置４０のメモリに格納しておいた音響波検出素子の感度分布の画像データを表示させる。そして、作業者が、表示された感度分布の画像からＰＣの入力デバイスを用いて任意の領域を選択する。そして、設定モジュール４１は、選択された任意の領域を、所定の感度領域として設定することができる。このとき、例えば、感度分布の画像を表示しながら、マウスによる認識や、タッチパネル上のセンサによる認識方法で始点から終点までを結んで任意の領域を選択することができる。
【００４８】
なお、設定モジュール４１は、選択された任意の領域の感度分布に基づき、所定の感度領域として設定してもよい。例えば、最も小さい感度を基準として所定の感度領域を設定したりすることができる。
【００４９】
また、それぞれの音響波検出素子について個別に所定の感度領域を設定してもよいし、あるいは、１つの音響波検出素子について所定の感度領域を設定し、その所定の感度領域と同様の感度領域を、他の音響波検出素子に対して設定してもよい。
【００５０】
（Ｓ２００：所定の感度領域に関心領域が含まれない音響波検出素子が取得した検出信号を用いずに、関心領域における初期音圧を取得する工程）
この工程では、Ｓ１００で設定した所定の感度領域に、関心領域が含まれない音響波検出素子について、その音響波検出素子が取得した関心領域に対応する検出信号を用いずに、関心領域における初期音圧を取得する。そして、この初期音圧のデータを信号処理装置４０のメモリに格納する。
【００５１】
例えば、上記で示したシミュレーションの例の場合、音響波検出素子ｅ１に対応する所定の感度領域に関心領域３３が含まれていなかった。そのため、初期音圧取得モジュール４２は、検出信号Ｐｄ１（ｒＴ）、Ｐｄ２（ｒＴ）、Ｐｄ３（ｒＴ）のうち、音響波検出素子ｅ１が取得した関心領域に対応する検出信号Ｐｄ１（ｒＴ）を用いずに、画像再構成を行うことにより、式（４）で示された初期音圧Ｐ０’（ｒＴ）を取得する。
【００５２】
このとき、初期音圧取得モジュール４２が行う画像再構成アルゴリズムとしては、例えば、トモグラフィー技術で通常に用いられるタイムドメインあるいはフーリエドメインでの逆投影などがある。
【００５３】
なお、本発明においては、関心領域の一部にでも所定の感度領域が含まれていれば、所定の感度領域に関心領域が含まれているとすることができる。
【００５４】
また、本発明において検出信号を用いないとは、初期音圧を取得する際に、検出信号を全く用いないことも、実質的に検出信号を用いていないことも含む概念である。
【００５５】
（Ｓ３００：初期音圧を取得する際に用いない検出信号に対応する光量値を用いずに、関心領域における積算光量値を取得する工程）
この工程では、Ｓ２００で用いなかった検出信号に対応する、関心領域における光量値を用いずに、関心領域における積算光量値を取得する。そして、この積算光量値のデータを、信号処理装置４０のメモリに格納する。
【００５６】
例えば、光量値取得モジュール４３は、光量値Φ１（ｒＴ）、Φ２（ｒＴ）、Φ３（ｒＴ）のうち、初期音圧取得モジュール４２が用いなかった検出信号Ｐｄ１（ｒＴ）に対応する関心領域における光量値Φ１（ｒＴ）を用いずに、式（６）で示される関心領域における積算光量値Φ’（ｒＴ）を取得する。
【００５７】
【数６】
0005871958.tif 000007
【００５８】
すなわち、光量値取得モジュール４３は、初期音圧取得モジュール４２が初期音圧を算出する際に用いた検出信号に対応する、関心領域における光量値を用いて、関心領域における積算光量値を取得している。
【００５９】
ここで、前述したパラメータを用いて、式（６）から求めた関心領域３３における積算光量値は、Φ’（ｒＴ）＝７５００ｍＪ／ｍ２となる。
【００６０】
なお、本発明において光量値を用いないとは、積算光量値を取得する際に、光量値を全く用いないことも、実質的に光量値を用いていないことも含む概念である。
【００６１】
また、本実施形態では、照射光１２の照射条件が一定であるため、関心領域３３に照射される光量は一定である。そのため、複数の音響波検出素子が取得した複数の検出信号に対応する光量値も一定となる。このような場合、光量値取得モジュール４３は、初期音圧取得モジュール４２が初期音圧を取得する際に用いた検出信号の数に、関心領域３３に到達した光量を掛けた値を、関心領域３３における積算光量値として取得してもよい。本発明において、このように取得した積算光量値も、光量値を用いずに取得した積算光量値、として取り扱う。
【００６２】
（Ｓ４００：関心領域における初期音圧および積算光量値を用いて、関心領域における光学特性値を取得する工程）
この工程では、Ｓ２００で取得した関心領域における初期音圧と、Ｓ３００で取得した関心領域における積算光量値とを用いて、関心領域における光学特性値としての吸収係数を取得する。
【００６３】
例えば、上記で示されたシミュレーションの例の場合、光学特性値取得モジュール４４は、式（３）に示す初期音圧Ｐ０’（ｒＴ）と、式（６）に示す積算光量値Φ’（ｒＴ）と、を式（１）に適用させる。そして、式（７）で示される関心領域３３における吸収係数μａ（ｒＴ）を取得する。ここで、グリュナイゼン係数Γ＝１としている。
【００６４】
【数７】
0005871958.tif 000008
【００６５】
例えば、前述したパラメータを用いて、式（７）より求められた関心領域３３における吸収係数は、μａ＝０．０８８／ｍｍである。一方、式（５）より求められた関心領域３３における吸収係数は、μａ＝０．０５９／ｍｍであった。そして、シミュレーションで設定した光吸収体３１の吸収係数は、μａ＝０．０８８／ｍｍである。すなわち、式（７）によれば、式（５）と比べ、より精度良く吸収係数を求めることができる。
【００６６】
以上のように、所定の感度領域に関心領域に含まれない音響波検出素子について、その音響波検出素子が取得した関心領域に対応する検出信号と、その検出信号に対応する関心領域における光量値と、を用いないことにより、ノイズによる誤差が少なく、定量性の高い吸収係数を取得することができる。
【００６７】
なお、以上の工程を複数の波長において行うことにより、波長毎における吸収係数を取得してもよい。そして、これらの吸収係数を用いて、光学特性値としての酸素飽和度などを取得してもよい。
【００６８】
また、以上の工程を含んだプログラムを、コンピュータとしての信号処理装置４０に実行させてもよい。
【００６９】
（第２の実施形態）
図３Ａ〜３Ｃは、本実施形態に係る被検体情報取得装置の模式図である。
【００７０】
本実施形態に係る被検体情報取得装置は、１つの音響波検出素子を備えた音響波検出器２０を有している。また、被検体３０と音響波検出器２０とを相対的に移動させるための検出器走査機構２１を有している。本実施形態において、検出器走査機構２１は、１つの音響波検出素子を備えた音響波検出器２０を紙面右方向に走査することにより、複数の位置で光音響波を検出可能としている。ここで、図３Ａ、３Ｂ、３Ｃに示すそれぞれの位置における音響波検出素子を、ｅ１，ｅ２，ｅ３とする。また、三角形で示した領域は、音響波検出素子に対応する所定の感度領域を示す。
【００７１】
本発明において、複数の音響波検出素子とは、音響波検出素子が複数の位置で光音響波を検出可能であることを指す。すなわち、本実施形態のように、音響波検出器２０を走査することにより、複数の位置で光音響波を検出可能とした音響波検出素子についても、複数の音響波検出素子とする。
【００７２】
また、本実施形態に係る被検体情報取得装置は、照射光１２を走査するために、光学系１１を走査する光学走査機構１３が設けられている。そして、本実施形態では、音響波検出器２０と照射光１２とが同期して走査されている。このように、音響波検出器２０と照射光１２とを同期して走査することによって、音響波検出素子に対応する所定の感度領域（三角形の領域）に照射光１２が常に照射されるため、Ｓ／Ｎの高い検出信号を常に取得することが可能となる。
【００７３】
本実施形態に係る被検体情報取得装置においても、第１の実施形態と同様に、音響波検出素子ｅ１に対応する所定の感度領域に関心領域３３が含まれていない。そのため、初期音圧取得モジュール４２は、音響波検出素子ｅ１が取得した関心領域３３に対応する検出信号を用いずに、関心領域３３における初期音圧を取得する。そして、光量値取得モジュール４３は、初期音圧を取得する際に用いない検出信号に対応する関心領域３３における光量値を用いずに、関心領域３３における積算光量値を取得する。そして、光学特性値取得モジュール４４は、この初期音圧と積算光量値とを用いて、式（７）で示される関心領域３３における吸収係数を取得する。このように吸収係数を取得することにより、本実施形態においても、高精度に吸収係数を取得することができる。
【００７４】
（第３の実施形態）
第１の実施形態および第２の実施形態では、所定の感度領域に関心領域が含まれていない音響波検出素子について、その音響波検出素子が取得した関心領域に対応する検出信号、および、その検出信号に対応する関心領域における光量値を用いずに吸収係数を取得していた。一方、本実施形態では、上記検出信号と、この検出信号に対応する関心領域における光量値とを低減して、吸収係数を取得する。
【００７５】
以下に、本実施形態に係る被検体情報取得方法を、図１に示す被検体情報取得装置を用いて説明する。
【００７６】
本実施形態において、初期音圧取得モジュール４２は、所定の感度領域に関心領域３３が含まれない音響波検出素子ｅ１が取得した、関心領域３３に対応する検出信号に、第１の低減係数を掛ける。そして、初期音圧取得モジュール４２は、第１の低減係数がかけられた検出信号も用いて、関心領域３３における初期音圧を取得する。
【００７７】
このように、所定の感度領域に関心領域が含まれない音響波検出素子について、この音響波検出素子が取得した関心領域に対応する検出信号に、第１の低減係数を掛けることにより、Ｓ／Ｎが低い検出信号を低減して、初期音圧を取得することができる。そのため、ノイズによる誤差の少ない初期音圧を取得することができる。
【００７８】
次に、光量値取得モジュール４３は、第１の低減係数をかけられた検出信号に対応する関心領域３３における光量値に、第２の低減係数を掛ける。そして、光量値取得モジュール４３は、第２の低減係数がかけられた光量値も用いて、関心領域３３における積算光量値を取得する。
【００７９】
そして、光学特性値取得モジュール４４は、初期音圧取得モジュール４２が取得した初期音圧と、光量値取得モジュール４３が取得した積算光量値とを用いて、関心領域３３における吸収係数を取得する。
【００８０】
このように、検出信号に第１の低減係数を掛けることに加え、この検出信号に対応する光量値にも第２の低減係数を掛けることにより、吸収係数を高精度に取得することができる。
【００８１】
なお、第１の低減係数および第２の低減係数は、１より小さい値である。また、関心領域によって別の低減係数を設定してもよい。また、第１の低減係数および第２の低減係数は同一の値であることが好ましい。ここで、同一の値とは、全く同一の値も、吸収係数を取得する際に実質的に同一となるような値も含む概念である。
【００８２】
（第４の実施形態）
本発明は、図４Ａ〜４Ｃに示される被検体情報取得装置や、図５Ａ〜５Ｃに示される被検体情報取得装置にも適用可能である。図４Ａ〜４Ｃに示す被検体情報取得装置は、検出器走査機構２１が音響波検出器２０を被検体３０の周囲を回転走査することにより、複数の位置で光音響波を検出可能としている。また、被検体３０と音響波検出器２０との間の音響インピーダンスマッチングを図るために、被検体３０は水槽８１に満たされた水８０に浸かっている。また、被検体３０を走査する被検体走査機構３４を有している。このような構成とすることで、保持板等で形状を規定出来ない部位でも測定可能となる。また、被検体に対して多くの方向に検出素子を設置可能となるため、情報量の多いデータの取得が可能となる。
【００８３】
図４Ａ〜４Ｃに示す被検体情報取得装置では、図４Ａの状態から、被検体走査機構３４が被検体を紙面下方向に走査することにより、図４Ｂの状態となる。そして、図４Ｂの状態から、検出器走査機構２１が音響波検出器２０を走査することにより、図４Ｃの状態となる。ここで、図４Ａ、４Ｂ、４Ｃのそれぞれの状態における音響波検出素子を、ｅ１，ｅ２，ｅ３とする。また、点線で示した三角形の領域は、音響波検出素子に対応する所定の感度領域を示す。
【００８４】
また、図５Ａ〜５Ｃで示された被検体情報取得装置は、音響波検出器２０と光学系１１とが１つのハウジング７０に収めて設けられている。また、このハウジング７０はハンドヘルド機構７１を備え、作業者がハンドヘルド機構７１を把持して、ハウジング７０を走査することを可能にしている。このようにハウジング７０を走査することによって、音響波検出素子が複数の位置で光音響波を検出可能としている。図５Ａ〜５Ｃにおいては、作業者がハンドヘルド機構７１を把持してハウジング７０を紙面右方向に走査することにより、音響波検出素子が光音響波を検出している。ここで、図５Ａ、５Ｂ、５Ｃのそれぞれの状態における音響波検出素子を、ｅ１，ｅ２，ｅ３とする。また、点線で示した三角形の領域は、音響波検出素子に対応する所定の感度領域を示す。
【００８５】
ただし、本実施形態は、他の実施形態と異なり、音響波検出器２０を機械的に走査するのではなく、ハンドヘルド機構７１を作業者が把持してハウジング７０を自由に走査している。そのため、光音響波３２を検出した時の音響波検出器２０と関心領域３３との位置関係を把握することができない。しかし、音響波検出器２０が取得した検出信号から関心領域に対応する検出信号を抽出するためには、音響波検出器２０と関心領域３３との位置関係を把握する必要がある。そこで、本実施形態においては、ハウジング７０の位置、すなわちハウジング７０に収められた音響波検出器２０及び光学系１１の位置を検出するための位置検出器７２をハウジング７０が備えていることが好ましい。
【００８６】
図４Ａ〜４Ｃおよび図５Ａ〜５Ｃに示す被検体情報取得装置において、音響波検出素子ｅ１に対応する所定の感度領域に関心領域３３が含まれていない。そのため、信号処理装置４０は、第１および第２の実施形態で説明した被検体情報取得方法や、第３の実施形態で説明した被検体情報取得方法を用いて、関心領域３３の吸収係数を取得することができる。このように吸収係数を取得することにより、本実施形態においても、高精度に吸収係数を取得することができる。
【００８７】
以下、主要な構成について説明する。
【００８８】
（光源１０）
光源１０は、５ナノ秒乃至５０ナノ秒のパルス光を発生可能な光源を備えている。光源としては大きな出力が得られるレーザーが好ましいが、レーザーのかわりに発光ダイオードなどを用いることも可能である。レーザーとしては、固体レーザー、ガスレーザー、色素レーザー、半導体レーザーなど様々なレーザーを使用することができる。理想的には、出力が強く連続的に波長を変えられる、Ｎｄ：ＹＡＧ励起のＴｉ：Ｓａレーザーや、アレキサンドライトレーザーがよい。異なる波長の単波長レーザーを複数で保有していてもよい。
【００８９】
（光学系１１）
光源１０から出射されたパルス光は、典型的にはレンズやミラーなどの光学部品により、所望の光分布形状に加工されながら被検体に導かれるが、光ファイバなどの光導波路などを用いて伝搬させることも可能である。光学系１１は、例えば、光を反射するミラーや、光を集光したり拡大したり形状を変化させるレンズ、光を拡散させる拡散板などである。このような光学部品は、光源から発せられたパルス光が被検体に所望の形状で照射されれば、どのようなものを用いてもかまわない。なお、光はレンズで集光させるより、ある程度の面積に広げる方が被検体への安全性ならびに診断領域を広げられるという観点で好ましい。なお、照射光を走査するために、光学系１１に光学走査機構を設けてもよい。
【００９０】
（音響波検出器２０）
光により被検体表面及び被検体内部で発生する光音響波を検出する検出器である音響波検出器２０は、音響波を検知し、アナログ信号である電気信号に変換するものである。以後、単に探触子あるいはトランスデューサということもある。圧電現象を用いたトランスデューサ、光の共振を用いたトランスデューサ、容量の変化を用いたトランスデューサなど音響波信号を検知できるものであれば、どのような音響波検出器を用いてもよい。
【００９１】
また、音響波検出器２０は、複数の音響波検出素子を備えている。この複数の音響波検出素子を１次元、または２次元に並べてアレイ状に配置することにより、複数の位置で光音響波を検出可能としている。このよう多次元配列素子を用いることで、同時に複数の位置で音響波を検出することができ、検出時間を短縮できると共に、被検体の振動などの影響を低減できる。
【００９２】
なお、複数の位置で光音響波を検出可能とするために、音響波検出器２０が検出器走査機構２１によって機械的に走査可能に構成されていてもよい。また、作業者が把持して音響波検出器２０を自由に走査できるハンドヘルド機構を備えることもできる。
【００９３】
（信号収集器４７）
音響波検出器２０より得られた電気信号を増幅し、その電気信号をアナログ信号からデジタル信号に変換する信号収集器４７を有することが好ましい。信号収集器４７は、典型的には増幅器、Ａ／Ｄ変換器、ＦＰＧＡ（Ｆｉｅｌｄ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ）チップなどで構成される。音響波検出器から得られる検出信号が複数の場合は、同時に複数の信号を処理できることが望ましい。それにより、画像を形成するまでの時間を短縮できる。なお、本明細書において「検出信号」とは、音響波検出器２０から出力されるアナログ信号も、信号収集器４７によりＡＤ変換されたデジタル信号も含む概念である。
【００９４】
（信号処理装置４０）
信号処理装置４０は、画像再構成などを行うことにより被検体内部の光学特性値を取得する。信号処理装置４０には、典型的にはワークステーションなどが用いられ、画像再構成処理などがあらかじめプログラミングされたソフトウェアにより行われる。例えば、ワークステーションで使われるソフトウェアとしては、設定モジュール４１、初期音圧取得モジュール４２、光量値取得モジュール４３、光学特性値取得モジュール４４などがある。
【００９５】
なお、それぞれのモジュールを、別々のハードウェアとして設けてもよい。この場合、それぞれのハードウェアを総じて信号処理装置４０としてもよい。
【００９６】
また、場合によっては、信号収集器４７、信号処理装置４０は一体化される場合もある。この場合、ワークステーションで行うようなソフトウェア処理ではなく、ハードウェア処理により被検体の光学特性値を生成することもできる。
【００９７】
本発明は上記実施の形態に制限されるものではなく、本発明の精神及び範囲から離脱することなく、様々な変更及び変形が可能である。従って、本発明の範囲を公にするために以下の請求項を添付する。
【符号の説明】
【００９８】
２０  音響波検出器
３０  被検体
４０  信号処理装置
４１  設定モジュール
４２  初期音圧取得モジュール
４３  光量値取得モジュール
４４  光学特性値取得モジュール
【図１】
0005871958.tif 000009
【図２】
0005871958.tif 000010
【図３Ａ】
0005871958.tif 000011
【図３Ｂ】
0005871958.tif 000012
【図３Ｃ】
0005871958.tif 000013
【図４Ａ】
0005871958.tif 000014
【図４Ｂ】
0005871958.tif 000015
【図４Ｃ】
0005871958.tif 000016
【図５Ａ】
0005871958.tif 000017
【図５Ｂ】
0005871958.tif 000018
【図５Ｃ】
0005871958.tif 000019