(57)【要約】      （修正有）【課題】被検体内部で変位が生じた場合であっても、複数のタイミングで取得した受信信号を用いて特性情報を精度良く再現するための技術を提供する。【解決手段】光源１１０と、光源から出射された光が照射されることにより被検体から発生する光音響波を受信して時系列の電気信号に変換する受信手段１３０と、時系列の電気信号を用いて前記被検体内部の特性情報を取得する処理手段１５０と、を有し、光源は、複数のタイミングで光を出射し、受信手段は、複数のタイミングで光音響波を受信し、処理手段は、複数のタイミングに対応する複数の時系列の電気信号のそれぞれに対して、複数のタイミング間の被検体内部の変位量に基づいて、時系列の電気信号が投影される投影位置座標を決定し、複数の時系列の電気信号を決定された投影位置座標へ投影して特性情報を取得する。【選択図】図１
2014224806.tif 000002
【特許請求の範囲】
【請求項１】
光源と、
前記光源から出射された光が照射されることにより被検体から発生する光音響波を受信して時系列の電気信号に変換する受信手段と、
前記時系列の電気信号を用いて前記被検体内部の特性情報を取得する処理手段と、
を有し、
前記光源は、複数のタイミングで光を出射し、
前記受信手段は、前記複数のタイミングで前記光音響波を受信し、
前記処理手段は、前記複数のタイミングに対応する複数の時系列の電気信号のそれぞれに対して、前記複数のタイミング間の前記被検体内部の変位量に基づいて、前記時系列の電気信号が投影される投影位置座標を決定し、
前記複数の時系列の電気信号を前記決定された投影位置座標へ投影して前記特性情報を取得する
ことを特徴とする被検体情報取得装置。
【請求項２】
前記複数のタイミング間の前記被検体内部の変位量を検出する第１の検出手段をさらに有する
ことを特徴とする請求項１に記載の被検体情報取得装置。
【請求項３】
前記第１の検出手段は、前記被検体に送信されたのち前記被検体内で反射した超音波に基づく超音波信号を用いて前記複数のタイミング間の前記被検体内部の変位量を検出することを特徴とする請求項２に記載の被検体情報取得装置。
【請求項４】
前記第１の検出手段は、前記複数のタイミングのそれぞれにおいて取得された複数の前記超音波信号に基づいて前記複数のタイミング間の前記被検体内部の変位量を検出する
ことを特徴とする請求項３に記載の被検体情報取得装置。
【請求項５】
前記第１の検出手段は、スペックルトラッキング法により、前記複数のタイミング間の前記被検体内部の変位量を検出する
ことを特徴とする請求項４に記載の被検体情報取得装置。
【請求項６】
前記受信手段は、前記被検体に超音波を送信する機能と、前記被検体内で反射した超音波を受信する機能を有する
ことを特徴とする請求項３に記載の被検体情報取得装置。
【請求項７】
前記複数のタイミングのそれぞれにおける前記受信手段と前記被検体との相対位置を検出する第２の検出手段をさらに有し、
前記処理手段は、前記複数のタイミングのそれぞれにおける前記受信手段と前記被検体との相対位置、および前記複数のタイミングのそれぞれにおける前記被検体内部の変位量に基づいて、前記投影位置座標を決定する
ことを特徴とする請求項１ないし６のいずれか１項に記載の被検体情報取得装置。
【請求項８】
前記第２の検出手段は、前記被検体に送信されたのち前記被検体内で反射した超音波に基づく超音波画像にブロックマッチング法を適用して、前記複数のタイミングのそれぞれにおける前記受信手段と前記被検体との相対位置を検出する
ことを特徴とする請求項７に記載の被検体情報取得装置。
【請求項９】
前記受信手段は、磁気センサまたは光学センサを有し、
前記第２の検出手段は、前記磁気センサまたは前記光学センサの出力に基づいて、前記
複数のタイミングのそれぞれにおける前記受信手段と前記被検体との相対位置を検出することを特徴とする請求項７に記載の被検体情報取得装置。
【請求項１０】
前記特性情報を用いて前記被検体内部の画像を表示する表示手段をさらに有する
ことを特徴とする請求項１から９のいずれか１項に記載の被検体情報取得装置。
【請求項１１】
前記受信手段は、前記複数のタイミングにおいて前記受信手段と前記被検体との相対位置を変更可能に構成されている
ことを特徴とする請求項７に記載の被検体情報取得装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被検体情報取得装置に関する。
【背景技術】
【０００２】
レーザなどの光源から生体等の被検体に光を照射し、入射した光に基づいて得られる被検体内の情報を画像化する光イメージングの研究が医療分野で進められている。この光イメージング技術の一つとして、Ｐｈｏｔｏ  Ａｃｏｕｓｔｉｃ  Ｉｍａｇｉｎｇ（ＰＡＩ：光音響イメージング）がある。光音響イメージングでは、光源から照射されたのち、被検体内で伝播・拡散したパルス光のエネルギーを吸収した組織から発生した弾性波が受信される。そして、受信信号に基づき被検体内の特性情報が画像化される。
【０００３】
光音響イメージングにおいて被検体に光が照射されると、被検部位（生体の腫瘍等）とそれ以外の組織では光エネルギーの吸収率に差があるため、被検部位が光エネルギーを吸収して瞬間的に膨張し、弾性波を発生させる。この弾性波を音響受信器で受信した信号を解析処理することにより、特性情報が得られる。ここで、特性情報とは、初期音圧分布、光吸収エネルギー密度分布、光吸収係数分布などの光学特性値分布である。また、これらの情報を複数の波長の光で計測することにより、被検体内の特定物質（血液中のヘモグロビン濃度や、血液の酸素飽和度など）の定量的計測にも利用できる。
【０００４】
これらの特性情報を正確に算出するためには、被検体内で三次元的に発生した弾性波を全周囲から受信して、それらの受信信号を使用して画像再構成することが望まれる。しかし多くの場合、音響受信器と測定部位の物理的な制約から、弾性波を全周囲から受信することは困難である。その結果、限られた受信条件から特性情報を推定することとなるため、画像化した特性情報にデータ欠落や再構成アーティファクトといったノイズが発生する。そこで、この課題を鑑みて特許文献１では、光音響波の取得を複数回、異なる位置で行い、受信信号に対して画像再構成して特性情報を取得することが提案されている。ここで画像再構成とは、任意に抽出した受信信号（もしくは受信信号を任意に重みづけなど処理した投影信号）を各再構成ピクセル（ボクセル）に割り当てる（投影する）処理のことである。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０１０−０８８４９７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、実際に生体組織などに対して多数回の測定を行うと拍動や血流などの動きを伴い、それぞれのタイミングの間に被検体内部で変位が生じてしまう可能性がある。この場合、それぞれのタイミングで取得した受信信号を用いて特性情報を正確に推定して再現することは困難である。
【０００７】
本発明は上記課題に鑑みてなされたものである。本発明の目的は、被検体内部で変位が生じた場合であっても、複数のタイミングで取得した受信信号を用いて特性情報を精度良く再現するための技術を提供することである。
【課題を解決するための手段】
【０００８】
本発明は、以下の構成を採用する。すなわち、光源と、前記光源から出射された光が照
射されることにより被検体から発生する光音響波を受信して時系列の電気信号に変換する受信手段と、前記時系列の電気信号を用いて前記被検体内部の特性情報を取得する処理手段と、を有し、前記光源は、複数のタイミングで光を出射し、前記受信手段は、前記複数のタイミングで前記光音響波を受信し、前記処理手段は、前記複数のタイミングに対応する複数の時系列の電気信号のそれぞれに対して、前記複数のタイミング間の前記被検体内部の変位量に基づいて、前記時系列の電気信号が投影される投影位置座標を決定し、前記複数の時系列の電気信号を前記決定された投影位置座標へ投影して前記特性情報を取得することを特徴とする被検体情報取得装置である。
【発明の効果】
【０００９】
本発明によれば、被検体内部で変位が生じた場合であっても、複数のタイミングで取得した受信信号を用いて特性情報を精度良く再現するための技術を提供できる。
【図面の簡単な説明】
【００１０】
【図１】本発明に係る被検体情報取得装置の模式図。
【図２】本発明に係る被検体情報取得方法のフローを示す図。
【図３】実施例１に係る被検体情報取得装置のタイムチャートを説明する図。
【図４】実施例１に係る音響受信器の操作範囲を示す図。
【図５】実施例２に係る被検体情報取得装置の処理フローを示す図。
【図６】実施例４に係る被検体情報取得装置の処理フローを示す図。
【図７】実施例１に係る被検体情報取得装置の表示画面を示す図。
【図８】実施例２に係る被検体情報取得装置の表示画面を示す図。
【図９】実施例４に係る被検体情報取得装置のメモリ保存領域を示す図。
【図１０Ａ】音響受信器と信号受信する領域の関係を説明する図。
【図１０Ｂ】音響受信器と信号受信する領域の関係を説明する別の図。
【図１０Ｃ】音響受信器と信号受信する領域の関係を説明する別の図。
【図１０Ｄ】音響受信器と信号受信する領域の関係を説明する別の図。
【図１１】実施例３に係る被検体情報取得装置の処理フローを示す図。
【図１２】実施例５に係る被検体情報取得装置の処理フローを示す図。
【図１３】実施例５に係る音響受信器を示す図。
【発明を実施するための形態】
【００１１】
以下に図面を参照しつつ、本発明の好適な実施の形態について説明する。ただし、以下に記載されている構成部品の寸法、材質、形状およびそれらの相対配置などは、発明が適用される装置の構成や各種条件により適宜変更されるべきものであり、この発明の範囲を以下の記載に限定する趣旨のものではない。
【００１２】
本発明の被検体情報取得装置は、被検体に光（電磁波）を照射することにより被検体内で発生し伝播した音響波を受信して、被検体内部の特性情報を画像データとして取得する光音響効果を利用した装置である。本発明の被検体情報取得装置は、超音波の送信機能と、被検体内からの反射波（エコー波）を受信する機能を有し、画像の改善などのためにエコー波から得られる情報を使用できる。
【００１３】
取得される特性情報とは、光照射によって生じた音響波の発生源分布、被検体内の初期音圧分布、あるいは初期音圧分布から導かれる光エネルギー吸収密度分布や吸収係数分布、組織を構成する物質の濃度分布を示す。組織を構成する物質とは、例えば、酸素飽和度分布や酸化・還元ヘモグロビン濃度分布などの血液成分、あるいは脂肪、コラーゲン、水分などである。また、エコー波に基づく電気信号に対して既知の情報処理を行うことにより得られる、被検体内部の音響インピーダンスの分布を示す情報のことも、一種の特性情報と捉えてもよい。
【００１４】
本発明でいう音響波とは、典型的には超音波であり、音波、音響波と呼ばれる弾性波を含む。光音響効果により発生した音響波のことを、光音響波または光超音波と呼ぶ。これ以降、音響波の中でも、音響受信器によって被検体内に送信されたものや、送信後、被検体内で反射したものを指すとき、「超音波」と記載する。後者は「エコー波」と記載する場合もある。また、音響波の中でも、光照射により被検体内で発生したものを指すとき、「光音響波」と記載する。ただしこれは両者を区別する便宜上のものであり、それぞれの音響波の波長等を限定するものではない。
【００１５】
本発明では音響受信器の被検体に対する方向や位置を変更させながら複数の光音響波を時系列にサンプリングし、時系列に得られたサンプリングデータを用いて再構成することが好ましい。方向等の変更は、例えば音響受信器としてのハンドヘルド探触子を、作業者が揺動させながら被検体に当てることで実現できる。これにより、一回の（あるいは特定方向のみからの）測定によるデータを用いた場合に比べて、より開口を広げた再構成が行われる。すなわち、音響受信器を被検体内部の関心領域に対して揺動させることで、三次元的に発生した音響波をより大きな立体角で捉えられる。その結果、複数の位置で得られた光音響波の時系列の受信信号を用いた再構成により、実際の特性情報をより正確に推定して再現できる。
【００１６】
しかしながら、上記方法では、音響受信器の被検体に対する方向や位置を変更させながら複数のタイミングで光音響波を受信する必要がある。そのため、実際に生体組織などに対して多数回の測定を行うと拍動や血流などの動きを伴い、それぞれのタイミングの間に被検体内部で変位が生じてしまう可能性がある。この場合に、被検体内部の変位を考慮せずに時系列の受信信号を用いた再構成により特性情報を取得すると、特性情報を正確に推定して再現することは困難である。
【００１７】
そこで、本発明は、複数のタイミング間の被検体内部の変位量に基づいて、複数のタイミングで得られた各時系列の受信信号が投影される投影位置座標を決定する。すなわち、本発明は、複数のタイミングで得られた複数の時系列の受信信号の各サンプリングデータに対して、複数のタイミング間の被検体内部の変位量に基づいて、各サンプリングデータが投影される投影位置座標を決定する。
【００１８】
これにより、複数のタイミング間で被検体内部に変位が生じてしまった場合でも、被検体内部の変位を考慮して決定された投影位置座標に基づいて実際の特性情報を正確に推定して再現することができる。
なお、画像再構成により複数のタイミングで得られた時系列の受信信号に基づいて特性情報を取得した後では、複数のタイミング間の被検体内部の変位による誤差を補正することは困難である。
【００１９】
（装置構成）
図１の模式図を用いて、本発明に係る被検体情報取得装置を説明する。装置は、光源１１０、光学系１２０、音響受信器１３０、制御装置１４０、信号処理装置１５０、表示装置１６０を有する。以下、各構成要素について説明する。
【００２０】
（被検体１００及び光吸収体１０１）
これらは本発明の被検体情報取得装置の一部を構成するものではないが、以下に説明する。本発明の被検体情報取得装置は、人や動物の悪性腫瘍や血管疾患などの診断や化学治療の経過観察などを主な目的とする。よって、被検体としては生体、具体的には人体や動物の乳房や頸部、腹部などの診断の対象部位が想定される。
また、被検体内部にある光吸収体とは、被検体内部で相対的に吸収係数が高いものを指
す。例えば、人体が測定対象であればオキシヘモグロビンあるいはデオキシヘモグロビンや、それらを含む多く含む血管、あるいは新生血管を多く含む悪性腫瘍が光吸収体となる。その他、頸動脈壁のプラークなども光吸収体となる。
【００２１】
（光源１１０）
光源としては、数ナノから数マイクロ秒オーダーのパルス光を発生可能なパルス光源が好ましい。具体的には効率的に光音響波を発生させるため、１０ナノ秒程度のパルス幅が使われる。光源としてはレーザのかわりに発光ダイオードなども使用できる。レーザとしては、固体レーザ、ガスレーザ、色素レーザ、半導体レーザなど様々なレーザを使用できる。照射光の波長は、被検体内部まで光が伝搬する波長が望ましい。具体的には、被検体が生体の場合、５００ｎｍ以上１２００ｎｍ以下である。
【００２２】
（光学系１２０）
光源から出射された光は、典型的にはレンズやミラーなどの光学部品により、所望の光分布形状に加工されながら被検体に導かれる。また、光ファイバなどの光導波路などを用いて光を伝搬させることも可能である。光学系は、例えば、光を反射するミラーや、光を集光したり拡大したり形状を変化させるレンズ、光を拡散させる拡散板などである。このような光学部品は、光源から発せられた光が被検体に所望の形状で照射されれば、どのようなものを用いてもかまわない。なお、光はレンズで集光させるより、ある程度の面積に広げる方が、生体への安全性ならびに診断領域を広げられるという観点で好ましい。
【００２３】
（音響受信器１３０）
音響受信器は、音響波（光音響波およびエコー波）を受信し、アナログの電気信号に変換する。圧電現象、光の共振、静電容量の変化等を用いたものなど、どのような音響受信器を用いてもよい。ここで、音響受信器は、音響波を送受信するトランスデューサが、アレイ状に複数個配列されたものが好ましい。これにより、複数の位置で音響波を受信して複数の信号を出力できるので、測定時間の短縮やＳＮ比の向上が期待できる。
音響受信器は通常、トランスデューサが筐体に保持された探触子の形態で提供される。なお、本明細書において音響受信器を探触子とも呼ぶ。
【００２４】
音響受信器は、被検体に対して位置を変更可能であることが好ましい。位置の変更方法は、ハンドヘルド探触子であれば作業者の手技による。ただし、音響受信器の走査機構を設置することも可能である。被検体情報取得装置が、作業者に対し、音声や画面表示の形式で揺動のガイド指示を与えることも可能である。
音響受信器は、被検体に超音波を送信する送信器の機能と、被検体の内部を伝搬したエコー波を受信する受信器の機能を備えることが好ましい。これにより、同一領域での信号検知や省スペース化が期待できる。ただし、送信器と受信器を別にしてもよい。また、光音響波とエコー波の受信器を別にしてもよい。音響受信器は、本発明の受信手段に相当する。
【００２５】
（制御装置１４０）
制御装置は、音響受信器１３０より出力されたアナログの電気信号に対して、増幅処理と、デジタル変換処理を施す。制御装置は、典型的には増幅器、Ａ／Ｄ変換器、ＦＰＧＡ（Ｆｉｅｌｄ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ）チップ、ＣＰＵなどで構成される。音響受信器１３０が複数のトランスデューサを有し、複数の受信信号を出力する場合、制御装置も、複数の信号を同時に処理できることが望ましい。それにより、処理時間を短縮できる。
また制御装置は、パルス光の照射と、それに引き続く光音響波の受信のタイミングを制御する。具体的には、パルス光の照射タイミングの制御や、パルス光をトリガ信号とした電気信号の送受信タイミングの制御を行う。
【００２６】
ここで、本発明において、光音響波信号とは、音響受信器１３０から出力されたアナログの時系列の電気信号も、制御装置１４０で処理された時系列の信号も含む概念である。また、超音波信号とは、エコー波を受信した音響受信器１３０から出力されるアナログの時系列の電気信号と、制御装置１４０での処理後の時系列の信号を含む概念である。
【００２７】
（信号処理装置１５０）
信号処理装置は、デジタル信号に基づき被検体内部の情報を生成する。信号処理装置としては、典型的にはワークステーションなどの情報処理装置が用いられる。後述する補正処理や画像再構成処理などは、あらかじめプログラミングされたソフトウェアにより行われる。
ソフトウェアは、本発明の特徴的な処理である補正処理を行う変位量算出モジュール１５１を含んでいる。ソフトウェアはまた、画像再構成モジュール１５２を含んでいる。なお、それぞれのモジュールを、信号処理装置１５０とは別の装置として設けてもよい。また信号処理装置１５０は２Ｄ空間、３Ｄ空間のいずれにも信号処理を適用できる。
【００２８】
変位量算出モジュール１５１は、音響受信器１３０の変位量、例えば作業者の手技によるハンドヘルド探触子の方向や位置の変化量を算出する。また、ハンドヘルド探触子の押し付けや体動などに起因する、被検体内部の変位量も算出できる。
【００２９】
画像再構成モジュール１５２は、光音響波信号を用いて画像再構成を行い、特性情報を形成する。また、超音波信号を用いて、被検体の音響インピーダンス分布などの特性情報を生成する。
画像再構成アルゴリズムとしては、トモグラフィー技術で既知の手法を用いる。例えば、タイムドメインあるいはフーリエドメインでの逆投影や、整相加算（ディレイ・アンド・サム）などである。再構成の時間を多く使える場合は、繰り返し処理による逆問題解析法などの画像再構成手法を使用してもよい。
また、画像再構成モジュール１５２は、超音波信号に対して、位相を整合するための遅延加算処理や、その後の加算処理を行う。これにより、被検体内での音響インピーダンス等の特性情報や、被検体内での散乱に起因するスペックルパターンデータを形成できる。
【００３０】
なお、光音響イメージングにおいて、フォーカスした音響受信器１３０を用いることで、画像再構成なしに生体内の特性情報画像を形成できる。そのような場合には、画像再構成アルゴリズムを用いた信号処理は必要ない。
また、変位量算出モジュール１５１および画像再構成モジュール１５２のそれぞれは、典型的にはＣＰＵ、ＧＰＵなどの素子や、ＦＰＧＡ、ＡＳＩＣなどの回路から構成される。なお、変位量算出モジュール１５１および画像再構成モジュール１５２のそれぞれは、１つの素子や回路から構成されるだけではなく、複数の素子や回路から構成されていてもよい。また、変位量算出モジュール１５１および画像再構成モジュール１５２が行う各処理をいずれの素子や回路が実行してもよい。
【００３１】
また、制御装置、信号処理装置は一体化される場合もある。この場合、ワークステーションで行うようなソフトウェア処理ではなく、ハードウェア処理により被検体の音響インピーダンスなどの特性情報や、光学特性値分布を生成できる。制御装置および信号処理装置の画像再構成モジュールは、本発明の処理手段に相当する。変位量算出モジュールは、本発明の検出手段に相当する。
【００３２】
（表示装置１６０）
表示装置は、信号処理装置から出力される光学特性値分布等の特性情報を表示する。表示装置としては例えば、液晶ディスプレイ、プラズマディスプレイ、有機ＥＬディスプレ
イ、ＦＥＤなどを使用できる。なお、表示装置は、被検体情報取得装置の本体とは別に提供されていても良い。例えば、取得した被検体情報を有線または無線で表示装置に表示させてもよい。表示装置は、本発明の表示手段に相当する。
【００３３】
（処理フロー）
次に、図２のフロー図を参照して、被検体情報取得装置を制御して行われる被検体情報取得方法の手順を説明する。なお、制御装置が、制御装置内のメモリに保存された被検体情報取得方法が記述されたプログラムを読み出し、被検体情報取得装置の各構成の作動を制御することにより以下のフローを実行させる。
【００３４】
（Ｓ１００：光音響波信号を取得する工程）
この工程では、音響受信器１３０が、被検体から発生した光音響波を受信し、光音響波信号を生成する。
まず、光源１１０から出射されたパルス光１２１が光学系１２０を介して被検体１００に照射される。パルス光１２１は光吸収体１０１に吸収され、光音響波１０３が発生する。制御装置１４０はパルス光の出射を検知して、音響受信器１３０に光音響波１０３の受信を開始させる。音響受信器１３０から出力された光音響波信号は、制御装置１４０での処理を経て、メモリに格納される。
なお、本実施形態においては、本工程を複数のタイミングで行う。
【００３５】
（Ｓ２００：超音波信号を取得する工程）
この工程では、被検体に対して超音波を送受信することにより、被検体の音響インピーダンスなどの特性情報や、被検体内での散乱に起因するスペックルパターンなどのデータを取得する。音響インピーダンスは、例えば、Ｓ５００における音響受信器１３０の変位量算出に利用できる。スペックルパターンは、例えば、ハンドヘルド探触子の押し付け等による被検体内部の変位量算出に利用できる。
音響受信器１３０は、被検体１００に対して超音波１０２ａを送信する。送信した超音波は被検体内で反射して、エコー波１０２ｂとなる。音響受信器１３０は、そのエコー波１０２ｂを受信し、超音波信号を出力する。そして、制御装置１４０が増幅およびＡ／Ｄ変換処理を行ったのち、メモリに格納する。
超音波１０２ａを送信するタイミングは、Ｓ１００におけるパルス光１２１と同時でも構わない。また、超音波を送受信する回数は１回に限らず、次のパルス光照射までの間に複数回行ってもよい。
【００３６】
（Ｓ３００：Ｓ１００、Ｓ２００の繰り返し回数を判定する工程）
この工程では、信号処理装置１５０において、Ｓ１００での光音響波信号の取得回数と、Ｓ２００での超音波信号の取得回数が、所定数に達しているかどうかが判定される。所定数に達していなければＳ１００、Ｓ２００を繰り返す。
繰り返し回数は予め所望の回数を設定しておく。あるいは、測定開始前に作業者がタッチパネル、キーボード等の入力手段から入力してもよい。あるいは、外部に設置した押しボタンを押し込んでいる間は測定を繰返すように設定するなど、測定をしながら回数を決定可能としてもよい。測定をしながら回数を決定する場合の一例として、ハンドヘルド探触子に接触センサを設けておき、作業者がハンドヘルド探触子を被検体に当てている間は測定が行われるようにしてもよい。
【００３７】
（Ｓ４００：光音響波の関心領域を設定する工程）
この工程では、画像再構成モジュール１５２が、関心領域（ＲＯＩ：Ｒｅｇｉｏｎ  ｏｆ  Ｉｎｔｅｒｅｓｔ）を設定する。本フローにおいては、ＲＯＩとは、光音響波信号に基づき再構成を行う領域（光学特性値分布を算出する領域）を指す。
ＲＯＩを設定する際には、Ｓ２００で繰り返し測定された各超音波信号から再構成され
た超音波画像のうち、いずれか１つの超音波画像の画像領域を用いればよい。例えば、複数の超音波画像のうち、最初の測定画像や、最後の測定画像、中間の回数での測定画像などを使用する。また、それらの測定画像を表示装置１６０に表示し、作業者がＲＯＩを任意に設定してもよい。この場合、マウスやタッチパネル等の入力機器や、座標を指示する入力方法を使用できる。
【００３８】
（Ｓ５００：光音響波受信時の音響受信器とＲＯＩとの位置を取得する工程）
この工程では、第２の検出手段としての変位量算出モジュール１５１が、Ｓ４００でＲＯＩを設定した超音波画像を取得した際の音響受信器１３０の位置と、Ｓ１００の工程で光音響波を受信した際の音響受信器１３０の位置との間の変位量を算出する。変位量は、並進（音響受信器１３０に沿った方向（Ｘ）、および音響受信器１３０に鉛直下面方向（Ｙ）の移動量）、回転値（θ）などで表現される。このような変位は、作業者がハンドヘルド探触子を手技で移動させることや、三次元的に配置された探触子を機械的に三次元でスキャンすることにより発生する。
【００３９】
変位量算出に当たっては、音響受信器１３０に磁気センサを配置し、磁気の変動から変位量を算出する方法がある。また、音響受信器１３０に光学式のトラックボール等の光学センサを配置し、それに用いる赤外線を測定する方法もある。また、既定の形状に沿ってステッピングモーターなどで音響受信器１３０を制御する場合、制御データを利用して変位量を算出してもよい。
【００４０】
また、Ｓ２００で取得された超音波信号や、それに基づく超音波画像を用いて変位量を算出してもよい。この方法は、外部のセンサや機器類を使用しない点で好ましい。超音波画像には音響インピーダンスやスペックルパターンによる特徴があり、それらを元に画像間の相関を算出できる。
この場合、変位量算出モジュール１５１は、各超音波画像にブロックマッチング法などの既知のアルゴリズムを適用する。ブロックマッチング法では、画像を一定の大きさの画像領域（ブロック）に分割し、各超音波画像において対応するブロックの位置を探索し、ブロックの移動および回転量を算出する。探索する領域の評価関数には、画像の相関値を示す既知の指標Ｓを用いる。探索範囲の中でもっとも相関性が高く算出される動きベクトルを、各ブロックの変位量として算出する。
【００４１】
相関値の指標Ｓとしては、ＳＡＤ（Ｓｕｍ  Ｏｆ  Ａｂｓｏｌｕｔｅ  Ｄｉｆｆｅｒｅｎｃｅ）、ＳＳＤ（Ｓｕｍ  Ｏｆ  Ｓｑｕａｒｅｄ  Ｄｉｆｆｅｒｅｎｃｅ）、ＣＣ（Ｃｒｏｓｓ−Ｃｏｒｒｅｌａｔｉｏｎ）などがある。指標Ｓとしては、また、ＮＣＣ（Ｎｏｒｍａｌｉｚｅｄ  Ｃｒｏｓｓ−Ｃｏｒｒｅｌａｔｉｏｎ）、ＺＮＣＣ（Ｚｅｒｏ−ｍｅａｎ  Ｎｏｒｍａｌｉｚｅｄ  Ｃｒｏｓｓ−Ｃｏｒｒｅｌａｔｉｏｎ）などもある。あるいは、これらの指標Ｓを、近傍のブロックの値により補完、補正してもよい。
【００４２】
例えば、ＳＡＤに関しては、両画像のブロック内の各画素をｆ（ｉ，ｊ）、ｇ（ｉ，ｊ）とすると、指標Ｓは下記の式（１）で示される。
【数１】
2014224806.tif 000003
また、ＳＳＤに関しては、指標Ｓは下記の式（２）で示される。
【数２】
2014224806.tif 000004
また、ＣＣに関しては、指標Ｓは下記の式（３）で示される。
【数３】
2014224806.tif 000005
また、ＮＣＣに関しては、指標Ｓは下記の式（４）で示される。
【数４】
2014224806.tif 000006
また、ＺＮＣＣに関しては、指標Ｓは下記の式（５）で示される。
【数５】
2014224806.tif 000007
【００４３】
各画像領域の相関値を算出する際には、変位量算出モジュール１５１が算出した中で、最も値の高い相関値の指標Ｓを抽出する。なお、領域分割の方法としては、画像領域を等面積の矩形などに分割する方法がある。このとき、画像の輝度値などによって分割範囲や分割数を変更してもよい。
【００４４】
ブロックマッチング法における音響受信器１３０の変位量算出には様々な方法がある。例えば、音響受信器１３０の近傍の矩形領域のみを対象として上記の探索を行い、発見された矩形領域と元の矩形領域の差をそのまま変位量としてもよい。また、音響受信器１３０の周囲の矩形領域の変位量を求め、それを用いて線形補間による近似演算を行い、音響受信器１３０自身の変位量を算出してもよい。あるいは、多数の矩形領域の変位量の平均値を利用してもよい。なお、画面端などで正確に変位が算出できない場合は、最頻値を利用して誤差を軽減してもよい。
【００４５】
なお、音響受信器１３０の位置を超音波画像から求めた場合、超音波信号受信時点での位置が算出される。したがって、光音響画像の補正にかかる位置情報を用いる場合、超音波信号受信のタイミングと光音響波信号受信のタイミングを対応づける必要がある。
かかる対応付けの一例として、超音波信号受信時の音響受信器１３０の位置を、近似的に光音響波受信時の音響受信器１３０の位置としてもよい。さらに、光音響波受信時の音響受信器１３０の位置を、その光音響波受信の前後に行われた超音波受信時の音響受信器
１３０の位置から推定してもよい。
【００４６】
（Ｓ６００：音響波受信器の変位量をもとに光音響波信号を回転・移動する工程）
この工程では、画像再構成モジュール１５２が、Ｓ５００で算出した音響受信器の変位量をもとに、光音響波信号受信時の音響受信器１３０の座標系を変換する。これにより、異なる方向から取得した光音響波に由来する光音響信号を、同じ基準に揃えた上で比較、合成できる。
例えば、Ｓ５００で算出した変位角度をθ、水平方向移動量をｘ’、それに垂直な向きの移動量をｙ’とする。このとき、音響受信器の位置ベクトルｘ０、ｙ０に関して、座標変換後の音響受信器の位置ベクトルｘ１、ｙ１は次の式（６）で示される。
【数６】
2014224806.tif 000008
この座標変換を変位量に基づいてそれぞれ実行し、信号処理装置１５０内部のメモリに格納する。音響受信器の位置や方向に対応する座標系の一例を図９に示す。この例では、鉛直下方を測定するときの音響受信器９２１が取得する信号はｘａ−ｙａ座標系で表現され、斜めに測定するときの音響受信器９２０による取得信号はｘａ−ｙａ座標系で表現される。
なお、Ｓ５００の結果から音響受信器の変位量が微量であり、その誤差を補正する必要がない場合は、Ｓ６００を省略してもよい。
【００４７】
（Ｓ７００：被検体内の変位量を取得する工程）
この工程では、第１の検出手段としての変位量算出モジュール１５１において、Ｓ２００あるいはＳ１００で取得した被検体情報を用いて、異なる画像間での被検体内の変位量を算出する。
被検体内の変位量の算出は、被検体内の画像領域の比較により行われる。まず、画像再構成モジュール１５２が、Ｓ６００の工程によって揃えられた座標系に基づいて被検体の画像再構成を行う。つぎに変位量算出モジュール１５１が、各画像領域についての変位量をアルゴリズムに従って算出し、信号処理装置１５０内部のメモリに格納する。アルゴリズムの一例として、超音波画像の分野で既知の、スペックルパターンを用いたスペックルトラッキング法を利用できる。
各画像領域の変位量は、格子状に分割した画像領域それぞれについて算出すればよい。より好ましくは、ＲＯＩ内で再構成したいピクセルのピッチにあわせた格子サイズで変位量を算出する。
【００４８】
（Ｓ８００：被検体内の変位量に基づき、投影位置座標を決定する工程）
この工程では、画像再構成モジュール１５２が、Ｓ７００で算出した変位量に基づき、光音響波信号データの投影位置座標を決定する。ここで、投影される空間が３次元の場合は投影位置座標を投影位置ボクセルと呼び、２次元の場合は投影位置ピクセルと呼ぶ。
また、光音響波信号データとは、タイムドメインで再構成をする際に、光音響信号を時系列に従って投影位置ボクセルに割り振るデータを指す。その他、光音響波信号データは各トランスデューサが時系列で受信した光音響信号を基準としていればどのような信号でもよい。例えば、光音響波信号データは、Ｓ１００で取得した光音響波信号（ＲＦ信号）や、包絡線検波したＲＦ信号、ＲＦ信号に一定の重みづけを行った信号、ＲＦ信号にある係数をコンボリューションした信号であってもよい。他にも、光音響波信号データは、Ｒ
Ｆ信号に音響受信器の法線ベクトルと音響受信器−投影位置ボクセル間ベクトルとの角度に応じた重みづけを行った信号などであってもよい。
【００４９】
図１０Ａ〜図１０Ｄを用いて投影データについて説明する。
図１０Ａに、音響受信器１０００の１つのトランスデューサと、被検体内を格子状に分割した各画像領域を示す。ここで、音響受信器１０００が光音響信号１０１０を受信した場合、光照射からの時間および被検体内の音速から、音響受信器１０００から等距離に存在する、同心円状のいずれかの画像領域（投影データ領域１０２０と呼ぶ）に音源があると推定される。またこのとき、図１０Ｂに示すように、Ｓ７００で算出された変位量データは、各分割領域における変位ベクトル１０３０で表わされるものとする。このとき、図１０Ｃのように図１０Ｂと同じ大きさで分割された領域に対して、変位ベクトル１０３０を用いた補正を行う。すると、図１０Ｃにおける各投影データ領域１０２１はベクトルに応じて移動し、図１０Ｄにおける変位量を考慮した投影データ領域１０２２となる。
【００５０】
例えばＵＢＰ法（Ｕｎｉｖｅｒｓａｌ  Ｂａｃｋ−Ｐｒｏｊｅｃｔｉｏｎ）であれば初期音圧分布ｐ０（ｒ）は次式（７）により時間空間上で推定される。
【数７】
2014224806.tif 000009
このとき投影データに相当する項ｂ（ｒ０，ｔ）を、式（８）に示す。
【数８】
2014224806.tif 000010
また、任意のデータ領域に対する受信器の立体角の項ｄΩ０は、式（９）となる。
【数９】
2014224806.tif 000011
ここで、ｐ（ｒ）はＳ１００で取得した光音響波信号、ｒは位置、ｔは時間であり、θは受信器と任意のデータ領域とがなす角度である。
また音源の大きさに比べて、音源と測定位置の距離が十分大きい場合、投影データ項ｂ（ｒ０，ｔ）は、式（１０）のようにしてもよい。
【数１０】
2014224806.tif 000012
【００５１】
このときＳ７００で求まった投影位置ボクセルの位置に応じて投影データのｔが変わるため、この工程で決定された投影位置ボクセルに対するｔを用いて投影データ項ｂ（ｒ０，ｔ）を抽出する。
例えばＤＡＳ法（Ｄｅｌａｙ  Ａｎｄ  Ｓｕｍ）もしくは整相加算法であれば初期音圧分布ｐ０（ｒ）は次式（１１）により時間空間上で推定される。
【数１１】
2014224806.tif 000013
ここでτは、次式（１２）で表される、i番目のトランスデューサからのｆ番目の投影
位置ボクセルのｘ座標ｘｆ、ｙ座標ｙｆ、音速ｃ０より求まる、Ｍ個のトランスデューサから投影位置ボクセルまでの遅延時間となる。
【数１２】
2014224806.tif 000014
ωは窓関数などの重みファクターなどでありトランスデューサの指向性や再構成画像の所望の解像度やＳＮに応じて変更される。
このとき被検体内部の変位が考慮された投影データとして、各トランスデューサが受信した時系列の信号から投影位置ボクセルに対して位相がそろうサンプリングデータＳｉを抽出する。
【００５２】
なお、Ｓ７００の結果から被検体内部の変位量が微量であり、その誤差を補正する必要がない場合は、Ｓ８００を省略してもよい。
また、本工程において投影位置座標を補正する代わりに、Ｓ６００と同様に各タイミングの被検体内部の変位量に基づいて各タイミングにおける投影データの座標系を補正してもよい。
また、Ｓ６００において光音響信号の座標系を補正する代わりに、本工程の工程と同様に音響受信器の変位量に基づいて光音響信号が投影される投影位置座標を補正してもよい。
【００５３】
（Ｓ９００：被検体内部の変位量に基づき決定された投影データ画像再構成を行う工程）
この工程では、画像再構成モジュール１５２が、Ｓ８００で抽出した投影データに基づいたすべてのフレームの投影データ信号から、ＵＢＰであれば式（７）、ＤＡＳであれば式（１１）などにしたがってＲＯＩの領域の画像再構成を行う。あるいは被検体内部の変位量を考慮しない場合は、元の光音響波信号から画像再構成を行う。
画像再構成された光音響波画像、超音波画像は、表示装置１６０が画像として表示する。それぞれの画像を重畳させたり並置させたりして単一画面に表示してもよいし、画像を切り替えて表示してもよい。
【００５４】
＜実施例１＞
本実施例では、規定回数受信した光音響波に由来する光音響信号、および、エコー波に由来する超音波信号から、補正した光音響画像を取得する方法を説明する。処理のフロー自体は図２で示したものと同じであり、ステップＳ３００における信号取得回数の設定は１００回とする。
音響受信器１３０には１Ｄリニアアレイ探触子を用いる。また、装置構成として、探触子を走査してデータを取得している最中は探触子走査に関する指示を画面上に表示し、既定のシーケンス終了後はデータ取得完了のメッセージを表示させるようにした。作業者に対する操作指示としては、被検体を複数の方向から測定するために探触子を移動させるようなものが好ましい。
【００５５】
シーケンス開始後、まず、Ｓ１００において光音響波信号を取得する。次に、Ｓ２００において超音波信号を取得する。
このときのＳ１００、Ｓ２００のタイミングについて図３のタイミングチャートを用いて説明する。図３において，光音響波信号は発光をトリガに所定時間（例えば３０μｓｅｃ）取得される。発光の間隔は，光源１１０の発光周波数から決まる。例えば発光周波数が１０Ｈｚとすると、発光の間隔は１００ｍｓｅｃである。１回の発光時間を５０ｎｓｅｃとすると、光音響の信号取得終了から次の発光までの時間は９９ｍｓｅｃ以上ある。その期間に超音波の送信を行い、送信後、所定時間（例えば６０μｓｅｃ）、エコー波を取得する。
【００５６】
本実施例では、Ｓ３００において前記工程を１００回繰り返す間、図４に示すように音響受信器を移動させる。まず、図４（ａ）のように、音響受信器４１０の受信面の法線が、鉛直下向きに対して斜め２０度程度になるように配置する。そして、被検体内の音源４００からの光音響波およびエコー波を取得する。そして徐々に音響受信器を移動させ、角度を変化させていく。途中の図４（ｂ）では、音響受信器４１１が、鉛直下向きにある音源４０１からの光音響波信号を受信する。最後の図４（ｃ）では、図４（ａ）と逆方向に音響受信器４１２を配置して、音源４０２に対して鉛直下向きに対して２０度程度傾いた方向からの信号を取得する。図４（ａ）〜図４（ｃ）にかけて、それぞれの位置で信号が受信され、各位置での座標系に従って情報が保存される。
【００５７】
Ｓ４００におけるＲＯＩの設定に関しては、５０回目に測定した超音波信号に基づく再構成画像から、音響受信器サイズに対応する幅４０ｍｍ、深さ３０ｍｍの領域を抽出する。
Ｓ５００の探触子の変位量算出に関しては、１００回分の超音波信号を画像再構成モジュール１５２で再構成した画像を分割して、各分割領域の移動を追跡する。具体的には、まず、探触子直下にあたる幅３０ｍｍ、深さ１０ｍｍの領域を２ｍｍ四方の領域に分割した。次にＲＯＩ設定に用いた５０回目の画像を基準として、各分割領域について、ＺＮＣＣを相関指標として、最も相関の高い並進量及び回転量を算出する。そして、それらの値の中で最頻値をもって探触子の変位量とする。
Ｓ６００における座標系変換に関しては、Ｓ５００で算出した並進量、回転量を用いて各超音波取得時の音響受信器の座標系を変換し、メモリ上に保管する。
【００５８】
Ｓ７００における被検体内部の変位量算出に関しては、画像再構成モジュール１５２が、Ｓ６００において算出された探触子の座標位置に従って、１００データ分の超音波信号に基づく再構成を行う。そして再構成された画像を２ｍｍ四方の領域に分割して、それぞれＲＯＩに設定した５０回目の画像の各分割領域に対して最も相関の高い並進量、回転量を算出してメモリ上に保管する。このとき、相関指標としてＺＮＣＣを使用する。
Ｓ８００の投影位置座標の補正に関しては、画像再構成モジュール１５２が、Ｓ７００でメモリに保管された被検体内部の変位量を用いて、投影位置座標の補正を行う。
そしてＳ９００において、画像再構成モジュール１５２が、ＲＯＩ内の投影データを積算し、重みづけをおこなうＵＢＰ法による再構成アルゴリズムを適用して画像再構成する。
【００５９】
この手法による視認性向上の効果を図７に示す。これらは柱状の光吸収体を再構成した
ものである。図７（ａ）に、従来の手法による再構成画像を示す。一方、図７（ｂ）には、本実施例の手法による再構成画像を示す。一般に、音響受信器に垂直方向に配置された光吸収体からの音響波の観測は、原理的に難しい。そのため図７（ａ）では、光吸収体の上下端しか画像化できていない。一方、本実施例の手法によれば、光音響波の受信位置（被検体に対する方向）が様々であるため、垂直方向以外から信号を取得できる。そのため、図７（ｂ）に示すように柱状の光吸収体を再現できる。
【００６０】
従来手法の画像および本発明の手法の画像をそれぞれ、実際の光吸収体の配置形状と比較した。すると、従来の手法での再現度が３％であるのに対して、本手法の再現度はほぼ１００％である。すなわち本手法により、光吸収体の形状をより正確に再現できることが分かる。これは、探触子を移動させて光音響波をより多くの方向から受信することで、探触子の開口が広がり、信号データを補完しあえるようになるためである。このように、本発明によれば探触子の受信方向に制約がある場合でも、探触子を移動させ、開口を実質的に広げてアーティファクトやノイズを低減した再構成画像を取得できる。その際、音響受信器の位置に応じて適切な座標系変換等の処理が行われる。また本実施例では、超音波画像による被検体内部の音源の変位を考慮した補正が行われるため、音響受信器の移動や体動などの影響を抑制して良好な画像を取得できる。
【００６１】
＜実施例２＞
実施例２では、光音響波取得時の音響受信器の位置を推定する方法を示す。図５は本実施例の処理フロー図である。なお、実施例１と共通の部分については説明を省略する。
【００６２】
音響受信器１３０は１Ｄリニアアレイ探触子とし、信号の取得回数の設定は１００とする。Ｓ４００工程までの処理は、実施例１と同様である。
次に、Ｓ５０１に示した方法で、超音波画像をもとにＲＯＩを設定した音響受信器に対しての相対変位（相対位置）として、並進量と回転量を算出する。次に、算出した値に基づいて、光音響波を取得した時間に相当する時点での並進量と回転量を線形補間により算出し、光音響波取得時の音響受信器の座標を推定する。
次に、Ｓ６００において、Ｓ５０１で算出した並進量と回転量を用いて座標系を変換し、メモリ上に保管する。
そして、Ｓ９００で、Ｓ６００で算出された座標系に基づき１００データから再構成を行って画像を算出する。
【００６３】
この手法による視認性向上の効果を図８に示す。これらは２ｍｍ径の点状の光吸収体を再構成したものである。図８（ａ）に従来の手法による再構成画像を示す。一方、図８（ｂ）には本実施例の手法によって画像再構成された結果を示す。図８（ａ）では再構成をするための信号情報が十分でないために、光吸収体を示す明るい点の左右に再構成アーティファクトが発生している。一方、本実施例の手法では、音響受信器の取得位置（被検体に対する方向）が様々であるため、垂直方向以外から信号を取得できる。その結果、図８（ｂ）に示すように、再構成アーティファクトを抑制し、光吸収体の形状を再現できている。
【００６４】
従来手法の画像および本発明の手法の画像をそれぞれ、実際の光吸収体の配置形状と比較した。すると、従来の手法である図８（ａ）では、光吸収体のサイズが実物の３５０％であるのに対して、本手法の図８（ｂ）では１１０％である。すなわち本手法により、実際の吸収体の形状をより正確に再現できることが分かる。つまり、従来法のアーティファクトを再現する画像データを仮に積算したとしても、アーティファクトのＳＮが上がるだけで現実の形状とは異なるものが再現される。一方、本手法によれば、アーティファクトを抑えて現実の形状を再現できる。
【００６５】
＜実施例３＞
本実施例では、図１１のフローを参照しつつ、過去に取得したデータを利用してリアルタイムに画像を出力する方法を説明する。
音響受信器は１Ｄリニアアレイ探触子とし、メモリに保存する信号の取得回数の設定は５０とした。
【００６６】
まず、Ｓ１００で、被検体に光を照射して光音響波信号を取得する。
次に、Ｓ２００において、超音波を送受信する。
次に、Ｓ３０００においてＳ１００およびＳ２００で取得した信号を最後に取得したデータから最大で５０回分までをメモリに順次保管する。
次に、Ｓ４０１において、最後に取得した超音波から得られる画像領域をＲＯＩに設定する。
次に、Ｓ５０００において、メモリに保管されている光音響波取得時の音響受信器の位置を、変位量算出モジュール１５１を用いて超音波画像から算出する。
次に、Ｓ７００で被検体内部の変位量算出を行い、Ｓ８００で投影データの投影位置座標の補正をＳ７００で算出された被検体内部の変位量を用いて行う。
そしてＳ９００において、Ｓ５００で作成された座標を元に、画像再構成モジュール１５２が、Ｓ８００で作成されたＲＯＩ内の投影データを重ね合わせる再構成アルゴリズムを適用して画像再構成する。
【００６７】
上記過程を任意の回数繰り返す際にプローブの位置を適時動かすことで、常にメモリに保管されるデータが更新され続けたままリアルタイムに過去データを利用しながら仮想的に探触子数を増やした再構成画像を表示し続けることが可能になる。すなわち本実施例によれば、被検体内の情報をより正確にリアルタイムに再現できる。
【００６８】
＜実施例４＞
本実施例では、図６のフローを参照しつつ、音響受信器の変位量の算出を画像データ取得ごとに行う方法を説明する。この方法であれば、計算の処理が分散され、メモリの負荷を抑制して画像を出力できる。
【００６９】
音響受信器１３０は１Ｄリニアアレイ探触子とし、信号の取得回数の設定は１００とした。また音響受信器１３０には磁気センサを一体化させてある。磁気センサは、光音響波受信時の音響受信器の位置や傾き状況を検知し、信号処理装置１５０に送信した。
【００７０】
まず、Ｓ１４００において、一回目に取得した超音波から得られる画像領域をＲＯＩに設定する。
次に、Ｓ２１００で、被検体に光を照射して光音響波信号を取得する。
次に、Ｓ３２００において、超音波を送受信する。
次に、Ｓ４５０１において、ＲＯＩに設定した超音波取得時の音響受信器の位置を基準として、直前の光音響波取得時の音響受信器の位置からの音響受信器の並進量と回転量を、磁気センサの出力値から求める。
次に、Ｓ５６００では、Ｓ４５０１で算出した並進量と回転量を用いて座標系を変換する。
【００７１】
次に、Ｓ６７０１において、Ｓ５６００で算出した座標系に基づき、超音波の画像再構成を行って画像を算出する。そして再構成画像の比較結果から、ＲＯＩ内の並進（音響受信器に沿った方向（Ｘ）の移動量、および、音響受信器に鉛直下面方向（Ｙ））量、すなわち被検体内部の変位量を算出する。
次に、Ｓ７８００において、光音響波のＲＯＩ内の投影データの投影位置座標をＳ６７０１で算出された値をもとに補正する。
次に、Ｓ８０００では、Ｓ７８００での補正値をメモリに保存する。
上記の処理をＳ９３０２に従って２０回繰り返す。このときメモリに保管する領域を図９に示す。ここでは、図面上斜め方向を向いているときの音響受信器９２０から再構成領域９１０が再構成される。また、変更後の位置にあるときの音響受信器９２１から画像化領域９１１が再構成される。この場合、メモリに保管する領域は共通の画像領域９１２とする。
最後に、Ｓ１０９００にて、２０回分の共通部分のデータを用いて再構成を行う。
【００７２】
本実施例の手法によれば、メモリの使用が抑制可能であり、処理負荷が分散されるため、処理速度を向上させることができ、リアルタイムの処理を好適に実施できる。そして、音響受信器が複数の方向から信号を取得することにより、信号の取れていない箇所を再構成することで生じる再構成アーティファクトを抑制できる。すなわち本実施例によれば、被検体内の情報をより正確にリアルタイムに再現できる。
【００７３】
＜実施例５＞
本実施例では、図１２のフローを参照しつつ、音響受信器としてトランスデューサが三次元配置された探触子を使用する方法を説明する。この方法によれば、三次元の再構成像で視認性を向上させることができる。
【００７４】
トランスデューサが三次元配置された探触子としては図１３に示すように半球の表面上に多数のトランスデューサを配置した音響受信器１３０３を使用した。半球の内部は生体との音響カップリングを保ち、光を透過するジェル材を封入し、半球の下面は平面に設置できる構成とした。半球上面部に光照明部（１３０１）を有し、その隣接部に磁気センサ（１３０４）を一体化して配置した。磁気センサは、光音響波受信時の音響受信器の位置や傾き状況を検知し、信号処理装置１５０に送信した。
またメモリに保存する信号の設定は５０とした。
【００７５】
まず、Ｓ１００で、被検体に光を照射して光音響波信号を取得する。
次に、Ｓ１００１において、Ｓ１００を実施した時点での磁気センサのデータを受信する。
次に、Ｓ３００１においてＳ１００およびＳ１００１で取得した信号を、最後に取得したデータから最大で５０回分までをメモリに順次保管する。
次に、Ｓ４００において、最後に取得した時点での探触子直下の領域２ｃｍ＾３のボリューム領域をＲＯＩに設定する。
次に、Ｓ５０００において、最後に取得した時点での音響受信器の位置座標を基準にして、メモリに保管されている光音響波取得時の音響受信器の３次元位置座標をメモリに保管されている磁気センサデータから算出する。
次に、Ｓ７００で被検体内部の変位量算出を行い、Ｓ８００で投影データの投影位置座標の補正をＳ７００で算出された被検体内部の変位量を用いて行う。
そしてＳ９００において、Ｓ５００で作成された３次元位置座標を元に、画像再構成モジュール１５２が、メモリに保管されている受信信号をもとに投影データを作成し、ＲＯＩ内の投影データを３次元的に重ね合わせる再構成アルゴリズムを適用し再構成する。
【００７６】
上記過程を任意の回数繰り返す際に音響受信器の位置を適時動かすことで、常にメモリに保管されるデータが更新され続けたままリアルタイムに過去データを利用しながら仮想的に音響波受信器の数を増やして３次元再構成像を表示し続けることが可能になる。すなわち本実施例によれば、被検体内の情報をより正確に３次元でリアルタイムに再現できる。
以上、好適な実施形態を説明したが、本発明はこれらの実施形態に限らず、特許請求の範囲を逸脱しない限りにおいて、種々の変形例、応用例も包含するものである。
【符号の説明】
【００７７】
１１０：光源、１３０：音響受信器、１５０：信号処理装置
【図１】
2014224806.tif 000015
【図２】
2014224806.tif 000016
【図３】
2014224806.tif 000017
【図４】
2014224806.tif 000018
【図５】
2014224806.tif 000019
【図６】
2014224806.tif 000020
【図７】
2014224806.tif 000021
【図８】
2014224806.tif 000022
【図９】
2014224806.tif 000023
【図１０Ａ】
2014224806.tif 000024
【図１０Ｂ】
2014224806.tif 000025
【図１０Ｃ】
2014224806.tif 000026
【図１０Ｄ】
2014224806.tif 000027
【図１１】
2014224806.tif 000028
【図１２】
2014224806.tif 000029
【図１３】
2014224806.tif 000030