(57)【特許請求の範囲】
【請求項１】
光を照射された測定対象から伝搬する音響波を受信して電気信号に変換する複数の振動子と、
前記複数の振動子の少なくとも一部の振動子の指向軸が集まるように前記複数の振動子が配置された探触子と、
前記複数の振動子の位置に関する位置情報を取得する位置情報取得部と、
前記電気信号と、前記位置情報とに基づいて、前記測定対象の特性情報を取得する特性情報取得部と、
を有し、
前記位置情報取得部は、
点音源が、当該点音源と前記探触子が所定の相対位置となるような第一の位置に設置された場合の、前記第一の位置と前記複数の振動子のそれぞれの間の距離である第一のデータ群を取得し、
前記測定対象が前記点音源であるときに、前記点音源から実際に伝搬した前記音響波に由来する前記電気信号に基づいて、前記点音源と前記複数の振動子のそれぞれの間の距離である第二のデータ群を算出し、
前記第一のデータ群および前記第二のデータ群に基づいて前記位置情報を算出する
ことを特徴とする被検体情報取得装置。
【請求項２】
前記特性情報取得部は、前記複数の振動子それぞれが出力した前記電気信号を、前記複数の振動子それぞれの位置の設計値からのずれに基づいて校正する
ことを特徴とする請求項１に記載の被検体情報取得装置。
【請求項３】
前記位置情報取得部は、前記位置情報に基づいて前記点音源を前記第一の位置に移動させ、前記第一の位置にある前記点音源に関する前記第二のデータ群を取得し、
前記特性情報取得部は、前記第一の位置にある前記点音源に関する前記第二のデータ群を前記校正に用いる
ことを特徴とする請求項２に記載の被検体情報取得装置。
【請求項４】
前記点音源を支持する治具をさらに有し、
前記位置情報取得部は、前記治具と前記探触子の相対位置を変化させることにより前記点音源を移動させる
ことを特徴とする請求項３に記載の被検体情報取得装置。
【請求項５】
前記特性情報取得部は、前記測定対象の単位領域ごとに、前記複数の振動子から出力された前記電気信号から、前記音響波の媒質の音速、および、前記複数の振動子と前記単位領域の距離に基づいて、前記単位領域に対応するデータを選択して加算することで前記単位領域の特性情報を取得するものであり、
前記複数の振動子と前記単位領域の距離は、前記第二のデータ群に基づいて校正されたものである
ことを特徴とする請求項２ないし４のいずれか１項に記載の被検体情報取得装置。
【請求項６】
前記第一の位置は、前記複数の振動子の指向軸が集まる位置である
ことを特徴とする請求項１ないし５のいずれか１項に記載の被検体情報取得装置。
【請求項７】
前記探触子は、内周面に前記複数の振動子が配置された半球状の部材であり、
前記第一の位置は、前記半球状の探触子の曲率中心点である
ことを特徴とする請求項６に記載の被検体情報取得装置。
【請求項８】
前記点音源は、光を照射されて音響波を発生させる球状の部材、または、音響波を送信する音響トランスデューサである
ことを特徴とする請求項１ないし７のいずれか１項に記載の被検体情報取得装置。
【請求項９】
光を照射された測定対象から伝搬する音響波を受信して電気信号に変換する複数の振動子と、
前記複数の振動子の少なくとも一部の振動子の指向軸が集まるように前記複数の振動子が配置された探触子と、
前記複数の振動子の位置に関する位置情報を取得する位置情報取得部と、
特性情報取得部と、
を有する被検体情報取得装置の制御方法であって、
前記位置情報取得部が、点音源が、当該点音源と前記探触子が所定の相対位置となるような第一の位置に設置された場合の、前記第一の位置と前記複数の振動子のそれぞれの間の距離である第一のデータ群を取得するステップと、
前記位置情報取得部が、前記測定対象が前記点音源であるときに、前記点音源から実際に伝搬した前記音響波に由来する前記電気信号に基づいて、前記点音源と前記複数の振動子のそれぞれの間の距離である第二のデータ群を算出するステップと、
前記位置情報取得部が、前記第一のデータ群および前記第二のデータ群に基づいて前記位置情報を算出するステップと、
前記特性情報取得部が、前記電気信号と、前記位置情報とに基づいて、前記測定対象の特性情報を取得するステップと、
を有することを特徴とする被検体情報取得装置の制御方法。
【請求項１０】
前記特性情報取得部が、前記複数の振動子それぞれが出力した前記電気信号を、前記複数の振動子それぞれの位置の設計値からのずれに基づいて校正するステップ
をさらに有することを特徴とする請求項９に記載の被検体情報取得装置の制御方法。
【請求項１１】
前記位置情報取得部が、前記位置情報に基づいて前記点音源を前記第一の位置に移動させ、前記第一の位置にある前記点音源に関する前記第二のデータ群を取得するステップと、
前記特性情報取得部が、前記第一の位置にある前記点音源に関する前記第二のデータ群を前記校正に用いるステップと、
をさらに有することを特徴とする請求項１０に記載の被検体情報取得装置の制御方法。
【請求項１２】
前記被検体情報取得装置は、前記点音源を支持する治具をさらに有し、
前記位置情報取得部が、前記治具と前記探触子の相対位置を変化させることにより前記点音源を移動させるステップ
をさらに有することを特徴とする請求項１１に記載の被検体情報取得装置の制御方法。
【請求項１３】
前記特性情報取得部は、前記測定対象の単位領域ごとに、前記複数の振動子から出力された前記電気信号から、前記音響波の媒質の音速、および、前記複数の振動子と前記単位領域の距離に基づいて、前記単位領域に対応するデータを選択して加算することで前記単位領域の特性情報を取得するものであり、
前記複数の振動子と前記単位領域の距離は、前記第二のデータ群に基づいて校正されたものである
ことを特徴とする請求項１０ないし１２のいずれか１項に記載の被検体情報取得装置の制御方法。
【請求項１４】
前記第一の位置は、前記複数の振動子の指向軸が集まる位置である
ことを特徴とする請求項９ないし１３のいずれか１項に記載の被検体情報取得装置の制御方法。
【請求項１５】
前記探触子は、内周面に前記複数の振動子が配置された半球状の部材であり、
前記第一の位置は、前記半球状の探触子の曲率中心点である
ことを特徴とする請求項１４に記載の被検体情報取得装置の制御方法。
【請求項１６】
前記点音源は、光を照射されて音響波を発生させる球状の部材、または、音響波を送信する音響トランスデューサである
ことを特徴とする請求項９ないし１５のいずれか１項に記載の被検体情報取得装置の制御方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被検体情報取得装置およびその制御方法に関する。
【背景技術】
【０００２】
近年、光音響イメージングを利用して被検体内を画像化する研究が行われている。光音響イメージングの１つにＰｈｏｔｏａｃｏｕｓｔｉｃ  Ｔｏｍｏｇｒａｐｈｙ（ＰＡＴ：光音響トモグラフィー）がある。
【０００３】
光源から発生したパルス光を被検体に照射すると、光エネルギーを吸収した被検体内部の組織から音響波が発生する。この現象を光音響効果と呼ぶ。この発生した音響波を被検体の周囲に複数配置された振動子によって検出する。そして、その受信信号を信号処理することで被検体内の情報を取得できる。これが光音響トモグラフィーによるイメージングの原理である。
【０００４】
例えば、ヘモグロビンの吸収が大きい近赤外光をパルス光に用いる場合、被検体内のヘモグロビンすなわち血液の存在する個所をイメージングできる。血管イメージング結果の悪性腫瘍の診断などへの利用が期待されている。
【０００５】
光音響トモグラフィーによるイメージングの例として、特許文献１には、光の照射領域と音響波を受信する振動子を移動させながら被検体からの音響波を受信し、被検体情報を再構成する手法が示されている。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】米国特許第５８４００２３号明細書
【発明の概要】
【発明が解決しようとする課題】
【０００７】
複数の振動子で受信された受信信号を用いて画像再構成をする際に、注目位置と各振動子との距離を基に処理が行われる。しかし、再構成画像の解像度に影響しないように複数の振動子を探触子上に精度よく機械的に配置することは難しい。そのため、機械的な配置精度による誤差で、再構成画像の解像度が劣化するという問題がある。
【０００８】
本発明は上記課題に鑑みてなされたものである。本発明の目的は、複数の振動子の位置を、振動子の受信信号を用いて精度よく校正することにある。
【課題を解決するための手段】
【０００９】
本発明は、以下の構成を採用する。すなわち、
光を照射された測定対象から伝搬する音響波を受信して電気信号に変換する複数の振動子と、
前記複数の振動子の少なくとも一部の振動子の指向軸が集まるように前記複数の振動子が配置された探触子と、
前記複数の振動子の位置に関する位置情報を取得する位置情報取得部と、
前記電気信号と、前記位置情報とに基づいて、前記測定対象の特性情報を取得する特性情報取得部と、
を有し、
前記位置情報取得部は、
点音源が、当該点音源と前記探触子が所定の相対位置となるような第一の位置に設置された場合の、前記第一の位置と前記複数の振動子のそれぞれの間の距離である第一のデータ群を取得し、
前記測定対象が前記点音源であるときに、前記点音源から実際に伝搬した前記音響波に由来する前記電気信号に基づいて、前記点音源と前記複数の振動子のそれぞれの間の距離である第二のデータ群を算出し、
前記第一のデータ群および前記第二のデータ群に基づいて前記位置情報を算出する
ことを特徴とする被検体情報取得装置である。
【００１０】
本発明はまた、以下の構成を採用する。すなわち、
光を照射された測定対象から伝搬する音響波を受信して電気信号に変換する複数の振動子と、
前記複数の振動子の少なくとも一部の振動子の指向軸が集まるように前記複数の振動子が配置された探触子と、
前記複数の振動子の位置に関する位置情報を取得する位置情報取得部と、
特性情報取得部と、
を有する被検体情報取得装置の制御方法であって、
前記位置情報取得部が、点音源が、当該点音源と前記探触子が所定の相対位置となるような第一の位置に設置された場合の、前記第一の位置と前記複数の振動子のそれぞれの間の距離である第一のデータ群を取得するステップと、
前記位置情報取得部が、前記測定対象が前記点音源であるときに、前記点音源から実際に伝搬した前記音響波に由来する前記電気信号に基づいて、前記点音源と前記複数の振動子のそれぞれの間の距離である第二のデータ群を算出するステップと、
前記位置情報取得部が、前記第一のデータ群および前記第二のデータ群に基づいて前記位置情報を算出するステップと、
前記特性情報取得部が、前記電気信号と、前記位置情報とに基づいて、前記測定対象の特性情報を取得するステップと、
を有することを特徴とする被検体情報取得装置の制御方法である。
【発明の効果】
【００１１】
本発明によれば、複数の振動子の位置を、振動子の受信信号を用いて精度よく校正することができる。
【図面の簡単な説明】
【００１２】
【図１】被検体情報取得装置を示す模式図
【図２】被検体情報取得装置における被検体情報取得のフロー図
【図３】校正データ取得のフロー図
【図４】校正データ取得時における仮想点音源と複数の振動子の概略図
【図５】校正データ取得時における点音源と複数の振動子の概略図
【図６】半球状の探触子の模式図
【図７】校正前と校正後の再構成画像の比較
【発明を実施するための形態】
【００１３】
以下に図面を参照しつつ、本発明の好適な実施の形態について説明する。ただし、以下に記載されている構成部品の寸法、材質、形状およびそれらの相対配置などは、発明が適用される装置の構成や各種条件により適宜変更されるべきものである。よって、この発明の範囲を以下の記載に限定する趣旨のものではない。
【００１４】
本発明は、被検体から伝搬する音響波を検出し、被検体内部の特性情報を生成し、取得する技術に関する。よって本発明は、被検体情報取得装置またはその制御方法、あるいは
被検体情報取得方法や信号処理方法として捉えられる。本発明はまた、これらの方法をＣＰＵやメモリ等のハードウェア資源を備える情報処理装置に実行させるプログラムや、そのプログラムを格納した記憶媒体としても捉えられる。
【００１５】
本発明の被検体情報取得装置には、被検体に光（電磁波）を照射することにより被検体内で発生した音響波を受信して、被検体の特性情報を画像データとして取得する光音響効果を利用した装置を含む。この場合、特性情報とは、光音響波を受信することにより得られる受信信号を用いて生成される、被検体内の複数位置のそれぞれに対応する特性値の情報である。
【００１６】
光音響測定により取得される特性情報は、光エネルギーの吸収率を反映した値である。例えば、光照射によって生じた音響波の発生源、被検体内の初期音圧、あるいは初期音圧から導かれる光エネルギー吸収密度や吸収係数、組織を構成する物質の濃度や量を含む。また、物質濃度として酸化ヘモグロビン濃度と還元ヘモグロビン濃度を求めることにより、酸素飽和度分布を算出できる。また、総ヘモグロビン濃度、グルコース濃度、コラーゲン濃度、メラニン濃度、脂肪や水の体積分率なども求められる。
【００１７】
被検体内の各位置の特性情報に基づいて、二次元または三次元の特性情報分布が得られる。分布データは画像データとして生成され得る。特性情報は、数値データとしてではなく、被検体内の各位置の分布情報として求めてもよい。すなわち、初期音圧分布、エネルギー吸収密度分布、吸収係数分布や酸素飽和度分布などの分布情報である。
【００１８】
本発明でいう音響波とは、典型的には超音波であり、音波、音響波と呼ばれる弾性波を含む。探触子等により音響波から変換された電気信号を音響信号とも呼ぶ。ただし、本明細書における超音波または音響波という記載は、それらの弾性波の波長を限定する意図ではない。光音響効果により発生した音響波は、光音響波または光超音波と呼ばれる。光音響波に由来する電気信号を光音響信号とも呼ぶ。
【００１９】
以下の記載や図面において、原則として同一の構成要素には同一の符号を付し、細かい説明は省略する。以下の記載においては、被検体情報取得装置の一例として、被検体内部の特性情報を取得し、画像データを生成可能な光音響装置について説明する。
【００２０】
（装置構成）
図１は本実施形態における被検体情報取得装置の装置構成概略図である。本実施形態の被検体情報取得装置は、探触子１０２、音響波を受信する複数の振動子１０３、信号受信部１０４、信号処理部１０５、位置制御部１０６、光源１０７、システム制御部１０９、特性情報取得部１１０、表示部１１１を備える。
【００２１】
被検体１０１は測定対象である。本発明の被検体情報取得装置は生体の悪性腫瘍や血管疾患を診断することを目的としている。したがって、測定対象としては人の乳房などの生体が想定される。また、装置性能の確認などでは生体特性を模擬したファントムも測定対象として想定される。生体の特性とは音響特性と光学特性のことを表わす。さらに後述する振動子位置の校正の際に用いられる点音源も測定対象に含まれる。
【００２２】
被検体１０１は保持部１１２によって保持されている。保持部１１２は取り付け部１１３に取り付けられている。保持部１１２は測定対象によってサイズや形状を変更してもよい。また保持部１１２を取り付けないで測定することも可能である。保持部１１２は光や音響波を透過する特性を持つことが望ましい。被検体１０１が乳房の場合、伏臥位の被検者を支持する寝台に乳房を挿入するための開口を設けて、開口から下方に垂らされた乳房を測定するとよい。
【００２３】
探触子１０２は、半球状の支持体の内周面に複数の振動子１０３が配置された構成である。振動子１０３の例としては、圧電現象を利用した圧電素子などの変換素子、ＣＭＵＴなどの静電容量型の変換素子がある。ただし変換素子の方式に制限はない。保持部１１２と被検体の間や、保持部１１２と探触子１０２の間には、音響インピーダンスを整合させるための音響整合材を配置することが好ましい。音響整合材としては水、超音波ジェル、ひまし油などがある。
【００２４】
複数の振動子１０３の受信感度が高い方向が集まる領域を高解像度化できる。このような領域を本明細書では高感度領域と呼ぶ。振動子１０３は、一般的に、受信面の法線方向で最も感度が高い。この方向を指向軸と呼ぶ。この指向軸が集まる領域が高感度領域であり、図１のような半球状の探触子では、半球の曲率中心周辺が高感度領域となる。
【００２５】
また、探触子１０２の形状は半球状に限定されない。全ての振動子１０３の指向軸が平行にはならない場合（すなわち、複数の振動子１０３の指向軸の少なくとも一部が集まる場合）、高感度領域が形成される。したがって、探触子１０２は、球冠状、球帯状、楕円体の一部、または、平面または曲面の組み合わせ、などの形状でもよい。
【００２６】
光源１０７はパルス光を発生させる。光源１０７はシステム制御部１０９からの制御信号を受けてパルス光を発生する。光音響効果によって音響波を発生させるには、パルス幅１００ｎｓｅｃ程度のパルスが用いられる。光の波長としては６００〜１０００ｎｍ程度が望ましい。レーザの種類としては、Ｎｄ：ＹＡＧレーザやアレクサンドライトレーザ、Ｔｉ：ｓａレーザなどが用いられる。また、半導体レーザでもよい。また、光源１０７としてフラッシュランプや発光ダイオードを用いてもよい。波長可変レーザを用いたり、波長が互いに異なる複数のレーザを用いたりすることで、物質ごとの波長吸収スペクトルの違いを利用した測定（例えば酸素飽和度測定）が可能になる。
【００２７】
光源１０７で発生させたレーザ光は光学系１０８を伝送して被検体１０１に照射される。光学系１０８としてはレンズやプリズム、ミラー、光ファイバなどが用いられる。
【００２８】
信号受信部１０４は複数の振動子１０３の受信信号を増幅する信号増幅機、アナログの電気信号をデジタル電気信号に変換するＡＤ変換機を含む。生成されたデジタル電気信号は特性情報取得部１１０へ出力される。信号受信部１０４は、光学系１０８からの同期信号をトリガーに動作を開始する。測定対象が光音響波を発生させるものである場合、光源制御に関する同期信号や、光センサが照射光を検知して出力する信号をトリガーにする。
【００２９】
信号処理部１０５は受信信号から振動子位置の校正処理を行い、その結果の校正データを特性情報取得部１１０へ出力する。もしくは保持している校正データを特性情報取得部１１０へ出力する。信号処理部１０５の処理フローについては、後の（校正データ取得のフロー）で詳細を述べる。信号処理部１０５や、後述する位置制御部１０６、特性情報算出部１１０などは、処理回路や情報処理装置で構成できる。情報処理装置としては、ＣＰＵやメモリ等の演算資源を備え、プログラムの指示に従って動作するＰＣやワークステーションが好適である。各ブロックは、同一の情報処理装置で動作するモジュールとして構成されてもよいし、物理的に別々に構成されてもよい。信号処理部は、本発明の位置情報取得部として機能する。
【００３０】
位置制御部１０６は探触子１０２を移動させることにより、複数の振動子１０３と被検体１０１の相対位置を変化させる。これにより、複数の振動子１０３が形成する高感度領域が被検体内部で移動する。その結果、取得した被検体情報画像内の感度のばらつきが軽減する。この探触子１０２の移動は、２次元方向でもよいし、３次元方向でもよい。探触
子１０２を移動させる構成として例えば、ボールねじやアクチュエータを備え、プログラムされた軌跡で移動するＸＹステージを利用できる。
【００３１】
システム制御部１０９は、被検体情報取得装置に含まれる各ブロックと情報をやり取りして、それぞれの動作タイミングや動作内容を統合的に制御する。
【００３２】
特性情報取得部１１０は、信号受信部１０４から出力された受信信号を用いて画像再構成を行い、特性情報を算出する。画像再構成の手法としてはＵｎｉｖｅｒｓａｌ  Ｂａｃｋ  Ｐｒｏｊｅｃｔｉｏｎ（ＵＢＰ）やＦｉｌｔｅｒｅｄ  Ｂａｃｋ  Ｐｒｏｊｅｃｔｉｏｎ（ＦＢＰ）、整相加算などの公知の再構成手法が用いられる。
【００３３】
表示部１１１は、特性情報取得部１１０によって生成された被検体情報を表示する。また、操作者が装置を操作するのに必要なＵＩも表示される。表示部１１１としては液晶ディスプレイ、プラズマディスプレイなど、任意の表示装置を利用できる。表示部１１１は、被検体情報取得装置と一体に提供されてもよいし、別体でもよい。
【００３４】
（点音源）
本発明の点音源は、所定の位置に設置され、任意のタイミングで音響波を発生させる。良好な校正のためには、等方的に音響波を送信できる部材が好ましい。例えば、光の照射を受けて光音響波を発生させる球状の部材を利用できる。材質としては、球の加工精度（真球度）を高めやすく、かつ音響波発生効率が高いものが好ましい。例えば金属球の表面を塗料で黒くコーティングした部材が好適である。他に、樹脂、ゴム、カーボンなども利用できる。なお、音響波を等方的に送信できるのであれば、点音源として音響トランスデューサを用いてもよい。
【００３５】
点音源を所定の位置に設置する際には、点音源を糸や針金でぶら下げて支持する治具を用いるとよい。治具または探触子の少なくともいずれか一方を移動させることで、点音源と探触子の相対位置を変化させられる。光音響効果を利用する点音源であれば、光を照射することで音響波が発生する。
【００３６】
（被検体情報取得のフロー）
以下、図２を参照して、被検体情報取得のフローを説明する。
ステップＳ２０１では、操作者が被検体情報の取得のために必要なレーザや探触子位置制御などのパラメータを設定する。
ステップＳ２０２では、ステップＳ２０１で設定された探触子位置制御に関するパラメータを基に、位置制御部１０６が探触子を指定の位置へ移動させる。複数位置での撮像を設定した場合には１つめの指定位置へ移動する。
【００３７】
ステップＳ２０３では、ステップＳ２０１で設定されたレーザに関するパラメータを基に、パルス光を照射する。パルス光は光学系１０８を伝送して被検体１０１へ照射され、被検体１０１から音響波が発生する。光学系１０８はパルス光の伝送と同時に信号受信部１０４へ同期信号を送信する。これにより信号受信部１０４は受信動作を開始するため、被検体１０１からの音響波に基づく電気信号が受信される。受信された音響波に由来するアナログ電気信号は、信号増幅機とＡＤ変換機により、増幅されたデジタル電気信号へと変換され、特性情報取得部１１０へ出力される。
【００３８】
ステップＳ２０４では、所定の範囲の画像データを生成するために必要な撮像がすべて完了したかどうかの判断がなされる。所定の範囲は、ユーザーの指定や予め設定された値により決定される。撮像が完了していない場合は、次の指定位置への移動を行い、音響波の取得を繰り返す。
【００３９】
ステップＳ２０５では、特性情報取得部１１０が、取得した受信信号、レーザや探触子位置制御の情報から画像再構成を行い、被検体情報を示す画像データを生成する。典型的な画像再構成においては、被検体内部の単位領域ごとに各振動子の受信信号が合計され、初期音圧が求められる。このとき媒体の音速と、単位領域から振動子までの距離に基づき、受信信号中から適切な検出時間のデジタル信号が選択される。よって、単位領域から振動子までの距離が正確であればあるほど、画像再構成の精度が上がる。逆に、振動子の位置が、製造時の精度、各振動子の電気的特性によって生じる信号の遅れ、経時変化、などによって設計値からずれている場合、校正データを用いて適切なデータを選択することが必要となる。そして各単位領域の初期音圧に基づいて初期音圧分布が求められる。初期音圧分布、吸収係数分布、酸素飽和度分布などの形で生成された被検体情報は、表示部１１１に表示される。
【００４０】
（点音源を用いる校正データ取得について）
振動子位置の校正は、装置の組み立て時や定期的なメンテナンス時に行うことが望ましい。本発明で校正とは、各振動子１０３の設計値からのズレを測定し、校正データとして保存しておき、画像再構成時の補正に用いることを言う。校正は例えば、装置の製造時、出荷時、定期点検時に行われる。
【００４１】
校正データ取得の際は、まず初めに、点音源を半球状の探触子１０２の特異相対位置に設置する。特異相対位置とは、あらかじめ複数の振動子１０３の極座標の偏角が分かっている位置や、複数の振動子１０３の指向軸が集まる位置である。半球状の探触子１０２の場合、特異相対位置の１つとして曲率中心点がある。
【００４２】
点音源から音響波が発生すると、複数の振動子１０３が、媒質内を伝搬する音響波を受信して受信信号を出力する。そして信号処理部１１０は、受信信号から点音源に由来する成分を検出することで、点音源から複数の振動子１０３それぞれまでの距離を算出する。距離算出方法としては、時系列的な受信信号の強度の立ち上がりを検出する方法や、所定の閾値を越えた強度を検出する方法がある。音響波の発生から点音源由来成分の検出までに経過した時間と、音響波の媒質（例えば音響整合材）の音速とを用いて距離を算出できる。
【００４３】
この各振動子の距離情報に基づいて校正データが取得できる。校正データは、受信信号に基づいて画像再構成を行うときに、振動子ごとに用いられる。校正データを用いることにより、再構成対象の単位領域（ピクセルまたはボクセル）と振動子の距離が正確に取得できるので、時系列のデジタル信号から最適なデータを選択可能になる。
【００４４】
校正データは、任意の形式でメモリ等に保持する。例えば、探触子がホームポジションにあるときの各振動子の設計上の位置情報が、装置に設定されたＸＹＺ座標系または極座標系における座標値として保持されているものとする。この場合、校正データは、設計値からのズレ量として保持される。またはメモリ上の座標値を上書きしてもよい。または校正日時と合わせてズレ量または座標値をバージョン管理してもよい。あるいは、実測された各振動子の、座標系における座標値を保持してもよい。
【００４５】
しかし、点音源を探触子１０２の曲率中心点に配置することは容易ではない。点音源の設置方法としては治具を用いるが、設置の精度を知ることは難しい。その結果、再構成画像の解像度低下につながる。そこで本発明では、上記のように複数の振動子１０３の受信信号を用いて点音源の位置を取得する。位置制御部１０６は、その位置情報を基に点音源と探触子１０２の相対位置を変化させ、曲率中心点まで移動させる。これにより点音源が望ましい位置に設置できる。
【００４６】
（振動子位置の校正データ取得のフロー）
以下に、図３を参照して、振動子位置の校正データ取得のフローを説明する。ここで、図４は、仮想点音源４０１が半球状の探触子１０２の曲率中心位置付近のある位置（ｘ、ｙ、ｚ）に配置されたと仮定したシミュレーションの様子を示す。
【００４７】
ステップＳ３０１では、図４のように仮定した場合の、仮想点音源４０１と複数の振動子１０３の距離Ｒｎ（ｘ、ｙ、ｚ）を計算する。ｎは１からＮまでとし、Ｎは複数の振動子１０３の数を表す。この距離を算出する際の複数の振動子１０３の位置には、機械設計値（すなわち校正前の振動子位置情報）を用いる。この（ｘ、ｙ、ｚ）を曲率中心位置付近で変更し、それぞれの距離Ｒｎ（ｘ、ｙ、ｚ）を計算する。複数位置（ｘ、ｙ、ｚ）のピッチを細かくし、範囲を広くすることで点音源位置の推定精度が向上する。
【００４８】
この距離Ｒｎ（ｘ、ｙ、ｚ）は計算して保持しておいてもよいし、後述するステップＳ３０４の計算を行う際に随時計算してもよい。各振動子と仮想点音源との距離は、本発明の第一のデータ群に相当する。第一のデータ群は、校正に先立って予めメモリに格納しておいてもよい。また、探触子の経時的な変形の影響を低減するための校正の場合は、前回の校正処理により取得された各振動子の位置に基づいて、新たな第一のデータ群を取得してもよい。
【００４９】
ステップＳ３０２では、図５に示されるように、実際の被検体情報取得装置において、点音源５０１を半球状の探触子１０２の曲率中心点付近に設置する。
【００５０】
ステップＳ３０３では、実際の点音源５０１から音響波を伝搬させる。図５に示されるように、複数の振動子１０３のそれぞれが音響波を受信して受信信号を出力する。特性情報取得部１１０は、受信信号を用いて点音源５０１から複数の振動子１０３までの距離ｒｎを算出する。ｎは１からＮまでとし、Ｎは複数の振動子１０３の数を表す。ここでは、受信信号の立ち上がりを検出することで、距離ｒｎを算出する。この際、点音源５０１の球径を考慮することが好ましい。また、複数の振動子１０３や信号受信部１０４などの受信に関する遅延量も考慮する。さらに、算出距離精度向上のために受信信号を補間処理してもよい。実際の点音源と各振動子の距離を示す情報は、本発明の第二のデータ群に相当する。
【００５１】
ステップＳ３０４では、特性情報取得部１１０は、仮想点音源４０１と各振動子１０３の間の距離Ｒｎ（ｘ、ｙ、ｚ）と、点音源５０１から各振動子１０３までの距離ｒｎの差の２乗が最小となる（ｘ、ｙ、ｚ）を算出する。そして、その位置を点音源５０１の設置位置の推定結果とする。すなわち以下の式（１）のｄ（ｘ、ｙ、ｚ）が最小となる（ｘ、ｙ、ｚ）を算出する。
【数１】
0006598667.tif 000002
【００５２】
ステップＳ３０５では、推定された点音源５０１の位置を基に、点音源５０１が曲率中心点から所定の誤差内に配置されているかを判定する。点音源５０１が中心から所定の誤差以上離れた位置にある場合は、ステップＳ３０６に移行する。誤差が所定よりも小さい場合はステップＳ３０７に移行する。
【００５３】
ステップＳ３０６では、図５に示されるように、ステップＳ３０４で推定された点音源
５０１の位置から曲率中心点までの誤差分だけ、点音源５０１と探触子１０２の相対位置を変化させる。この相対位置の変化は点音源５０１を移動させてもよいし、探触子１０２の位置を移動させてもよい。
【００５４】
ステップＳ３０７では、点音源５０１から複数の振動子１０３までの距離ｒｎを、探触子１０２の曲率中心点から複数の振動子１０３までの距離として振動子位置の校正データを作成し、保持する。
【００５５】
この校正データ取得に用いた点音源は、本発明の第一のデータ群および第二のデータ群を用いて算出された位置情報に基づいて、適切に配置されたものである。したがって、上記手法で得られた校正データを用いて振動子位置の誤差を補正し、画像再構成を行うことで、高解像度な被検体情報を生成できる。また、このように複数の振動子１０３の受信信号を用いて振動子位置を校正することで、振動子の物理的な位置だけでなく、受信回路の特性のばらつきなども補正することができる。
【００５６】
（実施例）
以下、より詳細な実施例を説明する。探触子１０２は、５１２個の複数の振動子１０３が半球状支持体の内周面に配置された部材である。複数の振動子１０３はφ１．５ｍｍの振動子とする。図６は、探触子の平面模式図である。複数の振動子１０３は半球上で３次元スパイラルを形成する。振動子配置の座標系は、曲率中心点を原点とした極座標系の半径ｒ、極角θ、方位角φと、直交座標系ｘ、ｙ、ｚで定義されている。
【００５７】
ここで、各振動子と曲率中心点との距離は、機械設定値から±０．１ｍｍの範囲でランダムに誤差を持っているものとする。したがって、極座標系の半径ｒの方向にも、±０．１ｍｍの範囲の誤差を持っている。これは、各振動子と曲率中心点の間の距離が、設計値から±０．１ｍｍの誤差を持つことに相当する。なお探触子１０２が走査可能な場合、探触子１０２が寝台の開口の直下にある場合をホームポジションとして、この位置で点音源を用いた校正を行う。
【００５８】
（比較例）
本実施例では、φ０．１ｍｍの球形をした測定対象を曲率中心点付近に設置した状態で、光出射部６０１からパルスを照射する。信号受信部のサンプリング周波数は４０ＭＨｚとする。各振動子が受信した光音響波の時間軸は、測定対象と各振動子間の媒質の音速から距離に換算される。そして特性情報取得部が、振動子配置の機械設計値（すなわち校正前の各振動子の位置座標）を基に、受信した音響波を逆投影するＵＢＰ法によって初期音圧分布を再構成する。
【００５９】
図７（ａ）は、再構成された点音源像である。点音源像の半値幅は０．２５ｍｍとなる。この結果には、各振動子と曲率中心点との距離が機械設定値から±０．１ｍｍの誤差を持っていることに起因する解像度の低下が含まれている。
【００６０】
（本発明の手法と効果）
この解像度の低下を、振動子（探触子）の位置情報に基づいて校正データを作成し、画像再構成に用いることで改善する。振動子位置の校正にはφ１．５ｍｍの点音源を用いる。そして、上記の振動子位置の校正データ取得のフローに従って、点音源を半球状の探触子の曲率中心点に移動させる。この際の信号受信部のサンプリング周波数も、上の比較例と同様に、４０ＭＨｚとする。また、Ｒｎ（ｘ、ｙ、ｚ）を算出する際の（ｘ、ｙ、ｚ）のピッチは０．０１ｍｍとする。
【００６１】
以上の条件で、各振動子の受信信号から曲率中心点と各振動子間の距離を算出する。具
体的には、受信信号の立ち上がり位置から取得できる伝搬時間と、点音源と各振動子間の媒質の音速から、曲率中心点と想定される位置から各振動子までの距離を算出する。距離算出の際には、点音源の球径、各振動子の応答の遅延特性、信号受信部のサンプリングに関するオフセット量や遅延量を考慮することが好ましい。これにより距離算出の精度が向上する。
【００６２】
この受信信号から算出した距離を各振動子の極座標系の半径ｒとし、極角θ、方位角φから直交座標系ｘ、ｙ、ｚも修正する。これにより振動子位置の校正データが取得できる。校正データを設計値からの差分として保存してもよい。
【００６３】
この振動子位置の校正データを用いて、比較例と同じφ０．１ｍｍの球形をした測定対象の初期音圧分布をＵＢＰ法にて再構成する。図７（ｂ）は、算出された点音源像である。この点音源像の半値幅は０．１９ｍｍとなる。この結果から解像度が改善し、校正の効果があることが確認された。
【００６４】
（その他の実施例）
本発明は、上述の実施形態の１以上の機能を実現するプログラムを、ネットワーク又は記憶媒体を介してシステム又は装置に供給し、そのシステム又は装置のコンピュータにおける１つ以上のプロセッサがプログラムを読出し実行する処理でも実現可能である。また、１以上の機能を実現する回路（例えば、ＡＳＩＣ）によっても実現可能である。
【符号の説明】
【００６５】
１０２：探触子、１０３：振動子、１０５：信号処理部、１１０：特性情報取得部
【図１】
0006598667.tif 000003
【図２】
0006598667.tif 000004
【図３】
0006598667.tif 000005
【図４】
0006598667.tif 000006
【図５】
0006598667.tif 000007
【図６】
0006598667.tif 000008
【図７】
0006598667.tif 000009