(57)【要約】【課題】音響波を受信する装置において、被検体の変位の影響を低減する。【解決手段】被検体から発生する音響波を受信して信号を出力する探触子と、被検体と探触子の相対位置を変化させることにより、探触子に被検体を走査させる走査手段と、被検体に対する音響波測定の実行を制御する制御手段と、信号を用いて被検体の特性情報を取得する情報取得手段を備え、制御手段は、被検体の所定の領域について、第一の音響波測定および第二の音響波測定を実行するものであり、第一の音響波測定は第二の音響波測定よりも短時間の間に実行されるものであり、情報取得手段は、第一の音響波測定で取得された第一の信号および第二の音響波測定で取得された第二の信号を用いて、被検体の変位を示す変位情報を取得する被検体情報取得装置を用いる。【選択図】図３
2020039809.tif 000002
【特許請求の範囲】
【請求項１】
被検体から発生する音響波を受信して信号を出力する探触子と、
前記被検体と前記探触子の相対位置を変化させることにより、前記探触子に前記被検体を走査させる走査手段と、
前記被検体に対する音響波測定の実行を制御する制御手段と、
前記信号を用いて前記被検体の特性情報を取得する情報取得手段と、
を備え、
前記制御手段は、前記被検体の所定の領域について、第一の音響波測定および第二の音響波測定を実行するものであり、前記第一の音響波測定は前記第二の音響波測定よりも短時間の間に実行されるものであり、
前記情報取得手段は、前記第一の音響波測定で取得された第一の信号および前記第二の音響波測定で取得された第二の信号を用いて、前記被検体の変位を示す変位情報を取得する
ことを特徴とする被検体情報取得装置。
【請求項２】
前記情報取得手段は、前記第一の信号を用いて第一の被検体形状情報を取得し、前記第二の信号を用いて第二の被検体形状情報を取得し、前記第一の被検体形状情報および前記第二の被検体形状情報を用いて前記変位情報を取得する
ことを特徴とする請求項１に記載の被検体情報取得装置。
【請求項３】
前記情報取得手段は、前記第二の信号を用いて前記特性情報を取得し、当該特性情報を、前記変位情報を用いて補正する
ことを特徴とする請求項２に記載の被検体情報取得装置。
【請求項４】
前記制御手段は、前記第一の音響波測定での測定位置の数が、前記第二の音響波測定での測定位置の数よりも少なくなるように制御を行う
ことを特徴とする請求項１ないし３のいずれか１項に記載の被検体情報取得装置。
【請求項５】
前記制御手段は、前記第一の音響波測定における複数の測定位置が、それぞれ、前記第二の音響波測定における複数の測定位置のいずれかに重なるような制御を行う
ことを特徴とする請求項４に記載の被検体情報取得装置。
【請求項６】
前記制御手段は、前記第一の音響波測定における複数の測定位置により囲まれる領域が、前記第二の音響波測定における複数の測定位置により囲まれる領域を包含するように制御を行う
ことを特徴とする請求項４または５に記載の被検体情報取得装置。
【請求項７】
前記走査手段は、
前記第二の音響波測定の際に、前記探触子をラスタ軌道または曲線を含む軌道にてさせるものであり、
前記第一の音響波測定の際に、前記探触子を前記ラスタ軌道または前記曲線を含む軌道と交差する方向に移動させる
ことを特徴とする請求項１ないし６のいずれか１項に記載の被検体情報取得装置。
【請求項８】
前記情報取得手段は、前記第一の音響波測定での複数の測定位置において前記被検体の表面の位置を取得し、前記第二の音響波測定での複数の測定位置において前記被検体の表面の位置を取得し、前記第一の音響波測定で取得された前記表面の位置と前記第二の音響波測定で取得された前記表面の位置に基づいて、前記変位情報を取得する
ことを特徴とする請求項２に記載の被検体情報取得装置。
【請求項９】
前記情報取得手段は、前記第一の音響波測定での測定位置と前記第二の音響波測定での測定位置が一致しない場合、補間処理を行ったのち前記変位情報を取得する
ことを特徴とする請求項１ないし８のいずれか１項に記載の被検体情報取得装置。
【請求項１０】
前記被検体にパルス光を照射して前記音響波を発生させる照射手段をさらに有する
ことを特徴とする請求項１ないし９のいずれか１項に記載の被検体情報取得装置。
【請求項１１】
前記第一の音響波測定における前記パルス光の繰り返し周波数は、前記第二の音響波測定における前記パルス光の繰り返し周波数よりも高い
ことを特徴とする請求項１０に記載の被検体情報取得装置。
【請求項１２】
前記第一の音響波測定における前記パルス光のエネルギーは、前記第二の音響波測定における前記パルス光のエネルギーよりも小さい
ことを特徴とする請求項１０または１１に記載の被検体情報取得装置。
【請求項１３】
前記音響波は、前記探触子から前記被検体に送信されたのち反射したものである
ことを特徴とする請求項１ないし９のいずれか１項に記載の被検体情報取得装置。
【請求項１４】
前記特性情報に基づく画像を表示する表示手段をさらに有し、
前記表示手段は、前記補正の前の画像と前記補正の後の画像を、並列に、または切り替え可能に表示する
ことを特徴とする請求項３に記載の被検体情報取得装置。
【請求項１５】
探触子と、走査手段と、制御手段と、情報取得手段と、を備える被検体情報取得装置の制御方法であって、
前記探触子が、被検体から発生する音響波を受信して信号を出力する出力ステップと、
前記走査手段が、前記被検体と前記探触子の相対位置を変化させることにより、前記探触子に前記被検体を走査させる走査ステップと、
前記制御手段が、前記被検体に対する音響波測定の実行を制御する制御ステップと、
前記情報取得手段が、前記信号を用いて前記被検体の特性情報を取得する情報取得ステップと、
を備え、
前記制御ステップは、前記被検体の所定の領域について、第一の音響波測定および第二の音響波測定を実行するものであり、前記第一の音響波測定は前記第二の音響波測定よりも短時間の間に実行されるものであり、
前記情報取得ステップは、前記第一の音響波測定で取得された第一の信号および前記第二の音響波測定で取得された第二の信号を用いて、前記被検体の変位を示す変位情報を取得する
ことを特徴とする被検体情報取得装置の制御方法。
【請求項１６】
前記情報取得ステップは、前記第一の信号を用いて第一の被検体形状情報を取得し、前記第二の信号を用いて第二の被検体形状情報を取得し、前記第一の被検体形状情報および前記第二の被検体形状情報を用いて前記変位情報を取得する
ことを特徴とする請求項１５に記載の被検体情報取得装置の制御方法。
【請求項１７】
前記情報取得ステップは、前記第二の信号を用いて前記特性情報を取得し、当該特性情報を、前記変位情報を用いて補正する
ことを特徴とする請求項１６に記載の被検体情報取得装置の制御方法。
【請求項１８】
前記制御ステップは、前記第一の音響波測定での測定位置の数が、前記第二の音響波測定での測定位置の数よりも少なくなるように制御を行う
ことを特徴とする請求項１５ないし１７のいずれか１項に記載の被検体情報取得装置の制御方法。
【請求項１９】
前記制御ステップは、前記第一の音響波測定における複数の測定位置が、それぞれ、前記第二の音響波測定における複数の測定位置のいずれかに重なるような制御を行う
ことを特徴とする請求項１８に記載の被検体情報取得装置の制御方法。
【請求項２０】
前記制御ステップは、前記第一の音響波測定における複数の測定位置により囲まれる領域が、前記第二の音響波測定における複数の測定位置により囲まれる領域を包含するように制御を行う
ことを特徴とする請求項１８または１９に記載の被検体情報取得装置の制御方法。
【請求項２１】
前記走査ステップは、
前記第二の音響波測定の際に、前記探触子をラスタ軌道または曲線を含む軌道にてさせるものであり、
前記第一の音響波測定の際に、前記探触子を前記ラスタ軌道または前記曲線を含む軌道と交差する方向に移動させる
ことを特徴とする請求項１５ないし２０のいずれか１項に記載の被検体情報取得装置の制御方法。
【請求項２２】
前記情報取得ステップは、前記第一の音響波測定での複数の測定位置において前記被検体の表面の位置を取得し、前記第二の音響波測定での複数の測定位置において前記被検体の表面の位置を取得し、前記第一の音響波測定で取得された前記表面の位置と前記第二の音響波測定で取得された前記表面の位置に基づいて、前記変位情報を取得する
ことを特徴とする請求項１６に記載の被検体情報取得装置の制御方法。
【請求項２３】
前記情報取得ステップは、前記第一の音響波測定での測定位置と前記第二の音響波測定での測定位置が一致しない場合、補間処理を行ったのち前記変位情報を取得する
ことを特徴とする請求項１５ないし２２のいずれか１項に記載の被検体情報取得装置の制御方法。
【請求項２４】
前記被検体情報取得装置は、照射手段をさらに備え、
前記照射手段が、前記被検体にパルス光を照射して前記音響波を発生させる照射ステップをさらに有する
ことを特徴とする請求項１５ないし２３のいずれか１項に記載の被検体情報取得装置の制御方法。
【請求項２５】
前記第一の音響波測定における前記パルス光の繰り返し周波数は、前記第二の音響波測定における前記パルス光の繰り返し周波数よりも高い
ことを特徴とする請求項２４に記載の被検体情報取得装置の制御方法。
【請求項２６】
前記第一の音響波測定における前記パルス光のエネルギーは、前記第二の音響波測定における前記パルス光のエネルギーよりも小さい
ことを特徴とする請求項２４または２５に記載の被検体情報取得装置の制御方法。
【請求項２７】
前記音響波は、前記探触子から前記被検体に送信されたのち反射したものである
ことを特徴とする請求項１５ないし２３のいずれか１項に記載の被検体情報取得装置の制御方法。
【請求項２８】
前記被検体情報取得装置は、表示手段をさらに備え、
前記表示手段が、前記特性情報に基づく画像を表示する表示ステップをさらに有し、
前記表示ステップは、前記補正の前の画像と前記補正の後の画像を、並列に、または切り替え可能に表示する
ことを特徴とする請求項１７に記載の被検体情報取得装置の制御方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被検体情報取得装置およびその制御方法に関する。
【背景技術】
【０００２】
音響波を受信して生体などの被検体の情報を取得する技術として、光音響イメージング装置や超音波エコーイメージング装置などの被検体情報取得装置が提案されている。
光音響イメージング装置は、特に皮膚疾患や乳がんの診断での有用性が示されている。可視光や近赤外光等を被検体組織に照射すると、被検体内部の光吸収物質（例えば血液中のヘモグロビン等）が光エネルギーを吸収して瞬間的に膨張し、音響波を発生させる。この現象を光音響効果といい、発生した音響波を光音響波とも呼ぶ。光音響イメージング装置は、この光音響波を計測することで被検体組織の情報を可視化する。
【０００３】
光音響効果を利用した断層撮影の技術を光音響イメージング（Ｐｈｏｔｏａｃｏｕｓｔｉｃ  Ｉｍａｇｉｎｇ：ＰＡＩ）と呼ぶ。光音響イメージングは、被検体内部の吸収係数に関連した情報を画像化できる。吸収係数とは被検体組織が光エネルギーを吸収する率である。吸収係数に関連した情報として、例えば、光音響波が発生した瞬間の音圧である初期音圧がある。初期音圧は光エネルギー（光強度）と吸収係数の積に比例する。さらに、吸収係数は被検体組織を構成する成分の濃度に依存する。したがって、吸収係数からそれら成分の濃度を取得できる。また、濃度情報から酸素飽和度等の情報を得ることもできる。これら情報を分析することで、被検体内の腫瘍組織や、その周辺組織を判別するなど、医療診断への応用が期待されている。
【０００４】
また、超音波エコーイメージング装置は、被検体内部の音響インピーダンスの差を可視化して形態情報を取得できるため、診断を含めた様々な分野で用いられている。
【０００５】
探触子やトランスデューサを走査するタイプの被検体情報取得装置では、走査軌道上の各走査点（測定位置）において取得された信号を使って画像再構成を行う。このような走査型の被検体情報取得装置の課題として、走査中の被検体の変位がある。すなわち、走査中に被検体の体動などにより変位が起こると、被検体情報の精度が低下する虞がある。
【０００６】
特許文献１では、超音波プローブの位置ずれを求める光音響計測装置が開示されている。また、非特許文献１では、光音響イメージング装置において被検体の体動を補正する方法が開示されている。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２０１４−０６１１２４号公報
【非特許文献】
【０００８】
【非特許文献１】“Ｍｏｔｉｏｎ  ｃｏｒｒｅｃｔｉｏｎ  ｉｎ  ｏｐｔｏａｃｏｕｓｔｉｃ  ｍｅｓｏｓｃｏｐｙ”，Ｓｃｉｅｎｔｉｆｉｃ  Ｒｅｐｏｒｔｓ  ７，  ａｒｔｉｃｌｅ  ｎｕｍｂｅｒ：１０３８６（２０１７）
【非特許文献２】“Ｕｎｉｖｅｒｓａｌ  ｂａｃｋ−ｐｒｏｊｅｃｔｉｏｎ  ａｌｇｏｒｉｔｈｍ  ｆｏｒ  ｐｈｏｔｏａｃｏｕｓｔｉｃ  ｃｏｍｐｕｔｅｄ  ｔｏｍｏｇｒａｐｈｙ”，  ＰＨＹＳＩＣＡＬ  ＲＥＶＩＥＷ  Ｅ  ７１，  ０１６７０６  （２００５）
【発明の概要】
【発明が解決しようとする課題】
【０００９】
特許文献１に記載の光音響計測装置は、プローブの位置ずれを求めるために、パターンが形成された部材を備えている。しかし、この部材を体動の検知に用いた場合、部材上のパターンによりアーティファクトが発生することや、部材の音速が被検体と異なることが原因となり、画質の劣化が生じる虞がある。また、非特許文献１に記載の方法では、全走査時間に対して周期が十分短い短周期の体動は補正できるが、全走査時間にわたって起きる長周期の体動を補正することは困難であった。
【００１０】
本発明は、上記課題に鑑みてなされたものである。本発明の目的は、音響波を受信する装置において、被検体の変位の影響を低減することにある。
【課題を解決するための手段】
【００１１】
本発明は、以下の構成を採用する。すなわち、
被検体から発生する音響波を受信して信号を出力する探触子と、
前記探触子を走査して前記被検体と前記探触子の相対位置を変化させる走査手段と、
前記探触子が前記音響波を受信するときの測定位置を制御することにより、前記被検体に対する音響波測定を実行する制御手段と、
前記信号を用いて前記被検体の特性情報を取得する情報取得手段と、
を備え、
前記制御手段は、第一の音響波測定および第二の音響波測定を実行するものであり、前記第一の音響波測定は前記第二の音響波測定よりも短時間の間に行われるものであり、
前記情報取得手段は、前記第一の音響波測定で取得された第一の信号および前記第二の音響波測定で取得された第二の信号を用いて、前記被検体の変位を示す変位情報を取得する
ことを特徴とする被検体情報取得装置である。
本発明はまた、以下の構成を採用する。すなわち、
探触子と、走査手段と、制御手段と、情報取得手段と、を備える被検体情報取得装置の制御方法であって、
前記探触子が、被検体から発生する音響波を受信して信号を出力する出力ステップと、
前記走査手段が、前記被検体と前記探触子の相対位置を変化させることにより、前記探触子に前記被検体を走査させる走査ステップと、
前記制御手段が、前記被検体に対する音響波測定の実行を制御する制御ステップと、
前記情報取得手段が、前記信号を用いて前記被検体の特性情報を取得する情報取得ステップと、
を備え、
前記制御ステップは、前記被検体の所定の領域について、第一の音響波測定および第二の音響波測定を実行するものであり、前記第一の音響波測定は前記第二の音響波測定よりも短時間の間に実行されるものであり、
前記情報取得ステップは、前記第一の音響波測定で取得された第一の信号および前記第二の音響波測定で取得された第二の信号を用いて、前記被検体の変位を示す変位情報を取得する
ことを特徴とする被検体情報取得装置の制御方法である。
【発明の効果】
【００１２】
本発明によれば、音響波を受信する装置において、被検体の変位の影響を低減できる。
【図面の簡単な説明】
【００１３】
【図１】第一の実施形態の被検体情報取得装置の模式図および処理部周辺模式図。
【図２Ａ】第一の実施形態の全体的な処理を示すフロー図。
【図２Ｂ】第一の実施形態の被検体情報取得方法のフロー図。
【図３】第一の実施形態の第一および第二の音響波測定の典型例を示した図。
【図４】第一の実施形態の体動補正量算出方法の模式図。
【図５】体動補正方法の模式図。
【図６】表示部の例を示した図。
【図７】体動量が許容閾値を超えた場合を示した図。
【図８】第一の音響波測定の軌道を例示した図。
【図９】第二の実施形態の被検体情報取得方法のフロー図。
【図１０】第二の実施形態の第一および第二の音響波測定の典型例を示した図。
【図１１】第二の実施形態の体動補正量算出方法の模式図。
【図１２】第三の実施形態の被検体情報取得装置の模式図および処理部周辺模式図。
【図１３】第三の実施形態の被検体情報取得方法のフロー図。
【図１４】第三の実施形態の第一および第二の音響波測定の典型例を示した図。
【図１５】被検体断面と到達光量の模式図。
【図１６】光源のパルス周波数とパルスエネルギーの関係の模式図。
【図１７】被検体情報取得装置の別の構成例を示す模式図。
【発明を実施するための形態】
【００１４】
以下に図面を参照しつつ、本発明の好適な実施の形態について説明する。ただし、以下に記載されている構成部品の寸法、材質、形状およびそれらの相対配置などは、発明が適用される装置の構成や各種条件により適宜変更されるべきものである。よって、この発明の範囲を以下の記載に限定する趣旨のものではない。
【００１５】
本発明は、被検体に光（電磁波）を照射し、被検体内部から発生した音響波を用いて被検体内部の特性情報（被検体情報）を生成し、取得する技術に関する。よって本発明は、光音響装置またはその制御方法、被検体情報取得装置またはその制御方法、あるいは被検体情報取得方法や信号処理方法として捉えられる。本発明はまた、これらの方法をＣＰＵやメモリ等のハードウェア資源を備える情報処理装置に実行させるプログラムや、そのプログラムを格納した、コンピュータにより読み取り可能な非一時的な記憶媒体としても捉えられる。
【００１６】
また、本発明の被検体情報取得装置は、被検体に音響波を照射し、被検体内の特定位置で反射、散乱されて伝搬してきた音響波を受信（検出）する、エコー技術を利用した装置を含む。このような被検体情報取得装置は、音響波の反射、散乱特性に基づく被検体内部の特性情報を画像データ等の形式で得ることから、超音波エコーイメージング装置とも呼べる。この場合、本発明は、超音波エコーイメージング装置またはその制御方法、被検体情報取得装置またはその制御方法、あるいは被検体情報取得方法や信号処理方法として捉えられる。本発明はまた、これらの方法をＣＰＵやメモリ等のハードウェア資源を備える情報処理装置に実行させるプログラムや、そのプログラムを格納した、コンピュータにより読み取り可能な非一時的な記憶媒体としても捉えられる。
【００１７】
光音響装置における特性情報とは、光音響波に由来する受信信号を用いて生成される、被検体内の複数位置のそれぞれに対応する、光エネルギーの吸収量や吸収率を反映した値である。この特性情報には例えば、単一の波長の光照射によって生じた音響波の発生源、被検体内の初期音圧、あるいは初期音圧から導かれる光エネルギー吸収密度や吸収係数を含む。また、互いに異なる複数の波長により得られる特性情報から、組織を構成する物質の濃度を取得できる。物質濃度として酸化ヘモグロビン濃度と還元ヘモグロビン濃度を求めることにより、酸素飽和度分布を算出できる。また、物質濃度としては、グルコース濃度、コラーゲン濃度、メラニン濃度、脂肪や水の体積分率なども求められる。
【００１８】
また、超音波エコーイメージング装置における特性情報とは、被検体内の音響インピーダンス差、音響インピーダンス差のある位置、音速、密度などを示す情報である。
【００１９】
被検体内の各位置の特性情報に基づいて、二次元または三次元の特性情報分布が得られる。分布データは画像データとして生成され得る。特性情報は、数値データとしてではなく、被検体内の各位置の分布情報として求めてもよい。すなわち、初期音圧分布、エネルギー吸収密度分布、吸収係数分布や酸素飽和度分布などの分布情報である。あるいは超音波エコーイメージングの場合、音響インピーダンスに関する分布情報である。
【００２０】
本発明でいう音響波とは、典型的には超音波であり、音波、音響波と呼ばれる弾性波を含む。トランスデューサ等により音響波から変換された電気信号を音響信号とも呼ぶ。ただし、本明細書における超音波または音響波という記載は、それらの弾性波の波長を限定する意図ではない。光音響効果により発生した音響波は、光音響波または光超音波と呼ばれる。光音響波に由来する電気信号を光音響信号とも呼ぶ。超音波エコーイメージングで発生した音響波は、超音波エコーまたはエコー波とも呼ばれる。超音波エコーに由来する電気信号を超音波信号やエコー信号とも呼ぶ。分布データは、光音響画像データ、超音波画像データ、再構成画像データなどとも呼ばれる。
【００２１】
本発明の被検体情報取得装置は、人や動物の血管疾患や悪性腫瘍などの診断や化学治療の経過観察に好適である。被検体の例として、被検者の乳房や手のような生体の一部、マウスなどヒト以外の動物、無生物、ファントムなどが挙げられる。
【００２２】
［第一の実施形態］
以下、図面を参照しながら、本発明の実施形態について詳細に説明する。なお、同一の構成要素には原則として同一の符号を付して説明を省略する。
【００２３】
＜被検体情報取得装置の構成＞
図１（ａ）は、本実施形態にかかる被検体情報取得装置の模式図である。以下、装置の各構成要素について説明する。装置は、探触子１１０、照射部１２０、走査部１３０、処理部１４０、表示部１５０を有する。測定対象は、被検体１００である。図１（ｂ）は、処理部１４０と周辺の構成との関係を示す模式図である。
【００２４】
処理部１４０は、バス２００を介して被検体情報取得装置の各構成要素の動作を制御している。また、処理部１４０は、後述する被検体情報取得方法が記述されたプログラムを保持しており、そのプログラムを読み出して被検体情報取得装置に被検体情報取得方法を実行させる。
【００２５】
照射部１２０は、光Ｌを被検体１００に照射する。すると、光音響効果により被検体内部や被検体表面から光音響波ＰＡが発生する。探触子１１０は、伝搬した音響波を受信して時系列の電気信号を取得し、受信信号とする。処理部１４０は受信信号に対する処理を行い、表示部１５０に表示させるための画像データを生成する。
【００２６】
＜構成要素の詳細な説明＞
以下、本実施形態に係る被検体情報取得装置の各構成の詳細を説明する。
【００２７】
（被検体１００）
被検体は本発明の被検体情報取得装置の一部を構成するものではないが、以下に説明する。本発明の被検体情報取得装置は、人や動物の悪性腫瘍や血管疾患などの診断や化学治療の経過観察などを主な目的とする。よって、被検体としては生体、具体的には人体や動物の四肢、指、乳房、頭部、頸部、腹部などの診断対象部位が想定される。
また、被検体内部にある光吸収物質は、被検体内部で相対的に光吸収係数が高いものとする。例えば、人体が測定対象であればオキシヘモグロビンあるいはデオキシヘモグロビンやそれらを含む多く含む血管あるいは新生血管を多く含む悪性腫瘍、血管の他にメラニン等の色素を含む正常皮膚や疾患を持つ皮膚等が測定対象たる光吸収物質となる。その他、頸動脈壁のプラークなども測定対象となる。
【００２８】
（探触子１１０）
探触子１１０は、音響波を検出可能な素子であるトランスデューサを含む。トランスデューサは、音響波を受信してアナログ信号である電気信号に変換する。探触子またはトランスデューサは、音響探触子、探触子、音響波探触子、音響波検出素子、音響波検出器、音響波受信器、などとも呼ばれる。トランスデューサ１１０は、圧電現象、光の共振、静電容量の変化を用いたものなど、音響波を受信できれば、どのようなものでもよい。本実施形態で利用される音響波は、典型的には数１００ＫＨｚから１００ＭＨｚの周波数成分で構成される。そのためトランスデューサとしては、これらの周波数を検出可能なものが好ましい。また探触子は、感度が高く、周波数帯域が広いことが望ましい。例えばＰＺＴ（チタン酸ジルコン酸鉛）などを用いた圧電素子、ＰＶＤＦ（ポリフッ化ビニリデン）などの高分子圧電膜材料、ＣＭＵＴ（容量性マイクロマシン超音波トランスデューサ）、ファブリペロー干渉計を用いたものなどが挙げられる。
【００２９】
探触子１１０が備えるトランスデューサの数は、単一でも複数でも良い。単一のトランスデューサの場合は、矩形、円形、平面、球面、楕円面等の形状をとることができる。単一素子のトランスデューサを用いる光音響顕微鏡も、本発明の適用対象である。
複数のトランスデューサの場合は、複数のトランスデューサが１Ｄ、１．５Ｄ、２Ｄ等に配置されたアレイトランスデューサを用いても良い。また、複数のトランスデューサを、各トランスデューサの指向軸が集まる高感度領域が形成されるように、お椀状や球冠状の支持体に配置しても良い。また、アレイトランスデューサの支持体又はお椀状の支持体と、光を照射する照射部とを一体化して、同時に移動可能としても良い。
【００３０】
（照射部１２０）
照射部１２０は、非図示の光源で発生させた光を被検体１００に照射する。照射部１２０は、処理部１４０による制御に従って動作しても良いし、処理部１４０と協調しながら、照射部自身が備える制御回路による制御に従って動作しても良い。前者の場合は処理部１４０が照射制御部を兼ねており、後者の場合は照射部が照射制御部を備えている。また、照射制御部を別途設けても良い。照射制御部は、ユーザにより指定された照射制御情報や、予めメモリに保存された照射制御情報を取得し、照射タイミング、照射光量、パルス光のパルス長やパルス間隔、照射光の波長などの照射条件を制御する。また、走査制御部の位置制御情報と連動することにより、照射位置を制御することもできる。照射部は、本実施形態での照射手段に相当する。
【００３１】
照射部１２０は、典型的にはレンズやミラーなどの光学部品を含む光学系から成る。照射部１２０は、光を所望の分布形状に成型して被検体１００に照射する。光学部品として、光ファイバなどの光を伝搬させる導波路、光を反射するミラー、光を集光したり拡大したり形状を変化させるレンズ、光を分散・屈折・反射するプリズム、光を拡散させる拡散板などを利用できる。光源から発せられた光を被検体に所望の形状で照射できるのであれば、どのような光学部品を用いても良い。なお、光源が発する光そのものを所望の光として被検体１００に照射できる場合、光源を照射部１２０とみなせる。
【００３２】
非図示の光源としては、数ナノから数マイクロ秒オーダーのパルス光を発生可能なパルス光源が好ましい。効率的に光音響波を発生させるためには、被検体の熱特性に応じて十分短い時間に光を照射させる必要がある。具体的には、光源が数百ナノ秒以下のパルス幅
の光を発生可能であることが好ましい。被検体が生体である場合は、光源から発生するパルス光のパルス幅は１０〜５０ナノ秒程度が好適である。
【００３３】
パルス光の波長は、被検体を構成する成分のうち特定の成分に吸収される特定の波長であって、被検体内部まで光が伝搬する波長であることが好ましい。被検体が生体の場合、好適な波長は５００ｎｍ以上、１２００ｎｍ以下であり、より好ましくは、７００ｎｍ以上、１１００ｎｍ以下である。ただし、比較的生体表面付近の生体組織の光学特性値分布を求める場合は、上記の波長領域よりも範囲の広い波長領域（例えば４００ｎｍ以上、１６００ｎｍ以下）も使用できる。光源としてはレーザ、フラッシュランプ、発光ダイオードを利用できる。レーザとして、固体レーザ、ガスレーザ、色素レーザ、半導体レーザなど様々なレーザを使用できる。例えば、アレキサンドライトレーザ、Ｙｔｔｒｉｕｍ−Ａｌｕｍｉｎｉｕｍ−Ｇａｒｎｅｔレーザ、Ｔｉｔａｎ−Ｓａｐｐｈｉｒｅレーザなどを利用できる。
【００３４】
ここでは光を発生する光源を挙げたが、電磁波を発生する手段にしてもよい。例えば、マイクロウェーブ源を用いても、光音響イメージングと同様の原理により被検体情報を取得できる。
【００３５】
（走査部１３０）
走査部１３０は、探触子１１０を走査して、被検体１００に対する相対的な位置を変化させる。また、走査部１３０は、照射部１２０を走査して、被検体１００に対する相対的な位置を変化させても良い。探触子１１０と照射部１２０を連動させて走査しても良いし、探触子１１０と照射部１２０を独立して走査しても良い。探触子１１０と照射部１２０の移動を連動させる場合、両者を結合して一体的に移動させても良い。走査部は、本実施形態での走査手段に相当する。
【００３６】
図１（ａ）における走査部１３０は、探触子１１０及び照射部１２０をｘｙ平面で一体的に走査する。走査部１３０は、処理部１４０による制御に従って動作しても良いし、処理部１４０と協調しながら、走査部自身が備える制御回路による制御に従って動作しても良い。前者の場合は処理部１４０が走査制御部を兼ねており、後者の場合は走査部が走査制御部を備える。また、走査制御部を別途設けても良い。走査制御部は、ユーザにより指定された走査制御情報や、予めメモリに保存された走査制御情報を取得し、走査の軌跡、走査開始と終了のタイミング、走査中の速度などの走査条件を制御する。また、照射制御部の照射制御情報と連動することにより、照射位置と照射タイミングを制御することもできる。
【００３７】
処理部１４０が走査制御を行う場合、処理部１４０は、電磁波照射と音響波の受信を行った位置（測定位置）を、エンコーダ等の位置情報取得手段を用いて座標値としてメモリに保存する。測定位置情報は、被検体の画像化において用いられる。本実施形態ではｘｙ平面での走査としたが、ｚ方向も含めた３次元的な走査であってもよい。
【００３８】
走査経路として、ラスタ軌道、螺旋軌道、円軌道等、様々な軌道を利用できる。音響波測定（すなわち、照射部１２０からの光照射と探触子１１０による音響波受信）は、これら軌道上の測定位置において行われる。音響波測定における「測定位置」や「走査点」という言葉は、走査軌道上において光を照射して音響波を受信する位置を指す。ここで、探触子の停止・音響波測定・探触子の再移動を繰り返す方式（ステップアンドリピート方式）の場合は、探触子が一時停止する位置を測定位置とする。また、探触子が連続的に移動する方式の場合は、光を照射する時の探触子位置を測定位置としても良いし、当該光照射により発生する光音響波を受信する期間における任意の位置（例えば、受信期間の中心時刻における探触子の位置）を測定位置としても良い。
【００３９】
また、ハンドヘルド型の探触子１１０を手動によって走査しても良い。この場合、探触子の位置に関する情報として、探触子の空間位置や姿勢を取得すると良い。探触子の空間位置や姿勢は、例えば、モーションキャプチャー用カメラや、磁気により位置を取得する装置などを用いて取得できる。
【００４０】
（処理部１４０）
処理部１４０は、受信信号を用いて被検体内部の被検体情報を取得するための演算を行う。典型的にはＣＰＵ、ＧＰＵなどの素子や、ＦＰＧＡ、ＡＳＩＣなどの回路から構成される。処理部１４０は、プログラムや、制御情報や、音響波測定の結果などを保存するメモリを備えることが好ましい。なお、処理部１４０は、１つの素子や回路から構成されるだけではなく、複数の素子や回路から構成されていてもよい。また、被検体情報取得方法で行われる各処理をいずれの素子や回路が実行してもよい。各処理を実行する装置を総称して本実施形態に係る処理部とする。処理部１４０として、典型的にはワークステーションやパーソナルコンピュータを利用できる。処理部１４０としてワークステーション等を用いる場合、当該ワークステーションのＵＩ（例えばキーボード、マウス、タッチパネル等）を用いてユーザからの指示情報の入力を受け付けても良い。処理部は、本実施形態での制御手段および情報取得手段に相当する。なお、図１（ｂ）に示すように、制御手段１４２および情報取得手段１４４を、処理部１４０を構成する機能ブロックとして実装しても良い。
【００４１】
処理部１４０は、Ａ／Ｄ変換器や、信号増幅器を含んでもよい。Ａ／Ｄ変換器は、探触子１１０により音響波から変換されたアナログの電気信号をデジタル信号に変換する。信号増幅器は受信信号に増幅処理を施す。また、Ａ／Ｄ変換器や信号増幅器を処理部１４０とは別の信号処理部として設けても良い。
【００４２】
また、処理部１４０は、同時に複数の信号をパイプライン処理できるように構成されていることが好ましい。これにより、被検体情報を取得するまでの時間を短縮することができる。また、処理部１４０は非一時的な記録媒体を有し、被検体情報取得方法で行われるそれぞれの処理を、自身が実行するプログラムとして保存しておくことができる。
また、処理部１４０は、探触子１１０と共通の筺体に収められた構成で提供されてもよい。ただし、筺体に収められた処理部で一部の信号処理を行い、残りの信号処理を筺体の外部に設けられた処理部で行ってもよい。この場合、筺体の内部および外部に設けられた処理部を総称して、本実施形態に係る処理部とすることができる。その他にも、図示された被検体情報取得装置の各構成要素の配置は一例であり、総体として本発明に必要な処理を実行できるのであれば、どのような配置でも良い。
【００４３】
（表示部１５０）
表示部１５０は、処理部１４０から出力される特性情報である被検体情報を表示する装置である。表示部１５０として、例えば、液晶ディスプレイ、プラズマディスプレイ、有機ＥＬディスプレイ、ＦＥＤなどを利用できる。なお、表示部１５０または処理部１４０は、被検体情報を表示する際に、輝度値の調整等の画像処理を行っても良い。また、処理部１４０は、被検体情報の他に、操作者や被検者への指示やメッセージを表示部１５０に表示させても良い。表示部は、本実施形態での表示手段に相当する。
【００４４】
＜被検体情報取得方法＞
次に、本実施形態に係る被検体情報取得方法の各工程を、図面を参照して説明する。なお、各工程は、処理部１４０が被検体情報取得装置の各構成の動作を制御することにより実行される。
【００４５】
図２Ａは、全体的な処理を示すフロー図である。
（ステップＳ１０）処理部１４０は、被検体情報取得の内容に関する情報を取得する。この情報には例えば、被検体の種類や被検体における所定の関心領域の広さや深さ、被検体情報の種類、被検体情報の所望の精度など、音響波測定に関する種々の情報が含まれる。処理部１４０は、ユーザにより入力された情報を取得したり、予めメモリに保存されている情報を読み出したりして音響波測定の内容に関する情報を取得する。
【００４６】
（ステップＳ２０）処理部１４０は、後述する第一の音響波測定および第二の音響波測定に関する情報を取得する。音響波測定に関する情報には、少なくとも走査軌道と、軌道上の音響波測定位置が含まれる。音響波測定位置は光照射位置および音響波受信位置を含む。音響波測定に関する情報も、ユーザ入力情報やメモリに保存された情報など任意の方法で取得できる。なお、ステップＳ１０とＳ２０とを一つのステップとしてもよい。
【００４７】
（ステップＳ３０）処理部１４０は、被検体情報取得の内容に関する情報と、音響波測定に関する情報に基づいて、装置の制御情報を設定する。制御情報には少なくとも、照射制御情報と走査制御情報が含まれる。処理部１４０は、第一および第二の音響波測定における軌道や測定位置に基づいて、測定開始後の各タイミングにおける、走査部１３０による移動方向や移動距離、光の照射タイミング、音響波受信タイミングを算出する。そして処理部１４０は、算出された情報に基づいて走査部１３０に関する走査制御情報、照射部１２０に関する照射制御情報、などのパラメータを算出して設定する。また、探触子のトランスデューサによる音響波の受信に関する受信制御情報のパラメータを算出して設定する。なお、予め保存済みの制御情報を読み出して設定するようにしても良い。その際、被検体の種類やユーザの指定に応じて複数の既定のパラメータから、今回の測定で利用するものを選択してもよい。
【００４８】
（ステップＳ４０）操作者は、被検体を所定の位置に設置する。
（ステップＳ５０）操作者は、制御情報の設定や被検体の設置が完了したことを確認した後、図２Ｂで述べる音響波測定を開始する。
【００４９】
図２Ｂは、本実施形態における被検体情報取得方法を示すフロー図である。ステップＳ１１０〜Ｓ１３０は第一の音響波測定を構成し、ステップＳ１４０〜Ｓ１６０は第二の音響波測定を構成する。なお、第一の音響波測定のうち、特に探触子の移動について検討する際に、第一の走査という名称を用いる。同様に、第二の音響波測定のうち、特に探触子の移動について検討する際に、第二の走査という名称を用いる。また、第一の音響波測定により得られた音響波や音響信号は、第一の信号と呼ぶことができる。同様に、第二の音響波測定により得られた音響波や音響信号は、第二の信号と呼ぶことができる。
【００５０】
（ステップＳ１１０：被検体内に光を照射し光音響波を発生させる工程）
照射部１２０が、被検体１００に光を照射する。すると、被検体内部や被検体表面の光吸収物質が光エネルギーを吸収して光音響波を発生させる。
（ステップＳ１２０：光音響波を受信して受信信号を取得する工程）
探触子１１０のトランスデューサが光音響波を受信（検出）して、処理部１４０に受信信号を出力する。
【００５１】
（ステップＳ１３０：第一の音響波測定の終了を判定する工程）
本工程では、処理部１４０が、第一の音響波測定を終了させるかどうかを判定する。具体的には処理部１４０は、第一の音響波測定で受信信号を取得すべき測定位置（図３（ａ）の白丸）全てで測定を終えるまで、探触子１１０および／または照射部１２０を走査してＳ１１０とＳ１２０を繰り返す。一方、未測定の測定位置が残っていない場合、第一の音響波測定を終了させて第二の音響波測定に移行する。
【００５２】
（ステップＳ１４０：被検体内に光を照射し光音響波を発生させる工程）
（ステップＳ１５０：光音響波を受信して受信信号を取得する工程）
第二の音響波測定においても、これらの処理は、ステップＳ１１０およびＳ１２０の処理と同様に実施される。
【００５３】
（ステップＳ１６０：第二の音響波測定の終了を判定する工程）
本工程では、処理部１４０が、第二の音響波測定を終了させるかどうかを判定する。具体的には処理部１４０は、第二の音響波測定で受信信号を取得すべき測定位置（図３（ｂ）の黒丸）全てで測定を終えるまで、探触子１１０および／または照射部１２０を走査してＳ１４０とＳ１５０を繰り返す。一方、未測定の測定位置が残っていない場合、第二の音響波測定を終了させる。
【００５４】
ここで、第一の音響波測定と第二の音響波測定について図３を用いて説明する。図３の破線は、ラスタ軌道で走査する場合の走査軌道を示す。図３（ｂ）は第二の音響波測定の様子を表し、黒丸は各測定位置を示す。処理部１４０は、第二の音響波測定の各測定位置で取得した受信信号を用いて被検体情報を生成する。
図３（ａ）は第一の音響波測定の様子を表し、白丸は各測定位置を示す。第一の音響波測定における測定位置の数は、第二の音響波測定における測定位置の数より少ない。そのため、第一の音響波測定は第二の音響波測定より短時間で完了する。これにより、第一の音響波測定では長周期の体動の影響を抑えた受信信号を取得することができる。図３の例では、第一の走査は、ラスタ走査である第二の走査の主走査方向と交差する方向に、直線状の経路を辿って行われる。
【００５５】
なお、長周期の体動とは、被検者の突発的で短期的な体動とは異なり、例えば被検者が長時間姿勢を維持することが難しいために少しずつ姿勢を変化させてしまうことを指す。例えば、被検部位が被検者の足であり、測定のために足を上げる姿勢を取らなければならない場合、測定が長期間に渡ると、被検者の疲労によって足の位置が徐々に変位（低下）してしまう。ただし、長周期の体動と短周期の体動の区別は相対的なものであり、時間や測定内容で決定されるものではない。音響波測定に要する時間全体のうち、被検体情報の精度に影響を与える程度の長期間に渡る体動であれば、本発明による補正の対象となる。
【００５６】
（ステップＳ１７０：第一の音響波測定における表面形状情報を取得する工程）
本工程では、Ｓ１１０からＳ１３０で取得した受信信号を用いて被検体表面位置（ｚ座標）を取得して、第一の表面形状情報とする。表面形状情報は、あるタイミングにおける被検体のｚ方向の変位を示す情報である。本実施形態では、第二の音響波測定と並行してＳ１７０の処理を行うので、総処理時間を短縮できる。ただし、第一の音響波測定の完了後であればいつ行ってもよく、第二の音響波測定の前や後に、Ｓ１７０の処理を行ってもよい。ここで、被検体表面という用語は、必ずしも厳密に被検体の最表面を意図するものではない。被検体の変位を検出するために有効な被検体情報が取得できる部位であれば、基準部位として利用可能である。かかる基準部位の情報を、Ｓ１７０やＳ１８０における情報取得の対象として良い。例えば、後述の図１５に示すように、被検体の最表面よりも深部にメラニン層が存在する。かかるメラニン層の被検体情報を取得して、体動量を算出するための基礎としても良い。Ｓ１７０やＳ１８０で取得する被検体表面位置のことを、基準部位位置と呼んでも良い。第一の音響波測定における表面形状情報は、第一の被検体形状情報と呼ぶことができる。
【００５７】
Ｓ１７０における被検体表面位置の取得には、例えば非特許文献１に記載の方法を利用できる。すなわち、処理部１４０は、時系列の受信信号から、被検体表面付近のメラニン層に由来する信号成分を検出する。そして、被検体表面に由来する信号が検出された時間
に音速をかけることで、被検体表面からトランスデューサまでの距離を取得する。この処理を被検体全体について行うことで、表面形状情報を取得する。得られた第一の表面形状情報は、測定位置ごとに測定のタイミングが異なる情報である。
上記の方法であれば、装置に新たな構成要素を設けずとも表面形状情報を取得できる。ただし、被検体表面の形状情報が得られるのであれば、方法はこれに限られない。例えば、撮像装置により光学像を取得して表面形状を取得してもよい。また、トランスデューサから被検体に超音波を送信し、エコー波がトランスデューサに戻るまでの時間に基づいて表面形状を取得しても良い。
【００５８】
（ステップＳ１８０：第二の音響波測定における表面形状情報を取得する工程）
本工程では、Ｓ１４０からＳ１６０で取得した受信信号を用いて被検体表面位置（ｚ座標）を取得して、第二の表面形状情報とする。被検体表面位置の取得には、Ｓ１７０と同様の方法を利用できる。第二の音響波測定における表面形状情報は、第二の被検体形状情報と呼ぶことができる。
【００５９】
処理部１４０は、第二の表面形状情報と第一の表面形状情報を組み合わせる。具体的には処理部１４０は、Ｓ１８０で得た第二の音響波測定における各測定位置の被検体表面位置を、第一の音響波測定が行われた時間ｔに対してプロットする。なお、第一の音響波測定は第二の音響波測定よりも短い時間で行われているため、図４のように時間を基準にしてプロットする場合は、第一の音響波測定の測定タイミングを第二の音響波測定の測定タイミングに適合させるような時間補正を行う。一方、測定位置または走査距離を基準にしてプロットを行う場合、第一の音響波測定と第二の音響波測定の測定位置または走査距離を合わせるようにする。
【００６０】
第一の音響波測定および第二の音響波測定における被検体表面位置をプロットした様子を図４（ａ）に示す。図４（ａ）の縦軸は被検体表面位置を示し、単位は例えば［ｍｍ］である。横軸は時間を示し、単位は例えば［ｓ］である。なお、横軸の時間は、第二の音響波測定の時間を基準としている。図４（ａ）の黒丸は、第二の音響波測定の各測定位置での被検体表面位置である。図４（ａ）の白丸は、第一の音響波測定の各測定位置での被検体表面位置を、測定位置が最も近い第二の音響波測定の測定位置の時間にプロットしたものである。なお、音響波測定の測定位置を設定する際には、第一の音響波測定の各測定位置が、第二の音響波測定の各測定位置のいずれかに重なるようにすることが好ましい。
【００６１】
（ステップＳ１９０：体動補正量を取得する工程）
本工程では、Ｓ１７０とＳ１８０で得た被検体表面位置を用いて、体動の補正量を取得する。図４（ｂ）の縦軸は、ｚ方向における被検体表面の変位量に対応する体動量Δｚを表し、単位は例えば［ｍｍ］である。体動量は、被検体の変位の量などを示す変位情報と考えられる。
【００６２】
Ｓ１９０において、まず処理部１４０は、第二の音響波測定の測定位置に対応する第一の音響波測定が存在する時間を抽出する。図４（ａ）からは、ｔ１〜ｔ５まで５つの測定タイミングが抽出される。なお、時間ではなく測定位置や走査距離を基準にしてプロットを行う場合にも同様に対応位置の抽出を行う。
続いて処理部１４０は、第二の音響波測定における被検体表面位置から第一の音響波測定における被検体表面位置を減算することで、ｚ方向の体動量Δｚを取得する。これを図４（ｂ）に示す。ここで得られた体動量は、被検体の変位を示す変位情報と呼ぶことができる。
以上の処理によって、従来は想定されていなかった、長周期に渡る被検体の変位を示す情報を取得できるという効果が得られる。
【００６３】
続いて処理部１４０は、図４（ｂ）にプロットされた体動量を、第二の音響波測定における各測定位置の時間に対して補間処理を行う。補間後の様子を図４（ｃ）に示す。図４（ｃ）のように、第二の音響波測定における全ての測定位置での補間値を求めることが好ましい。
続いて処理部１４０は、図４（ｃ）の符号を反転して図４（ｄ）に示す体動補正量Δｚ’を得る。
【００６４】
以上の処理によって、長期的な被検体の変位を補正して変位の影響を低減するための補正情報を取得できるという効果が得られる。
【００６５】
（ステップＳ２００：被検体情報を取得する工程）
本工程では、Ｓ１９０で取得した体動補正量を用いて被検体１００内の初期音圧分布ｐ０（ｒ）を算出する。
ここで、ｒは画像再構成する位置の位置ベクトル、ｒ０は測定位置の位置ベクトルを表す。ｄΩ０はトランスデューサが位置ｒを見込む立体角、Ω０は全測定位置の立体角の合計を表す。
【数１】
2020039809.tif 000003
は、時間ｔに関する時系列受信信号を表す。また、
【数２】
2020039809.tif 000004
は、時間ｔに音速をかけた距離換算時間を表す。
【００６６】
このとき、非特許文献２によれば、式（１）の逆投影法により受信信号を用いて被検体１００内の初期音圧分布ｐ０（ｒ）を算出できる。この処理を画像再構成と言う。
【数３】
2020039809.tif 000005
【００６７】
（補正の原理）
ここでＳ２００では、式（１）に体動補正量Δｚ’を含めた式（２）を用いる。
【数４】
2020039809.tif 000006
図５を用いて式（２）の説明をする。図５において、縦軸はｚ方向の位置を示す。横軸のうち上側は、被検体の様子を示す。簡便化のために、被検体はｚ方向には変位するが、ｘｙ方向には変位しないものとする。横軸のうち下側は、被検体に対する探触子の相対位置が走査によって変化する様子を示す。
【００６８】
図５（ａ）から図５（ｂ）にかけて、第二の音響波測定中に被検体１００が体動を起こして、体動量Δｚの変位が発生した場合を示す。この間、探触子は、ｉ番目の測定位置からｉ＋１番目の測定位置に移動している。また本実施形態では、図５（ｂ）に示すように、被検体内部の音源（光吸収物質）の変位量は、Ｓ１９０で得た体動量Δｚと同じとしている。この状況下で画像再構成を行うと、図５（ａ）と図５（ｂ）とで音源位置が異なることに起因して、ｉ番測定位置の受信信号に由来する音源像と、ｉ＋１番測定位置の受信信号に由来する音源像が、異なる位置に逆投影される。その結果、解像度やコントラストの低下、音源像の変形を引き起こす。
【００６９】
このような画像劣化を補正する様子を図５（ｃ）に示す。処理部１４０は、画像再構成上の探触子１１０の位置を体動量と反対方向に、体動補正量Δｚ’だけずらす。これにより、図５（ｃ）にあるように、ｉ番測定位置とｉ＋１番測定位置とで見かけの音源位置が一致する。その結果、体動が補正される。本処理により、被検体の変位を補正した被検体情報を取得するという効果が得られる。
【００７０】
（表示例）
Ｓ２００の後、表示部１５０は、式（２）により体動が補正された初期音圧分布ｐ０（ｒ）を表示してもよい。表示の例を図６に示す。図６（ａ）は、体動補正処理のない場合のウィンドウ１５１ａと、補正がある場合のウィンドウ１５１ｂとを共に画面に表示する例である。補正前はコントラスト低下が発生しているが、補正後は体動の影響が低減された高コントラストな画像が表示される。操作者は、並列に表示された双方の画面を比較して補正処理の効果や精度を確認できる。
【００７１】
図６（ｂ）では、検出した体動量をウィンドウ１５１ｃに表示することで効果や精度の確認をし易くしている。ここでは横軸に時間を、縦軸に体動量を取っているが、横軸には測定位置や走査距離を取ってもよい。
【００７２】
図６（ｃ）と図６（ｄ）は、補正処理のありなしを切り替え可能にするボタン１５１ｄを有する例である。操作者はＧＵＩとして表示されるボタンをクリックすることで、図６（ｃ）のような補正前の状態と、図６（ｄ）のような補正後の状態とを切り替えることができる。この場合、補正前と補正後の像が画面上の同じ位置に表示されるので、比較が容易になる。例えば、補正処理はバックグラウンドで実行されていてボタンは画像を切り替えるだけとすることで、素早く表示を切り替えることができる。
【００７３】
図６（ｅ）は、補正のパラメータを入力するスライドバー１５１ｅがある場合を示す。例えば、平滑化のｘ方向、ｙ方向それぞれのカットオフ周波数を設定することができる。その他、補正に関するいかなるパラメータを追加しても良い。表示の切り替えを行うためのボタンやスライドバーの調整は、コンピュータのＵＩを操作者用の入力手段として利用することで実行できる。
【００７４】
また、測定中に大きな体動が生じた場合、表示部１５０に注意表示をしても良い。図７を用いて、第一の音響波測定を終えた後、第二の音響波測定の最中に許容閾値を超える体動が発生した場合について説明する。
【００７５】
図７（ａ）〜図７（ｃ）はそれぞれ、第二の音響波測定途中における図４（ａ）〜図４（ｃ）に対応する。図７の例では、第二の音響波測定中、測定と同タイミングで体動量Δｚを計算しながら、予め設定された体動量の許容閾値Δｚｔｈを超えたかどうかを判定する。体動量が許容閾値以下であれば走査を継続し、体動量が許容閾値を超えた場合、図７（ｄ）に示すように、表示部１５０に注意表示８１０を表示する。この時、キャンセルボタン８２０を備えておいて、操作者が測定を中止できることが好ましい。また、体動が許
容閾値を超えた場合、装置側で測定を中止するようにしても良い。許容閾値は、シミュレーションで体動補正処理の補正限界を求めておく等により設定することができる。
【００７６】
（走査軌道及び測定位置の例）
図８は、第一の音響波測定と第二の音響波測定の様々な態様の例を示す。図８（ａ）〜図８（ｅ）は、第二の音響波測定の際に図３（ｂ）と同様の走査軌跡及び測定位置を採用する場合の、第一の音響波測定における測定位置を示す。なお、ラスタ走査に限らず、主走査方向への主走査と、主走査方向に交差する副走査方向への副走査を含む走査方式であれば構わない。図８（ａ）の場合、第二の音響波測定における測定位置を４つおきに抽出した測定位置群を、第一の音響波測定における測定位置とする。この場合、第一の走査を直線状にできる点も好ましい。図８（ｂ）の場合、第二の音響波測定における測定位置を５つおきに抽出した測定位置群を、第一の音響波測定における測定位置とする。
【００７７】
図８（ａ）、図８（ｂ）の例では、第二の音響波測定の測定位置の数に関して、第一の音響波測定の測定位置が等間隔に配置されている。また、仮に第二の音響波測定の測定位置が走査距離に関して等間隔に配置されているとすると、第一の音響波測定の測定位置もまた、走査距離に関して等間隔に配置される。また、仮に第二の音響波測定の測定位置が走査時間に関して等間隔に配置されているとすると、第一の音響波測定の測定位置もまた、走査時間に関して等間隔に配置される。ただし、時間に関して等間隔な配置を採用する場合、ラスタ走査において主走査と副走査の切り替えなどのために走査速度が一定にならない場合もあることを考慮することが好ましい。
【００７８】
図８（ｃ）、図８（ｄ）、図８（ｅ）は、第一の音響波測定における測定位置を直線上に配置した例である。これらの場合、第一の走査の経路を直線状に設定できる。
【００７９】
図８（ｆ）は、第二の走査の軌道が任意の曲線であるような場合に、曲線と直線の交点を、第一の音響波測定での測定位置とする例である。
図８（ｇ）、図８（ｈ）は、第二の走査が螺旋軌道の場合の例である。図８（ｇ）では、螺旋軌道と直線の交点に、第一の音響波測定の測定位置を配置する。図８（ｈ）では、螺旋軌道上の第二の音響波測定における測定位置（黒丸）から３つおきに抽出した測定位置群を、第一の音響波測定における測定位置（白丸）とする。なお図８（ｈ）では、黒丸と白丸が重なる位置は、白丸を優先して表示している。
【００８０】
図９のいずれの場合でも、第一の音響波測定における測定位置と、第二の音響波測定における測定位置は同じ位置であることが好ましい。第一の音響波測定と第二の音響波測定の間で測定位置を一致させることにより、図４（ｂ）における体動量取得の精度が向上する。なお、ここで言う位置の一致とは、体動量算出や被検体情報算出における所望の精度との関係で一致していると見なせる程度の一致度合いであれば良い。
【００８１】
なお、第一の音響波測定と第二の音響波測定の間で測定位置が一致しない場合は、最近傍の測定位置同士を対応させたり、空間的な補間処理を行ったりすると良い。Ｓ１９０の体動量の補間においては、スプライン補間や多項式補間等を用いることができる。さらに、図８（ａ）、図８（ｂ）、図８（ｈ）のように、測定位置が走査時間（または走査距離）に対して等間隔となる場合は、標本化定理に基づくｓｉｎｃ補間やＬａｎｃｚｏｓ補間といった高精度な補間処理を使用できる。また、等間隔でない場合であっても、不等間隔標本化定理に基づく高精度な補間を実施することができる。補間の精度を高めるためには、測定位置間の間隔が体動の周期の半分より短くなるような第二の音響波測定における測定位置を、第一の音響波測定における測定位置とするのが好ましい。
【００８２】
本発明による長周期の体動の補正方法を、非特許文献１に記載されている短周期の体動
の補正方法と組み合わせても良い。非特許文献１では、抽出した被検体表面形状を平滑化し、元の被検体表面形状と平滑化形状の差分をとることで短周期の体動量を取得している。この方法を本発明の長周期体動補正と組み合わせる場合、平滑化に用いるフィルタ特性が排他的になることが好ましい。非特許文献１の短周期体動量補正は式（３）で表される。Δzsmallは短周期の体動量、ｚ（ｘ，ｙ）は被検体表面ｚ座標、ｚ’（ｘ，ｙ）は平滑化被検体表面ｚ座標、ｈ（ｘ，ｙ）は平滑化フィルタ、大文字は周波数領域であることを表す。Ｈの最大値が１に規格化してある場合、式（３）は平滑化フィルタを通過した成分が短周期の体動量になることを示す。
【数５】
2020039809.tif 000007
【００８３】
よって、長周期体動補正では、処理部１４０は、Ｓ１７０とＳ１８０において抽出したｚ（ｘ，ｙ）にフィルタｈ（ｘ，ｙ）をかけたｚ’（ｘ，ｙ）に対して、Ｓ１９０以降の処理を実施することで、短周期体動補正と排他的な成分を補正できる。例えば図３（ａ）においては、ｈ（ｘ，ｙ）のｙ方向のみのｈｙ（ｙ）を平滑化に用いることができる。以上により、長周期体動補正と短周期体動補正とで二重に補正をかけることを抑制でき、高精度な体動補正を実行することができる。
【００８４】
以上、本実施形態によれば、第一の音響波測定により短期間の間に被検体表面形状を取得することで体動を補正することができ、精度の高い被検体情報を生成することができる。
【００８５】
［第二の実施形態］
本実施形態では、被検体表面形状を空間的に補間して体動補正量を取得する。なお、第一の実施形態と同一の構成要素には、原則として同一の符号を付して説明を省略する。
【００８６】
＜被検体情報取得装置の構成＞
本実施形態にかかる被検体情報取得装置の模式図および処理部の周辺構成は、上述の第一の実施形態と同じである。本実施形態では、図１（ａ）における処理部１４０の実行する処理の内容が異なる。
【００８７】
＜被検体情報取得方法＞
次に、本実施形態に係る被検体情報取得方法の各工程を、図９を参照して説明する。なお、各工程は、処理部１４０が被検体情報取得装置の各構成の動作を制御することにより実行される。Ｓ１０１０、Ｓ１０２０、Ｓ１０４０からＳ１０８０、Ｓ１１００は、それぞれ第一の実施形態のＳ１１０、Ｓ１２０，Ｓ１４０からＳ１０８０、Ｓ２００と同一であるため、説明を省略する。
【００８８】
（Ｓ１０３０：第一の音響波測定の終了を判定する工程）
被検体情報取得装置は、図１０（ａ）の全ての測定位置で測定を終えるまで、探触子１１０および、または照射部１２０を走査して、ステップＳ１０１０とＳ１０２０の処理を繰り返す。そしてステップＳ１０３０において、第一の音響波測定の終了を判定する。
【００８９】
本実施形態の第一の音響波測定における測定位置を図１０（ａ）に白丸で示す。本実施形態の第二の音響波測定における測定位置を図１０（ｂ）に黒丸で示す。第一の音響波測定で得た被検体表面位置を空間的に補間するためには、図１０（ａ）に示す第一の音響波測定における測定位置群が、第二の音響波測定における走査領域（図１０（ｂ）に示す黒丸で囲まれる領域）を包含するようにすることが望ましい。例えば第一の走査と第二の走査が同じ軌道を描く場合、第一の音響波測定の開始測定位置と終了測定値が、それぞれ、第二の音響波測定の開始測定位置と終了測定位置と重なるようにすると良い。その場合、第一の音響波測定時の測定領域と第二の音響波測定時の測定領域は、互いに互いを包含していると言える。
また包含関係にならない場合でも、少なくとも、第一の音響波測定における測定位置群により囲まれる測定領域が、第二の音響波測定における測定位置群により囲まれる測定領域と重なる範囲が、可及的に広くなるようにすることが好ましい。
【００９０】
第一の実施形態同様、第一の音響波測定における測定位置の数は、第二の音響波測定における測定位置の数より少ない。そのため、第一の音響波測定は第二の音響波測定より短時間で完了するため、第一の音響波測定において取得した受信信号に基づいて、長周期の体動の影響を抑えた被検体表面情報（被検体形状情報）を取得できるようになる。
【００９１】
（Ｓ１０９０：体動補正量を取得する工程）
本工程では、ステップＳ１０７０とＳ１０８０で得た被検体表面位置を用いて、体動の補正量を取得する。図１１（ａ）は、第一の音響波測定により取得した被検体表面位置である。白丸で示された点は、図１０（ａ）の各測定位置に対応する。図１１（ｂ）は、第二の音響波測定により取得した被検体表面位置である。黒丸で示された点は、図１０（ｂ）の各測定位置に対応する。図を比較して分かるように、白丸で示された点は、黒丸で示された点よりも不足している。
【００９２】
そこで処理部１４０は、図１１（ａ）の各測定位置におけるデータを用いて、不足している測定位置の補間を行う。これにより、図１１（ｃ）に示されるように、被検体表面位置が補間された第一の表面形状情報が得られる。続いて処理部１４０は、図１１（ｂ）の第二の表面形状情報から図１１（ｃ）の第一の表面形状情報を減算し差分をとる。これにより、図１１（ｄ）に示す体動量Δｚを得る。本実施形態の体動量は、体動量分布情報として取得される。続いて処理部１４０は、図１１（ｄ）の符号を反転して図１１の（ｅ）に示す体動補正量Δｚ’を算出する。
【００９３】
以上、本実施形態によれば、第一および第二の音響波測定により取得された被検体表面形状に基づいて、体動を補正するための情報を得ることができる。本実施形態の体動補正量は被検体の形状に沿った分布情報として取得されるため、補正の精度が向上し、精度の高い被検体情報を生成可能となる。
【００９４】
［第三の実施形態］
本実施形態では、光源の繰り返し周波数（Ｐｕｌｓｅ  Ｒｅｐｅｔｉｔｉｏｎ  Ｒａｔｅ、ＰＲＲ）と放出する光パルスのエネルギー（Ｐｕｌｓｅ  Ｅｎｅｒｇｙ、ＰＥ）の関係を利用して高速な第一の音響波測定を行い、体動補正量を取得する。なお、第一の実施形態および第二の実施形態と同一の構成要素には、原則として同一の符号を付して説明を
省略する。
【００９５】
＜被検体情報取得装置の構成＞
図１２（ａ）は、本実施形態にかかる被検体情報取得装置の模式図である。図１２（ｂ）は、処理部１４０と周辺の構成との関係を示す模式図である。本実施形態では、処理部１４０の処理内容の一部が第一の実施形態とは異なる。本実施形態は光源の制御に特徴があるため、図中に光源１８０を明示している。光源１８０は、処理部１４０の制御に応じて照射部１２０から光を出射させる。ただし本実施形態に必要な光制御を実施できるのであれば、外部の光源から光を導く構成でも構わない。
【００９６】
（光源１８０）
光源１８０は、第一の実施形態で説明した光源の特徴に加えて、ＰＲＲを変更することができる。本実施形態では、第一の音響波測定と第二の音響波測定とでＰＲＲを変更する。
【００９７】
＜被検体情報取得方法＞
次に、本実施形態に係る被検体情報取得方法の各工程を、図１３を参照して説明する。なお、各工程は処理部１４０が被検体情報取得装置の各構成の動作を制御することにより実行される。Ｓ１００２０、Ｓ１００３０、Ｓ１００６０からＳ１０１２０はそれぞれ、第一の実施形態におけるＳ１１０、Ｓ１２０、Ｓ１４０からＳ２００と同一であるため、説明を省略する。
【００９８】
（Ｓ１００１０：光源のＰＲＲを第一の音響波測定で用いる値に設定する工程）
本工程では、光源のＰＲＲをＳ１００５０で設定される値（第二の音響波測定の走査で用いる値）より大きな値に設定する。すなわち、第一の音響波測定の繰り返し周波数は、第二の音響波測定の繰り返し周波数よりも高い。その結果、Ｓ１００２０におけるパルス光照射とＳ１００３０における音響波受信は、比較的早い周期で実行される。
【００９９】
（Ｓ１００４０：第一の音響波測定の終了を判定する工程）
本工程では、第一の音響波測定の終了を判定する。第一の音響波測定で受信信号を取得すべき測定位置（図１４（ａ）の白丸）全てで測定を終えるまで、探触子１１０および、または照射部１２０を走査してＳ１００２０とＳ１００３０を繰り返す。本実施形態では、第二の音響波測定より第一の音響波測定のＰＲＲが大きい。そのため、図１４のように第一の音響波測定における測定位置群と第二の音響波測定における測定位置群が同じであっても、第一の音響波測定を第二の音響波測定より短時間で完了できる。そのため、第一の音響波測定では長周期の体動の影響を抑えた受信信号を取得できる。さらに、測定位置数を第一の実施形態や第二の実施形態より増やすことができるため、体動補正量の精度を高めることができる。なお、必ずしも第一の音響波測定と第二の音響波測定との間で測定位置の数は同じでなくても良い。上記実施形態と同様に、第一の音響波測定時の測定位置数を減らして、第一の音響波測定にかかる時間を更に短縮しても良い。
【０１００】
（Ｓ１００５０：光源のＰＲＲを第二の音響波測定で用いる値に設定する工程）
本工程では、光源のＰＲＲをＳ１００１０で設定される値（第一の音響波測定の走査で用いる値）より小さな値に設定する。すなわち、第二の音響波測定の繰り返し周波数は、第一の音響波測定の繰り返し周波数よりも低い。その結果、Ｓ１００６０におけるパルス光照射とＳ１００７０における音響波受信は、比較的遅い周期で実行される。
【０１０１】
ここで本実施形態でのＰＲＲおよびＰＥの設定の特徴を説明する。図１５（ａ）は、被検体断面の模式図である。被検体表面は表皮の角質層などで構成され、その下部（深部）にはメラニン色素を含むメラニン層が存在する。更にその下部には真皮は皮下組織が存在
し、血管を含む層も存在する。図１５（ｂ）は、第一の音響波測定において被検体内部に到達する光量を示す到達光量模式図である。図１５（ｃ）は、第二の音響波測定において被検体内部に到達する光量を示す到達光量模式図である。図１５（ｂ）および図１５（ｃ）においては、ハッチングが濃いほど光量が大きい。
【０１０２】
ここで、第一の音響波測定では被検体表面の形状情報さえ取得できれば良い。そのため図１５（ｂ）では、ＰＥの設定値を、被検体表面近傍のメラニン層まで光が到達する程度に小さくしている。一方、第二の音響波測定では、診断対象となることが多い血管等の層（被検体深部）まで光を到達させる必要がある。そのため図１５（ｃ）では、ＰＥの設定値を図１５（ｂ）の場合よりも大きくしている。さらに、このように第一の音響波測定におけるＰＥを第二の音響波測定よりも小さくすることにより、第一の音響波測定におけるＰＲＲの設定値を大きくしても、被検体の安全のために定められる規格（ＭＰＥ等）を満足させることが可能になる。また図１６に示すように、光源１８０は一般に、ＰＲＲとＰＥが単調減少の傾向にある。よって、第一の音響波測定では第二の音響波測定よりＰＲＲを高く設定できるため、走査にかかる時間を短くすることができる。
【０１０３】
本実施形態のステップＳ１０１１０で体動補正量を算出する際には、第一の実施形態で説明した時間的または距離的な補間を行っても良いし、第二の実施形態で説明した空間的な補間を行っても良い。
【０１０４】
以上、本実施形態によれば、第一の音響波測定の内容に応じた光源ＰＲＲを設定することで、より高精度な被検体表面形状情報を取得できる。その結果、体動補正情報の算出精度も向上するため、精度の高い被検体情報を生成できる。
【０１０５】
［第四の実施形態］
＜被検体情報取得装置の別の構成例＞
図１７は、本発明を適用可能な被検体情報取得装置の別の構成例を示す模式図である。図１７の被検体情報取得装置は、探触子ユニット１７０１、探触子ユニット保持機構１７１３、信号取得部１７１９、光源１７２０、装置制御部１７２２、表示装置１７２１を有する。
【０１０６】
探触子ユニット１７０１は、被検体に対して光を照射し、被検体から発生した音響波を受信するユニットである。探触子ユニット１７０１は、被検体への光の照射を行う光照射部１７０３と、音響波の受信を行う音響探触子１７０２と、走査機構１７０４を有する。光照射部１７０３および音響探触子１７０２は、走査機構１７０４によって一体的に移動可能に構成される。探触子ユニット１７０１と被検体１７０９は、生体接触膜１７０６を介して接触する。以下の説明では、生体接触膜１７０６を「（探触子ユニットと被検体との）接触面」と称する。光照射部は、本実施形態の照射手段に相当する。走査機構は、本実施形態の走査手段に相当する。
【０１０７】
生体接触膜１７０６は、ポリエチレンテレフタラート等によって構成された膜である。生体接触膜１７０６は、被検体によって変形しない強度と、光と音響波を透過させる特性を有する材質であることが好ましい。本実施形態では、生体接触膜の開口部分は、３０ｍｍ×３０ｍｍである。生体接触膜１７０６と音響探触子１７０２との間には、音響マッチング剤（音響伝播媒質）である水１７０５が貯留されている。なお、生体接触膜１７０６は、膜の中での音響波の多重反射を避けるために、１００ミクロン程度の厚みとすることが好ましい。
【０１０８】
探触子ユニット保持機構１７１３は、探触子ユニット１７０１を保持および移動するための機構である。探触子ユニット保持機構１７１３は、Ｚ軸方向への移動を行うためのＺ
軸ステージ１７１１と、Ｘ軸方向への移動を行うためのＸ軸ステージ１７１６を含む。Ｚ軸ステージ１７１１は、Ｚ軸ハンドル１７１２によって移動可能に構成される。これにより、探触子ユニット１７０１を被検体１７０９に対してＺ軸方向に移動させられる。Ｚ軸ステージの位置はＺ軸エンコーダ１７１４によって検出される。これにより、探触子ユニットのＺ軸方向における位置を算出できる。また、Ｘ軸ステージ１７１６は、Ｘ軸ハンドル１７１７によって移動可能である。これにより、探触子ユニット１７０１を被検体１７０９に対してＸ軸方向に移動させられる。Ｘ軸ステージの位置はＸ軸エンコーダ１７１８によって検出される。これにより、探触子ユニットのＸ軸方向における位置を算出できる。
【０１０９】
光源１７２０は、被写体に照射するパルス光を発生させる装置である。光源１７２０として、上述の実施形態と同様の光源装置を利用できる。波長、パルス長などについても上述の実施形態と同様に設定できる。なお、光照射のタイミング、波形、強度等は、後述する装置制御部１７２２によって制御される。本構成例では、パルス幅を１０ナノ秒とし、繰り返し周波数を２００Ｈｚとする。また、５３２ｎｍと１０６４ｎｍの波長を切り替えることができるＹＡＧレーザを使用する。５３２ｎｍは、生体における吸収が大きい波長であるが、本実施形態における光音響装置は、被検体表面から５ｍｍ程度までを測定の対象とするため、当該波長が利用できる。なお、１０６４ｎｍの波長を用いることによって、血管とメラニンを識別することもできる。
【０１１０】
光源１７２０から出射した光は、光照射部１７０３である光ファイバを用いて被検体１７０９に照射される。なお、光ファイバは、音響探触子１７０２を中心にリング状に配置してもよい。また、光はレンズで集光させるより、ある程度の面積に広げる方が、生体への安全性ならびに診断領域を広げられるという観点で好ましい。
【０１１１】
音響探触子１７０２は、被検部の内部から到来する音響波を受信して、時系列の電気信号に変換する手段である。音響探触子１７０２として、上述の実施形態と同様の探触子またはトランスデューサを利用できる。受信する波長や受信方式についても上述の実施形態と同様に設定できる。
【０１１２】
本実施形態における音響探触子１７０２は、ＰＺＴと音響レンズからなる音響フォーカス型の探触子であり、所定の焦点から発生した音響波を効率よく受信することができる。直径は６ｍｍで、中心周波数は５０ＭＨｚである。探触子の先端には石英ガラスによる音響レンズが組みつけられ、その開口数は０．６である。ＸＹ平面内の解像度は、音響探触子１７０２の性能によって決まり、本実施形態では６０μｍ程度である。また、深さ方向の解像度は、検出できる波長の８割程度（３０μｍ程度）となる。焦点は、探触子から４ｍｍ離れた位置にあり、生体接触膜１７０６の位置と一致する。なお、焦点の位置は、探触子側に近づけて配置したほうがよい場合もあり、その場合は、例えば０．５ｍｍ近づける。
【０１１３】
信号取得部１７１９は、音響探触子１７０２が取得したアナログの電気信号を増幅してデジタル信号に変換する手段である。信号取得部１７１９は、受信信号を増幅する増幅器、アナログ信号をデジタル変換するＡ／Ｄ変換器を用いて構成されてもよい。また、信号取得部１７１９は、複数のプロセッサや演算回路から構成されていてもよい。
【０１１４】
本実施形態では、サンプリング周波数を５００ＭＨｚ、サンプリング数を８１９２とする。サンプリングは、光照射のタイミングを表すトリガ信号の発生から所定の時間経過後に開始される。なお、信号取得部１７１９は、受信信号を記憶するＦＩＦＯ等のメモリと、ＦＰＧＡチップ等の演算回路をさらに有していてもよい。また、装置制御部１７２２は、汎用コンピュータや、専用に設計されたワークステーションによって実現されてもよい
。
【０１１５】
装置制御部１７２２は、デジタル変換された信号（光音響信号）に基づいて、再構成処理を行うことで、被検体の内部の光吸収係数や酸素飽和度等といった被検体情報を取得する手段（本発明における画像生成手段）である。具体的には、収集された電気信号から三次元の被検体内における初期音圧分布を生成する。
また、装置制御部１７２２は、被検体に照射される光量に関する情報に基づいて、被検体内における三次元の光強度分布を生成する。三次元の光強度分布は、二次元の光強度分布に関する情報から光拡散方程式を解くことで取得できる。また、光音響信号から生成された被検体内の初期音圧分布と、三次元の光強度分布とを用いて、被検体内の吸収係数分布を得ることができる。また、複数の波長における吸収係数分布を演算することで、被検体内の酸素飽和度分布を得ることができる。装置制御部は、本実施形態の制御手段および情報取得手段に相当する。
【０１１６】
なお、装置制御部１７２２は、光量分布の計算や背景の光学係数取得に必要な情報処理、信号補正など、所望の処理を実施する機能を有していてもよい。また、装置制御部１７２２は、後述する表示装置や入力インタフェースを介して、測定パラメータの変更、測定の開始・終了、画像の処理方法の選択、患者情報や画像の保存、データの解析などに関する指示を取得してもよい。また、装置制御部１７２２は、光音響装置が有する各構成要素の制御を行う手段でもある。例えば、被検体に対する光照射の制御、音響波や光音響信号の受信制御、探触子ユニットの移動制御など、装置全体の制御に関する指令を行う。
【０１１７】
装置制御部１７２２は、ＣＰＵとＲＡＭ、不揮発メモリ、制御ポートを有するコンピュータで構成してもよい。不揮発メモリに格納されたプログラムがＣＰＵで実行されることにより制御が行われる。装置制御部１７２２は、汎用コンピュータや、専用に設計されたワークステーションによって実現されてもよい。また、装置制御部１２２の演算機能を担うユニットは、ＣＰＵやＧＰＵ等のプロセッサ、ＦＰＧＡチップ等の演算回路で構成されていてもよい。これらのユニットは、単一のプロセッサや演算回路から構成されるだけでなく、複数のプロセッサや演算回路から構成されていてもよい。
【０１１８】
また、装置制御部１７２２の記憶機能を担うユニットは、ＲＯＭ、磁気ディスクやフラッシュメモリなどの非一時記憶媒体や、ＲＡＭなどの揮発性の媒体であってもよい。なお、プログラムが格納される記憶媒体は、非一時記憶媒体である。なお、これらのユニットは、１つの記憶媒体から構成されるだけでなく、複数の記憶媒体から構成されていてもよい。装置制御部１７２２の制御機能を担うユニットは、ＣＰＵなどの演算素子で構成される。
【０１１９】
表示装置１７２１は、装置制御部１７２２が取得した情報およびその加工情報を表示する手段であり、典型的にはディスプレイ装置である。表示装置１７２１は、複数の装置であってもよいし、単一の装置に複数の表示部を備え、並列表示が可能な装置であってもよい。なお、表示装置１７２１には、高解像度でカラー表示が可能な３０型以上、コントラスト比１０００：１以上の装置を用いることが望ましい。表示装置は、本実施形態の表示手段に相当する。
【０１２０】
上記のように、本構成例の被検体情報取得装置は、探触子ユニットを保持及び移動させるための探触子ユニット保持機構とは別に、探触子ユニットの内部に配置され、音響探触子と光照射部とを移動させる走査機構を備える。このような構成においても、走査機構の移動制御に基づいて音響探触子および光照射部を走査して、測定時間が比較的短い第一の音響波測定と、測定時間が比較的長い第二の音響波測定を行うことができる。その結果、第二の音響波測定により得られた被検体情報を補正して、体動などの被検体の変位の影響
が低減された高精度な光音響画像データを取得することが可能となる。
【０１２１】
［第五の実施形態］
本発明は、光音響イメージング装置に限られず、被検体表面位置を取得できるいかなる装置にも適用することができる。例えば、超音波エコーイメージング装置や超音波トモグラフィー装置であれば、被検体表面と音響マッチング剤の音響インピーダンス差によって被検体表面位置を取得できる。
【０１２２】
例えば超音波エコーイメージング装置を例として説明する。超音波エコーイメージング装置は、第一の超音波エコー測定として比較的に短期間の測定を行い、第二の超音波エコー測定として比較的に長期間の測定を行う。そして、第一の超音波エコー測定の結果に基づいて体動補正情報を算出し、第二の超音波エコー測定により得られた音響インピーダンス分布を補正する。このとき、第一の超音波エコー測定を短時間化する方法としては、測定位置の数を減らす方法や、エコー送信の周波数を高くする方法などがある。
【０１２３】
また、上記各実施形態において図１のｚ方向の補正について説明してきたが、ｘ方向およびｙ方向を対象として補正を行っても良い。例えば、音響波の集束点を形成できる探触子において、集束点位置に対応する受信信号が第一の音響波測定における測定位置と、対応する第二の音響波測定における測定位置とで一致するようにすることで、ｘ方向、ｙ方向の体動量を取得できる。音響波の集束点を形成できる探触子として例えば、アレイトランスデューサ、集束型トランスデューサ等がある。
【０１２４】
以上、特定の実施形態を参照しながら、本発明について詳解してきた。しかしながら、本発明は上記特定の形態に限定されず、本発明の技術思想を逸脱しない範囲で実施形態の修正をすることができる。以上述べたように、本発明によれば、光音響イメージングや超音波エコーイメージングの画質を改善しつつ、少ない計算時間で画像を提供することができる。
【０１２５】
（その他の実施例）
本発明は、上述の実施例の１以上の機能を実現するプログラムを、ネットワーク又は記憶媒体を介してシステム又は装置に供給し、そのシステム又は装置のコンピュータにおける１つ以上のプロセッサがプログラムを読出し実行する処理でも実現可能である。また、１以上の機能を実現する回路（例えば、ＡＳＩＣ）によっても実現可能である。
【符号の説明】
【０１２６】
１１０：探触子、１２０：照射部、１３０：走査部、１４０：処理部
【図１】
2020039809.tif 000008
【図２Ａ】
2020039809.tif 000009
【図２Ｂ】
2020039809.tif 000010
【図３】
2020039809.tif 000011
【図４】
2020039809.tif 000012
【図５】
2020039809.tif 000013
【図６】
2020039809.tif 000014
【図７】
2020039809.tif 000015
【図８】
2020039809.tif 000016
【図９】
2020039809.tif 000017
【図１０】
2020039809.tif 000018
【図１１】
2020039809.tif 000019
【図１２】
2020039809.tif 000020
【図１３】
2020039809.tif 000021
【図１４】
2020039809.tif 000022
【図１５】
2020039809.tif 000023
【図１６】
2020039809.tif 000024
【図１７】
2020039809.tif 000025