(57)【特許請求の範囲】
【請求項１】
測定対象から伝播した弾性波が受信機により受信されて得られた圧力信号の初期値を取得する信号初期化手段と、
前記受信機のインパルス応答の初期値を取得するインパルス応答初期化手段と、
前記圧力信号の初期値および前記インパルス応答の初期値に基づいて修正圧力信号を出力する出力手段と、
を有する情報処理装置であって、
前記出力手段は、
インパルス応答を用いて圧力信号を修正する信号修正手段と、
高コントラスト性、正値性、および周波数帯域特性を拘束条件として圧力信号を修正する信号拘束手段と、
圧力信号を用いてインパルス応答を修正するインパルス応答修正手段と、
インパルス応答の時間帯域特性を拘束条件としてインパルス応答を修正するインパルス応答拘束手段とを有し、
前記信号修正手段、前記信号拘束手段、前記インパルス応答修正手段および前記インパルス応答拘束手段の処理を、所定の条件に従って反復して実行することにより、前記修正圧力信号を生成する
ことを特徴とする情報処理装置。
【請求項２】
前記インパルス応答初期化手段は、前記受信機のインパルス応答の実測または近似値の推算により、前記インパルス応答の初期値を取得する
ことを特徴とする請求項１に記載の情報処理装置。
【請求項３】
ｘｋはｋ回目の反復における圧力信号の推定値、γは緩和係数、ｈはインパルス応答、ｙは実際に測定された圧力信号、＊は複素共役、積記号は畳みこみ積であるとしたときに、
前記信号修正手段は、式（１）により圧力信号を修正する
ことを特徴とする請求項１または２に記載の情報処理装置。
【数１】
0006053265.tif 000004
【請求項４】
前記式（１）における緩和係数は、０より大きく２より小さい実数を前記インパルス応答のパワースペクトルの最大値で除したものを用いる
ことを特徴とする請求項３に記載の情報処理装置。
【請求項５】
ｈｋはｋ回目の反復におけるインパルス応答、γは緩和係数、ｘは圧力信号、ｙは実際に測定された圧力信号、＊は複素共役、積記号は畳みこみ積であるとしたときに、
前記インパルス応答修正手段は、式（２）によりインパルス応答を修正する
ことを特徴とする請求項１ないし４のいずれか１項に記載の情報処理装置。
【数２】
0006053265.tif 000005
【請求項６】
前記式（２）における緩和係数は、０より大きく２より小さい実数を、前記式（２）における圧力信号のパワースペクトルの最大値で除したものを用いる
ことを特徴とする請求項５に記載の情報処理装置。
【請求項７】
前記信号拘束手段は、前記高コントラスト性に対応する基線除去、前記正値性に対応する負値除去、前記周波数帯域特性に対応する高周波成分除去の順に圧力信号の修正を行う
ことを特徴とする請求項１ないし６のいずれか１項に記載の情報処理装置。
【請求項８】
前記信号拘束手段は、圧力信号の平均値、中央値および低周波成分のいずれかに所定の係数を乗じたのち当該圧力信号から減ずることにより、前記基線除去を行う
ことを特徴とする請求項７に記載の情報処理装置。
【請求項９】
前記基線除去における前記所定の係数は、０以上の実数である
ことを特徴とする請求項８に記載の情報処理装置。
【請求項１０】
前記信号拘束手段は、所定のカットオフ周波数を持つローパスフィルタにより、前記高周波成分除去を行う
ことを特徴とする請求項７ないし９のいずれか１項に記載の情報処理装置。
【請求項１１】
前記高周波成分除去における前記所定のカットオフ周波数は、前記受信機の帯域中心の周波数の１倍以上１０倍以下の値である
ことを特徴とする請求項１０に記載の情報処理装置。
【請求項１２】
前記インパルス応答拘束手段は、所定のカットオフ時間を持つ減衰関数をインパルス応答に乗じることにより、当該インパルス応答を修正する
ことを特徴とする請求項１ないし１１のいずれか１項に記載の情報処理装置。
【請求項１３】
前記インパルス応答拘束手段における前記所定のカットオフ時間は、前記受信機の帯域中心の周波数と、前記受信機の帯域幅を表すＱ値とから算出される
ことを特徴とする請求項１２に記載の情報処理装置。
【請求項１４】
前記出力手段は、
前記信号修正手段、前記信号拘束手段の順に圧力信号の修正を繰り返す第一のループ処理と、
前記インパルス応答修正手段、前記インパルス応答拘束手段の順にインパルス応答の修正を繰り返す第二のループ処理と、
前記第一のループ処理および前記第二のループ処理を繰り返す第三のループ処理の、回数及び順番を制御する
ことを特徴とする請求項１ないし１３のいずれか１項に記載の情報処理装置。
【請求項１５】
前記出力手段は、前記第一のループ処理、前記第二のループ処理および前記第三のループ処理を、それぞれ所定の回数だけ実行する
ことを特徴とする請求項１４に記載の情報処理装置。
【請求項１６】
前記出力手段は、前記第一のループ処理における圧力信号の修正量が所定の基準値を下回るまで前記第一のループ処理を実行する
ことを特徴とする請求項１４に記載の情報処理装置。
【請求項１７】
前記出力手段は、前記第二のループ処理におけるインパルス応答の修正量が所定の基準値を下回るまで前記第二のループ処理を実行する
ことを特徴とする請求項１４に記載の情報処理装置。
【請求項１８】
請求項１ないし１７のいずれか１項に記載の情報処理装置と、
当該情報処理装置によって修正された圧力信号を用いて、バックプロジェクション法、ｋ−ｓｐａｃｅ法、最小二乗法、制約付き最小二乗法および非線形最適化法のいずれかの手法により前記測定対象内部の画像データを作成する作成手段と、
を有することを特徴とする光音響装置。
【請求項１９】
測定対象から伝播した弾性波が受信機により受信されて得られた圧力信号の初期値を取得する信号初期化ステップと、
前記受信機のインパルス応答の初期値を取得するインパルス応答初期化ステップと、
前記圧力信号の初期値および前記インパルス応答の初期値に基づいて修正圧力信号を出力するステップと、
を有する情報処理方法であって、
前記出力するステップでは、
インパルス応答を用いて圧力信号を修正する信号修正ステップと、
高コントラスト性、正値性、および周波数帯域特性を拘束条件として圧力信号を修正する信号拘束ステップと、
圧力信号を用いてインパルス応答を修正するインパルス応答修正ステップと、
インパルス応答の時間帯域特性を拘束条件としてインパルス応答を修正するインパルス応答拘束ステップとが、所定の条件に従って反復して実行されることにより、前記修正圧力信号が生成される
ことを特徴とする情報処理方法。
【請求項２０】
請求項１９に記載の情報処理方法の各ステップを情報処理装置に実行させるプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、情報処理装置および情報処理方法に関する。
【背景技術】
【０００２】
近年、光音響イメージングの技術が開発されている。これは、測定対象に光を照射したときに光音響効果により発生する音響波（光音響波）を、受信素子を有する探触子などの受信機で受信し、音響波の圧力信号を求める技術である。これにより測定対象内部の三次元圧力分布や光学特性値分布を求めること、すなわち再構成が可能になる。よって、生体を測定対象とすれば、生体内部を画像化して腫瘍の観察など医療分野への応用ができる。
【０００３】
光音響イメージングで測定される圧力信号は、受信機のインパルス応答の影響を受けて劣化する。そのため再構成で得られる三次元圧力分布画像の解像度及びコントラストが低下する。上記解像度及びコントラストを改善するには、受信機のインパルス応答の影響を除去する必要がある。かかる除去操作は、数学的には逆畳み込み演算と等価であり、「デコンボリューション」と呼ばれる。
【０００４】
光音響イメージングにおけるデコンボリューション手法として非特許文献１の例がある。非特許文献１では、受信機のインパルス応答を実測して圧力信号のデコンボリューションを行う手法が開示されている。具体的には、励起光を微小領域に集中させ、その領域から発せられる圧力信号をインパルス応答の近似値として用いる手法である。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】Wang et al., “Photoacoustic imaging with deconvolution algorithm,” Phys. Med. Biol. 49 (2004) 3117-24
【発明の概要】
【発明が解決しようとする課題】
【０００６】
非特許文献１の手法では、インパルス応答の実測が必要であった。一般に、受信機のインパルス応答を正確に測定するには、十分な高周波成分を含む既知の入力信号を生成し、これに対する受信機の応答を正確に測定する必要がある。しかし、ノイズ存在下でかかる測定を行うのは困難である。
また、実際の応用の場面では、測定環境の違いや機器の経時変化に伴い、インパルス応答が変化する。上記の手法では、これらの変化に応じて何度もインパルス応答を測定する必要が生じるため、現実的な運用は困難であった。
【０００７】
本発明は、上記の問題点に鑑みてなされたものであり、光音響イメージングの際に、インパルス応答の影響を除去するデコンボリューション手法を好適に行うための技術を提供することを目的とする。
【課題を解決するための手段】
【０００８】
上記目的を解決するために、本発明は以下の構成を採用する。
すなわち、測定対象から伝播した弾性波が受信機により受信されて得られた圧力信号の初期値を取得する信号初期化手段と、前記受信機のインパルス応答の初期値を取得するインパルス応答初期化手段と、前記圧力信号の初期値および前記インパルス応答の初期値に基づいて修正圧力信号を出力する出力手段と、有する情報処理装置であって、前記出力手段は、インパルス応答を用いて圧力信号を修正する信号修正手段と、高コントラスト性、正値性、および周波数帯域特性を拘束条件として圧力信号を修正する信号拘束手段と、圧力信号を用いてインパルス応答を修正するインパルス応答修正手段と、インパルス応答の時間帯域特性を拘束条件としてインパルス応答を修正するインパルス応答拘束手段とを有し、前記信号修正手段、前記信号拘束手段、前記インパルス応答修正手段および前記インパルス応答拘束手段の処理を、所定の条件に従って反復して実行することにより、前記修正圧力信号を生成することを特徴とする情報処理装置である。
本発明はまた、以下の構成を採用する。
すなわち、測定対象から伝播した弾性波が受信機により受信されて得られた圧力信号の
初期値を取得する信号初期化ステップと、前記受信機のインパルス応答の初期値を取得するインパルス応答初期化ステップと、前記圧力信号の初期値および前記インパルス応答の初期値に基づいて修正圧力信号を出力するステップと、を有する情報処理方法であって、前記出力するステップでは、インパルス応答を用いて圧力信号を修正する信号修正ステップと、高コントラスト性、正値性、および周波数帯域特性を拘束条件として圧力信号を修正する信号拘束ステップと、圧力信号を用いてインパルス応答を修正するインパルス応答修正ステップと、インパルス応答の時間帯域特性を拘束条件としてインパルス応答を修正するインパルス応答拘束ステップとが、所定の条件に従って反復して実行されることにより、前記修正圧力信号が生成されることを特徴とする情報処理方法である。
【発明の効果】
【００１１】
本発明によれば、光音響イメージングの際に、インパルス応答の影響を除去するデコンボリューション手法を好適に行うための技術を提供できる。
【図面の簡単な説明】
【００１２】
【図１】本発明を適用した光音響イメージング装置の模式図。
【図２】本発明の実施の形態における処理の各工程を示したフローチャート。
【図３】本発明の適用前後で圧力信号を比較した図。
【図４】本発明の適用前後で再構成結果を比較した図。
【図５】本発明の前処理手段の構成を示した模式図。
【発明を実施するための形態】
【００１３】
以下、図面を参照しながら、本発明の実施形態について詳細に説明する。
本発明の技術は、光音響イメージングにおいて、測定対象内部の三次元圧力分布あるいは光学特性値分布を求める（再構成する）前の、圧力信号を整形する段階に適用される。従って本発明は、光音響イメージング装置としてはもちろん、光音響イメージング装置に組み込まれたり、接続されたりして使用される画像処理装置あるいは情報処理装置として実施することもできる。本発明はまた、コンピュータ等の画像処理装置や情報処理装置に
本発明の各工程における処理を実行させる画像処理プログラムまたは情報処理プログラムとしても実現できる。本発明はまた、プログラムに従い画像処理装置や情報処理装置が行う画像処理方法や情報処理方法としても実現できる。
【００１４】
光音響イメージング装置においては、光照射によって生じた弾性波の圧力信号の発生源分布（三次元圧力分布）、あるいはそれから導かれる吸収係数分布などの光学特性値分布が求められる。さらに光学特性値分布から、酸素飽和度分布や酸化・還元ヘモグロビン濃度分布のような物質濃度分布も求められる。さらにこれら各種の値の分布に基づき、例えば医療診断の目的で、表示手段に表示するための測定対象内部の画像データを作成し、画像化（再構成）を行うことができる。
【００１５】
＜光音響イメージング装置の構成と機能＞
図１は、本発明の実施の形態を含む光音響イメージング装置全体の模式図である。以下、装置の各構成要素の機能を説明する。
【００１６】
光源１はレーザ光源であり、測定対象２にパルス状のレーザ光１１を入射させる。測定対象が生体の場合、光源からは生体を構成する成分のうち特定の成分に吸収される波長の光を照射する。光源としては数ナノから数百ナノ秒オーダーのパルス光を発生させ得るパルス光源が好ましい。光源としてはレーザが好ましいが、レーザのかわりに発光ダイオードなどを用いることも可能である。レーザとしては、固体レーザ、ガスレーザ、色素レーザ、半導体レーザなど様々なレーザを使用することができる。光源が使用する波長に関しては、生体内において吸収が少ない７００ｎｍから１１００ｎｍの波長領域が好ましい。ただし、比較的生体表面付近の生体組織の光学特性値分布を求める場合は、上記の波長領域よりも範囲の広い、例えば４００ｎｍから１６００ｎｍの波長領域を使用することも可能である。
【００１７】
測定対象２は内部に吸光係数分布を持つ物体であり、光源１より入射した光を物体内の各点における吸光係数に応じて吸収する。吸収された光のエネルギーが物体内の各点のグリュナイゼン係数に応じて弾性波を生じさせることにより、測定対象内部に圧力分布を生じさせる。この圧力分布を本明細書では三次元圧力分布と呼ぶ。
【００１８】
測定対象２の内部に発生した三次元圧力分布は弾性波２１として周囲に伝播する。超音波受信機３に達した前記弾性波は、時系列の圧力データとして測定及び記録される。これを本明細書では圧力信号と呼称する。この弾性波は、音波、超音波、音響波、光音響波、光超音波とも呼ばれる。
超音波受信機は、弾性波を受信して電気信号に変換する検出器である。電気信号には、必要に応じてＡＤ変換や増幅処理が施される。測定対象から主に発生する１００ＫＨｚから１００ＭＨｚの弾性波を受信できる検出器が好適に用いられる。超音波受信機は、圧電現象を用いるもの、光の共振を用いるもの、容量の変化を用いるものなど、弾性波を受信できるものであればどのような検出器でもよい。また、複数の受信素子がリニア状あるいは平面状に配列された超音波受信機を用いれば、検出時間を短縮できると共に、測定対象の振動の影響を低減できる。超音波受信機は、本発明の受信機に相当する。
【００１９】
前記圧力信号に対して前処理手段４にてノイズ除去やデコンボリューション等の前処理工程を行う。本発明はこの前処理工程を特徴づけるものである。
前処理の済んだ圧力信号に対して再構成手段５にて再構成工程を行う。再構成とは圧力信号から前記三次元圧力分布を復元推定する工程である。再構成手段は、本発明の作成手段に相当する。
再構成された三次元圧力分布に対して、後処理手段６にて特徴抽出やコントラスト調整等の後処理工程を行う。
これら前処理手段、再構成手段、後処理手段は、例えばコンピュータ等の情報処理装置で実行されるプログラムのモジュールとして実現することもできるし、専用の回路を用いても構わない。
【００２０】
後処理の結果として得られた三次元圧力分布は、表示手段７に表示される。これらの工程により測定対象２内部の吸光係数分布に関する情報を三次元的に可視化できる。
以上が、本発明の実施の形態に係る光音響イメージング装置及びその動作の概要である。
【００２１】
＜デコンボリューションの各工程＞
次に、図２を参照しながら、本発明の実施の形態における各工程を説明する。なお、以下の説明では圧力信号の推定値及びインパルス応答の推定値を単に信号、インパルス応答と記し、特に区別が必要な場合以外、推定値である旨を書かない。各工程は、図５に示すように、前処理手段４の各構成要素によって行われる。
【００２２】
はじめに、ステップＳ２１において、信号初期化手段２１にて受信機より圧力信号の実測値を取得する（後述の式（１）、式（２）におけるｙ）。また、信号の初期値（式（１）におけるｘ０）を０とする。
【００２３】
次に、ステップＳ２２において、インパルス応答初期化手段２２にてインパルス応答の初期値（式（２）におけるｈ０）を算出もしくは取得する。算出する場合は、例えば、周波数特性において受信機と等価な減衰振動子モデルを用いてインパルス応答の近似値を推算してもよい。取得する場合は、例えば、ノイズを含むインパルス応答の実測値を近似値として用いてもよい。要するに真のインパルス応答を定性的に近似する時系列データを用いればよいのであって、例えば、機器の経時変化等でインパルス応答が若干変化しても、初期値を変化させる必要はない。
【００２４】
次に、ステップＳ２３において、反復制御手段２３にて信号のデコンボリューションを行う反復計算を開始する。図２では、前記反復計算を二重ループとして示してある。二重ループの中には、まず、後述の信号修正手段と信号拘束手段とによる処理をこの順で適当な回数反復して適用する内部ループ（第一のループ処理）がある。さらに、インパルス応答修正手段とインパルス応答拘束手段とによる処理をこの順で適当な回数反復して適用する内部ループ（第二のループ処理）がある。そして、これらの内部ループを順に繰り返し行う工程を外部ループ（第三のループ処理）としてある。以下、この二重ループによる処理を説明する。
しかし、これはあくまで実施形態の一例を示したに過ぎず、各修正手段および各拘束手段による処理の本質を損なわない限りにおいて、適宜、反復の順番や反復の回数を変更することができる。
【００２５】
次に、ステップＳ２４において、信号修正手段２４にて、この時点でのインパルス応答を真値とみなして信号の修正を行う。修正の方法は、例えば、式（１）に示すＬａｎｄｗｅｂｅｒ反復により行うことができる。なお、式（１）において、ｘｋはｋ回目の反復における信号の推定値、γは緩和係数、ｈはインパルス応答、ｙは実際に測定された圧力信号、＊は複素共役、積記号（○に×）は畳みこみ積である。
【数１】
0006053265.tif 000002
【００２６】
式（１）における緩和係数γは、反復の収束する速さを調整する係数であって、０より
大きく２より小さい実数をインパルス応答のパワースペクトル（フーリエ変換の絶対値の２乗）の最大値で除した値を使うことができる。これは反復写像のスペクトル半径が１より小さいこと（収束条件）から定まる。さらに言えば、反復が安定して行える限り、緩和係数は大きい値であるほど好ましい。
【００２７】
次に、ステップＳ２５において、信号拘束手段２５にて、信号の正値性、高コントラスト性、周波数帯域特性を拘束条件として信号を再修正する。
具体的には、例えば、基線除去（高コントラスト性に対応）、負値除去（正値性に対応）、高周波成分除去（周波数帯域特性に対応）をこの順で行うことにより信号を拘束してもよい。
【００２８】
また、前記基線除去は、例えば、信号の平均値、中央値、低周波成分のいずれかに所定の係数を乗じたのち元の信号から減ずることにより、基線を除去してもよい。この処理は信号のバックグラウンドが反復計算中に増加するのを抑制し、以って信号の高コントラスト性を保証するためのものである。したがってその所定の係数は０以上の実数であればよい。
【００２９】
また、前記負値除去は、例えば、単純に信号の負の部分を０に置換することで実現できる。この処理は圧力信号の時間積分（速度ポテンシャル）が正値であるという物理的要請を保証するためのものである。
【００３０】
また、前記高周波成分除去は、例えば、適当なカットオフ周波数によって特徴づけられるローパスフィルタによって、信号の高周波成分を除去してもよい。この処理は信号に含まれる高周波数のノイズが反復計算中に増加するのを抑制するためのものであり、修正後の信号の帯域上限を定める。したがって、そのカットオフ周波数は、受信機の帯域中心を特徴づける周波数の１倍以上１０倍以下の値とするのが好ましい。
【００３１】
次に、ステップＳ２６において、反復制御手段２３が、Ｓ２３で説明したような反復制御を行う。圧力信号修正処理を終了する条件が満たされていれば、本フローを終了し、再構成処理に移行する。
次に、ステップＳ２７において、インパルス応答修正手段２７にて、この時点での信号を真値とみなしてインパルス応答の修正を行う。修正の方法は、例えば、式（２）に示すＬａｎｄｗｅｂｅｒ反復により行うことができる。なお、式（２）において、ｈｋはｋ回目の反復におけるインパルス応答、γは緩和係数、ｘは信号、ｙは実際に測定された圧力信号、＊は複素共役、積記号（○に×）は畳みこみ積である。
【数２】
0006053265.tif 000003
【００３２】
式（２）における緩和係数に関する条件は、前記ステップＳ２４における条件と同様に、反復の収束する速さを調整する係数であって、０より大きく２より小さい実数を信号のパワースペクトルの最大値で除したものを用いてもよい。
【００３３】
次に、ステップＳ２８において、インパルス応答拘束手段２８にて、インパルス応答の時間帯域特性を拘束条件としてインパルス応答を再修正する。
具体的には、適当なカットオフ時間により特徴づけられる減衰関数をインパルス応答に乗じることにより、インパルス応答を拘束してもよい。この処理はインパルス応答が有限の時間で減衰するという物理的要請を保証するためのものである。したがってそのカットオフ時間は、インパルス応答の典型的な減衰時間を示す指標であればよく、例えば、受信
機の帯域中心を特徴づける周波数と、受信機の帯域幅を特徴づけるＱ値とから算出してもよい。
【００３４】
最後に、ステップＳ２３に戻り、反復制御手段２３にて、反復計算を終了するか否かの判定を行う。判定は、例えば、所定の反復回数に達したことで終了としてもよいし、前記信号修正手段もしくは前記インパルス応答修正手段における各修正量が所定の基準値を下回ったことで終了としてもよいし、またそれらの条件を組み合わせて用いてもよい。終了しない場合は上記の反復計算を継続する。
【００３５】
本フローにおける、信号初期化手段２１、インパルス応答初期化手段２２、反復制御手段２３、信号修正手段２４、信号拘束手段２５、インパルス応答修正手段２７、インパルス応答拘束手段２８は、前処理手段４の一部を構成するものである。この様子を図５に示す。反復制御手段２３は、不図示の制御線により他の構成要素を制御している。それぞれの手段は例えば、前処理手段を実現するプログラムの一部として構成することも可能であり、専用の回路を用いても構わない。また、信号修正手段、信号拘束手段、インパルス応答修正手段およびインパルス応答拘束手段は、インパルス応答の初期値と圧力信号の初期値に基づいて、最終的に修正された圧力信号である修正圧力信号を生成し出力する出力手段と捉えることもできる。
【００３６】
以上の各工程により、信号とインパルス応答とが同時に拘束条件下で最も有り得る値に収束する。なお、信号は圧力信号の速度ポテンシャルに収束する。
【００３７】
＜実施例＞
以下、実施例として、本発明の情報処理装置により上記フローに従って算出したデコンボリューションの結果について述べる。
【００３８】
測定データ１は下記のように測定したものである。光源として、Ｔｉ：Ｓａレーザ（波長７９７ｎｍ）を用いた。測定対象として、ウレタンゲルのブロックに直径０．３ｍｍのゴムワイヤを５ｍｍ間隔で包埋したものを用いた。超音波受信機として、ＣＭＵＴ（Capacitive Micromachined Ultrasonic Transducer）を用いた。
【００３９】
測定データ２は下記のように測定したものである。光源として、Ｔｉ：Ｓａレーザ（波長７９７ｎｍ）を用いた。測定対象として、ウレタンゲルのブロックに直径２ｍｍのチューブ状光吸収体（ウレタンゲルに光を吸収するインクを加えたもの）を１０ｍｍ間隔で包埋したものを用いた。超音波受信機として、ピエゾ素子を用いた超音波トランスデューサを用いた。
【００４０】
再構成には、既知のバックプロジェクション法を用いた。なお、本発明は再構成以前の圧力信号の整形に係るものであるから、再構成手法の如何に依らずその効果が発揮される。例えばバックプロジェクション法の代わりに、ｋ−ｓｐａｃｅ法、最小二乗法、制約付き最小二乗法、非線形最適化法のいずれかの手法により再構成してもよい。
【００４１】
各設定値は次のとおりである。インパルス応答の初期値は、減衰振動子モデルから算出したものを用いた。信号修正手段とインパルス応答修正手段における緩和係数γは、「１．９」をそれぞれインパルス応答、信号のパワースペクトルの最大値で割った値を用いた。信号拘束手段における基線除去は、信号の中央値に１をかけて元の信号から減ずる方法を採った。信号拘束手段におけるカットオフ周波数は、受信機の帯域中心周波数の３倍とした。インパルス応答拘束手段におけるカットオフ時間は、６０μｓを用いた。反復回数は、信号修正及び拘束ループを６回、インパルス応答修正及び拘束ループを６回とし、これら一連の処理を順に３００回反復した。
上記の設定値を用いて、図２に示したフローチャートに従ってデコンボリューションを行った。
【００４２】
図３は、測定データ１に対してデコンボリューションを行った結果を圧力信号として表示したものである。図３において、縦軸は時間を、横軸は超音波受信機の受信素子の位置を表す。また、右側のスケールは色の濃さに応じて圧力（Ｐ）を表す。図３（ａ）に示すように、デコンボリューション前には主信号（矢印に「Ｍ」の文字で示す）の後の時刻に、複数の本来存在しないピークが見られる。一方、図３（ｂ）に示すように、デコンボリューション後にはそれらのピークは見られなかった。
【００４３】
図４は、測定データ２に対してデコンボリューションを行い、再構成した結果である。図４（ａ）に示すように、デコンボリューション前には解像度もコントラストも低い。一方、図４（ｂ）に示すように、デコンボリューション後には、解像度、コントラストともに向上した。
【００４４】
なお、上記計算に要した時間は１．４ＧＦｌｏｐｓの計算機環境において９秒間であった。
【００４５】
上記の結果から、本発明のデコンボリューション処理によって、光音響イメージングにおける解像度、コントラストが改善し、本発明の課題が解決されたことが示された。
【００４６】
以上、特定の実施形態を参照しながら、本発明について詳解してきた。しかしながら、本発明は上記特定の形態に限定されず、本発明の技術思想を逸脱しない範囲で実施形態の修正をすることができる。
【符号の説明】
【００４７】
３：超音波受信機、４：前処理手段、２１：信号初期化手段、２２：インパルス応答初期化手段、２３：反復制御手段、２４：信号修正手段、２５：信号拘束手段、２７：インパルス応答修正手段、２８：インパルス応答拘束手段
【図１】
0006053265.tif 000006
【図２】
0006053265.tif 000007
【図３】
0006053265.tif 000008
【図４】
0006053265.tif 000009
【図５】
0006053265.tif 000010