(57)【特許請求の範囲】
【請求項１】
光源と、
前記光源から光を照射された被検体より発生する音響波を検出し、時系列の信号に変換する受信素子を含む検出器と、
光が照射されたときの光量分布を表わす情報を格納するメモリと、
前記光量分布を表わす情報を用いて、前記時系列の信号を補正する信号処理部と、
前記補正された時系列の信号を用いて被検体情報を形成する情報処理部と、
を有し、
前記信号処理部は、前記受信素子の中心から所定の角度範囲内の領域における前記光量分布を時間データに変換し、前記時間データを用いて前記時系列の信号の補正を行う
ことを特徴とする被検体情報取得装置。
【請求項２】
前記受信素子により変換された前記時系列の信号を、時系列のデジタル信号に変換する変換器を更に有し、
前記信号処理部は、前記時間データを用いて前記時系列のデジタル信号の補正を行う
ことを特徴とする請求項１に記載の被検体情報取得装置。
【請求項３】
前記信号処理部は、前記所定の角度範囲内の領域における前記光量分布を、各受信時刻に対応する、前記受信素子から等距離の位置における光量の合計値を表す前記時間データに変換する
ことを特徴とする請求項１または２に記載の被検体情報取得装置。
【請求項４】
前記信号処理部は、前記所定の角度範囲内の領域における前記光量分布を、前記受信素子の指向性に基づいて重みづけし、前記重みづけされた光量分布を前記時間データに変換する
ことを特徴とする請求項１から３のいずれか１項に記載の被検体情報取得装置。
【請求項５】
前記所定の角度範囲は、前記受信素子の指向性から決定された角度範囲である
ことを特徴とする請求項１から４のいずれか１項に記載の被検体情報取得装置。
【請求項６】
光源と、
前記光源から光を照射された被検体より発生する音響波を検出し、時系列の信号に変換する受信素子を含む検出器と、
光が照射されたときの光量分布を表わす情報を格納するメモリと、
前記光量分布を表わす情報を用いて、前記時系列の信号を補正する信号処理部と、
前記補正された時系列の信号を用いて被検体情報を形成する情報処理部と、
を有し、
前記信号処理部は、前記光量分布を前記受信素子の指向性に基づいて重みづけし、前記重みづけされた光量分布を時間データに変換し、前記時間データを用いて前記時系列の信号の補正を行う
ことを特徴とする被検体情報取得装置。
【請求項７】
前記受信素子により変換された前記時系列の信号を、時系列のデジタル信号に変換する変換器を更に有し、
前記信号処理部は、前記時間データを用いて前記時系列のデジタル信号の補正を行う
ことを特徴とする請求項６に記載の被検体情報取得装置。
【請求項８】
前記信号処理部は、前記重みづけされた光量分布を、各受信時刻に対応する、前記受信素子から等距離の位置における光量の合計値を表す前記時間データに変換する
ことを特徴とする請求項６または７に記載の被検体情報取得装置。
【請求項９】
前記時系列の信号を前記時間データで除算する
ことを特徴とする請求項１から８のいずれか１項に記載の被検体情報取得装置。  
【請求項１０】
前記信号処理部は、前記被検体内の各位置に照射される光量の差を低減するように、前記時系列の信号を補正する
ことを特徴とする請求項１から９のいずれか１項に記載の被検体情報取得装置。
【請求項１１】
前記信号処理部は、前記光量分布を算出する演算部を有する
ことを特徴とする請求項１から１０のいずれか１項に記載の被検体情報取得装置。
【請求項１２】
前記演算部は、前記被検体の光学特性値および形状を用いて、前記被検体内の光量分布を算出する
ことを特徴とする請求項１１に記載の被検体情報取得装置。
【請求項１３】
前記演算部は、前記時系列の信号を解析することにより、前記被検体の光学特性値を得る
ことを特徴とする請求項１２に記載の被検体情報取得装置。
【請求項１４】
前記信号処理部は、前記光量分布に基づいて、前記受信素子の受信面から垂直な方向に、前記被検体の深さと光量の関係を示す光量分布プロファイルを作成し、前記光量分布プロファイルと前記被検体内での音速とを用いて前記時間データを作成する
ことを特徴とする請求項１から１３のいずれか１項に記載の被検体情報取得装置。
【請求項１５】
前記情報処理部は、前記被検体情報に基づいて画像データを生成する
ことを特徴とする請求項１から１４のいずれか１項に記載の被検体情報取得装置。
【請求項１６】
光を照射された被検体より発生する音響波を検出し、時系列の信号に変換する工程と、
光が照射されたときの光量分布を表わす情報を用いて、前記時系列の信号を補正する工
程と、
前記補正された時系列の信号を用いて被検体情報を形成する工程と、
を有し、
前記時系列の信号を補正する工程は、
前記受信素子の中心から所定の角度範囲内の領域における前記光量分布を時間データに変換する工程と、
前記時間データを用いて前記時系列の信号の補正を行う工程と、を含む
ことを特徴とする被検体情報取得方法。
【請求項１７】
光を照射された被検体より発生する音響波を検出し、時系列の信号に変換する工程と、
光が照射されたときの光量分布を表わす情報を用いて、前記時系列の信号を補正する工程と、
前記補正された時系列の信号を用いて被検体情報を形成する工程と、
を有し、
前記時系列の信号を補正する工程は、
前記光量分布を前記受信素子の指向性に基づいて重みづけする工程と、
前記重みづけされた光量分布を時間データに変換する工程と、
前記時間データを用いて前記時系列の信号の補正を行う工程と、を含む
ことを特徴とする被検体情報取得方法。
【請求項１８】
前記時系列の信号を、時系列のデジタル信号に変換する工程を更に有し、
前記時系列の信号を補正する工程では、前記光量分布を表わす情報を用いて、前記時系列のデジタル信号を補正する
ことを特徴とする請求項１６または１７に記載の被検体情報取得方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被検体情報取得装置および被検体情報取得方法に関する。
【背景技術】
【０００２】
レーザーなどの光源から生体に光を照射し、入射した光に基づいて得られる生体内の情報を画像化する光イメージング装置の研究が医療分野で積極的に進められている。この光イメージング技術の一つとして、Photo Acoustic Tomography（ＰＡＴ：光音響トモグラ
フィー）がある。
【０００３】
光音響トモグラフィーでは、光源から発生したパルス光を生体に照射し、生体内で伝播・拡散したパルス光のエネルギーを吸収した生体組織から発生した音響波（典型的に超音波である）を検出する。すなわち、腫瘍などの被検部位とそれ以外の組織との光エネルギーの吸収率の差を利用し、被検部位が照射された光エネルギーを吸収して瞬間的に膨張する際に発生する弾性波を音響波検出器（探触子やトランスデューサーとも言われる）で受信する。この受信信号を解析処理することにより初期圧力分布あるいは光吸収エネルギー密度分布（吸収係数分布と光量分布の積）に比例した画像を得ることができる。
【０００４】
また、これらの画像情報は、様々な波長の光で計測することにより、被検体内の特定物質、例えば血液中に含まれるヘモグロビン濃度や血液の酸素飽和度などの定量的計測にも利用できる。近年、この光音響トモグラフィーを用いて、小動物の血管像をイメージングする前臨床研究や、この原理を乳がんなどの診断に応用する臨床研究が積極的に進められている。
【０００５】
光音響トモグラフィー技術を用いた装置では、画像化する領域に対して音響波の測定領域が不十分な場合がある。例えば、被検体の全周囲方向からではなく、ある特定の方向からしか音響波を受信できない場合である。このような場合、ストリーク（放射線）状のアーティファクトが画像上に出現することが知られている。
【０００６】
一方、特許文献１には、ボクセルあるいはピクセルごとに光量分布の画像データで除算し、吸収係数分布の画像データを求めるやり方が開示されている。
【０００７】
ここで、上記のような、音響波の測定領域が不十分な測定系のＰＡＴ装置で得られた光吸収エネルギー分布の画像データを、特許文献１のやり方で求めると、光量の弱い領域においてアーティファクトが強調された画像が出力される課題があった。
これは、画像では初期音圧像とアーティファクトの区別がつかないため、初期音圧像及びアーティファクトともに光量補正が行われるためである。
【０００８】
このような課題を解決するために、特許文献２では、音響波探触子から取得されるアナログ受信信号に対して、バリアブルゲインアンプを用いて光量補正を行い、その後に画像再構成により吸収係数分の画像データを算出する方法が開示されている。
【０００９】
このようにアナログ受信信号に対して、被検体の光量分布に基づいたゲインで補正すれば、吸収係数分布像のコントラストに寄与する受信信号のみを光量補正できる。そのため、初期音圧分布の画像データ自体を光量分布の画像データで補正する従来方法に比べて、アーティファクトの少ない画像を得ることができる。
【先行技術文献】
【特許文献】
【００１０】
【特許文献１】特開２０１０−０８８６２７号公報
【特許文献２】国際公開第２０１０／０２４２９０号
【発明の概要】
【発明が解決しようとする課題】
【００１１】
しかしながら，特許文献２ではアナログ受信信号をアンプで被検体の光量分布に対応した増幅を行うため、音圧の時間分解信号に対して任意の増幅を行うことができるバリアブルゲインアンプをすべての受信素子に対して用意する必要がある。その結果、装置のコストが増加する課題があった。
【００１２】
また、アナログ受信信号自体を増幅するため、光量補正していない受信信号を解析して、被検体の光量分布データを求めるようなアダプティブ処理を行う場合、つまり、光量分布データが測定前に求められない場合、受信信号の光量補正ができない課題があった。
【００１３】
本発明は上記の課題に鑑みてなされたものであり、その目的は、光音響トモグラフィーにおいて、装置のコストを抑制しつつ、光量補正によるアーティファクトの増加を低減するための技術を提供することである。
【課題を解決するための手段】
【００１４】
本発明は以下の構成を採用する。すなわち、光源と、前記光源から光を照射された被検体より発生する音響波を検出し、時系列の信号に変換する受信素子を含む検出器と、光が照射されたときの光量分布を表わす情報を格納するメモリと、前記光量分布を表わす情報を用いて、前記時系列の信号を補正する信号処理部と、前記補正された時系列の信号を用いて被検体情報を形成する情報処理部と、を有し、前記信号処理部は、前記受信素子の中心から所定の角度範囲内の領域における前記光量分布を時間データに変換し、前記時間データを用いて前記時系列の信号の補正を行うことを特徴とする被検体情報取得装置である。
本発明はまた、以下の構成を採用する。すなわち、光源と、前記光源から光を照射された被検体より発生する音響波を検出し、時系列の信号に変換する受信素子を含む検出器と、光が照射されたときの光量分布を表わす情報を格納するメモリと、前記光量分布を表わす情報を用いて、前記時系列の信号を補正する信号処理部と、前記補正された時系列の信号を用いて被検体情報を形成する情報処理部と、を有し、前記信号処理部は、前記光量分布を前記受信素子の指向性に基づいて重みづけし、前記重みづけされた光量分布を時間データに変換し、前記時間データを用いて前記時系列の信号の補正を行うことを特徴とする被検体情報取得装置である。
【００１５】
本発明はまた、以下の構成を採用する。すなわち、光を照射された被検体より発生する音響波を検出し、時系列の信号に変換する工程と、光が照射されたときの光量分布を表わす情報を用いて、前記時系列の信号を補正する工程と、前記補正された時系列の信号を用いて被検体情報を形成する工程と、を有し、前記時系列の信号を補正する工程は、前記受信素子の中心から所定の角度範囲内の領域における前記光量分布を時間データに変換する工程と、前記時間データを用いて前記時系列の信号の補正を行う工程と、を含むことを特徴とする被検体情報取得方法である。
本発明はまた、以下の構成を採用する。すなわち、光を照射された被検体より発生する音響波を検出し、時系列の信号に変換する工程と、光が照射されたときの光量分布を表わす情報を用いて、前記時系列の信号を補正する工程と、前記補正された時系列の信号を用いて被検体情報を形成する工程と、を有し、前記時系列の信号を補正する工程は、前記光量分布を前記受信素子の指向性に基づいて重みづけする工程と、前記重みづけされた光量分布を時間データに変換する工程と、前記時間データを用いて前記時系列の信号の補正を行う工程と、を含むことを特徴とする被検体情報取得方法である。
【発明の効果】
【００１６】
本発明によれば、光音響トモグラフィーにおいて、装置のコストを抑制しつつ、光量補正によるアーティファクトの増加を低減することが可能になる。
【図面の簡単な説明】
【００１７】
【図１】被検体情報取得装置の構成を模式的に示した図。
【図２】受信信号の処理の一例を説明するフロー図。
【図３】光量分布、光量分布と光量補正の一例を示す模式図。
【図４】第１の実施例における測定対象と画像の補正を示す模式図。
【発明を実施するための形態】
【００１８】
以下、図面を参照しつつ本発明をより詳細に説明する。なお、同一の構成要素には原則として同一の参照番号を付して、説明を省略する。
【００１９】
（被検体情報取得装置）
図１を参照しながら本実施形態にかかる被検体情報取得装置の構成を説明する。本実施形態の被検体情報取得装置は、被検体の内部の光学特性値に関する被検体情報を画像データとして取得する装置である。本発明において、音響波とは、典型的には超音波であり、音波、超音波、音響波、光音響波、光超音波と呼ばれる弾性波を含む。本発明の被検体情報取得装置とは、被検体に光（電磁波）を照射することにより被検体内で発生した音響波（典型的には超音波）を受信して、被検体情報を画像データとして取得する光音響効果を利用した装置を含む。取得される被検体情報とは、光照射によって生じた音響波の発生源分布や、被検体内の初期音圧分布、あるいは初期音圧分布から導かれる光エネルギー吸収密度分布や、吸収係数分布、組織を構成する物質の濃度情報分布を示す。物質の濃度情報分布とは、例えば、酸素飽和度分布や酸化・還元ヘモグロビン濃度分布などである。
【００２０】
本実施形態の被検体情報取得装置は、基本的なハード構成として、光源１１、音響波の検出器としての音響波探触子１７、信号処理部１９を有する。光源１１から発せられたパルス光１２は例えばレンズ、ミラー、光ファイバ、拡散板などの光学系１３により所望の光分布形状に加工されながら導かれ、生体などの被検体１５に照射される。被検体１５の内部を伝播した光のエネルギーの一部が血管などの光吸収体（結果的に音源となる）１４に吸収されると、その光吸収体１４の熱膨張により音響波（典型的には超音波）１６が発生する。これは「光音響波」と呼ばれることもある。音響波１６は音響波探触子１７により受信されてアナログの受信信号（電気信号）に変換される。そして、アナログ受信信号は、信号収集器１８で増幅やデジタル変換された後、信号処理部１９で被検体の画像データに変換される。また、その画像データは表示装置２０上で画像として表示される。
【００２１】
（光源１１）
被検体が生体の場合、光源１１からは生体を構成する成分のうち特定の成分に吸収される特定の波長の光を照射する。光源は、本実施形態の被検体情報取得装置と一体として設けられていても良いし、光源を分離して別体として設けられていても良い。光源としては数ナノから数百ナノ秒オーダーのパルス光を発生可能なパルス光源が好ましい。具体的には、効率的に光音響波を発生させるため、１０ナノ秒程度のパルス幅が使われる。光源としては大出力が得られるためレーザーが好ましいが、レーザーのかわりに発光ダイオードなどを用いることも可能である。レーザーとしては、固体レーザー、ガスレーザー、色素レーザー、半導体レーザーなど様々なレーザーを使用することができる。照射のタイミング、光強度などは不図示の光源制御部によって制御される。なお、この光源制御部は通常、光源と一体化されている。本発明において、使用する光源の波長は、被検体内部まで光が伝搬する波長を使うことが望ましい。具体的には、被検体が生体の場合、５００ｎｍ以上１２００ｎｍ以下である。
【００２２】
（光学系１３）
光源１１から照射された光１２は、光学系１３により、所望の光分布形状に加工されながら被検体に導かれる。光１２を、光ファイバなどの光導波路などを用いて伝搬させることも可能である。光学系１３は、例えば、光を反射するミラーや、光を集光したり拡大したり形状を変化させるレンズ、光を拡散させる拡散板、光ファイバなどである。このような光学部品は、光源から発せられた光１２が被検体１５に所望の形状で照射されれば、どのようなものを用いてもかまわない。なお、光はレンズで集光させるより、ある程度の面積に広げる方が生体への安全性ならびに診断領域を広げられるという観点で好ましい。
【００２３】
（被検体１５及び光吸収体１４）
これらは本発明の被検体情報取得装置の一部を構成するものではないが、以下に説明する。本発明の被検体情報取得装置は、血管の造影、人や動物の悪性腫瘍や血管疾患などの
診断や化学治療の経過観察などを主な目的とする。よって、被検体１５としては生体、具体的には人体や動物の乳房や指、手足などの診断の対象部位が想定される。動物では、ねずみなど小動物の場合は特定の部位だけではなく、小動物全体が対象となることもある。
【００２４】
被検体内部の光吸収体１４は、被検体内で相対的に吸収係数が高いものを示す。光吸収体１４には、使用する光の波長にもよるが、例えば、人体が測定対象であれば酸化あるいは還元ヘモグロビンやそれらを含む多く含む血管あるいは新生血管を多く含む悪性腫瘍が該当する。
【００２５】
（音響波探触子１７）
音響波探触子１７は、パルス光により被検体表面及び被検体内部などで発生する音響波を検出する検出器である。音響波探触子１７は、音響波を検知し、アナログ受信信号（以下、単に「アナログ信号」ともいう）である電気信号に変換する。以後、単に探触子あるいはトランスデューサということもある。圧電現象を用いたトランスデューサ、光の共振を用いたトランスデューサ、容量の変化を用いたトランスデューサなど音響波を検知できるものであれば、どのような音響波検出器を用いてもよい。本実施形態の音響波探触子１７は、典型的には複数の受信素子が１次元あるいは２次元に配置されたものが良い。このような多次元配列素子を用いることで、同時に複数の場所で音響波を検出することができ、計測時間を短縮できる。その結果、被検体の振動などの影響を低減できる。
なお、本発明において音響波探触子１７は、被検体１５に対して、特定の方向のみの音響波を受信できるように設置されることを想定している。
【００２６】
（信号収集器１８）
本実施形態の被検体情報取得装置は、探触子１７より得られたアナログ受信信号を増幅し、デジタル受信信号（以下、単に「デジタル信号」ともいう）に変換する信号収集器１８を有することが好ましい。信号収集器１８は、典型的には増幅器、Ａ／Ｄ変換器、ＦＰＧＡ（Field Programmable Gate Array）チップなどで構成される。探触子から得られる
受信信号が複数の場合は、同時に複数の信号を処理できることが望ましい。それにより、画像を形成するまでの時間を短縮できる。なお、本明細書において「受信信号」とは、探触子１７から取得されるアナログ信号も、その後ＡＤ変換されたデジタル信号も含む概念である。そして、光吸収体からの光音響波に起因する受信信号は「光音響信号」ともいう。
【００２７】
（信号処理部１９）
信号処理部１９は、主な役割として、信号収集器１８から得られたデジタル受信信号を用いて、画像再構成を行い、被検体の内部の光学特性値情報を画像データとして取得する。また、本発明の特徴的処理として、信号処理部１９は、信号収集器１８から得られたデジタル受信信号に対して光量分布に応じた補正（以下、「光量補正」ともいう）を行う。その結果、測定領域が制限された場合に得られた初期圧力分布の画像あるいは光吸収エネルギー吸収分布の画像に対して、光量補正を行わない場合と比較して、相対的にアーティファクトの少ない画像を出力することができる。なお、ここで言う測定領域が制限された場合とは、被検体を取り囲むすべての方向からではなく、図１の被検体１５と音響波探触子１７の関係のように、ある特定の方向のみからしか音響波を受信できないことを示している。
【００２８】
通常、信号処理部１９には典型的にはワークステーションなどのコンピュータが用いられ、あらかじめプログラミングされたソフトウェアにより、受信信号処理や画像再構成処理などが行われる。そのため、バリアブルゲインアンプなどのハードウェア構成を追加する場合に比べて、信号処理部１９で行われる処理は、低いコストで目的を達成できる。
【００２９】
ワークステーションで使われるソフトウェアは、例えば、受信信号に対してノイズ処理や光量補正などを行う信号処理モジュール１９ａと画像再構成を行う画像再構成モジュール１９ｂ、被検体の光量分布を算出する光量分布算出モジュール１９ｄからなる。また、本発明の信号処理部１９は、被検体の光量分布を保存するメモリ１９ｃを有し、信号処理モジュールではこのメモリ１９ｃに保存されたデータを参照しながら受信信号を補正する。ただし、信号処理モジュール１９ａ、画像再構成モジュール１９ｂ、メモリ１９ｃ及び光量分布算出モジュール１９ｄは、通常、ワークステーションなどのコンピュータ及び付属のソフトウェアとして扱われる。そのため、３つのモジュール及びメモリは、図１のように一つの信号処理装置１９として扱われることが多い。
【００３０】
信号処理モジュール１９ａでは、基本機能として、信号収集器１８から得られたデジタル受信信号に対して、被検体１５の光量分布データに応じた光量補正を行う。なお、被検体１５の光量分布データはメモリ１９ｃに保存されているものを用いる。保存されている光量分布データは光量分布算出モジュール１９ｄで計算されたもの、あるいは、予め他の場所で算出された光量分布データでも良い。ここで、光量分布データとは、光を照射された被検体内の各位置での光量を集積したデータを言う。
【００３１】
光量分布算出モジュール１９ｄでは、有限要素法やモンテ・カルロ法など生体などの被検体の光量分布を算出できる演算方法を用いて、被検体内の光量分布を算出する。メモリ１９ｃは光量分布算出モジュール１９ｄで算出された光量分布データや予め算出された光量分布データを格納する。光量分布算出モジュールは、本発明の演算部に相当する。
【００３２】
画像再構成モジュール１９ｂでは、基本機能として、信号処理モジュール１９ａから得られる光量補正後の受信信号データを用いて、画像再構成による画像データの形成が行われる。画像再構成アルゴリズムとして、例えば、トモグラフィー技術で通常に用いられるタイムドメインあるいはフーリエドメインでの逆投影などが使われる。なお、再構成に多くの時間をかけることが可能な場合は、繰り返し処理を行う逐次再構成法（iterative method）などの画像再構成手法も利用することができる。ＰＡＴの画像再構成手法には、代表的なものとして、フーリエ変換法、ユニバーサルバックプロジェクション法やフィルタードバックプロジェクション法、逐次再構成法などがあり、本発明においてはどのような画像再構成手法を用いても構わない。画像再構成モジュールは、本発明の画像処理部に相当する。
【００３３】
（表示装置２０）
表示装置２０は信号処理部１９で出力される画像データを表示する装置であり、典型的には液晶ディスプレイなどが利用される。なお、本発明の画像診断装置とは別に提供されていても良い。
【００３４】
（受信信号の処理）
次に、本発明の特徴的な処理である、信号処理部１９で行うデジタル受信信号の光量補正処理と画像再構成処理の一例について、図２と図３も参照しつつ説明する。下記の番号は図２のフロー図における処理の番号と一致する。
【００３５】
処理（１）（Ｓ２００）：被検体内の光量分布を算出する工程
【００３６】
光量分布算出モジュール１９ｄにて、被検体の吸収係数分布、等価散乱係数分布、光照射強度分布や被検体形状に応じて、被検体内部の光量分布を算出する。
【００３７】
被検体の吸収係数分布、等価散乱係数分布は、生体模擬ファントムなど予め物性値が既知の場合、その値を利用する。一方、吸収係数分布、等価散乱係数分布の値が未知の場合
、既知の方法を用いて受信信号から推定することができる。その場合、信号収集器１８から得られたデジタル受信信号を解析処理し、被検体の光学特性値である平均的な吸収係数や等価散乱係数を推定する。また、受信信号からの推定が困難な場合、文献などに記載されている統計データなどから値を推定することもできる。
光照射強度分布は、被検体に照射される領域での光強度の分布であり、装置構成及び被検体形状が決まれば予め求めることができる。
被検体形状はカメラなどの撮影データから求めることもできるし、圧迫板などを用いて、被検体の形状を特定の形状に整形することも可能である。
【００３８】
このようなデータ群が既知であれば、光輸送方程式あるいは光拡散方程式を解くことで、被検体の光量分布を計算により求めることができる。なお、光輸送方程式あるいは光拡散方程式はモンテ・カルロ法、差分法、有限要素法など数値計算分野ですでに知られている方法を用いることで、解くことが可能である。
【００３９】
上述したように被検体の光量分布算出方法には様々な方法を考えることができるが、ここでは一例をあげて説明する。例えば、被検体の形状、平均的な吸収係数及び等価散乱係数、光照射強度は予め求められていると仮定する。次に、それらの情報を用いて、光拡散方程式を有限要素法により解いた光量分布データを求めると、図３（ａ）のような光量分布データが求められる。この例では、図１のように光１２を音響波探触子１７とは反対方向から入射した場合の例を示している。
【００４０】
図３（ａ）においては、光照射近傍の光量が大きく、そこから離れるに従って、光が減衰し光量が小さくなっている。上記のような光量分布算出処理は、演算に必要な設定値である被検体の吸収係数分布、等価散乱係数分布、光照射強度分布、被検体形状が決まっていれば予め算出することができる。したがって、演算時間を低減するためにあらかじめ計算しておくこともできる。この場合、図１の信号処理器１９に示したように、光量分布算出モジュール１９ｄを装置に組み込む必要はない。ただし、被検体の吸収係数分布や散乱係数分布などが受信信号から推定される場合などは、信号収集器１８で信号が受信された後に算出する必要がある。その場合、図１に示したように、光量分布算出モジュール１９ｄを必要とする。
【００４１】
処理（２）（Ｓ２０１）：処理（１）で算出された光量分布データをメモリに保存する工程
【００４２】
信号処理器１９がワークステーションなどのコンピュータである場合、上記の処理で算出された光量分布データをハードディスク内のファイルに保存し、使用するときにそのデータをワークステーション内のメモリ１９ｃに格納する。あるいは、算出と同時に、その結果をメモリ１９ｃに格納する。このように光量分布データをメモリ１９ｃに格納できれば、どのような方法を用いても構わない。
【００４３】
処理（３）（Ｓ２０２）：メモリに保存された光量分布データを用いて、受信信号に対して光量分布（光量の大きさ）に応じた補正を行う工程
【００４４】
メモリに保存された光量分布データから音響波探触子１７の受信素子に対応する、前記受信素子の受信面に垂直な方向の光量分布プロファイルを算出する。次に、被検体の平均的な音速に応じて、前述の光量分布プロファイルの距離データを時間データに変換する。その時間データに変換された光量分布プロファイルを受信信号データに対して徐算することで、光量分布の影響を補正する。
【００４５】
光量分布プロファイルは例えば、対応する受信素子の指向性の範囲の中心に作成しても
良い。この場合、受信素子を２次元配列した音響波検出器であれば、受信面に垂直な方向に光量分布プロファイルが作成される。
【００４６】
上記の処理を、図３を参照しつつ詳しく説明する。図３（ｂ）は、図３（ａ）の点線ａの領域の光量分布プロファイルを、横軸を距離（超音波探触子１７の受信素子からの距離）、縦軸を光量としてプロットしたものである。なお、点線ａは、音響波探触子１７に含まれる一つの受信素子の受信面に対して垂直な方向である。ここでは、被検体を介して受信素子と反対側から光が照射されるため、受信素子から距離が離れるに従って、光量が増加している。受信素子側に近い位置は、照射された光が減衰しているため光量が小さい。
【００４７】
この図３（ｂ）のデータを被検体の平均音速を用いて、横軸を時間に変更する。つまり、距離を音速で割ればよい。そこで得られたデータが図３（ｃ）である。図３（ｃ）は図３（ｂ）のＸ軸が距離から時間（光照射から受信素子で受信されるまでの時間）に変更されただけで、グラフの形状などは変化しない。
すなわち、光量分布プロファイルは、ある受信素子から被検体に垂直な方向に線を引いた時に、その線上における光量を、被検体の深さを表わす値（距離や時間）と対応付けて示すデータである。したがって光量分布プロファイルは光量分布データから作成可能である。
【００４８】
なお、上記では各受信信号を補正する光量プロファイルを受信素子の受信面に対して垂直な領域とする例を示したが、光量プロファイル作成方法はこれに限られない。図３（ａ）の点線ｂは、一つの受信素子の中心とした、所定の角度θの範囲を示している。このとき、受信素子ごとに算出する光量プロファイルとして、この受信素子の中心から角度θの範囲にあり、かつ、受信素子の中心から等距離にある光量値を指向性で重み付けして積分した光量プロファイルを用いても良い。なお、この場合、角度範囲は受信素子が音響波を受信できる角度範囲であり、通常は、指向性から決定される。
【００４９】
図３（ｄ）は音響波探触子の受信素子で得られたデジタル受信信号の一例である。図３（ｄ）のグラフにおいて横軸は、右にいくほど、受信された時刻が早いことを意味する。通常、被検体１５内部にある光吸収体１４から発生する光音響波１６は、図３（ｄ）に示されているようなＮ型ライクな形状をしており、探触子に近い領域にある光吸収体から発生した光音響波ほど早く受信される。また、同じ形状、かつ、同じ吸収係数を持つ光吸収体から発生した光音響波ならば、受信素子に近いものほど、また、吸収する光量が大きいほど、大きな音圧で受信される。
【００５０】
通常、受信信号の音波の伝搬による音圧の減衰（回折効果）は、画像再構成において補正され、画像に反映される。一方、画像再構成においては光量の補正は行われないため、画像再構成後の画像は光量分布の影響を受けた画像となる。具体的に説明すると、図３（ｄ）のような受信信号データを再構成した後に得られる画像は、通常、初期音圧分布ｐ０（ｒ）である。初期音圧分布ｐ０（ｒ）は、下記の式（１）のように吸収係数μａ（ｒ）と光量φ（ｒ）の積に比例する。また、μａ（ｒ）とφ（ｒ）の積は光吸収エネルギー密度Ｈ（ｒ）と呼ばれる。
【数１】
0005645637.tif 000002
【００５１】
この式（１）から明らかなように、初期音圧分布ｐ０（ｒ）あるいは光吸収エネルギー密度Ｈ（ｒ）の値は、光量φ（ｒ）に影響される。従来の方法では、下の式（２）に示すように、再構成画像から吸収係数分布の画像データを求めるために、得られた画像データ
であるｐ０（ｒ）あるいはＨ（ｒ）をφ（ｒ）で除算し、μａ（ｒ）を求め、光量の影響を排除する。この場合、Ｈ（ｒ）あるいはｐ０（ｒ）の画像データにアーティファクトが存在すると、それも光量で補正されるため、光量が小さい領域ではそこに存在するアーティファクトが強調される課題がある。
【数２】
0005645637.tif 000003
【００５２】
一方、本発明では、この被検体内の各位置（各部位）に照射される光量の大きさに応じてデジタル受信信号を補正することが特徴である。つまり、被検体内の光量の差を低減するように、デジタル受信信号を補正する。具体的には、図３（ｄ）のデジタル受信信号に対して、図３（ｃ）の時間データに変換した光量分布データを受信時間ごとに除算する。例えば、この補正処理で得られるデータが図３（ｅ）である。この図から分かるように、光量の低い受信時間帯の信号が強調され、光量の影響が排除されている。このようにデジタル受信信号に対して光補正を行えば、光吸収体からの光音響波に起因する光音響信号のみに光補正を行うことが可能になる。
【００５３】
式を用いて上記を説明する。ここでは議論を単純化するために、データを最初から離散的に扱う。まず、連続値である初期音圧分布ｐ０（ｒ）の画像を離散化して、行列として扱い、初期音圧分布の行列表現をＰ０とする。また、初期音圧分布の離散化データから発生する光音響波を各受信素子で検出される音響波データに変換する感度行列をＹとする。なお、このＹはヤコビ行列である。また、各受信素子で検出される音圧データの行列表現をＤとすると、以下の式（３）が成り立つ。
Ｄ＝Ｙ・Ｐ０  …（３）
【００５４】
また、Ｐ０は上記で述べたように吸収係数分布と光量分布の積で表現されるため、Ｐ０＝Ｍａ・Φと表現することができる。ここでΦは光量分布を表現する行列であり、Ｍａは吸収係数分布を表現する行列である。つまり、式（３）は、以下の式（４）となる。
Ｄ＝Ｙ・Ｍａ・Φ  …（４）
両辺にΦの逆行列Φ−１をかければ、式（４）は、以下の式（５）となる。
Ｄ・Φ−１＝Ｙ・Ｍａ  …（５）
【００５５】
式（５）は、Φの逆行列を求め、それを各受信素子で検出される音圧データ行列Ｄにかけ、そのデータ（Ｄ・Φ−１）を用いて画像再構成することで、吸収係数分布の画像データＭａを得ることができることを示している。このとき、式（２）のように画像データに光量分布補正をする必要はない。その結果、受信信号に光量分布補正を行うことができるため、画像再構成後に発生するアーティファクトを増加させることはなくなる。
なお、Φの逆行列は、例えば、受信素子の中心から等距離にある光量値を指向性で重み付けして積分したもの、近似的には、受信素子の受信面に対して垂直な方向の光量分布プロファイルなどで表現できる。また、ここではΦの逆行列として一例をあげて示したが、各受信信号が被検体内の光量分布の影響を低減するように補正されれば、デジタル受信信号に対してどのような補正手段を用いてもかまわない。つまり被検体内の光量の差が低減されるように各受信信号が補正されれば、どのような補正手段を用いてもよい。
【００５６】
処理（４）（Ｓ２０３）：光量補正されたデジタル受信信号を用いて画像再構成処理を行い、被検体の光学特性値分布に関連した画像データを形成する工程
【００５７】
図３（ｅ）に示された光量補正されたデジタル受信信号を用いて画像再構成を行い、被検体１５の吸収係数分布に関連した画像データを形成する。この処理に関しては、通常の
光音響トモグラフィーで用いられるどのような画像再構成処理を用いることが可能である。例えば、タイムドメインあるいはフーリエドメインでの逆投影などである。
【００５８】
例えば、演算にタイムドメイン法の一つであるユニバーサルバックプロジェクション法を利用した場合、下記の式（６）に従って画像再構成を行う。
【数３】
0005645637.tif 000004
ここで、μａ（ｒ）は吸収係数分布、Ａは比例定数、ｂ（ｒ０，ｔ）は投影データ、d
Ω0は任意の観測点Pに対する検出器dS0の立体角である。この投影データを式（６）の積
分に従って逆投影することで吸収係数分布μａ（ｒ）の画像データを得ることができる。
【００５９】
また、b(r0,t)は、以下の式（７）である。
【数４】
0005645637.tif 000005
また、dΩ0は、以下の式（８）である。
【数５】
0005645637.tif 000006
ここで、ｐ’は光量分布補正された検出音波データ、θは検出器と任意の観測点Pとが
なす角度である。これらの式で表したように、本発明では、通常の画像再構成とは異なり、光量分布補正された検出音波データを投影データとして用いることが特徴である。
【００６０】
図１の例では、画像形成において、音響波の受信範囲が制限され、アーティファクトが発生する。そのような場合であっても、以上の工程を行うことで、初期音圧分布画像あるいは光吸収エネルギー密度分布画像に対して、光量補正を行う従来技術よりも、コストアップすることなく、アーティファクトによる画像劣化が少ない画像を提供できる。
【００６１】
＜実施例１＞
本実施形態を適用した光音響トモグラフィーを用いた被検体情報取得装置の一例について説明する。図１の装置概略図を用いて説明する。本実施例においては、光源１１として波長１０６４ｎｍで約１０ナノ秒のパルス光を発生するＱスイッチＹＡＧレーザーを用いた。パルス光１２から発せられる光パルスのエネルギーは０．６Ｊである。パルス光１２をミラーとビームエキスパンダーなどの光学システム１３を用いて、半径約２ｃｍ程度まで広げたあと、光を探触子１７とは反対側の被検体に照射できるように光学系１３をセッティングした。
【００６２】
被検体１５としては生体を模擬した図４（ａ）のような長方形ファントムを用いた。フ
ァントムは、１％に希釈したイントラリピッドを寒天で固めたものを利用した。このファントムは光学的に均質であり、その平均的な吸収係数は約０．０１ｍｍ−１、散乱係数は約０．７ｍｍ−１であった。また、このファントムのサイズは幅：５ｃｍ、高さ：５ｃｍ、奥行き：４ｃｍとした。
ファントム内には図４（ａ）に示されているように、直径０．３ｍｍの円柱状ゴムワイヤが光吸収体１４として８個、図のように埋め込まれている。なお、それぞれの光吸収体１４の吸収係数値はどれも同じである。また、ファントムや音響波探触子は音響マッチングのため、脱気された水で充たされた水槽の中に置かれている。
【００６３】
このようにセッティングされたファントムに対して、図１のように探触子１７とは反対側のファントム表面にパルス光１２を照射した。なお、音響波探触子１７としてはＰＺＴ（ジルコン酸チタン酸鉛）で作られた超音波トランスデューサを用いた。このトランスデューサは２次元アレイ型で、素子数は３４５（１５×２３）の素子、素子ピッチは２ｍｍである。また、素子の幅は約２ｍｍである。
【００６４】
図４（ａ）のように、パルス光１２をファントム表面に照射すると、光が円柱状光吸収体１４で吸収されることにより生じる光音響波が発生する。それらの光音響波を超音波トランスデューサ１７で、３４５チャンネル同時に受信した。その受信信号をアンプ、ＡＤコンバーター、ＦＰＧＡからなる信号収集器１８を用いて処理して、全チャンネルでの受信信号のデジタルデータを取得した。なお、信号のＳ／Ｎ比を向上させるために、３２回レーザーを照射し、得られたすべての受信信号を平均化した。その後、得られたデジタルデータを信号処理器１９であるワークステーション（ＷＳ）へ転送し、ＷＳ内に保存した。
【００６５】
次に、ＷＳ内のソフトウェアプログラムである信号処理モジュールにて、ファントムの光学特性値、ファントム形状、光照射強度分布を用いて、ファントム内の光量分布データを有限要素法により求めた。さらに、この光量分布データから各受信素子の受信面に対して垂直方向の光量分布データを算出し、各受信データを光量で補正した。ここでは、時間データに直した光量プロファイルで受信データを除算する補正方法を用いた。
【００６６】
さらに、この補正データを用いてＷＳ内のソフトウェアプログラムである再構成モジュール１９ａにて、画像再構成を行った。ここでは複数の画像再構成手法の中でタイムドメイン方式であるユニバーサルバックプロジェクション法を用いて３次元のボリュームデータを形成した。このとき使用したボクセル間隔は０．０２５ｃｍとした。画像化範囲は３．０ｃｍ×４．６ｃｍ×４．０ｃｍである。そのときに得られた画像の一例を図４（ｂ）に示す。なお、この図は３次元画像データにおいて、すべての吸収体が画像化できる方向の最大輝度を投影したＭＩＰ（Maximum Intensity Projection）像を示している。
【００６７】
次に、ＷＳ内に保存した光量補正していない受信信号を用いて、上記で採用された画像再構成手法を用いて、上記と同じ範囲である３．０ｃｍ×４．６ｃｍ×５．０ｃｍの画像を算出した。さらに、従来技術と同様に、その画像を光量分布で補正した。ここでは再構成画像データの各ボクセル値を光量分布データのボクセル値で除算する方法を用いた。そのときに得られた画像の一例を図４（ｃ）に示す。図４（ｃ）も３次元画像データから、すべての吸収体が画像化できる方向の最大輝度を投影したＭＩＰ像である。
【００６８】
図４（ｂ）と（ｃ）を比較する。図４（ｃ）では画像再構成後の光量補正のため、図中の点線で示した領域、つまり、光量の弱い領域のアーティファクトが目立っている。これは光量補正により、光量の少ない領域にあるすべての像が強調された結果である。一方、図４（ｂ）においては、受信信号に対して光量補正を行ったため、アーティファクトは強調されない。つまり、受信信号のみが光量補正されているため、アーティファクトを強調
することはない。このように、音響波の計測領域が制限された場合において、デジタル受信信号を光量補正することで、従来法よりもアーティファクトの少ない画像が得られる。
【００６９】
＜実施例２＞
本実施形態を適用した光音響トモグラフィーを用いた被検体情報取得装置の一例について説明する。本実施例においては、実施例１と同様なファントム及び測定系を用いた。ただし、本実施例では受信信号の解析により、被検体の光学特性値を求める。
【００７０】
まず、図１と同じように探触子１７とは反対方向から光を照射して光音響波を受信する。次に、図１には図示していないが探触子方向からも光を照射して光音響波を受信する。それぞれの受信信号に対して画像再構成を行い、２つの初期音圧分布画像を取得する。次に、それらの比をとり、その比のデータと計算により求めた光分布の比が最小になるように被検体の平均的な吸収係数と等価散乱係数を求めた。次に、求めた被検体の平均的な吸収係数と等価散乱係数、被検体形状、照射強度分布を用いて、光量分布算出モジュール１９ｃで被検体の光量分布を求めた。次に、この光量分布データで実施例１と同様な方法で、受信信号を光量補正し、その補正された受信信号を用いて画像再構成を行い、画像を得た。
このような方法で得られた画像は図４（ｂ）と同様なものであり、従来技術で得られた画像である図４（ｃ）よりもアーティファクトが低減されていた。
【００７１】
以上のことから、アダプティブに受信信号を解析して被検体の光学定数を求め、その値を利用して受信信号を光補正しても、従来法よりもアーティファクトの少ない画像が得られる。
【００７２】
＜実施例３＞
本実施形態を適用した光音響トモグラフィーを用いた被検体情報取得装置の一例について説明する。本実施例においては、実施例１と同様なファントム及び測定系を用いた。
【００７３】
まず、実施例１と同様にファントム内の光量分布を算出した。次に、各受信素子の中心から全角で３０度以内、かつ、受信素子の中心から等距離にある光量分布の値を検出素子の指向性で重み付けしたものを積分した光量プロファイルを算出した。この光量プロファイルを用いて、実施例１と同様な方法で各受信信号データを補正し、補正した受信信号データを作成した。この補正した受信信号データを用いて、画像再構成を行い、相対的な吸収係数分布画像を算出した。
その結果、従来技術よりもアーティファクトが低減され、かつ、図４（ｂ）よりも定量性の向上した吸収係数分布画像を得ることができた。
【００７４】
以上のことから、探触子の指向性を考慮した光量プロファイルで受信信号を補正することで、定量性も保ちつつ、従来技術よりもアーティファクトが低減された光学特性値分布画像を得ることができる。
【符号の説明】
【００７５】
１１：光源，１７：音響波探触子，１９：信号処理部，１９ａ：信号処理モジュール，１９ｂ：画像再構成モジュール，１９ｃ：メモリ，１９ｄ：光量分布算出モジュール
【図１】
0005645637.tif 000007
【図２】
0005645637.tif 000008
【図３】
0005645637.tif 000009
【図４】
0005645637.tif 000010