複数の剪断波を用いる撮像方法および装置
本発明の撮像方法は、剪断成分と圧縮成分とを含む力学波を粘弾性媒体内で発生させ、粘弾性媒体の運動パラメータを力学波の伝搬中にさまざまな点で求めることにある。本方法は、力学波の圧縮成分による誤差を取り除くように運動パラメータが処理される修正段階を含む。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、複数の剪断波を用いる撮像方法および装置に関する。
【０００２】
特に、本発明は、剪断成分および圧縮成分を含む力学波（onde mecanique） が粘弾性媒体内で発生され、かつ、この力学波の伝搬が続いている間における、 様々な点での粘弾性媒体の運動を記述する少なくとも一つのパラメータを求める 撮像方法に関する。
【０００３】
したがって、特に、粘弾性媒体の他の部分と硬度が異なる複数の領域、また は、粘弾性媒体の他の部分と力学緩和時間が異なる複数の領域を特定するために 、定性的および／または定量的な分析を実施することができる。
【背景技術】
【０００４】
ＰＣＴ特許出願公開明細書第０４／２１０３８号は、そのような方法の例を 示しており、対象とする力学波、本方法では剪断波、を生成する局所的なパルス 状の推力の離れた位置での発生に続く、平面内の運動（特に変位）を記述してい るパラメータが求められる。
【発明の開示】
【発明が解決しようとする課題】
【０００５】
この方法はすでに完全に満足できるが、本発明の目的は、その精度を改良し 、必要であれば、３次元の撮像に使用できるように、さらにこれを完全にするこ とである。
【課題を解決するための手段】
【０００６】
この目的のために、対象となる方法は、運動パラメータが、上記の力学波の 圧縮成分による影響を取り除くように処理される補正ステップを含むことを特徴 とする。
【０００７】
したがって、パルス状の推力が使用されるときでも、それによって発生する 弾性力学波の圧縮成分が運動パラメータの計測を妨害する場合、剪断成分の効果 だけが考慮される。例えば３次元撮像などの比較的長い計測を実施するために、 推力がなくても持続する力学波が使用される場合、力学波の圧縮成分は、運動パ ラメータの計測中も依然として存在し、その結果、この計測を大きく妨害する。
本発明による処理は、この妨害を取り除くことができ、必要であれば、粘弾性媒 体の複数の粘弾性パラメータを２または３次元にマッピングすることができる。
【０００８】
本発明による方法の好ましい実施形態では、以下の複数の構成の一つまたは ２つ以上を、さらに任意で用いてもよい。
【０００９】
−上記の運動パラメータの回転が補正ステップ中に計算される。
【００１０】
−粘弾性媒体内の複数の点における力学波の伝搬が、エコーグラフで観測さ れる少なくとも１つの観測ステップ（ｂ）を含む。
【００１１】
−観測ステップ（ｂ）は、以下のサブステップを含む。
【００１２】
複数の圧縮超音波（ondes ultrasonores de compression）の連続するショ ットが、複数の変換器の列によって、（焦点が合わされているかどうかにかかわ らず）少なくとも１秒あたり１０ショットの割合で、粘弾性媒体内に伝送される サブステップ（ｂ１）。
【００１３】
粘弾性媒体に相互作用する複数の圧縮超音波によって発生した複数のエコー を含む粘弾性媒体から受信した複数の音響信号がリアルタイムで検出および記録 されるサブステップ（ｂ２）。
【００１４】
−以下のステップが実施される少なくとも１つの処理ステップ（ｃ）をさら に含む。
【００１５】
力学波の連続する複数の伝搬画像を求めるために、サブステップ（ｂ２）の 間に粘弾性媒体から受信した複数の連続する音響信号が処理されるステップ（ｃ １）。
【００１６】
連続する複数の伝搬画像から粘弾性媒体の運動パラメータが様々な点で求め られるステップ（ｃ２）。
【００１７】
−粘弾性媒体内での力学波の連続する複数の伝搬画像は、いくつかの視点に 沿って生成され、これらの連続する画像は、運動パラメータのベクトル値を求め るために組み合わされる。
【００１８】
−運動パラメータの３次元ベクトル値が求められる。
【００１９】
−いくつかの視点に沿った力学波の連続する複数の伝搬画像を生成するため に、解析平面内の伝送／検出方向に沿って複数の圧縮超音波を伝送し検出する複
数の超音波変換器の列（複数の圧電変換器の直線状の列、または位相列、１．５
Ｄ列、１．７５Ｄ列、または２Ｄ配列などの任意の他の列）が使用され、直線状
の列は、粘弾性媒体の外側で変位することで、解析平面を回転させる（これらの
構成は、前述の構成とは独立して任意に使用することが可能で、そのため、本発
明の主題は、これらの構成自体でもある）。
【００２０】
−列は、解析平面を回転しない並進運動によって変位し、その後、列は、回 転することで、解析平面を回転させる。
【００２１】
−複数の変換器の列の伝送／検出方向は、解析平面を並進も回転もさせずに 変化し、その後、列は、並進および回転することで、解析平面を回転させる。
【００２２】
−運動パラメータは粘弾性媒体の変位、変位速度および応力から選択される 。
【００２３】
−粘弾性媒体内の伝搬パラメータのマップを求めるために、粘弾性媒体内の 少なくとも特定の複数の点での、力学波の剪断成分の伝搬を記述する少なくとも １つのパラメータが、運動パラメータの時間の経過に伴う変化から計算されるマ ッピングステップ（ｄ）をさらに含む。
【００２４】
−マッピングステップ中に計算される剪断波の伝搬パラメータは、複数の剪 断波の速度、剪断弾性率、ヤング率、複数の剪断波の減衰、剪断弾性、剪断粘性 、および力学緩和時間、または、異方性媒体の場合これらのパラメータうちのい くつかのパラメータから選択される。
【００２５】
また、本発明の主題は、剪断成分と圧縮成分とを含む力学波を粘弾性媒体内 で発生させる手段と、力学波の伝搬中に様々な点で粘弾性媒体の運動を記述する 少なくとも１つのパラメータを求める手段とを含み、力学波の圧縮成分による影 響を取り除くように運動パラメータを処理する補正手段を含むことを特徴とする 撮像装置である。
【００２６】
本発明の装置の好ましい実施形態では、以下の複数の構成の１つまたは２つ 以上をさらに任意で用いてもよい。
【００２７】
−補正手段は、運動パラメータの回転を計算するように構成されている。
【００２８】
−力学波の伝搬を粘弾性媒体内の複数の点でエコグラフィーによって観測す る複数の超音波変換器の制御されている列を含む。
【００２９】
−装置は、以下の複数の手段を含む。
【００３０】
複数の変換器の列によって、少なくとも１秒あたり１０ショットの割合で（ 焦点が合わされているかどうかにかかわらず）複数の圧縮超音波の連続するショ ットを粘弾性媒体内に伝送する手段（ｂ１）。
【００３１】
粘弾性媒体の複数の反射粒子に相互作用する複数の焦点が合っていない圧縮 超音波によって発生した複数のエコーを含む、粘弾性媒体から受信した複数の音 響信号をリアルタイムで検出および記録する手段（ｂ２）。
【００３２】
−装置は、以下のために構成されている複数の処理手段を含む。
【００３３】
力学波の連続する複数の伝搬画像を求めるために、粘弾性媒体から受信した 複数の連続する音響信号を処理する（ｃ１）。
【００３４】
連続する複数の伝搬画像から粘弾性媒体の運動パラメータを観測フィールド 内の様々な点で求める（ｃ２）。
【００３５】
−複数の処理手段は、粘弾性媒体内でいくつかの視点に沿って力学波のいく つかの連続する伝搬画像を生成し、連続する画像を運動パラメータのベクトル値 を求めるために組み合わせるように構成されている。
【００３６】
−複数の処理手段は運動パラメータの３次元ベクトル値を求めるように構成 されている。
【００３７】
−複数の超音波変換器の列はいくつかの視点に沿った力学波の連続する複数 の伝搬画像を生成するために、解析平面内の伝送／検出方向に沿って複数の圧縮 超音波を伝送して検出し、列（複数の圧電変換器の直線状の列、または位相列、 １．５Ｄ列、１．７５Ｄ列、または２Ｄ配列などの任意の他の列）を粘弾性媒体 の外側で変位させる複数の変位手段を含み、複数の手段は解析平面を回転させる ように構成されている（これらの構成は、前述の構成とは独立して任意に使用す ることが可能で、そのため、本発明の主題は、これらの構成自体でもある）。
【００３８】
−装置は、複数の変位手段を制御するように構成されている制御手段を含み 、制御手段は、解析平面を回転させずに列を並進運動によって変位させ、その後 、列を回転させることで、解析平面を回転させる。
【００３９】
−装置は、複数の変位手段を制御するように構成されている制御手段を含み 、制御手段は、解析平面を並進も回転もさせずに、複数の変換器の列の伝送／検 出方向を変化させ、その後、列を並進および回転させることで、解析平面を並進 および回転させるように構成されている。
【００４０】
装置は、粘弾性媒体内の少なくとも特定の複数の点での力学波の剪断成分の 伝搬を記述する少なくとも１つのパラメータを計算するように構成されているマ ッピング手段を含む。
【００４１】
本発明の他の特徴および利点は、添付の図面と共に非限定的な例として示し ているその実施形態の１つの以下の説明を通して明らかになるであろう。
【発明を実施するための最良の形態】
【００４２】
さまざまな図において、同一の符号は、同一のまたは同様な要素を示してい る。
【００４３】
図１に示している撮像装置１は、圧縮超音波を散乱させる粘弾性媒体２内の 複数の弾性剪断波の伝搬を研究することを目的とし、粘弾性媒体２は、例えば、 以下のようなものである。
【００４４】
−不活性体、特に、産業用途のための品質管理の場合、または
−生体、例えば、医療用途における患者の体の一部
図に示している例では、撮像装置１は、特に、女性患者３の乳房の医療用撮 像のために設計され、そのため、乳房は、撮像される媒体２を構成している。
【００４５】
この目的のために、撮像装置１は、例えば、患者がうつぶせに横たわる検査 テーブル４を含んでいてもよい。このテーブルは、患者の少なくとも１つの乳房 ２を貫通させることができる開口５を備えており、貫通した乳房は、ジェルまた は液体、好ましくは水などの超音波カプリング剤の槽６に垂れ下がる。この液体 槽６は、容器７によって境界が定められおり、容器７は超音波変換器の列９と乳 房に接触している振動発生器８も含む。容器７は、特に振動発生器８が乳房２に 接触して適切に配置する、および、変換器列９が乳房に沿って適切に配置されて いることを確認するために透明であってもよい。
【００４６】
検査テーブル４は、必要に応じて、患者の両方の乳房が通過し、同じ液体槽 ６または２つの別個の槽に垂れ下がるように構成されていてもよい。つまり、撮 像装置は、患者の両方の乳房を同時に撮像するように任意に構成されていてもよ い。
【００４７】
振動発生器８は、複数の力学（弾性）波、特に圧縮成分と剪断成分とを含む 複数の持続波（ondes entretenues）を患者の乳房２の内部に伝送するために、 例えば、患者の乳房２の端部に接触するように配置される高さ調整プレート１１ 内で、低周波の力学振動を発生する音響変換器１０を含んでもよい。対象として いる複数の力学波は、例えば、０．１Ｈｚと１０ｋＨｚとの間の周波数を含んで もよい。
【００４８】
ただし、上記の複数の力学波は、例えば、以下のようなもので、前述の振動 発生器以外のもので発生させることができる。
【００４９】
−複数の低周波を発生し、乳房２に対して異なる方法で配置されている変換 器。
【００５０】
−例えば、上述の文献ＰＣＴ特許出願公開明細書第０４／２１０３８号、ま たは、米国特許出願第５ ６０６９７１号で説明されたように、放射圧力を使用 して時間変調されており焦点が合わされている複数の超音波によって複数の低周 波力学的信号を遠隔発生する超音波変換器の列。
【００５１】
図１に示している例では、プレート１１の高さは、乳房２に接触するように 調整可能である。この調整は、任意で、例えば、密封電気ステッピングモータな どで作動するねじシリンダ１３で行ってもよい。
【００５２】
シリンダ１３は、垂直軸Ｚに沿って、両方向矢印１４で示された両方向にプ レート１１および音響変換器１０を移動させるように、例えば、支持ロッド１２ を介して、プレート１１に接続されていてもよい。
【００５３】
図１および図２に示しているように、複数の変換器の列９は、例えば、直線 状列の形態でもよく、この場合、その複数の変換器は、垂直平面内に揃えられ、 複数の超音波を送受信するように、垂直解析平面（Ｄ−Ｚ）内の同一の送信／検 出方向Ｄに対して全て平行に構成される。
【００５４】
変換器列９は、例えば、固定プラットフォーム１６によって支持される回転 プラットフォーム１５に取り付けられていてもよい。２つのプラットフォーム１ ５、１６は水辺平面（Ｘ−Ｙ）内に配置され、患者の乳房２が中央通路１７を自 由に通過できるように開口していることが有利である。回転プラットフォーム１ ５は、例えば、プラットフォーム１５および１６の間に組み込まれている密閉電 気ステッピングモータ（図１および２では、図示しない）で動作する調整システ ムによって、固定プラットフォーム１６に対して垂直中央軸Ｚ０を中心に角度方 向に移動可能である。回転プラットフォーム１５を動かすことによって、解析平 面（Ｄ−Ｚ）は回転し、それによって、特に以下で説明するように３次元画像を 再構成する目的で、乳房２をいくつかの視点で撮像することが可能になる。
【００５５】
変換器列９は、解析平面が回転する前にいくつかの解析平面（例えば約３０ の解析平面）を撮像するように、さらに解析平面（Ｄ−Ｚ）に垂直に両方向矢印 １９によって示された両方向に並進運動することができてもよい。この目的のた めに、変換器の列は、例えば中心軸Ｚ０に対して半径方向に垂直な方向に沿って 延びており、回転プラットフォーム１５に固定されている横断方向ビーム２０上 を両方向１９にスライドするように取り付けられていてもよい。変換器列９は、 例えば、ビーム２０に内蔵されているナットやねじの組み立て品を制御する密閉 ステッピングモータによって、ビーム２０に沿って動作する調整システムによっ て移動させる。
【００５６】
変換器列９は、必要な場合、ビームに沿ってスライドし、列９が回転軸Ｒを 中心に前述の方向１９に平行にピボット運動可能な回り継ぎ手２１を介してビー ム２０に接続されていてもよい。したがって、立体図面において、任意で、方向 Ｄを変化させる、つまり、解析平面（Ｄ−Ｚ）を変化させずに方向Ｄの水平平面 （Ｘ−Ｙ）に対する傾斜角度を変化させてもよい。以下で説明するように、この 任意の動きによって、乳房２の画像が３次元で再構成されるように、乳房２が撮 像される視点も変化させることができる。
【００５７】
変換器列９のピボット運動は、例えば、密閉ステッピングモータ２２によっ て制御され、列９を両方向矢印２３によって定められる方向に変位させることが できる。
【００５８】
特に変換器の列９が乳房２の観測領域の所望の高さをカバーするのに十分に 長くない場合、固定プラットフォーム１６は、任意で、いくつかの高さで乳房を 撮像するために、垂直方向に移動させてもよい。
【００５９】
図示の例では、固定プラットフォーム１６の高さは、垂直軸Ｚに沿って両方 向矢印２６にて定義される方向にプラットフォーム１６を移動させるように、例 えば、支持ロッド２５を介してプラットフォーム１６に接続されている密閉電気 ステッピングモータ等によって動作するねじシリンダ２４で調整できる。
【００６０】
図３に示しているように、撮像デバイスは、さらに（キーボードなどの少な くとも１つの入力インターフェイス２８とスクリーンなどの出力インターフェイ ス２９とを含む）マイクロコンピュータ２７または任意の他の電子中央処理装置 を含んでもよく、以下を制御する。
【００６１】
−音響変換器１０（ＴＡ）
−シリンダ１３のモータ３０（ＭＯＴ．１）
−回転プラットフォーム１５の回転を制御するモータ３１（ＭＯＴ．２）
−変換器列９の方向１９のスライドを制御するモータ３２（ＭＯＴ．３）
−変換器列９の方向２３のピボット運動を制御するモータ２２（ＭＯＴ．４ ）
−シリンダ２４のモータ３３（ＭＯＴ．５）
−電子中央処理装置（ＣＰＵ）３４
−電子信号処理ラック３５のマイクロプロセッサ等
中央処理装置３４は、変換器列９のｎ個の変換器３６（Ｔ１、Ｔ２、．．． Ｔｎ）も制御する（ｎは、例えば、１２８以上）。
【００６２】
ラック３５は、任意で以下のものを含んでもよい。
【００６３】
−中央処理装置３４に接続されているメモリ３７（Ｍ）
−同様に中央処理装置３４に接続されている信号処理用の専用回路３８（Ｄ ＳＰ）
−中央処理装置３４に接続されている加算機３９（Ｓ）
−加算器に接続されており、各々が変換器３６に割り当てられているｎ個の メモリ４０（Ｍ１、Ｍ２、．．．Ｍｎ）（各メモリ４０は例えば１２８Ｍｂ以上 の容量を含んでもよい）
−各変換器３６と対応しているメモリ４０との間に介在しているｎ個のサン プラ４１（Ｅ１、Ｅ２、．．．Ｅｎ）
前述のデバイスは、以下のように動作する。
【００６４】
マイクロコンピュータ２７は、励起ステップ（ａ）の間に、音響変換器１０ によって粘弾性媒体２内に複数の持続弾性力学波を発生する。
【００６５】
さらに、粘弾性媒体２内の力学波の伝搬は、励起ステップ（ａ）に伴う観測 ステップ（ｂ）の間に、エコーグラフとして観測される。この目的のために、マ イクロコンピュータ２７は、列９の複数の変換器３６を介して複数の圧縮超音波 （例えば０．５と１００ＭＨｚとの間の、好ましくは０．５と１５ＭＨｚとの間 の、例えば４ＭＨｚ付近の周波数で）を伝送させる。
【００６６】
これらの圧縮波は、カプリングジェルつまり液体６を通過した後、媒体２を 貫通し、そのとき媒体２内に含まれている複数の散乱粒子によって反射され、そ れによって媒体２の運動を監視することができる。これらの対象となる散乱粒子 は、媒体２内の任意の異質成分からなり、特に、医療用途の場合、人体組織に存 在している複数のコラーゲン粒子からなる。
【００６７】
したがって、伝送される複数の圧縮波は、例えば、「平面」圧縮超音波（つ まり、この場合、波面が解析平面（Ｄ−Ｚ）内の直線である波）または、例えば 、さまざまな変換器３６を介して複数の音響信号を伝送することによって発生さ れた波の、媒体２内の観察領域全体を明らかにする、任意の他の種類の焦点波ま たは非焦点波（focused or unfocused wave）である。
【００６８】
前述の観測ステップ（ｂ）は、変換器列９の各位置で以下のサブステップを 含んでもよい。
【００６９】
（ｂ１）マイクロコンピュータ２７は、列９に、圧縮超音波の連続するショ ットを粘弾性媒体２内に伝送させる。
【００７０】
（ｂ２）マイクロコンピュータ２７は、列９に、粘弾性媒体２から受け付け た、粘弾性媒体２内の複数の散乱粒子と相互作用する複数の圧縮超音波によって 発生された複数のエコーを含む複数の音響信号を、リアルタイムで検出および記 録させる。これらのエコーは、粘弾性媒体の運動の連続する画像に（直接または 間接的に）対応する。
【００７１】
ステップ（ｂ１）間に、圧縮超音波は、１秒あたり１０および１００００シ ョットの間、好ましくは１秒あたり１０および５０００ショットの間の割合（こ の割合は、圧縮波が媒体２内を往復する時間、つまり、媒体２の方向Ｄの厚さに よって制限される。これは、新しい圧縮波が送信される前に、プローブ３６にて 受信された圧縮波によって発生された全てのエコーが、プローブ３６にて受信さ れなければならないためである）で伝送される。
【００７２】
各圧縮超音波は、複数の剪断波よりもはるかに高い伝搬速度（例えば、人体 内を約１５００ｍ／ｓ）で媒体内２を伝搬し、媒体２内の複数の散乱粒子と相互 作用することによって、複数のエコー、または、エコーグラフの分野において「 スペックル（speckle）雑音」という用語でそれ自体が知られている、他の同等 の乱れを発生する。
【００７３】
このスペックル雑音は、各圧縮超音波の発射後に、サブステップ（ｂ２）間 に、変換器Ｔ１．．．Ｔｎによって検出される。したがって、ｊ番目の発射後に 各変換器Ｔｉによって検出された信号ｓｉｊ（ｔ）は、対応するサンプラＥｉに よって、先ず、高周波（例えば３０から１００ＭＨｚの範囲の周波数）でサンプ リングされ、そして、リアルタイムでデジタル化（例えば１２ビットを超える） される。
【００７４】
このようにサンプリングおよびデジタル化された信号ｓｉｊ（ｔ）は、再び リアルタイムで、対応するメモリＭｉに保存される。
【００７５】
また、変換器列９の各位置に対して、マイクロコンピュータ２７は、通常は 以降に、以下の内容の処理ステップ（ｃ）を実行する。
【００７６】
（ｃ１）マイクロコンピュータ２７は、サブステップ（ｂ２）間で、連続す る複数の伝搬画像を得るために、粘弾性媒体２から受信した連続する複数の音響 信号を処理する。
【００７７】
（ｃ２）マイクロコンピュータ２７は、観察領域の様々な点で、粘弾性媒体 ２の少なくとも１つの運動パラメータ（例えば、変位または変位速度）を求める 。
【００７８】
サブステップ（ｂ２）間で、変換器列９の同じ位置での計測に対応している すべての信号ｓｉｉ（ｔ）が保存された後に、中央処理装置３４は、これらの信 号を加算回路Ｓによって、サブステップ（ｃ１）に対応している従来のビーム構 成手順によって再処理させる（または、中央処理装置自体がこの処理を実行した り、この処理がマイクロコンピュータ２７内で実行されてもよい）。
【００７９】
したがって、各々がｊ番目の発射後の平面（Ｄ−Ｚ）内の観測領域の画像に 対応している信号Ｓｊ（ｒ，ｚ）が生成される。
【００８０】
例えば、信号Ｓｊ（ｔ）は以下の式によって求めることができる。
【００８１】
【数１】
【００８２】
ここで、
−ｓｉｊは、ｊ番目の圧縮波の発射後にｉ番目の変換器によって検出された 未処理信号である
−ｔ（ｒ、ｚ）は、圧縮超音波が観測領域内の座標（ｒ、ｚ）の点に到達す るのにかかった時間である。ここで、ｊ番目の発射の開始時がｔ＝０で、ｒは、 Ｄ方向に原点から計測した横座標であり、ｚは、対象の点の垂直座標である。
【００８３】
−ｄｉ（ｒ，ｚ）は、観測領域内の座標（ｒ、ｚ）の点と、ｉ番目の変換器 との間の距離、または、その距離の概算値である。
【００８４】
−Ｖは、観測される粘弾性媒体内の、複数の音響圧縮超音波の平均伝搬速度 である。
【００８５】
−αｉ（ｒ、ｚ）は、アポダイゼーションの原理を考慮した、重み付け係数 である（実際には、多くの場合、αｉ（ｒ、ｚ）＝１と見なすことができる）。
【００８６】
任意事項であるビーム構成ステップ後、中央処理装置３４は、複数の画像信 号ｓｊを、（ｒ、ｚ）をラック７が含む中央メモリＭに保存する。ここで、ｊは 圧縮波発射の番号である。これらの信号は、自体が画像処理を行うマイクロコン ピュータ４内にも保存されてもよい。
【００８７】
サブステップ（ｃ２）の直後または後に、相関計算（correlation）と有利 には相互相関計算(intercorrelation)によって、１対で、好ましくは基準画像と １対で処理される、これらの画像は、以下のものである。
【００８８】
−前述のようにして以前に求められ、以降の複数の変位画像（または、例え ば３０の変位画像など、限定された数の以降の変位画像）の基準画像として使用 される変位画像。
【００８９】
−または、予備の初期観測ステップ（ａ０）中に、剪断波を発生させる励起 ステップ（ａ）の前に、変換器列９に１つまたは２つ以上の超音波を伝送させて 前述の複数の連続変位画像と同様にして得られた画像。（そのため、励起段階の 前にいくつかの圧縮超音波が伝送され、粘弾性媒体の複数の反射粒子と相互作用 する各圧縮超音波によって発生される複数のエコーが記録され、これらのエコー は粘弾性媒体と連続するいくつかの予備画像に対応しており、粘弾性媒体の初期 画像は、連続する複数の初期画像を組み合わせることによって、具体的には、復数の初期画像の複数の画素値を平均することによって求められる）。
【００９０】
前述の相関計算は、例えば、ラック３５の一部を構成しているＤＳＰ回路内 で実行したり、中央処理装置３４またはマイクロプロセッサ２７内でプログラム したりすることができる。
【００９１】
この相互相関処理中に、例えば、相互相関関数＜Ｓｊ（ｘ，ｙ），Ｓｊ（ｘ ，ｙ）＞は、超音波エコーを発生させる各粒子５によって受ける変位を求めるた めに最大化される。
【００９２】
そのような相互相関計算の例は、Ｏ‘Ｄｏｎｎｅｌ他、（“Ｉｎｔｅｒｎａ ｌ ｄｉｓｐｌａｃｅｍｅｎｔ ａｎｄ ｓｔｒａｉｎ ｉｍａｇｉｎｇ ｕｓ ｉｎｇ ｓｐｅｃｋｌｅ ｔｒａｃｋｉｎｇ（スペックルトラッキングを使用し た内部変位と歪みの撮像）”、ＩＥＥＥ Ｔｒａｎｓａｃｔｉｏｎ ｏｎ Ｕｌ ｔｒａｓｏｎｉｃｓ， Ｆｅｒｒｏｅｌｅｃｔｒｉｃｓ ａｎｄ Ｆｒｅｑｕｅ ｎｃｙ Ｃｏｎｔｒｏｌ Ｖｏｌ．４１， Ｎｏ．３，Ｍａｙ １９９４，ｐｐ ３１４−３２５）とＯｐｈｉｒ他、（“Ｅｌａｓｔｇｒａｐｈｙ：ａ ｑｕａｎ ｔｉｔａｔｉｖｅ ｍｅｔｈｏｄ ｆｏｒ ｉｍａｇｉｎｇ ｔｈｅ ｅｌａｓ ｔｉｃｉｔｙ ｏｆ ｂｉｏｌｏｇｉｃａｌ ｔｉｓｓｕｅｓ（エラストグラフ ィー：生物組織の弾性の定量的撮像方法）”、超音波撮像、Ｖｏｌ．１３、ｐｐ １１１−１３４、１９９１）による従来技術で示されている。
【００９３】
Ｔａｎｔｅｒ他（Ｍ．Ｔａｎｔｅｒ， Ｊ．Ｂｅｒｃｏｆｆ， Ｌ．Ｓａｎ ｄｒｉｎ ａｎｄ Ｍ．Ｆｉｎｋ、“Ｕｌｔｒａｆａｓｔ ｃｏｍｐｏｕｎｄ ｉｍａｇｉｎｇ ｆｏｒ ２Ｄ ｍｏｔｉｏｎ ｖｅｃｔｏｒ ｅｓｔｉｍａｔ ｉｏｎ： Ａｐｐｌｉｃａｔｉｏｎ ｔｏ ｔｒａｎｓｅｎｔ ｅｌａｓｔｇｒ ａｐｈｙ（２Ｄ運動ベクトル予測のための超高速複合撮像：過渡エラストグラフ ィーへの応用）”、ＩＥＥＥ、Ｔｒａｎｓａｃｔｉｏｎ Ｕｌｔｒａｓｏｎｉｃ ｓ， Ｆｅｒｒｏｅｌｅｃｔｒｉｃｓ ａｎｄ Ｆｒｅｑｕｅｎｃｙ Ｃｏｎｔ ｒｏｌ ４９（１０）、ｐｐ１３６３−１３７４、２００２）によって提案され ている相互相関技法を一般化すると、変換器列９の各位置について、剪断波の影 響下にある解析平面（Ｄ−Ｚ）の観測領域の各位置
【００９４】
【数２】
【００９５】
で複数の剪断波によって発生した一式の変位ベクトル
【００９６】
【数３】
【００９７】
が得られる（これらの変位ベクトルは２次元であり、方向Ｄに沿った成分と方 向Ｚに沿った成分とを含む）。
【００９８】
また、前述の軸Ｒを中心に列をピボット運動させることによって変換器列９ の場所を修正しながら、一連の複数の超音波の発射をもう１回行うことによって 、解析平面（Ｄ−Ｚ）内で、計測を実行する、または、複数の変位を計測するた めの前述の技法を洗練する（つまり、観測領域を拡大する）ことも可能である。
【００９９】
それから、複数の変位の計測は、ビーム２０に沿って連続するいくつかの位 置に変換器列９を変位させることによって、いくつかの平行な解析平面（Ｄ−Ｚ ）内で繰り返される。したがって、例えば、すべての平行な平面について約１か ら２秒の合計時間で、乳房の撮像の場合では約３０の平行な平面内で、変位を計 測することも可能になる。
【０１００】
それから、プラットフォーム１５は、予め定められた角度で回転し、それに よって解析平面（Ｄ−Ｚ）も回転させ、新しい方向Ｄに平行ないくつかの解析平 面内で複数の変位を計測する前述の動作が前述のように反復される。したがって 、新しい方向Ｄに沿った成分と方向Ｚに沿った成分とを含む２次元の複数の変位 場
【０１０１】
【数４】
【０１０２】
が、求められる。
【０１０３】
方向Ｄの少なくとも２つの向き（例えば、互いに９０°の２つの向き）、ま たは必要に応じて２つより多い向き（例えば、３６までの角度の向き）を含む複
数の変位の計測後に、粘弾性媒体２内の３次元の変位の場
【０１０４】
【数５】
【０１０５】
を再構成することができる。
【０１０６】
そのためには、例えば、時間の関数として、軸Ｘ、Ｙ、Ｚのそれぞれに沿っ た複数の変位の値が収集され、必要に応じて、さまざまな計測値を時間的に再設 定後（そのために、信号は、周期的な変調を伴って変換器１０によって送信され 、そして、時間の原点は、超音波の発射の各シーケンスの前の周期的な変調にお ける全く同一の点とされる）、これらの値が平均と補間とによって組み合わされ て、３次元の変位ベクトルの１つの場
【０１０７】
【数６】
【０１０８】
となる。
【０１０９】
次に、この変位の場
【０１１０】
【数７】
【０１１１】
（または他の移動パラメータ）は、音響変換器１０にて発生された粘性力学波 の圧縮成分の影響を除去するように補正される。このために、回転演算子が、ベ クトル場
【０１１２】
【数８】
【０１１３】
を求めるために、場
【０１１４】
【数９】
【０１１５】
に適用される。
【０１１６】
したがって、音響変換器１０にて発生した粘性力学波の剪断成分の影響だけ が考慮される。
【０１１７】
この一式のベクトル
【０１１８】
【数１０】
【０１１９】
は、メモリＭまたはマイクロコンピュータ４に保存され、そして、例えば、具 体的にはマイクロコンピュータ４の画面４ａを使用して、ｑの成分の値がグレー レベルまたは色レベルなどの光学的パラメータとして表示されるゆっくりしたフ ィルムの形態で表示することができる。
【０１２０】
したがって、医療用途での健康な組織とガン性の組織など、媒体２内の特性 が異なる領域間の剪断波の伝搬の違いが完全に表示することができる。
【０１２１】
この弾性波の伝搬のフィルムは、さらに、前述の装置１によって同様に生成 可能な通常のエコーグラフの画像にさらに重ねて表示することが可能である。
【０１２２】
さらに、観測されている媒体２の各点の複数の変位だけでなく、媒体２の歪 み、つまり、各変位ベクトルの両成分の複数の空間変数についての微分を成分と するベクトルも計算することもできる。これらの歪みベクトルを、剪断波の伝搬 を明確に表示するために、変位ベクトルのようにフィルムの形態で使用すること が可能で、これらの歪みベクトルは、観測している媒体２に対するプローブ３６ の変位を除外することができるという利点も含む。これらの歪みベクトルは、変 位ベクトルの場合に前述したように、音響変換器１０の粘性力学波の圧縮要素の 影響を取り除くために回転演算子を適用することによって処理することができる 。
【０１２３】
次に、マイクロコンピュータ２７は、媒体２内の剪断波の伝搬を記述してい る少なくとも１つのパラメータが、回転演算子を適用することによって処理され た運動パラメータの変化（変位または歪み速度）から計算される、変位場または 歪み場に基づいてマッピングステップ（ｄ）を実行することが有利である。
【０１２４】
マッピングステップ間に計算される剪断波の伝搬パラメータは、剪断波の速 度ＣT、剪断弾性率μまたはヤング率Ｅ（＝３μ）、複数の剪断波の減衰、剪断 弾性、剪断粘性ζ、または組織の力学緩和時間、または異方性の媒体の場合これ らのパラメータうちのいくつかのパラメータから選択される。
【０１２５】
例えば、観測領域の様々な点で剪断弾性率μと剪断粘性ζを計算することが できる。
【０１２６】
そのために、以下の伝搬方程式（１）を使用することができ（媒体２は、１ 次近似として局所的に一様かつ等方な粘弾性媒体と見なせると仮定する）、複数 の粘性力学波によって媒体２内の各位置
【０１２７】
【数１１】
【０１２８】
で発生する変位
【０１２９】
【数１２】
【０１３０】
は以下の式に従う。
【０１３１】
【数１３】
【０１３２】
ここでρは媒体の密度、μは剪断弾性率（再構成の対象）、λは第２ラメ係数 、
【０１３３】
【数１４】
【０１３４】
は媒体内への拡散を考慮した剪断粘性、ζは圧縮波の粘性である。剪断弾性 率μと第２ラメ係数λとは、ヤング率Ｅとポアソン比σに関して以下の式で与え られる。
【０１３５】
【数１５】
【０１３６】
各ベクトル場
【０１３７】
【数１６】
【０１３８】
は、非回転部分
【０１３９】
【数１７】
【０１４０】
、圧縮不能部分
【０１４１】
【数１８】
【０１４２】
、および非回転で圧縮不能な調和部分
【０１４３】
【数１９】
【０１４４】
とに分解することができる。調和部分
【０１４５】
【数２０】
【０１４６】
は、大局的な複数の変位を考慮し、そのため無視でき、ここでは、この分解を 使用して複数の波の場合を研究しているため、式（１）を複数の剪断波の伝搬を 表している式と複数の圧縮波の伝搬を表している式との２つの式に分離できるこ とを示すことが可能であり、各々の伝搬速度はｃTとｃLとであって（粘性の影 響は無視）、つまり以下のようになる。
【０１４７】
【数２１】
【０１４８】
式（１）に回転演算子を適用することによって、変換器１０にて発生された 弾性力学波の圧縮成分の影響をもはや有していない簡略化された式が得られる。
【０１４９】
【数２２】
【０１５０】
それから、式（４）からパラメータμとζと３次空間微分の現在の値を計算 することができる。式（４）は、例えば、系の転置行列（疑似逆行列）が乗算さ れることによって簡単に解くことができる。式（５）は、２つの未知数（密度ρ は水の密度に等しいとみなされる）に対して実際には６つの方程式（複素数表記 で３つが
【０１５１】
【数２３】
【０１５２】
の実部で３つが虚部）を示していることがわかる。
【０１５３】
したがって、定常状態では、以下の表記のもとで、
【０１５４】
【数２４】
【０１５５】
使用可能な方程式は、以下のように表すことができる。
【０１５６】
【数２５】
【０１５７】
転置行列を乗じることによって、以下の式が得られる。
【０１５８】
【数２６】
【０１５９】
ここでｉ＝ｘ，ｙ，ｚであって、複数の同一の添え字に対するアインシュタイ ン規約が適用されている。定常状態の場合、μとζとの間には相関はなく、その ため互いに独立しているこれらの２つの計算されたパラメータをマッピングする のがより容易になる。
【０１６０】
ただし、局所的に一様で等方性の媒体という近似ができない場合にも、所望 の複数の伝搬パラメータを求めるために同様な計算を使用することができる。こ の場合、より複雑な伝搬方程式が使用される、この伝搬方程式は、以下のように なる。
【０１６１】
【数２７】
【０１６２】
ここでμiは変位ベクトル
【０１６３】
【数２８】
【０１６４】
のｉ番目の成分であり、ρは媒体２の密度であり、λiklmは階数４の弾性テン ソルであり、ηiklmは階数４の粘性テンソルである（複数の同一の添え字に対す るアインシュタイン規約を使用）。
【０１６５】
この式（１′）は、媒体が横方向に等方性を含む、つまり、媒体が一方向に だけ揃えられている繊維を含むと仮定できる場合に簡単になる。この場合、物質 を記述している複数のパラメータを以下の２つのグループに分割することができ る。
【０１６６】
−２つの剪断係数（μ□，μ⊥）を含む１つのグループは複数の繊維の方向 に水平（□）および垂直（⊥）な剪断波の伝搬を記述している。
【０１６７】
−他方のグループ（λ□，λ⊥，λM）は様々な方向の複数の縦波（つまり 複数の圧縮波）の伝搬を記述している
歪みテンソルと応力テンソルとの間の関係は次の式で示される。
【０１６８】
【数２９】
【０１６９】
ここでｕik＝１／２（∂xkｕi＋∂xiuk）は歪みテンソルであってσikは応 力テンソルである。
【０１７０】
変位ベクトル
【０１７１】
【数３０】
【０１７２】
の対応する伝搬方程式は次の式で与えられる。
【０１７３】
【数３１】
【０１７４】
ここでζは剪断粘性である。複数の縦波の粘性（ξ）は、比較的低い周波数（ ＨｚまたはｋＨｚの範囲）では無視できる。さらに、
【０１７５】
【数３２】
【０１７６】
はポアソン比の値のために小さい値を含む。したがって、ξ、ζおよび
【０１７７】
【数３３】
【０１７８】
は値が小さく、そのため式９の右辺の第３項は合理的に無視することができる 。式８を式９に代入すると以下の式を得ることができる。
【０１７９】
【数３４】
【０１８０】
ここで表記τ＝μ□−μ⊥であって、縦波のすべての弾性率は等しい、つまり λ＝λ□＝λ⊥＝λMと仮定している。
【０１８１】
回転演算子の適用後に以下の式が得られる。
【０１８２】
【数３５】
【０１８３】
この式（１１）によって、特にμ⊥、τおよびζを計算するために使用可能 な６つの方程式が実際に得られる。
【０１８４】
ただし、本発明は媒体２の運動パラメータ（変位等）を計測するどのような 方法にも適用可能である。したがってこのパラメータは必要に応じてＭＲＩによ って求められるであろう。
【０１８５】
説明した方法は、例えばＰＣＴ特許出願公開明細書第０４／８５９７８号に 記述されているような媒体２の音響光学撮像と組み合わせることもできる。
【図面の簡単な説明】
【０１８６】
【図１】本発明の一実施形態の複数の剪断波によって撮像する撮像装置の概略 図である。
【図２】図１の装置の一部の詳細な透視図である。
【図３】図１の装置のブロック図である。
【技術分野】
【０００１】
本発明は、複数の剪断波を用いる撮像方法および装置に関する。
【０００２】
特に、本発明は、剪断成分および圧縮成分を含む力学波（onde mecanique） が粘弾性媒体内で発生され、かつ、この力学波の伝搬が続いている間における、 様々な点での粘弾性媒体の運動を記述する少なくとも一つのパラメータを求める 撮像方法に関する。
【０００３】
したがって、特に、粘弾性媒体の他の部分と硬度が異なる複数の領域、また は、粘弾性媒体の他の部分と力学緩和時間が異なる複数の領域を特定するために 、定性的および／または定量的な分析を実施することができる。
【背景技術】
【０００４】
ＰＣＴ特許出願公開明細書第０４／２１０３８号は、そのような方法の例を 示しており、対象とする力学波、本方法では剪断波、を生成する局所的なパルス 状の推力の離れた位置での発生に続く、平面内の運動（特に変位）を記述してい るパラメータが求められる。
【発明の開示】
【発明が解決しようとする課題】
【０００５】
この方法はすでに完全に満足できるが、本発明の目的は、その精度を改良し 、必要であれば、３次元の撮像に使用できるように、さらにこれを完全にするこ とである。
【課題を解決するための手段】
【０００６】
この目的のために、対象となる方法は、運動パラメータが、上記の力学波の 圧縮成分による影響を取り除くように処理される補正ステップを含むことを特徴 とする。
【０００７】
したがって、パルス状の推力が使用されるときでも、それによって発生する 弾性力学波の圧縮成分が運動パラメータの計測を妨害する場合、剪断成分の効果 だけが考慮される。例えば３次元撮像などの比較的長い計測を実施するために、 推力がなくても持続する力学波が使用される場合、力学波の圧縮成分は、運動パ ラメータの計測中も依然として存在し、その結果、この計測を大きく妨害する。
本発明による処理は、この妨害を取り除くことができ、必要であれば、粘弾性媒 体の複数の粘弾性パラメータを２または３次元にマッピングすることができる。
【０００８】
本発明による方法の好ましい実施形態では、以下の複数の構成の一つまたは ２つ以上を、さらに任意で用いてもよい。
【０００９】
−上記の運動パラメータの回転が補正ステップ中に計算される。
【００１０】
−粘弾性媒体内の複数の点における力学波の伝搬が、エコーグラフで観測さ れる少なくとも１つの観測ステップ（ｂ）を含む。
【００１１】
−観測ステップ（ｂ）は、以下のサブステップを含む。
【００１２】
複数の圧縮超音波（ondes ultrasonores de compression）の連続するショ ットが、複数の変換器の列によって、（焦点が合わされているかどうかにかかわ らず）少なくとも１秒あたり１０ショットの割合で、粘弾性媒体内に伝送される サブステップ（ｂ１）。
【００１３】
粘弾性媒体に相互作用する複数の圧縮超音波によって発生した複数のエコー を含む粘弾性媒体から受信した複数の音響信号がリアルタイムで検出および記録 されるサブステップ（ｂ２）。
【００１４】
−以下のステップが実施される少なくとも１つの処理ステップ（ｃ）をさら に含む。
【００１５】
力学波の連続する複数の伝搬画像を求めるために、サブステップ（ｂ２）の 間に粘弾性媒体から受信した複数の連続する音響信号が処理されるステップ（ｃ １）。
【００１６】
連続する複数の伝搬画像から粘弾性媒体の運動パラメータが様々な点で求め られるステップ（ｃ２）。
【００１７】
−粘弾性媒体内での力学波の連続する複数の伝搬画像は、いくつかの視点に 沿って生成され、これらの連続する画像は、運動パラメータのベクトル値を求め るために組み合わされる。
【００１８】
−運動パラメータの３次元ベクトル値が求められる。
【００１９】
−いくつかの視点に沿った力学波の連続する複数の伝搬画像を生成するため に、解析平面内の伝送／検出方向に沿って複数の圧縮超音波を伝送し検出する複
数の超音波変換器の列（複数の圧電変換器の直線状の列、または位相列、１．５
Ｄ列、１．７５Ｄ列、または２Ｄ配列などの任意の他の列）が使用され、直線状
の列は、粘弾性媒体の外側で変位することで、解析平面を回転させる（これらの
構成は、前述の構成とは独立して任意に使用することが可能で、そのため、本発
明の主題は、これらの構成自体でもある）。
【００２０】
−列は、解析平面を回転しない並進運動によって変位し、その後、列は、回 転することで、解析平面を回転させる。
【００２１】
−複数の変換器の列の伝送／検出方向は、解析平面を並進も回転もさせずに 変化し、その後、列は、並進および回転することで、解析平面を回転させる。
【００２２】
−運動パラメータは粘弾性媒体の変位、変位速度および応力から選択される 。
【００２３】
−粘弾性媒体内の伝搬パラメータのマップを求めるために、粘弾性媒体内の 少なくとも特定の複数の点での、力学波の剪断成分の伝搬を記述する少なくとも １つのパラメータが、運動パラメータの時間の経過に伴う変化から計算されるマ ッピングステップ（ｄ）をさらに含む。
【００２４】
−マッピングステップ中に計算される剪断波の伝搬パラメータは、複数の剪 断波の速度、剪断弾性率、ヤング率、複数の剪断波の減衰、剪断弾性、剪断粘性 、および力学緩和時間、または、異方性媒体の場合これらのパラメータうちのい くつかのパラメータから選択される。
【００２５】
また、本発明の主題は、剪断成分と圧縮成分とを含む力学波を粘弾性媒体内 で発生させる手段と、力学波の伝搬中に様々な点で粘弾性媒体の運動を記述する 少なくとも１つのパラメータを求める手段とを含み、力学波の圧縮成分による影 響を取り除くように運動パラメータを処理する補正手段を含むことを特徴とする 撮像装置である。
【００２６】
本発明の装置の好ましい実施形態では、以下の複数の構成の１つまたは２つ 以上をさらに任意で用いてもよい。
【００２７】
−補正手段は、運動パラメータの回転を計算するように構成されている。
【００２８】
−力学波の伝搬を粘弾性媒体内の複数の点でエコグラフィーによって観測す る複数の超音波変換器の制御されている列を含む。
【００２９】
−装置は、以下の複数の手段を含む。
【００３０】
複数の変換器の列によって、少なくとも１秒あたり１０ショットの割合で（ 焦点が合わされているかどうかにかかわらず）複数の圧縮超音波の連続するショ ットを粘弾性媒体内に伝送する手段（ｂ１）。
【００３１】
粘弾性媒体の複数の反射粒子に相互作用する複数の焦点が合っていない圧縮 超音波によって発生した複数のエコーを含む、粘弾性媒体から受信した複数の音 響信号をリアルタイムで検出および記録する手段（ｂ２）。
【００３２】
−装置は、以下のために構成されている複数の処理手段を含む。
【００３３】
力学波の連続する複数の伝搬画像を求めるために、粘弾性媒体から受信した 複数の連続する音響信号を処理する（ｃ１）。
【００３４】
連続する複数の伝搬画像から粘弾性媒体の運動パラメータを観測フィールド 内の様々な点で求める（ｃ２）。
【００３５】
−複数の処理手段は、粘弾性媒体内でいくつかの視点に沿って力学波のいく つかの連続する伝搬画像を生成し、連続する画像を運動パラメータのベクトル値 を求めるために組み合わせるように構成されている。
【００３６】
−複数の処理手段は運動パラメータの３次元ベクトル値を求めるように構成 されている。
【００３７】
−複数の超音波変換器の列はいくつかの視点に沿った力学波の連続する複数 の伝搬画像を生成するために、解析平面内の伝送／検出方向に沿って複数の圧縮 超音波を伝送して検出し、列（複数の圧電変換器の直線状の列、または位相列、 １．５Ｄ列、１．７５Ｄ列、または２Ｄ配列などの任意の他の列）を粘弾性媒体 の外側で変位させる複数の変位手段を含み、複数の手段は解析平面を回転させる ように構成されている（これらの構成は、前述の構成とは独立して任意に使用す ることが可能で、そのため、本発明の主題は、これらの構成自体でもある）。
【００３８】
−装置は、複数の変位手段を制御するように構成されている制御手段を含み 、制御手段は、解析平面を回転させずに列を並進運動によって変位させ、その後 、列を回転させることで、解析平面を回転させる。
【００３９】
−装置は、複数の変位手段を制御するように構成されている制御手段を含み 、制御手段は、解析平面を並進も回転もさせずに、複数の変換器の列の伝送／検 出方向を変化させ、その後、列を並進および回転させることで、解析平面を並進 および回転させるように構成されている。
【００４０】
装置は、粘弾性媒体内の少なくとも特定の複数の点での力学波の剪断成分の 伝搬を記述する少なくとも１つのパラメータを計算するように構成されているマ ッピング手段を含む。
【００４１】
本発明の他の特徴および利点は、添付の図面と共に非限定的な例として示し ているその実施形態の１つの以下の説明を通して明らかになるであろう。
【発明を実施するための最良の形態】
【００４２】
さまざまな図において、同一の符号は、同一のまたは同様な要素を示してい る。
【００４３】
図１に示している撮像装置１は、圧縮超音波を散乱させる粘弾性媒体２内の 複数の弾性剪断波の伝搬を研究することを目的とし、粘弾性媒体２は、例えば、 以下のようなものである。
【００４４】
−不活性体、特に、産業用途のための品質管理の場合、または
−生体、例えば、医療用途における患者の体の一部
図に示している例では、撮像装置１は、特に、女性患者３の乳房の医療用撮 像のために設計され、そのため、乳房は、撮像される媒体２を構成している。
【００４５】
この目的のために、撮像装置１は、例えば、患者がうつぶせに横たわる検査 テーブル４を含んでいてもよい。このテーブルは、患者の少なくとも１つの乳房 ２を貫通させることができる開口５を備えており、貫通した乳房は、ジェルまた は液体、好ましくは水などの超音波カプリング剤の槽６に垂れ下がる。この液体 槽６は、容器７によって境界が定められおり、容器７は超音波変換器の列９と乳 房に接触している振動発生器８も含む。容器７は、特に振動発生器８が乳房２に 接触して適切に配置する、および、変換器列９が乳房に沿って適切に配置されて いることを確認するために透明であってもよい。
【００４６】
検査テーブル４は、必要に応じて、患者の両方の乳房が通過し、同じ液体槽 ６または２つの別個の槽に垂れ下がるように構成されていてもよい。つまり、撮 像装置は、患者の両方の乳房を同時に撮像するように任意に構成されていてもよ い。
【００４７】
振動発生器８は、複数の力学（弾性）波、特に圧縮成分と剪断成分とを含む 複数の持続波（ondes entretenues）を患者の乳房２の内部に伝送するために、 例えば、患者の乳房２の端部に接触するように配置される高さ調整プレート１１ 内で、低周波の力学振動を発生する音響変換器１０を含んでもよい。対象として いる複数の力学波は、例えば、０．１Ｈｚと１０ｋＨｚとの間の周波数を含んで もよい。
【００４８】
ただし、上記の複数の力学波は、例えば、以下のようなもので、前述の振動 発生器以外のもので発生させることができる。
【００４９】
−複数の低周波を発生し、乳房２に対して異なる方法で配置されている変換 器。
【００５０】
−例えば、上述の文献ＰＣＴ特許出願公開明細書第０４／２１０３８号、ま たは、米国特許出願第５ ６０６９７１号で説明されたように、放射圧力を使用 して時間変調されており焦点が合わされている複数の超音波によって複数の低周 波力学的信号を遠隔発生する超音波変換器の列。
【００５１】
図１に示している例では、プレート１１の高さは、乳房２に接触するように 調整可能である。この調整は、任意で、例えば、密封電気ステッピングモータな どで作動するねじシリンダ１３で行ってもよい。
【００５２】
シリンダ１３は、垂直軸Ｚに沿って、両方向矢印１４で示された両方向にプ レート１１および音響変換器１０を移動させるように、例えば、支持ロッド１２ を介して、プレート１１に接続されていてもよい。
【００５３】
図１および図２に示しているように、複数の変換器の列９は、例えば、直線 状列の形態でもよく、この場合、その複数の変換器は、垂直平面内に揃えられ、 複数の超音波を送受信するように、垂直解析平面（Ｄ−Ｚ）内の同一の送信／検 出方向Ｄに対して全て平行に構成される。
【００５４】
変換器列９は、例えば、固定プラットフォーム１６によって支持される回転 プラットフォーム１５に取り付けられていてもよい。２つのプラットフォーム１ ５、１６は水辺平面（Ｘ−Ｙ）内に配置され、患者の乳房２が中央通路１７を自 由に通過できるように開口していることが有利である。回転プラットフォーム１ ５は、例えば、プラットフォーム１５および１６の間に組み込まれている密閉電 気ステッピングモータ（図１および２では、図示しない）で動作する調整システ ムによって、固定プラットフォーム１６に対して垂直中央軸Ｚ０を中心に角度方 向に移動可能である。回転プラットフォーム１５を動かすことによって、解析平 面（Ｄ−Ｚ）は回転し、それによって、特に以下で説明するように３次元画像を 再構成する目的で、乳房２をいくつかの視点で撮像することが可能になる。
【００５５】
変換器列９は、解析平面が回転する前にいくつかの解析平面（例えば約３０ の解析平面）を撮像するように、さらに解析平面（Ｄ−Ｚ）に垂直に両方向矢印 １９によって示された両方向に並進運動することができてもよい。この目的のた めに、変換器の列は、例えば中心軸Ｚ０に対して半径方向に垂直な方向に沿って 延びており、回転プラットフォーム１５に固定されている横断方向ビーム２０上 を両方向１９にスライドするように取り付けられていてもよい。変換器列９は、 例えば、ビーム２０に内蔵されているナットやねじの組み立て品を制御する密閉 ステッピングモータによって、ビーム２０に沿って動作する調整システムによっ て移動させる。
【００５６】
変換器列９は、必要な場合、ビームに沿ってスライドし、列９が回転軸Ｒを 中心に前述の方向１９に平行にピボット運動可能な回り継ぎ手２１を介してビー ム２０に接続されていてもよい。したがって、立体図面において、任意で、方向 Ｄを変化させる、つまり、解析平面（Ｄ−Ｚ）を変化させずに方向Ｄの水平平面 （Ｘ−Ｙ）に対する傾斜角度を変化させてもよい。以下で説明するように、この 任意の動きによって、乳房２の画像が３次元で再構成されるように、乳房２が撮 像される視点も変化させることができる。
【００５７】
変換器列９のピボット運動は、例えば、密閉ステッピングモータ２２によっ て制御され、列９を両方向矢印２３によって定められる方向に変位させることが できる。
【００５８】
特に変換器の列９が乳房２の観測領域の所望の高さをカバーするのに十分に 長くない場合、固定プラットフォーム１６は、任意で、いくつかの高さで乳房を 撮像するために、垂直方向に移動させてもよい。
【００５９】
図示の例では、固定プラットフォーム１６の高さは、垂直軸Ｚに沿って両方 向矢印２６にて定義される方向にプラットフォーム１６を移動させるように、例 えば、支持ロッド２５を介してプラットフォーム１６に接続されている密閉電気 ステッピングモータ等によって動作するねじシリンダ２４で調整できる。
【００６０】
図３に示しているように、撮像デバイスは、さらに（キーボードなどの少な くとも１つの入力インターフェイス２８とスクリーンなどの出力インターフェイ ス２９とを含む）マイクロコンピュータ２７または任意の他の電子中央処理装置 を含んでもよく、以下を制御する。
【００６１】
−音響変換器１０（ＴＡ）
−シリンダ１３のモータ３０（ＭＯＴ．１）
−回転プラットフォーム１５の回転を制御するモータ３１（ＭＯＴ．２）
−変換器列９の方向１９のスライドを制御するモータ３２（ＭＯＴ．３）
−変換器列９の方向２３のピボット運動を制御するモータ２２（ＭＯＴ．４ ）
−シリンダ２４のモータ３３（ＭＯＴ．５）
−電子中央処理装置（ＣＰＵ）３４
−電子信号処理ラック３５のマイクロプロセッサ等
中央処理装置３４は、変換器列９のｎ個の変換器３６（Ｔ１、Ｔ２、．．． Ｔｎ）も制御する（ｎは、例えば、１２８以上）。
【００６２】
ラック３５は、任意で以下のものを含んでもよい。
【００６３】
−中央処理装置３４に接続されているメモリ３７（Ｍ）
−同様に中央処理装置３４に接続されている信号処理用の専用回路３８（Ｄ ＳＰ）
−中央処理装置３４に接続されている加算機３９（Ｓ）
−加算器に接続されており、各々が変換器３６に割り当てられているｎ個の メモリ４０（Ｍ１、Ｍ２、．．．Ｍｎ）（各メモリ４０は例えば１２８Ｍｂ以上 の容量を含んでもよい）
−各変換器３６と対応しているメモリ４０との間に介在しているｎ個のサン プラ４１（Ｅ１、Ｅ２、．．．Ｅｎ）
前述のデバイスは、以下のように動作する。
【００６４】
マイクロコンピュータ２７は、励起ステップ（ａ）の間に、音響変換器１０ によって粘弾性媒体２内に複数の持続弾性力学波を発生する。
【００６５】
さらに、粘弾性媒体２内の力学波の伝搬は、励起ステップ（ａ）に伴う観測 ステップ（ｂ）の間に、エコーグラフとして観測される。この目的のために、マ イクロコンピュータ２７は、列９の複数の変換器３６を介して複数の圧縮超音波 （例えば０．５と１００ＭＨｚとの間の、好ましくは０．５と１５ＭＨｚとの間 の、例えば４ＭＨｚ付近の周波数で）を伝送させる。
【００６６】
これらの圧縮波は、カプリングジェルつまり液体６を通過した後、媒体２を 貫通し、そのとき媒体２内に含まれている複数の散乱粒子によって反射され、そ れによって媒体２の運動を監視することができる。これらの対象となる散乱粒子 は、媒体２内の任意の異質成分からなり、特に、医療用途の場合、人体組織に存 在している複数のコラーゲン粒子からなる。
【００６７】
したがって、伝送される複数の圧縮波は、例えば、「平面」圧縮超音波（つ まり、この場合、波面が解析平面（Ｄ−Ｚ）内の直線である波）または、例えば 、さまざまな変換器３６を介して複数の音響信号を伝送することによって発生さ れた波の、媒体２内の観察領域全体を明らかにする、任意の他の種類の焦点波ま たは非焦点波（focused or unfocused wave）である。
【００６８】
前述の観測ステップ（ｂ）は、変換器列９の各位置で以下のサブステップを 含んでもよい。
【００６９】
（ｂ１）マイクロコンピュータ２７は、列９に、圧縮超音波の連続するショ ットを粘弾性媒体２内に伝送させる。
【００７０】
（ｂ２）マイクロコンピュータ２７は、列９に、粘弾性媒体２から受け付け た、粘弾性媒体２内の複数の散乱粒子と相互作用する複数の圧縮超音波によって 発生された複数のエコーを含む複数の音響信号を、リアルタイムで検出および記 録させる。これらのエコーは、粘弾性媒体の運動の連続する画像に（直接または 間接的に）対応する。
【００７１】
ステップ（ｂ１）間に、圧縮超音波は、１秒あたり１０および１００００シ ョットの間、好ましくは１秒あたり１０および５０００ショットの間の割合（こ の割合は、圧縮波が媒体２内を往復する時間、つまり、媒体２の方向Ｄの厚さに よって制限される。これは、新しい圧縮波が送信される前に、プローブ３６にて 受信された圧縮波によって発生された全てのエコーが、プローブ３６にて受信さ れなければならないためである）で伝送される。
【００７２】
各圧縮超音波は、複数の剪断波よりもはるかに高い伝搬速度（例えば、人体 内を約１５００ｍ／ｓ）で媒体内２を伝搬し、媒体２内の複数の散乱粒子と相互 作用することによって、複数のエコー、または、エコーグラフの分野において「 スペックル（speckle）雑音」という用語でそれ自体が知られている、他の同等 の乱れを発生する。
【００７３】
このスペックル雑音は、各圧縮超音波の発射後に、サブステップ（ｂ２）間 に、変換器Ｔ１．．．Ｔｎによって検出される。したがって、ｊ番目の発射後に 各変換器Ｔｉによって検出された信号ｓｉｊ（ｔ）は、対応するサンプラＥｉに よって、先ず、高周波（例えば３０から１００ＭＨｚの範囲の周波数）でサンプ リングされ、そして、リアルタイムでデジタル化（例えば１２ビットを超える） される。
【００７４】
このようにサンプリングおよびデジタル化された信号ｓｉｊ（ｔ）は、再び リアルタイムで、対応するメモリＭｉに保存される。
【００７５】
また、変換器列９の各位置に対して、マイクロコンピュータ２７は、通常は 以降に、以下の内容の処理ステップ（ｃ）を実行する。
【００７６】
（ｃ１）マイクロコンピュータ２７は、サブステップ（ｂ２）間で、連続す る複数の伝搬画像を得るために、粘弾性媒体２から受信した連続する複数の音響 信号を処理する。
【００７７】
（ｃ２）マイクロコンピュータ２７は、観察領域の様々な点で、粘弾性媒体 ２の少なくとも１つの運動パラメータ（例えば、変位または変位速度）を求める 。
【００７８】
サブステップ（ｂ２）間で、変換器列９の同じ位置での計測に対応している すべての信号ｓｉｉ（ｔ）が保存された後に、中央処理装置３４は、これらの信 号を加算回路Ｓによって、サブステップ（ｃ１）に対応している従来のビーム構 成手順によって再処理させる（または、中央処理装置自体がこの処理を実行した り、この処理がマイクロコンピュータ２７内で実行されてもよい）。
【００７９】
したがって、各々がｊ番目の発射後の平面（Ｄ−Ｚ）内の観測領域の画像に 対応している信号Ｓｊ（ｒ，ｚ）が生成される。
【００８０】
例えば、信号Ｓｊ（ｔ）は以下の式によって求めることができる。
【００８１】
【数１】
【００８２】
ここで、
−ｓｉｊは、ｊ番目の圧縮波の発射後にｉ番目の変換器によって検出された 未処理信号である
−ｔ（ｒ、ｚ）は、圧縮超音波が観測領域内の座標（ｒ、ｚ）の点に到達す るのにかかった時間である。ここで、ｊ番目の発射の開始時がｔ＝０で、ｒは、 Ｄ方向に原点から計測した横座標であり、ｚは、対象の点の垂直座標である。
【００８３】
−ｄｉ（ｒ，ｚ）は、観測領域内の座標（ｒ、ｚ）の点と、ｉ番目の変換器 との間の距離、または、その距離の概算値である。
【００８４】
−Ｖは、観測される粘弾性媒体内の、複数の音響圧縮超音波の平均伝搬速度 である。
【００８５】
−αｉ（ｒ、ｚ）は、アポダイゼーションの原理を考慮した、重み付け係数 である（実際には、多くの場合、αｉ（ｒ、ｚ）＝１と見なすことができる）。
【００８６】
任意事項であるビーム構成ステップ後、中央処理装置３４は、複数の画像信 号ｓｊを、（ｒ、ｚ）をラック７が含む中央メモリＭに保存する。ここで、ｊは 圧縮波発射の番号である。これらの信号は、自体が画像処理を行うマイクロコン ピュータ４内にも保存されてもよい。
【００８７】
サブステップ（ｃ２）の直後または後に、相関計算（correlation）と有利 には相互相関計算(intercorrelation)によって、１対で、好ましくは基準画像と １対で処理される、これらの画像は、以下のものである。
【００８８】
−前述のようにして以前に求められ、以降の複数の変位画像（または、例え ば３０の変位画像など、限定された数の以降の変位画像）の基準画像として使用 される変位画像。
【００８９】
−または、予備の初期観測ステップ（ａ０）中に、剪断波を発生させる励起 ステップ（ａ）の前に、変換器列９に１つまたは２つ以上の超音波を伝送させて 前述の複数の連続変位画像と同様にして得られた画像。（そのため、励起段階の 前にいくつかの圧縮超音波が伝送され、粘弾性媒体の複数の反射粒子と相互作用 する各圧縮超音波によって発生される複数のエコーが記録され、これらのエコー は粘弾性媒体と連続するいくつかの予備画像に対応しており、粘弾性媒体の初期 画像は、連続する複数の初期画像を組み合わせることによって、具体的には、復数の初期画像の複数の画素値を平均することによって求められる）。
【００９０】
前述の相関計算は、例えば、ラック３５の一部を構成しているＤＳＰ回路内 で実行したり、中央処理装置３４またはマイクロプロセッサ２７内でプログラム したりすることができる。
【００９１】
この相互相関処理中に、例えば、相互相関関数＜Ｓｊ（ｘ，ｙ），Ｓｊ（ｘ ，ｙ）＞は、超音波エコーを発生させる各粒子５によって受ける変位を求めるた めに最大化される。
【００９２】
そのような相互相関計算の例は、Ｏ‘Ｄｏｎｎｅｌ他、（“Ｉｎｔｅｒｎａ ｌ ｄｉｓｐｌａｃｅｍｅｎｔ ａｎｄ ｓｔｒａｉｎ ｉｍａｇｉｎｇ ｕｓ ｉｎｇ ｓｐｅｃｋｌｅ ｔｒａｃｋｉｎｇ（スペックルトラッキングを使用し た内部変位と歪みの撮像）”、ＩＥＥＥ Ｔｒａｎｓａｃｔｉｏｎ ｏｎ Ｕｌ ｔｒａｓｏｎｉｃｓ， Ｆｅｒｒｏｅｌｅｃｔｒｉｃｓ ａｎｄ Ｆｒｅｑｕｅ ｎｃｙ Ｃｏｎｔｒｏｌ Ｖｏｌ．４１， Ｎｏ．３，Ｍａｙ １９９４，ｐｐ ３１４−３２５）とＯｐｈｉｒ他、（“Ｅｌａｓｔｇｒａｐｈｙ：ａ ｑｕａｎ ｔｉｔａｔｉｖｅ ｍｅｔｈｏｄ ｆｏｒ ｉｍａｇｉｎｇ ｔｈｅ ｅｌａｓ ｔｉｃｉｔｙ ｏｆ ｂｉｏｌｏｇｉｃａｌ ｔｉｓｓｕｅｓ（エラストグラフ ィー：生物組織の弾性の定量的撮像方法）”、超音波撮像、Ｖｏｌ．１３、ｐｐ １１１−１３４、１９９１）による従来技術で示されている。
【００９３】
Ｔａｎｔｅｒ他（Ｍ．Ｔａｎｔｅｒ， Ｊ．Ｂｅｒｃｏｆｆ， Ｌ．Ｓａｎ ｄｒｉｎ ａｎｄ Ｍ．Ｆｉｎｋ、“Ｕｌｔｒａｆａｓｔ ｃｏｍｐｏｕｎｄ ｉｍａｇｉｎｇ ｆｏｒ ２Ｄ ｍｏｔｉｏｎ ｖｅｃｔｏｒ ｅｓｔｉｍａｔ ｉｏｎ： Ａｐｐｌｉｃａｔｉｏｎ ｔｏ ｔｒａｎｓｅｎｔ ｅｌａｓｔｇｒ ａｐｈｙ（２Ｄ運動ベクトル予測のための超高速複合撮像：過渡エラストグラフ ィーへの応用）”、ＩＥＥＥ、Ｔｒａｎｓａｃｔｉｏｎ Ｕｌｔｒａｓｏｎｉｃ ｓ， Ｆｅｒｒｏｅｌｅｃｔｒｉｃｓ ａｎｄ Ｆｒｅｑｕｅｎｃｙ Ｃｏｎｔ ｒｏｌ ４９（１０）、ｐｐ１３６３−１３７４、２００２）によって提案され ている相互相関技法を一般化すると、変換器列９の各位置について、剪断波の影 響下にある解析平面（Ｄ−Ｚ）の観測領域の各位置
【００９４】
【数２】
【００９５】
で複数の剪断波によって発生した一式の変位ベクトル
【００９６】
【数３】
【００９７】
が得られる（これらの変位ベクトルは２次元であり、方向Ｄに沿った成分と方 向Ｚに沿った成分とを含む）。
【００９８】
また、前述の軸Ｒを中心に列をピボット運動させることによって変換器列９ の場所を修正しながら、一連の複数の超音波の発射をもう１回行うことによって 、解析平面（Ｄ−Ｚ）内で、計測を実行する、または、複数の変位を計測するた めの前述の技法を洗練する（つまり、観測領域を拡大する）ことも可能である。
【００９９】
それから、複数の変位の計測は、ビーム２０に沿って連続するいくつかの位 置に変換器列９を変位させることによって、いくつかの平行な解析平面（Ｄ−Ｚ ）内で繰り返される。したがって、例えば、すべての平行な平面について約１か ら２秒の合計時間で、乳房の撮像の場合では約３０の平行な平面内で、変位を計 測することも可能になる。
【０１００】
それから、プラットフォーム１５は、予め定められた角度で回転し、それに よって解析平面（Ｄ−Ｚ）も回転させ、新しい方向Ｄに平行ないくつかの解析平 面内で複数の変位を計測する前述の動作が前述のように反復される。したがって 、新しい方向Ｄに沿った成分と方向Ｚに沿った成分とを含む２次元の複数の変位 場
【０１０１】
【数４】
【０１０２】
が、求められる。
【０１０３】
方向Ｄの少なくとも２つの向き（例えば、互いに９０°の２つの向き）、ま たは必要に応じて２つより多い向き（例えば、３６までの角度の向き）を含む複
数の変位の計測後に、粘弾性媒体２内の３次元の変位の場
【０１０４】
【数５】
【０１０５】
を再構成することができる。
【０１０６】
そのためには、例えば、時間の関数として、軸Ｘ、Ｙ、Ｚのそれぞれに沿っ た複数の変位の値が収集され、必要に応じて、さまざまな計測値を時間的に再設 定後（そのために、信号は、周期的な変調を伴って変換器１０によって送信され 、そして、時間の原点は、超音波の発射の各シーケンスの前の周期的な変調にお ける全く同一の点とされる）、これらの値が平均と補間とによって組み合わされ て、３次元の変位ベクトルの１つの場
【０１０７】
【数６】
【０１０８】
となる。
【０１０９】
次に、この変位の場
【０１１０】
【数７】
【０１１１】
（または他の移動パラメータ）は、音響変換器１０にて発生された粘性力学波 の圧縮成分の影響を除去するように補正される。このために、回転演算子が、ベ クトル場
【０１１２】
【数８】
【０１１３】
を求めるために、場
【０１１４】
【数９】
【０１１５】
に適用される。
【０１１６】
したがって、音響変換器１０にて発生した粘性力学波の剪断成分の影響だけ が考慮される。
【０１１７】
この一式のベクトル
【０１１８】
【数１０】
【０１１９】
は、メモリＭまたはマイクロコンピュータ４に保存され、そして、例えば、具 体的にはマイクロコンピュータ４の画面４ａを使用して、ｑの成分の値がグレー レベルまたは色レベルなどの光学的パラメータとして表示されるゆっくりしたフ ィルムの形態で表示することができる。
【０１２０】
したがって、医療用途での健康な組織とガン性の組織など、媒体２内の特性 が異なる領域間の剪断波の伝搬の違いが完全に表示することができる。
【０１２１】
この弾性波の伝搬のフィルムは、さらに、前述の装置１によって同様に生成 可能な通常のエコーグラフの画像にさらに重ねて表示することが可能である。
【０１２２】
さらに、観測されている媒体２の各点の複数の変位だけでなく、媒体２の歪 み、つまり、各変位ベクトルの両成分の複数の空間変数についての微分を成分と するベクトルも計算することもできる。これらの歪みベクトルを、剪断波の伝搬 を明確に表示するために、変位ベクトルのようにフィルムの形態で使用すること が可能で、これらの歪みベクトルは、観測している媒体２に対するプローブ３６ の変位を除外することができるという利点も含む。これらの歪みベクトルは、変 位ベクトルの場合に前述したように、音響変換器１０の粘性力学波の圧縮要素の 影響を取り除くために回転演算子を適用することによって処理することができる 。
【０１２３】
次に、マイクロコンピュータ２７は、媒体２内の剪断波の伝搬を記述してい る少なくとも１つのパラメータが、回転演算子を適用することによって処理され た運動パラメータの変化（変位または歪み速度）から計算される、変位場または 歪み場に基づいてマッピングステップ（ｄ）を実行することが有利である。
【０１２４】
マッピングステップ間に計算される剪断波の伝搬パラメータは、剪断波の速 度ＣT、剪断弾性率μまたはヤング率Ｅ（＝３μ）、複数の剪断波の減衰、剪断 弾性、剪断粘性ζ、または組織の力学緩和時間、または異方性の媒体の場合これ らのパラメータうちのいくつかのパラメータから選択される。
【０１２５】
例えば、観測領域の様々な点で剪断弾性率μと剪断粘性ζを計算することが できる。
【０１２６】
そのために、以下の伝搬方程式（１）を使用することができ（媒体２は、１ 次近似として局所的に一様かつ等方な粘弾性媒体と見なせると仮定する）、複数 の粘性力学波によって媒体２内の各位置
【０１２７】
【数１１】
【０１２８】
で発生する変位
【０１２９】
【数１２】
【０１３０】
は以下の式に従う。
【０１３１】
【数１３】
【０１３２】
ここでρは媒体の密度、μは剪断弾性率（再構成の対象）、λは第２ラメ係数 、
【０１３３】
【数１４】
【０１３４】
は媒体内への拡散を考慮した剪断粘性、ζは圧縮波の粘性である。剪断弾性 率μと第２ラメ係数λとは、ヤング率Ｅとポアソン比σに関して以下の式で与え られる。
【０１３５】
【数１５】
【０１３６】
各ベクトル場
【０１３７】
【数１６】
【０１３８】
は、非回転部分
【０１３９】
【数１７】
【０１４０】
、圧縮不能部分
【０１４１】
【数１８】
【０１４２】
、および非回転で圧縮不能な調和部分
【０１４３】
【数１９】
【０１４４】
とに分解することができる。調和部分
【０１４５】
【数２０】
【０１４６】
は、大局的な複数の変位を考慮し、そのため無視でき、ここでは、この分解を 使用して複数の波の場合を研究しているため、式（１）を複数の剪断波の伝搬を 表している式と複数の圧縮波の伝搬を表している式との２つの式に分離できるこ とを示すことが可能であり、各々の伝搬速度はｃTとｃLとであって（粘性の影 響は無視）、つまり以下のようになる。
【０１４７】
【数２１】
【０１４８】
式（１）に回転演算子を適用することによって、変換器１０にて発生された 弾性力学波の圧縮成分の影響をもはや有していない簡略化された式が得られる。
【０１４９】
【数２２】
【０１５０】
それから、式（４）からパラメータμとζと３次空間微分の現在の値を計算 することができる。式（４）は、例えば、系の転置行列（疑似逆行列）が乗算さ れることによって簡単に解くことができる。式（５）は、２つの未知数（密度ρ は水の密度に等しいとみなされる）に対して実際には６つの方程式（複素数表記 で３つが
【０１５１】
【数２３】
【０１５２】
の実部で３つが虚部）を示していることがわかる。
【０１５３】
したがって、定常状態では、以下の表記のもとで、
【０１５４】
【数２４】
【０１５５】
使用可能な方程式は、以下のように表すことができる。
【０１５６】
【数２５】
【０１５７】
転置行列を乗じることによって、以下の式が得られる。
【０１５８】
【数２６】
【０１５９】
ここでｉ＝ｘ，ｙ，ｚであって、複数の同一の添え字に対するアインシュタイ ン規約が適用されている。定常状態の場合、μとζとの間には相関はなく、その ため互いに独立しているこれらの２つの計算されたパラメータをマッピングする のがより容易になる。
【０１６０】
ただし、局所的に一様で等方性の媒体という近似ができない場合にも、所望 の複数の伝搬パラメータを求めるために同様な計算を使用することができる。こ の場合、より複雑な伝搬方程式が使用される、この伝搬方程式は、以下のように なる。
【０１６１】
【数２７】
【０１６２】
ここでμiは変位ベクトル
【０１６３】
【数２８】
【０１６４】
のｉ番目の成分であり、ρは媒体２の密度であり、λiklmは階数４の弾性テン ソルであり、ηiklmは階数４の粘性テンソルである（複数の同一の添え字に対す るアインシュタイン規約を使用）。
【０１６５】
この式（１′）は、媒体が横方向に等方性を含む、つまり、媒体が一方向に だけ揃えられている繊維を含むと仮定できる場合に簡単になる。この場合、物質 を記述している複数のパラメータを以下の２つのグループに分割することができ る。
【０１６６】
−２つの剪断係数（μ□，μ⊥）を含む１つのグループは複数の繊維の方向 に水平（□）および垂直（⊥）な剪断波の伝搬を記述している。
【０１６７】
−他方のグループ（λ□，λ⊥，λM）は様々な方向の複数の縦波（つまり 複数の圧縮波）の伝搬を記述している
歪みテンソルと応力テンソルとの間の関係は次の式で示される。
【０１６８】
【数２９】
【０１６９】
ここでｕik＝１／２（∂xkｕi＋∂xiuk）は歪みテンソルであってσikは応 力テンソルである。
【０１７０】
変位ベクトル
【０１７１】
【数３０】
【０１７２】
の対応する伝搬方程式は次の式で与えられる。
【０１７３】
【数３１】
【０１７４】
ここでζは剪断粘性である。複数の縦波の粘性（ξ）は、比較的低い周波数（ ＨｚまたはｋＨｚの範囲）では無視できる。さらに、
【０１７５】
【数３２】
【０１７６】
はポアソン比の値のために小さい値を含む。したがって、ξ、ζおよび
【０１７７】
【数３３】
【０１７８】
は値が小さく、そのため式９の右辺の第３項は合理的に無視することができる 。式８を式９に代入すると以下の式を得ることができる。
【０１７９】
【数３４】
【０１８０】
ここで表記τ＝μ□−μ⊥であって、縦波のすべての弾性率は等しい、つまり λ＝λ□＝λ⊥＝λMと仮定している。
【０１８１】
回転演算子の適用後に以下の式が得られる。
【０１８２】
【数３５】
【０１８３】
この式（１１）によって、特にμ⊥、τおよびζを計算するために使用可能 な６つの方程式が実際に得られる。
【０１８４】
ただし、本発明は媒体２の運動パラメータ（変位等）を計測するどのような 方法にも適用可能である。したがってこのパラメータは必要に応じてＭＲＩによ って求められるであろう。
【０１８５】
説明した方法は、例えばＰＣＴ特許出願公開明細書第０４／８５９７８号に 記述されているような媒体２の音響光学撮像と組み合わせることもできる。
【図面の簡単な説明】
【０１８６】
【図１】本発明の一実施形態の複数の剪断波によって撮像する撮像装置の概略 図である。
【図２】図１の装置の一部の詳細な透視図である。
【図３】図１の装置のブロック図である。
【特許請求の範囲】
【請求項１】
剪断成分および圧縮成分を含む力学波が粘弾性媒体内（２）で発生され、か つ、前記力学波の伝搬における様々な点で前記粘弾性媒体（２）の運動を記述す る少なくとも１つのパラメータが求められる撮像方法において、
前記粘弾性媒体（２）内の複数の点での前記力学波の伝搬が、エコーグラフ で観測される少なくとも１つの観測ステップ（ｂ）を含み、
前記力学波の前記圧縮成分による影響を取り除くように、前記運動パラメー タが処理される補正ステップを含むことを特徴とする撮像方法。
【請求項２】
前記補正ステップの間に、前記運動パラメータの回転が計算される、請求項 １に記載の撮像方法。
【請求項３】
前記観測ステップ（ｂ）は、
圧縮超音波の連続するショットが、互いに独立して制御される変換器の列（ ９）によって、少なくとも一秒当たり１０ショットの割合で前記粘弾性媒体（２ ）中に伝送されるサブステップ（ｂ１）と、
前記粘弾性媒体に相互作用する前記複数の圧縮超音波によって生成される複 数のエコーを含む、前記粘弾性媒体（２）から受信した複数の音響信号が、リア ルタイムで検出および記録されるサブステップ（ｂ２）と、を含む、請求項１に 記載の撮像方法。
【請求項４】
前記サブステップ（ｂ２）の間に前記粘弾性媒体から受信した前記複数の連 続する音響信号が、前記力学波の連続する複数の伝搬画像を求めるために処理さ れるステップ（ｃ１）と、
前記粘弾性媒体の前記運動パラメータが、前記連続する複数の伝搬画像から 様々な点で求められるステップ（ｃ２）と、を含む少なくとも１つの処理ステッ プ（ｃ）をさらに含む、請求項３に記載の撮像方法。
【請求項５】
前記粘弾性媒体内の前記力学波の前記連続する複数の伝搬画像は、いくつか の視点に沿って生成され、該連続する画像は、前記運動パラメータのベクトル値 を求めるために組み合わされる、先行する請求項のいずれか１項に記載の撮像方 法。
【請求項６】
前記運動パラメータの３次元ベクトル値が求められる、請求項５に記載の撮 像方法。
【請求項７】
前記いくつかの視点に沿った力学波の連続する複数の伝搬画像を生成するた めに、解析平面（Ｄ−Ｚ）内の伝送／検出方向（Ｄ）に沿って、前記複数の圧縮 超音波を伝送および検出する複数の超音波変換器の列（９）が使用され、前記直 線状の列（９）は、前記粘弾性媒体の外側で変位することで、前記解析平面（Ｄ −Ｚ）を回転させる、請求項５または請求項６に記載の撮像方法。
【請求項８】
前記直線状の列（９）は、前記解析平面（Ｄ−Ｚ）が回転しない並進運動に よって変位し、その後、前記直線状の列（９）は、回転することで、前記解析平 面を回転させる、請求項７に記載の撮像方法。
【請求項９】
前記列の変換器の伝送／検出方向（Ｄ）は、前記解析平面（Ｄ−Ｚ）を並進 も回転もさせずに変化し、その後、前記直線状の列（９）は、並進および回転す ることで、前記解析平面を並進および回転させる、請求項８に記載の撮像方法。
【請求項１０】
前記運動パラメータは、前記粘弾性媒体（２）の変位、変位速度および応力 から選択される、先行する請求項のいずれか１項に記載の撮像方法。
【請求項１１】
前記粘弾性媒体内の伝搬パラメータのマップを求めるために、前記粘弾性媒 体（２）内の少なくとも特定の複数の点での、前記力学波の剪断成分の伝搬を記 述する少なくとも１つのパラメータが、前記運動パラメータの時間の経過に伴う 変化から計算されるマッピングステップ（ｄ）をさらに含む、先行する請求項の いずれか１項に記載の撮像方法。
【請求項１２】
前記マッピングステップ中に計算される複数の剪断波の伝搬パラメータは、 前記複数の剪断波の速度、剪断弾性率、ヤング率、前記複数の剪断波の減衰、剪 断弾性、剪断粘性および力学緩和時間、または、異方性媒体の場合には、これら のパラメータうちのいくつかのパラメータ、から選択される、請求項１１に記載 の撮像方法。
【請求項１３】
剪断成分と圧縮成分とを含む力学波を粘弾性媒体内（２）で発生させる手段 と、前記力学波の伝搬中に様々な点で前記粘弾性媒体の運動を記述する少なくと も１つのパラメータを求める手段と、を含む撮像装置において、
前記粘弾性媒体（２）内の複数の点での前記力学波の伝搬をエコーグラフで 観測するために制御された複数の超音波変換器の列（９）を含み、
前記力学波の前記圧縮成分による影響を取り除くように前記運動パラメータ を処理する補正手段（２７）と、を含むことを特徴とする撮像装置。
【請求項１４】
前記補正手段（２７）は、前記運動パラメータの回転を計算するように構成 される、請求項１３に記載の撮像装置。
【請求項１５】
前記複数の変換器の列（９）によって、少なくとも１秒あたり１０ショット の割合で、前記複数の圧縮超音波の連続するショットを前記粘弾性媒体（２）内 に伝送し（ｂ１）、
前記粘弾性媒体の複数の反射粒子に相互作用する前記複数の圧縮超音波によ って発生した複数のエコーを含む、前記粘弾性媒体（２）から受信した複数の音 響信号をリアルタイムで検出および記録する（ｂ２）、複数の手段（２７、３４ ）を含む請求項１３に記載の撮像装置。
【請求項１６】
前記力学波の連続する複数の伝搬画像を求めるために、前記粘弾性媒体から 受信した、前記複数の連続する音響信号を処理し（ｃ１）、
前記粘弾性媒体の前記運動パラメータを、前記連続する複数の伝搬画像から 、観測場内の様々な点で求める（ｃ２）ように構成される、複数の処理手段（３ ４、３８、３７）を含む請求項１５に記載の撮像装置。
【請求項１７】
前記複数の処理手段（３４、３８、３７）は、いくつかの視点に沿って、前 記粘弾性媒体（２）内の力学波のいくつかの連続する伝搬画像を生成し、そして 、前記連続する画像を、前記運動パラメータのベクトル値を求めるために組み合 わせるように構成される、請求項１６に記載の撮像装置。
【請求項１８】
前記複数の処理手段は、前記運動パラメータの３次元ベクトル値を求めるよ うに構成されている、請求項１７に記載の撮像装置。
【請求項１９】
前記複数の変換器の列（９）は、前記いくつかの視点に沿った力学波の連続 する複数の伝搬画像を生成するために、解析平面（Ｄ−Ｚ）内の伝送／検出方向 （Ｄ）に沿って前記複数の圧縮超音波を伝送および検出し、前記装置は、前記粘 弾性媒体（２）の外側で、前記直線状の列（９）を変位させる複数の変位手段（ ３０、３１、３２、２２、３３）を含み、前記複数の手段は前記解析平面（Ｄ− Ｚ）を回転させるように構成されている、請求項１７または請求項１８に記載の 撮像装置。
【請求項２０】
前記複数の変位手段（３０、３１、３２、３３）を制御するように構成され る制御手段（２７）を含み、前記制御手段は、前記解析平面（Ｄ−Ｚ）を回転さ せずに前記列（９）を並進運動によって変位させ、その後、前記列（９）を回転 させることで、前記解析平面を回転させる、請求項１９に記載の撮像装置。
【請求項２１】
前記複数の変位手段（３０、３１、３２、２２、３３）を制御するように構 成されている制御手段（２７）を含み、前記制御手段（２７）は、前記解析平面 （Ｄ−Ｚ）を並進も回転もさせずに、前記複数の変換器の列（９）の前記伝送／ 検出方向（Ｄ）を変化させ、その後、前記直線状の列（９）を並進および回転さ せることで、前記解析平面（Ｄ−Ｚ）を並進および回転させるように構成されて いる、請求項１９または請求項２０に記載の撮像装置。
【請求項２２】
前記粘弾性媒体内の前記伝搬パラメータのマップを求めるために、前記粘弾 性媒体（２）内の少なくとも特定の複数の点での前記力学波の前記剪断成分の伝 搬を記述する少なくとも１つのパラメータを計算するように構成されているマッ ピング手段（２７）を含む、請求項１９から請求項２１のいずれか１項に記載の 撮像装置。
【請求項１】
剪断成分および圧縮成分を含む力学波が粘弾性媒体内（２）で発生され、か つ、前記力学波の伝搬における様々な点で前記粘弾性媒体（２）の運動を記述す る少なくとも１つのパラメータが求められる撮像方法において、
前記粘弾性媒体（２）内の複数の点での前記力学波の伝搬が、エコーグラフ で観測される少なくとも１つの観測ステップ（ｂ）を含み、
前記力学波の前記圧縮成分による影響を取り除くように、前記運動パラメー タが処理される補正ステップを含むことを特徴とする撮像方法。
【請求項２】
前記補正ステップの間に、前記運動パラメータの回転が計算される、請求項 １に記載の撮像方法。
【請求項３】
前記観測ステップ（ｂ）は、
圧縮超音波の連続するショットが、互いに独立して制御される変換器の列（ ９）によって、少なくとも一秒当たり１０ショットの割合で前記粘弾性媒体（２ ）中に伝送されるサブステップ（ｂ１）と、
前記粘弾性媒体に相互作用する前記複数の圧縮超音波によって生成される複 数のエコーを含む、前記粘弾性媒体（２）から受信した複数の音響信号が、リア ルタイムで検出および記録されるサブステップ（ｂ２）と、を含む、請求項１に 記載の撮像方法。
【請求項４】
前記サブステップ（ｂ２）の間に前記粘弾性媒体から受信した前記複数の連 続する音響信号が、前記力学波の連続する複数の伝搬画像を求めるために処理さ れるステップ（ｃ１）と、
前記粘弾性媒体の前記運動パラメータが、前記連続する複数の伝搬画像から 様々な点で求められるステップ（ｃ２）と、を含む少なくとも１つの処理ステッ プ（ｃ）をさらに含む、請求項３に記載の撮像方法。
【請求項５】
前記粘弾性媒体内の前記力学波の前記連続する複数の伝搬画像は、いくつか の視点に沿って生成され、該連続する画像は、前記運動パラメータのベクトル値 を求めるために組み合わされる、先行する請求項のいずれか１項に記載の撮像方 法。
【請求項６】
前記運動パラメータの３次元ベクトル値が求められる、請求項５に記載の撮 像方法。
【請求項７】
前記いくつかの視点に沿った力学波の連続する複数の伝搬画像を生成するた めに、解析平面（Ｄ−Ｚ）内の伝送／検出方向（Ｄ）に沿って、前記複数の圧縮 超音波を伝送および検出する複数の超音波変換器の列（９）が使用され、前記直 線状の列（９）は、前記粘弾性媒体の外側で変位することで、前記解析平面（Ｄ −Ｚ）を回転させる、請求項５または請求項６に記載の撮像方法。
【請求項８】
前記直線状の列（９）は、前記解析平面（Ｄ−Ｚ）が回転しない並進運動に よって変位し、その後、前記直線状の列（９）は、回転することで、前記解析平 面を回転させる、請求項７に記載の撮像方法。
【請求項９】
前記列の変換器の伝送／検出方向（Ｄ）は、前記解析平面（Ｄ−Ｚ）を並進 も回転もさせずに変化し、その後、前記直線状の列（９）は、並進および回転す ることで、前記解析平面を並進および回転させる、請求項８に記載の撮像方法。
【請求項１０】
前記運動パラメータは、前記粘弾性媒体（２）の変位、変位速度および応力 から選択される、先行する請求項のいずれか１項に記載の撮像方法。
【請求項１１】
前記粘弾性媒体内の伝搬パラメータのマップを求めるために、前記粘弾性媒 体（２）内の少なくとも特定の複数の点での、前記力学波の剪断成分の伝搬を記 述する少なくとも１つのパラメータが、前記運動パラメータの時間の経過に伴う 変化から計算されるマッピングステップ（ｄ）をさらに含む、先行する請求項の いずれか１項に記載の撮像方法。
【請求項１２】
前記マッピングステップ中に計算される複数の剪断波の伝搬パラメータは、 前記複数の剪断波の速度、剪断弾性率、ヤング率、前記複数の剪断波の減衰、剪 断弾性、剪断粘性および力学緩和時間、または、異方性媒体の場合には、これら のパラメータうちのいくつかのパラメータ、から選択される、請求項１１に記載 の撮像方法。
【請求項１３】
剪断成分と圧縮成分とを含む力学波を粘弾性媒体内（２）で発生させる手段 と、前記力学波の伝搬中に様々な点で前記粘弾性媒体の運動を記述する少なくと も１つのパラメータを求める手段と、を含む撮像装置において、
前記粘弾性媒体（２）内の複数の点での前記力学波の伝搬をエコーグラフで 観測するために制御された複数の超音波変換器の列（９）を含み、
前記力学波の前記圧縮成分による影響を取り除くように前記運動パラメータ を処理する補正手段（２７）と、を含むことを特徴とする撮像装置。
【請求項１４】
前記補正手段（２７）は、前記運動パラメータの回転を計算するように構成 される、請求項１３に記載の撮像装置。
【請求項１５】
前記複数の変換器の列（９）によって、少なくとも１秒あたり１０ショット の割合で、前記複数の圧縮超音波の連続するショットを前記粘弾性媒体（２）内 に伝送し（ｂ１）、
前記粘弾性媒体の複数の反射粒子に相互作用する前記複数の圧縮超音波によ って発生した複数のエコーを含む、前記粘弾性媒体（２）から受信した複数の音 響信号をリアルタイムで検出および記録する（ｂ２）、複数の手段（２７、３４ ）を含む請求項１３に記載の撮像装置。
【請求項１６】
前記力学波の連続する複数の伝搬画像を求めるために、前記粘弾性媒体から 受信した、前記複数の連続する音響信号を処理し（ｃ１）、
前記粘弾性媒体の前記運動パラメータを、前記連続する複数の伝搬画像から 、観測場内の様々な点で求める（ｃ２）ように構成される、複数の処理手段（３ ４、３８、３７）を含む請求項１５に記載の撮像装置。
【請求項１７】
前記複数の処理手段（３４、３８、３７）は、いくつかの視点に沿って、前 記粘弾性媒体（２）内の力学波のいくつかの連続する伝搬画像を生成し、そして 、前記連続する画像を、前記運動パラメータのベクトル値を求めるために組み合 わせるように構成される、請求項１６に記載の撮像装置。
【請求項１８】
前記複数の処理手段は、前記運動パラメータの３次元ベクトル値を求めるよ うに構成されている、請求項１７に記載の撮像装置。
【請求項１９】
前記複数の変換器の列（９）は、前記いくつかの視点に沿った力学波の連続 する複数の伝搬画像を生成するために、解析平面（Ｄ−Ｚ）内の伝送／検出方向 （Ｄ）に沿って前記複数の圧縮超音波を伝送および検出し、前記装置は、前記粘 弾性媒体（２）の外側で、前記直線状の列（９）を変位させる複数の変位手段（ ３０、３１、３２、２２、３３）を含み、前記複数の手段は前記解析平面（Ｄ− Ｚ）を回転させるように構成されている、請求項１７または請求項１８に記載の 撮像装置。
【請求項２０】
前記複数の変位手段（３０、３１、３２、３３）を制御するように構成され る制御手段（２７）を含み、前記制御手段は、前記解析平面（Ｄ−Ｚ）を回転さ せずに前記列（９）を並進運動によって変位させ、その後、前記列（９）を回転 させることで、前記解析平面を回転させる、請求項１９に記載の撮像装置。
【請求項２１】
前記複数の変位手段（３０、３１、３２、２２、３３）を制御するように構 成されている制御手段（２７）を含み、前記制御手段（２７）は、前記解析平面 （Ｄ−Ｚ）を並進も回転もさせずに、前記複数の変換器の列（９）の前記伝送／ 検出方向（Ｄ）を変化させ、その後、前記直線状の列（９）を並進および回転さ せることで、前記解析平面（Ｄ−Ｚ）を並進および回転させるように構成されて いる、請求項１９または請求項２０に記載の撮像装置。
【請求項２２】
前記粘弾性媒体内の前記伝搬パラメータのマップを求めるために、前記粘弾 性媒体（２）内の少なくとも特定の複数の点での前記力学波の前記剪断成分の伝 搬を記述する少なくとも１つのパラメータを計算するように構成されているマッ ピング手段（２７）を含む、請求項１９から請求項２１のいずれか１項に記載の 撮像装置。
【図１】
【図２】
【図３】
【図２】
【図３】
【公表番号】特表２００８−５３４１９８（Ｐ２００８−５３４１９８Ａ）
【公表日】平成２０年８月２８日（２００８．８．２８）
【国際特許分類】
【出願番号】特願２００８−５０４７９６（Ｐ２００８−５０４７９６）
【出願日】平成１８年３月３０日（２００６．３．３０）
【国際出願番号】ＰＣＴ／ＦＲ２００６／０００７０２
【国際公開番号】ＷＯ２００６／１０６２１３
【国際公開日】平成１８年１０月１２日（２００６．１０．１２）
【出願人】（５０２１５５３５６）サントル ナショナル デ ラ ルシェルシュ シィアンティフィク （セ．エヌ．エール．エス．） (10)
【氏名又は名称原語表記】ＣＥＮＴＲＥ ＮＡＴＩＯＮＡＬ ＤＥ ＬＡ ＲＥＣＨＥＲＣＨＥ ＳＣＩＥＮＴＩＦＩＱＵＥ （Ｃ．Ｎ．Ｒ．Ｓ．）
【出願人】（５０７３３２７６７）ユニヴェルシテ パリ ７ − デニ デデロ (1)
【氏名又は名称原語表記】ＵＮＩＶＥＲＳＩＴＥ ＰＡＲＩＳ ７ − ＤＥＮＩＳ ＤＩＤＥＲＯＴ
【Ｆターム（参考）】
【公表日】平成２０年８月２８日（２００８．８．２８）
【国際特許分類】
【出願日】平成１８年３月３０日（２００６．３．３０）
【国際出願番号】ＰＣＴ／ＦＲ２００６／０００７０２
【国際公開番号】ＷＯ２００６／１０６２１３
【国際公開日】平成１８年１０月１２日（２００６．１０．１２）
【出願人】（５０２１５５３５６）サントル ナショナル デ ラ ルシェルシュ シィアンティフィク （セ．エヌ．エール．エス．） (10)
【氏名又は名称原語表記】ＣＥＮＴＲＥ ＮＡＴＩＯＮＡＬ ＤＥ ＬＡ ＲＥＣＨＥＲＣＨＥ ＳＣＩＥＮＴＩＦＩＱＵＥ （Ｃ．Ｎ．Ｒ．Ｓ．）
【出願人】（５０７３３２７６７）ユニヴェルシテ パリ ７ − デニ デデロ (1)
【氏名又は名称原語表記】ＵＮＩＶＥＲＳＩＴＥ ＰＡＲＩＳ ７ − ＤＥＮＩＳ ＤＩＤＥＲＯＴ
【Ｆターム（参考）】
[ Back to top ]