(57)【要約】 【課題】  本発明は、音響波が固体を透過する際の波形の歪みに起因する、音響波の受信信号の歪みを補正することができる処理システムを提供することを目的とする。【解決手段】  本発明に係る処理システムは、音響波が第１の媒質を透過するときの、複数の周波数に対応する複素透過率を表す透過率フィルタを取得する透過率フィルタ取得手段と、透過率フィルタを用いて電気信号を補正することにより、補正された電気信号を取得する補正手段と、を有する。【選択図】  図１
2017184972.tif 000002
【特許請求の範囲】
【請求項１】
被検体で発生し、第１の媒質を透過した音響波を受信手段が受信することにより得られた電気信号を処理するための処理システムであって、
音響波が前記第１の媒質を透過するときの、複数の周波数に対応する複素透過率を表す透過率フィルタを取得する透過率フィルタ取得手段と、
前記透過率フィルタを用いて前記電気信号を補正することにより、補正された電気信号を取得する補正手段と、を有する
ことを特徴とする処理システム。
【請求項２】
前記補正手段は、
前記電気信号を前記透過率フィルタでデコンボリューションすることにより、前記補正された電気信号を取得する
ことを特徴とする請求項１に記載の処理システム。
【請求項３】
前記被検体内の関心位置で発生し、前記受信手段に到達する音響波の前記第１の媒質への入射角度情報を取得する角度取得手段を有し、
前記透過率フィルタ取得手段は、前記入射角度情報に基づいて前記透過率フィルタを取得する
ことを特徴とする請求項２に記載の処理システム。
【請求項４】
前記透過率フィルタ取得手段は、前記入射角度情報に基づいて前記透過率フィルタを計算する
ことを特徴とする請求項３に記載の処理システム。
【請求項５】
入射角度情報と透過率フィルタとの関係を表す関係テーブルまたは関係式が保存された記憶手段を更に有し、
前記透過率フィルタ取得手段は、前記入射角度情報に基づいて前記関係テーブルまたは前記関係式を参照することにより前記入射角度に対応する前記透過率フィルタを取得する
ことを特徴とする請求項４に記載の処理システム。
【請求項６】
前記被検体内での音速と、前記第１の媒質と前記被検体との間に配置された第２の媒質内での音速とを取得する音速取得手段を有し、
前記角度取得手段は、前記被検体内での音速と前記第２の媒質内での音速とを用いて、前記入射角度情報を取得する
ことを特徴とする請求項３から５のいずれか１項に記載の処理システム。
【請求項７】
前記第２の媒質内での音速が保存された記憶手段を更に有し、
前記音速取得手段は、
前記第２の媒質内での音速を前記記憶手段から読み出すことにより取得し、
前記電気信号と前記第２の媒質内での音速とを用いて、前記被検体内での音速を取得する
ことを特徴とする請求項６に記載の処理システム。
【請求項８】
前記被検体内の関心位置と前記受信手段との位置関係及び前記透過率フィルタの関係を表す関係テーブルまたは関係式が保存された記憶手段と、
前記関心位置と前記受信手段との位置関係に関する情報を取得する位置取得手段と、を有し、
前記透過率フィルタ取得手段は、前記位置関係に関する情報に基づいて前記関係テーブルまたは前記関係式を参照することにより、前記位置関係に対応する前記透過率フィルタを取得する
ことを特徴とする請求項１または２に記載の処理システム。
【請求項９】
前記受信手段の位置情報を測定する位置測定手段を更に有し、
前記位置取得手段は、前記測定手段により測定された前記受信手段の位置情報を用いて、前記位置関係を取得する
ことを特徴とする請求項８に記載の処理システム。
【請求項１０】
前記位置取得手段は、
前記補正された電気信号を用いて前記関心位置の被検体情報を取得する
ことを特徴とする請求項３から９のいずれか１項に記載の処理システム。
【請求項１１】
前記受信手段の位置情報を測定する位置測定手段と、
ずれ量推定手段と、を更に有し、
前記受信手段は、少なくとも１つの受信素子を含み、
前記少なくとも１つの受信素子は、所定の音源を含む前記被検体で発生した音響波を受信することにより前記電気信号を出力し、
前記位置取得手段は、
前記位置測定手段により測定された前記受信手段の位置情報を用いて、前記少なくとも１つの受信素子が所定の位置に配置されていると仮定したときの、前記少なくとも１つの受信素子と前記所定の音源との距離を取得し、
前記補正された電気信号に基づいて、前記所定の音源と前記少なくとも１つの受信素子との距離を推定し、
前記ずれ量推定手段は、前記少なくとも１つの受信素子が前記所定の位置に配置されていると仮定したときの前記距離と、前記補正された電気信号に基づいて推定された前記距離とに基づいて、前記少なくとも１つの受信素子の、前記所定の位置からのずれ量を推定する
ことを特徴とする請求項１から１０のいずれか１項に記載の処理システム。
【請求項１２】
前記補正手段は、
前記電気信号のＳＮ比を評価し、
前記ＳＮ比が閾値よりも小さいときに前記透過率フィルタに基づいて前記電気信号に対する位相補正を行い、振幅補正を行わないデコンボリューションを行うことにより前記補正された電気信号を取得し、
前記ＳＮ比が閾値よりも大きいときに前記透過率フィルタに基づいて前記電気信号に対する位相補正も振幅補正も行うデコンボリューションを行うことにより前記補正された電気信号を取得する
ことを特徴とする請求項１から１１のいずれか１項に記載の処理システム。
【請求項１３】
前記透過率フィルタ取得手段は、
前記第１の媒質の厚み情報に基づいて前記透過率フィルタを変更する
ことを特徴とする請求項１から１２のいずれか１項に記載の処理システム。
【請求項１４】
請求項１から１３のいずれか１項に記載の処理システムと、
光を前記被検体に照射する光照射手段と、
前記光照射手段からの光が前記被検体に照射されることにより発生する音響波を受信することにより前記電気信号を出力する前記受信手段と、を有する
ことを特徴とする情報取得装置。
【請求項１５】
請求項１から１３のいずれか１項に記載の処理システムと、
前記第１の媒質を介して超音波を前記被検体に送受信する送受信手段を更に有し、
前記送受信手段は、前記被検体内で反射した前記超音波を前記音響波として受信することにより前記電気信号を出力し、
前記処理手段は、
超音波の送信時および受信時における、超音波が前記第１の媒質を透過するときの前記透過率フィルタを取得する
ことを特徴とする情報取得装置。
【請求項１６】
前記受信手段の受信帯域の中心周波数は、
当該中心周波数を有する音響波の水中での波長が前記所定の音源のサイズよりも大きくなるような中心周波数である
ことを特徴とする請求項１４または１５に記載の情報取得装置。
【請求項１７】
前記受信手段の受信帯域の中心周波数は、
当該中心周波数を有する音響波の水中での波長が前記所定の音源のサイズの１／５となるような中心周波数である
ことを特徴とする請求項１４または１５に記載の情報取得装置。
【請求項１８】
前記第１の媒質は、前記被検体を保持する保持手段である
ことを特徴とする請求項１４から１７のいずれか１項に記載の情報取得装置。
【請求項１９】
被検体で発生し、第１の媒質を透過した音響波に由来する電気信号を処理する信号処理方法であって、
音響波が前記第１の媒質を透過するときの、複数の周波数に対応する複素透過率を表す透過率フィルタを取得する工程と、
前記透過率フィルタを用いて前記電気信号を補正することにより、補正された電気信号を取得する工程と、
を有する
ことを特徴とする信号処理方法。
【請求項２０】
請求項１９に記載の信号処理方法をコンピュータに実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、音響波に由来する電気信号を処理する処理システムに関する。
【背景技術】
【０００２】
音響波を受信することにより得られた受信信号を用いて画像化する手法として、光音響トモグラフィ（ＰＡＴ：Ｐｈｏｔｏａｃｏｕｓｔｉｃ  Ｔｏｍｏｇｒａｐｈｙ）や超音波エコートモグラフィなどが提案されている。
特許文献１では、被検体と音響波の探触子との間に配置された圧迫板に起因する音線の屈折を考慮して超音波画像を取得する超音波エコー装置が開示されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】米国特許第６６０７４８９号明細書
【発明の概要】
【発明が解決しようとする課題】
【０００４】
一般に、乳房などの生体中あるいは液体中では、音響波は縦波として伝搬する。一方、固体中では、音響波は縦波だけでなく横波でも伝搬可能である。生体中あるいは液体中を伝搬した音響波（縦波）は、固体中を、縦波および横波が共存した状態で伝搬する。さらに、固体中を伝搬し、探触子がある側の液体の媒質へ到達すると、音響波は再び縦波となって媒質中を伝搬する。縦波と横波とでは固体中を伝搬する際の音速が異なっている。
【０００５】
そのため、液体、固体、液体を縦波−縦波−縦波で伝搬した音響波（音響波１）と、縦波−横波−縦波で伝搬した音響波（音響波２）では位相や強度が異なる。媒質中の音響波は音響波１と音響波２とが干渉した状態で伝搬するが、干渉の程度は固体に入射する音響波の角度に依存する。また、固体の厚さが音響波の波長と同程度、あるいはそれ以下の場合は、スネルの法則で規定される臨界角以上で音響波が入射した場合でもその一部は固体（保持手段）を透過することができる。臨界角は、スネルの法則により出射角が９０度となる入射角のことであり、全反射角とも呼ぶ。
【０００６】
これらの現象などに起因して、音響波は、固体に入射する音響波の角度に依存して、波形が歪んだ状態で探触子に到達する。
【０００７】
特許文献１の方法では、圧迫板に入射する音響波の角度に依存した音響波の波形の変化を考慮していない。このため、特許文献１の方法を用いた場合、音響波の波形の変化に伴う受信信号の歪みによって、得られる情報の精度が低下する。
【０００８】
そこで、本発明は、音響波が固体を透過する際の波形の歪みに起因する、音響波の受信信号の歪みを補正することができる処理システムを提供することを目的とする。
【課題を解決するための手段】
【０００９】
本発明に係る処理システムは、音響波が第１の媒質を透過するときの、複数の周波数に対応する複素透過率を表す透過率フィルタを取得する透過率フィルタ取得手段と、透過率フィルタを用いて電気信号を補正することにより、補正された電気信号を取得する補正手段と、を有する。
【発明の効果】
【００１０】
本発明に係る処理システムによれば、音響波が固体を透過する際の波形の歪みに起因する、音響波の受信信号の歪みを補正することができる。
【図面の簡単な説明】
【００１１】
【図１】実施例１に係る光音響装置の模式図
【図２】実施例１に係るコンピュータ周辺の構成の模式図
【図３】実施例１に係る被検体情報取得フローを示す図
【図４】実施例１に係る音線の計算方法を説明するための図
【図５】実施例１に係る受信信号に対するデコンボリューション処理を説明するための図
【図６】実施例１に係る伝搬時間の計算方法を説明するための図
【図７】実施例１に係る光音響画像情報を示す図
【図８】実施例２に係る光音響装置の模式図
【図９】実施例３に係る超音波エコー装置の模式図
【図１０】実施例３に係る超音波の伝搬を示す模式図
【図１１】実施例３に係る被検体情報取得フローを示す図
【図１２】実施例３に係る超音波の入射角度を示す模式図
【図１３】実施例４に係る光音響装置の模式図
【図１４】実施例４に係る被検体情報取得フローを示す図
【図１５】実施例４に係る補正された受信信号の例
【図１６】実施例４に係る別の被検体情報取得フローを示す図
【図１７】実施例４によって得られた位置補正量を示す図
【図１８】実施例５によって得られた保持カップの厚さを示す図
【図１９】実施例５における保持カップの幾何学情報を示す図
【発明を実施するための形態】
【００１２】
以下、図面を参照しつつ本発明をより詳細に説明する。なお、同一の構成要素には原則として同一の参照番号を付して、説明を省略する。
【００１３】
＜光音響装置＞
図１は、本実施形態にかかる光音響装置の模式図である。以下、装置の各構成要素について説明する。図１に示す装置は、光照射部１１０、保持カップ１２０、探触子１３０、信号データ収集部１４０、コンピュータ１５０、表示部１６０、入力部１７０を有する。測定対象は、被検体１００である。
【００１４】
光照射部１１０がパルス光１１３を被検体１００に照射し、被検体１００内で音響波が発生する。光に起因して光音響効果により発生する音響波を光音響波とも呼ぶ。探触子１３０は、光音響波を受信することによりアナログ信号としての電気信号を出力する。信号データ収集部１４０は、探触子１３０から出力されたアナログ信号としての電気信号をデジタル信号に変換し、コンピュータ１５０に出力する。コンピュータ１５０は、信号データ収集部１４０から出力されたデジタル信号を、光音響波に由来する信号データとして記憶する。
【００１５】
コンピュータ１５０は、記憶されたデジタル信号に対して信号処理を行うことにより、被検体１００に関する情報（被検体情報）を表す画像データを生成する。また、コンピュータ１５０は、得られた画像データに対して画像処理を施した後に、画像データを表示部１６０に出力する。表示部１６０は、被検体１００に関する情報の画像を表示する。ユーザーとしての医師は、表示部１６０に表示された被検体に関する情報の画像を確認することにより、診断を行うことができる。
【００１６】
本実施形態に係る光音響装置により得られる被検体情報は、光音響波の発生音圧（初期音圧）、光エネルギー吸収密度、光吸収係数、及び被検体を構成する物質の濃度に関する情報などの少なくとも１つである。物質の濃度に関する情報とは、オキシヘモグロビン濃度、デオキシヘモグロビン濃度、総ヘモグロビン濃度、または酸素飽和度等である。総ヘモグロビン濃度とは、オキシヘモグロビン濃度およびデオキシヘモグロビン濃度の和のことである。酸素飽和度とは、全ヘモグロビンに対するオキシヘモグロビンの割合のことである。本実施形態に係る光音響装置は、被検体内の各位置（２次元または３次元の空間の各位置）における上記情報の値を表す画像データを取得する。すなわち、本実施形態に係る光音響装置は、被検体情報を取得する被検体情報取得装置ともいえる。
【００１７】
以下、本実施形態に係る光音響装置の各構成の詳細を説明する。
【００１８】
（光照射部１１０）
光照射部１１０は、パルス光１１３を発する光源１１１と、光源１１１から射出されたパルス光１１３を被検体１００へ導く光学系１１２とを含む。
【００１９】
光源１１１が発する光のパルス幅としては、１ｎｓ以上、１００ｎｓ以下のパルス幅であってもよい。また、光の波長としては４００ｎｍから１６００ｎｍ程度の範囲の波長であってもよい。生体表面近傍の血管を高解像度でイメージングする場合は、血管での吸収が大きい波長（４００ｎｍ以上、７００ｎｍ以下）としてもよい。一方、生体の深部をイメージングする場合には、生体の背景組織（水や脂肪など）において典型的に吸収が少ない波長（７００ｎｍ以上、１１００ｎｍ以下）の光を用いてもよい。
【００２０】
光源１１１としては、レーザーや発光ダイオードを用いることができる。また、複数波長の光を用いて測定する際には、波長の変換が可能な光源であってもよい。なお、複数波長を被検体に照射する場合、互いに異なる波長の光を発生する複数台の光源を用意し、それぞれの光源から交互に照射することも可能である。複数台の光源を用いた場合もそれらをまとめて光源として表現する。レーザーとしては、固体レーザー、ガスレーザー、色素レーザー、半導体レーザーなど様々なレーザーを使用することができる。例えば、Ｎｄ：ＹＡＧレーザーやアレキサンドライトレーザーなどのパルスレーザーを光源１１１として用いてもよい。また、Ｎｄ：ＹＡＧレーザー光を励起光とするＴｉ：ｓａレーザーやＯＰＯ（Ｏｐｔｉｃａｌ  Ｐａｒａｍｅｔｒｉｃ  Ｏｓｃｉｌｌａｔｏｒｓ）レーザーを光源１１１として用いてもよい。また、光源１１１としてマイクロウェーブ源を用いてもよい。
【００２１】
光学系１１２には、レンズ、ミラー、光ファイバ等の光学素子を用いることができる。乳房等を被検体１００とする場合、パルス光のビーム径を広げて照射することが好ましいため、光学系１１２の光出射部は光を拡散させる拡散板等で構成されていてもよい。一方、光音響顕微鏡においては、解像度を上げるために、光学系１１２の光出射部はレンズ等で構成し、ビームをフォーカスして照射してもよい。また、光伝送手段としての光学系１１２に複数の光ファイバを束ねたファイバーバンドルを採用してもよい。また、光伝送手段としての光学系１１２に、ミラーを内包した関節によって複数の中空の導波管を接続した光学系（多関節アームとも呼ぶ）を採用してもよい。
【００２２】
なお、光照射部１１０が光学系１１２を備えずに、光源１１１から直接被検体１００にパルス光１１３を照射してもよい。
【００２３】
（保持カップ１２０）
第１の媒質である保持手段としての保持カップ１２０は被検体１００の形状を測定中に保持するために使用される。保持カップ１２０により被検体１００を保持することによって、被検体１００の動きの抑制および被検体１００の位置を保持カップ１２０内に留めることができる。保持カップ１２０の材料には、ポリカーボネートやポリエチレン、ポリエチレンテレフタレート等、樹脂材料を用いることができる。
【００２４】
保持カップ１２０は、被検体１００を保持できる硬度を有する材料であることが好ましい。保持カップ１２０は、測定に用いる光を透過する材料であってもよい。保持カップ１２０は、インピーダンスが被検体１００と同程度の材料で構成されていてもよい。乳房等の曲面を有するものを被検体１００とする場合、凹型に成型した保持カップ１２０であってもよい。この場合、保持カップ１２０の凹部分に被検体１００を挿入することができる。
【００２５】
保持カップ１２０は、取り付け部１２１に取り付けられている。取り付け部１２１は、被検体の大きさに合わせて複数種類の保持カップ１２０を交換可能に構成されていてもよい。例えば、取り付け部１２１は、曲率半径や曲率中心などの異なる保持カップに交換できるように構成されていてもよい。
【００２６】
また、保持カップ１２０には保持カップ１２０の緒元が登録されたタグ１２２が設置されていてもよい。例えば、タグ１２２には、保持カップ１２０の曲率半径、曲率中心、縦波音速、横波音速、識別ＩＤ等の緒元を登録することができる。タグ１２２に登録された緒元は、読み取り部１２３により読み出され、コンピュータ１５０に転送される。保持カップ１２０が取り付け部１２１に取り付けられたときに容易にタグ１２２を読み取るために、読み取り部１２３は取り付け部１２１に設置されていてもよい。例えば、タグ１２２はバーコードであり、読み取り部１２３はバーコードリーダである。
【００２７】
（探触子１３０）
受信手段としての探触子１３０は、音響波を受信することにより電気信号を出力するトランスデューサ１３１と、トランスデューサ１３１を支持する支持体１３２とを含む。
【００２８】
トランスデューサ１３１を構成する部材としては、ＰＺＴ（チタン酸ジルコン酸鉛）に代表される圧電セラミック材料や、ＰＶＤＦ（ポリフッ化ビニリデン）に代表される高分子圧電膜材料などを用いることができる。また、圧電素子以外の素子を用いてもよい。例えば、静電容量型トランスデューサ（ＣＭＵＴ：Ｃａｐａｃｉｔｉｖｅ  Ｍｉｃｒｏ−ｍａｃｈｉｎｅｄ  Ｕｌｔｒａｓｏｎｉｃ  Ｔｒａｎｓｄｕｃｅｒｓ）、ファブリペロー干渉計を用いたトランスデューサなどを用いることができる。なお、音響波を受信することにより電気信号を出力できる限り、いかなるトランスデューサを採用してもよい。また、トランスデューサにより得られる信号は時間分解信号である。つまり、受信素子により得られる信号の振幅は、各時刻にトランスデューサで受信される音圧に基づく値（例えば、音圧に比例した値）を表したものである。
【００２９】
光音響波を構成する周波数成分は、典型的には１００ＫＨｚから１００ＭＨｚであるので、トランスデューサ１３１として、これらの周波数を検出することのできるものを採用してもよい。
【００３０】
支持体１３２は、機械的強度が高い金属材料やプラスチック材料などから構成されていてもよい。本実施例において支持体１３２は半球殻形状であり、半球殻上に複数のトランスデューサ１３１を支持できるように構成されている。この場合、支持体１３２に配置されたトランスデューサ１３１の指向軸は半球の曲率中心付近に集まる。そして、複数のトランスデューサ１３１から出力された電気信号群を用いて画像化したときに曲率中心付近の画質が高くなる。なお、支持体１３２はトランスデューサ１３１を支持できる限り、いかなる構成であってもよい。支持体１３２は、１Ｄアレイ、１．５Ｄアレイ、１．７５Ｄアレイ、２Ｄアレイと呼ばれるような平面又は曲面内に、複数のトランスデューサを並べて配置してもよい。
【００３１】
また、支持体１３２は音響マッチング材１９０を貯留する容器としても機能してもよい。すなわち、支持体１３２をトランスデューサ１３１と被検体１００との間に音響マッチング材１９０を配置するための容器としてもよい。
【００３２】
また、探触子１３０が、トランスデューサ１３１から出力される時系列のアナログ信号を増幅する増幅器を備えてもよい。また、探触子１３０が、トランスデューサ１３１から出力される時系列のアナログ信号を時系列のデジタル信号に変換するＡ／Ｄ変換器を備えてもよい。すなわち、探触子１３０が後述する信号データ収集部１４０を備えてもよい。
【００３３】
なお、音響波を様々な角度で検出できるようにするために、理想的には被検体１００を全周囲から囲むようにトランスデューサ１３１を配置することが好ましい。ただし、被検体１００が大きく全周囲を囲むようにトランスデューサを配置できない場合は、図１に示したように半球状の支持体上にトランスデューサを配置して全周囲を囲む状態に近づけてもよい。なお、トランスデューサの配置や数及び支持体の形状は被検体に応じて最適化すればよく、本発明に関してはどのような探触子１３０を採用してもよい。
【００３４】
（信号データ収集部１４０）
信号データ収集部１４０は、トランスデューサ１３１から出力されたアナログ信号である電気信号を増幅するアンプと、アンプから出力されたアナログ信号をデジタル信号に変換するＡ／Ｄ変換器とを含む。信号データ収集部１４０は、ＦＰＧＡ（Ｆｉｅｌｄ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ）チップなどで構成されてもよい。信号データ収集部１４０から出力されるデジタル信号は、コンピュータ１５０内の記憶手段に記憶される。信号データ収集部１４０は、Ｄａｔａ  Ａｃｑｕｉｓｉｔｉｏｎ  Ｓｙｓｔｅｍ（ＤＡＳ）とも呼ばれる。本明細書において電気信号は、アナログ信号もデジタル信号も含む概念である。なお、信号データ収集部１４０は、光照射部１１０の光射出部に取り付けられた光検出センサと接続されており、パルス光１１３が光照射部１１０から射出されたことをトリガーに、同期して処理を開始してもよい。
【００３５】
（コンピュータ１５０）
コンピュータ１５０は、処理手段、記憶手段、制御手段を含む。各構成の機能については処理フローの説明の際に説明する。
【００３６】
記憶手段は、ＲＯＭ（Ｒｅａｄ  ｏｎｌｙ  ｍｅｍｏｒｙ）、磁気ディスクやフラッシュメモリなどの非一時記憶媒体で構成することができる。また、記憶手段は、ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）などの揮発性の媒体であってもよい。なお、プログラムが格納される記憶媒体は、非一時記憶媒体である。
【００３７】
処理手段としての演算機能を担うユニットは、ＣＰＵ、ＧＰＵ（Ｇｒａｐｈｉｃｓ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）、ＤＳＰ（Ｄｉｇｉｔａｌ  Ｓｉｇｎａｌ  Ｐｒｏｃｅｓｓｏｒ）等のプロセッサ、ＦＰＧＡ（Ｆｉｅｌｄ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ）チップ等の演算回路で構成されることができる。これらのユニットは、単一のプロセッサや演算回路から構成されるだけでなく、複数のプロセッサや演算回路から構成されていてもよい。
【００３８】
制御手段は、ＣＰＵなどの演算素子で構成される。制御手段は、光音響装置の各構成の動作を制御する。制御手段は、入力部１７０からの測定開始などの各種操作による指示信号を受けて、光音響装置の各構成を制御してもよい。また、制御手段は、記憶手段に格納されたプログラムコードを読み出し、光音響装置の各構成の作動を制御する。
【００３９】
コンピュータ１５０は専用に設計されたワークステーションであってもよい。また、コンピュータ１５０の各構成は異なるハードウェアによって構成されてもよい。また、コンピュータ１５０の少なくとも一部の構成は単一のハードウェアで構成されてもよい。
【００４０】
図２は、本実施例に係るコンピュータ１５０の具体的な構成を示す。本実施例に係るコンピュータ１５０は、ＣＰＵ１５１、ＧＰＵ１５２、ＲＡＭ１５３、ＲＯＭ１５４、外部記憶装置１５５から構成される。また、コンピュータ１５０には、表示部１６０としての液晶ディスプレイ１６１、入力部１７０としてのマウス１７１、キーボード１７２が接続されている。
【００４１】
また、コンピュータ１５０および複数のトランスデューサ１３１は、共通の筺体に収められた構成で提供されてもよい。ただし、筺体に収められたコンピュータで一部の信号処理を行い、残りの信号処理を筺体の外部に設けられたコンピュータで行ってもよい。この場合、筺体の内部および外部に設けられたコンピュータを総称して、本実施形態に係るコンピュータとすることができる。
【００４２】
トランスデューサから出力された信号を処理するシステムをまとめて処理システムと称する。処理システムは複数の処理装置から構成されていてもよい。
【００４３】
（表示部１６０）
表示部１６０は、液晶ディスプレイや有機ＥＬ（Ｅｌｅｃｔｒｏ  Ｌｕｍｉｎｅｓｃｅｎｃｅ）などのディスプレイである。コンピュータ１５０により得られた被検体情報等に基づく画像や特定位置の数値等を表示する装置である。表示部１６０は、画像や装置を操作するためのＧＵＩを表示してもよい。なお、被検体情報の表示にあたっては、表示部１６０またはコンピュータ１５０において画像処理（輝度値の調整等）を行った上で表示することもできる。
【００４４】
（入力部１７０）
入力部１７０は、ユーザーが操作可能な、マウスやキーボードなどで構成されることができる。また、表示部１６０をタッチパネルで構成し、表示部１６０を入力部１７０としてもよい。
【００４５】
なお、光音響装置の各構成はそれぞれ別の装置として構成されてもよいし、一体となった１つの装置として構成されてもよい。また、光音響装置の少なくとも一部の構成が一体となった１つの装置として構成されてもよい。
【００４６】
（音響マッチング材１９０）
第２の媒質である音響マッチング材１９０は光音響装置の構成ではないが説明する。音響マッチング材１９０は、被検体１００とトランスデューサ１３１との間を音響波伝搬させるためのものである。音響マッチング材１９０は、変形可能な部材から構成され、被検体１００と接触したときに変形してもよい。すなわち、音響マッチング材１９０は、被検体１００とトランスデューサ１３１との間の隙間をできるだけ低減するために、被検体に合わせて変形することができる材料から構成されていてもよい。音響マッチング材１９０は音響波減衰が少ないものであってもよい。音響マッチング材１９０には、被検体１００の音響インピーダンスとトランスデューサ１３１の音響インピーダンスとの間の音響インピーダンスを有する部材を採用してもよい。特に、被検体の音響インピーダンスに近い材料を選択してもよい。照射光が音響マッチング材１９０を透過する場合、照射光に対して透明であってもよい。音響マッチング材１９０には、水、超音波ジェルなどを用いてもよい。
【００４７】
本実施形態においては、支持体１３２が音響マッチング材１９０を貯留する容器として機能する。なお、光音響装置は、支持体１３２とは別に、トランスデューサ１３１と被検体１００との間に音響マッチング材１９０を貯留することができる容器を備えていてもよい。容器としては、プラスチック容器や金属容器などを採用してもよい。
【００４８】
（被検体１００）
被検体１００は光音響装置を構成するものではないが、以下に説明する。本実施例に係る光音響装置は、人や動物の悪性腫瘍や血管疾患などの診断や化学治療の経過観察などを目的としてもよい。よって、被検体１００としては生体、具体的には人体や動物の乳房や頸部、腹部などの診断の対象部位が想定される。例えば、人体が測定対象であれば、オキシヘモグロビンあるいはデオキシヘモグロビンやそれらを含む多く含む血管あるいは腫瘍の近傍に形成される新生血管などを光吸収体の対象としてもよい。また、頸動脈壁のプラークなどを光吸収体の対象としてもよい。また、メチレンブルー（ＭＢ）、インドシニアングリーン（ＩＣＧ）などの色素、金微粒子、またはそれらを集積あるいは化学的に修飾した外部から導入した物質を光吸収体としてもよい。
【実施例１】
【００４９】
図１に示す光音響装置を用いて本実施例を説明する。光源１１１は、波長８００ｎｍ、パルス幅２０ｎｓｅｃ、繰り返し周波数１０Ｈｚ、パルスエネルギー１５０ｍＪのパルス光を発生するチタンサファイアレーザである。光学系１１２は、光伝送手段としてのファイバーバンドルと、光照射手段として光ビームを拡大するための凹レンズを内包している。ファイバーバンドルの終端部は凸レンズに接続されている。
【００５０】
支持体１３２は、光照射手段としての凸レンズとトランスデューサ１３１とを支持している。支持体１３２において、トランスデューサ１３１を支持している部分の内面が半径１２７ｍｍの半球状となっており、半球面に沿って５１２個のトランスデューサ１３１が配置されている。支持体１３２の内側には、音響マッチング材１９０としての水が貯留されている。
【００５１】
保持カップ１２０は、０．５ｍｍのポリエチレンテレフタレート樹脂で構成されている。被検体１００としては例えば人の乳房などが用いられる。保持カップ１２０において、被検体１００と接する部分は半径２３０ｍｍの球面となっている。被検体１００と保持カップ１２０との間は音響マッチングのため音響マッチング材としての超音波ジェル（図中不指示）で埋められている。
【００５２】
本実施例では、素子サイズ３ｍｍ角、中心検出周波数２ＭＨｚの圧電素子からなるトランスデューサ１３１を５１２個、半球面上に並べている。
【００５３】
以下、図３に示すフロー図に沿って、本実施例に係る被検体情報の取得方法を説明する。
【００５４】
（Ｓ１１０：光を照射する工程）
光源１１１で発生させた光は、光学系１１２を介してパルス光１１３として被検体１００に照射される。そして、被検体１００内部でパルス光１１３が吸収され、光音響効果により光音響波が生じる。
【００５５】
（Ｓ１２０：光音響波を受信する工程）
本工程では、探触子１３０が光音響波を受信して、トランスデューサ１３１から電気信号を出力する。出力された受信信号は、コンピュータ１５０に渡される。本実施例において信号データ収集部１４０は、トランスデューサ１３１から出力されたアナログ信号をサンプリングレート２０ＭＨｚでＡＤ変換し、信号データを取得する。また、信号データ収集部１４０は、冗長な情報を除くため、被検体１００にパルス光を照射した瞬間の信号をサンプル０番として、サンプル０〜９９９番の信号を除いて保存する。すなわち、信号データ収集部１４０は、サンプル１０００番から３０４７番まで２０４８サンプルの信号データを取得する。
【００５６】
なお、光音響装置は、光音響波を受信したときの探触子１３０の位置を測定する位置測定部を備えていてもよい。例えば、位置測定部としては、探触子１３０に搭載された磁気センサや、探触子１３０を機械的に走査するステージのエンコーダなどを用いることができる。例えば、光照射部１１０から光が射出されたことをトリガーとして、位置測定部がそのときの探触子１３０の位置情報をコンピュータ１５０に出力してもよい。コンピュータ１５０は、このようにして得られた探触子１３０の位置情報を後述するＳ１４０−１７０での処理に用いてもよい。
【００５７】
（Ｓ１３０：音速情報を取得する工程）
音速取得手段としてのコンピュータ１５０は、被検体１００内での縦波音速ｃ１、保持カップ１２０内（第１の媒質内）での縦波音速ｃ２ｌ、保持カップ１２０内での横波音速ｃ２ｔ、音響マッチング材１９０内（第２の媒質内）での縦波音速ｃ３を取得する。
【００５８】
保持カップ１２０内での縦波音速及び横波音速については予め測定したデータを記憶手段に保存しておき、本工程においてコンピュータ１５０が記憶手段から読み出すことにより取得してもよい。また、保持カップ１２０の温度に対する保持カップ１２０内での縦波音速及び横波音速の関係式または関係テーブルが記憶手段に予め保存されていてもよい。そして、本工程において温度測定部が保持カップ１２０の温度を測定し、コンピュータ１５０が測定された温度に対応する音速を関係式または関係テーブルに従って取得してもよい。
【００５９】
また、ユーザーが入力部１７０により保持カップ１２０内での音速を入力し、コンピュータ１５０がその情報を受け取ることにより取得してもよい。
【００６０】
また、コンピュータ１５０は、取り付け部１２１に取り付けられた保持カップ１２０に対応する保持カップ１２０内での音速を取得してもよい。例えば、読み取り部１２３が保持カップ１２０に付されたタグ１２２に登録された保持カップ１２０内での音速の情報を読み取り、コンピュータ１５０に転送してもよい。コンピュータ１５０は、読み取り部１２３により読み出された保持カップ１２０内での音速の情報を取得してもよい。また、ユーザーが入力部１７０により保持カップ１２０に付された識別ＩＤを入力し、コンピュータ１５０が入力された識別ＩＤに対応する保持カップ１２０内の音速を記憶手段から読み出すことにより取得してもよい。
【００６１】
被検体１００内での縦波音速についても保持カップ１２０内での縦波音速及び横波音速と同様の方法で取得してもよい。ただし、被検体１００内での縦波音速は被検体１００毎に異なるため、被検体１００毎に新たなデータを取得することが好ましい。コンピュータ１５０は、被検体１００から発生した音響波に由来する信号を用いて被検体１００内での音速を取得してもよい。例えば、コンピュータ１５０は、仮の音速を設定し、仮の音速を用いて信号から再構成画像を生成し、再構成画像の画質（コントラストや解像度など）を評価する。コンピュータ１５０はこの処理を異なる仮の音速に繰り返し適用し、再構成画像の画質が所定の閾値より高くなったときの仮の音速を被検体１００の音速として設定する。この方法以外にも、コンピュータ１５０は、仮の音速を用いて特定の位置で発生した音響波に由来する、複数のトランスデューサから出力された信号のばらつきを評価することにより被検体１００の音速を取得してもよい。この場合、コンピュータ１５０は、信号のばらつきが所定の閾値よりも小さくなるときの仮の音速を被検体１００の音速として取得する。なお、被検体１００から発生した音響波に由来する信号を用いて被検体１００の音速を取得できる上記以外の方法を採用してもよい。また、被検体１００以外の構成の音速については記憶手段に保存された既知の値を用いてもよい。また、保持カップ１２０が十分に薄い場合（音響波の波長よりも小さい場合）、保持カップ１２０を無視し、音響マッチング材１９０の音速を既知として被検体１００の音速を取得してもよい。この方法によれば、装置規模を大きくすることなく、被検体１００毎の固有の音速を取得することができる。コンピュータ１５０は、その他公知の手法により被検体１００内での音速を取得してもよい。
【００６２】
音響マッチング材１９０内での縦波音速についても保持カップ１２０内での縦波音速及び横波音速と同様の方法で取得してもよい。コンピュータ１５０は、被検体１００から発生した音響波に由来する信号を用いて音響マッチング材１９０内での音速を取得してもよい。また、温度測定部が音響マッチング材１９０の温度を測定し、コンピュータ１５０が測定された温度に対応する音速を関係式または関係テーブルに従って取得してもよい。コンピュータ１５０は、その他公知の手法により音響マッチング材１９０内での音速を取得してもよい。
【００６３】
上記ではコンピュータ１５０が音速そのものを取得する例を説明したが、本工程では音速を推定することのできるパラメータであればいかなるパラメータを取得してもよい。例えば、密度ρと体積弾性率Ｋとから音速は求められるため、本工程において密度ρと体積弾性率Ｋとを取得し、それらのパラメータから音速を推定してもよい。本明細書では、縦波または横波の伝搬速度（音速）の他に、これらの音速を推定することのできるパラメータも音速情報に含む。
【００６４】
（Ｓ１４０：入射角度情報を取得する工程）
角度取得手段としてのコンピュータ１５０は、図４（ａ）に示すように保持カップ１２０、トランスデューサ１３１を配置したときに、トランスデューサ１３１に到達する音響波が保持カップ１２０に入射したときの入射角度情報を取得する。本実施例では、コンピュータ１５０が、後に画像化されるボクセル１０１で発生した音響波が保持カップ１２０に入射するときの入射角度情報を取得する例を説明する。
【００６５】
本工程において、コンピュータ１５０は、保持カップ１２０が存在しないものとして音響波の音線を算出する。そして、コンピュータ１５０は、算出された音線が保持カップ１２０に入射する角度を入射角度情報として取得する。
【００６６】
図４（ｂ）に示すようにボクセル１０１で発生した音響波の音線を直線近似した場合の例を説明する。被検体１００と音響マッチング材１９０との音速差が小さい場合や、ボクセル１０１のサイズが比較的大きい場合は、図４（ｂ）に示すように音線を直線で近似してもよい。コンピュータ１５０は、屈折を考慮した場合と直線近似した場合の到達時間の差に相当する距離（到達時間の差×音速）が、ボクセルサイズよりも十分小さいときに直線近似を採用してもよい。例えば、その値がボクセルサイズの１／１０よりも小さいときに、到達時間の差に相当する距離が十分小さいといえる。
【００６７】
例えば、被検体１００の音速が１４６０ｍ／ｓ、音響マッチング材１９０の音速が１４８０ｍ／ｓ、ボクセルサイズが０．２５ｍｍであるとする。この場合、入射角５０度において、屈折を考慮した場合と直線近似した場合の到達時間の差は２．８ｎｓ、距離に換算して４．２μｍ程度である。この値はボクセルサイズの１／１０（２５μｍ）よりも小さいので、この例では直線近似してもよい。なお、到達時間の差は、音速差の二乗に比例し、ｔａｎ（入射角）の二乗に比例するので、これを考慮して、屈折を考慮するか直進近似とするかを選択すればよい。
【００６８】
コンピュータ１５０は、ボクセル１０１とトランスデューサ１３１とを直線で結び、この直線と保持カップ１２０とのなす角度θ１を入射角度情報として取得する。このとき、位置取得部としてのコンピュータ１５０は、トランスデューサ１３１の位置および保持カップ１２０の形状情報はコンピュータ１５０の記憶手段に予め保存されたものを読み出して取得してもよい。また、コンピュータ１５０は、保持カップ１２０の形状情報を取得する形状取得手段を用いて得られた形状情報を使用してもよい。
【００６９】
一方、被検体１００と音響マッチング材１９０との音速差が大きい場合や、ボクセル１０１のサイズが比較的小さい場合は、図４（ｃ）に示すように屈折を考慮した音線を算出する。コンピュータ１５０は、屈折を考慮した場合と直線近似した場合の到達時間の差に相当する距離（到達時間の差×音速）が、ボクセルサイズよりも十分大きいときに屈折を考慮した計算を採用してもよい。例えば、その値がボクセルサイズの１／１０よりも大きいときに、到達時間の差に相当する距離が十分大きいといえる。
【００７０】
例えば、被検体１００の音速が１４００ｍ／ｓ、音響マッチング材１９０の音速が１４８０ｍ／ｓ、ボクセルサイズが０．２５ｍｍであるとする。この場合、入射角５０度において、屈折を考慮した場合と直線近似した場合の到達時間の差は４８ｎｓ、距離に換算して７０μｍ程度である。この値はボクセルサイズの１／１０（２５μｍ）よりも大きいので、この例では屈折を考慮してもよい。
【００７１】
コンピュータ１５０は、式（１）で表すスネルの法則にしたがってボクセル１０１とトランスデューサ１３１とを結ぶ音線を算出し、音線と保持カップ１２０とのなす角度θ１を入射角度情報として取得する。
【００７２】
【数１】
2017184972.tif 000003
ここで、θ３は、音響マッチング材１９０中の音線と保持カップ１２０とのなす角度であり、以下、伝搬角と表す場合もある。
【００７３】
コンピュータ１５０は、本工程の処理を処理対象とする全ボクセルおよび全トランスデューサに対して実行する。
【００７４】
なお、位置取得部としてのコンピュータ１５０は、あるボクセルとあるトランスデューサとの位置関係を取得し、この位置関係により決定された入射角度情報を、近接するボクセルまたはトランスデューサに対応づけてもよい。また、コンピュータ１５０は、あるボクセルとあるトランスデューサとの位置関係により決定された入射角度情報に対して補間処理を行うことにより、その他のボクセルとトランスデューサとの組み合わせに対応する入射角度情報を取得してもよい。これらの場合、コンピュータ１５０は、全てのボクセルとトランスデューサとの組み合わせについて音線を算出しなくてもよい。
【００７５】
（Ｓ１５０：入射角度情報に基づいた透過率フィルタを取得する工程）
透過率フィルタ取得手段としてのコンピュータ１５０は、Ｓ１４０で取得した入射角度情報に基づいて、受信信号に対する透過率フィルタを取得する。本実施例において、透過率フィルタは、音響波が保持カップ１２０の被検体側からトランスデューサ側に透過する際の透過率を表す。透過率は、音響波の周波数に依存するため、周波数ごとに取得される。また、透過率は位相情報を持つため、複素数で表現される。以下、複数の周波数に対応する複素透過率をまとめて透過率フィルタと呼ぶ。なお、本発明で取り扱う音響波に含まれる周波数成分に応じて、透過率フィルタの周波数帯域を決定してもよい。
【００７６】
音響波が保持カップ１２０を透過する際、保持カップ１２０の第１面（被検体側の面）と第２面（探触子側の面）による多重反射に起因する透過波の歪みが発生する。本実施例のように保持カップ１２０が曲面の場合、この多重反射の影響が小さいことを本発明者は実験的に見出した。この場合、多重反射を考慮せずに透過率フィルタの計算を行ってもよい。
【００７７】
コンピュータ１５０は、保持カップの曲率半径に基づいて多重反射を考慮した透過率フィルタを適用するのか、多重反射を考慮しない透過率フィルタを適用するのかを選択してもよい。例えば、保持カップ１２０の曲率半径が３００ｍｍより小さい場合は、多重反射を考慮しない透過率フィルタを適用してもよい。一方、コンピュータ１５０は、保持カップ１２０の曲率半径が３００ｍｍより大きい場合は、多重反射を考慮した透過率フィルタを適用してもよい。また、コンピュータ１５０は、保持カップ１２０の曲率半径が１００ｍｍより小さい場合に、多重反射を考慮しない透過率フィルタを適用するようにしてもよい。また、コンピュータ１５０は、保持カップ１２０の曲率半径が１００ｍｍより大きい場合に、多重反射を考慮した透過率フィルタを適用するようにしてもよい。
【００７８】
また、コンピュータ１５０は、入射角度情報に基づいて多重反射を考慮した透過率フィルタを適用するのか、多重反射を考慮しない透過率フィルタを適用するのかを選択してもよい。例えば、入射角度が臨界角付近（例えば、臨界角の±５度）である場合、多重反射を考慮した透過率フィルタを適用し、それ以外の場合、多重反射を考慮しない透過率フィルタを適用してもよい。
【００７９】
また、コンピュータ１５０は、保持カップ１２０の厚さに基づいて多重反射を考慮した透過率フィルタを適用するのか、多重反射を考慮しない透過率フィルタを適用するのかを選択してもよい。
【００８０】
すなわち、コンピュータ１５０は、多重反射を考慮した透過率フィルタと、多重反射を考慮しない透過率フィルタという互いに異なる透過率フィルタを選択できるように構成されていてもよい。
【００８１】
本実施例では多重反射を考慮しない場合の透過率フィルタの計算について説明する。多重反射を考慮した場合の透過率フィルタの計算については実施例２で述べる。
【００８２】
被検体１００の音速（縦波）をｃ１、保持カップ１２０の縦波、横波に対する音速をそれぞれｃ２Ｌ、ｃ２Ｔ、音響マッチング材１９０の音速（縦波）をｃ３、入射角をθ１、保持カップ１２０の厚さをＴとする。また、被検体１００の音響インピーダンスをＺ１、保持カップ１２０の縦波、横波に対する音響インピーダンスをそれぞれＺ２Ｌ、Ｚ２Ｔ、音響マッチング材１９０の音響インピーダンスをＺ３とする。入射角θ１はＳ１４０計算されたものである。
【００８３】
保持カップ１２０中の縦波、横波に対する伝搬角θ２Ｌ、θ２Ｔ、および音響マッチング材１９０中の伝搬角θ３は、スネルの法則により以下の関係がある。
【００８４】
【数２】
2017184972.tif 000004
例えば、ｃ２Ｌがｃ１より大きい場合、ｓｉｎθ２Ｌ＞１となる場合がある。このとき、ｃｏｓθ２Ｌは虚数となる。ここで、保持カップ１２０の面に垂直な方向への伝搬に伴って減衰するように符号を決める。すなわち、
【００８５】
【数３】
2017184972.tif 000005
とする。ｃｏｓｑ２Ｔについても同様に、虚数となる場合は、
【００８６】
【数４】
2017184972.tif 000006
とする。
【００８７】
入射角θ１で保持カップ１２０の第１面に被検体１００側から入射した音響波（縦波）が縦波として保持カップ１２０側に透過する波の複素振幅透過率をｔ１２Ｌとする。また、このときに横波として透過する波の複素振幅透過率をｔ１２Ｔとする。この場合、ｔ１２Ｌ、ｔ１２Ｔは以下のように求められる。
【００８８】
【数５】
2017184972.tif 000007
保持カップ１２０中を第１面から第２面まで音響波が伝搬する際、伝搬に伴って音響波の位相がシフトする。位相シフトは音響波の周波数に依存する。また、ｃｏｓθ２Ｌあるいはｃｏｓθ２Ｔが虚数となる場合は減衰となる。音響波の角周波数をω、横波、縦波に対する位相シフト量（あるいは減衰量）をそれぞれφ２Ｌ、φ２Ｔとすると、φ２Ｌ、φ２Ｔは以下のように求められる。
【００８９】
【数６】
2017184972.tif 000008
保持カップ１２０の第２面に入射した縦波の音響波または横波の音響波が縦波として音響マッチング材１９０側に透過する波の複素振幅透過率をそれぞれｔ２３Ｌ、ｔ２３Ｔとすると、ｔ２３Ｌ、ｔ２３Ｔは以下のように求められる。
【００９０】
【数７】
2017184972.tif 000009
角周波数ωに対する保持カップ１２０の透過率フィルタｆ（ω）は次のように表すことができる。
【００９１】
【数８】
2017184972.tif 000010
ここで、後段のｅｘｐ項は、保持カップ１２０の領域に音響マッチング材１９０が存在した場合に発生する位相シフト量を補正するために設けている。このｅｘｐ項があることにより、Ｓ１６０で説明する受信信号補正を行ったときに、保持カップ１２０の領域を音響マッチング材１９０で置き換えた場合と等価な受信信号を得ることができる。
【００９２】
以上、コンピュータ１５０は、Ｓ１４０で取得した入射角度情報を用いて、透過率フィルタを取得できる。本実施例のおけるコンピュータ１５０は、さらに保持カップ１２０の厚さ、保持カップ１２０の縦波音速及び横波音速、被検体の縦波音速、音響マッチング材の縦波音速を用いて透過率フィルタを取得できる。
【００９３】
後述するＳ１６０でのデコンボリューションでは、受信信号を離散フーリエ変換して得た周波数空間での受信信号に対して処理を行う。本実施例の場合、受信信号を離散フーリエ変換すると、−９９９０２３４．３７５Ｈｚから９７６５．６２５Ｈｚごとに１０００００００Ｈｚまで２０４８点の周波数情報が得られる。したがって、本工程では、後述するＳ１６０での処理に必要な２０４８点の角周波数に対応して、式１１で示された透過率フィルタは計算される。
【００９４】
以上、コンピュータ１５０が、すべての音響波の伝搬経路における透過率フィルタを毎回計算する例を示したが、透過率フィルタの取得方法はこれに限ったものではない。例えば、コンピュータ１５０内の記憶手段に入射角度と透過率フィルタとの関係を表す関係式または関係テーブルを記憶しておいてもよい。そして、コンピュータ１５０が、Ｓ１４０で取得された入射角度情報をもとに記憶手段に保存された関係式または関係テーブルを参照して、その入射角度に対する透過率フィルタを求めてもよい。なお、関係テーブルに含まれない入射角度に対する透過率フィルタは、前後のデータから補間することで生成してもよい。このような構成をとれば、透過率フィルタを取得するための計算時間を短縮できる。
【００９５】
なお、すべてのトランスデューサ１３１に対して本工程を実施する。本工程の処理を、各トランスデューサに対して時系列に適用してもよいし、並列処理で複数のトランスデューサ１３１に対して同時に行ってもよい。
【００９６】
（Ｓ１６０：信号データに透過率フィルタをデコンボリューションする工程）
補正手段としてのコンピュータ１５０は、Ｓ１５０により得られた透過率フィルタを応答関数として、Ｓ１２０で得られた信号データをデコンボリューションすることで補正された信号データを取得する。この処理は、透過率フィルタから求めたデコンボリューションフィルタを信号データに重畳することで補正された信号データを得る処理、と言い換えることができる。デコンボリューションフィルタによる重畳処理は、保持カップ１２０を透過する際の音響波の変化を補正する、という効果がある。すなわち、このデコンボリューション処理により、トランスデューサが保持カップなしで音響波を受信した場合と等価な信号データを得ることができる。受信信号補正手段は、保持カップを透過することによって発生する信号データの歪みをなくすことができる。
【００９７】
さらに、コンピュータ１５０は、トランスデューサの受信特性（インパルス応答）をあらかじめ取得しておき、それを応答関数としてデコンボリューションする機能を設けてもよい。この処理により、インパルス応答による信号データの劣化を補正することができる。なお、インパルス応答を応答関数としたデコンボリューションと透過率フィルタを応答関数としたデコンボリューションとはどちらが先でもよく、また、同時であってもよい。
【００９８】
コンピュータ１５０の行うこの計算を以下で説明する。
【００９９】
あるトランスデューサ１３１における受信信号をＳ０（ｔ）、そのフーリエ変換信号をＦ［Ｓ０］（ω）、透過率フィルタをｆ（ω）とする。また、逆フーリエ変換をＦ−−１［］で表す。
【０１００】
ノイズがなく設計からのずれがない理想的な状態では、補正された受信信号Ｓ（ｔ）は以下のように表すことができる。
【０１０１】
【数９】
2017184972.tif 000011
ここで、ｒｅａｌは実部のみを抽出する関数である。
【０１０２】
実際には、ノイズがあったり、計算上ｆ（ω）が０となる入射角においてはゼロ割が発生するなど、式１２は現実的でない。したがって、一般的に知られるウィナーフィルターをデコンボリューションフィルタとして用いる。
【０１０３】
デコンボリューションフィルタをＤ（ω）として、Ｄ（ω）を以下のように求める。
【０１０４】
【数１０】
2017184972.tif 000012
ここで、ｃｏｎｊは複素共役、Ｃは定数である。Ｃについては信号を著しく劣化させないように経験的に設定するものである。一般的には｜ｆ（ω）｜の最大値の数％に設定するとよい。
【０１０５】
デコンボリューションフィルタＤ（ω）を用いて、以下のような計算で補正された受信信号Ｓ（ｔ）を求める。
【０１０６】
【数１１】
2017184972.tif 000013
このようにして得られた補正された受信信号Ｓ（ｔ）は保持カップ１２０の領域を音響マッチング材１９０で置き換えた場合とほぼ等価な受信信号となっている。
【０１０７】
実測された受信信号を用いて信号の補正を行った場合の例を図５に示す。ここでは、音源（光吸収体）として０．３ｍｍのゴム球をトランスデューサ１３１が配置された球面の中心近傍に置いている。音源の周りは水を満たしている。音速は、ｃ１＝１４８６ｍ／ｓ、ｃ３＝１４８６ｍ／ｓ、ｃ２Ｌ＝２３４０ｍ／ｓ、ｃ２Ｔ＝９１２ｍ／ｓを用いた。また、音響インピーダンスは、Ｚ１＝１４８５６００ｋｇ／ｍ２ｓ、Ｚ３＝１４８５６００ｋｇ／ｍ２ｓ、Ｚ２Ｌ＝２９７１８００ｋｇ／ｍ２ｓ、Ｚ２Ｔ＝１１５８２４０ｋｇ／ｍ２ｓを用いた。
【０１０８】
図５（ａ）は、保持カップ１２０への入射角が４７度となる位置のトランスデューサで受信された受信信号、図５（ｂ）は式１４を用いて補正された受信信号である。また、図５（ｃ）は比較例であり、保持カップ１２０をはずした状態で、同じ音源からの信号を図５（ａ）と同じトランスデューサで受信された受信信号である。図５中、横軸はサンプリングナンバー、縦軸は信号強度を示している。図５（ｃ）の波形に比べ、図５（ａ）では波形が歪んでいるが、図５（ｂ）の補正された受信信号では図５（ｃ）の波形に近くなっており、波形の補正が行われていることがわかる。
【０１０９】
また、図５（ａ）から図５（ｂ）の変換により信号の振幅も増加していることがわかる。これは、式１４の処理が保持カップ１２０の透過の際の損失（反射による損失）を補う効果があることを意味する。
【０１１０】
保持カップ１２０の透過の際の損失を補うことは、通常は画像の質を向上する方向に働く。ただ、信号強度が小さくＳＮ比が小さいといった場合にはノイズを持ち上げてしまい、画像の質を低下させる場合もある。その場合は、位相シフトによる形状の補正は行うが振幅の補正は行わないようにデコンボリューションフィルタを変えてもよい。具体的には、
【０１１１】
【数１２】
2017184972.tif 000014
が一定値となるようにＤ（ω）を規格化すればよい。
【０１１２】
その他、透過率フィルタｆ（ω）の複素共役をデコンボリューションフィルタＤ（ω）とする、という方法もある。これも、振幅の補正は行わずに位相シフトによる形状の補正は行う、という効果がある。
【０１１３】
コンピュータ１５０は、信号データのＳＮ比を評価し、信号データのＳＮ比が閾値より小さい場合に振幅の補正は行わずに、位相シフトによる形状補正（位相補正）を行うデコンボリューションを行ってもよい。また、コンピュータ１５０は、信号データのＳＮ比が閾値より大きい場合に位相補正も振幅補正も含めた、式１４に記載されたようなデコンボリューションを行ってもよい。
【０１１４】
保持カップ１２０の前後の音速を一定とみなしてよい場合、補正された受信信号Ｓ（ｔ）は、図６（ａ）のように、指定されたボクセル１０１から発生した音響波が境界のない空間を伝搬してトランスデューサ１３１で受信された信号とみなすことができる。また、保持カップ１２０の前後の音速差を考慮する場合、補正された受信信号Ｓ（ｔ）は、図６（ｂ）に示すように伝搬した音響波の受信信号とみなすことができる。すなわち、補正された受信信号Ｓ（ｔ）は、指定されたボクセル１０１から発生した音響波が、保持カップ１２０の第１面（被検体側の面）で１回の屈折を伴って伝搬してトランスデューサ１３１で受信された信号とみなすことができる。図６（ｂ）において保持カップ１２０を点線で示した。
【０１１５】
なお、保持カップ１２０の厚さが十分小さい場合、本実施例のように保持カップの一面を音速の境界として透過率フィルタを取得しても誤差が小さい。この場合、計算が簡易であるため、計算量が少ない。一方、保持カップ１２０の厚さが大きくなると、実際の伝搬経路との計算上の伝搬経路との差が大きくなる。そこで、コンピュータ１５０は、保持カップ１２０の厚さが閾値より大きい場合、図６（ｃ）に示すように、保持カップ１２０の内部を音速の境界として透過率フィルタを取得してもよい。特に保持カップ１２０の中心を音速の境界とすることで誤差を小さくすることができる。図６（ｃ）において、保持カップ１２０を点線で表し、音速の境界面を実線で表した。コンピュータ１５０は、保持カップの厚さに基づいてどの面を音速の境界とするのかを決定してもよい。例えば、コンピュータ１５０は、保持カップ１２０の厚さが１ｍｍより大きい場合に、保持カップ１２０の内部を音速の境界として計算するようにしてもよい。一方、コンピュータ１５０は、保持カップ１２０の厚さが１ｍｍより小さい場合に、保持カップ１２０の表面を音速の境界として計算するようにしてもよい。図６（ｃ）のように保持カップ１２０の内部（例えば中央）を音速の境界として計算する場合、（１１）の後段のexp項は、保持カップ１２０の領域の半分に被検体１００が存在するものとして計算すべきである。また、保持カップ１２０の領域のもう半分に音響マッチング材１９０が存在するものとして計算すべきである。
【０１１６】
なお、すべてのトランスデューサ１３１に対して本工程を実施する。本工程の処理を、各トランスデューサに対して時系列に適用してもよいし、並列処理で複数のトランスデューサ１３１に対して同時に行ってもよい。
【０１１７】
（Ｓ１７０：補正された信号データを用いて被検体情報を取得する工程）
画像再構成手段としてのコンピュータ１５０は、Ｓ１６０で得られた補正された受信信号Ｓ（ｔ）を再構成処理して指定されたボクセル１０１の被検体情報を得る。再構成アルゴリズムはユニバーサルバックプロジェクション（Ｕｎｉｖｅｒｓａｌ  Ｂａｃｋ−ｐｒｏｊｅｃｔｉｏｎ、以下ＵＢＰ）アルゴリズムを用いる。Ｓ１６０で説明した受信信号の補正により、保持カップ１２０を透過することに伴う波形の歪みは補正されている。そのため、本工程では、図６（ａ）、図６（ｂ）、または図６（ｃ）で表された系で受信されたものとして再構成処理を行えばよい。なお、補正された信号データから被検体情報を得られる限り、いかなる再構成処理を適用してもよい。
【０１１８】
図６（ａ）のように保持カップ１２０の前後の音速を一定とみなしてよい場合、指定されたボクセル１０１とトランスデューサ１３１の距離を、被検体１００または音響マッチング材１９０の縦波音速で割った時間が到達時間となる。補正された受信信号Ｓ（ｔ）中のその到達時間に相当する情報（振幅）が指定されたボクセル１０１の被検体情報のもととなる情報となる。
【０１１９】
また、保持カップ１２０の前後の音速差を考慮する場合の到達時間は以下のようになる。すなわち、図６（ｂ）において、指定されたボクセル１０１と屈折する境界までの距離を被検体１００の音速で割った時間と、屈折する境界からトランスデューサ１３１−ｉまでの距離を音響マッチング材１９０の音速で割った時間との和が到達時間となる。補正された受信信号Ｓ（ｔ）中のその到達時間に相当する情報（振幅）が指定されたボクセル１０１の被検体情報のもととなる情報となる。これは図６（ｃ）においても同様である。
【０１２０】
なお、コンピュータ１５０は、被検体情報を取得した全てのボクセルに対して、Ｓ１４０〜Ｓ１７０の処理を行うことにより、二次元あるいは三次元の被検体情報の空間分布を取得することができる。
【０１２１】
コンピュータ１５０は上記処理により、被検体情報として二次元あるいは三次元の初期音圧分布を取得してもよい。また、コンピュータ１５０は、被検体１００に照射された光の被検体１００内の光フルエンス分布を取得し、初期音圧分布と光フルエンス分布とを用いて光吸収係数分布を取得してもよい。また、コンピュータ１５０は、光吸収係数分布を用いて酸素飽和度分布等の濃度分布を取得してもよい。例えば、複数の波長の光に対応する光吸収係数分布を用いて濃度分布を取得することができる。コンピュータ１５０は、本工程で得られた初期音圧分布、光吸収係数分布、または濃度分布等の被検体情報を表示部１６０に出力する。
【０１２２】
（Ｓ１８０：被検体情報を取得する工程）
コンピュータ１５０は、Ｓ１７０で得られた被検体情報を用いて、画像化領域の被検体情報を表示部１６０に表示させる。表示部１６０は、被検体情報として、初期音圧分布、吸収係数分布、または濃度分布（酸素飽和度分布）等を表示することができる。表示部１６０に表示された被検体情報は、音響波が保持カップ１２０を透過することに伴う波形の歪みが抑制された情報であるため、医師等の作業者が診断等に用いる上で好適な情報となっている。
【０１２３】
ここで、実測された受信信号を用いて光音響画像情報（三次元の被検体情報）を得た場合の例を図７に示す。ここでは、音源（光吸収体）として０．３ｍｍのゴム球をトランスデューサ１３１が配置された球面の中心近傍に置いている。音源の周りは水を満たしている。
【０１２４】
図７（ａ）、（ｂ）は、補正された受信信号を用いた場合に得られた光音響画像情報である。ここで、図１において、光照射部１１０から被検体１００の方向（紙面上方向）を＋Ｚ方向とし、それに垂直な面をＸＹ面としている。図７（ａ）はＸＹ断面、図７（ｂ）はＹＺ断面を示している。図７（ｃ）、（ｄ）は、受信信号を補正しないでそのまま用いた場合に得られる光音響画像情報である。図７（ｃ）はＸＹ断面、図７（ｄ）はＹＺ断面を示している。
【０１２５】
図７からわかるように、補正された受信信号を用いた場合は、点像として再現できている。一方、受信信号をそのまま用いた場合は音響波が保持カップ１２０を透過することに伴う波形の歪みにより、点がｚ方向に２つに分離してしまっていることがわかる。
【０１２６】
以上説明したように、本実施例の光音響装置によれば、音響波が保持カップを透過する際の波形の歪み（振幅、位相）を補正することで、被検体情報を高精度に得ることが可能となる。
【０１２７】
なお、図には示していないが、支持体１３２が走査手段としてのステージに接続されており走査可能な構成となっていてもよい。この場合、より体積の大きい被検体情報を取得することができる。
【実施例２】
【０１２８】
図８は、実施例２に係る光音響装置を示す模式図である。なお、同一の構成要素には原則として同一の参照番号を付して、説明を省略する。実施例１では、乳房のような形状の被検体を測定対象として想定していたため保持カップ１２０を球面状としていた。本実施例では、手や足など比較的平たい形状の被検体を測定対象として想定し、保持面が平坦な保持カップ１２０を用いる。
【０１２９】
探触子１３０は、トランスデューサを２次元アレイ状に配置した受信手段である。支持体１３２は、光学系１１２、探触子１３０を支持している。また、支持体１３２は、音響マッチング材１９０である水を貯留する貯留手段としても機能する。
【０１３０】
保持手段としての保持板１２０は、０．５ｍｍのポリエチレンテレフタレート樹脂からなる。保持板１２０において、被検体１００と接する部分は平面となっている。被検体１００と保持板１２０との間は音響マッチングのため超音波ジェル（図中不指示）で埋められている。保持板１２０と探触子１３０との間隔は３０ｍｍになるように設定されている。
【０１３１】
本実施例では、探触子１３０として、素子サイズ２ｍｍ角、中心検出周波数１ＭＨｚの圧電素子からなるトランスデューサを２次元アレイ状に１０×２０個並べたものを用いている。
【０１３２】
本実施例では、信号データ収集部１４０が、被検体１００にパルス光を照射した瞬間の信号をサンプル０番として、サンプル０番から２０４７番までの２０４８サンプルの信号を取得する。
【０１３３】
本実施例のように保持板１２０の保持面が平面の場合、多重反射を考慮した透過率フィルタを用いたほうがよいことを本発明者は実験的に見出した。すなわち、本実施例においては、音響波が保持板１２０を透過する際、保持板１２０の第１面（被検体側の面）と第２面（探触子１３０側の面）による多重反射に起因した透過波の歪みが発生する。
【０１３４】
被検体１００の音速（縦波）をｃ１、保持板１２０の縦波、横波に対する音速をそれぞれｃ２Ｌ、ｃ２Ｔ、音響マッチング材１９０の音速（縦波）をｃ３、入射角をθ１、保持板１２０の厚さをＴとする。また、被検体１００の音響インピーダンスをＺ１、保持板１２０の縦波、横波に対する音響インピーダンスをそれぞれＺ２Ｌ、Ｚ２Ｔ、音響マッチング材１９０の音響インピーダンスをＺ３とする。入射角θ１はコンピュータ１５０によって計算されたものである。
【０１３５】
保持板１２０中の縦波、横波に対する伝搬角をθ２Ｌ、θ２Ｔ、および音響マッチング材１９０中の伝搬角をθ３とする。これらの関係は式１―４で説明したとおりである。
【０１３６】
入射角θ１で保持板１２０の第１面に被検体１００側から入射した音響波（縦波）が保持板１２０を透過して音響マッチング材１９０側に縦波として伝搬する波の透過率フィルタｆ（ω）は以下のように求められる。これは、保持板１２０中をある方向に伝搬する縦波、横波と逆方向に伝搬する縦波、横波を仮定し、保持板１２０の第１面と第２面での波の連続性を考慮して方程式を解くことで導くことができる。
【０１３７】
式１６―１９で表すように、Ｆ，Ｇ，Ｎ，Ｍを定義すると、多重反射を考慮した透過率フィルタｆ（ω）は式２０で表すことができる。
【０１３８】
【数１３】
2017184972.tif 000015
ここで、式１６−１９のＦ、Ｇ、Ｍ、Ｎは、式２０を簡略に表現するため便宜上導入した関数であり、角周波数ωの関数である。また、式２０の後段のｅｘｐ項は、保持板１２０の領域に音響マッチング材１９０が存在した場合に発生する位相シフト量を補正するために設けている。このｅｘｐ項があることにより、Ｓ１６０で説明するデコンボリューション処理を行ったときに、保持板１２０の領域を音響マッチング材１９０で置き換えた場合と等価な受信信号を得ることができる。
【０１３９】
本実施例においても、音響波が保持手段を透過する際の波形の歪み（位相、時間）を補正することで、高精細な光音響画像（被検体情報の空間分布）を得ることが可能となる。
【０１４０】
実施例１での透過率フィルタ（式１１に相当）は保持手段内（第１の媒質内）の多重反射がないものとして導いたフィルタである。一方、実施例２での透過率フィルタ（式２０に相当）は保持手段内の多重反射が理想的に発生するものとして導いたフィルタである。保持手段の形状によっては不完全な多重反射を考慮したほうが良い場合もある。その場合は、例えば、式１１と式２０の平均値を透過率フィルタとして用いたり、重みを付けた平均値を周波数フィルタとして用いたりしてもよい。重みについては、光音響画像の画質（例えばコントラスト）が最良になるように調整すればよい。なお、ユーザーが入力部１７０を用いて多重反射を考慮したフィルタと多重反射を考慮しないフィルタの重みを変更してもよい。
【実施例３】
【０１４１】
図９は、本実施例に係る被検体情報取得装置としての超音波エコー装置の模式図である。本実施例において探触子１３０は、被検体１００に超音波を送信し、被検体１００の内部で反射した超音波（エコー）を受信する。すなわち、本実施例において探触子１３０は、超音波を送受信する送受信手段として機能する。
【０１４２】
保持カップ１２０は、０．５ｍｍのポリエチレンテレフタレート樹脂からなる。保持カップ１２０において、被検体１００と接する部分は半径２３０ｍｍの球面となっている。被検体１００と保持カップ１２０との間は音響マッチングのため超音波ジェル（図中不指示）で埋められている。
【０１４３】
本実施例では、探触子１３０は、焦点距離３０ｍｍ、中心検出周波数６ＭＨｚの圧電素子のトランスデューサ１３１を備えている。
【０１４４】
送信制御手段としてのコンピュータ１５０は、トランスデューサ１３１にパルス信号を送って所望のパルス形状の超音波を送信させる。送信制御手段は一般にはパルサーとも呼ばれる。信号データ収集部１４０は、トランスデューサ１３１で受信された電気信号をサンプリングレート５０ＭＨｚでサンプリングする。信号データ収集部１４０は、被検体１００にパルス光を照射した瞬間の信号をサンプル０番として、サンプル１００番から２１４７番まで２０４８サンプルの信号を取得する。
【０１４５】
図１０（ａ）はトランスデューサ１３１から被検体１００側に超音波を送信している状態を示している。被検体１００中に存在する反射体である。これまで説明したように、送信された超音波には、保持カップ１２０を透過する際に形状の変化や時間ずれといった歪み（振幅および位相の変化）が発生する。そのように歪んだ超音波の一部は、被検体１００内に存在する反射体１０２にぶつかり反射する。
【０１４６】
図１０（ｂ）は反射体１０２で反射した歪みを伴った超音波がトランスデューサ１３１に到達して受信される状態を示している。図１０（ａ）と同様、反射された超音波には、保持カップ１２０を透過する際に再度波形の歪みや時間ずれが発生する。そのようにしてさらに歪んだ超音波がトランスデューサ１３１によって受信される。このようにして得られた受信信号から超音波エコー像を得た場合、歪みを伴った画像となる。
【０１４７】
本実施例は、そのような受信信号の歪みを補正することが可能な超音波エコー装置を提供するものである。本実施例において得られる被検体情報とは、Ｂモード画像データ、ドップラー画像データ（血流速度分布）、エラストグラフィ画像データ（歪み分布）などである。
【０１４８】
以下、図１１に示す本実施例に係る超音波エコー画像を得るためのフローチャートを用いて本実施例を説明する。
【０１４９】
（Ｓ２１０：超音波を送信する工程）
送信制御手段としてのコンピュータ１５０が、トランスデューサ１３１にパルス信号を送って所望のパルス形状の超音波を送信させる。なお、送信手段としての探触子１３０は、送信された超音波が保持カップ１２０を透過して被検体１００に照射されるように配置されている。
【０１５０】
（Ｓ２２０：超音波の反射波（エコー）を受信する工程）
探触子１３０は、Ｓ２１０で送信された超音波が反射体１０２で反射されることにより発生した超音波エコーを受信することにより電気信号を出力する。信号データ収集部１４０は、探触子１３０から出力されたアナログ信号としての電気信号をデジタル信号に変換し、コンピュータ１５０に出力する。
【０１５１】
（Ｓ２３０：音速情報を取得する工程）
音速取得手段としてのコンピュータ１５０は、Ｓ１３０と同様の方法で、被検体１００内での縦波音速ｃ１、保持カップ１２０内での縦波音速ｃ２ｌ、保持カップ１２０内での横波音速ｃ２ｔ、音響マッチング材１９０内での縦波音速ｃ３を取得する。
【０１５２】
（Ｓ２４０：入射角度情報を取得する工程）
角度取得手段としてのコンピュータ１５０は、図１２（ａ）、（ｂ）に示すように保持カップ１２０、トランスデューサ１３１を配置したときに、トランスデューサ１３１に到達する音響波が保持カップ１２０に入射したときの入射角度情報を取得する。本実施例では、コンピュータ１５０が、後に画像化されるボクセル１０１で発生した音響波が保持カップ１２０に入射するときの入射角度情報を取得する例を説明する。本実施例では、反射体１０２の存在する位置を、ボクセル１０１の位置（関心位置）として説明する。
【０１５３】
本実施例では、トランスデューサ１３１はフォーカス型であるため実際には球面波が送信されるが、計算上は、その主軸の方向にトランスデューサ１３１の中央から超音波が送信されたものと仮定する。図１２はこの状態を示したものである。図１２（ａ）は超音波の送信時を表し、図１２（ｂ）は関心位置で発生した超音波エコーの受信時を表す。この仮定は、トランスデューサの焦点距離が極端に小さいものでなければ実質的に問題ない。
【０１５４】
コンピュータ１５０は、図１２（ａ）に示す送信時について、式１で表されるスネルの法則を満たすような保持カップ１２０の第１面（被検体側の面）への入射角度θ３を取得する。図１２（ａ）では、図４（ｃ）の場合とは音響波の伝搬方向が逆であるため、入射角度はθ３となり、伝搬角はθ１となる。
【０１５５】
また、コンピュータ１５０は、図１２（ｂ）に示す受信時について、式１で示されたスネルの法則を満たすような保持カップの第１面への入射角度θ１を取得する。図１２（ｂ）では、図４（ｃ）の場合と音響波の伝搬方向が同じであるため、入射角度も同様にθ１となる。
【０１５６】
（Ｓ２５０：入射角度情報を用いて透過率フィルタを取得する工程）
透過率フィルタ取得手段としてのコンピュータ１５０は、Ｓ２４０で取得した送信時及び受信時の入射角度情報に基づいて透過率フィルタを取得する。
【０１５７】
まずコンピュータ１５０は、送信時の入射角度θ３に基づいて、送信時に対応する透過率フィルタを取得する。被検体１００の音速（縦波）をｃ１、保持カップ１２０の縦波、横波に対する音速をそれぞれｃ２Ｌ、ｃ２Ｔ、音響マッチング材１９０の音速（縦波）をｃ３、保持カップ１２０の厚さをＴとする。また、被検体１００の音響インピーダンスをＺ１、保持カップ１２０の縦波、横波に対する音響インピーダンスをそれぞれＺ２Ｌ、Ｚ２Ｔ、音響マッチング材１９０の音響インピーダンスをＺ３とする。ここで、添え字１、２、３は実施例１と同様、被検体を１、保持カップを２、音響マッチング材を３としている。
【０１５８】
保持カップ１２０中の縦波、横波に対する伝搬角をθ２Ｌ、θ２Ｔとする。これらの関係は式１−４で説明したとおりである。
【０１５９】
多重反射を考慮しない場合、角周波数ωに対する保持カップ１２０の送信方向の透過率フィルタｆ１（ω）は次のように表すことができる。ここでφ２Ｌ、φ２Ｔは式７，８と同一である。
【０１６０】
【数１４】
2017184972.tif 000016
【０１６１】
ここで、式２５後段のｅｘｐ項は、保持カップ１２０の領域に音響マッチング材１９０が存在した場合に発生する位相シフト量を補正するために設けている。
【０１６２】
多重反射を考慮する場合、角周波数ωに対する保持カップ１２０の送信方向の透過率フィルタｆ１（ω）は次のように表すことができる。ここで、Ｎ、Ｍは式１８、１９と同一である。
【０１６３】
【数１５】
2017184972.tif 000017
次に、コンピュータ１５０は、受信時の入射角度情報を用いて、受信時に対応する透過率フィルタを取得する。受信時における音響波の伝搬方向は実施例１あるいは実施例２と同じ方向である。したがって、多重反射を考慮しない場合、角周波数ωに対する保持カップ１２０の受信方向の透過率フィルタｆ２（ω）は式１１のｆ（ω）と同一である。また、多重反射を考慮する場合、角周波数ωに対する保持カップ１２０の受信方向の透過率フィルタｆ２（ω）は式２０のｆ（ω）と同一である。
【０１６４】
後述するＳ２６０でのデコンボリューション処理では、受信信号を離散フーリエ変換して得た周波数空間での受信信号に対して処理を行う。本実施例の場合、受信信号を離散フーリエ変換すると、−９９９０２３４．３７５Ｈｚから９７６５．６２５Ｈｚごとに１０００００００Ｈｚまで２０４８点の周波数情報が得られる。したがって、Ｓ２６０での処理に必要な２０４８点の角周波数に対応して透過率フィルタは計算されている。
【０１６５】
なお、本実施例では、送信時と受信時とで音速の境界が同じであると仮定して透過率フィルタを取得したが、送信時と受信時とで音速の境界を変更して透過率フィルタを取得してもよい。これにより、透過率フィルタの精度を高くすることができる。ただし、送信時と受信時とで音速の境界を同じとする場合、計算コストを小さくすることができる。
【０１６６】
（Ｓ２６０：信号データに透過率フィルタをデコンボリューションする工程）
補正手段としてのコンピュータ１５０は、Ｓ２５０で得られた送信時および受信時の透過率フィルタｆ１（ω）、ｆ２（ω）を応答関数として、Ｓ２２０で得られた信号データをデコンボリューションする。これにより補正された信号データが得られる。この計算を以下で説明する。
【０１６７】
まず、デコンボリューションフィルタＤ１（ω）、  Ｄ２（ω）を以下のように求める。
【０１６８】
【数１６】
2017184972.tif 000018
ここで、Ｃは定数であり、信号を著しく劣化させないように経験的に設定するものである。ユーザーが入力部１７０を用いて定数Ｃの値を入力してもよい。
【０１６９】
受信信号をＳ０（ｔ）として、補正された受信信号Ｓ（ｔ）は以下のように求めることができる。
【０１７０】
【数１７】
2017184972.tif 000019
このようにして得られた補正された受信信号Ｓ（ｔ）は保持カップ１２０の領域を音響マッチング材１９０で置き換えた場合とほぼ等価な受信信号となっている。
【０１７１】
（Ｓ２７０：補正された信号データを用いて被検体情報を取得する工程）
画像再構成手段としてのコンピュータ１５０は、Ｓ２６０で得られた補正された受信信号Ｓ（ｔ）を再構成処理して指定されたボクセル１０１の被検体情報を得る。本工程では、コンピュータ１５０が補正された受信信号を包絡線処理して１次元の超音波エコー画像情報を得る。なお、補正された信号データから被検体情報を得られる限り、いかなる再構成処理を適用してもよい。
【０１７２】
また、探触子１３０を機械的または電子的に走査し、Ｓ２１０−Ｓ２７０の工程を繰り返し実行することにより２次元あるいは３次元の超音波エコー画像情報を得る。
【０１７３】
コンピュータ１５０は、本工程で得られたＢモード画像データ、ドップラー画像データ、またはエラストグラフィ画像データ等の被検体情報を表示部１６０に出力する。
【０１７４】
（Ｓ２８０：被検体情報を取得する工程）
コンピュータ１５０は、Ｓ２７０で得られた被検体情報を用いて、画像化領域の被検体情報を表示部１６０に表示させる。表示部１６０は、被検体情報として、Ｂモード画像データ、ドップラー画像データ、またはエラストグラフィ画像データ等を表示することができる。表示部１６０に表示された被検体情報は、音響波が保持カップ１２０を透過することに伴う波形の歪みが抑制された情報であるため、医師等の作業者が診断等に用いる上で好適な情報となっている。
【０１７５】
本実施例では、超音波が保持カップを透過する際の波形の歪み（振幅、位相）を補正することで、高精細な超音波エコー画像を得ることが可能となる。
【実施例４】
【０１７６】
図１３に示す光音響装置を用いて本実施例を説明する。図１３において既出の構成と同様の構成については同じ符号を付し、詳細な説明を省略する。被検体情報取得装置において、被検体情報を高精度に得るためにはトランスデューサ１３１の取り付け位置を正確に把握しておくことが重要である。本実施例では、上述のＳ１６０で取得した補正された電気信号を用いて、トランスデューサ１３１の取り付け位置の情報を取得する方法について説明する。
【０１７７】
本実施例では、音源５１３を保持カップ１２０中の所定の位置に配置する。音源５１３には０．３ｍｍのゴム球を用いた。音源５１３は、トランスデューサ１３１が配置された球面状の支持体１３２の球面中心に配置されることが望ましい。また、音源５１３の位置は、球面中心でなくてもよいがその座標が既知であることが望ましい。音源５１３の周りは水５１４で満たしておく。音源５１３の材質はレーザー光の波長に対して吸収するものであれば何でもよい。また、音源５１３は理想的には球体であることが望ましい。球体でなくても、その大きさ（縦幅、横幅、高さの最大値）がトランスデューサ１３１の検出帯域の中心周波数を有する音響波の水中での波長以下であればよい。この場合、音源から発生する音響波の波形は球体から発生する音響波の波形とほぼ等しくなるため、実質的には音源を球体とみなすことができる。すなわち、トランスデューサ１３１の受信帯域の中心周波数は、当該中心周波数を有する音響波の水中での波長が音源５１３のサイズの最大値よりも大きくなるような中心周波数としてもよい。
【０１７８】
ところで、音源が大きい場合、照射された光は音源の表面で吸収され、音源の内部には届かない場合がある。したがって、音源は、その表面だけ光吸収作用を有する中空の光吸収体として扱うことができる。音源から発生した音響波はインパルス関数に近い形状を有する。例えば、鉄球などの金属球を用いた場合がこの場合に相当する。このとき、音源から発生して所望の方向に伝搬する音響波と、それとは逆方向に伝搬して音源の裏面で反射した音響波とが干渉して探触子で検出された場合、時間的な誤差が生じる。したがって、音源の裏面で反射した音響波を時間的に分離できるよう、音源の直径は、トランスデューサの受信帯域の中心周波数を有する音響波の水中での波長の５倍以上であることが望ましい。トランスデューサ１３１の受信帯域の中心周波数は、当該中心周波数を有する音響波の水中での波長が音源５１３のサイズの最大値の１／５となるような中心周波数としてもよい。
【０１７９】
図１４は、配置位置が既知である音源５１３からの音響波を利用して、トランスデューサ１３１の取り付け位置情報を取得するフローチャートである。図１４は、上述したＳ１６０で補正された電気信号を取得した後の処理フローを示したものである。
【０１８０】
（Ｓ１７１：補正された電気信号を用いて音響波の到達時間を取得する工程）
コンピュータ１５０は、Ｓ１６０で得られた補正された電気信号を用いて、音源５１３から発生した音響波がトランスデューサ１３１に到達した到達時間を計算する。この計算を以下で説明する。
【０１８１】
図１５は、Ｓ１６０で得られる、補正された受信信号の例である。横軸はサンプリングナンバー、縦軸は信号強度を表す。図１５のように、補正された受信信号は、信号強度の最大点と最小点の中心位置を中心として点対称な形状をしている。そこで、最大点と最小点の間で値がゼロとなるサンプリングナンバーを到達時間の基準とする。サンプリングナンバーは離散的な整数であるが、より精度を高めるため、Ｘ軸と交わる前後の値から線形補間して小数点以下の数値も算出している。もちろん、補間の方法はこれに限ったものではない。ここで求めたサンプリングナンバーをＮｉとする。ここで、ｉはトランスデューサ１３１の素子ナンバーである。
【０１８２】
実施例１で述べたように、信号データ収集部１４０では、音源５１３にパルス光を照射した瞬間の信号をサンプル０番として、サンプル１０００番から３０４７番まで２０４８サンプルの信号を取得している。サンプリング周波数をｆｓａｍｐとすると、音響波の到達時間ｔｉは次式で表すことができる。
【０１８３】
【数１８】
2017184972.tif 000020
（Ｓ１７２：音響波の到達時間に基づいてトランスデューサの位置補正量を取得する工程）
ずれ量推定部としてのコンピュータ１５０は、Ｓ１７１で取得した音響波の到達時間ｔｉに基づいてトランスデューサ１３１の位置補正量を計算する。すなわち、コンピュータ１５０は、トランスデューサ１３１の設計位置（所定の位置）からのずれ量を位置補正量として計算する。この計算を以下で説明する。
【０１８４】
音源５１３の位置を（Ｘ０，Ｙ０，Ｚ０）とし、トランスデューサ１３１−ｉの設計上の位置を（ｘｉ，ｙｉ，ｚｉ）とする。位置補正量ΔＲｉは次式で表すことができる。
【０１８５】
【数１９】
2017184972.tif 000021
ここで、ｃは水の音速である。位置補正量ΔＲｉは設計上の距離と実際の距離との差を意味する。位置補正量ΔＲｉはコンピュータ１５０の記憶手段に記憶しておき、後述する画像再構成を実行するときに用いる。
【０１８６】
（Ｓ１７３：トランスデューサの位置補正量を用いて被検体情報を取得する工程）
コンピュータ１５０は、Ｓ１７２で取得したトランスデューサ１３１の位置補正量を用いてトランスデューサ１３１の位置情報を取得する。実際に、画像再構成によって画像を生成するとき、位置補正量ΔＲｉは次のように用いられる。
【０１８７】
トランスデューサ１３１が配置されている半球面の支持体１３２の設計上の半径（実施例１では１２７ｍｍ）をＲとし、トランスデューサ１３１−ｉの設計上の位置は、半球面の中心を原点として、球座標系で、
【０１８８】
【数２０】
2017184972.tif 000022
で表されているものとする。このとき、コンピュータ１５０は、トランスデューサ１３１−ｉの実際の位置を、
【０１８９】
【数２１】
2017184972.tif 000023
であるとして置き換える。すなわち、距離の誤差を半径方向のずれとして補正を行う。本実施例では、支持体１３２にトランスデューサ１３１を取り付ける穴をあけてそこにトランスデューサ１３１を差し込んで固定する、といった作製工程を考慮している。式３３が表す誤差補正は、この作製工程において、半径方向が最も発生しやすいずれの方向であることに対応させたものである。
【０１９０】
以上が音源５１３の位置が既知の場合のトランスデューサ１３１の位置情報を取得する方法である。また、上記方法によれば、保持カップと透過する際に発生する音響波の歪みが補正された信号を用いてトランスデューサの位置情報を精度よく取得することができる。
【０１９１】
ところで、上記方法によれば、音源５１３の位置決めの精度がトランスデューサ１３１の位置情報の精度に依存する。そこで、音源５１３の位置が正確に把握できていない場合のトランスデューサ１３１の位置情報を取得する方法について以下に説明する。すなわち、音源５１３の位置が既知ではなく、未知である場合にトランスデューサ１３１の位置情報を取得する方法を説明する。図１６は、音源５１３の位置が未知である場合のトランスデューサ１３１の位置情報取得のフローを示す。
【０１９２】
（Ｓ４１０：音源の位置を仮設定する工程）
コンピュータ１５０は、音源５１３の位置を仮に定める。コンピュータ１５０は、（ｘｐ、ｙｐ、ｚｐ）を音源５１３の位置として仮設定する。Ｓ１４０、Ｓ１５０、Ｓ１６０、Ｓ１７１では、本工程で仮設定された位置（ｘｐ、ｙｐ、ｚｐ）に音源５１３が配置されていると仮定して処理が実行される。
【０１９３】
（Ｓ１７５：補正された信号から求められる距離と、幾何学的に求められる距離との差を表すエラー関数を取得する工程）
コンピュータ１５０は、Ｓ１７１で得られた到達時間に基づいて、音源５１３とトランスデューサ１３１との距離を取得する。また、コンピュータ１５０は、Ｓ４１０で仮設定された音源５１３の位置とトランスデューサ１３１の位置の設計値とに基づいて、音源５１３とトランスデューサ１３１との距離を取得する。そして、コンピュータ１５０は、これらの距離の差を表すエラー関数を取得する。具体的には、次式によってエラー関数Ｅ（ｖ，ｘ，ｙ，ｚ）を定義する。
【０１９４】
【数２２】
2017184972.tif 000024
ここで、ｖは水の音速を示す変数であり、定数ｃと区別するため文字を変えている。ｔｉは式３０で表される音響波の到達時間である。また、ｘ、ｙ、ｚは音源５１３の位置を示す変数である。
【０１９５】
（Ｓ１７６：エラー関数を最小とする音源の位置を取得する工程）
式３４のエラー関数Ｅを最小にするｖ、ｘ、ｙ、ｚの値を求める。これは、
【０１９６】
【数２３】
2017184972.tif 000025
を満たすｖ、ｘ、ｙ、ｚを求めればよい。このとき得られた解を、ｖａ、ｘａ、ｙａ、ｚａとする。
【０１９７】
（Ｓ１７７：収束したか否かを判定する工程）
コンピュータ１５０は、Ｓ１７５で得られた音源５１３の位置（ｘａ、ｙａ、ｚａ）が収束値であるか否かを判定する。
【０１９８】
収束したと判定された場合はＳ１７２に進む。収束したと判定されたときの解を、Ｖａ、Ｘａ、Ｙａ、Ｚａとする。Ｓ１７２において、コンピュータ１５０は、収束したときにＳ１７１で取得した音響波の到達時間ｔｉに基づいてトランスデューサ１３１の位置補正量を計算する。位置補正量ΔＲｉは次式で表すことができる。
【０１９９】
【数２４】
2017184972.tif 000026
ここで、音源５１３の位置を（Ｘａ，Ｙａ，Ｚａ）とし、トランスデューサ１３１−ｉの設計上の位置を（ｘｉ，ｙｉ，ｚｉ）とする。
【０２００】
収束していないと判定された場合は、Ｓ４１０に進む。２回目以降のＳ４１０では、Ｓ１７６で取得された音源５１３の位置（ｘａ、ｙａ、ｚａ）を音源の位置として仮設定する。そして、Ｓ１４０、Ｓ１５０、Ｓ１６０、Ｓ１７１では、仮設定された位置（ｘａ、ｙａ、ｚａ）に音源５１３が配置されていると仮定して処理が実行される。
【０２０１】
例えば、コンピュータ１５０は、Ｓ１７６で得られた音源５１３の位置（ｘａ、ｙａ、ｚａ）と、Ｓ４１０で仮設定された位置（ｘｐ、ｙｐ、ｚｐ）との距離が閾値以下となるときに収束したと判定してもよい。また、コンピュータ１５０は、Ｓ１７６で得られた音源５１３の位置（ｘａ、ｙａ、ｚａ）と、ひとつ前のループで得られた音源５１３の位置（ｘａ、ｙａ、ｚａ）との距離が閾値以下となるときに収束したと判定してもよい。また、コンピュータ１５０は、Ｓ４１０〜Ｓ１７６の工程を所定の回数繰り返したときに収束したと判定してもよい。距離に対する閾値や繰り返し回数については、ユーザーが入力部１７０を用いて指定してもよい。
【０２０２】
なお、本実施例においては、音速ｖを変数としたが、変数ではなく、一定としてもよい。また、発光タイミング（音響波の発生タイミング）とサンプリング開始タイミングとの時間的なずれを調整するために、Ｔを変数とした次式のエラー関数を定義してもよい。
【０２０３】
【数２５】
2017184972.tif 000027
また、球体の比較的大きな音源５１３を用いた場合、音源５１３の中央ではなく表面から音響波は発生する。よって、位置補正量ΔＲｉの算出の際は、その半径分の伝搬時間を考慮する必要がある。音源５１３の半径をｒ０、音源５１３の中心位置を（Ｘ０，Ｙ０，Ｚ０）とすると、位置補正量ΔＲｉは次式で表すことができる。
【０２０４】
【数２６】
2017184972.tif 000028
以上が、音源５１３の位置が未知である場合のトランスデューサ１３１の位置情報を取得する方法である。本実施例によれば、保持カップと透過する際に発生する音響波の歪みが補正された信号を用いてトランスデューサの位置情報を精度よく取得することができる。
【０２０５】
実際に本実施例の方法により位置補正量を取得した例を図１７に示す。図１７は、音源５１３として直径０．３ｍｍのゴム球を用いた場合に算出された位置補正量を示している。横軸は、トランスデューサ１３１の素子番号であり、５１２個あるうちの１００番から２００番までを示している。縦軸は、その素子番号での位置補正量を示している。
【０２０６】
本実施例では、音源５１３として０．３ｍｍのゴム球を用いた場合を示したがそれに限ったものではない。例えば、音源としては鉄球などの金属球を用いることもできる。
【実施例５】
【０２０７】
実施例１〜４で説明した情報取得装置において、高精細な被検体情報を得るためには保持カップ１２０の厚さを正確に把握しておくことが重要である。例えば、樹脂材料を保持カップに適した形状に成形する方法として、真空成形や圧空成形がある。どちらもシート状の樹脂材料を加熱して柔らかくして型に当てて冷却する、という方法である。このような成形プロセスによって保持カップの膜厚に分布ができる場合がある。これは、加熱温度、型に当てる速度など様々な条件に依存し、周辺よりも中央が厚くなるケースも、逆に周辺よりも中央が薄くなるケースもある。均一な膜厚分布とするには成形方法の条件の最適化が必要でありコストが高くなる要因となる。不均一な膜厚分布を有する保持カップを用いた場合でも、予めその膜厚分布を把握しておけば、実施例１から４で述べた受信信号の補正を精度よく実施することが可能となる。しかし、保持カップ１２０のような曲面の厚さを非破壊で測定するには特殊な装置が必要となる。
【０２０８】
そこで、本実施例では、図１３に示す光音響装置を用いて、破壊することなく保持カップ１２０の厚さに関する情報を取得する方法について説明する。
【０２０９】
まず、図１３に示すように保持カップ１２０が設置された状態でＳ１１０、Ｓ１２０の工程を実行し、受信信号データをコンピュータ１５０に保存する。保持カップ１２０を設置した状態で得られた受信信号データを第１の受信信号データと呼ぶ。
【０２１０】
続いて、保持カップ１２０を外した状態でＳ１１０、Ｓ１２０の工程を実行し、受信信号データをコンピュータ１５０に保存する。保持カップ１２０を外した状態で得られた受信信号データを第２の受信信号データと呼ぶ。なお、保持カップ１２０を外した測定を先に行ってもよい。
【０２１１】
続いて、コンピュータ１５０は、保持カップ１２０の厚さに関するデータを仮で設定する。厚さに関するデータは均一の厚さを示すデータでもよいし、厚さに分布を持たせたデータであってもよい。
【０２１２】
続いて、コンピュータ１５０は、仮設定された保持カップ１２０の厚さに関するデータを用いて、式１１や式２０にしたがって透過率フィルタを取得する。そして、コンピュータ１５０は、Ｓ１６０で説明した方法と同様に、第１の受信信号データに対して、得られた透過率フィルタをデコンボリューションすることにより、補正された受信信号データを取得する。ここで補正された受信信号データを第３の受信信号データと呼ぶ。
【０２１３】
続いて、コンピュータ１５０は、Ｓ１７１で説明した方法と同様に、第２の受信信号データを用いて保持カップ１２０を外したときの音響波の到達時刻ｔ２を取得する。また、コンピュータ１５０は、同様に、第３の受信信号データを用いて保持カップ１２０を外したと仮定したときの音響波の到達時刻ｔ３を取得する。
【０２１４】
続いて、コンピュータ１５０は、第２の受信信号データに基づいて得られた到達時刻ｔ２と第３の受信信号データに基づいて得られた到達時刻ｔ３との差Δｔを算出する。そして、コンピュータ１５０は、保持カップ１２０の厚さに関するデータをアップデートし、上記の処理を再度行い、Δｔが０または閾値よりも小さくなるまで繰り返す。コンピュータ１５０は、Δｔが０または閾値よりも小さくなったときの保持カップ１２０の厚さに関するデータを真値として取得する。なお、保持カップ１２０の厚さに関するデータのアップデートを所定の回数行って得られたデータを真値としてもよい。
【０２１５】
実際に本実施例の方法により保持カップ１２０の測定した結果を図１８に示す。横軸は保持カップ１２０の中心からの距離、縦軸は保持カップ１２０の厚さである。本実施例では、保持カップ１２０として、厚さが中心からの距離に依存して分布があるような保持カップを用いている。図中、○は実測点、実線はフィッティングカーブである。ここで、図１９に示す保持カップ１２０の底部ＯからＲ方向（Ｒ＝√（Ｘ２＋Ｙ２））の距離ｄを、保持カップ１２０の中心からの距離としている。また、保持カップ１２０の表面に対する法線方向への厚さを、当該距離ｄにおける保持カップ１２０の厚さＴ（ｄ）としている。なお、本実施例において保持カップ１２０の位置については既知であるとして計算している。すなわち、本実施例においては保持カップ１２０の音源５１３側の面の位置が既知であるとして計算している。
【０２１６】
この実施例では、保持カップ１２０の中心からの距離をｄ（単位ｍｍ）として、厚さＴ（ｄ）（単位ｍｍ）は、下記のようになった。
【０２１７】
【数２７】
2017184972.tif 000029
この厚さに関する情報をコンピュータ１５０内の記憶手段に記憶する。なお、本実施例ではＸＹ方向に向かうにつれて保持カップの厚さが変化する場合を想定して計算したが、いかなる保持カップの厚さの変化に対応してもよい。
【０２１８】
本実施例によれば、保持カップと透過する際に音響波の歪みが発生する場合であっても、保持カップの厚さに関する情報を精度良く取得することができる。また、本実施例によれば、保持カップの厚さが未知である場合でも、精度良く透過率フィルタを取得することができる。
【符号の説明】
【０２１９】
１２０  保持カップ
１３０  探触子
１５０  コンピュータ
【図１】
2017184972.tif 000030
【図２】
2017184972.tif 000031
【図３】
2017184972.tif 000032
【図４】
2017184972.tif 000033
【図５】
2017184972.tif 000034
【図６】
2017184972.tif 000035
【図７】
2017184972.tif 000036
【図８】
2017184972.tif 000037
【図９】
2017184972.tif 000038
【図１０】
2017184972.tif 000039
【図１１】
2017184972.tif 000040
【図１２】
2017184972.tif 000041
【図１３】
2017184972.tif 000042
【図１４】
2017184972.tif 000043
【図１５】
2017184972.tif 000044
【図１６】
2017184972.tif 000045
【図１７】
2017184972.tif 000046
【図１８】
2017184972.tif 000047
【図１９】
2017184972.tif 000048