JP6686066B2 - 光音響装置 - Google Patents

Description

本発明は、光を被検体に照射し発生する光音響波を受信する光音響装置に関する。

レーザーなどの光源から生体に光を照射し、入射した光に基づいて得られる生体内の情報を画像化する光イメージング装置の研究が医療分野で積極的に進められている。この光イメージング技術の一つとして、Ｐｈｏｔｏ Ａｃｏｕｓｔｉｃ Ｔｏｍｏｇｒａｐｈｙ（ＰＡＴ：光音響トモグラフィー）がある。光音響トモグラフィーでは、光源から発生したパルス光を生体に照射し、生体内で伝播・拡散したパルス光のエネルギーを吸収した生体組織から発生した音響波を検出する（特許文献１）。すなわち、腫瘍などの被検部位とそれ以外の組織との光エネルギーの吸収率の差を利用し、被検部位が照射された光エネルギーを吸収して瞬間的に膨張する際に発生する弾性波、すなわち光音響波をトランスデューサで受信する。この検出信号を解析処理することにより、生体内の光学特性分布、特に、光エネルギー吸収密度分布を得ることができる。これらの情報は、被検体内の特定物質、例えば血液中に含まれるグルコースやヘモグロビンなどの定量的計測にも利用できる。その結果、新生血管の増殖を伴う悪性腫瘍が存在する場所の特定などに利用できる。

また、非特許文献１には、光音響顕微鏡としての適用例が開示されている。非特許文献１によれば、パルス光を被検体に照射し発生する超音波をトランスデューサにて受信し画像化している。さらにパルス光の波長を変化させることによって被検体の分光特性を画像化している。

米国特許第５，８４０，０２３号明細書

ＩＥＥＥ Ｊｏｕｒｎａｌ ｏｆ Ｓｅｌｅｃｔｅｄ Ｔｏｐｉｃｓ ｉｎ Ｑｕａｎｔｕｍ Ｅｌｅｃｔｒｏｎｉｃｓ，Ｖｏｌ．１４，Ｎｏ．１，１７１−１７９（２００８）

ＰＡＴは、局所的な被検部位で吸収されて発生する音響波を測定することで、局所的な光の吸収情報を得ることができる。被検部位で発生する音響波の初期音圧Ｐは、光照射点から被検部位までの距離ｒを用いて、下の式（１）のように表される。

Ｐ（ｄ）＝Γμ_ａ（ｒ）Φ（ｒ） ・・・式（１）
ここで、
Γ：グリュナイゼン係数（熱−音響変換効率）
μ_ａ（ｒ）：距離ｒにおける位置での吸収係数
Φ（ｒ）：距離ｒにおける位置での光強度
である。弾性特性値であるグリュナイゼン係数Γは、熱膨張係数βと音速ｃの二乗の積を定圧比熱Ｃｐで割ったものである。Γは生体組織が決まればほぼ一定の値をとることが知られているので、音響波の大きさである音圧Ｐの変化を時分割で測定することによりμ_ａとΦの積、すなわち、光エネルギー吸収密度分布Ｈを得ることができる。そして、Ｈから光強度Φ（ｒ）を除算することでμ_ａ（ｒ）を求めることができる。

ここで、光音響波を発生させるために用いられるパルスレーザーは、その原理上常に一定光量のパルス光を発生させることはできず、ある程度の経時的な出力変動があることが知られている。具体的には、その光量変動が１０％以上に及ぶこともある。発生するパルス光量が変動すると、被検体内の局所領域における光量Φ（ｒ）も変動する。上記のように、光量Φと光音響波の強度Ｐは比例関係にあるため、光音響波もレーザーパルスごとに同様の変動を生じる。したがって、このような変動を放置して、光エネルギー吸収密度Ｈ分布や吸収係数μ_ａ分布を画像化した場合には、同じく再構成後の画面上に強度むらを生じることとなり、測定の定量性が失われる原因となる。

しかしながら上記特許文献１では、光源からの経時的な出力変動に関する記載はない。また、非特許文献１にはパルス光量を測定するためのセンサを用いて光量を補正することが記載されているものの、その測定、使用法、補正対象については明確に記されていない。特に非特許文献１は光音響顕微鏡としての使用が前提となっているために、厚みのある被検体を測定するときに重要となる深さ方向の光減衰については明確に記されていなかった。

本発明は、このような背景技術及び課題認識に基づいてなされたものである。本発明の目的は、光源の出力光量が経時的に変化する場合でも、それにより生じる画像への悪影響を低減させることができる光音響装置を提供することである。

本発明に係る処理装置は、光が被検体に照射されることにより発生する音響波を検出して得られる検出信号であって、複数回の光照射に対応した前記検出信号に基づく前記被検体内部の画像再構成データを取得する画像再構成データ取得部と、前記被検体に照射される光照射毎の照射光の光量と、予め記憶されている前記光の照射面における面内強度分布に基づいて前記被検体内部の光強度分布を取得し、前記画像再構成データを前記被検体内部の前記光強度分布に基づいて前記画像再構成データを補正する補正処理部と、を備える。
また、別の本発明に係る処理装置は、被検体内部の情報をボクセルデータとして取得するための処理装置であって、被検体への光の照射により生じる音響波を受信して得られる、複数回の光照射に対応する信号を取得する信号取得部、前記被検体へ複数回照射される光の光量変動に関する情報を取得するセンサ部、前記被検体へ一度に照射される照射光の照射面内における強度分布情報を記憶しているメモリ部、及び前記信号取得部により取得される前記信号、前記センサ部により取得される前記光量変動に関する情報、及び前記メモリ部に予め記憶されている前記強度分布情報を用いて、前記被検体に複数回照射される照射光の照射毎の光量変動が前記被検体内部の情報に対応するボクセルデータに与える影響を低減するための補正処理を行う補正処理部、を備える。
更にまた、別の本発明に係る処理方法は、被検体内部の情報をボクセルデータとして取得するための処理方法であって、被検体への光の照射により該被検体から音響波が生じる光照射を複数回行い、該光照射により生じた該音響波を受信して得られる信号を取得する工程、
前記被検体へ複数回照射される光の光量変動に関する情報を取得する工程、及び前記信号、前記光量変動に関する情報、及び予め取得されている前記被検体へ一度に照射される照射光の照射面内における強度分布情報を用いて、前記被検体に複数回照射される照射光の照射毎の光量変動が前記被検体内部の情報に対応するボクセルデータに与える影響を低減するための補正処理を行う工程を有する。

本発明によれば、光音響装置に光量測定器を設け、測定した出力光量の経時変化を考慮して光音響信号の強度を補正するため、光源の出力に経時的な出力変動があっても、画像に悪影響の少ない光音響装置を提供することができる。

本発明の実施形態１による光音響装置の構成を模式的に示した図である。 レーザーパルスの出力変動を説明するための図である。 本発明の実施形態１において、検出信号の強度の補正量を決定する一例を説明するためのフロー図である。 本発明のホトセンサの配置について説明するための図である。 本発明の実施形態２による光音響装置の構成を模式的に示した図である。 本発明の実施形態２において、検出信号の強度の補正量を決定する一例を説明するためのフロー図である。 本発明の実施形態３による光音響装置の構成を模式的に示した図である。

以下、図面を参照しつつ本発明をより詳細に説明する。なお、同一の構成要素には原則として同一の参照番号を付して、説明を省略する。

（実施形態１：光音響装置）
まず、図１を参照しながら本実施形態にかかる光音響装置の構成を説明する。
本実施形態の光音響装置は、被検体の内部の情報を画像化する光音響イメージング装置である。被検体が生体の場合、光音響装置は、悪性腫瘍や血管疾患などの診断や化学治療の経過観察などを目的として、生体情報の画像化を可能とする。本発明における「被検体内部の情報」とは、光照射によって生じた音響波の発生源分布の情報であり、生体内の初期音圧分布の情報、それから導かれる光エネルギー吸収密度分布の情報、それらの情報から得られる生体組織を構成する物質の濃度分布の情報などである。例えば、物質の濃度分布とは酸素飽和度などである。

本実施形態の光音響装置は、基本的なハード構成として、パルスレーザー２ａ、検出器５、ホトセンサ８ａを有する。パルスレーザー２ａは、被検体にパルス光を照射するための光源である。

生体などの被検体３は必要に応じて、これを両側から圧迫固定するプレート４ａ、４ｂに固定される。光源からの光は、例えばレンズ、ミラー、光ファイバなどの光学系（不図示）によって、プレート４ｂ表面に導かれ、被検体に照射される。被検体３の内部を伝播した光のエネルギーの一部が血管などの光吸収体に吸収されると、その光吸収体から熱膨張により音響波（典型的には超音波）が発生する。これは「光音響波」ということもある。すなわち、パルス光の吸収により、光吸収体の温度が上昇し、その温度上昇により体積膨張が起こり、音響波が発生する。このときの光パルスの時間幅は、光吸収体に吸収エネルギーを効率に閉じ込めるために、熱・ストレス閉じ込め条件が当てはまる程度にすることが好ましい。典型的には１ナノ秒から０．２秒程度である。

音響波を検出するための検出器５は、被検体内で発生した音響波を検出し、アナログ信号である電気信号に変換する。この検出器から取得される検出信号は「光音響信号」ともいう。

この光音響信号を処理して被検体内部の情報を取得するための信号処理部１５は、本実施形態においては、受信アンプ６、Ａ／Ｄコンバータ７、信号補正部１１、画像再構成処理部１２、光減衰補正部１６から構成される。検出器５から取得した光音響信号は、受信アンプ６によって増幅され、Ａ／Ｄ変換器７によってデジタル信号としての光音響信号に変換される。信号補正部１１は本実施形態の特徴的構成の一つであり、このデジタル信号に対して強度の補正を行う。そして、画像再構成処理部１２によって三次元情報に演算処理なされた後、得られたボクセルデータに対して、光減衰補正部１６にて被検体内での光減衰を考慮した補正を行う。そして、必要に応じて画像表示部１３にて被検体の光音響像が表示される。また、全ての要素はシステム制御部１によって制御されている。ここで、「光音響像」とは、得られた被検体内部の情報を３次元空間の座標で示すと共に、その情報を輝度情報に変換して表示した像である。

ここで、本実施形態の特徴的部分を簡潔に説明する。レーザー２ａの出力光量は、光量測定器であるホトセンサ８ａによって測定される。レーザー２ａの出力光量に経時的な変動が生じる場合、この変動もホトセンサ８ａにより測定される。そして、信号補正部１１では、出力光量の経時変化に基づいて、光音響信号の強度変動を抑制するように、光音響信号の強度を補正する。すなわち、出力光量の経時変化（時間変動）に起因する光音響信号の強度変動を低減することができる。

（光源及び光源からの出力光量の変動）
パルスレーザー２ａにより発生するレーザー光はパルス毎に変動する。光量変動の一例を図２に示す。図はＹＡＧレーザーで約５Ｗ（５００ｍＪ）にて１０Ｈｚのパルス光を６０秒発生させた場合の測定出力の経時変化を測定したものである。図より、出力光量は１０％程度の光量変動を確認することができる。

被検体が生体の場合、光源からは生体を構成する成分のうち特定の成分に吸収される特定の波長の光を照射する。光源としては１ナノから０．２ナノ秒オーダーのパルス光を発生可能なパルス光源が好ましい。光源としてはレーザーが好ましいが、レーザーのかわりに発光ダイオードなどを用いることも可能である。レーザーとしては、固体レーザー、ガスレーザー、色素レーザー、半導体レーザーなど様々なレーザーを使用することができる。

なお、図２のようなレーザーの出力光量の変動は、レーザー励起光源であるフラッシュランプの光量変動が主要因であると推測される。よって、本発明の光源として、フラッシュランプそのもの、あるいは、これを励起光源として含むレーザーを使用すると、本発明の効果をより得ることができる。しかし、勿論本発明の光源がこれらに限定されることはなく、フラッシュランプを含まない半導体レーザーや発光ダイオードであっても光量変動を生じうる光源であれば、本発明を用いることができる。

なお、本実施の形態においては、単一の光源の例を示しているが、複数の光源を用いても良い。複数光源の場合は、生体に照射する光の照射強度を上げるため、同じ波長を発振する光源を複数用いても良いし、光学特性値分布の波長による違いを測定するために、発振波長の異なる光源を複数個用いても良い。なお、光源として、発振する波長の変換可能な色素やＯＰＯ（Ｏｐｔｉｃａｌ Ｐａｒａｍｅｔｒｉｃ Ｏｓｃｉｌｌａｔｏｒｓ）を用いることができれば、光学特性値分布の波長による違いを測定することも可能になる。使用する波長としては、生体内において吸収が少ない７００ｎｍから１１００ｎｍの波長帯域から選択される。ただし、比較的生体表面付近の生体組織の光学特性値分布を求める場合は、上記の波長帯域よりも範囲の広い、例えば４００ｎｍから１６００ｎｍの波長帯域から選択して使用する。

光源から照射される光は必要に応じて光導波路などを用いて伝搬させることも可能である。図１で示してはいないが、光導波路としては、光ファイバが好ましい。光ファイバを用いる場合は、それぞれの光源に対して、複数の光ファイバを使用して、生体表面に光を導くことも可能であるし、複数の光源からの光を一本の光ファイバに導き、一本の光ファイバのみを用いて、すべての光を生体に導いても良い。また、主に光を反射するミラーや、光を集光したり拡大したり形状を変化させるレンズなどの光学部品で光を導いてもよい。このような光学部品は、光源から発せられた光が被検体表面の光照射領域に所望の形状で照射されれば、どのようなものを用いてもかまわない。

（検出信号の補正１）
以下に、本実施形態における検出信号の補正について詳細に説明する。
図１のように被検体が平板上に固定されており、レーザー２ａによる光照射領域が被検体表面に２次元的に設定されており、光照射領域が画像化範囲に対して十分大きい場合をモデルに説明する。被検体表面に照射されるパルス光の光量をΦ_０とする。被検体内では光は吸収、散乱によって表面から遠ざかるに従い指数関数的に減衰する。すなわち、
Φ（ｒ）＝Φ_０・ｅｘｐ（−μ_ｅｆｆ・ｒ） ・・・式（２）と表現できる。μ_ｅｆｆは被検体の平均的な等価減衰係数である。式（２）と式（１）から、
Ｐ（ｒ）＝Γμ_ａ（ｄ）Φ_０・ｅｘｐ（−μ_ｅｆｆ・ｒ） ・・・式（３）となる。本発明においては、Φ_０がパルスごとに変動することが問題となる。例えば、あるパルス１の出力光量Φ_０１に対して、別のパルス２の出力光量Φ_０２が０．９Φ_０１であった場合、パルス２による光音響波の音圧Ｐ_２（ｒ）＝０．９Ｐ_１（ｒ）となる。

すると、パルス１から作成した被検体内の画像（μ_ａ分布）とパルス２から作成した画像とで、音圧に対応して輝度信号が相違することになり、画像の再現性が得られない。このように、同じ場所を複数回測定する場合に、画像の再現性低下により被検体内の情報の誤認の原因となる。また、レーザー及び検出器を被検体表面に沿って走査しながら測定する場合、複数のパルスに応答した音圧から１つの画像を作成する。この場合、上記の光量変動が原因で画像内に輝度ムラが生じ、これも被検体内の情報の誤認の原因となる。

そこで、本実施形態においては、ホトセンサ８ａで各パルスの出力光量Φ_０ｎを測定する。そして、音圧Ｐ_ｎ（ｒ）が常に一定の初期光量Φ_０のような基準光量によって得られたとみなせるように、それぞれの光音響波による検出信号の強度を補正する。上記の例では、Ｐ_１（ｒ）を基準としたときに、Ｐ_２（ｒ）を１／０．９倍した強度が得られたとして、検出信号を補正する。このようにすれば、光源の出力に経時的な出力変動があっても、その影響を低減して被検体内部の音源の位置と音圧に関する情報を取得することができる。

なお、上記のパルス１の出力光量Φ_０１に対する、別のパルス２の出力光量Φ_０２の比の逆数、つまりΦ_０１／Φ_０２を本明細書において「補正係数」と呼ぶ。また、上記の検出信号の補正はアナログ信号、デジタル信号いずれに対しても行うことは出来るが、本実施形態では、Ａ／Ｄコンバータ７によってデジタル信号に変換された量に対して補正を行っている。デジタル信号はＡ／Ｄコンバータ７よりサンプリング周波数毎に音圧Ｐ（ｒ）の値が出力されるので、あるパルスに対応するＰ（ｒ）を示すデジタル信号に対して、当該パルスの補正係数を乗算処理することによって補正を行う。

以下、さらに詳細に説明する。本実施形態ではレーザー２ａのパルス光量をホトセンサ８ａにてパルス毎に検出した後、光量メモリ９ａに記憶する。信号処理の確実性の観点から、このように出力光量を記憶するメモリを有することが好ましい。補正量決定部１０は、メモリ９ａに記憶された出力光量の経時変化に関するデータを読み出し、それぞれの検出信号に対する補正量（補正係数）を決定する。決定された補正量によって、信号補正部１１でそれぞれの検出信号の強度を補正する。

図３に補正量演算の処理フローを示す。補正量決定部１０は、光量メモリ９ａからのデータを読み込む（Ｓ３０１）。光量メモリ９ａによって得られたパルスの光量に基づいて、それぞれのパルスに対応する補正係数が計算される（Ｓ３０２）。補正係数はあらかじめ測定された光量を基準値として用いる。ここで本明細書において、複数のパルスごとの補正係数セットを「補正量テーブル」という。そして、補正量テーブルは信号補正部１１に転送され（Ｓ３０４）、得られた光音響信号に対して演算される。つまり、サンプリング周波数毎に取得された音圧Ｐ（ｄ）のデジタル信号に対して補正係数が乗算される。

（画像再構成及び光減衰補正）
上記のように補正されたデジタル信号に対して、画像再構成処理部１２により画像再構成を行う。ＰＡＴの画像再構成とは、検出器で受信された音圧Ｐ_ｄ（ｒ_ｄ，ｔ）から被検体内の発生したときの初期音圧の分布Ｐ_０（ｒ）を導き出すことであり、数学的には逆問題と呼ばれる。ＰＡＴの画像再構成手法で代表的に使われているＵｎｉｖｅｒｓａｌ Ｂａｃｋ Ｐｒｏｊｅｃｔｉｏｎ（ＵＢＰ）法などはＰＨＹＳＩＣＡＬ ＲＥＶＩＥＷ Ｅ７１，０１６７０６（２００５）やＲＥＶＩＥＷ ＯＦ ＳＣＩＥＮＴＩＦＩＣ ＩＮＳＴＲＵＭＥＮＴＳ，７７，０４２２０１（２００６）に記載されている。

ここまでで、被検体内部の情報として初期音圧分布、あるいはμ_ａとΦの積、すなわち、光エネルギー吸収密度分布Ｈを得ることができた。Φを一定として、ＨをΦで除算すれば、被検体内の吸収係数μ_ａ（ｒ）分布を求めることができる。しかし、被検体内の局所領域に照射される光量は、上記のように指数関数的に減衰するため、同じ吸収係数を有する組織が２つあった場合、被検体表面からより遠い組織から発生する音響波の音圧は、近い組織からのそれより小さくなる。そこで、正確な吸収係数分布を得るため、このような光減衰の影響も補正することが好ましい。この補正を本明細書では「光減衰補正」という。

具体的には、画像再構成処理部１２から出力される光エネルギー吸収密度分布Ｈを示すボクセルデータのそれぞれに対して、該当するボクセルの位置における光量を除算する処理を、光減衰補正部１６で行う。該当するボクセルの位置における光量は、上記式（２）より算出することができる。

このような構成とすれば、光減衰の影響も考慮して被検体内の正確な吸収係数分布を画像化することができる。

（信号補正１の態様）
なお、本実施形態ではパルスごとの光量変動に対する補正を、Ａ／Ｄコンバータ７によってデジタル信号になった後、画像再構成処理を行う前のデータに対して行う例を示したが、本発明はこれに限定されない。この補正は、画像再構成後のボクセルデータに対して行うこともできる。つまり補正量決定部１０によって計算された補正係数を光減衰補正部１６に入力することによって補正演算を行っても同様に光音響信号の補正が可能である。また、Ａ／Ｄ変換する前のアナログデータに対して、補正を行うこともできる。この場合、補正量決定部１０の出力によって、受信アンプ６のゲインを制御する方法を取ることができる。すなわち、本明細書における「検出信号」には、アナログ信号、Ａ／Ｄ変換後のデジタル信号、このデジタル信号を画像再構成した後の輝度データのいずれも含まれる。

また、ある光が照射されたときの被検体内の光量分布が設定できれば、プレート４の両側より光を入射した場合、またさらに多方向からの光照射がなされた場合においても同様に光音響信号の補正が可能である。

また、上記の形態においてはトランスデューサ５及びレーザー２ａは固定された状態にて光音響信号の取得を行ったが、トランスデューサ５及びレーザー２ａを走査しながら光音響信号を取得する場合においても、各走査位置における光量測定データを取得することによって同様の光音響信号の補正が可能である。

また、使用するレーザーパルス２ａはある太さを持つレーザービームであるが、その断面の強度に空間的なむらが存在する場合には、その光量分布も考慮した上で三次元空間内にて補正光量の演算を行うことによって同様の光音響信号の補正が可能である。

（各構成の詳細な説明）
検出器（探触子）５は、音波や超音波のような音響波を検知し、電気信号に変換するものである。圧電現象を用いたトランスデューサ、光の共振を用いたトランスデューサ、容量の変化を用いたトランスデューサなど音響波信号を検知できるものであれば、どのような音響波検出器を用いてもよい。本実施形態の検出器５は、検出器５は複数のトランスデューサ素子を持つアレイ型であることが好ましい。このような２次元に配列されたトランスデューサ素子を用いることで、同時に複数の場所で音響波を検出することができる。これにより検出時間を短縮できると共に、被検体の振動などの影響を低減できる。また、検出器５とプレート４ｂ、被検体３との間には、音波の反射を抑えるためのジェルや水などの音響インピーダンスマッチング剤を使うことが望ましい。

光量測定器としては、ホトダイオードに代表されるホトセンサやパイロエレクトリックセンサが代表的なものであるが、１次元又は２次元のホトセンサアレイが必要な場合は、ＣＣＤのイメージセンサや、ＣＭＯＳイメージセンサ、また光依存性抵抗（ＬＤＲ）等を使用しても同様の効果を得ることが出来る。

光量測定器としてのホトセンサ８ａの好ましい配置について、図４を説明して説明する。１８は反射ミラー、１９はレーザーを示す。

図４（ａ）は、プレート４ｂに到達する前の光学系の途中において、反射ミラーの漏れ光を検出する場合である。このように、反射ミラーの後方にホトセンサ８ａを配置することで、光源からの光の一部を検出することができる。漏れ光の比率を既知としておけば、光源の光量変動を算出することは可能である。

図４（ｂ）は、プレート４ｂから反射した一部の光を検出する場合である。このように、プレート４ｂ近傍にホトセンサ８ａを配置することも好ましい。

図４（ｃ）は、プレート４ｂ内部を伝播した一部も光を検出する場合である。このように、プレート４ｂの終端部にホトセンサ８ａを設けることもできる。特に、プレートに対して斜めに光が入射する場合、このようなプレート内部の伝播光が増加するため、好ましい形態である。

メモリとしては、ＰＣもしくは制御ボード上のメモリが考えられるが、ホトセンサユニットに附属のメモリ、もしくはレーザーパルス周期以上の速度があればハードディスク等を用いても同様の効果を得ることが出来る。

（実施例１）
以下、実施例１では本発明における光音響装置を乳房検査に使用した場合について具体的に述べる。本実施例における乳房検査では、一般に用いられるＸ線マンモグラフィーと同様の乳房圧迫を行うものとする。つまり乳房においては平均的な乳房圧迫厚さである深さ４ｃｍまでの光音響信号を取得することができれば良い。

本実施例では、光源として波長１０６４ｎｍ、１０ヘルツ駆動、パルス幅５ナノ秒、１パルスあたりの出力が１．６ＪのＱスイッチＹＡＧレーザーを使用した。この条件で人体に照射できるレーザー光強度は１００ｍＪ／ｃｍ２以下であることがＪＩＳに定められているため、射出されるレーザー光を４ｃｍ角に広げるよう照明光学系を設計した。

これにより乳房を圧迫し両側から光を照射した場合に、光音響信号を発生させる範囲を深さ４ｃｍ、幅４ｃｍとした。またこの範囲内の光音響信号を取得するため、超音波トランスデューサの１辺は４ｃｍとした。また、素子ピッチは２ｍｍとし４００素子の二次元探触子を構成した。そして使用周波数１ＭＨｚとした。ホトセンサは浜松ホトニクス社製Ｓ５９７３ ＰＩＮフォトダイオードを用いた。

以上の条件にて光音響信号を取得するにあたっては、レーザー照射の際の照射光量は常に光量メモリに蓄えられる。各パルスに対応する補正係数は、測定された光量変動の最大値を１として計算される。図１の構成で説明した補正方法に従い検出信号の補正を行った。画像再構成を行い生成された光音響画像をボリュームデータとして蓄積し画面に表示した。

本方法を用いることによって光音響像の強度むらを改善することが可能であるが、電気ノイズに基づく光音響信号自体の再現性が２乃至３％程度、照射光の分布むらが２％程度存在するために同程度の画像むらは残る。しかしながら本方法を用いることによって従来８乃至１０％程度あった画像むらを３乃至４％程度にまで改善することが可能となった。

なお、照射光の分布は光路内にＣＣＤ等を別途配置することによって測定することが出来、光音響信号の再現性は本システムを、レーザー照射を行わない状態にて動作させることによって測定することが出来る。画像むらは複数回画像取得を行った際の同一ピクセルにおける輝度値変化の標準偏差の３倍と定義する。

なお、本実施例においては光音響装置を乳房検査に使用した場合について具体的に述べたが、人体のその他の部位さらに人体以外の被検体の測定についても同様の処理をすることによって同様の効果を得ることが可能である。

（実施形態２）
実施形態１においてはレーザー光を片側からのみ照射したが、本実施形態ではレーザー光をプレート４の両側から照射した場合の補正方法について述べる。図５に本実施形態を説明するための光音響装置の一例を示す。レーザーパルスの照射、光量データ及び光音響信号の取得は実施形態１と同様である。実施形態１と異なる点はレーザー光をレーザー２ａ，２ｂの両側から照射し、レーザー２ａ，２ｂの出力を２つのホトセンサ８ａ，８ｂにて検出する点である。検出された光量は光量メモリ９ａ，９ｂに蓄えられ、補正量決定部１０ａ、１０ｂに転送される。そして実施形態１と同様の方法を用いて２つの補正係数が算出され、信号補正部１１に補正係数が転送される。

図６を用いて、補正量決定部１０ａ，ｂ及び信号補正部１６の動作を詳述する。光量メモリ９ａ，９ｂによって得られた（Ｓ６０１）各レーザーからの光量はあらかじめ測定された各々の光量変動の最大値を１として補正係数が計算される（Ｓ６０２）。そして規格化された各光量と、各係数を用い深さ方向の相対減衰量を計算し、２つの光量データより２つの補正量テーブルが作成される（Ｓ６０３）。これら深さ方向の補正量テーブルはそれぞれ同じ深さに対して加算処理され、合成された一つの補正量テーブルが生成される（Ｓ６０４）。

つまり、プレート４の一方の側から照射された光による音圧は前述の式（３）と同様
Ｐ（ｒ）＝Γμ_ａ（ｒ）Φ_０Ａ・ｅｘｐ（−μ_ｅｆｆ・ｒ） ・・・式（３）であるが、もう一方の側から照射された光による音圧は次の式で表される。
Ｐ（ｒ）＝Γμ_ａ（ｒ）Φ_０Ｂ・ｅｘｐ（−μ_ｅｆｆ・（Ｄ−ｒ）） ・・・式（４）

ここでＤは圧迫プレート４ａ、４ｂ間の距離である。また、Φ_０Ａは、当該音響波を生じさせたパルス光のうち、レーザー２ａが発したパルス光について、光量変動の補正係数を掛けた後の初期光量である。Φ_０Ｂは、当該音響波を生じさせたパルス光のうち、レーザー２ｂが発したパルス光について、光量変動の補正係数を掛けた後の初期光量である。１つの補正量テーブルは上式の光量と光の減衰分をサンプリング周波数に対して数値化し、加算されたものであるため次の式で表される。
Ｃ（ｒ）＝１／（Φ_０Ａ・（ｅｘｐ（−μ_ｅｆｆ・ｒ）＋Φ_０Ｂｅｘｐ（−μ_ｅｆｆ・（Ｄ−ｒ）））） ・・・式（５）

ここでＣ（ｒ）は補正量テーブルの各値であり、Φ_０Ａ、Φ_０Ｂ、ｒ、Ｄ、μ_ｅｆｆは既知であるため求めることが出来る。そして合成された補正量テーブルＣ（ｄ）は、得られた画像再構成データに対して乗算される（Ｓ６０５）。この際、補正量テーブルは深さ方向（被検体表面からの直進距離）の関数であるのに対して、画像再構成データは３次元データであるが、深さ方向に対してのみ乗算処理を行い、広がり方向（被検体表面の面内方向）にはすべて同様の処理を行う。その後、光音響画像をボリュームデータとして蓄積し画面に表示する。この手法は、画像再構成後のデジタルデータに対して、光源からの光量変動を考慮する補正と、深さ方向の光減衰を考慮する補正をまとめて行う手法である。この場合の検出信号は、輝度信号に変換された音圧信号でありうる。

本実施形態を用いることによって、板４の両側から光照射を行う場合であっても実施形態１と同様に画像むらを改善することができる。

（実施形態３）
実施形態１，２においては被検体表面には場所に依存せず一定の光量Φ_０が照射されているとして光音響信号の強度を補正したが、実施形態３ではレーザー照射面（被検体表面）において初期光量に強度分布がある場合の補正方法について述べる。

図７に本実施形態を説明するための光音響装置の一例を示す。レーザーパルスの照射、光量データ及び光音響信号の取得は実施形態２と同様である。被検体３に照射するレーザーはトランスデューサ５と同じ４ｃｍ角に拡大され照射されるが、その面内には強度分布が存在する。これは、マルチモードレーザーを使用する場合に特に見られる。その面内強度分布は光量分布メモリ１４ａ，１４ｂにそれぞれあらかじめ測定され、蓄えられている。

ここで補正される補正量は光の強度分布を反映したものとなる。つまり補正量は光の強度分布、グリューナイセン係数、吸収係数、補正係数を掛けた初期光量、減衰係数の関数となるが、光の強度分布は解析的に解くことは出来ないため光伝播シミュレーションによって計算される。計算によって三次元方向の光量分布が求められるが、式を用いて説明すると、式（３）は今回の場合は
Ｐ（ｒ）＝Γμ_ａ（ｄ）Φ_０（ｘ、ｙ、ｒ、ｔ）・ｅｘｐ（−μ_ｅｆｆ・ｒ） ・・・式（６）
であり、光量Φ_０（ｘ、ｙ、ｒ、ｔ）は空間（ｘ、ｙ、ｒ）分布と時間ｔ毎に発振される各パルスの関数である。この光量Φ_０がシミュレーションによって算出され、各係数を乗算した後、逆数を取ることによって三次元空間内の各場所における補正係数が計算される。

本実施形態による補正量演算は、時間的な光量変動に加えて、被検体表面での空間的な照射光量分布も考慮して信号補正できる点が優れている。光音響信号の再現性は実施例１と同様である。これらの要因を考慮して光音響像の取得、補正を行った所、光量変動補正、光量分布補正の無い形態で８乃至１０％あった光音響像の強度むらを３％程度、すなわち光音響信号の再現性と同程度にまで改善することができた。

（実施形態４）
実施形態１〜３においてはトランスデューサ５及びレーザー２ａは固定された状態で光音響信号の取得を行ったが、本実施形態ではこれらをプレートに沿って走査しながら光音響信号を取得する場合について述べる。

システム構成は図４に示された通りであるが、本実施形態ではレーザー２ａ，２ｂ、トランスデューサ５を被検体に対して走査する可動機構を有する。なお、光源から被検体への入射部位が走査されればよく、必ずしもレーザーそのものを走査する必要はない。この場合、実施形態１〜３では考慮されていなかった走査中の光量変動を補正する必要が生ずる。

以上の光源及び超音波トランスデューサを用いて乳房を走査することを考えると、走査領域は２０ｃｍ×２０ｃｍであり、ストライプ幅は４ｃｍであるために４ｃｍ×２０ｃｍのストライプを５回形成することとした。また、トランスデューサの走査はいわゆるステップアンドリピートと呼ばれる方式、すなわちトランスデューサの移動と停止を繰り返し、トランスデューサが停止している際にレーザー照射を行う方式を採用し、光音響信号を取得した。一方、レーザー照射に際してはその照射光量は常に光量メモリに蓄えられる。

光量の規格化、演算処理、及び光音響信号の補正は実施形態２と同様であるが、補正された光音響信号は走査終了まで一旦メモリに蓄えられる。そして全ての走査が終了した後に補正された光音響信号を用い画像再構成を行い、生成された光音響画像をボリュームデータとして蓄積し画面に表示する。

本方法を用いることによって光音響像の強度むらを改善することが可能であるが、光音響信号自体の再現性が２乃至３％程度存在し、一方走査による光音響信号の平均化効果はＳ／Ｎを４乃至５倍高める効果がある。これらの影響を考慮すると従来の走査、光量補正の無い形態で８乃至１０％程度あった画像むらを１％以下に改善することが可能となった。

（実施形態５）
また、本発明は、以下の処理を実行することによっても実現される。即ち、上述した実施形態の機能を実現するソフトウェア（プログラム）を、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータ（またはＣＰＵやＭＰＵ等）がプログラムを読み出して実行する処理である。

１ システム制御部
２ａ，２ｂ レーザー
３ 被検体
４ａ、４ｂ プレート
５ 検出器
６ 受信アンプ
７ アナログ−デジタルコンバータ
８ａ，８ｂ ホトセンサ
９ａ，９ｂ 光量メモリ
１０ 補正量演算部
１１ 信号補正部
１２ 画像再構成処理部
１５ 信号処理部
１６ 光減衰補正部

Claims (14)

1. 光が被検体に照射されることにより発生する音響波を検出して得られる検出信号であって、複数回の光照射に対応した前記検出信号に基づく前記被検体内部の画像再構成データを取得する画像再構成データ取得部と、
   前記被検体に照射される光照射毎の照射光の光量と、予め記憶されている前記光の照射面における面内強度分布に基づいて前記被検体内部の光強度分布を取得し、前記被検体内部の前記光強度分布に基づいて前記画像再構成データを補正する補正処理部と、を備える
   ことを特徴とする処理装置。
2. 前記補正処理部は、前記被検体に照射される光の該被検体の深さ方向における光の減衰に基づいて前記被検体内部の前記光強度分布を取得することを特徴とする請求項１に記載の処理装置。
3. アナログ信号である前記検出信号をデジタル信号に変換するＡ／Ｄ変換器を更に有し、前記画像再構成データ取得部は、前記デジタル信号に基づいて前記画像再構成データを取得することを特徴とする請求項１または２に記載の処理装置。
4. 前記画像再構成データ取得部は、前記検出信号を前記被検体内部の初期音圧分布あるいは光エネルギー吸収密度分布に変換することにより前記画像再構成データを取得することを特徴とする請求項１から３のいずれか１項に記載の処理装置。
5. 被検体内部の情報をボクセルデータとして取得するための処理装置であって、
   被検体への光の照射により生じる音響波を受信して得られる、複数回の光照射に対応する信号を取得する信号取得部、
   前記被検体へ複数回照射される光の光量変動に関する情報を取得するためのセンサ部、
   前記被検体へ一度に照射される照射光の照射面内における強度分布情報を記憶しているメモリ部、及び
   前記信号取得部により取得される前記信号、前記センサ部により取得される前記光量変動に関する情報、及び前記メモリ部に予め記憶されている前記強度分布情報を用いて、前記被検体に複数回照射される照射光の照射毎の光量変動が前記被検体内部の情報に対応するボクセルデータに与える影響を低減するための補正処理を行う補正処理部、を備えることを特徴とする処理装置。
6. 前記補正処理部は、前記被検体に照射された光の該被検体の深さ方向の光減衰の影響も考慮した補正を行うことを特徴とする請求項５に記載の処理装置。
7. 前記補正処理部は、前記信号をデジタル信号に変換する前のアナログ信号に対して、前記光量変動に関する情報を用いて、前記補正処理を行う請求項５あるいは６に記載の処理装置。
8. 前記補正処理部は、前記信号をデジタル信号に変換した後、前記光量変動に関する情報を用いて、前記デジタル信号を補正する処理を行うことを特徴とする請求項５あるいは６に記載の処理装置。
9. 前記補正処理部は、前記信号をデジタル信号に変換した後、画像再構成処理を行うことで得られるボクセルデータに対して、前記光量変動に関する情報を用いた補正を行うことを特徴とする請求項５あるいは６に記載の処理装置。
10. 被検体内部の情報をボクセルデータとして取得するための処理方法であって、
    被検体への光の照射により該被検体から音響波が生じる光照射を複数回行い、該光照射により生じた該音響波を受信して得られる信号を取得する工程、
    前記被検体へ複数回照射される光の光量変動に関する情報を取得する工程、及び
    前記信号、前記光量変動に関する情報、及び予め取得されている前記被検体へ一度に照射される照射光の照射面内における強度分布情報を用いて、前記被検体に複数回照射される照射光の照射毎の光量変動が前記被検体内部の情報に対応するボクセルデータに与える影響を低減するための補正処理を行う工程を有することを特徴とする処理方法。
11. 前記補正処理を行う工程が、前記被検体に照射された光の該被検体の深さ方向の光減衰の影響も考慮した補正を行う工程であることを特徴とする請求項１０に記載の処理方法。
12. 前記補正処理を行う工程は、前記信号をデジタル信号に変換する前のアナログ信号に対して、前記光量変動に関する情報を用いて、前記補正処理を行う請求項１０あるいは１１に記載の処理方法。
13. 前記補正処理を行う工程が、前記信号をデジタル信号に変換した後、前記光量変動に関する情報を用いて、前記デジタル信号を補正する処理を含むことを特徴とする請求項１０あるいは１１に記載の処理方法。
14. 前記補正処理を行う工程が、前記信号をデジタル信号に変換した後、画像再構成処理を行うことで得られるボクセルデータに対して、前記光量変動に関する情報を用いた補正を行うことを特徴とする請求項１０あるいは１１に記載の処理方法。

Images