JP6411667B2 - 光音響計測装置および光音響計測装置の信号処理方法 - Google Patents

Description

本発明は、被検体内で発生した光音響信号に基づいて光音響画像を生成する光音響計測装置、並びにそれに利用される信号処理方法に関する。

近年、光音響効果を利用した非侵襲の計測法が注目されている。この計測法は、ある適宜の波長（例えば、可視光、近赤外光または中間赤外光の波長帯域）を有するパルス光を被検体に向けて出射し、被検体内の吸収物質がこのパルス光のエネルギーを吸収した結果生じる弾性波である光音響波を検出して、その吸収物質の濃度を定量的に計測するものである。被検体内の吸収物質とは、例えば血液中に含まれるグルコースやヘモグロビンなどである。また、このような光音響波を検出しその検出信号に基づいて光音響画像を生成する技術は、光音響イメージング（ＰＡＩ：Photoacoustic Imaging）あるいは光音響トモグラフィー（ＰＡＴ：Photo Acoustic Tomography）と呼ばれている。

光音響イメージングでは、例えば特許文献１および２に示されるように、被検体に向けて測定光を出射する光出射部と、この光出射部と並設され、測定光の出射によって被検体内に生じた光音響波を検出する音響波検出部とを有するプローブが多く使用される。

特開２０１３−１５８５３１号公報 特表２０１５−５１９１８３号公報

本願発明者らは、上述のように光出射部と音響波検出部とが並設されてなるプローブを使用した場合、測定光が入射する被検体の表面部分で生じた光音響波が、被検体表面を伝って他の信号と一緒に音響波検出部により検出され、アーチファクト（偽画像）の発生原因となることを見出した。このようなアーチファクトは、表示された光音響画像の診断性能を損なうものである。

特許文献１には、被検体の深さ方向に相異なる領域ごとに、光音響波検出信号に対する処理を変えてアーチファクトを抑制する技術が記載されている。また特許文献２には、光音響波検出信号に対して低周波カット処理を施すことにより、アーチファクトを除去する技術が記載されている。

しかし、特許文献１に記載された技術は、エラストグラフィー測定やＢモード測定を前提とするものであり、それらの測定では、測定光が入射する被検体の表面部分で生じた光音響波に起因するアーチファクト発生領域を検出することができない。したがって特許文献１に記載された技術では、測定光が入射する被検体の表面部分で生じた光音響波に起因するアーチファクトを精度良く低減あるいは除去することは不可能である。

また、特許文献２に記載された技術は、アーチファクトが発生する領域を認識しないで低周波カット処理を行っているので、この技術でも、測定光が入射する被検体の表面部分で生じた光音響波に起因するアーチファクトを精度良く低減あるいは除去することは不可能である。

本発明は上記問題に鑑みてなされたものであり、測定光が入射する被検体の表面部分で生じた光音響波に起因するアーチファクトを精度良く低減あるいは除去できる光音響計測装置、並びにそれに利用される信号処理方法を提供することを目的とする。

本発明による光音響計測装置は、
被検体に測定光を出射する光出射部と、光出射部に並設され、測定光の出射によって被検体内に生じた光音響波を検出する音響波検出部とを有するプローブ、
音響波検出部が出力した光音響波検出信号に基づいて光音響画像を生成する画像生成部、
少なくとも光出射部と音響波検出部との位置関係に基づいて、光音響画像中のアーチファクト発生領域とアーチファクト非発生領域とを判別する領域判別手段、および、
上記光音響波検出信号のうち、領域判別手段が判別したアーチファクト非発生領域の光音響画像に対応する第一の光音響波検出信号に対して第一のフィルタリング処理を行い、かつ、上記光音響波検出信号のうち、領域判別手段が判別したアーチファクト発生領域の光音響画像に対応する第二の光音響波検出信号に対して第二のフィルタリング処理を行うフィルタ手段を有し、
第二のフィルタリング処理は、第一のフィルタリング処理に比べて、所定の周波数よりも低い周波数範囲の光音響波検出信号をより多く低減することを特徴とするものである。

なお、上記「光音響波検出信号」、「第一の光音響波検出信号」および「第二の光音響波検出信号」は、光音響波を検出して音響波検出部から出力された信号そのものだけでなく、その信号が何らかの処理を受けた後の信号（デジタル化されたデータ等も含む）も含んで意味するものである。

また、上記「所定の周波数」とは、或る値の一定の周波数を意味するものではなく、予め定められた周波数を意味するものであって、その値自体は任意に設定され得るものである。

本発明の光音響計測装置においては、
フィルタ手段がバンドパスフィルタリング処理を行うバンドパスフィルタを含み、
第二のフィルタリング処理としてのバンドパスフィルタリング処理が、第一のフィルタリング処理としてのバンドパスフィルタリング処理よりも、低周波数側のカットオフ周波数をより高周波数側に設定したものであることが望ましい。

また、本発明の光音響計測装置においては、
フィルタ手段がハイパスフィルタリング処理を行うハイパスフィルタを含み、
第二のフィルタリング処理としてのハイパスフィルタリング処理が、第一のフィルタリング処理としてのハイパスフィルタリング処理よりも、カットオフ周波数をより高周波数側に設定したものであってもよい。

また、本発明の光音響計測装置においては、領域判別手段が、光出射部と音響波検出部との位置関係と、アーチファクト発生領域およびアーチファクト非発生領域の境界とを対応させた対応表に基づいて、アーチファクト発生領域とアーチファクト非発生領域とを判別するものであることが望ましい。

また、本発明の光音響計測装置においては、領域判別手段が、光出射部と音響波検出部との位置関係から、アーチファクト発生領域およびアーチファクト非発生領域の境界を算出する計算式により、アーチファクト発生領域とアーチファクト非発生領域とを判別するものであってもよい。

また、本発明の光音響計測装置においては、領域判別手段が、光音響波が伝搬する被検体の部位における音速に応じて上記判別を補正するものであることが望ましい。

また、本発明の光音響計測装置においては、フィルタ手段が、アーチファクト非発生領域およびアーチファクト発生領域の境界を含む範囲内に設定された少なくとも１つの境界領域に関する光音響波検出信号に対して、第一のフィルタリング処理および第二のフィルタリング処理とは異なる境界領域用フィルタリング処理を行うものであることが望ましい。

上記の境界領域用フィルタリング処理は、第一のフィルタリング処理と第二のフィルタリング処理との中間的なフィルタリング処理であることが望ましい。

さらに、上記境界領域用フィルタリング処理の通過特性は、第一のフィルタリング処理の通過特性と、第二のフィルタリング処理の通過特性とに基づいて、境界領域の被検体深さ方向位置に応じて定められたものであることが望ましい。

また、本発明の光音響計測装置は、アーチファクト非発生領域およびアーチファクト発生領域の境界を含む範囲内に設定された少なくとも１つの境界領域に関する画像を、第一のフィルタリング処理を受けた第一の光音響波検出信号と、第二のフィルタリング処理を受けた第二の光音響波検出信号とを、境界領域の被検体深さ方向位置に応じて重み付け加算して得られた信号によって構成する手段をさらに備えていることが望ましい。

一方、本発明による光音響計測装置の信号処理方法は、
被検体に測定光を出射する光出射部と、光出射部に並設され、測定光の出射によって被検体内に生じた光音響波を検出する音響波検出部とを有するプローブ、および、
音響波検出部が出力した光音響波検出信号に基づいて光音響画像を生成する画像生成部を備えた光音響計測装置において、
少なくとも光出射部と音響波検出部との位置関係に基づいて、光音響画像中のアーチファクト発生領域とアーチファクト非発生領域とを判別し、
上記光音響波検出信号のうち、アーチファクト非発生領域の光音響画像に対応する第一の光音響波検出信号に対して第一のフィルタリング処理を行い、
上記光音響波検出信号のうち、アーチファクト発生領域の光音響画像に対応する第二の光音響波検出信号に対して第二のフィルタリング処理を行い、
第二のフィルタリング処理を、第一のフィルタリング処理に比べて、所定の周波数よりも低い周波数範囲の光音響波検出信号をより多く低減するものとすることを特徴とするものである。

なお、本発明による光音響計測装置の信号処理方法においては、
フィルタリング処理がバンドパスフィルタリング処理であり、
第二のフィルタリング処理としてのバンドパスフィルタリング処理が、第一のフィルタリング処理としてのバンドパスフィルタリング処理よりも、低周波数側のカットオフ周波数をより高周波数側に設定したものであることが望ましい。

また、本発明による光音響計測装置の信号処理方法においては、
フィルタリング処理がハイパスフィルタリング処理であり、
第二のフィルタリング処理としてのハイパスフィルタリング処理が、第一のフィルタリング処理としてのハイパスフィルタリング処理よりも、カットオフ周波数をより高周波数側に設定したものであってもよい。

また、本発明による光音響計測装置の信号処理方法においては、光出射部と音響波検出部との位置関係と、アーチファクト発生領域およびアーチファクト非発生領域の境界とを対応させた対応表に基づいて、アーチファクト発生領域とアーチファクト非発生領域とを判別することが望ましい。

あるいは、本発明による光音響計測装置の信号処理方法においては、光出射部と音響波検出部との位置関係から、アーチファクト発生領域およびアーチファクト非発生領域の境界を算出する計算式により、アーチファクト発生領域とアーチファクト非発生領域とを判別するようにしてもよい。

また、本発明による光音響計測装置の信号処理方法においては、光音響波が伝搬する被検体の部位における音速に応じて上記判別を補正することが望ましい。

また、本発明による光音響計測装置の信号処理方法においては、アーチファクト非発生領域およびアーチファクト発生領域の境界を含む範囲内に設定された少なくとも１つの境界領域に関する光音響波検出信号に対して、第一のフィルタリング処理および第二のフィルタリング処理とは異なる境界領域用フィルタリング処理を行うことが望ましい。

上記の境界領域用フィルタリング処理は、第一のフィルタリング処理と第二のフィルタリング処理との中間的なフィルタリング処理であることが望ましい。

そして上記境界領域用フィルタリング処理の通過特性は、第一のフィルタリング処理の通過特性と、第二のフィルタリング処理の通過特性とに基づいて、境界領域の被検体深さ方向位置に応じて定められたものであることが望ましい。

また、本発明による光音響計測装置の信号処理方法においては、アーチファクト非発生領域およびアーチファクト発生領域の境界を含む範囲内に設定された少なくとも１つの境界領域に関する画像を、第一のフィルタリング処理を受けた第一の光音響波検出信号と、第二のフィルタリング処理を受けた第二の光音響波検出信号とを、境界領域の被検体深さ方向位置に応じて重み付け加算して得られた信号によって構成することが望ましい。

本発明による光音響計測装置は、
被検体に測定光を出射する光出射部と、光出射部に並設され、測定光の出射によって被検体内に生じた光音響波を検出する音響波検出部とを有するプローブ、
音響波検出部が出力した光音響波検出信号に基づいて光音響画像を生成する画像生成部、
少なくとも光出射部と音響波検出部との位置関係に基づいて、光音響画像中のアーチファクト発生領域とアーチファクト非発生領域とを判別する領域判別手段、および、
上記光音響波検出信号のうち、領域判別手段が判別したアーチファクト非発生領域の光音響画像に対応する第一の光音響波検出信号に対して第一のフィルタリング処理を行い、かつ、上記光音響波検出信号のうち、領域判別手段が判別したアーチファクト発生領域の光音響画像に対応する第二の光音響波検出信号に対して第二のフィルタリング処理を行うフィルタ手段を有し、
第二のフィルタリング処理は、第一のフィルタリング処理に比べて、所定の周波数よりも低い周波数範囲の光音響波検出信号をより多く低減するように構成されているので、
測定光が入射する被検体の表面部分で生じた光音響波に起因するアーチファクトを精度良く低減あるいは除去可能となる。

本発明の一実施形態による光音響計測装置の全体構成を示す概略図 偽画像の発生を説明する図 （Ａ）〜（Ｄ）は、被検体表面で生じた光音響波に起因してアーチファクト領域が発生した光音響画像を示す図、（Ｅ）および（Ｆ）はそれぞれ本図の（Ａ）および（Ｄ）の光音響画像を生成した際の光出射部と振動子アレイの位置関係を示す概略図 光出射部と振動子アレイとの距離と、観察可能深度との関係例を示すグラフ フィルタの通過特性の一例を示す概略図 フィルタの通過特性の別の例を示す概略図 フィルタの通過特性のさらに別の例を示す概略図 フィルタの通過特性のさらに別の例を示す概略図 境界領域を説明する概略図 フィルタの通過特性のさらに別の例を示す概略図 光出射部と振動子アレイとの距離を説明する図 本発明の別の実施形態による光音響計測装置の全体構成を示す概略図

以下、図面を参照して、本発明の実施形態について詳しく説明する。

図１は、本発明の一実施形態による光音響計測装置１０の全体構成を示す概略図である。なお図１において、プローブ１１の形状は概略的に示してある。本実施形態の光音響計測装置１０は、一例として、光音響信号に基づいて光音響画像を生成する機能を有するものであり、図１に概略的に示すように、超音波探触子からなるプローブ１１、超音波ユニット１２、レーザユニット１３、表示部１４および入力部１５等を備えている。以下、それらの構成要素について順次説明する。

プローブ１１は、例えば生体である被検体Ｍに向けて測定光および超音波を出射する機能と、被検体Ｍ内を伝搬する音響波Ｕを検出する機能とを有する。すなわちプローブ１１は、被検体Ｍに対する超音波の出射（送信）、および被検体Ｍで反射して戻って来た反射超音波（反射音響波）の検出（受信）を行うことができる。

本明細書において「音響波」とは、超音波および光音響波を含む用語である。ここで、「超音波」とはプローブにより送信された弾性波およびその反射波を意味し、「光音響波」とは吸収体６５が測定光を吸収することにより発する弾性波を意味する。また、プローブ１１が発する音響波は超音波に限定されるものでは無く、被検対象や測定条件等に応じて適切な周波数を選択してさえいれば、可聴周波数の音響波を用いてもよい。なお被検体Ｍ内の吸収体６５としては、例えば血管、金属部材等が挙げられる。

プローブ１１は一般に、セクタ走査対応のもの、リニア走査対応のもの、コンベックス走査対応のもの等が用意され、それらの中から適宜のものが撮像部位に応じて選択使用される。またプローブ１１には、後述するレーザユニット１３から発せられた測定光であるレーザ光Ｌを、光出射部４０まで導光させる接続部としての光ファイバ６０が接続されている。

プローブ１１は、音響波検出器である振動子アレイ２０と、この振動子アレイ２０を間に置いて、振動子アレイ２０の両側に各々１つずつ配設された合計２つの光出射部４０と、振動子アレイ２０および２つの光出射部４０等を内部に収容した筐体５０とを備えている。

本実施形態において振動子アレイ２０は、超音波送信素子としても機能する。振動子アレイ２０は、図示外の配線を介して、超音波送信用の回路および音響波受信用の回路等と接続されている。

振動子アレイ２０は、電気音響変換素子である超音波振動子が複数、一方向に並設されてなるものである。超音波振動子は、例えば圧電セラミクス、またはポリフッ化ビニリデン（ＰＶＤＦ）のような高分子フィルムから構成された圧電素子である。超音波振動子は、受信した音響波Ｕを電気信号に変換する機能を有している。なお、振動子アレイ２０は音響レンズを含んでもよい。

本実施形態における振動子アレイ２０は、上述の通り、複数の超音波振動子が一次元に配列されてなるものであるが、複数の超音波振動子が二次元に配列されてなる振動子アレイが用いられてもよい。

上記超音波振動子は、上述した通り超音波を送信する機能も有する。すなわち、この超音波振動子に交番電圧が印加されると、超音波振動子は交番電圧の周波数に対応した周波数の超音波を発生させる。なお、超音波の送信と受信は互いに分離させてもよい。つまり、例えばプローブ１１とは異なる位置から超音波の送信を行い、その送信された超音波に対する反射超音波をプローブ１１で受信するようにしてもよい。

光出射部４０は、光ファイバ６０によって導光されたレーザ光Ｌを被検体Ｍに向けて出射させる部分である。本実施形態において光出射部４０は、光ファイバ６０の先端部、つまり測定光の光源であるレーザユニット１３から遠い方の端部によって構成されている。図１に示されるように、本実施形態では２つの光出射部４０が、振動子アレイ２０を間に置いて、振動子アレイ２０の例えばエレベーション方向の両側に配置されている。このエレベーション方向とは、複数の超音波振動子が一次元に配列された場合、その配列方向に対して直角で、振動子アレイ２０の検出面に平行な方向である。

なお光出射部は、光ファイバ６０の先端に光学的に結合させた導光板および拡散板から構成されてもよい。そのような導光板は、例えばアクリル板や石英板から構成することができる。また拡散板としては、マイクロレンズが基板上にランダムに配置されているレンズ拡散板や、例えば拡散微粒子が分散された石英板等を使用することができる。レンズ拡散板としてはホログラフィック拡散板やエンジニアリング拡散板を用いてもよい。

図１に示されるレーザユニット１３は、例えばＱスイッチアレキサンドライトレーザ等のフラッシュランプ励起Ｑスイッチ固体レーザを有し、測定光としてのレーザ光Ｌを発する。レーザユニット１３は、例えば超音波ユニット１２の制御部３４からのトリガ信号を受けてレーザ光Ｌを出力するように構成されている。レーザユニット１３は、１〜１００ｎｓｅｃ（ナノ秒）のパルス幅を有するパルスレーザ光Ｌを出力するものであることが好ましい。

レーザ光Ｌの波長は、計測の対象となる被検体Ｍ内の吸収体６５の光吸収特性に応じて適宜選択される。例えば計測対象が生体内のヘモグロビンである場合、つまり血管を撮像する場合、一般的にその波長は、近赤外波長域に属する波長であることが好ましい。近赤外波長域とはおよそ７００〜８５０ｎｍの波長域を意味する。しかし、レーザ光Ｌの波長は当然これに限られるものではない。またレーザ光Ｌは、単波長のものでもよいし、例えば７５０ｎｍおよび８００ｎｍ等の複数波長を含むものでもよい。レーザ光Ｌが複数の波長を含む場合、これらの波長の光は、同時に出射されてもよいし、交互に切り替えながら出射されてもよい。

なおレーザユニット１３は、上に述べたアレキサンドライトレーザの他、同様に近赤外波長域のレーザ光を出力可能なＹＡＧ−ＳＨＧ（Second harmonic generation：第二次高調波発生）−ＯＰＯ（Optical Parametric Oscillation：光パラメトリック発振）レーザや、Ｔｉ−Ｓａｐｐｈｉｒｅ（チタン−サファイア）レーザ等を用いて構成することもできる。

光ファイバ６０は、レーザユニット１３から出射されたレーザ光Ｌを、２つの光出射部４０まで導く。光ファイバ６０は特に限定されず、石英ファイバ等の公知のものを使用することができる。例えば１本の太い光ファイバが用いられてもよいし、あるいは複数の光ファイバが束ねられてなるバンドルファイバが用いられてもよい。一例としてバンドルファイバが用いられる場合、１つにまとめられたファイバ部分の光入射端面から上記レーザ光Ｌが入射するようにバンドルファイバが配置され、そしてバンドルファイバの２つに分岐されたファイバ部分の各先端部が前述した通り光出射部４０を構成する。

超音波ユニット１２は、受信回路２１、受信メモリ２２、データ分離手段２３、光音響画像生成部２４、超音波画像生成部２９、表示制御部３０、送信制御回路３３および制御部３４を有する。

制御部３４は、光音響計測装置１０の各部を制御するものであり、本実施形態では図示外のトリガ制御回路を備える。このトリガ制御回路は、例えば光音響画像を取得する場合には、レーザユニット１３に光トリガ信号を送る。これによりレーザユニット１３のＱスイッチ固体レーザにおいて励起源のフラッシュランプが点灯し、レーザロッドの励起が開始される。このレーザロッドの励起状態が維持されている間、レーザユニット１３はレーザ光Ｌを出力可能な状態となる。

制御部３４は、その後トリガ制御回路からレーザユニット１３へＱスイッチトリガ信号を送信する。つまり制御部３４は、このＱスイッチトリガ信号によって、レーザユニット１３からのレーザ光Ｌの出力タイミングを制御する。また制御部３４は、Ｑスイッチトリガ信号の送信と同時に、サンプリングトリガ信号を受信回路２１に送信する。このサンプリングトリガ信号は、受信回路２１のＡＤ変換器（Analog to Digital convertor）における光音響信号のサンプリングの開始タイミングを規定する。このように、サンプリングトリガ信号を使用することにより、レーザ光Ｌの出力と同期して光音響信号をサンプリングすることが可能となる。

制御部３４は、超音波画像を取得する場合は、送信制御回路３３に超音波送信を指示する超音波トリガ信号を送信する。送信制御回路３３は、超音波送信トリガ信号を受けると、プローブ１１から超音波を送信させる。制御部３４は、超音波送信のタイミングに合わせて受信回路２１にサンプリングトリガ信号を送信し、反射超音波信号のサンプリングを開始させる。

以上述べた光音響画像あるいは超音波画像を取得する際、プローブ１１は、被検体Ｍに対して例えば前述したエレベーション方向に少しずつ位置が変えられて、レーザ光Ｌあるいは超音波により被検体Ｍが走査される。そこで上記光音響信号あるいは反射超音波信号のサンプリングはこの走査と同期して、音響波検出ラインを一ラインずつずらしながらなされる。なお上記走査は、術者がプローブ１１を手操作で動かすことによってなされてもよいし、あるいは自動走査機構を用いてなされてもよい。

受信回路２１は、プローブ１１の振動子アレイ２０が出力する光音響波検出信号を受信し、受信した検出信号を受信メモリ２２に格納する。受信回路２１は典型的には、低ノイズアンプ、可変ゲインアンプ、ローパスフィルタ、およびＡＤ変換器を含んで構成される。プローブ１１が出力する光音響波検出信号は、低ノイズアンプで増幅された後に、可変ゲインアンプで深度に応じたゲイン調整がなされ、ローパスフィルタで高周波成分がカットされた後にＡＤ変換器でデジタル信号に変換され、受信メモリ２２に格納される。受信回路２１は、例えば１つのＩＣ（Integrated Circuit）で構成される。なお、上記ローパスフィルタによる高周波成分のカットは、ＡＤ変換時に折り返し雑音が発生することを防止するためのものであり、一般にそのカットオフ周波数は、ＡＤ変換のサンプリング周波数の半分程度の１０ＭＨｚ〜３０ＭＨｚ程度とされる。

本実施形態においてプローブ１１は、光音響波の検出信号と反射超音波の検出信号とを出力する。そこで受信メモリ２２には、デジタル化された光音響波検出信号および反射超音波検出信号が格納される。データ分離手段２３は、受信メモリ２２から光音響画像用のデータ、つまりデジタル化された光音響波検出信号を読み出して、光音響画像生成部２４に送信する。またデータ分離手段２３は、受信メモリ２２から反射超音波画像用のデータ、つまりデジタル化された反射超音波検出信号を読み出して、超音波画像生成部２９に送信する。

光音響画像生成部２４は、受信メモリ２３から受信した光音響波検出信号を再構成して、光音響画像を生成する。具体的には、光音響画像生成部２４は、振動子アレイ２０の各超音波振動子の位置に応じた遅延時間で各超音波振動子からの信号に基づく光音響波検出信号を加算し、１ライン分の光音響波検出信号を生成する（遅延加算法）。光音響画像生成部２４は、遅延加算法に代えて、ＣＢＰ法（Circular Back Projection）により再構成を行ってもよい。あるいは光音響画像生成部２４は、ハフ変換法又はフーリエ変換法を用いて再構成を行ってもよい。再構成された複数ライン分の光音響波検出信号は、検波処理および対数変換処理等の信号処理を経て、光音響画像（断層画像）を表示するための信号として表示制御部３０に送られる。

なお上記複数ライン分の光音響波検出信号は、デジタル化や再構成等の処理が施されたもので、プローブ１１の振動子アレイ２０が出力した光音響波検出信号そのものではないが、この振動子アレイ２０が出力した光音響波検出信号に基づくものであるので、「光音響波検出信号」と称することとする。また本実施形態では、光音響画像生成部２４において前述したフィルタリング処理がなされるが、その処理については後に詳しく説明する。

超音波画像生成部２９は、受信メモリ２２に格納された反射超音波検出信号に対して、基本的に上記光音響波検出信号に対するのと同様の処理を施して、超音波画像（断層画像）を示す複数ライン分の反射超音波検出信号を生成する。超音波画像生成部２９は、そのようにして生成された反射超音波検出信号を表示制御部３０に出力する。

表示制御部３０は、上記光音響波検出信号に基づいて光音響画像を、また上記反射超音波検出信号に基づいて超音波画像を、それぞれ表示部１４に表示させる。これら２つの画像は別々に、あるいは合成されて合成画像として表示部１４に表示される。後者の場合、表示制御部３０は、例えば光音響画像と超音波画像とを重畳させて画像合成を行う。このように、光音響画像に加えて超音波画像を生成し表示させれば、光音響画像では画像化することができない部分を超音波画像において観察可能となる。

次に、以上の述べた通りの基本構成を有する光音響計測装置１０において、前述したアーチファクトの発生を防止する点について説明する。まず図２および図３を参照して、このアーチファクトについて詳しく説明する。

図２は、図１に示す吸収体６５の一例である血管Ｖから生じた光音響波と、測定光が入射する被検体表面で生じた光音響波とが検出される様子を示す概略図である。測定光Ｌを被検体Ｍに照射した際、理想的には血管Ｖからの光音響波Ｕ１のみを検出することが望まれるが、実際には測定光Ｌが入射する被検体表面部分４４で発生した光音響波Ｕ２も検出される。この光音響波Ｕ２がアーチファクト（偽画像）の原因となる。光音響波Ｕ２が発生してから検出されるまでの時間は、光出射部４０と振動子アレイ２０（または個々の超音波振動子。以下同じ）との間隔、より詳細には、測定光Ｌの当接平面上の到達領域と音響波検出部との間隔に依存する。つまり、光出射部４０と振動子アレイ２０との間隔が大きい程、光音響波Ｕ２が被検体中を伝わる距離が長いため、光音響波Ｕ２が発生してから検出されるまでの時間は長くなる。

なお、上記「当接平面」とは、プローブ１１の先端（つまり、被検体Ｍと接触するプローブ表面と音響波検出部の中心軸との交点）を通り音響波検出部の検出面に平行な平面を意味する。また、「到達領域」とは、当接平面と測定光Ｌが交わる領域を意味する。

そして、光音響波Ｕ２が発生してから検出されるまでの時間は、光音響画像中のアーチファクト領域の現れる範囲に影響する。図３の（Ａ）〜（Ｆ）は、光出射部４０と振動子アレイ２０との間隔Ｗ１に依存して、観察領域（アーチファクト非発生領域）Ｒ１およびアーチファクト発生領域Ｒ２の大きさが変化する様子を示す図である。具体的に、図３の（Ａ）〜（Ｄ）は、被検体表面で生じた光音響波Ｕ２に起因してアーチファクト発生領域Ｒ２が生じた光音響画像Ｐを示す図である。また、図３の（Ｅ）および（Ｆ）はそれぞれ、図３の（Ａ）および（Ｄ）の光音響画像Ｐを生成した際の光出射部４０と振動子アレイ２０の位置関係を示す概略図である。

この図３から、光出射部４０と振動子アレイ２０との間隔Ｗ１が大きい程、アーチファクト非発生領域Ｒ１とアーチファクト発生領域Ｒ２との境界Ｂが下がる、つまり被検体Ｍのより深い方に移動することが分かる。これは、光音響画像Ｐの上下方向は、時間軸に対応し、光出射部４０と振動子アレイ２０との間隔Ｗ１が大きい程、光音響波Ｕ２の信号が検出されるまでの時間が遅くなるためである。以上述べた光出射部４０と振動子アレイ２０との間隔Ｗ１に対する、アーチファクト非発生領域Ｒ１とアーチファクト発生領域Ｒ２との境界Ｂの関係の一例を図４に示す。この図４の横軸が上記間隔Ｗ１であり、また縦軸の「観察可能深度」がすなわちアーチファクト非発生領域Ｒ１の深度、つまりは上記境界Ｂの深さ位置である。

以上の知見に基づいて本実施形態では、一例として図１の光音響画像生成部２４において、バンドパスフィルタにより以下のフィルタリング処理がなされる。なお、このフィルタリング処理は、予め与えられたプログラムに基づいて、演算処理によってなされるものである。

このバンドパスフィルタリング処理は、前述したように再構成されて光音響画像を示す光音響波検出信号のうち、アーチファクト非発生領域Ｒ１に関する光音響波検出信号（第一の光音響波検出信号）と、アーチファクト発生領域Ｒ２に関する光音響波検出信号（第二の光音響波検出信号）とでフィルタの通過特性を変えて行われる。ここで、上記第一の光音響波検出信号に対するバンドパスフィルタリング処理を第一のバンドパスフィルタリング処理（本発明における第一のフィルタリング処理）と称し、上記第二の光音響波検出信号に対するバンドパスフィルタリング処理を第二のバンドパスフィルタリング処理（本発明における第二のフィルタリング処理）と称する。

上記バンドパスフィルタリング処理の基本的な通過特性を図５に示す。この図５中、破線で示すのが上記第一のバンドパスフィルタリング処理における通過特性であり、実線で示すのが上記第二のバンドパスフィルタリング処理における通過特性である。つまり、第一の光音響波検出信号に対する低周波数側のカットオフ周波数ｆＬ１よりも、第二の光音響波検出信号に対する低周波数側のカットオフ周波数ｆＬ２の方が、より高周波数側に設定されている。より詳しくは、例えばプローブ１１の振動子アレイ２０の検出中心周波数が６．５ＭＨｚの場合、第一の光音響波検出信号に対する低周波数側のカットオフ周波数ｆＬ１は１ＭＨｚ程度、高周波数側のカットオフ周波数ｆＨは１０ＭＨｚ程度とされ、第二の光音響波検出信号に対する低周波数側のカットオフ周波数ｆＬ２は２ＭＨｚ程度、高周波数側のカットオフ周波数ｆＨは１０ＭＨｚ程度とされる。

以上の通り本実施形態では、第二の光音響波検出信号に対する第二のバンドパスフィルタリング処理は、第一の光音響波検出信号に対する第一のバンドパスフィルタリング処理に比べて、所定の周波数（２ＭＨｚ）よりも低い周波数範囲の光音響波検出信号をより多く低減するものとなっている。ここで、「所定の周波数よりも低い周波数範囲」とは、２ＭＨｚよりも低い周波数範囲に限られるものではなく、光音響波検出信号の最も低周波数側に存在する、或る幅を持った周波数範囲であり、「所定の周波数」の数値は、適宜設定され得るものである。

先に図２を参照して説明したように、被検体表面部分４４で発生してアーチファクトの原因となる光音響波Ｕ２は、レーザ光Ｌが入射する被検体Ｍの表面部分全体から発生するため、この光音響波Ｕ２を検出した信号は比較的低周波数域のものとなる。そこで、上述したように、第一の光音響波検出信号と比べて、第二の光音響波検出信号の低周波数範囲の信号成分をより多く除去すれば、表示部１４に表示される光音響画像のアーチファクト発生領域Ｒ２においてアーチファクトが発生しなくなるか、あるいは低減される。その一方、アーチファクト非発生領域Ｒ１において、被検体Ｍの診断上有意な低周波数範囲の信号成分が不用意に除去されることもなくなる。そこで表示部１４において、全体として診断性能の高い光音響画像を表示可能となる。それに対して、先に述べた特許文献２に記載された技術では、診断上有意な低周波数範囲の信号成分が除去される可能性がある。

なお、上記のアーチファクト非発生領域Ｒ１とアーチファクト発生領域Ｒ２とを判別するためには、それら両領域の境界Ｂを知ることが必要である。先に図３および図４を参照して説明した通り、上記境界Ｂの被検体深さ方向位置は基本的に、光出射部４０と振動子アレイ２０との間隔Ｗ１と対応している。そこで本実施形態では、図１に示した光音響画像生成部２４の内部メモリに、プローブ１１毎の上記境界Ｂの被検体深さ方向位置と上記間隔Ｗ１との対応表をルックアップテーブルの形で記憶させておき、光音響画像生成部２４はこのルックアップテーブルを参照して境界Ｂを、つまりはアーチファクト非発生領域Ｒ１およびアーチファクト発生領域Ｒ２を判別するようにしている。こうすることにより、光音響計測装置１０に対して複数種類のプローブ１１が交換可能に取り付けられる場合でも、使用されるプローブ１１の認識情報を例えば入力部１５から入力して、そのプローブ１１に関して正しくアーチファクト非発生領域Ｒ１およびアーチファクト発生領域Ｒ２を判別可能となる。以上の説明から明らかな通り、本実施形態では光音響画像生成部２４が、本発明における領域判別手段およびフィルタ手段を構成している。

ここで、フィルタ手段としては、ハミング窓、ハニング窓、ブラックマン窓、あるいはカイザー窓等の窓関数によりサイドローブの影響が低減されているものを適用することが望ましい。

なお、上記境界Ｂの被検体深さ方向位置と上記間隔Ｗ１との関係は、観察部位での音速、つまり光音響波が伝搬する被検体Ｍの部位における音速に応じて変化する。そこで、上記対応表における境界Ｂの被検体深さ方向位置と間隔Ｗ１との関係を、上記音速に応じて補正してもよい。この補正により、境界Ｂの被検体深さ方向位置が補正され、つまりはアーチファクト非発生領域Ｒ１とアーチファクト発生領域Ｒ２の判別が補正されることになる。

上記音速は、図１に示す入力部１５から入力するようにしてもよいし、あるいは、この音速と観察部位との対応関係を図示外の記憶手段に記憶させておき、入力部１５から入力された観察部位の情報から上記対応関係を参照して音速を知るようにしてもよい。

さらに、境界Ｂの被検体深さ方向位置と間隔Ｗ１との関係を、上記音速に応じて補正する他、測定光であるレーザＬの強度、レーザＬのパルス幅、および振動子アレイ２０の検出面上における音響部材の有無、のうちの少なくとも１つに基づいて上記関係を補正するようにしてもよい。

また、境界Ｂの被検体深さ方向位置と間隔Ｗ１との対応表から境界Ｂを求める代わりに、間隔Ｗ１から境界Ｂの被検体深さ方向位置を算出する計算式を記憶させておき、図１に示す入力部１５から入力された間隔Ｗ１から上記計算式によって境界Ｂを求めるようにしてもよい。

さらには、操作者が表示部１４に示されるアーチファクト非発生領域Ｒ１とアーチファクト発生領域Ｒ２の表示を確認しながら、境界Ｂを指定する情報を入力部１５から入力するようにしてもよい。

上記境界Ｂを指定する情報とは、境界Ｂの被検体深さ方向位置そのものを示す情報でなくてもよい。例えば図１１に示されるようなｘ（光出射部４０の出射端における光軸の位置から振動子アレイ２０の中心軸までの距離）、ｙ（光出射部４０の出射端における光軸の位置から当接平面Ｓ１までの距離）およびθ（当接平面Ｓ１と測定光の光軸とが成す角度）でもよい。ｘ、ｙおよびθが得られれば、光出射部４０と振動子アレイ２０の間隔Ｗ１が得られ、これからさらに境界Ｂの被検体深さ方向位置を得ることができる。

また、光出射部４０の出射端における光軸が当接平面に対して傾いている場合には、次のようにして「測定光の当接平面上の到達領域と音響波検出部との間隔」を定義することができる。例えば図１１に示すように検出面に何も装着されていない場合には、検出面と振動子アレイ２０の中心軸との交点をＰ１とし、測定光の到達領域におけるエネルギープロファイルＥＰをガウス分布で近似したときの最大点をＰ２とする。この際、上記点Ｐ１が、当接平面が通るプローブの先端となる。そして、この点Ｐ１および点Ｐ２の距離Ｗ２を「到達領域と音響波検出部との間隔」とすることができる。

また、先に述べた対応表や計算式に基づいて定められたアーチファクト非発生領域Ｒ１とアーチファクト発生領域Ｒ２の表示を操作者が表示部１４で確認し、それよりも適正な領域分けが考えられる場合は、その領域分けの境界Ｂを指定する情報を操作者が入力部１５から入力して、アーチファクト非発生領域Ｒ１とアーチファクト発生領域Ｒ２を再設定するようにしてもよい。

なお、以上述べたバンドパスフィルタリング処理により、第一の光音響波検出信号および第二の光音響波検出信号の低周波数範囲の信号成分を除去する代わりに、ハイパスフィルタを用いて上記低周波数範囲の信号成分を除去するようにしてもよい。図６に、そのようなハイパスフィルタの基本的な通過特性例を示す。この図６中、破線で示すのが第一の光音響波検出信号に対する第一のハイパスフィルタリング処理（本発明における第一のフィルタリング処理）の通過特性であり、実線で示すのが第二の光音響波検出信号に対する第二のハイパスフィルタリング処理（本発明における第二のフィルタリング処理）の通過特性である。つまりこの場合は、第一の光音響波検出信号に対するカットオフ周波数ｆ１よりも、第二の光音響波検出信号に対するカットオフ周波数ｆ２の方が、より高周波数側に設定されている。

例えば上記カットオフ周波数ｆ１は１ＭＨｚ程度、カットオフ周波数ｆ２は２ＭＨｚ程度とされる。以上のようにして、この場合も、第一の光音響波検出信号と比べて、第二の光音響波検出信号の低周波数範囲の信号成分をより多く除去すれば、図５に通過特性を示したバンドパスフィルタを用いる場合と、基本的に同様の効果が得られる。

また、バンドパスフィルタを用いる場合、図７に示す通過特性や、あるいは図８に示す通過特性のバンドパスフィルタを適用することもできる。これらの図においても、破線で示すのが第一の光音響波検出信号に対する第一のバンドパスフィルタリング処理の通過特性であり、実線で示すのが第二の光音響波検出信号に対する第二のバンドパスフィルタリング処理の通過特性である。

図７に示す通過特性は、図５に示した通過特性と比べると、第一の光音響波検出信号に対する低周波数側のカットオフ周波数がｆＬ１とｆＬ２との間で、ロールオフ状に連続的に変化している点で異なるものである。この場合も、第一の光音響波検出信号と比べて、第二の光音響波検出信号の低周波数範囲の信号成分をより多く除去できるので、図５に通過特性を示したバンドパスフィルタを用いる場合と、基本的に同様の効果が得られる。

なお図７に示す例では、第一の光音響波検出信号に対する第一のバンドパスフィルタリング処理の低周波数側のカットオフ周波数がｆＬ１とｆＬ２との間で連続的に変化しているが、このカットオフ周波数がｆＬ１とｆＬ２との間で段階的に変化するようにしてもよい。

図８に示す通過特性は、図５に示した通過特性と比べると、第二の光音響波検出信号に対するフィルタリングの利得が、第一の光音響波検出信号に対するフィルタリングの利得よりも高く設定されている点で異なるものである。この場合も、第一の光音響波検出信号と比べて、第二の光音響波検出信号の低周波数範囲の信号成分をより多く除去できるので、図５に通過特性を示したバンドパスフィルタを用いる場合と、基本的に同様の効果が得られる。またこの場合は、第二の光音響波検出信号の低周波数範囲の信号成分をより多く除去することにより、第二の光音響波検出信号の全体的な信号強度が低減する分を、上述のように第二の光音響波検出信号に対するフィルタリングの利得を上げることによって補償可能となる。

なお、以上述べたようにアーチファクト非発生領域Ｒ１に関する光音響波検出信号に対するフィルタリングの特性と、アーチファクト発生領域Ｒ２に関する光音響波検出信号に対するフィルタリングの特性を互いに異なるものとすると、フィルタリング後の光音響波検出信号に基づいて作成かつ表示された光音響画像において、アーチファクト非発生領域Ｒ１とアーチファクト発生領域Ｒ２の境界Ｂで画像の周波数特性が急激に変化しやすい。以下、この周波数特性の急激な変化を抑制するようにした実施形態について説明する。

なお、以下に説明する処理は、例えば図１の光音響画像生成部２４において行われる。この実施形態では、図９に概略的に示すアーチファクト非発生領域Ｒ１とアーチファクト発生領域Ｒ２の境界Ｂを含む範囲内に、一例として２つの境界領域ＢＲ１および境界領域ＢＲ２を設定する。本例では、境界領域ＢＲ１はアーチファクト非発生領域Ｒ１内に属し、境界領域ＢＲ２はアーチファクト発生領域Ｒ２内に属している。そしてこの場合、境界領域ＢＲ１および境界領域ＢＲ２の各光音響画像に対応する光音響波検出信号に対して、前述したバンドパスフィルタによりバンドパスフィルタリング処理が施されるが、それらのバンドパスフィルタリング処理は、アーチファクト非発生領域Ｒ１の光音響画像に対応する第一の光音響波検出信号に対する第一のバンドパスフィルタリング処理とも、またアーチファクト発生領域Ｒ２の光音響画像に対応する第二の光音響波検出信号に対する第二のバンドパスフィルタリング処理とも異なる、境界領域用バンドパスフィルタリング処理とされている。

ここで、上記各バンドパスフィルタリングにおける基本的な通過特性を図１０に示す。この図１０中、破線で示すのが上記第一の光音響波検出信号に対する通過特性、実線で示すのが上記第二の光音響波検出信号に対する通過特性、２点鎖線で示すのが境界領域ＢＲ１の光音響画像に対応する光音響波検出信号に対する通過特性、そして１点鎖線で示すのが境界領域ＢＲ２の光音響画像に対応する光音響波検出信号に対する通過特性である。

以上のように、境界領域ＢＲ１の光音響画像に対応する光音響波検出信号に対する境界領域用バンドパスフィルタリング処理および境界領域ＢＲ２の光音響画像に対応する光音響波検出信号に対する境界領域用バンドパスフィルタリング処理は、アーチファクト非発生領域Ｒ１の光音響画像に対応する第一の光音響波検出信号に対する第一のバンドパスフィルタリング処理と、アーチファクト発生領域Ｒ２の光音響画像に対応する第二の光音響波検出信号に対する第二のバンドパスフィルタリング処理との中間的な処理とされている。

より具体的には、境界領域ＢＲ１用のバンドフィルタリング処理の通過特性の一つである低周波数側のカットオフ周波数ｆＬ３（＝１．３ＭＨｚ）と、境界領域ＢＲ２用のバンドフィルタリング処理の低周波数側のカットオフ周波数ｆＬ４（＝１．７ＭＨｚ）が各々、アーチファクト非発生領域Ｒ１用の第一のバンドパスフィルタリング処理の低周波数側のカットオフ周波数ｆＬ１（＝１ＭＨｚ）と、アーチファクト発生領域Ｒ２用の第二のバンドパスフィルタリング処理の低周波数側のカットオフ周波数ｆＬ２（＝２ＭＨｚ）とに基づいて、境界領域ＢＲ１、境界領域ＢＲ２の被検体深さ方向位置に応じて定められている。

以上のように各バンドパスフィルタリング処理の低周波数側のカットオフ周波数ｆＬ１〜ｆＬ４を設定することにより、アーチファクト非発生領域Ｒ１とアーチファクト発生領域Ｒ２の境界Ｂにおいて、光音響画像の周波数特性が急激に変化することが無くなる。つまり表示部１４に表示される光音響画像においては、アーチファクト非発生領域Ｒ１、境界領域ＢＲ１、境界領域ＢＲ２、およびアーチファクト発生領域Ｒ２の各境界において周波数特性が少しずつ変化するようになり、光音響画像の診断性能が高められる。

なお、以上説明した実施形態においては、アーチファクト非発生領域Ｒ１とアーチファクト発生領域Ｒ２の境界Ｂを含む範囲に、２つの境界領域ＢＲ１および境界領域ＢＲ２を設定しているが、境界領域は１つ以上、いくつ設定しても構わない。そして境界領域を１つだけ設定する場合は、その境界領域の端縁が境界Ｂと揃った状態となっていてもよい。境界領域を多く設定するほど、各領域の境界における周波数特性の変化をより滑らかにすることができるが、それに伴ってフィルタリング処理に要する時間も長くなるので、それらの点を双方考慮して適切な数に設定するのがよい。

次に、アーチファクト非発生領域Ｒ１とアーチファクト発生領域Ｒ２の境界Ｂで画像の周波数特性が急激に変化することを抑制可能とした、本発明の別の実施形態について説明する。なお、以下に説明する処理は、例えば図１の光音響画像生成部２４において行われる。

この実施形態でも、図９に示すように、アーチファクト非発生領域Ｒ１およびアーチファクト発生領域Ｒ２の境界を含む範囲内に、少なくとも１つの境界領域を設定する。本実施形態では図９の例と同じように、２つの境界領域ＢＲ１および境界領域ＢＲ２を設定する。そしてアーチファクト非発生領域Ｒ１の光音響画像に対応する第一の光音響波検出信号に対しては、低周波数側のカットオフ周波数ｆＬ１を１ＭＨｚとした第一のバンドパスフィルタリング処理を施して、第一の処理済み光音響波検出信号を得る。一方、アーチファクト発生領域Ｒ２の光音響画像に対応する第二の光音響波検出信号に対しては、低周波数側のカットオフ周波数ｆＬ２を２ＭＨｚとした第二のバンドパスフィルタリング処理を施して、第二の処理済み光音響波検出信号を得る。

そして上記第一の処理済み光音響波検出信号および第二の処理済み光音響波検出信号を、境界領域ＢＲ１および境界領域ＢＲ２の被検体深さ方向位置に応じて重み付け加算して、境界領域ＢＲ１用の光音響波検出信号、並びに境界領域ＢＲ２用の光音響波検出信号を得る。この場合は、例えば第一の処理済み光音響波検出信号に重み付け係数０．７を与え、第二の処理済み光音響波検出信号に重み付け係数０．３を与えて重み付け加算を行い、境界領域ＢＲ１用の光音響波検出信号を得る。また、例えば第一の処理済み光音響波検出信号に重み付け係数０．３を与え、第二の処理済み光音響波検出信号に重み付け係数０．７を与えて重み付け加算を行い、境界領域ＢＲ２用の光音響波検出信号を得る。

以上の通りにして得た境界領域ＢＲ１用の光音響波検出信号により境界領域ＢＲ１の光音響画像を生成し、境界領域ＢＲ２用の光音響波検出信号により境界領域ＢＲ２の光音響画像を生成し、またアーチファクト非発生領域Ｒ１の境界領域ＢＲ１以外の部分については上記第一の処理済み光音響波検出信号により光音響画像を生成し、そしてアーチファクト発生領域Ｒ２の境界領域ＢＲ２以外の部分については上記第二の処理済み光音響波検出信号により光音響画像を生成し、それらの光音響画像を合成して１枚の光音響画像を生成し表示する。

そのようにして表示部１４に表示される光音響画像は、アーチファクト非発生領域Ｒ１、境界領域ＢＲ１、境界領域ＢＲ２、およびアーチファクト発生領域Ｒ２の各境界において周波数特性が少しずつ変化するようになり、光音響画像の診断性能が高められる。

また本実施形態は、周波数特性の急激な変化を抑制するために、前述したようにアーチファクト非発生領域Ｒ１、アーチファクト発生領域Ｒ２、境界領域ＢＲ１および境界領域ＢＲ２の各光音響波検出信号に対するフィルタリング処理の特性を変える場合と比べれば、計算処理が軽くなる点で優れている。

以上、光音響画像と共に反射超音波画像も生成かつ表示可能とされた光音響計測装置１０について説明したが、本発明は、反射超音波画像の生成および表示は行わず、光音響画像の生成および表示のみを行うように構成された光音響計測装置に適用することも勿論可能である。図１２は、そのように構成された光音響計測装置１０の一例を示すものである。この図１２に示す光音響計測装置１０は、図１に示したものと比べると、データ分離手段２３、超音波画像生成部２９および送信制御回路３３が除かれた形のものとなっている。

１０ 光音響計測装置
１１ プローブ
１２ 超音波ユニット
１３ レーザユニット
１４ 表示部
１５ 入力部
２０ 振動子アレイ
２１ 受信回路
２２ 受信メモリ
２３ データ分離手段
２４ 光音響画像生成部
２９ 超音波画像生成部
３０ 表示制御部
３３ 送信制御回路
３４ 制御部
４０ 光出射部
５０ 筺体
６０ 光ファイバ
６５ 吸収体
Ｂ アーチファクト非発生領域とアーチファクト発生領域との境界
Ｌ レーザ光（測定光）
Ｍ 被検体
Ｒ１ アーチファクト非発生領域
Ｒ２ アーチファクト発生領域
Ｕ 音響波
Ｕ１、Ｕ２ 光音響波

Claims (20)

1. 被検体に測定光を出射する光出射部と、該光出射部に並設され、前記測定光の出射によって前記被検体内に生じた光音響波を検出する音響波検出部とを有するプローブ、
   前記音響波検出部が出力した光音響波検出信号に基づいて光音響画像を生成する画像生成部、
   少なくとも前記光出射部と前記音響波検出部との位置関係に基づいて、前記光音響画像中のアーチファクト発生領域とアーチファクト非発生領域とを判別する領域判別手段、および、
   前記光音響波検出信号のうち、前記領域判別手段が判別したアーチファクト非発生領域の光音響画像に対応する第一の光音響波検出信号に対して第一のフィルタリング処理を行い、かつ、前記光音響波検出信号のうち、前記領域判別手段が判別したアーチファクト発生領域の光音響画像に対応する第二の光音響波検出信号に対して第二のフィルタリング処理を行うフィルタ手段を有し、
   前記第二のフィルタリング処理は、前記第一のフィルタリング処理に比べて、所定の周波数よりも低い周波数範囲の光音響波検出信号をより多く低減することを特徴とする光音響計測装置。
2. 前記フィルタ手段がバンドパスフィルタリング処理を行うバンドパスフィルタを含み、
   前記第二のフィルタリング処理としての前記バンドパスフィルタリング処理が、前記第一のフィルタリング処理としての前記バンドパスフィルタリング処理よりも、低周波数側のカットオフ周波数をより高周波数側に設定したものである請求項１記載の光音響計測装置。
3. 前記フィルタ手段がハイパスフィルタリング処理を行うハイパスフィルタを含み、
   前記第二のフィルタリング処理としての前記ハイパスフィルタリング処理が、前記第一のフィルタリング処理としての前記ハイパスフィルタリング処理よりも、カットオフ周波数をより高周波数側に設定したものである請求項１記載の光音響計測装置。
4. 前記領域判別手段が、前記光出射部と前記音響波検出部との位置関係と、前記アーチファクト発生領域およびアーチファクト非発生領域の境界とを対応させた対応表に基づいて、前記アーチファクト発生領域とアーチファクト非発生領域とを判別する請求項１から３いずれか１項記載の光音響計測装置。
5. 前記領域判別手段が、前記光出射部と前記音響波検出部との位置関係から、前記アーチファクト発生領域およびアーチファクト非発生領域の境界を算出する計算式により、前記アーチファクト発生領域とアーチファクト非発生領域とを判別する請求項１から３いずれか１項記載の光音響計測装置。
6. 前記領域判別手段が、前記光音響波が伝搬する被検体の部位における音速に応じて前記判別を補正する請求項１から５いずれか１項記載の光音響計測装置。
7. 前記フィルタ手段が、前記アーチファクト非発生領域およびアーチファクト発生領域の境界を含む範囲内に設定された少なくとも１つの境界領域に関する光音響波検出信号に対して、前記第一のフィルタリング処理および前記第二のフィルタリング処理とは異なる境界領域用フィルタリング処理を行う請求項１から６いずれか１項記載の光音響計測装置。
8. 前記境界領域用フィルタリング処理が、前記第一のフィルタリング処理と前記第二のフィルタリング処理との中間的なフィルタリング処理である請求項７記載の光音響計測装置。
9. 前記境界領域用フィルタリング処理の通過特性が、前記第一のフィルタリング処理の通過特性と、前記第二のフィルタリング処理の通過特性とに基づいて、前記境界領域の被検体深さ方向位置に応じて定められたものである請求項８記載の光音響計測装置。
10. 前記アーチファクト非発生領域およびアーチファクト発生領域の境界を含む範囲内に設定された少なくとも１つの境界領域に関する画像を、前記第一のフィルタリング処理を受けた前記第一の光音響波検出信号と、前記第二のフィルタリング処理を受けた前記第二の光音響波検出信号とを、前記境界領域の被検体深さ方向位置に応じて重み付け加算して得られた信号によって構成する手段をさらに備えた請求項１から６いずれか１項記載の光音響計測装置。
11. 被検体に測定光を出射する光出射部と、該光出射部に並設され、前記測定光の出射によって前記被検体内に生じた光音響波を検出する音響波検出部とを有するプローブ、および、
    前記音響波検出部が出力した光音響波検出信号に基づいて光音響画像を生成する画像生成部を備えた光音響計測装置において、
    少なくとも前記光出射部と前記音響波検出部との位置関係に基づいて、前記光音響画像中のアーチファクト発生領域とアーチファクト非発生領域とを判別し、
    前記光音響波検出信号のうち、前記アーチファクト非発生領域の光音響画像に対応する第一の光音響波検出信号に対して第一のフィルタリング処理を行い、
    前記光音響波検出信号のうち、前記アーチファクト発生領域の光音響画像に対応する第二の光音響波検出信号に対して第二のフィルタリング処理を行い、
    前記第二のフィルタリング処理を、前記第一のフィルタリング処理に比べて、所定の周波数よりも低い周波数範囲の光音響波検出信号をより多く低減するものとすることを特徴とする光音響計測装置の信号処理方法。
12. 前記フィルタリング処理がバンドパスフィルタリング処理であり、
    前記第二のフィルタリング処理としての前記バンドパスフィルタリング処理が、前記第一のフィルタリング処理としての前記バンドパスフィルタリング処理よりも、低周波数側のカットオフ周波数をより高周波数側に設定したものである請求項１１記載の光音響計測装置の信号処理方法。
13. 前記フィルタリング処理がハイパスフィルタリング処理であり、
    前記第二のフィルタリング処理としての前記ハイパスフィルタリング処理が、前記第一のフィルタリング処理としての前記ハイパスフィルタリング処理よりも、カットオフ周波数をより高周波数側に設定したものである請求項１１記載の光音響計測装置の信号処理方法。
14. 前記光出射部と前記音響波検出部との位置関係と、前記アーチファクト発生領域およびアーチファクト非発生領域の境界とを対応させた対応表に基づいて、前記アーチファクト発生領域とアーチファクト非発生領域とを判別する請求項１１から１３いずれか１項記載の光音響計測装置の信号処理方法。
15. 前記光出射部と前記音響波検出部との位置関係から、前記アーチファクト発生領域およびアーチファクト非発生領域の境界を算出する計算式により、前記アーチファクト発生領域とアーチファクト非発生領域とを判別する請求項１１から１３いずれか１項記載の光音響計測装置の信号処理方法。
16. 前記光音響波が伝搬する被検体の部位における音速に応じて前記判別を補正する請求項１１から１５いずれか１項記載の光音響計測装置の信号処理方法。
17. 前記アーチファクト非発生領域およびアーチファクト発生領域の境界を含む範囲内に設定された少なくとも１つの境界領域に関する光音響波検出信号に対して、前記第一のフィルタリング処理および前記第二のフィルタリング処理とは異なる境界領域用フィルタリング処理を行う請求項１１から１６いずれか１項記載の光音響計測装置の信号処理方法。
18. 前記境界領域用フィルタリング処理が、前記第一のフィルタリング処理と前記第二のフィルタリング処理との中間的なフィルタリング処理である請求項１７記載の光音響計測装置の信号処理方法。
19. 前記境界領域用フィルタリング処理の通過特性が、前記第一のフィルタリング処理の通過特性と、前記第二のフィルタリング処理の通過特性とに基づいて、前記境界領域の被検体深さ方向位置に応じて定められたものである請求項１８記載の光音響計測装置の信号処理方法。
20. 前記アーチファクト非発生領域およびアーチファクト発生領域の境界を含む範囲内に設定された少なくとも１つの境界領域に関する画像を、前記第一のフィルタリング処理を受けた前記第一の光音響波検出信号と、前記第二のフィルタリング処理を受けた前記第二の光音響波検出信号とを、前記境界領域の被検体深さ方向位置に応じて重み付け加算して得られた信号によって構成する請求項１１から１６いずれか１項記載の光音響計測装置の信号処理方法。

Images