JP4679507B2 - 音響光学撮像方法および音響光学撮像装置 - Google Patents

Description

本発明は、音響光学撮像方法および音響光学撮像装置に関する。

音響光学撮像では、撮像物によって散乱された信号波の音響成分で局所振動子のうなりを検出する。この信号波は、光学特性に関する情報を得ようとする撮像物のある点における音響周波数の振動によって、周波数がシフトされたものである。

単一画素の音響光学撮像では、周波数のシフトを伴わずに散乱された局所振動子としての信号波の成分と、情報を伝達する信号波の音響成分とのうなりの測定に大きなノイズを伴う。なぜならば、このうなりの測定は検出面のたった１点でおこなわれる一方、これら２つの成分は、この検出面においてランダムに変化するからである。より良い信号／ノイズ比を示す情報を得るには、このうなりの振幅の２乗の、時間に関する総和をとらねばならない。

この欠点を解決するために、非特許文献１に記載されているように、多画素検出装置を用いて、時間よりも、検出装置の画素に関する総和をとることが可能である。これにおいては、情報を伝達する信号波の音響成分と、変調側波帯（局所振動子）と、の間のうなりが、一般的に低い取得周波数を有する多画素検出装置で検出するに十分低い周波数であるように、レーザの出力を音響周波数に近い周波数に変調する必要がある。
「Ultrasonic tagging of photon paths in scattering media: parallel speckle modulation processing」，Optics Letters，Vol.24, No.3, 1st, February 1999, page 181.

しかしながら、局所振動子のウェイトは大抵あまりにも低いという重大な問題が残る。したがって、最適なノイズでヘテロダイン検出を行うには、ヘテロダインのゲインが低くなりすぎてしまう。

本発明はこれらの欠点を解消することを目的としている。

この目的のために、以下の工程から構成される撮像物の音響光学撮像方法が、本発明にしたがって提供される。
（ａ）周波数ｆ_Ｉの入射光波と周波数ｆ_Ｒの基準光波とを生成し、当該基準波は入射波とコヒーレントであるとともに、既知の位相差φ_ｉ（ｔ）を示し、
（ｂ）振動生成装置によって、第１の撮像物方向に沿って、音響周波数ｆ_Ａで、撮像物の領域を振動させ、
（ｃ）前記入射波を撮像物に照射して、散乱信号波を生成し、
（ｄ）散乱信号波の少なくとも一部を検出装置に照射し、
（ｅ）基準波を撮像物に透過させずに検出装置に照射し、検出装置の点ｒに、時間ｔで変化する干渉画像Ｉ（ｒ，ｔ）を生成し、
（ｆ）干渉画像Ｉ（ｒ，ｔ）からデジタル情報を抽出し、
（ｇ）デジタル情報が関連する撮像物の測定点の座標を得る。

このように、局所振動子として働く基準波は、撮像物を透過する必要はない。これにより、充分なレベルの局所振動子が得られ、よりよい信号／ノイズ比で撮像のための測定点、特に医療のための撮像の測定点に関する有用な情報を得ることができる。更に、この撮像方法は、例えば、医療用画像に関連した撮像した組織の安全基準にかなった低い音響又は光学出力であっても有用な信号を得られる。

本発明の好ましい形態において、以下の条件によるものである。

工程（ｆ）では、検出装置に照射された散乱信号波の一部の音響成分が検出され、当該音響成分は、入射波の周波数ｆ_Ｉと音響周波数ｆ_Ａの倍音との合計（ｆ_Ｉ±Ｈ・ｆ_Ａ、Ｈは０でない整数）に対応する周波数にある。

工程（ａ）では、基準波は、入射波の周波数ｆ_Ｉと音響周波数ｆ_Ａの倍音との合計に等しいか、または略等しい周波数ｆ_Ｒ（ｆ_Ｒ〜ｆ_Ｉ±Ｈ・ｆ_Ａ、Ｈは０でない整数）で生成される。

工程（ｂ）では、撮像物に位置する焦点に集中した音響波が生成され、工程（ｇ）では、測定点の座標が、前記焦点の座標として得られる。

工程（ａ）から（ｇ）は、撮像物における音響波の複数の焦点に対して繰り返し行なわれ、これら複数の焦点は、第１の撮像物方向に沿って配列される。

第１サイクルでは、工程（ａ）から（ｆ）は、音響波の第１周波数ｆ_Ａと基準波の第１周波数ｆ_Ｒとで行なわれ、少なくとも第２サイクルでは、工程（ａ）から（ｆ）は、音響波の第２周波数ｆ’_Ａと基準波の第２周波数ｆ’_Ｒとで繰り返し行なわれ、これら第２周波数は、それぞれ、第１周波数でコード化され、さらに、次の工程を含む。

（ｆ‘）各サイクルの工程（ｆ）において、使用された周波数の関数として得られたデジタル情報をデコードすることによって、少なくとも１つのデジタル情報が得られ、工程（ｇ）では、各サイクルの工程（ｆ）において、使用された周波数の関数として得られたデジタル情報をデコードすることによって、工程（ｆ’）で得られたデジタル情報に関連する、撮像物の少なくとも１つの測定点の座標が得られる。

以下のような一連の操作が実行される。
座標（Ｕ，Ｖ，Ｗ）を持つ点の周囲に、第１の撮像物方向に沿って延びる座標（［Ｕ−Ｄｘ，Ｕ＋Ｄｘ］，Ｖ，Ｗ）を有する点の間に集中する音響波の周波数のスキャンが行なわれ、
ｆ_Ｒがｆ_Ｉ±Ｈ・ｆ_Ａ、（Ｈは０でない整数）に略等しいか、または等しく保ちながら、基準波の周波数ｆ_Ｒでスキャンを共同で行ない、
広がった間隔の点集合（［Ｕ−Ｄｘ，Ｕ＋Ｄｘ］，Ｖ，Ｗ）に対応する干渉画像Ｉ（ｆ_Ａ，Ｖ，Ｗ，ｒ）が、各画素ｒおよび各周波数ｆ_Ａに対して記録され、
１次元周波数→時間フーリエ変換が、干渉画像Ｉ（ｆ_Ａ，Ｖ，Ｗ，ｒ）の周波数ｆ_Ａにしたがって、各画素ｒに対して行われ、
座標（Ｕ’，Ｖ，Ｗ，）の少なくとも１つの測定点に対応する少なくとも１つの干渉画像Ｉ（ｒ）が、フーリエ変換後に得られる時間と、第１の撮像物方向に沿った長さＵ’（Ｕ’はほぼＵに等しい）とを、撮像物における音響波の伝播速度を用いて置き換えることによって得られる。

少なくとも工程（ａ）から（ｇ）を、撮像物の第１の方向（ｘ_０）と平行にはならないように、撮像物に対する振動生成装置の位置を変えた後に繰り返す。

工程（ｆ）では、干渉画像Ｉ（ｒ，ｔ）をもとに、音響成分の複素振幅Ｅ_Ｓ（ｒ）を推定する。

検出装置は単一画素検出装置であり、工程（ｆ）では、対象物によって散乱された複素振幅Ｅ_Ｓ（ｒ）の場の強度として得られる。

検出装置は多画素検出装置であり、工程（ｆ）では、検出装置の画素ｒの、対象物によって散乱された複素振幅場Ｅ_Ｓ（ｒ）の強度の少なくとも一部にわたる合計として、デジタル情報が抽出される。

工程（ｄ）では、検出装置の各画素によって見られる散乱信号波の一部の角度範囲を、少なくとも１つの方向に沿って制限するように、空間フィルタリング装置が使用される（検出装置の各画素によって見られる散乱信号波の一部の、平均角度方向を明確にすることが可能である）。

第１絞り方向に沿った長さがＸで、第３絞り方向に沿った長さがＹの絞りを備える空間フィルタリング装置と、撮像物の下流部に位置する焦点を持つ焦点距離がＬのレンズと、が用いられ、検出装置の各画素によって見られる散乱信号波の一部の角度範囲を制限し、検出装置に入射した基準波が全て平面波である（音響波の入射方向と絞りの方向とは、関連している必要はない）。

撮像物と検出装置との間に、後者から距離Ｌを離して配置され、第１絞り方向に沿った長さがＸで、第３絞り方向に沿った長さがＹの絞りを備える空間フィルタリング装置が用いられ、検出装置の各画素によって見られる散乱信号波の一部の角度範囲を制限し、検出装置に照射される基準波は、絞りの平面に位置する原点から発せられる球面波である（音響波の入射方向と絞りの方向とは、関連している必要はない）。

基準波と散乱信号波とは、ゼロでない角度θ_Ｙをなしながら検出装置上で干渉し、角度θ_Ｙは、検出装置上の２つの波の入射面で測定される。

使用する検出装置は多画素検出装置であり、検出装置の平面で空間的に急激に変化する複素振幅Ｅ_Ｓ（ｒ）の音響成分の一部は、分離される（検出装置は、基準波の方向に対し準直交しており、分離しようとする部分は、干渉画像Ｉ（ｒ）の、空間的には激しく、時間的には緩やかに変化する成分に対応する）。

検出装置は、第１検出装置方向に沿った行、第３検出装置方向に沿った列から成るマトリックスとして配列された画素を備え、
工程（ｆ）では、
（ｆ１）少なくとも１つの行または列に対する１次元フーリエ変換が、検出装置の行または列に沿って、波動ベクトル空間へ、複素振幅Ｅ_Ｓ（ｒ）に対して実行され、この行または列に対して場ＴＦ_１ Ｅ_Ｓ（ｋ）が得られ、
（ｆ２）いくつかの総和領域が、波数ベクトル空間において定義され、
（ｆ３）少なくとも１つの領域の各点ｋでの場ＴＦ_１ Ｅ_Ｓ（ｋ）の強度が、この領域中で総和され、
（ｆ４）デジタル情報は、各領域でこのように得られた総和の線形結合として抽出される（この線形結合は、たった１項のみを含む可能性がある）。

検出装置は、第１検出装置方向に沿った行、第３検出装置方向に沿った列から成るマトリックスとして配列された画素を備え、
工程（ｆ）では、
（ｆ１）検出装置の平面から波数ベクトル空間への、複素振幅Ｅ_Ｓ（ｒ）の２次元フーリエ変換が行われ、場ＴＦ_２ Ｅ_Ｓ（ｋ）が得られ、
（ｆ２）いくつかの総和領域が、波数ベクトル空間において定義され、
（ｆ３）少なくとも１つの領域の各点ｋでの場ＴＦ_２ Ｅ_Ｓ（ｋ）の強度が、この領域中で総和され、
（ｆ４）デジタル情報は、各領域でこのように得られた総和の線形結合として抽出される（この線形結合は、たった１項のみを含む可能性がある）。

角度θ_Ｙは３Ｙ／２Ｌに略等しく、工程（ｆ２）では、いわゆる中央領域の第１総和領域、いわゆる左の領域の第２総和領域、いわゆる右の領域の第３総和領域が定義され、工程（ｆ４）では、デジタル情報は、左の領域の総和と右の領域の総和との線形結合として抽出される（この線形結合は、たった１項のみを含む可能性がある）。

工程（ａ）では、
波長λのレーザ光源は、周波数ｆ_Ｌの放射波を発生し、
放射波の振幅変調手段は、入射周波数ｆ_Ｉの搬送波と、周波数ｆ_Ｒの少なくとも１つの振幅変調側波帯と、を生成し、
半反射鏡は、横帯域波の一部と入射波を形成する搬送波の一部とを伝送し、搬送波の一部と基準波を形成する側波帯の一部とを反射する

工程（ａ）では、
波長λのレーザ光源は、周波数ｆ_Ｌの放射波を発生し、
第１の音響光変調器は、放射波の一部を伝送して、撮像物に対する入射波を形成し、さらに、第１の周波数がシフトされた波を生成し、当該シフト波の周波数は、放射波に対しδｆ_１（負でありうる）だけシフトされたものであり、
第２の音響光変調器は、第１の周波数がシフトされた波を受けて、第２の周波数がシフトされた波を生成し、当該シフト波の周波数は、第１の周波数がシフトされた波に対しδｆ_２（負でありうる）だけシフトされたものであり、基準波（ＲＥＦ）を形成する第２の周波数がシフトされた波の周波数は、入射波（ＩＮＣ）に対しδｆ＝δｆ_１＋δｆ_２だけシフトしたものであり、このようにして２つの波の位相差φ_ｉ（ｔ）を決定する（大抵、δｆ_１とδｆ_２は反対の符号となる）。

工程（ａ）では、電子的スレービングによって位相が固定された２つの独立なレーザ光源は、これらの間で既知の位相差φ_ｉ（ｔ）を示す入射波と基準波とを生成する。

工程（ａ）では、
波長λのレーザ光源は、周波数ｆ_Ｌの放射波を発生し、
半反射鏡は、放射波の一部を伝送して、撮像物に対する入射波を形成し、放射波の第２の部分を伝送し、
第１の音響光変調器は、放射波の第２の部分を受けて、当該放射波に対しδｆ_１（負でありうる）だけシフトした第１の周波数がシフトされた波を生成し、
第２の音響光変調器は、第１の周波数がシフトされた波を受けて、シフト波に対しδｆ_２（負でありうる）だけシフトされた第２の周波数がシフトされた波を生成し、基準波を形成する第２の周波数がシフトされた波の周波数は、入射波に対しδｆ＝δｆ_１＋δｆ_２だけシフトされたものであり、このようにして２つの波の位相差φ_ｉ（ｔ）を決定する（大抵、δｆ_１とδｆ_２は反対の符号となる）。

撮像物は、生物学的組織である。

振動生成装置は、撮像物の領域の音響キューを得るために用いられ、工程（ｆ）で抽出されたデジタル情報は、この音響キューと共同して使用される。

他の態様によれば、本発明は、撮像物に対する音響光学撮像装置に関連し、この音響光学撮像装置は、
周波数ｆ_Ｉの入射光波、および、この入射光波とコヒーレントでかつ既知の位相差φ_ｉ（ｔ）を示す周波数ｆ_Ｒの基準光波（ＲＥＦ）を生成する手段と、
第１の撮像物方向に沿って、音響周波数ｆ_Ａで撮像物の領域の振動する振動生成装置と、
撮像物に入射光波を照射し、これにより散乱信号波を生成する手段と、
検出装置と、
散乱信号波の少なくとも一部を検出装置に照射する手段と、
基準波を撮像物に透過させずに検出装置（ＤＥＴ）に照射し、検出装置の点ｒに時間ｔで変化する干渉画像Ｉ（ｒ，ｔ）生成する手段と、
デジタル情報、および、当該デジタル情報が関連する撮像物の測定点の座標を干渉画像から抽出する手段と、を備える。

本発明の好ましい形態において、以下の条件によるものである、

音響光学撮像装置は、以下の要素をさらに備える。
前記撮像物の測定点に関連するデジタル情報を可視化する手段。
撮像物を動かす手段。

音響光学撮像装置は、さらに、撮像物の下流部に位置する空間フィルタリング装置を備える。

本発明の、他の側面、目的および効果は、限定しない例によって与えられた実施形態の記載を読むことで明らかになるであろう。

本発明によれば、局所振動子として働く基準波は、撮像物を透過する必要はない。これにより、充分なレベルの局所振動子が得られ、よりよい信号／ノイズ比で撮像のための測定点、特に医療のための撮像の測定点に関する有用な情報を得ることができる。更に、この撮像方法は、例えば、医療用画像に関連した撮像した組織の安全基準にかなった低い音響又は光学出力であっても有用な信号を得られる。

図１は、波を生成する装置ＧＥＮを示し、この装置は、
撮像物ＯＢＪに照射される周波数ｆ_Ｉ、波長λの入射光波ＩＮＣと、
周波数ｆ_Ｒの基準光波ＲＥＦと、を生成する。

この入射波ＩＮＣおよび基準波ＲＥＦは、互いにコヒーレントであり、既知の位相差φ_ｉ（ｔ）を示す。これらの光波は、可視領域、あるいは赤外領域や紫外領域に放射される。

生成装置ＧＥＮは、基準光波ＲＥＦの周波数が、入射波に対してδｆに等しい値だけシフトされるように調整される。

このような生成装置を実現する可能性は多数存在しており、例えば、図２に示すような生成装置を用いることが可能である。この生成装置は、
周波数ｆ_Ｌで放射光波ＥＭＩを放射するレーザ源ＬＡＳと、
放射波ＥＭＩの一部を伝送して、周波数ｆ_Ｉで入射波ＩＮＣを形成し、さらに、周波数がδｆ_１だけシフトされた波ＤＥＣを生成する第１音響光学変調器ＭＡＯ１と、
周波数がｆ_Ｌ＋δｆ_１にシフトされた波ＤＥＣを受けて、再度この波を周波数δｆ_２だけシフトし、したがって、周波数がｆ_Ｒ＝ｆ_Ｌ＋δｆ_１＋δｆ_２（δｆ_１および／またはδｆ_２は、負でありえ、また、δｆ_１およびδｆ_２は、通常反対の符号である）の基準波ＲＥＦを生成する第２音響光学変調器ＭＡＯ２と、を備える。

これらの音響光学変調器ＭＡＯ１およびＭＡＯ２は、例えば、これらに照射される波、すなわち放射波ＥＭＩおよび位相波ＤＥＣに対し所定の角度で向けられた二酸化テルル（ＴｅＯ_２）から成る音響光学素子によって構成され、高周波生成器の下で周波数δｆ_１およびδｆ_２でそれぞれ振動し、非回折光と、周波数がシフトされた回折光と、の両方を伝送する。これら２つの波が互いに、数ヘルツから数十メガヘルツの間に固定された周波数δｆ＝δｆ_１＋δｆ_２にシフトされるように、レーザ波に対して基準波ＲＥＦの周波数をシフトさせる２つの音響光学変調器ＭＡＯ１とＭＡＯ２は、８０ＭＨｚの周波数を有する。

他の変形例としては、第１の音響光学変調器ＭＡＯ１の上流に置かれた半反射鏡（又はいかなる他の装置）によって、レーザ源ＬＡＳから入射波ＩＮＣが抽出され、その後、この半反射鏡（レーザ源ＬＡＳと第１の音響光学変調器との間に配置した他の同様な装置）によって伝送された入射波の部分に基づいて、シフトされた波ＤＥＣを生成する、その他の生成装置（図示せず）を用いることが可能である。

また、電子的スレービングによって位相が固定され、入射波ＩＮＣと基準波ＲＥＦとの間で既知の位相差φ_ｉ（ｔ）を示す２つの独立なレーザを用いることも可能である。

生成装置のこれら実施形態において、レーザ光源または光源ＬＡＳは、例えば、回折格子によって安定化され、波長λ＝８５０ｎｍのレーザ波ＥＭＩを放出し、６５ｍＡの電流に対し２０ｍＷの最大出力を供給するレーザダイオードを備える。略１５ｍＷの出力を有するレーザを用いることが可能であり、したがってこの例では、光学損失を考慮すると、略７ｍＷの出力の、撮像物に対する入射波ＩＮＣが与えられる。

したがって、周波数ｆ_Ｉの入射波ＩＮＣおよび周波数ｆ_Ｒの基準波ＲＥＦの、取得された２つの波は、周波数がδｆ＝δｆ_１＋δｆ_２だけシフトされた２つのコヒーレントな波である。それゆえ、これら２つの波の間の位相差φ_ｉ（ｔ）を制御できる。

図５は、図２の発生装置ＧＥＮを備える本発明の第１実施形態を示す。

入射波ＩＮＣが照射される撮像物ＯＢＪは、光波に対する散乱物であり、例えば、生物学的組織のサンプルである。この例では、サンプルは、入射波ＩＮＣの伝播方向に沿って略２０ｍｍの厚さを示す。

特に、このサンプルは、入射波ＩＮＣの伝播方向に対し、垂直である前面および下流部平板の間で圧縮され、これら２つの平板は、サンプルホルダ（図示せず）の一部を形成する。

この上流部平板は、例えば、全く透明であり、より詳細にはＤＭＭＡ（プレキシガラス（登録商標））から作られ、下流部平板は、不透明であり、より詳細には黒いベークライトから作られる。この後部平板は、例えば直径が略Ｘ＝２０ｍｍの円形穴が開けられる。

音響的な結合は、撮像物ＯＢＪと、この撮像物の外に配置した振動生成装置ＴＲＡＮＳと、の間で生成される。この目的のため、撮像物ＯＢＪとそのサンプルホルダとは、例えば、直径８０ｍｍで高さ１５０ｍｍのタンク１の中央に設置される。このタンクは、直径６０ｍｍの平面ガラス窓を、２４０ｍｍ離して具備されてもよく、また、撮像物ＯＢＪと振動生成装置ＴＲＡＮＳとの間の音響結合の実現を果たすための液体で満たされている。

何らかの理由により、撮像物ＯＢＪが浸されない場合は、この音響結合を、当業者に知られる他の手段で保証してもよい。例えば、振動生成装置を、直接撮像物ＯＢＪに搭載できる。

振動生成装置ＴＲＡＮＳは、音響振動の周波数がｆ_Ａ＝２．２ＭＨｚで、例えば７５ｍｍに設定された変化可能な焦点長を有するトランスデューサＰＺＴである。これは、例えば、５０Ωに２２ｄＢｍ、すなわち６３Ｖｃｃの最大出力の正弦波信号によって励起し、これは生物学的組織に対し良性を示す。例えば、３４ｄＢｍ程度の、異なる最大出力で使用することが可能である。

この振動生成装置は、撮像物の第１の方向ｘ_ｏに沿って向けられ、ｘ_ｏに沿って周波数ｆ_Ａの音響波を放出する。情報を得ようとする座標点（Ｕ，Ｖ，Ｗ）を中心にした撮像物の領域（Ｄｘ，Ｄｙ，Ｄｚ）は、振動生成装置ＴＲＡＮＳによって、周波数ｆ_Ａで振動される。振動する撮像物の領域の範囲（Ｄｙ，Ｄｚ）は、概ね、横方向ｙおよびｚ、すなわち、トランスデューサによって放出される音響波の伝播方向ｘ_ｏと垂直な平面に沿った音響波の焦点領域のサイズに対応する。この平面において、焦点領域は、（Ｖ，Ｗ）を中央に位置付けられる。Ｕは、音響トランスデューサと音響波の伝播方向に沿う後者の焦点との間の距離に対応する。この焦点の周囲において、トランスデューサの種類および撮像物の特性に依存する範囲Ｄｘの領域は、特に、周波数ｆ_Ａで振動する。したがって、トランスデューサの位置および方向とその焦点の位置とが、撮像物ＯＢＪの測定点を、座標（Ｕ，Ｖ，Ｗ）に決定する。

撮像物ＯＢＪによって全方向に散乱される散乱信号波ＤＩＦを形成するために、第２の撮像物方向ｚ_ｏ、できれば第１の撮像物方向ｘ_ｏと同一方向に沿って、撮像物ＯＢＪに入射波が照射される。撮像物ＯＢＪの内部では、波の一部が、音響周波数ｆ_Ａで振動する撮像物の領域（Ｄｘ，Ｄｙ，Ｄｚ）を通過する。散乱しそうな撮像物の点の動きが、音響周波数ｆ_Ａで散乱された波の位相の変調を発生する。この振動は、さらに、媒質の光学定数（周波数ｆ_Ａで）の変調を生成する。これらの２つの効果により、媒質を通る周波数ｆ_Ｉの入射波に対し、周波数ｆ_Ａによって周波数が変調された音響側波帯を生成する。したがって、散乱信号波ＤＩＦは、周波数ｆ_ＡＯ＝ｆ_Ｉ＋Ｈ・ｆ_Ａ（ここで、Ｈ＝１，２，．．．は、倍音のランクであり、一般にＨ＝１である）の、シフトされた周波数の音響成分を有する。音響光撮像の原理は、それから、この散乱信号波の音響成分を、ほぼ等しい周波数の局所振動子と干渉させることによって、撮像物によって散乱された信号波の音響成分を検出装置ＤＥＴで検出することから成る。これは、本発明に基づいて、基準波ＲＥＦを、撮像物に透過させずに、検出装置の局所振動子として用いることにより実行される。

検出装置は、観測方向ｚ_Ｄ（図における例では、ｘ_Ｄ，ｙ_Ｄ，ｚ_Ｄがｘ_Ｏ，ｙ_Ｏ，ｚ_Ｏに対応する）と略直行する平面ｘ_Ｄ，ｙ_Ｄ上で、少なくとも１つの検出セルを示す。

例えば、ｄｘ＝ｄｙ＝６．７μｍのサイズの１２８０×１０２４の四角の画素から成るデジタルＣＣＤカメラ（１２ビット）のような多画素検出装置ＤＥＴを使用できる。このＣＣＤカメラは、例えば、８５０ｎｍの波に対して、充分な量（例えば、５％より大きい）を有するように選択される。このカメラは、検出周波数が例えばｆ_Ｃ＝１２．５Ｈｚの“フルフレーム”型（ノンインターレース）、またはフレーム転送型でもよい。

撮像物ＯＢＪから放射する散乱信号波ＤＩＦを、検出装置ＤＥＴ上で基準波ＲＥＦと干渉させる装置が使用される。有限サイズの画素上で有効な検出を行なうために、２つの波がほぼ同一線（形成する角度はせいぜい５°）である必要がある。例えば、検出装置ＤＥＴへ向けて基準波ＲＥＦを導くために、半反射鏡板やスプリッタプリズムを使用できる。

基準波ＲＥＦの方向を、散乱信号波ＤＩＦの平均方向に対し、角度θ_Ｙだけずらすことが有効である。以下の記載と図１、５、７および９とにおいて、角度のずれは方向ｙに沿って行われ、音響波の方向はｘで、観測方向は入射波ＩＮＣ（方向ｚ）方向と同じで、ＣＣＤ検出装置の行列はｘおよびｙに沿って方向付けられるが、これらの方向は一例として定義されたものである。

散乱信号波ＤＩＦおよび基準波ＲＥＦは、検出装置上で干渉するようになされ、この後者によって、干渉画像Ｉ（Ｕ，Ｖ，Ｗ，ｒ，ｔ）が、検出装置上の点ｒの瞬間ｔにおいて撮影され、これが、音響周波数ｆ_Ａで振動する撮像物の座標（Ｕ，Ｖ，Ｗ）における測定点に対応する。

この装置は、記録された時間干渉画像から、座標（Ｕ，Ｖ，Ｗ）の測定点に関するデジタル情報を抽出できる処理装置ＣＡＬＣに接続される。このデジタル情報は、その後、撮像物の画像に表示される。この処理は、検出装置上で音響振動により周波数がシフトされた散乱された領域の音響成分の複素振幅Ｅｓ（Ｕ，Ｖ，Ｗ，ｒ）の計算を含む。

第１実施形態に限らず、複素振幅Ｅ_Ｓ（Ｕ，Ｖ，Ｗ，ｒ）の決定において、固定周波数ｆ_Ａの音響波は、座標（Ｕ，Ｖ，Ｗ）の測定点に集中され、基準波ＲＥＦの周波数ｆ_Ｒは、以下に説明されるように、Ｎ相復調、特に４相復調（Ｎ＝４）を実行するように選択される。

音響周波数ｆ_Ａで振動する撮像物ＯＢＪの座標（Ｕ，Ｖ，Ｗ）の測定点において、散乱信号波ＤＩＦは、周波数ｆ_Ｉ＋ｆ_Ａの音響成分を含む。処理ユニットＣＡＬＣは、時間変化する干渉画像Ｉ（Ｕ，Ｖ，Ｗ，ｒ，ｔ）を、４相復調によって、以下の方法で処理する。まず、周波数のズレδｆは、入射波ＩＮＣと基準波ＲＥＦとの間で、δｆ＝ｆ_Ａ＋ｆ_Ｃ／４＝ｆ_Ａ＋ｆ_Ｃ＝ｆ_Ａ＋３．１２５Ｈｚになるように選択される。しかしながら、この検出はまた、少なくとも２以上の任意のＮの相に対して実行でき、基準波の周波数は、δｆ＝ｆ_Ａ＋ｆ_Ｃ／Ｎとすることによって合わせられる。

その後、Ｎに対する干渉画像を、それぞれの時間Ｔ_ｉｎｔの間、また、入射波ＩＮＣと基準波ＲＥＦとの既知の位相差φ_ｉに対応するそれぞれの干渉画像の測定をする。基準波が周波数ｆ_Ｒ＝ｆ_Ｉ＋ｆ_Ａ＋ｆ_Ｃ／Ｎであり、信号波の音響成分が周波数ｆ_Ｉ＋ｆ_Ａである場合、基準周波数ＲＥＦと信号波の音響成分との位相差は、一般的な方式で図３に示されているが、Ｎが４に対する図４では、時間の関数として区間ごとに線形であり、Ｎ／ｆ_Ｃに等しい時間間隔にわたり、０から２πになる。

そのため、Ｎ個の異なる位相差の値φ_ｉに対応するＮ個の干渉画像Ｉ_１ ^１（Ｕ，Ｖ，Ｗ，ｒ，ｔ），．．．，Ｉ_ｉ ^１（Ｕ，Ｖ，Ｗ，ｒ，ｔ），．．．，Ｉ_Ｎ ^１（Ｕ，Ｖ，Ｗ，ｒ，ｔ）を検出する。Ｎは少なくとも２以上であり、図４の場合はＮが４に等しい。

この操作は、それぞれの既知の位相差φ_ｉに対し、ｎ個の干渉画像Ｉ_ｉ ^ｋ（Ｕ，Ｖ，Ｗ，ｒ，ｔ）（ｉ＝１．．．Ｎ，ｋ＝１．．．ｎ）を求めるために、ｎ（少なくとも回数１、例えば、ｎ＝３）回実行されてもよい。
１２個の干渉画像、すなわち、４相、３サイクル（Ｎ＝４，ｎ＝３）に対しては、測定の全持続時間は、例えば、数秒のオーダーである。

４相復調の計算は、検出装置の各画素ｒに対し、別個に実行される。処理装置ＣＡＬＣは、実際には、検出装置の各画素ｒに対し、以下の操作を実施する。
ＣＣＤカメラ（１／１２．５秒）のフレームの持続期間にわたり、検出装置によって測定された強度の積分 Ｉ_ｉ ^ｋ（Ｕ，Ｖ，Ｗ，ｒ）＝（１／Ｔ_ｉｎｔ）∫Ｉ（Ｕ，Ｖ，Ｗ，ｒ，ｔ）ｄｔ（時間ｔ＝ｔ_ｉｋとｔ＝ｔ_ｉｋ＋Ｔ_ｉｎｔとの間にとられる）。
Ｎ個の平均干渉画像Ｉ_ｉ（Ｕ，Ｖ，Ｗ，ｒ）を得るために、既知の各位相φ_ｉに対し、この位相に対して検出されたｎ個の干渉画像Ｉ_ｉ ^ｋ（Ｕ，Ｖ，Ｗ，ｒ） （ｋ＝１．．．ｎ）の平均をとる。
撮像物により散乱された信号波ＤＩＦの音響成分の複素振幅Ｅ_ｓ（Ｕ，Ｖ，Ｗ，ｒ）に関する情報を得るように測定されたＮ個の干渉画像のＮ相復調。例えば、ｎが１に等しく、また、特に、Ｎが４に等しく、φ_１＝π／２，φ_２＝π，φ_３＝３π／２，φ_４＝２πとした図４の場合、基準波は空間と時間において一定の複素振幅を示すと考えられので、音響成分の複素振幅Ｅ_Ｓ（Ｕ，Ｖ，Ｗ，ｒ）は、（Ｉ_４−Ｉ_２）＋ｊ（Ｉ_１−Ｉ_３）（ここで、ｊは、ｊ^２＝−１のような複素数である）と、撮像物ＯＢＪに伴う強度と、に比例する。

散乱信号波ＤＩＦの音響成分の複素振幅Ｅ_Ｓ（Ｕ，Ｖ，Ｗ，ｒ）を決定するため、当業者に知られた、いくつかの他の技術が、一般に使用される。

音響成分の複素振幅Ｅｓ（Ｕ，Ｖ，Ｗ，ｒ）の決定の第２の形態は、いわゆる“周波数チャーピング”処理の使用から構成される。この場合、周波数ｆ_Ａに固定された音響波を、方向ｘに沿った長さＵに集中させる代わりに、次の一連の動作を実行する。
Ｕの周囲に広がる［Ｕ−Ｄｘ，Ｕ＋Ｄｘ］の区間に集中された周波数ｆ_Ａの音響波のスキャンが行ない、
条件ｆ_Ｒ＝ｆ_Ｉ±Ｈ・ｆ_Ａ（ここで、Ｈは倍音のランクであり、一般に１である）を保持するように、周波数ｆ_Ｒの基準波のスキャンを、第１のスキャンと共同で行ない（後に注意するように、基準波は、検出しようとする音響成分の周波数と同じ周波数を有する）、
広がった区間の点集合（［Ｕ−Ｄｘ，Ｕ＋Ｄｘ］，Ｖ，Ｗ）に対応する干渉画像Ｉ（ｆ_Ａ，Ｖ，Ｗ，ｒ）が、位置ｒの各画素、および各周波数ｆ_Ａに対して記録され、
周波数→時間の１次元フーリエ変換が、位置ｒの各画素に対して、干渉画像Ｉ（ｆ_Ａ，Ｖ，Ｗ，ｒ）の周波数ｆ_Ａによって行われ、
方向ｘに沿った長さＵの、様々な値に対する複素干渉画像Ｉ（Ｕ，Ｖ，Ｗ，ｒ）が、撮像物における音響波の伝播速度を用いて、フーリエ変換後に得られる時間をｘ方向に沿った長さＵに置き換えることにより得られる。

基準波は、ほぼ一定の振幅の平面波であり、デコードされた複素干渉画像Ｉ（Ｕ，Ｖ，Ｗ，ｒ）は、決定しようとする散乱信号波ＤＩＦの音響成分の複素振幅Ｅｓ（Ｕ，Ｖ，Ｗ，ｒ）に比例する。

ｘ方向に沿ったコーディングを、ｆ_Ａとｆ_Ｒにしたがってコーディングした周波数に置き換えることからなる“周波数チャーピング”技術には、いくつもの変形例がある。これらの技術や他の技術は、本発明の記載において、Ｅｓ（Ｕ，Ｖ，Ｗ，ｒ）を決定するために使用する。

この第１実施形態では、さらに、検出装置の各画素で見られる信号波の一部の角度範囲を制限することを可能にする空間フィルタリングＣＯＬ（図１、５）の設置がなされる。この装置は、検出装置の平面において、散乱信号波ＤＩＦに伴うスペックル粒のサイズを制御するために有用である。

このような検出を行うには、スペックル粒のサイズが、カメラの画素のサイズに合わせられることを確認しておくことは有益である。この第１条件は、いわゆる“反エーリアシング”条件に対応する。

さらに、この空間フィルタリング装置の構造を慎重に選択することで、各種のノイズ成分から有用な信号を、本発明にしたがって分離することが可能である。

このフィルタリング装置は、例えば、矩形状の絞り２によって構成される。この絞り２は、ある方向に沿って伸びており、例えば、観測方向に対し垂直で、撮像物ＯＢＪの下流部に置かれ（したがって検出装置に対し略平行である）、例えば、撮像物ＯＢＪとサンプルホルダの下流部の板との間に置かれる。この絞り２は、例えば、０．５ｍｍの厚さのアルミニュームシートで、略４ｍｍで分けられた２枚の薄板から成る。観測可能な領域の上流部は、したがって、観測方向に対し垂直な２つの方向に沿って、略ＸｍｍとＹｍｍの形状の準矩形である。例えば、観測方向としてｚを用い、絞りの軸としてｘおよびｙを用いる。

図５に示された空間フィルタリング装置ＣＯＬの実施形態において、この装置ＣＯＬは、さらに、タンク１と検出装置ＤＥＴとの間に配置されるレンズ３を備える。このレンズ３の対象物への焦点は、絞りの平面上に設定される（散乱信号波ＤＩＦが、振動生成装置ＴＲＡＮＳの音響結合を用いられる液体を通過する場合は、この液体の光学定数も考慮することが好ましい）。ここに記載された実施形態においては、焦点距離Ｌ＝２５０ｍｍのレンズが使用されるが、他の焦点距離であっても適用できる。

この空間フィルタリング装置は、検出装置ＤＥＴに到達する散乱信号波ＤＩＦの一部ＳＩＧの角度の範囲を減少させる。これは、スペックル粒のサイズを、カメラの画素のサイズに合わせるのに有用である。

特に、空間フィルタリング装置を使用しない場合、撮像物ＯＢＪで散乱された信号波ＤＩＦは、πのオーダーの広域な立体角に拡散し、大きく異なった波数ベクトルＫ _Ｓ を有する平面波の重ね合わせに分解される。各波数ベクトルＫ _Ｓ は、検出装置（ｘ_Ｄ，ｙ_Ｄ）の平面上で、２つの座標Ｋ_ｘおよびＫ_ｙを有する。有用な情報を得るためには、散乱信号波の基準波による干渉の信号の空間変調に対応する縞の密度が、複数の基礎的な検出装置で構成されたマトリックス検出装置の解像度を超さないようにする必要がある（いわゆる“反エーリアシング”条件）。特に、検出装置の平面上の座標Ｋ_ｘｏ，Ｋ_ｙｏで、波数ベクトルＫ _ｏ を有する平面基準波ＲＥＦに対しては、この縞の密度の条件は、以下に定義される“反エーリアシング”条件を導く。

ＳＫ０＝ＳＩＮＣ（ｄ_ｘ・（Ｋ_ｘ−Ｋ_ｘ０））・ＳＩＮＣ（ｄ_ｙ・（Ｋ_ｙ−Ｋ_ｙ０））〜１−Ｅａ。ここで、ダミー変数ｘｘの関数ＳＩＮＣは、ｘｘ＝０に対して１に等しく、それ以外はｓｉｎ ｘｘ／ｘｘに等しく、Ｅａは、縞のコントラストの損失を数量化した減衰因子であり、検出装置の空間積分と関連する。また、ｄｘとｄｙは、検出装置を構成する基礎的な検出装置の、方向ｘ_Ｄおよびｙ_Ｄに沿った特徴的なサイズを示す。したがって、この測定は、基準波ＲＥＦの波数ベクトルＫ _ｏ の方向を取り囲む角度（±α_ｘ；±α_ｙ）の円錐に対応する、信号波ＳＩＧの基本角の領域に制限しなければならない。基本角の領域のサイズα_ｘおよびα_ｙは、前記“反エーリアシング”条件を満たすように、それぞれ実質的に、λ／２ｄ_ｘおよびλ／２ｄ_ｙに、少ないか等しくされる。ここで、λは基準波長ＲＥＦである。

撮像物ＯＢＪの下流部の絞りを使用することにより、特に、“反エーリアシング”条件を満たさない信号波の成分を除去することができる。

この空間フィルタリング装置ＣＯＬと検出装置との構造を慎重に選択することで、干渉画像Ｉ（Ｕ，Ｖ，Ｗ，ｒ，ｔ）の解析から得られる信号に現れる様々な項から、有用な信号をさらに分離できる。この問題は、４相復調の場合のみ論じられているが、“周波数チャーピング”や他の同種の技法による場合においても同様に論ずることができる。

干渉画像Ｉ（Ｕ，Ｖ，Ｗ，ｒ，ｔ）は、検出装置で見られる全強度Ｉ_Ｔに対応する。すなわち、複素振幅Ｅの絶対値の２乗に相当する（Ｉ_Ｔ＝｜Ｅ｜^２Ｅ・Ｅ^＊。ここで、Ｅ^＊はＥの複素共役である）。議論を簡潔にするために、周波数ｆ_ＡＯ＝ｆ_Ｉ＋Ｈ・ｆ_Ａ（Ｈ＝１）の、単一音響成分のみを検討する。複素振幅Ｅは、基準波ＲＥＦの振幅Ｅ_Ｒと、入射波の周波数ｆ_Ｉで散乱された信号波ＳＩＧの一部の振幅Ｅ_Ｉと、周波数ｆ_ＡＯの信号波の音響成分の振幅Ｅ_Ｓと、の合計から得られる。したがって、Ｅ＝Ｅ_Ｉ＋Ｅ_Ｒ＋Ｅ_Ｓを得る。

求めたい情報を検出するのに用いられる全強度Ｉ_Ｔは、６つの項の合計に対応する（Ｉ＝Ｅ_Ｉ・Ｅ_Ｉ ^＊＝（Ｅ_Ｉ＋Ｅ_Ｒ＋Ｅ_ｓ）・（Ｅ_Ｉ＋Ｅ_Ｒ＋Ｅ_ｓ）^＊）。
Ｅ_Ｉ・Ｅ_Ｉ ^＊は、周波数ｆ_Ｉで散乱された信号波と自身との干渉に対応する。すなわち、通常スペックルと通常スペックルとの干渉に対応し、
Ｅ_Ｓ・Ｅ_Ｓ ^＊は、音響光学周波数ｆ_ＡＯ＝ｆ_Ｉ＋ｆ_Ａで散乱された信号波の音響成分と自身との干渉に対応する。すなわち、音響光学スペックルと音響光学スペックルとの干渉に対応し、
Ｅ_Ｒ・Ｅ_Ｒ ^＊は、基準波と自体との干渉に対応し、
Ｅ_Ｒ・Ｅ_Ｉ ^＊は、基準波と通常スペックルとの干渉に対応し、
Ｅ_Ｉ・Ｅ_Ｓ ^＊は、通常スペックルと音響光学スペックルとの干渉に対応し、
Ｅ_Ｒ・Ｅ_Ｓ ^＊は、基準波と音響光学スペックルとの干渉に対応し、関連情報を伝送する項で構成される。

空間フィルタリング装置ＣＯＬは、撮像物から放射される波の角度の範囲を狭くし、これにより、この波は、検出装置のレベルでは準平面波として振舞う。これは、特に、複素振幅Ｅ_ＩおよびＥ_Ｓを有する散乱信号波の一部に対する場合である。また、本実施形態において、基準波ＲＥＦは平面波である。空間フィルタリング装置によって、対角項Ｅ_Ｒ・Ｅ_Ｒ ^＊、Ｅ_Ｉ・Ｅ_Ｉ ^＊、およびＥ_Ｓ・Ｅ_Ｓ ^＊（また、Ｅ_Ｉ・Ｅ_Ｓ ^＊も同様）は、検出装置の方向ｘおよびｙに沿って、空間的にゆっくりと変化する。さらに、Ｅ_Ｒ・Ｅ_Ｉ ^＊は、時間的に速く変化し（ほぼｆ_Ａの周波数で）、検出装置の低取得周波数のため、Ｅ_Ｒ・Ｅ_Ｉ ^＊は０に平均化される。したがって、関連した項Ｅ_Ｒ・Ｅ_Ｓ ^＊は、適当なデジタル処理で容易に分離できる（これにより、Ｅ_Ｓを決定できる）。基準波と散乱信号波の間に、充分な角度のずれθ_Ｙを選択する場合、この項は、時間的にゆっくり変化し、Ｙ方向に沿って空間的に速く変化する唯一の項となる。このＹ方向は、絞りの幅方向と、検出装置の平面の方向Ｙ_Ｄと、に対応する。

関連情報を抽出する１つの方法は、上記で計算される複素振幅Ｅ_Ｓ（Ｕ，Ｖ，Ｗ，ｒ）の、検出装置の平面ｘおよびｙ（あるいはｙ方向のみ）に沿ったフーリエ変換を実行することからなる。その結果、信号ＴＦ Ｅ_Ｓ（Ｕ，Ｖ，Ｗ，ｋ）を得る。ｋは波数ベクトル空間の座標である。フーリエ変換後得られた信号ＴＦ Ｅ_Ｓ（Ｕ，Ｖ，Ｗ，ｋ）のマッピングは、図６に示されており、波数ベクトル空間で角表示されたものである。ここに示された干渉画像に寄与する各種の項は、以下で説明される。

全時間にわたって完全に安定して試験が行われるならば、Ｎ相復調（または“周波数チャーピング”）により、理論上は、基準波と基準波との干渉項（Ｅ_Ｒ・Ｅ_Ｒ ^＊）を、完全に取り除くことができる。これは、あり得ない状況であり、極めて深刻な不快な成分が残る。しかし、この項は、検出装置の平面ｘおよびｙの方向に沿ってゆっくり変化するので、したがって、２次元フーリエ変換により、波数ベクトル空間において座標軸の原点に位置する狭いピークになる（図６の領域４）。

全時間にわたって完全に安定して試験が行われ、かつ、スペックルが静的で、自己無相関でないならば、Ｎ相復調（または“周波数チャーピング”）により、通常スペックルと通常スペックルとの干渉項（Ｅ_Ｉ・Ｅ_Ｉ ^＊）を、取り除くことができる。これは、あり得ない状況であり、したがって、極めて深刻な不快な成分が残る（図６の領域２）。この項は、撮像物への主なノイズ項であり、スペックルは静的ではない（例えば、生物学組織に対する）。通常スペックルと通常スペックルとの干渉項は、基準波と基準波との干渉項と同様に、波数ベクトル空間の原点に集中する。先に定義した空間フィルタリング装置によって、通常スペックルの領域の振幅は、［−Ｙ／２Ｌ；Ｙ／２Ｌ］の間隔に対応する有限の角度範囲を有する。通常スペックルと通常スペックルとの干渉項は、場の強度、すなわち、複素振幅と共役複素振幅との積に相当するので、波数ベクトル空間のその範囲を評価するには、波数ベクトル空間においてスペックル領域を自身でかき回す必要がある。したがって、この干渉項は、それ自身の領域の２倍の角度範囲を有する（間隔［−Ｙ_２＝−Ｙ／Ｌ；Ｙ_２＝＋Ｙ／Ｌ］）。同様に、最適な（Ｘ／Ｌ＝λ／２ｄ_ｘ）において、後者がエーリアシング条件を満たされれば、波数ベクトル空間の高さ方向には、この項は間隔［−Ｘ_２＝−Ｘ／Ｌ；Ｘ_２＝＋Ｘ／Ｌ］に相当する。このノイズは、したがって、範囲が定められた錐体状の空間を示す。空間は、波数ベクトル空間の座標の原点の中央に位置する。

音響光学スペックルと音響光学スペックルとの干渉（Ｅ_Ｓ・Ｅ_Ｓ ^＊）は、２次の対角項である。その低強度は別として、この項は、上述の通常スペックルと通常スペックルとの干渉項と異ならない。

基準波は、周波数ｆ_Ｒであり、通常スペックルは、周波数ｆ_Ｉである。基準波と通常スペックルとの干渉項（Ｅ_Ｒ・Ｅ_Ｉ ^＊）は、周波数が略ｆ_Ｒ−ｆ_Ｉであり、ｆ_Ａに略等しく、すなわち、略２．２ＭＨｚである。この干渉項は、検出装置の低取得周波数のため、取得時間にわたってゼロに平均化され、したがって省略してもよい。

同様に、通常スペックルと音響光学スペックルとの干渉項（Ｅ_Ｉ・Ｅ_Ｓ ^＊）は、２次の項であることに加えて、音響波の周波数ｆ_Ａと略等しい周波数を備える。したがって、各画像の取得時間にわたってゼロに平均化してもよい。したがって省略してもよい。

信号波ＳＩＧと基準波ＲＥＦとの干渉から抽出される有効な項は、したがって、基準波と音響光学スペックルとの干渉項（Ｅ_Ｒ・Ｅ_Ｓ ^＊）である。この項は、領域３と図６の間隔［Ｙ_３−；Ｙ_３＋］に対応する。この項は、絞り２から放射するスペックルと平面波との干渉に対応しているので、間隔［Ｙ_３−；Ｙ_３＋］は、絞りの画角、すなわちＹ_３−−Ｙ_３＋＝Ｙ／Ｌ、と同じ角度幅を有する。この項は、検出周波数ｆ_Ｒ−ｆ_ＡＯ（＝ｆ_Ｃ／４（４位相復調の場合）、＝０（“周波数チャーピング“の場合））で時間発展する。したがって、それは電子検出装置によって計測されてもよい。さらに、基準波は、角度的にずらされるので、領域３の中央は、θ_Ｙだけ領域３の中央がずれる（ここで、先に定義したように、角度θ_Ｙは、基準波ＲＥＦと信号波ＳＩＧとの間の角度である）。幅Ｙの絞りを使用する時、角度θ_Ｙの選択は、領域３の中心の位置を定義する。また、幅がＹの絞りを用いた場合、領域３の範囲は、いずれにせよＹ／Ｌとなる。したがって、領域３の画素において、役立つ情報だけが得られるように、領域２と３は重複しないように工夫するのが適当である。

例えば、θ_Ｙ＝３Ｙ／２Ｌを選択すると、Ｙ_３＋＝２Ｙ／ＬとＹ_３−＝Ｙ／Ｌ＝Ｙ_２が得られる。したがって、これら２つの領域が重複しないで、領域２の外側端と領域３の内側端とが重ね合わせられている。それから、図６の領域２および３の間には、ブランクの部分がない。したがって、基準波と音響光学スペックルとの干渉信号は、波数ベクトル空間において、この波数ベクトル空間の中心に位置する基準波の強度およびスペックルの強度（図６の領域４および２）から分離される。

領域３の範囲が、Ｙに直接比例しているので、“エーリアシング”条件Ｙ_３＋＜λ／２ｄ_ｙを満たされるまで、絞りの幅を広げようとする誘惑が存在する。しかしながら、一方で絞りの幅Ｙを広げることで領域３を最大にすることと、もう一方でこの面積と有効性との積の、この領域にわたる積分を最大にすることと、の間でなされる妥協が存在する。この有効性は、有限サイズの画素にわたる干渉画像の積分による、干渉画像のコントラスト損失に対応する。隣接する画素に対し、有効因子は、前に記載した“エーリアシング”の条件と同様のＳＩＮＣ則に従って見られる。この場合、“エーリアシング”角λ／２ｄ_ｙに達すると、有効性がなくなる。例えばＹ_３＋は、“エーリアシング”の極限よりもかなり小さく選択される。有効性が十分なように、各側で検出装置の略第６領域を活動しないように“エーリアシング”の極限、例えば、未満でＹ_３＋＝（２／３）・λ／２ｄ_ｙが選択される。Ｙ_３＋が２Ｙ／Ｌに等しく、それにより、使用したレーザと使用した検出装置の機能として空間フィルタリング装置の特性を応用できる。

図６は、したがって、２次元フーリエ変換の後に、本発明によって波数ベクトル空間で得られた複数の領域のマッピングを示す。

このマッピングは、範囲［−Ｙ／Ｌ，Ｙ／Ｌ］の中央列または領域２と、左列または領域１と、右列または領域３と、に分解される。領域２および領域４の内部は、基準波とそれ自身との干渉項（Ｅ_Ｒ・Ｅ_Ｒ ^＊）を示す。１次のオーダーで音響信号がない場合は、領域２および４の外側の画像で観測されるノイズは、基準波に伴う“ショットノイズ”に対応する。音響信号が存在する場合は、領域３は波数ベクトル空間の領域を示し、本発明に係る有用な信号が検出される。入射角θ_Ｙの符号によっては、この領域３は、当然ながら図６中左に位置する。

音響成分の複素振幅Ｅ_Ｓの計算を可能にする有効な情報は、基準波と音響光学スペックルとの干渉（Ｅ_Ｒ・Ｅ_Ｓ ^＊）に対応する。撮像物の測定点（Ｕ，Ｖ，Ｗ）に関連するデジタル情報は、例えば、領域３の画素にわたって計算された強度（｜ＴＦ Ｅ_Ｓ｜^２）を合計して算出される。図６の縦軸に関して領域３と対称な領域１の画素にわたる強度の総和は、制御として機能する。また、領域３の画素の強度の総和から、領域１の画素の強度の総和を引くことができ、これにより、どんな系統的な測定エラーも減少できる。領域２にわたって測定された強度の総和もまた、情報を含んでいる。なぜならば、この総和は、スペックルの無相関性に伴うノイズの度合を特徴付けるからである。

また、従来の画像処理の手段を用いることも可能である。すなわち、座標（Ｕ，Ｖ，Ｗ）の測定点で得られた信号の制御マッピングを表示するように（例えば、様々な領域の位置を確かめる）、４相復調、フーリエ変換、および強度の計算後に、例えば８×８、又は１６×１６、又は３２×３２の画素の強度を平均化する。

また、検出装置の各行に対してのみ、１次元フーリエ変換（ｙに沿って）を実行してもよく、この場合、信号ＴＦ_１ Ｅ_Ｓ（Ｕ，Ｖ，Ｗ，ｋ）は、図６に示されたものにはならず、領域４は、この図においてＸ−Ｙ全域に広がったものとなる。残る処理は正しいままである。

したがって、４相復調の場合、撮像物ＯＢＪの座標（Ｕ，Ｖ，Ｗ）の測定点に関連する、求めていたデジタル情報を得る。“周波数チャーピング”等に対し、このようにして得られたデジタル情報は、使用された周波数ｆ_Ａに関連しており、そして、上述の手段によって、［Ｕ−Ｄｘ；Ｕ＋Ｄｘ］間の複数の点と関連した座標およびデジタル情報を得る。

角度θ_Ｙの選択により、得られた信号を適切に分けることができるが、与えられた角度θ_Ｙの使用を可能にする装置の適切な配置には、制御工程が必要である。この制御工程においては、例えば、音響波を抑制し、周波数ｆ_Ｉで散乱した場の成分を検出できるように基準波の周波数ｆ_Ｒを調整することにより、画像を取得する。例えば、ｆ_Ｒ＝ｆ_Ｉ＋ｆ_Ｃ／４を選択し、４相復調を実行する。絞り４を通して見られる領域の鮮明な画像が、検出装置の平面において検出される信号のフーリエ変換によって得られるように、正確にレンズが配置される。

先に記載したように、フーリエ変換の計算において、計算行列の端が“エーリアシング”条件に対応する。この制御工程は、例えば、有効な信号領域の外側端が計算マトリックス端にそれ程近くなく、かつ／または、ノイズがある領域と有効な信号を有する領域が、この領域の内側端のレベルで接触しているが重複していないようにすることができる。

図７は、本発明の第２実施形態に係る方法を示しており、ここでは、空間フィルタリング装置は使用されていない。特に、低強度のレーザに対しては、無相関なスペックルの項が大きすぎない場合は、“ショットノイズ”が主要なノイズになる。第１実施形態の実現（図５）が、そのような場合であり、領域１と２にわたって測定した合計値が、ほぼ同じである場合である。また、この構成は、音響波の出力が十分であり、十分な音響光変換を達成する他の形態に対しても有効である。その結果、音響成分に伴う信号は、無相関なスペックルのノイズより大きくなる。

空間フィルタリング装置を省くことが可能である。なぜなら、通常スペックルと通常スペックルとの干渉項（Ｅ_Ｉ・Ｅ_Ｉ ^＊）が、有効な項（Ｅ_Ｒ・Ｅ_Ｓ ^＊）から空間的に分離できなくとも、この項は、基準波による“ショットノイズ”より十分に小さく在るからである。また、この場合、例えば半反射鏡によって、ゼロの角度のθ_Ｙを選ぶことが可能である。図６の領域４の数画素上に集中する、基準波とそれ自身との干渉項（Ｅ_Ｒ・Ｅ_Ｒ ^＊）のみが、フィルタリングするノイズ項を含む。このフィルタリングは、単純に、（０，０）周辺の波数ベクトル空間の成分をデジタル処理で取り除き、信号の強度の画素にわたる総和から領域４を排除することによって実行できる。このノイズ項を省略してもよく、これは、波数ベクトル空間よりも、実空間の画素にわたって強度の総和を計算できる（そして、この場合は、フーリエ変換を実行する必要も、もはやない）。

図８は、本発明の、特に図７の第２の例に係る方法を実現するために使用する他の生成装置を示す。レーザＬＡＳから放射した波長λの放射波ＥＭＩは、振幅変調装置ＭＡで受けられる。この装置は、周波数ｆ_Ｉの搬送波ＰＯＲと、周波数ｆ_Ｒの２つの振幅変調側波帯ＬＯＴＭＯＤおよびＬＯＴＭＯＤ’を生成する。これらの３つの波は、半透明な鏡板に照射され、この鏡板は、照射されたこれらの波のそれぞれの一部を、撮像物ＯＢＪに伝送する。搬送波の伝送された部分が、入射波ＩＮＣを構成する。さらに、この半透明な鏡板は、これらの波のそれぞれの一部を、検出装置に向けて反射し、ＬＡＴＭＯＤの反射された部分が、基準波ＲＥＦを構成する。

この装置は、空間フィルタリング装置が使用されない第２実施形態でのみ適用される。なぜならば、強度変調側波帯の場合、空間フィルタリングは、各種の干渉項を分離できないからである。

図９は、本発明の第３実施形態に係る空間フィルタリング装置の詳細を示す。多数の特徴を、第１実施形態を記載した時にすでに示したので、ここでは再度記載しない。この第３実施形態においては、基準波ＲＥＦは平面波ではなく周波数ｆ_Ｒの球面波である。このような基準球面波は、例えば、上述の生成装置ＧＥＮによって生成された基準平面波を基にし、この基準平面波をレンズ５によって、例えば絞り２の面に置かれた小さな鏡４に集光させることによって得られる。この基準波ＲＥＦは、さらに、鏡４に角度θ_Ｙで入射する。このとき、この基準波ＲＥＦは、反射されて信号波ＳＩＧと角度θ_Ｙを形成しながら、検出装置ＤＥＴに到達するように、入射される。この実施例においては、検出装置ＤＥＴに向かう信号波を平行にする集束レンズを使用する必要がない。具体的には、基準波の発散的な性質が同様な役割を果たす。

さらに、絞り２を撮像物ＯＢＪのサンプルホルダから取り外すことができ、また、この絞り２を、タンクと検出装置ＤＥＴとの間に配置するか、または、タンク１の前面に固定することができる。音響周波数ｆ_ＡのトランスデューサＴＲＡＮＳと撮像物ＯＢＪとの音響結合が、水１で満たされたタンクによって達成され、そしてこの後者が、基準波ＲＥＦが伝送する光路上に、鏡４に到達する前に置かれる場合には、この装置は有益である（また、絞りが使用される他の実施形態にも前述した記載が有効である）。

測定された干渉画像の処理は、上記の処理と同一である。ここで、Ｌは絞り２と検出装置ＤＥＴとの間の距離を示す。

図１０は、鏡４を備えた絞り２を示しており、検出装置ＤＥＴから視覚化したものである。絞り２は、例えば中央に、幅Ｙで高さＸのスリットを示す。図６に示すような画像を得るために、例えば円形の鏡４が、スリットの高さの中央の上部に配置される。さらに、鏡４は、絞り２のスリットに対して横にずれる。この横方向へのシフトは、処理装置ＣＡＬＣによって計算が実行されると、その効果が要求される図６の領域２および３のシフトに関連する。例えば、検出装置から見られる反射された基準波の湧き出し点と鏡に最も近い絞り２の縁との間の横方向のオフセット角をθ_Ｙに等しく選択することによって、重なり合わずに並列した図６の領域２および３が得られる。この場合、鏡が絞り２のスリットに侵入しないように、鏡４のサイズを選択するように注意が払われねばならない。したがって、図６の領域２と３とを分けるのは容易である。

本実施形態は、検出装置ＤＥＴの最適な位置を見出すことが、さらに容易にでき、これにより、撮像物ＯＢＪの調査された点（Ｕ，Ｖ，Ｗ）に対する良品質な情報を得られる。具体的には、Ｌを増大することによって、すなわち、絞り２から検出装置ＤＥＴを離すことによって、より良く“エーリアシング”条件を満たすように、信号波ＳＩＧと基準波ＲＥＦとの角度が減少する。一方、信号の強度は、この距離Ｌで減少する。検出の有効性と、すでに述べた図６の領域３のサイズと、の間でなされる妥協は、単純に、絞り２に対する検出装置ＤＥＴの位置によって影響される。

本実施形態において、検出装置の最適な位置付けは、図６の領域２および３の相対的な位置とは無関係になされ、自身が基準波の集光する点を適切に位置付けることによって実現される（絞り２上の鏡４の中心）。

図１〜１０に示された実施形態のいずれか１つに係る、上記の方法を実現することによって、撮像物の座標（Ｕ，Ｖ，Ｗ）の測定点に関連するデジタル情報（Ｎ相複調の場合）や、撮像物上の間隔（［Ｕ−Ｄｘ，Ｕ＋Ｄｘ］，Ｖ，Ｗ）の各点に関連するデジタル情報（“周波数チャーピング”や他の同様処理の場合）、すなわち、点状もしくは１次元の情報を得る。

２次元や３次元の被写体の画像を得るためには、他の方向（Ｎ相復調に対するｘやｙやｚ）（“周波数チャーピング”に対するｙやｚ）に沿った音響波の焦点の位置を置き換えた後に、上述した一連の動作全部を繰り返す必要がある。この焦点の位置の置き換えは、撮像物とトランスデューサとの音響結合を保持したまま、トランスデューサを置き換えたり（方向ｘやｙやｚに沿って）、また、例えば、異なる焦点を有する複数のトランスデューサを用いたりすることによって達成される。

空間フィルタリング装置ＣＯＬ（これを用いる場合）、検出装置ＤＥＴ、基準波を検出装置上で干渉させるための装置（場合によって）や、焦点の位置すなわち測定点の座標（Ｕ，Ｖ，Ｗ）の移動と共同して生成装置ＧＥＮの全てまたは一部を移動させることは有効である。あるいは、装置の残り全てを同様に動かさずにすむように、撮像物のみを反対方向に移すことが好ましい。

したがって、これら２つの装置の結合は、超音波映像の技法および本発明に基づく装置による光学情報によって、音響トランスデューサによって供給された純粋な音響情報を撮像物から得ることができる。

本発明に基づく方法の第１実施形態を示す図である。 本発明に基づく２つのコヒーレントな波を生成する装置の一例の詳細を示す図である。 ２つのコヒーレントな波の間の位相差を、時間の関数として、一般的な形式で示した図である。 ２つのコヒーレントな波の間の位相差を、時間の関数として、特定の場合で示した図である。 図２の生成装置を備える本発明に基づく第１実施形態を示す図である。 取得した信号のマッピングを示す図である。 本発明に基づく方法の第２実施形態を示す図である。 本発明に基づく２つのコヒーレントな波を生成する装置の他の例の詳細を示す図である。 本発明の第３実施形態の詳細を示す図である。 第３実施形態で用いられる絞りの後面を示す図である。

Claims (27)

1. 撮像物（ＯＢＪ）の音響光学撮像方法であって、
   （ａ）周波数ｆ_Ｉの入射光波（ＩＮＣ）と周波数ｆ_Ｒの基準光波（ＲＥＦ）とを生成する工程であって、当該基準光波は、入射光波（ＩＮＣ）とコヒーレントであるとともに、既知の位相差φ_ｉ（ｔ）を示す工程と、
   （ｂ）振動生成装置（ＴＲＡＮＳ）によって、第１の撮像物方向（ｘ_０）に沿って、音響周波数ｆ_Ａで、撮像物（ＯＢＪ）の領域（Ｄｘ，Ｄｙ，Ｄｚ）を振動させる工程と、
   （ｃ）前記入射光波（ＩＮＣ）を撮像物（ＯＢＪ）に照射して、散乱信号波（ＤＩＦ）を生成する工程と、
   （ｄ）散乱信号波の少なくとも一部を検出装置（ＤＥＴ）に照射する工程と、
   （ｅ）基準光波（ＲＥＦ）を撮像物（ＯＢＪ）に透過させずに前記検出装置（ＤＥＴ）に照射し、前記検出装置（ＤＥＴ）の点ｒに、時間ｔで変化する干渉画像Ｉ（ｒ，ｔ）を生成する工程と、
   （ｆ）干渉画像Ｉ（ｒ，ｔ）から撮像物（ＯＢＪ）の測定点に関するデジタル情報を抽出する工程、
   （ｇ）音響波の第１周波数ｆ_Ａと基準波（ＲＥＦ）の第１周波数ｆ_Ｒとで行なわれる工程（ａ）から工程（ｆ）の第１サイクルと、音響波の第１周波数ｆ_Ａでコード化された音響波の第２周波数ｆ’_Ａと基準波（ＲＥＦ）の第１周波数ｆ_Ｒでコード化された基準光波（ＲＥＦ）の第２周波数ｆ’_Ｒとで繰り返し行われる工程（ａ）から工程（ｆ）の第２サイクルと、を実行する工程と、
   （ｈ）デジタル情報に関連する音響周波数ｆ_Ａ，ｆ_Ａ’で振動する撮像物（ＯＢＪ）の測定点の３次元の座標（Ｕ，Ｖ，Ｗ）を得る工程であって、当該測定点の３次元の座標（Ｕ，Ｖ，Ｗ）は、使用された周波数の関数として各サイクルの工程（ｆ）において得られた撮像物（ＯＢＪ）の測定点に関するデジタル情報をデコードすることによって得られる工程と、から構成されることを特徴とする撮像物（ＯＢＪ）の音響光学撮像方法。
2. 工程（ｆ）では、前記検出装置（ＤＥＴ）に照射された散乱信号波の一部の音響成分が検出され、
   当該音響成分は、入射光波（ＩＮＣ）の周波数ｆ_Ｉと音響周波数ｆ_Ａの倍音との合計に対応する周波数にあることを特徴とする請求項１に記載の音響光学撮像方法。
3. 工程（ａ）では、前記基準光波（ＲＥＦ）は、入射光波（ＩＮＣ）の周波数ｆ_Ｉと音響周波数ｆ_Ａの倍音との合計に等しい周波数ｆ_Ｒで生成されることを特徴とする請求項２に記載の音響光学撮像方法。
4. 工程（ｂ）では、撮像物（ＯＢＪ）に位置する焦点に集束した音響波が生成され、
   工程（ｈ）では、測定点の３次元の座標（Ｕ，Ｖ，Ｗ）が、前記焦点の座標として得られることを特徴とする請求項３に記載の音響光学撮像方法。
5. 工程（ａ）から工程（ｈ）は、撮像物における音響波の複数の焦点に対して繰り返し行なわれ、
   これら複数の焦点は、第１の撮像物方向（ｘ_０）に沿って配列されることを特徴とする請求項４に記載の音響光学撮像方法。
6. 第１の撮像物方向（ｘ _０ ）に沿い、座標（Ｕ，Ｖ，Ｗ）を持つ点の周囲に広がる座標（［Ｕ−Ｄｘ，Ｕ＋Ｄｘ］，Ｖ，Ｗ）の区間に集束する音響波の周波数のスキャンが行なわれ、
   ｆ_Ｒがｆ_Ｉ±Ｈ・ｆ_Ａ、（Ｈは０でない整数）に等しく保ちながら、基準光波の周波数ｆ_Ｒでスキャンを共同で行ない、
   広がった間隔の点集合（［Ｕ−Ｄｘ，Ｕ＋Ｄｘ］，Ｖ，Ｗ）に対応する干渉画像Ｉ（ｆ_Ａ，Ｖ，Ｗ，ｒ）が、各画素ｒおよび各周波数ｆ_Ａに対して記録され、
   干渉画像Ｉ（ｆ_Ａ，Ｖ，Ｗ，ｒ）の周波数から時間への１次元フーリエ変換が、周波数ｆ_Ａにしたがって、各画素ｒに対して行われ、
   座標（Ｕ’，Ｖ，Ｗ，）の少なくとも１つの測定点に対応する少なくとも１つの干渉画像Ｉ（ｒ）が、フーリエ変換後に得られる時間と、第１の撮像物方向に沿った長さＵ’（Ｕ’はほぼＵに等しい）とを、撮像物における音響波の伝播速度を用いて置き換えることによって得られる操作が実行されることを特徴とする請求項１に記載の方法。
7. 少なくとも工程（ａ）から工程（ｈ）を、撮像物（ＯＢＪ）の第１の方向（ｘ_０）と平行にはならないように、撮像物（ＯＢＪ）に対する前記振動生成装置の位置を変えた後に繰り返すことを特徴とする請求項１から６のいずれかに記載の音響光学撮像方法。
8. 工程（ｆ）では、干渉画像Ｉ（ｒ，ｔ）をもとに、デジタル情報を得るために、音響成分の複素振幅Ｅ_Ｓ（ｒ）を推定することを特徴とする請求項１から７のいずれかに記載の音響光学撮像方法。
9. 前記検出装置（ＤＥＴ）は単一画素検出装置であり、
   工程（ｆ）では、対象物によって散乱された複素振幅Ｅ_Ｓ（ｒ）の場の強度として、デジタル情報が得られることを特徴とする請求項８に記載の音響光学撮像方法。
10. 前記検出装置は多画素検出装置であり、
    工程（ｆ）では、前記検出装置の画素ｒの、対象物によって散乱された複素振幅Ｅ_Ｓ（ｒ）の場の強度の少なくとも一部にわたる合計として、デジタル情報が抽出されることを特徴とする請求項８に記載の音響光学撮像方法。
11. 工程（ｄ）では、前記検出装置の各画素によって見られる散乱信号波（ＤＩＦ）の一部の角度範囲を、少なくとも１つの方向に沿って制限するように、空間フィルタリング装置（ＣＯＬ）が使用されることを特徴とする請求項８に記載の音響光学撮像方法。
12. 第１絞り方向に沿った長さがＸで、前記第１絞り方向に直交する第３絞り方向に沿った長さがＹの矩形状の絞りを備える空間フィルタリング装置（ＣＯＬ）と、撮像物（ＯＢＪ）の入射光波（ＩＮＣ）の伝播方向の下流部に位置する焦点を持つ焦点距離がＬのレンズと、が用いられ、
    前記検出装置の各画素によって見られる散乱信号波（ＤＩＦ）の一部の角度範囲を制限し、前記検出装置（ＤＥＴ）に入射した基準光波（ＲＥＦ）が全て平面波であることを特徴とする請求項１１に記載の音響光学撮像方法。
13. 撮像物（ＯＢＪ）と前記検出装置（ＤＥＴ）との間に、後者から距離Ｌを離して配置され、第１絞り方向に沿った長さがＸで、前記第１絞り方向に直交する第３絞り方向に沿った長さがＹの矩形状の絞りを備える空間フィルタリング装置（ＣＯＬ）が用いられ、
    前記検出装置の各画素によって見られる散乱信号波の一部の角度範囲を制限し、
    前記検出装置（ＤＥＴ）に照射される基準光波（ＲＥＦ）は、絞りの平面に位置する原点から発せられる球面波であることを特徴とする請求項１１に記載の音響光学撮像方法。
14. 基準光波（ＲＥＦ）と散乱信号波（ＤＩＦ）とは、ゼロでない角度θ_Ｙをなしながら前記検出装置上で干渉し、
    角度θ_Ｙは前記検出装置（ＤＥＴ）上の２つの波の入射面で測定されることを特徴とする請求項１２又は１３のいずれかに記載の音響光学撮像方法。
15. 使用する前記検出装置は多画素検出装置であり、
    前記検出装置の平面で空間的に急激に変化する複素振幅Ｅ_Ｓ（ｒ）の音響成分の一部は、分離されることを特徴とする請求項１２から１４のいずれかに記載の音響光学撮像方法。
16. 前記検出装置（ＤＥＴ）は、第１検出装置方向（ｘＤ）に沿った行および第３検出装置方向（ｙＤ）に沿った列から成るマトリックスとして配列された画素を備え、
    工程（ｆ）では、
    （ｆ１）少なくとも１つの行または列に対する１次元フーリエ変換が、前記検出装置（ＤＥＴ）の行（ｌ）または列（ｃ）に沿って、波動ベクトル空間へ、複素振幅Ｅ_Ｓ（ｒ）に対して実行され、この行または列に対して、対象物によって散乱された複素振幅ＥＳ（ｒ）の場として、場ＴＦ_１ Ｅ_Ｓ（ｋ）が得られ、
    （ｆ２）いくつかの総和領域が、波数ベクトル空間において定義され、
    （ｆ３）少なくとも１つの領域の各点ｋでの場ＴＦ_１ Ｅ_Ｓ（ｋ）の強度が、この領域中で総和され、
    （ｆ４）デジタル情報は、各領域でこのように得られた総和の線形結合として抽出されることを特徴とする請求項１２から１５のいずれかに記載の音響光学撮像方法。
17. 前記検出装置（ＤＥＴ）は、第１検出装置方向（ｘＤ）に沿った行および第３検出装置方向（ｙＤ）に沿った列から成るマトリックスとして配列された画素を備え、
    工程（ｆ）では、
    （ｆ１）前記検出装置（ＤＥＴ）の平面から波数ベクトル空間への、複素振幅Ｅ_Ｓ（ｒ）の２次元フーリエ変換が行われ、対象物によって散乱された複素振幅ＥＳ（ｒ）の場として、場ＴＦ_２ Ｅ_Ｓ（ｋ）が得られ、
    （ｆ２）いくつかの総和領域が、波数ベクトル空間において定義され、
    （ｆ３）少なくとも１つの領域の各点ｋでの場ＴＦ_２ Ｅ_Ｓ（ｋ）の強度が、この領域中で総和され、
    （ｆ４）デジタル情報は、各領域でこのように得られた総和の線形結合として抽出されることを特徴とする請求項１２から１５のいずれかに記載の音響光学撮像方法。
18. 基準波（ＲＥＦ）と散乱信号波（ＤＩＦ）との間の角度である角度θ_Ｙは３Ｙ／２Ｌに等しく、
    工程（ｆ２）では、いわゆる中央領域の第１総和領域、いわゆる左の領域の第２総和領域、いわゆる右の領域の第３総和領域が定義され、
    工程（ｆ４）では、デジタル情報は、左の領域の総和と右の領域の総和との線形結合として抽出されることを特徴とする請求項１６又は１７のいずれかに記載の音響光学撮像方法。
19. 工程（ａ）では、
    波長λのレーザー光源は、周波数ｆ_Ｌの放射波（ＥＭＩ）を発生し、
    放射波（ＥＭＩ）の振幅変調手段（ＭＡ）は、入射周波数ｆ_Ｉの搬送波（ＰＯＲ）と、周波数ｆ_Ｒの少なくとも１つの波に対応する振幅変調された側帯波（ＬＡＴＭＯＤ）と、を生成し、
    半反射鏡は、前記側帯波（ＬＡＴＭＯＤ）の一部と入射光波（ＩＮＣ）を形成する搬送波（ＰＯＲ）の一部とを伝送し、搬送波（ＰＯＲ）の一部と基準光波（ＲＥＦ）を形成する前記側帯波（ＬＡＴＭＯＤ）の一部とを反射することを特徴とする請求項１から１０のいずれかに記載の音響光学撮像方法。
20. 工程（ａ）では、
    波長λのレーザー光源は、周波数ｆ_Ｌの放射波（ＥＭＩ）を発生し、
    第１の音響光変調器（ＭＡＯ１）は、撮像物（ＯＢＪ）に対する入射光波（ＩＮＣ）を形成するために当該放射波（ＥＭＩ）の一部を伝送し、さらに、周波数が、当該放射波に対して、第１の周波数δｆ_１だけシフトされた波（ＤＥＣ）を生成し、
    第２の音響光変調器（ＭＡＯ２）は、第１の周波数だけシフトされた波（ＤＥＣ）を受けて、周波数が、第２の周波数だけシフトされた波を生成し、当該シフト波の周波数は、第１の周波数だけシフトされた波（ＤＥＣ）に対して、第２の周波数δｆ_２だけシフトされたものであり、基準光波（ＲＥＦ）を形成する第２の周波数だけシフトされた波の周波数は、入射光波（ＩＮＣ）に対し周波数δｆ＝δｆ_１＋δｆ_２だけシフトしたものであり、このようにして、これら２つの波間の既知の位相差φ_ｉ（ｔ）を決定することを特徴とする請求項１から１８のいずれかに記載の音響光学撮像方法。
21. 工程（ａ）では、
    電子的スレービングによって位相が固定された２つの独立なレーザー光源は、既知の位相差φ_ｉ（ｔ）を示す入射光波（ＩＮＣ）と基準光波（ＲＥＦ）を生成することを特徴とする請求項１から１８のいずれかに記載の音響光学撮像方法。
22. 工程（ａ）では、
    波長λのレーザー光源は、周波数ｆ_Ｌの放射波（ＥＭＩ）を発生し、
    半反射鏡は、撮像物（ＯＢＪ）に対する入射光波（ＩＮＣ）を形成するために当該放射波（ＥＭＩ）の一部を伝送し、当該放射波（ＥＭＩ）の第２の部分を伝送し、
    第１の音響光変調器（ＭＡＯ１）は、当該放射波の第２の部分を受けて、周波数が、当該放射波に対して、第１の周波数δｆ_１だけシフトした波（ＤＥＣ）を生成し、
    第２の音響光変調器（ＭＡＯ２）は、第１の周波数だけシフトされた波（ＤＥＣ）を受けて、シフト波（ＤＥＣ）に対して、周波数が、第２の周波数δｆ_２だけシフトされた波を生成し、基準光波（ＲＥＦ）を形成する第２の周波数だけシフトされた波の周波数は、入射光波（ＩＮＣ）に対し周波数δｆ＝δｆ_１＋δｆ_２だけシフトされたものであり、このようにして、これら２つの波間の既知の位相差φ_ｉ（ｔ）を決定することを特徴とする請求項１から１８のいずれかに記載の音響光学撮像方法。
23. 撮像物（ＯＢＪ）は、生物学的組織であることを特徴とすることを特徴とする請求項１から２２のいずれかに記載の音響光撮像方法。
24. 前記振動生成装置（ＴＲＡＮＳ）は、撮像物（ＯＢＪ）の領域（Ｄｘ，Ｄｙ，Ｄｚ）の音響キューを得るために用いられ、
    工程（ｆ）で抽出されたデジタル情報は、前記音響キューと共同して使用されることを特徴とする請求項１から２３のいずれかに記載の音響光撮像方法。
25. 撮像物（ＯＢＪ）に対する音響光学撮像装置であって、
    周波数ｆ_Ｉの入射光波（ＩＮＣ）、および、この入射光波（ＩＮＣ）とコヒーレントでかつ既知の位相差φ_ｉ（ｔ）を示す周波数ｆ_Ｒの基準光波（ＲＥＦ）を生成する手段と、
    第１の撮像物方向（ｘ_０）に沿って、音響波の第１周波数ｆ_Ａおよび音響波の第２周波数ｆ’_Ａで撮像物（ＯＢＪ）の領域（Ｄｘ，Ｄｙ，Ｄｚ）を振動する振動生成装置と、
    撮像物（ＯＢＪ）に入射光波（ＩＮＣ）を照射し、これにより散乱信号波（ＤＩＦ）を生成する手段と、
    検出装置（ＤＥＴ）と、
    前記散乱信号波（ＤＩＦ）の少なくとも一部（ＳＩＧ）を前記検出装置（ＤＥＴ）に照射する手段と、
    基準光波（ＲＥＦ）を撮像物（ＯＢＪ）に透過させずに前記検出装置（ＤＥＴ）に照射し、前記検出装置（ＤＥＴ）の点ｒに時間ｔで変化する干渉画像Ｉ（ｒ，ｔ）を生成する手段と、
    干渉画像Ｉ（ｒ，ｔ）から撮像物（ＯＢＪ）の測定点に関するデジタル情報を抽出するとともに、当該デジタル情報に関連する音響周波数ｆ_Ａ，ｆ_Ａ’で振動する撮像物の測定点の３次元の座標（Ｕ，Ｖ，Ｗ）を、前記振動生成装置により使用された周波数ｆ_Ａ，ｆ’_Ａの関数として撮像物（ＯＢＪ）の測定点に関するデジタル情報から抽出する手段（ＣＡＬＣ）と、を備えることを特徴する音響光撮像装置。
26. 前記撮像物の測定点に関連するデジタル情報を可視化する手段と、
    当該撮像物（ＯＢＪ）をずらす手段と、を備えることを特徴とする請求項２５記載の音響光撮像装置。
27. さらに、撮像物（ＯＢＪ）の入射光波（ＩＮＣ）の伝播方向の下流部に位置する空間フィルタリング装置を備えることを特徴とする請求項２５または２６に記載の音響光撮像装置。

Images