WO2013128846A1

WO2013128846A1 - 光音響画像生成装置及び方法

Info

Publication number: WO2013128846A1
Application number: PCT/JP2013/000929
Authority: WO
Inventors: 辻田　和宏
Original assignee: 富士フイルム株式会社
Priority date: 2012-02-28
Filing date: 2013-02-20
Publication date: 2013-09-06
Also published as: US9974440B2; JP2013176414A; JP5762995B2; US20140371571A1

Abstract

　光音響イメージングにおいて、所望の構造物を選択的に画像化する。光音響信号を検出し、検出した光音響信号を再構成して光音響画像データを生成する。光音響画像データを２次元以上の空間でフーリエ変換し、空間周波数領域の光音響画像データを生成する。空間周波数領域の光音響画像データから所定の空間周波数成分を抽出し、抽出された空間周波数成分をフーリエ逆変換して、空間周波数処理後の光音響画像データを生成する。

Description

光音響画像生成装置及び方法

　本発明は、光音響画像生成装置及び方法に関し、更に詳しくは、被検体に光を照射し、光照射により被検体内で生じた音響波を検出して光音響画像を生成する光音響画像生成装置及び方法に関する。

　生体内部の状態を非侵襲で検査できる画像検査法の一種として、超音波検査法が知られている。超音波検査では、超音波の送信及び受信が可能な超音波探触子を用いる。超音波探触子から被検体（生体）に超音波を送信させると、その超音波は生体内部を進んでいき、組織界面で反射する。超音波探触子でその反射音波を受信し、反射超音波が超音波探触子に戻ってくるまでの時間に基づいて距離を計算することで、内部の様子を画像化することができる。

　また、光音響効果を利用して生体の内部を画像化する光音響イメージングが知られている。一般に光音響イメージングでは、パルスレーザ光を生体内に照射する。生体内部では、生体組織がパルスレーザ光のエネルギーを吸収し、そのエネルギーによる断熱膨張により超音波（光音響信号）が発生する。この光音響信号を超音波プローブなどで検出し、検出信号に基づいて光音響画像を構成することで、光音響信号に基づく生体内の可視化が可能である。

　ここで、光音響イメージングに関し、空間周波数処理を行うことが特許文献１に記載されている。特許文献１では、光源からの照射光を被検体に照射し、複数の検出器素子で光音響信号を検出する。検出された光音響信号は、被検体の表面で発生した光音響信号と、被検体内部で発生した光音響信号とが含まれる。特許文献１では、同一受信時間における複数の検出器素子による検出信号を空間方向にフーリエ変換し、空間周波数信号を取得する。その空間周波数信号のうち所定の周波数以下の成分を低減し、低域成分が低減された空間周波数信号に対してフーリエ逆変換を行い、逆変換により得られた信号に基づいて光音響画像を生成する。被検体の表面で発生する光音響信号は、複数の検出器素子で同時に取得されることから、被検体表面で発生する光音響波の影響を軽減したイメージングが可能となる。

　また、空間周波数処理に関し、特許文献２には、超音波イメージングにおいて、イメージデータをフーリエ変換して空間周波数領域のイメージデータを生成し、フーリエ変換空間周波数領域のイメージデータの所定の低域成分を除去し、低域成分が除去された空間周波数領域のイメージデータを逆変換して実データ領域のイメージデータを生成することが記載されている。特許文献２には、低域成分はシェーディングに相当する部分であり、また、診断に必要なイメージ部分は低域に存在していないので、空間周波数の低い成分を除去することはシェーディングを除去することに相当すると記載されている。つまり、特許文献２では、イメージデータにおいて所望の見たい領域は空間周波数領域では高域に存在しているため、その必要な領域以外を除去するように空間周波数領域で低域除去を実行している。

特開２０１１－２１７７６７号公報特開平７－２８２２４７号公報

　ところで、光音響で画像化する部分として、例えば血管などが考えられる。血管には、太さの太いものから細いものまで種々のサイズの血管があり、その中から所望のサイズの血管を画像化できれば、画像診断に役立つと考えられる。しかしながら特許文献１では、光音響波の検出信号そのものに対し、同一時間の信号を空間方向にフーリエ変換しているため、所望のサイズの血管を抽出して画像化することはできない。また、特許文献２は、対象が超音波画像であると共に、単に見たい領域が空間周波数領域の高域に存在しているために、低域成分を除去する空間周波数処理を行っているだけであり、所望の構造物を抽出して画像化することはできない。

　本発明は、上記に鑑み、光音響イメージングにおいて、所望の構造物を選択的に画像化並びに強調することができる光音響画像生成装置及び方法を提供することを目的とする。

　上記目的を達成するために、本発明は、被検体に照射すべき光を出射する光源と、被検体に対する光照射により被検体内で生じた光音響波を検出する音響波検出手段と、光音響波の検出信号である光音響信号を再構成し画像データを生成する光音響画像生成手段と、画像データを２次元以上の空間でフーリエ変換し、空間周波数領域の画像データを生成するフーリエ変換手段と、空間周波数領域の画像データから所定の空間周波数成分を抽出する空間周波数処理手段と、抽出された空間周波数成分をフーリエ逆変換し、空間周波数処理後の画像データを生成するフーリエ逆変換手段とを備えたことを特徴とする光音響画像生成装置を提供する。

　空間周波数処理手段が、第１の空間周波数以上で、かつ、第１の空間周波数よりも高い第２の周波数領域以下の空間周波数成分を抽出することとしてもよい。

　上記の第１の空間周波数及び第２の空間周波数は、空間周波数領域（波数空間）における位置に応じて変化してもよい。

　本発明では、空間周波数処理手段が、ユーザが指定した観察対象条件に従って、抽出する所定の空間周波数成分を決定する構成とすることができる。

　３次元的に検出された光音響波に基づく３次元画像データから、３次元空間を構成する軸のうちの１つの軸に垂直な面に沿って断面を切り出し、該切り出した断面を含む１つの軸に沿った方向の所定範囲の画像データを統合した断面画像データを生成する断面画像生成手段を更に備え、フーリエ変換手段が断面画像データを２次元フーリエ変換することとしてもよい。

　断面画像生成手段が、１つの軸に沿った複数の位置で断面を切り出して複数の断面画像データを生成し、フーリエ変換手段が、生成された複数の断面画像データのそれぞれをフーリエ変換して複数の空間周波数領域の画像データを生成し、空間周波数処理手段が、生成された複数の空間周波数領域の画像データのそれぞれから所定の空間周波数成分を抽出し、フーリエ逆変換手段が、複数の空間周波数領域の画像データのそれぞれから抽出された所定の空間周波数成分をフーリエ逆変換し、複数の断面画像データのそれぞれについて空間周波数処理後の断面画像データを生成するようにしてもよい。

　断面画像生成手段は、断面を等間隔で切り出してもよい。

　断面画像生成手段は、被検体に深さ方向に対応する軸に沿って、超音響波検出手段の音響波検出面と平行な方向に断面を切り出してもよい。

　断面画像生成手段が、所定範囲内の画像データの最大値を投影し、又は所定範囲内の画像データを積分することで、所定範囲内の画像データを統合することとしてもよい。

　空間周波数領域の断面画像データから抽出された所定の空間周波数成分をフーリエ逆変換することで得られた空間周波数処理後の断面画像データを表示装置に表示させる画像表示制御手段を更に備える構成としてもよい。

　ユーザ操作に応じて、表示される断面画像データの断面位置を決定する断面位置選択手段を更に備え、表示制御手段が、ユーザが指定した断面位置に対応する空間周波数処理された断面画像データを表示装置に表示させてもよい。

　表示制御手段が、断面位置を指定するためのスライドバーを表示装置に表示させ、断面位置選択手段が、ユーザのスライドバーの操作に応じて断面位置を決定するようにすることもできる。

　表示制御手段が、空間周波数処理前の断面画像データと、空間周波数処理後の断面画像データとを並べて表示装置に表示させることとしてもよい。

　表示制御手段が、空間周波数処理後の断面画像データを２値化して表示装置に表示してもよい。

　フーリエ変換手段が、３次元的に検出された光音響波に基づく３次元画像データを３次元フーリエ変換するようにしてもよい。

　光源が相互に異なる複数の波長の光を出射し、音響波検出手段が、複数の波長の光のそれぞれが照射された後に被検体内で発生した複数の光音響波を検出し、光音響画像生成手段が、複数の波長の光に対応して検出された光音響信号をそれぞれ再構成して画像データし、フーリエ変換手段が、複数の波長の光に対応する画像データのそれぞれをフーリエ変換し、空間周波数処理手段が、フーリエ変換された、複数の波長の光に対応する空間周波数領域の画像データから、それぞれ各波長に対応した所定の空間周波数成分を抽出し、フーリエ逆変換手段が、抽出された処理空間周波数成分のそれぞれをフーリエ逆変換してもよい。

　複数の波長の光が第１の波長の光と第２の波長の光とを含み、空間周波数処理手段が、第１の波長の光に対応する空間周波数領域の画像データから抽出する空間周波数成分と、第２の波長の光に対応する空間周波数領域の画像データから抽出する空間周波数とが相互に異なるようにしてもよい。

　光音響信号から、被検体に照射された光の微分波形をデコンボリューションするデコンボリューション手段を更に備える構成としてもよい。

　本発明は、また、光源からの光を被検体に照射するステップと、被検体に対する光照射により被検体内で生じた光音響波を検出するステップと、光音響波の検出信号である光音響信号を再構成して画像データを生成するステップと、画像データを２次元以上の空間でフーリエ変換し、空間周波数領域の画像データを生成するステップと、空間周波数領域の画像データから所定の空間周波数成分を抽出するステップと、抽出された処理空間周波数成分をフーリエ逆変換し、空間周波数処理後の画像データを生成するステップとを有することを特徴とする光音響画像生成方法を提供する。

　本発明の光音響画像生成装置及び方法は、光音響信号を再構成した画像データをフーリエ変換して空間周波数領域の画像データに変換し、その画像データから所定の空間周波数成分を抽出する。抽出した空間周波数成分をフーリエ逆変換し、空間周波数処理後の画像データを生成する。観察対象のサイズなどに応じて抽出する空間周波数成分を適宜選択することで、所望の構造物を選択的に画像化することができる。

本発明の第１実施形態の光音響画像生成装置を示すブロック図。光音響信号の検出空間を示す図。ＸＺ断面の断層画像を示す図。断面画像生成手段が生成する断面画像を示す図。断面画像データの画像例を示す図。空間周波数領域の断面画像データを示す図。空間周波数フィルタの一例を示す図。抽出された所定の空間周波数成分の画像データを示す図。フーリエ逆変換することで得られた断面画像データを示す図。光音響画像生成装置の動作手順を示すフローチャート。楕円形の空間周波数フィルタを示す図。低周波領域を除去する空間周波数フィルタを示す図。本発明の第２実施形態の光音響画像生成装置を示すブロック図。表示画面例を示す図。本発明の第３実施形態の光音響画像生成装置を示すブロック図。

　以下、図面を参照し、本発明の実施の形態を詳細に説明する。図１は、本発明の第１実施形態の光音響画像生成装置を示す。光音響画像生成装置（光音響画像診断装置）１０は、超音波探触子（プローブ）１１、超音波ユニット１２、及び光源（レーザユニット）１３を備える。

　レーザユニット１３は、光源であり、被検体に照射する光（レーザ光）を生成する。レーザ光の波長は、観察対象物に応じて適宜設定すればよい。レーザユニット１３は、例えばヘモグロビンの吸収が大きい波長、具体的には７５０ｎｍや８００ｎｍの波長の光を出射する。レーザユニット１３が出射するレーザ光は、例えば光ファイバなどの導光手段を用いてプローブ１１まで導光され、プローブ１１から被検体に照射される。あるいは、プローブ１１以外の場所から光照射を行うこととしてもよい。

　プローブ１１は、被検体内からの音響波（超音波）を検出する音響波検出手段を有する。音響波検出手段の音響波検出器素子は、プローブ１１は、例えば一次元的に配列された複数の超音波振動子を有している。プローブ１１は、被検体内の測定対象物がレーザユニット１３からの光を吸収することで生じた光音響波（以下、光音響信号とも呼ぶ）を複数の超音波振動子により検出する。

　超音波ユニット１２は、受信回路２１、ＡＤ変換手段２２、受信メモリ２３、光音響画像再構成手段２４、検波・対数変換手段２５、画像構築手段２６、断面画像生成手段２７、フーリエ変換手段２８、空間周波数処理手段２９、フーリエ逆変換手段３０、表示制御手段３１、トリガ制御回路３２、及び制御手段３３を有する。制御手段３３は、超音波ユニット１２内の各部を制御する。受信回路２１は、プローブ１１が検出した光音響波の検出信号（光音響信号）を受信する。ＡＤ変換手段２２は、受信回路２１が受信した光音響信号をサンプリングしてデジタル信号に変換する。ＡＤ変換手段２２は、例えばＡＤクロック信号に同期して、所定のサンプリング周期で光音響信号をサンプリングする。

　トリガ制御回路３２は、レーザユニット１３に対して光出射を指示する光トリガ信号を出力する。レーザユニット１３は、図示しないＹＡＧやチタン－サファイアなどのレーザ媒質を励起するフラッシュランプ４１と、レーザ発振を制御するＱスイッチ４２とを含む。レーザユニット１３は、トリガ制御回路３２がフラッシュランプトリガ信号を出力すると、フラッシュランプ４１を点灯し、レーザ媒質を励起する。トリガ制御回路３２は、例えばフラッシュランプ４１がレーザ媒質を十分に励起させると、Ｑスイッチトリガ信号を出力する。Ｑスイッチ４２は、Ｑスイッチトリガ信号を受けるとオンし、レーザユニット１３からレーザ光を出射させる。フラッシュランプ４１の点灯からレーザ媒質が十分な励起状態となるまでに要する時間は、レーザ媒質の特性などから見積もることができる。

　なお、トリガ制御回路３２からＱスイッチを制御するのに代えて、レーザユニット１３内において、レーザ媒質を十分に励起させた後にＱスイッチ４２をオンにしてもよい。その場合は、Ｑスイッチ４２をオンにした旨を示す信号を超音波ユニット１２側に通知してもよい。ここで、光トリガ信号とは、フラッシュランプトリガ信号とＱスイッチトリガ信号の少なくとも一方を含む概念である。トリガ制御回路３２からＱスイッチトリガ信号を出力する場合はＱスイッチトリガ信号が光トリガ信号に対応し、レーザユニット１３にてＱスイッチトリガのタイミングを生成する場合はフラッシュランプトリガ信号が光トリガ信号に対応していてもよい。光トリガ信号が出力されることで、被検体に対するレーザ光の照射及び光音響信号の検出が行われる。

　また、トリガ制御回路３２は、ＡＤ変換手段２２に対して、サンプリング開始を指示するサンプリングトリガ信号を出力する。トリガ制御回路３２は、光トリガ信号の出力後、所定のタイミングでサンプリングトリガ信号を出力する。トリガ制御回路３２は、光トリガ信号の出力後、好ましくは、被検体に実際にレーザ光が照射されるタイミングで、サンプリングトリガ信号を出力する。例えばトリガ制御回路３２は、Ｑスイッチトリガ信号の出力と同期してサンプリングトリガ信号を出力する。ＡＤ変換手段２２は、サンプリングトリガ信号を受けると、プローブ１１にて検出された光音響信号のサンプリングを開始する。

　ＡＤ変換手段２２は、サンプリングした光音響信号を、受信メモリ２３に格納する。受信メモリ２３には、例えば半導体記憶装置を用いることができる。あるいは、受信メモリ２３に、その他の記憶装置、例えば磁気記憶装置を用いてもよい。受信メモリ２３には、光音響信号のサンプリングデータ（光音響データ）が格納される。

　光音響画像再構成手段２４は、受信メモリ２３から光音響信号を読み出し、読み出した光音響信号を再構成する。光音響画像再構成手段２４は、光音響信号に基づいて、断層画像である光音響画像の各ラインの画像データを生成する。ここで、再構成された光音響信号は、光音響画像とみなすことができる。光音響画像再構成手段２４は、遅延加算法（Delay and Sum、位相整合加算、整相加算と同義）により、光音響信号を再構成する。光音響画像再構成手段２４は、例えば６４素子分の光音響信号を、各素子（各超音波振動子）の位置に応じた遅延時間で加算し、１ライン分のデータを生成する。その際、被検体内の音速は一定であると仮定してよいし、音速分布を考慮して各素子の遅延時間を補正してもよい。遅延加算方向に代えて、ハフ変換法又はフーリエ変換法を用いて再構成を行ってもよい。

　検波・対数変換手段２５は、光音響画像再構成手段２４が出力する各ラインのデータの包絡線を生成し、その包絡線を対数変換してダイナミックレンジを広げる。画像構築手段２６は、対数変換が施された各ラインのデータに基づいて、光音響画像を生成する。画像構築手段２６は、例えば光音響信号（ピーク部分）の時間軸方向の位置を、断層画像における深さ方向の位置に変換して光音響画像を生成する。光音響画像再構成手段２４、検波・対数変換手段２５、及び画像構築手段２６は、光音響画像生成手段に相当する。光音響画像生成手段は、例えば３次元的に検出された光音響波に基づく３次元画像データを生成する。

　断面画像生成手段２７は、光音響の３次元画像データから所定の面の断面を切り出す。断面画像生成手段２７は、例えば３次元空間を構成する軸のうちの１つの軸に垂直な面に沿って、画像データを切り出す。断面画像生成手段２７は、切り出した断面を含む、その断面に垂直な方向の所定範囲の画像データを統合した断面画像データを生成する。例えば断面画像生成手段２７が、被検体の深さ方向に相当する軸のある位置で、プローブ１１の音響波検出面に平行な面に沿って断面を切り出す場合、その断面位置の前後（浅い方向と深い方向）の所定枚数分の断面画像を、１つの画像に統合する。断面画像生成手段２７は、例えば所定範囲内の画像データの最大値を投影することで、所定範囲内の画像データを１つの画像データに統合する。あるいは、所定範囲内の画像データを積分（平均）することで、所定範囲内の画像データを統合してもよい。

　フーリエ変換手段２８は、光音響の画像データをフーリエ変換し、空間周波数領域の画像データを生成する。本実施形態では、フーリエ変換手段２８は、断面画像生成手段２７が生成した断面画像データを、２次元フーリエ変換する。空間周波数処理手段２９は、フーリエ変換された空間周波数領域の画像データから、所定の空間周波数成分を抽出する。空間周波数処理手段２９は、例えば空間周波数領域の画像データから、第１の空間周波数以上で、かつ、第１の空間周波数よりも高い第２の周波数領域以下の空間周波数成分を選択的に抽出する。

　フーリエ逆変換手段３０は、空間周波数処理手段２９で抽出された空間周波数成分をフーリエ逆変換し、空間周波数処理後の画像データを生成する。表示制御手段３１は、フーリエ逆変換手段３０が逆変換を行うことで得られた空間周波数処理後の断面画像データを、ディスプレイ装置などの画像表示手段１４の表示画面上に表示させる。表示制御手段３１は、例えば階調データで表わされている空間周波数処理後の断面画像データを、２値化して表示装置に表示させてもよい。

　図２は、光音響信号の検出空間を示す。光音響信号の時間方向は、光音響画像の深さ方向（Ｚ方向）に対応している。プローブ１１は、例えばＸ方向に一次元配列された複数の検出器素子（超音波振動子）を有している。そのようなプローブ１１をＹ方向に走査することにより、光音響信号を３次元的に取得できる。複数の検出器素子が一次元配列されたプローブを走査するのに代えて、複数の検出器素子がＸ方向及びＹ方向に２次元配列されたプローブを用いてもよい。この場合は、プローブを走査しなくても、光音響信号を３次元的に取得できる。

　図３は、ＸＺ断面の断層画像（光音響画像）を示す。例えばプローブ１１がＸ方向に配列された複数の超音波振動子を有する場合、光音響画像生成手段は、ＸＺ断面の光音響画像を生成する。例えば血管がＹ方向に横行しているとすると、ＸＺ断面の光音響画像には円形の血管断面が現れる。プローブ１１をＹ方向に走査し、各走査位置においてＸＺ断面の光音響画像を生成し、それら断面画像を複数枚つなげることで、３次元の光音響画像データが得られる。

　図４は、断面画像生成手段２７が生成する断面画像を示す。断面画像生成手段２７は、例えばＺ軸に沿って、プローブ１１の音響波検出面と平行な面内（ＸＹ平面）で、３次元の光音響像データから断面を切り出し、所定枚数分の画像データを統合した断面画像データを生成する。断面画像生成手段２７は、例えば深さ方向（Ｚ方向）２ｍｍ分の厚みに相当する複数枚の画像データの最大値を投影することで、１つの断面画像データを生成する。断面画像生成手段２７は、例えばＺ軸に沿って等間隔に並ぶ複数の位置で断面を切り出し、２ｍｍごとに断面画像データを生成する。断面画像生成手段２７が生成する断面画像はＸＹ平面に平行な断面画像には限られず、断面画像生成手段２７がＸＺ平面又はＹＺ平面に平行な断面画像を生成するようにしてもよい。

　フーリエ変換手段２８は、図４に示す複数の断面画像データのそれぞれに対してフーリエ変換を行い、複数の空間周波数領域の断面画像データを生成する。空間周波数処理手段２９は、生成された複数の空間周波数領域の断面画像データのそれぞれから所定の空間周波数成分を抽出する。フーリエ逆変換手段３０は、複数の空間周波数領域の断面画像データのそれぞれから抽出された所定の空間周波数成分をフーリエ逆変換する。逆変換を行うことで、図４に示す複数の断面画像データから、空間周波数処理後の断面画像データを生成できる。複数の空間周波数処理後の断面画像データをつなげることで、空間周波数処理後の３次元光音響画像データが得られる。

　図５は、断面画像データの画像例を示す。図５に示す画像は、図４に示す複数の断面画像データのうちの１つに対応する。図６は、フーリエ変換された空間周波数領域（波数空間）の断面画像データを示す。図５に示す断面画像データを空間周波数領域のデータに変換すると、図６に示すデータが得られる。図６において、横軸は波数ｋｘを表し、縦軸は波数ｋｙを表す。ｋｘとｋｙの交点は原点であり、原点に近いほど空間周波数成分が低い。

　図７は、空間周波数処理手段２９が適用する空間周波数フィルタの一例を示す。図７において、白色の部分は空間周波数フィルタで遮断される部分を表し、グレーの部分は空間周波数フィルタを透過する部分を表している。空間周波数処理手段２９は、例えば図７に示すような空間周波数フィルタを、図６に示す波数空間の画像データに適用する。その結果、空間周波数処理手段２９は、波数空間における原点から距離が一定距離以内の低周波成分と、波数空間における原点からの距離が別の一定距離以上の高周波成分とを除去する。

　図８は、空間周波数処理手段２９が抽出する所定の空間周波数成分の画像データを示す。図８に示すように、空間周波数処理手段２９は、図６の空間周波数領域に画像データにおける、図７にグレーで示す輪帯状の領域内のデータを、所定の空間周波数成分のデータとして抽出する。

　図９は、図８に示す空間周波数成分の画像データをフーリエ逆変換することで得られた断面画像データを示す。波数空間において所定の空間周波数成分のみを選択的に抽出し、その抽出されたデータをフーリエ逆変換することで、所定の空間周波数成分に対応する部分を画像化することができる。図５に示す空間周波数処理前の断面画像データと、図９に示す断面画像データとを比較すると、画像ノイズを抑制して、血管部分をより認識しやすくなっていることが確認できる。

　図７に示す空間周波数フィルタは、観察対象物の大きさなどに応じて適宜決定すればよい。例えばユーザ（オペレータ）は、観察対象条件として、抽出したい血管の太さを指定する。空間周波数処理手段２９は、指定された血管の太さに対応する空間周波数フィルタを設定する。一例として、抽出対象として０．５ｍｍの血管が指定されたとき、空間周波数処理手段２９は、画像上で１～２cycle/mmに相当する周波数帯域（輪帯状）のみを選択的に透過させるような空間周波数フィルタを設定する。あるいは、２cycle/mmを中心として、そこから所定の範囲（例えば±α）の範囲を選択的に透過させるような空間周波数フィルタを設定してもよい。抽出対象の血管の太さが１ｍｍであれば、０．３～０．８cycle/mmの範囲の周波数成分を選択的に透過させるような空間周波数フィルタを用いればよい。

　上記では、ユーザが抽出したい血管の太さを指定する例を説明したが、ユーザが血管の太さを指定するのに代えて、観察部位を指定するようにしてもよい。空間周波数処理手段２９には、観察部位に対応する空間周波数フィルタの情報があらかじめ記憶されている。空間周波数処理手段２９は、ユーザが観察部位を指定すると、その観察部位に対応する空間周波数フィルタの情報を読み出し、読み出した空間周波数フィルタの情報に従って空間周波数フィルタを設定する。

　続いて動作手順を説明する。図１０は光音響画像生成装置１０の動作手順を示す。トリガ制御回路３２は、フラッシュランプトリガ信号をレーザユニット１３に出力する。レーザユニット１３では、フラッシュランプトリガ信号に応答してフラッシュランプ４１が点灯し、レーザ媒質の励起が開始される。トリガ制御回路３２は、Ｑスイッチトリガ信号をレーザユニット１３に送り、Ｑスイッチ４２をオンさせることで、レーザユニット１３からパルスレーザ光を出射させる（ステップＳ１）。トリガ制御回路３２は、例えばフラッシュランプトリガ信号を出力するタイミングと所定の時間関係にあるタイミングでＱスイッチトリガ信号を出力する。例えばトリガ制御回路３２は、フラッシュランプ発光から１５０μｓ秒後に、Ｑスイッチトリガ信号を出力する。

　レーザユニット１３から出射したレーザ光は、被検体に照射される。被検体内では、照射されたパルスレーザ光による光音響信号が発生する。プローブ１１は、被検体内で発生した光音響信号を検出する（ステップＳ２）。プローブが検出した光音響信号は、受信回路２１を介してＡＤ変換手段２２に入力される。ＡＤ変換手段２２は、光音響信号をサンプリングしてデジタルデータに変換し、受信メモリ２３に格納する。例えば複数の超音波振動子が一次元配列されたプローブ１１を走査しつつ、複数の走査位置で光照射及び光音響信号の検出を行うことで、光音響信号の３次元データが得られる。

　光音響画像再構成手段２４は、受信メモリ２３から光音響信号を読み出し、読み出した光音響信号を再構成する。検波・対数変換手段２５は、再構成された光音響信号に対して検波・対数変換を行い、画像構築手段２６は、検波・対数変換後の光音響信号に基づいて光音響画像を生成する（ステップＳ３）。光音響画像再構成手段２４、検波・対数変換手段２５、及び画像構築手段２６は、光音響信号の３次元データから、光音響の３次元画像データを生成する。

　断面画像生成手段２７は、光音響の３次元画像データから、所定の枚数分だけ所定面内の断面を統合した断面画像データを生成する（ステップＳ４）。断面画像生成手段２７は、例えばプローブ１１の音響波検出面に平行な断面を、所定の厚み分ずつ統合した断面画像データを生成する。複数枚の画像を１つの断面画像データに統合することで、例えば音響波検出面に垂直な方向に血管の位置が変動する場合でも、１つの断面画像データ内に変動する血管を納めることができる。フーリエ変換手段２８は、ステップＳ４で生成された断面画像をフーリエ変換する（ステップＳ５）。空間周波数処理手段２９は、フーリエ変換された空間周波数領域の断面画像データから所定の周波数成分を抽出する（ステップＳ６）。フーリエ逆変換手段３０は、抽出された所定の周波数成分のデータをフーリエ逆変換する（ステップＳ７）。

　表示制御手段３１は、画像表示手段１４の表示画面上に、逆変換により得られた断面画像データを表示する（ステップＳ８）。表示制御手段３１は、断面画像生成手段２７が生成した複数の断面画像データに対応した、空間周波数処理後の複数の断面画像データを並べて表示してもよい。また、表示制御手段３１は、断面画像生成手段２７が生成した空間周波数処理前の断面画像データ（図５）と、空間周波数処理後の断面画像データ（図９）とを並べて表示させてもよい。並べて表示するのに代えて、空間周波数処理前の断面画像データと空間周波数処理後の断面画像データとを、表示色を変えて重ねて表示させてもよい。

　本実施形態では、光音響信号を再構成した画像データをフーリエ変換し、空間周波数領域の画像データから所定の空間周波数成分を抽出してフーリエ逆変換する。例えば図７に示すような空間周波数フィルタを用いることで、不要な低周波成分を除去できると共に、組織由来の高周波成分（画像ノイズ）を除去することができる。

　ここで、被検体の内部に存在する光吸収体から発生した光音響信号は、光吸収体の位置と検出器素子との間の距離に応じた時間差を持って各検出器素子で検出され、複数の検出器素子で同時には検出されない。特許文献１では、検出信号そのものを空間方向にフーリエ変換しており、空間周波数領域の検出信号から所定の低域成分を低減したとしても、ある太さの血管を抽出することはできない。これに対し、本願発明では、再構成後のイメージデータをフーリエ変換する。本実施形態では、実数法に直した波数空間でフィルタリングするので、血管径情報に対してフィルタを選定することができ、観察対象に応じて使用するフィルタを適宜設定することで、所望の太さの血管を画像化することができる。

　なお、図７では原点からの距離が一定範囲の周波数成分を透過するフィルタを用いる例を示したが、抽出する空間周波数成分はこれには限定されない。透過帯域の下限となる空間周波数と、上限となる空間周波数とが、波数空間における位置に応じて変化するようにしてもよい。図１１は、楕円形の空間周波数フィルタを示す。波数空間において強度の強い領域があれば、その領域を長軸方向とする楕円状の空間周波数フィルタを用いてもよい。図１１に示す空間周波数フィルタを用いた場合、ｋｙ方向に比して、ｋｘ方向において高周波の周波数成分が抽出される。

　また、空間周波数フィルタは、所定帯域の空間周波数成分を通過させるものには限られない。図１２は、低周波領域を除去する空間周波数フィルタを示す。図１２に示す空間周波数フィルタを用いた場合、ある空間周波数以上の空間周波数成分を抽出することができる。空間周波数フィルタのフィルタ特性は円形や楕円形に限られず、多角形でもよい。

　本実施形態において、レーザユニット１３から相互に異なる複数の波長の光を被検体に照射し、プローブ１１にて、複数の波長の光のそれぞれが照射された後に被検体内で発生した複数の光音響波を検出してもよい。その場合、光音響画像生成手段は、複数の波長の光に対応して検出された光音響信号をそれぞれ再構成して光音響画像を生成する。フーリエ変換手段２８は、複数の波長の光に対応する光音響画像のそれぞれをフーリエ変換する。空間周波数処理手段２９は、フーリエ変換された、複数の波長の光に対応する空間周波数領域の画像データから、それぞれ各波長に対応した所定の空間周波数成分を抽出し、フーリエ逆変換手段３０は、抽出された処理空間周波数成分のそれぞれをフーリエ逆変換する。

　例えば、レーザユニット１３から、第１の波長の光と第２の波長の光とを出射し、被検体に照射し、各波長の光が照射されたときに検出された光音響信号に基づいて光音響画像を生成してもよい。例えば、第１の波長は血液（血管）の画像化に対応した波長であり、第２の波長は薬剤（造影剤）の画像化に対応した波長であるとする。あるいは、第１の波長は動脈の画像化に対応した波長で、第２の波長は静脈の画像化に対応した波長であってもよい。

　上記のように、第１の波長と第２の波長とでは画像化する対象物が異なる場合、空間周波数処理手段２９は、第１の波長に対応した光音響画像データと第２の波長に対応した光音響画像データとに対し、互いに異なる周波数特性の空間周波数フィルタを用いることができる。つまり、第１の波長の光に対応する空間周波数領域の光音響画像データと、第２音波長の光に対応する空間周波数領域の光音響画像データとから、それぞれ異なる空間周波数成分を抽出してもよい。例えば動脈と静脈とで血管の太さの違いがあらかじめ分かっている場合は、第１の波長に対応する光音響画像と第２の波長に対応する光音響画像とから異なる空間周波数成分を抽出する。各波長の光で画像化する対象物に応じた空間周波数成分を抽出することで、各波長に対応した光音響画像から観察したい対象物を抽出することができる。

　次いで、本発明の第２実施形態を説明する。図１３は、本発明の第２実施形態の光音響画像生成装置を示す。本実施形態の光音響画像生成装置１０ａは、超音波ユニット１２ａ内に断面位置選択手段３５を更に有する点で、図１に示す第１実施形態の光音響画像生成装置１０と相違する。その他の点は、第１実施形態と同様でよい。

　断面位置選択手段３５は、ユーザ操作に応じて、画像表示手段１４の表示画面上に表示される断面画像データの断面位置を決定する。表示制御手段３１は、ユーザが指定した断面位置に対応する空間周波数処理された断面画像データを画像表示手段１４に表示させる。断面位置選択手段３５は、例えば図４に示す複数の断面画像データのうちで、どの位置の断面画像データを表示対象とするのかをユーザに選択させる。

　図１４は、表示画面例を示す。表示制御手段３１は、例えば断面位置を指定するためのスライドバーを画像表示手段１４に表示させる。断面位置選択手段が、ユーザのスライドバーの操作に応じて表示する断面位置を決定する。例えばスライドバーの紙面向って左側が被検体の浅い位置に対応し、右側が深い位置に対応する。ユーザがスライドバーを操作すると、その操作に応じて、スライドバーの上部に表示される空間周波数処理後の断面画像データの断面位置が変化する。スライドバーは断面画像データの下部に設けられている必要はなく、その上部や左右に設けられていてもよい。また、スライドバーは、断面画像データに付随して設けられている必要はなく、表示画面上のどこかに表示されていればよい。

　本実施形態では、断面位置選択手段３５は、ユーザの操作に応じて表示されるべき断面画像データの断面位置を選択し、表示制御手段３１は、選択された断面位置の空間周波数処理後の断面画像データを表示する。ユーザは、断面位置選択手段３５を操作することで、任意の断面位置における断面画像を表示させることができる。また、例えばスライドバーなどを用いて断面位置を連続的に変化させるように操作することで、ユーザは、光音響画像を動画的に観察することができる。その他の効果は第１実施形態と同様である。

　続いて、本発明の第３実施形態を説明する。図１５は、本発明の第３実施形態の光音響画像生成装置を示す。本実施形態の光音響画像生成装置１０ｂは、超音波ユニット１２ｂ内にデコンボリューション手段３６を更に有する点で、図１に示す第１実施形態の光音響画像生成装置１０と相違する。その他の点は、第１実施形態と同様でよい。

　デコンボリューション手段３６は、光音響画像再構成手段２４で再構成された光音響信号から、被検体に照射された光の光強度の時間波形の微分波形である光微分波形をデコンボリューションした信号を生成する。デコンボリューション手段３６は、例えば離散フーリエ変換により、再構成された光音響信号を時間領域の信号から周波数領域の信号へと変換する。また、光微分波形についても、離散フーリエ変換により時間領域の信号から周波数領域の信号へと変換する。デコンボリューション手段３６は、フーリエ変換された光微分波形の逆数を逆フィルタとして求め、フーリエ変換された周波数領域の光音響信号に逆フィルタを適用する。逆フィルタが適用されることで、周波数領域の信号において、光微分波形がデコンボリューションされる。その後、フーリエ逆変換により、逆フィルタが適用された光音響信号を、周波数領域の信号から時間領域の信号へと変換する。再構成後の光音響信号から光微分波形をデコンボリューションするのに代えて、再構成前の光音響信号から光微分波形をデコンボリューションするようにしてもよい。

　光微分波形のデコンボリューションについて説明する。光吸収体であるミクロ吸収粒子を考え、このミクロ吸収粒子がパルスレーザ光を吸収して圧力波（光音響圧力波）が生じることを考える。時刻をｔとして、位置ｒにあるあるミクロ吸収粒子から発生する光音響圧力波を、位置Ｒで観測した場合の圧力波形ｐ_{ｍｉｃｒｏ}（Ｒ，ｔ）は、[Phys. Rev. Lett. 86(2001)3550.]より、以下の球面波となる。

　ここで、Ｉ（ｔ）は励起光の光強度の時間波形であり、係数ｋは、粒子が光を吸収して音響波を出力する際の変換係数であり、ｖ_ｓは被検体の音速である。また、位置ｒ、Ｒは、空間上の位置を示すベクトルである。ミクロ吸収粒子から発生する圧力は、上記式に示すように、光パルス微分波形に比例した球面波となる。

　実際にイメージングする対象から得られる圧力波形は、よりマクロな吸収体のサイズを有しているため、上記のミクロ吸収波形を重ね合わせた波形になると考える（重ね合わせの原理）。ここで、マクロな光音響波を発する粒子の吸収分布をＡ（ｒ－Ｒ）とし、そのマクロな吸収体からの圧力の観測波形をｐ_{ｍａｃｒｏ}（Ｒ，ｔ）とする。観測位置Ｒでは、各時刻において、観測位置Ｒから半径ｖ_ｓｔに位置する吸収粒子からの光音響波が観測されることになるため、観測波形ｐ_{ｍａｃｒｏ}（Ｒ，ｔ）は、以下の圧力波形の式で示される。

　上記式（１）からわかるように、観測波形は、光パルス微分のコンボリューション型を示す。観測波形から光パルス微分波形をデコンボリューションすることで、吸収体分布が得られる。

　本実施形態では、検出された光音響信号から被検体に照射された光の微分波形をデコンボリューションする。光微分波形をデコンボリューションすることで、光吸収体の分布を得ることができ、吸収分布画像を生成することができる。吸収体の分布を画像化できることで、空間周波数処理後の断面画像において、血管等の位置をより確認しやすくなる。その他の効果は第１実施形態と同様である。

　なお、空間周波数処理後の光音響画像データのユーザへの提供の仕方は任意であり、図９に示すような断面画像を表示するものには限られない。例えば３次元の光音響画像データに対して所定の周波数区間成分を抽出する処理を行い、空間周波数処理後の３次元光音響画像データに対し、最大値投影法などを用いて投影画像を生成し、画像表示手段１４に表示するようにしてもよい。あるいは、３次元光音響画像データに対してボリュームレンダリングなどを行って３次元表示を行ってもよい。

　上記各実施形態では、断面画像生成手段２７を用い、フーリエ変換手段２８が断面画像に対して２次元フーリエ変換する例について説明したが、これには限定されない。断面画像生成手段２７を用いずに、フーリエ変換手段２８で３次元の光音響画像データに対して３次元フーリエ変換を行うようにしてもよい。この場合も、所定空間周波数領域を抽出することで、所望の太さの血管などを画像化できる。

　上記各実施形態では、フーリエ変換手段２８が、画像構築手段２６が生成した画像に対してフーリエ変換を行うものとして説明したが、フーリエ変換手段２８がフーリエ変換する信号は少なくとも光音響信号再構成手段２４で再構成されていればよく、画像構築手段２６が生成した画像には限定されない。フーリエ変換手段２８は、光音響画像再構成手段２４の出力信号に対してフーリエ変換を行ってもよい。その場合は、フーリエ逆変換手段３０が逆変換することで得られた信号に対して、検波・対数変換以後の処理を行うようにすればよい。

　上記各実施形態では、主に血管部分を画像化する例について説明したが、これには限定されない。例えば光音響画像で、神経やリンパ管などの管状の構造物を画像化する場合についても、空間周波数領域で所定の周波数成分を抽出することで、所望のサイズの構造物を画像化することができる。

　以上、本発明をその好適な実施形態に基づいて説明したが、本発明の光音響画像生成装置は、上記実施形態にのみ限定されるものではなく、上記実施形態の構成から種々の修正及び変更を施したものも、本発明の範囲に含まれる。

Claims

　被検体に照射すべき光を出射する光源と、
　被検体に対する光照射により被検体内で生じた光音響波を検出する音響波検出部と、
　前記光音響波の検出信号である光音響信号を再構成し画像データを生成する光音響画像生成部と、
　前記画像データを２次元以上の空間でフーリエ変換し、空間周波数領域の画像データを生成するフーリエ変換部と、
　前記空間周波数領域の画像データから所定の空間周波数成分を抽出する空間周波数処理部と、
　前記抽出された空間周波数成分をフーリエ逆変換し、空間周波数処理後の画像データを生成するフーリエ逆変換部とを備えたことを特徴とする光音響画像生成装置。
　前記空間周波数処理部が、第１の空間周波数以上で、かつ、第１の空間周波数よりも高い第２の周波数領域以下の空間周波数成分を抽出することを特徴とする請求項１に記載の光音響画像生成装置。
　前記第１の空間周波数及び第２の空間周波数が、空間周波数領域における位置に応じて変化することを特徴とする請求項２に記載の光音響画像生成装置。
　前記空間周波数処理部が、ユーザが指定した観察対象条件に従って、前記抽出する所定の空間周波数成分を決定することを特徴とする請求項１に記載の光音響画像生成装置。
　３次元的に検出された光音響波に基づく３次元画像データから、３次元空間を構成する軸のうちの１つの軸に垂直な面に沿って断面を切り出し、該切り出した断面を含む前記１つの軸に沿った方向の所定範囲の画像データを統合した断面画像データを生成する断面画像生成部を更に備え、
　前記フーリエ変換部が前記断面画像データを２次元フーリエ変換することを特徴とする請求項１に記載の光音響画像生成装置。
　前記断面画像生成部が、前記１つの軸に沿った複数の位置で断面を切り出して複数の断面画像データを生成し、
　前記フーリエ変換部が、前記生成された複数の断面画像データのそれぞれをフーリエ変換して複数の空間周波数領域の画像データを生成し、
　前記空間周波数処理部が、前記生成された複数の空間周波数領域の画像データのそれぞれから所定の空間周波数成分を抽出し、
　前記フーリエ逆変換部が、複数の空間周波数領域の画像データのそれぞれから抽出された所定の空間周波数成分をフーリエ逆変換し、前記複数の断面画像データのそれぞれについて空間周波数処理後の断面画像データを生成することを特徴とする請求項５に記載の光音響画像生成装置。
　前記断面画像生成部が、前記断面を等間隔で切り出すことを特徴とする請求項６に記載の光音響画像生成装置。
　前記断面画像生成部が、被検体に深さ方向に対応する軸に沿って、前記超音響波検出部の音響波検出面と平行な方向に前記断面を切り出すことを特徴とする請求項５に記載の光音響画像生成装置。
　前記断面画像生成部が、前記所定範囲内の画像データの最大値を投影し、又は前記所定範囲内の画像データを積分することで、前記所定範囲内の画像データを統合することを特徴とする請求項５に記載の光音響画像生成装置。
　前記空間周波数領域の断面画像データから抽出された所定の空間周波数成分をフーリエ逆変換することで得られた空間周波数処理後の断面画像データを表示装置に表示させる画像表示制御部を更に備えたことを特徴とする請求項５に記載の光音響画像生成装置。
　ユーザ操作に応じて、表示される断面画像データの断面位置を決定する断面位置選択部を更に備え、前記表示制御部が、ユーザが指定した断面位置に対応する空間周波数処理された断面画像データを表示装置に表示させることを特徴とする請求項１０に記載の光音響画像生成装置。
　前記表示制御部が、前記断面位置を指定するためのスライドバーを表示装置に表示させ、前記断面位置選択部が、ユーザのスライドバーの操作に応じて前記断面位置を決定することを特徴とする請求項１１に記載の光音響画像生成装置。
　前記表示制御部が、空間周波数処理前の断面画像データと、前記空間周波数処理後の断面画像データとを並べて表示装置に表示させることを特徴とする請求項１０に記載の光音響画像生成装置。
　前記表示制御部が、前記空間周波数処理後の断面画像データを２値化して表示装置に表示させることを特徴とする請求項１０に記載の光音響画像生成装置。
　前記フーリエ変換部が、３次元的に検出された光音響波に基づく３次元画像データを３次元フーリエ変換することを特徴とする請求項１に記載の光音響画像生成装置。
　前記光源が相互に異なる複数の波長の光を出射し、
　前記音響波検出部が、前記複数の波長の光のそれぞれが照射された後に被検体内で発生した複数の光音響波を検出し、
　前記光音響画像生成部が、前記複数の波長の光に対応して検出された光音響信号をそれぞれ再構成して画像データし、
　フーリエ変換部が、前記複数の波長の光に対応する画像データのそれぞれをフーリエ変換し、
　前記空間周波数処理部が、前記フーリエ変換された、複数の波長の光に対応する空間周波数領域の画像データから、それぞれ各波長に対応した所定の空間周波数成分を抽出し、
　前記フーリエ逆変換部が、前記抽出された処理空間周波数成分のそれぞれをフーリエ逆変換することを特徴とする請求項１に記載の光音響画像生成装置。
　前記複数の波長の光が第１の波長の光と第２の波長の光とを含み、前記空間周波数処理部が前記第１の波長の光に対応する空間周波数領域の画像データから抽出する空間周波数成分と、前記第２の波長の光に対応する空間周波数領域の画像データから抽出する空間周波数とが相互に異なることを特徴とする請求項１６に記載の光音響画像生成装置。
　前記光音響信号から、被検体に照射された光の微分波形をデコンボリューションするデコンボリューション部を更に備えたことを特徴とする請求項１に記載の光音響画像生成装置。
　被検体に対する光照射後に、該光照射により被検体内で生じた光音響波を検出するステップと、
　前記光音響波の検出信号である光音響信号を再構成して画像データを生成するステップと、
　画像データを２次元以上の空間でフーリエ変換し、空間周波数領域の画像データを生成するステップと、
　前記空間周波数領域の画像データから所定の空間周波数成分を抽出するステップと、
　前記抽出された処理空間周波数成分をフーリエ逆変換し、空間周波数処理後の画像データを生成するステップとを有することを特徴とする光音響画像生成方法。