(57)【要約】  超音波処理装置であって、集束超音波を生成するように構成された超音波治療トランスデューサ（ＡＴＡ）と、超音波治療トランスデューサに接続された超音波撮像トランスデューサ（ＵＩＤ）と、電子システムとを備え、電子システムは、キャビテーション気泡（ＢＣ）クラウドを生成する超音波パルス列を発するように超音波治療トランスデューサを制御することと、処置領域の少なくとも１つの像を取得するように超音波撮像トランスデューサを制御することと、超音波撮像トランスデューサにより捕捉された、超音波治療トランスデューサにより発せられた超音波パルスの複数のエコー信号を取得し、処理して、キャビテーション気泡クラウドの像を再構築することであって、処理は時空間フィルタリングを含む、再構築することと、処置領域の上記像に重ねてキャビテーション気泡クラウドの上記像を表示することとを行うように構成された、超音波処理装置。
2021523776.tif 000002
【特許請求の範囲】
【請求項１】
超音波を用いて治療を提供する装置であって、
集束超音波の生成に適した治療超音波トランスデューサ（ＡＴＡ）と、
前記治療超音波トランスデューサに関連する撮像超音波トランスデューサ（ＵＩＤ）と、
電子システムと、
を備え、前記電子システムは、
焦点スポット（ＴＦ）が人間又は動物の体の処置領域内部に位置しているとき、前記トランスデューサの前記焦点スポットにキャビテーション気泡（ＢＣ）クラウドを生成させるのに適した持続時間でエネルギーの超音波パルス列を発するように前記治療超音波トランスデューサを制御することと、
前記治療超音波トランスデューサにより発せされる２つの超音波パルス間で、前記処置領域に向けられた超音波を発して、前記超音波のエコーを取得し、前記エコーを処理して、前記処置領域の少なくとも１つの像を再構築するように前記撮像超音波トランスデューサを制御することと、
前記治療超音波トランスデューサにより発せられたＮ＞１個の超音波パルスの複数のエコー信号を取得することであって、前記信号は前記撮像超音波トランスデューサにより捕捉される、取得することと、
ビーム形成アルゴリズムを使用して前記治療超音波トランスデューサにより発せられた前記Ｎ個の超音波パルスの前記エコー信号を処理し、各エコー像（Ｆ１〜ＦＮ）を形成し、時空間フィルタリングを使用して、前記キャビテーション気泡からの前記超音波パルスの後方散乱を表す前記エコー像の成分を抽出できるようにし、前記処置領域の組織からの後方散乱を表す成分から分離して、前記キャビテーション気泡クラウドの像（ＩＢＣ）を再構築することと、
前記処置領域の前記像に重ねて、前記キャビテーション気泡クラウドの前記像を表示することと、
を行うように構成される、超音波により治療を提供する装置。
【請求項２】
前記電子システムは、前記キャビテーション気泡クラウドの前記像の前記再構築中、前記治療超音波トランスデューサと前記焦点スポットとの間及び前記焦点スポットと前記撮像超音波トランスデューサとの間の前記超音波の伝播時間差の補償を導入するようにも構成される、請求項１に記載の超音波により治療を提供する装置。
【請求項３】
前記電子システムは、０．１μｓ〜５０μｓ、好ましくは０．５μｓ〜２０μｓの持続時間を有する超音波パルス列を発するよう前記治療超音波トランスデューサを制御するように構成される、請求項１又は２に記載の超音波により治療を提供する装置。
【請求項４】
Ｎは２〜１０，０００、好ましくは２〜１，０００を含む、請求項１〜３の何れか１項に記載の超音波により治療を提供する装置。
【請求項５】
前記電子システムは、並列ビーム形成アルゴリズムにより前記キャビテーション気泡クラウドの前記像を再構築するように構成される、請求項１〜４の何れか１項に記載の超音波により治療を提供する装置。
【請求項６】
前記電子システム及び前記治療超音波トランスデューサは、持続時間０．１μｓ〜５０μｓ、中心周波数１００ｋＨｚ〜５ＭＨｚ、反復率１〜１０００Ｈｚで超音波パルスを発するように構成される、請求項１〜５の何れか１項に記載の超音波により治療を提供する装置。
【請求項７】
前記電子システム及び前記撮像超音波トランスデューサは、中心周波数２〜１５ＭＨｚを有する超音波パルスを発するように構成される、請求項１〜６の何れか１項に記載の超音波により治療を提供する装置。
【請求項８】
前記電子システムは、特異値分解を使用して前記時空間フィルタリングを実行するように構成される、請求項１〜７の何れか１項に記載の超音波により治療を提供する装置。
【請求項９】
前記治療超音波トランスデューサ及び前記撮像超音波トランスデューサは、同軸に配置される、請求項１〜８の何れか１項に記載の超音波により治療を提供する装置。
【請求項１０】
前記治療超音波トランスデューサにより発せられた前記超音波パルスの前記焦点スポットを前記処置領域において移動させる手段（ＧＦＦ、ＣＰＭ）も備える、請求項１〜９の何れか１項に記載の超音波により治療を提供する装置。
【請求項１１】
前記電子システムは、前記キャビテーション気泡クラウドの前記像に応じて前記治療超音波トランスデューサにより発せられる前記超音波パルスの電力レベルを調整するようにも構成される、請求項１〜１０の何れか１項に記載の超音波により治療を提供する装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、人間又は動物の体の超音波による治療を提供する装置に関する。本発明は特に、動きのある臓器、すなわち、典型的には腹部臓器及び胸部臓器の高強度集束超音波（ＨＩＦＵ）、組織破砕、結石破砕、血栓破砕等による治療に適用可能である。
【背景技術】
【０００２】
組織破砕は、キャビテーション気泡クラウド（簡潔にするために「キャビテーションクラウド」とも呼ばれる）を生成する集束超音波パルスにより組織を機械的に破砕する技法である。種々の病理の治療への組織破砕の使用は、高リスク手術にとって代わることができる有望な手法である。キャビテーションは、高度に制御される焦点領域における望ましくない組織を非侵襲的に熱的影響なしで破壊できるようにする。しかしながら、人体は非常に不均質な媒質であり、治療が、直接的な視覚的フィードバックなしで体内で行われることは、キャビテーション領域を精密にモニタすることを極めて重要なものにする。特に、心臓用途では、胸郭が超音波の治療パスに大きな収差を生じさせる恐れがある。
【０００３】
同様の考慮事項が、結石破砕（腎臓結石の破砕）及び血栓破砕（血栓の破砕）等の概念的に組織破砕と同様の技法、更には高強度集束超音波による治療の分野にも該当する。
【０００４】
後方散乱超音波エコーの周波数の高調波の解析により、キャビテーションを受動的に検出することが可能であるが、気泡クラウドを精密に位置特定することはできない。
【０００５】
キャビテーション効果及び治療中の解剖学的構造をリアルタイムで表示することができる従来のＢモード超音波撮像技法を使用することが可能である。しかしながら、この手法は全体的に満足のいくものではない。特に、気泡は不均質な生体組織内で生成されるため、気泡のエコーは、組織からの全てのエコーの中から検出することが難しいことがある。したがって、気泡の識別は非常に主観的で非定量的なままであり、キャビテーションクラウドの輪郭の画定には問題が多い。
【０００６】
記事［１］は、時空間フィルタリングに関連する並列（又は「超高速」）超音波撮像技法を利用することにより、キャビテーション気泡のよりよい区別を得ることができることを教示している。しかしながら、この手法は、治療される組織が常に動く場合、例えば、心臓、肝臓、又は腎臓の場合等にあまり適さないことが証明されている。特に、撮像に使用される超音波は、関心のある全ての領域（又は関心のある領域よりも大きな領域）に適用され、エコー信号は、キャビテーションに起因するもののみではなく、この領域で生じるあらゆる時空間変動の影響を非常に受けやすい。これは、組織が静的であるか、又は準静的である場合、許容可能であるが、動いている組織が存在する場合、時空間フィルタリングは、キャビテーション気泡からのエコー信号を識別するはるかに効率が劣る方法になる。しかしながら、キャビテーションクラウド及び組織に対する位置を精密に見ることが可能なことが特に重要なのはまさしく後者の場合である。
【０００７】
受動的な検出方法は、キャビテーション気泡の位置推定にも使用されてもいる。例えば：
−［２］は、熱アブレーション治療トランスデューサにより生成された気泡を音響的にマッピングする受動的方法を記載している。この方法では、連続波放射線又は長持続時間（数百回の振動）の波列放射線により生成された気泡をマッピングすることができる。波列の持続時間により、気泡クラウドの位置を精密に識別することは可能ではない。［２］は、使用される波長と均等な解像度を有する後方散乱エネルギーのマップを計算する方法を提案しているが、これは、クラウドの境界を精密に識別することができない。
−［３］は、連続放射線をモニタするという同じ問題を解決する［２］に対する代替の方法を記載している。
−［４］は、特定の変更なしでの組織破砕に向けた方法［２］の適合を記載している。
−［５］は、時間反転を超音波パルスに適用する方法を記載している。引用文献［２］、［３］、及び［４］と同様に、この方法も長パルスを使用し、キャビテーション気泡クラウドの境界を精密に識別することができない。
【発明の概要】
【発明が解決しようとする課題】
【０００８】
本発明は、従来技術の上記欠点を解消することを目的とし、より詳細には、移動する組織の治療中を含め、キャビテーションクラウドを精密に位置特定しセグメント化することができる撮像技法を提供することを目的とする。
【課題を解決するための手段】
【０００９】
本発明によれば、この目的は、組織の能動的な超音波撮像、同期受動キャビテーション撮像、及び時空間フィルタリング（例えば、特異値分解）を関連付けることにより達成される。同期受動撮像は、治療パルスの放射シーケンスと同期された、媒質との治療超音波ビームの相互作用のエコーの再構築に基づく。組織の動きにロバストであり、したがって、撮像される組織が動いている場合であっても良好な性能を得ることができる。
【００１０】
本発明の技法は、撮像に使用される超音波パルス、すなわち、超音波治療パルスの持続時間が低減するにつれて増大する空間解像度でキャビテーション気泡クラウドを位置特定することができる。したがって、好ましくは、これらのパルスは０．１μｓ〜５０μｓの持続時間を有する（これは、何れの場合でも、集束治療パルスに望ましい）。引用文献［２］〜［５］に記載された撮像方法は、短パルスの使用から何も得るものがない。
【００１１】
時空間フィルタリングを受動キャビテーション撮像に適用することは、文献［６］において既に示唆されていたが、超音波により血液脳関門を開くという完全に異なる状況においてであった。この技法は、注入微小気泡及びキャビテーション気泡の生成には不十分な強度の超音波を使用する。さらに、エコーは頭蓋骨を通して取得され、それにより、治療される組織の撮像に能動超音波技法を使用することが略不可能になる。
【００１２】
したがって、本発明の一趣旨は、超音波を用いて治療を提供する装置であって、
−集束超音波の生成に適した治療超音波トランスデューサと、
−治療超音波トランスデューサに関連する撮像超音波トランスデューサと、
−電子システムと、
を備え、電子システムは、
−焦点スポットが人間又は動物の体の処置領域内部に位置しているとき、トランスデューサの上記焦点スポットにキャビテーション気泡クラウドを生成させるのに適した持続時間でエネルギーの超音波パルス列を発するように超音波トランスデューサを制御することと、
−治療超音波トランスデューサにより発せられる２つの超音波パルス間で、処置領域に向けられた超音波を発して、上記超音波のエコーを取得し、エコーを処理して、処置領域の少なくとも１つの像を再構築するように撮像超音波トランスデューサを制御することと、
−治療超音波トランスデューサにより発せられたＮ＞１個の超音波パルスの複数のエコー信号を取得することであって、上記信号は撮像超音波トランスデューサにより捕捉される、取得することと、
−ビーム形成アルゴリズムを使用して、取得された信号を処理して、各エコー像を形成し、時空間フィルタリングを使用して、キャビテーション気泡からの超音波パルスの後方散乱を表す上記エコー像の成分を抽出できるようにし、処置領域の組織からの後方散乱を表す成分からそれらを分離して、キャビテーション気泡クラウドの像を再構築する（これらの動作はこの順序で又は逆順で実行し得る）ことと、
−処置領域の上記像に重ねて、キャビテーション気泡クラウドの像を表示することと、
を行うように構成される、超音波により治療を提供する装置である。
【００１３】
そのような装置の特定の実施形態によれば、
−上記電子システムは、キャビテーション気泡クラウドの像の再構築中、治療超音波トランスデューサと焦点スポットとの間及び焦点スポットと撮像超音波トランスデューサとの間の超音波の伝播時間差の補償を導入するように構成することもできる。
−上記電子システムは、０．１μｓ〜５０μｓ、好ましくは０．５μｓ〜２０μｓの持続時間を有する超音波パルス列を発するよう治療超音波トランスデューサを制御するように構成することもできる。
−Ｎは２〜１０，０００、好ましくは２〜１，０００を含み得る。
−上記電子システムは、並列ビーム形成アルゴリズムによりキャビテーション気泡クラウドの像を再構築するように構成することもできる。
−上記電子システム及び治療超音波トランスデューサは、持続時間０．１μｓ〜５０μｓ、中心周波数１００ｋＨｚ〜５ＭＨｚ、反復率１〜１０００Ｈｚで超音波パルスを発するように構成し得る。
−上記電子システム及び撮像超音波トランスデューサは、中心周波数２〜１５ＭＨｚを有する超音波パルスを発するように構成し得る。
−上記電子システムは、特異値分解を使用して上記時空間フィルタリングを実行するように構成することもできる。
−治療超音波トランスデューサ及び撮像超音波トランスデューサは、同軸に配置し得る。
−本装置は、治療超音波トランスデューサにより発せられた超音波パルスの焦点スポットを処置領域において移動させる手段を備えることもできる。
−上記電子システムは、キャビテーション気泡クラウドの像に応じて治療超音波トランスデューサにより発せられる超音波パルスの電力レベルを調整するように構成することもできる。
【００１４】
本発明の他の特徴、詳細、及び利点は、例として与えられる添付図面を参照して与えられる説明を読むことで明らかになろう。
【図面の簡単な説明】
【００１５】
【図１】本発明の一実施形態による装置である。
【図２】そのような装置の機能概略図である。
【図３】本発明の一実施形態による受動撮像シーケンスである。
【図４】特異値分解アルゴリズムの図である。
【図５】本発明の一実施形態による装置を使用して取得されたキャビテーション像である。
【図６】本発明の技術的効果を示すグラフである。
【発明を実施するための形態】
【００１６】
図１の装置は、環形状の１組の、好ましくは５〜２０の数の基本トランスデューサで構成された治療トランスデューサＡＴＡの中心に配置された撮像超音波トランスデューサＵＩＤを備え、上記基本トランスデューサは、超音波を集束する球形ボウルＳＢにおいて位置合わせされ同心に配置される。撮像トランスデューサＵＩＤは、治療トランスデューサの最も内側のリングの中心に配置される。
【００１７】
トランスデューサＵＩＤは、周期性又は非周期性アレイに一次元（通常、線形）又は二次元に配置された複数の超音波検出要素を備える。それは、例えば、６４個の要素のアレイを含む二次元超音波検査プローブの問題であり得る。
【００１８】
撮像トランスデューサＵＩＤは足場Ｆに留められ、Ｏｘ軸（参照符号ＰＪはピボットジョイントを示す）の回りを相対的に回転運動できるようにする機械的リンクＭＬにより、治療トランスデューサに接続された機械的アーム（図示せず）に接続できるようにする。したがって、治療トランスデューサの向き（Ｏｚ’軸）を変更しながら−電気モータ（図示せず）により−、撮像トランスデューサの向き（Ｏｚ軸に沿った）を一定に保つことが可能であり、それにより、治療トランスデューサの焦点スポットの、Ｏｚ軸に直交する移動を可能にする。さらに、種々の基本トランスデューサの制御信号間のオフセットを細かく制御することにより、治療トランスデューサの焦点距離、ひいてはＯｚ’軸に沿った焦点スポットの位置を変更することが可能である。これにより、機械的手段及び電子的手段のハイブリッドにより処置領域の二次元スキャンを実行することができる。
【００１９】
変更も可能である、例えば、機械的リンクは、Ｏｘ及びＯｚの両方に直交するＯｙ軸である第２の軸の回りで治療トランスデューサが回転できるようにすることもできる。
【００２０】
図１の組立体は、欧州特許出願公開第３，２３６，４６７Ａ１号明細書により詳細に記載されている。多くの変形により本発明の実施が可能であるため、これは単なる例として与えられている。例えば、治療トランスデューサの焦点スポットによる処置領域のスキャンは、複数の機械的手段を介して達成することもでき、又は逆に、純粋に電子手段を介して（基本トランスデューサの行列アレイの使用により）達成することもできる。さらに、撮像トランスデューサが治療トランスデューサと同軸であることも重要ではなく、また治療トランスデューサに機械的にリンクされることも重要ではなく、撮像トランスデューサ及び治療トランスデューサの相対位置が制御可能であり、超音波治療パルスの焦点スポットが撮像トランスデューサの観測領域内部に配置されることで十分である。
【００２１】
さらに、撮像トランスデューサは、異なる数の要素を含んでもよく、又は２平面若しくは行列アレイタイプであってもよい。
【００２２】
何れの場合でも、電子システムＳＥＬは、
−キャビテーション気泡クラウドの形成の誘導に適した集束超音波パルスを生成し、
−処置領域をスキャンするようにこれらの超音波の焦点スポットを移動させ、
−処置領域及びキャビテーション気泡クラウドの超音波検査像を取得、表示し、画面Ｅに表示する
ように撮像トランスデューサ及び治療トランスデューサを制御するのに提供されなければならない。
【００２３】
電子システムＳＥＬは、メモリに記憶されたプログラムを実行する１つ又は複数のプロセッサ及びこの又はこれらのプロセッサの制御下で動作するアナログ及び／又はデジタル電子回路の両方を備える。図２は機能概略図を示す。この概略図では、参照符号ＰＴＤはデータプロセッサを示し、一般に、マイクロプロセッサ若しくはマイクロコントローラを構成する基板、コンピュータ、又はより複雑な組のプログラマブルデジタル電子回路の問題であり得る。このデータプロセッサは、インターフェースデバイスＩＣ（キーボード、コンピュータマウス等）を介してコマンド及び／又はパラメータをユーザ又はオペレータから受信し、関数生成器ＧＦ、撮像トランスデューサを制御する回路ＣＰＴＩ、ビーム形成回路ＧＦＦ、及び電気モータを制御する回路ＣＰＭの動作を制御し、撮像トランスデューサの取得回路から出力された信号を処理して像を再構築し、像は画面Ｅに表示される。
【００２４】
関数生成器ＧＦは、電子システムの他の機能ユニットに送られる複数の電子信号を生成する。第１の信号Ｓ１、例えば、周波数１ＭＨｚの正弦波信号はビーム形成回路ＧＦＦに送られ、ビーム形成回路ＧＦＦは、種々の位相シフトを有する同じ周波数の複数の個々の信号に分割し、これらの種々の信号は電力増幅器ＡＰにより増幅されてから、撮像トランスデューサの基本トランスデューサに適用される。先に引用した欧州特許出願公開第３，２３６，４６７Ａ１号明細書に詳細に説明されるように、ビーム形成回路ＧＦＦにより導入される位相シフトは、治療トランスデューサの焦点距離を変更させ、ひいては、超音波の焦点スポットＴＦの位置をＯｚ’軸に沿って変更させることができる。関数生成器は、はるかに低い周波数、例えば１００Ｈｚの方形波信号Ｓ２も生成し、これは電力増幅器をアクティブ化及び非アクティブ化し、このようにして、治療トランスデューサは、反復率１００Ｈｚで超音波パルス（「治療パルス」）を発し、パルスは、例えば持続時間８μｓを有する。
【００２５】
関数生成器ＧＦは、高周波数、例えば２ＭＨｚの第３の信号Ｓ３も生成し、これは信号Ｓ２と共に、撮像トランスデューサを制御する回路ＣＰＴＩに送られる。この回路は、治療トランスデューサにより発せられるより強力なパルス間の間隔で低強度超音波パルスを発する目的で撮像トランスデューサを制御する。これらのパルスの伝播は観測領域ＲＯを定義し、観測領域ＲＯは、二次元撮像トランスデューサの場合、高さがＯｚ軸と一致する台形の形状を有する。
【００２６】
撮像トランスデューサの取得回路ＣＡＴＩは、撮像トランスデューサＵＩＤにより検出されるエコー信号を取得し、エコー信号をデジタルフォーマットに変換してプロセッサＰＴＤに送信し、プロセッサＰＴＤはそれらを処理する。超音波治療パルスの放射に続き、持続時間（例えば、２５０μｓ）が焦点スポットＴＦの最大深度に依存する第１の時間窓において、撮像トランスデューサは治療パルスのエコーを検出し、したがって、受動モードで動作し、以下に詳細に説明するように、これらのエコー信号により、プロセッサＰＴＤはキャビテーション気泡クラウドの像を再構築することができる。第１の窓の終わりから続く治療パルスの放射までの第２の時間窓において、撮像トランスデューサは、それ自体が発したパルスのエコーを検出し、これにより、プロセッサＰＴＤは、従来の能動的超音波撮像技法を使用して処置領域の解剖学的構造の像を再構築することができる。
【００２７】
能動的超音波撮像に使用されるパルスの場合、同じトランスデューサが超音波のソース及び検出器として機能する。そのようなことは、超音波パルスが治療トランスデューサにより発せられ、そのエコーが撮像トランスデューサにより検出される受動的撮像では該当しない。したがって、受動モードで取得されるキャビテーション気泡クラウドの像を再構築できるようにするには、「行き」移動時間（治療トランスデューサから焦点スポットまで）と「帰り」移動時間（焦点スポットから撮像トランスデューサまで）との間の差を知る必要がある。装置の機械的構成が十分な精度で知られている場合、これは計算により取得し得るが、一般に、較正を実行することが好ましい。このために、幾つかの異なる方法で進むことが可能である。
−第１に、焦点スポットに対応してハイドロホンを配置し、２つのトランスデューサにより発せられたパルスの到着時間を測定することが可能である。
−第２に、焦点スポットに対応してリフレクタを配置し、撮像トランスデューサを受信モードで使用して、撮像トランスデューサ自体及び治療トランスデューサにより発せられたパルスのエコーの到着時間を測定することが可能である。
【００２８】
ハイブリッド手法（ハイドロホンを使用して「行き」移動時間を測定し、リフレクタを使用して「戻り」移動時間を測定する）も可能である。
【００２９】
どの方法が使用されても、処置領域全体をサンプリングする複数のポイントに対応して移動時間差を計算又は測定することが必要である。
【００３０】
図３は、キャビテーション気泡ＢＣのクラウドを生成する反復率１００Ｈｚ（周期１０ｍｓ）で発せられた持続時間８μｓを有する治療パルスＩＵＴ、各治療パルスの開始から２５０μｓに及ぶエコー取得の第１の時間窓ＦＴ１、及び第１の窓のそれぞれの終わりから続く治療パルスの放射まで及ぶ能動的撮像の第２の時間窓ＦＴ２を示すタイミング図である。
【００３１】
モータ制御回路ＣＰＭは、治療トランスデューサを旋回させる電気モータを作動させ、電気モータはビーム形成回路ＧＦＦと相互作用して、超音波パルスの焦点スポットＴＦを移動させ、処置領域をスキャンする。モータ制御回路ＣＰＭは、焦点スポットの移動が純粋に電子手段により達成される実施形態ではないことがある（逆に、他の実施形態では、ビーム形成回路は省き得る）。
【００３２】
上述した種々の機能ユニット（ＧＦ、ＣＰＴＩ、ＣＡＴＩ、ＣＰＭ、ＧＦＦ、ＡＰ）は必ずしも物理的に別個の構成要素に対応しない。例えば、１つの集積回路又は回路基板が複数のこれらのユニットの機能の全て又は幾つかを実行してもよい。逆に、１つのブロックの機能は複数の集積回路及び／又は回路基板により実行してもよい。
【００３３】
パルスの特徴は、非限定的な例として与えられた。より一般には、撮像トランスデューサは通常、持続時間０．１μｓ〜５０μｓ（好ましくは０．５μｓ〜２０μｓ）、中心周波数１００ｋＨｚ〜５ＭＨｚ、反復率１Ｈｚ〜１０００Ｈｚを有するパルス列を発し得、これらのパルスは、焦点において、ピーク正圧５０ＭＰａ〜１００ＭＰａ及びピーク負圧−２．５ＭＰａ〜−３０ＭＰａを生成するのに適する。撮像トランスデューサは一般に、治療トランスデューサよりも高い周波数、例えば２〜１５ＭＨｚの周波数で動作し、例えば、２ＭＨｚで超音波を発し得、受信モードでサンプリング周波数８ＭＨｚを有し得る。
【００３４】
先に示したように、本発明による装置は、処置領域の組織を表示する目的の能動モード及びキャビテーション気泡クラウドを表示する目的の受動モードの両方で超音波検査像を取得する。像は次に融合されて、画面Ｅに表示され、これにより、組織に対するキャビテーションクラウドの位置の視認が可能になる。好ましくは、これらの処理動作はリアルタイムで実行される。能動モードで取得された信号により、キャビテーションにより吸収されたエネルギーを推定することもでき、したがって、行っている治療の有効性（及び／又は危険）を評価することができる。
【００３５】
したがって、プロセッサＰＴＤは、治療パルスの強度を自動的に制御し、且つ／又は危険な場合に処理を停止し得る。例えば、プロセッサは、所望の属性（形状、寸法、エネルギー等）を有するキャビテーションクラウドが検出されるまで超音波パルスの強度を徐々に増大させ得、又は逆に、キャビテーションの持続を保証しながら強度を下げ得る。プロセッサは、所望の線量が付与された場合又はキャビテーションクラウドが所望の場所で検出されず、したがって、処置すべきではない組織を破損するリスクがある場合、処理を停止することもできる。
【００３６】
能動的撮像信号は、完全に従来通りに処理し得、したがって、これ以上詳述する必要はない。逆に、キャビテーション気泡の像を効率的に抽出できるように、プロセッサＰＴＤは、受動的に取得された信号に特定の処理を受けさせなければならない。この処理は、
−Ｎ個の一般に連続した治療パルスのエコー信号を取得することであって、通常、Ｎは２〜１０，０００、好ましくは２〜１，０００である、取得すること、
−これらのエコー信号にビーム形成再構築アルゴリズムを適用して、Ｎ個の像を再構築すること、及び
−時空間フィルタリングをそうして再構築された像に適用して、キャビテーション気泡に起因した寄与をそこから抽出すること
を含む。
【００３７】
Ｎの値は、エコー信号の時間的フィルタリングを可能にするために厳密に２以上でなければならない。許容可能な上限は超音波パルスの反復率に依存し、毎秒１０００パルスという高反復率及びＮ＝１０，０００の場合、１０秒毎に１つのフィルタリング像が取得され、これは大半の手術用途で許容可能な最小反復率である。
【００３８】
像再構築は、好ましくは、それ自体既知であるが、行きパス及び帰りパスでの伝播時間差の補償（上述）を含めるように変更された並列−又は超高速−型のビーム形成アルゴリズム（治療パルスに関して実行されるビーム形成と混同すべきではない）を使用して達成される。このアルゴリズムは基本的に、撮像プローブの種々の要素により検出された信号のコヒーレント合算からなり、これらの信号は、行きパス及び帰りパスでの超音波の伝播遅延を補償する（電子集束の原理）ように時間シフトされる。有効な集束のために、超音波源と検出器との間の空間オフセットを考慮に入れる必要があることが理解されよう。スペクトル領域では、オフセットの補償は信号の位相シフトを介して達成される。
【００３９】
時空間フィルタリングは、キャビテーション気泡に起因するエコーを特に組織からのものと区別できるようにする任意のブラインドソース分離アルゴリズムを使用し得る。これは例えば、特異値分解、（粗又は独立した）主成分分析、又は非負行列因数分解等の問題であり得る。以下、特異ベクトル分解を考慮する。
【００４０】
特異値分解（ＳＶＤ）は、行列を代数的に分解する技法である。像の局所統計量への適用に適し、最大エネルギーを少数の固有値に濃縮する。その本質は、Ｕがｍｘｍサイズの直交行列であり、Ｖがｎ×ｎ直交行列であり、Ｓが、対角線要素σｉがＸの特異値であり、その他の要素がゼロであるｍ×ｎ行列である形態に行列Ｘｍｘｎを因数分解することにある。ｎ＜ｍの場合、
【数１】
2021523776.tif 000003
であり、式中、Ｖ＊はＶの共役転置行列であり、Ｓの特異値の数はＸのランクに等しい。
【００４１】
ＳＶＤ法が像処理に適用される場合、最も低い特異値はノイズに関連し、像のエネルギーの大半は高値の特異値にコンパクト化される。本発明の場合、ＳＶＤ時空間フィルタングを使用して、特異ベクトルの賢明な選択を介して、望ましくない低周波成分をキャビテーション気泡クラウドに関連する高周波成分から分離する。これを図４に示す。
【００４２】
各治療パルスに関連し、それぞれが（ｎｘ，ｎｙ）サイズであり、バッファメモリＭＴに記憶された、並列ビーム形成アルゴリズムを使用して再構築されたＮ個１組の「受動」像Ｆ１、Ｆ２、・・・、ＦＮを考慮する。これらのデータは、カソラチ行列と呼ばれる（ｎｘ・ｎｙ，Ｎ）サイズの二次元時空間行列Ｘに再配置し得る。この行列の特異値分解（図４の参照符号ＳＶＤ）の本質は、行列Ｕ及びＶのそれぞれの列を形成する時間及び空間特異値並びに行列Ｓを形成する対応する特異値を見つけることにある。有利なことには、これらのベクトルはエネルギーの低下により並べられる。
【００４３】
実際の時空間フィルタリングの本質は、気泡クラウドを記述する特異ベクトルのみを使用して像を再構築することにある。これらのベクトルにｐ〜ｑ＞ｐのインデックスの連続特異値が関連付けられることが保証される。したがって、フィルタリングされた像ＸＢＣは
【数２】
2021523776.tif 000004
により与えられる（図４における参照符号ＲＩ）。
【００４４】
フィルタリングされたカソラチ行列ＸＢＣ（ｎｘｎＺ，Ｎ）の項は三次元行列ＸＢＣ（ｎｘ，ｎｙ，Ｎ）に再配置され、像の大半のエネルギー領域を位置特定できるようにする積算電力又はキャビテーションマップＣＭを計算することが可能である
【数３】
2021523776.tif 000005
（図４における参照符号ＩＰ）。
【００４５】
Ｎ個の像が１つのキャビテーションマップを与え、換言すれば、取得率がＮで除算されることに留意することが重要である。
【００４６】
キャビテーションマップは、能動的超音波検査により取得された組織の像に重ねられることが意図されるキャビテーション気泡クラウドの最終像として使用される（図４における参照符号ＩＢＣ）。マップＣＭの値を表すために、疑似色を使用し得る。
【００４７】
マップＣＭの空間解像度は、持続時間を伝播速度で乗算した積により与えられる、伝播方向における超音波パルスの長さによって制限される。したがって、可能な限り短く、最短において持続時間が超音波の１サイクルに等しいパルスが好ましくは使用される。パルスの持続時間により意味されるのは、治療トランスデューサの制御に使用される電子パルスの持続時間（すなわち、時間的サポート）であり、その有限帯域幅は必然的に、実際に発せられる超音波パルスの伸長に繋がる。実際には、パルスの持続時間は撮像解像度のみに関して選ばれず、治療プロトコールの要件にも応じて又は何よりも要件に応じて選ばれる。
【００４８】
パラメータＮ、ｐ、及びｑの最適値の選択は、問題となっている特定の用途に依存する。本発明者らは、組織からのエコー信号（なくすことが望ましい）が主に第１の特異ベクトルに集中し、Ｎ−１個の続く特異ベクトルが主にキャビテーション気泡（分離することが望ましい）の寄与を含むことを観測した。したがって、ｐは２（ｐ＝２）に等しく設定し得、ｑはＮ（ｑ＝Ｎ）に等しく設定し得る。パラメータＮに関して、本発明者らは、コントラスト対ノイズ比（以下に定義）が、このパラメータの値が約１０以下に留まる場合、Ｎに伴って増大し、それから横ばいになることを観測した。以下、３つの事例を考慮した：Ｎ＝６、Ｎ＝１０、及びＮ＝１４。
【００４９】
像（より厳密にはキャビテーションマップ）の品質を定義するコントラスト対ノイズ比は、
【数４】
2021523776.tif 000006
により与えられ、式中、＜＞ｉは空間平均を示し、σｉは領域ｉ（ｉ＝１：気泡クラウド、ｉ＝２：背景）におけるこれら２つの領域での標準偏差を示す。領域１及び２は手動で識別される。
【００５０】
変形として、時空間フィルタリングステップは、ビーム形成前に実施し得る。この場合、フィルタリングの空間成分は撮像プローブの要素に関連する。
【００５１】
本発明をテストするために、模型を１．７５Ｌ容量の水及び８％ポリビニルアルコール（ＰＶＡ）と共に準備した。超音波散乱体を１％セルロース（Ｓｉｇｍａｃｅｌｌ、２０μｍ、ＵＳＡ）を使用して添加した。研究所加熱装置を使用して水道水を９０℃まで加熱し、磁性ステアリング旋回ピンを用いて、必要とされる容量をビーカーに注いだ。次に、ＰＶＡを水中に溶解した。混合物を４０℃まで冷まし、セルロースを添加し、溶液をプラスチック製の正方形容器に注ぎ、冷凍庫に８時間配置した。次に、模型を溶かし、再び冷凍庫に８時間配置した。手順の終わりで、模型を離型し、キャビテーション実験のために水槽中に配置した。
【００５２】
本発明による装置のトランスデューサの組立体を２軸ステージ（ＰＩ，Ｍｉｃｏｓ，Ｇｅｒｍａｎｙ）に取り付け、模型を含む水槽内部に配置した。処置する体の領域の生理学的運動をモデリングするために、軸方向運動では最高速度１０ｍｍ／ｓ、運動振幅１０ｍｍ及び側方運動では最高速度１０ｍｍ／ｓ、運動振幅８ｍｍという異なる速度を有するトランスデューサの軸方向及び側方の運動を誘導した。
【００５３】
図５は、このようにして取得されたキャビテーション気泡クラウドの像を示す。クラウドは白色枠により識別される。アーチファクトは像の上下に存在し、水と模型との境界線からの超音波パルスの反射に起因した可能性が最も高い。
【００５４】
Ｎ＝６、１０、及び１４（ここでもｐ＝２及びｑ＝Ｎ）で得られたキャビテーション気泡クラウドの「受動」像のコントラスト対ノイズ比ＣＮＲを、両方とも静的条件下且つトランスデューサが動く状態で測定した。能動的撮像技法を使用しても同様の測定を実行した。本発明者らは、静止条件下でキャビテーション気泡の能動的撮像が受動的撮像よりもわずかによい結果を与えたが、トランスデューサが動く状態では逆であったことを観測した。これは、本発明における受動的方法の使用を確証する。
【００５５】
トランスデューサが動く状態で取得された結果を図６に示し、図６では、参照符号Ｐは受動像のＣＮＲを示し、Ａ１〜Ａ５は能動的方法を使用して取得された像のＣＮＲを示す。
−Ａ１：各治療パルス後の平面超音波の放射を用いて得られた能動的撮像。
−Ａ２：１つの像に結合された、各治療パルス後に有向発散波の１１回の放射のブロックがある超高速シーケンス。
−Ａ３：１０ミリ秒毎に１０の像のブロックを与える、各治療パルス後に発散波の１０回の放射のブロック（結合なし）がある超高速シーケンス。
−Ａ４：１０の像に結合される、３つの治療パルスにわたり分散した有向発散波の１１０回の放射のブロックがある超高速シーケンス。各像は、１１個の有向発散波の結合の結果である。このシーケンスは約３５ミリ秒毎に１０の像のブロックを与える。
【００５６】
実際に、観測されたのは、トランスデューサの動きが、能動的撮像技法の性能の非常に大きな低下（１桁）を誘導し、一方、本発明により得られた像のＣＮＲが約１／２だけしか低下しないことである。これは以下のように説明し得る。処置領域が動く（又は同等に、トランスデューサが動く）場合、組織からの反射は、キャビテーション気泡から反射された有用成分と干渉する高周波非コヒーレント時空間成分を含む。しかしながら、能動的撮像では、観測領域全体における組織は超音波に曝され、したがって、エコーを生成する。これとは対照的に、受動的撮像では、集束超音波が使用され、これは組織に起因した寄生寄与を低減する。
【００５７】
参照文献：
［１］Ｂ．Ａｒｎａｌ，Ｊ．Ｂａｒａｎｇｅｒ，Ｃ．Ｄｅｍｅｎｅ，Ｍ．Ｔａｎｔｅｒ  ａｎｄ  Ｍ．Ｐｅｒｎｏｔ（２０１６）．Ｉｎ  ｖｉｖｏ  ｒｅａｌ−ｔｉｍｅ  ｃａｖｉｔａｔｉｏｎ  ｉｍａｇｉｎｇ  ｉｎ  ｍｏｖｉｎｇ  ｏｒｇａｎｓ．Ｐｈｙｓ  Ｍｅｄ  Ｂｉｏｌ．２０１７  Ｆｅｂ  ７；６２（３）：８４３−８５７．
［２］Ｃｈｒｉｓｔｉａｎ  Ｃｏｖｉｅｌｌｏ，Ｒｉｃｈａｒｄ  Ｋｏｚｉｃｋ，Ｊａｍｅｓ  Ｃｈｏｉ，Ｍｉｋｌｏｓ  Ｇｙｏｅｎｇｙ，Ｃａｒｌ  Ｊｅｎｓｅｎ，Ｐｅｎｎｙ  Ｐｒｏｂｅｒｔ  Ｓｍｉｔｈ，ａｎｄ  Ｃｏｎｓｔａｎｔｉｎ−Ｃ．Ｃｏｕｓｓｉｏｓ（２０１５）．Ｐａｓｓｉｖｅ  ａｃｏｕｓｔｉｃ  ｍａｐｐｉｎｇ  ｕｔｉｌｉｚｉｎｇ  ｏｐｔｉｍａｌ  ｂｅａｍｆｏｒｍｉｎｇ  ｉｎ  ｕｌｔｒａｓｏｕｎｄ  ｔｈｅｒａｐｙ  ｍｏｎｉｔｏｒｉｎｇ．Ｔｈｅ  Ｊｏｕｒｎａｌ  ｏｆ  ｔｈｅ  Ａｃｏｕｓｔｉｃａｌ  Ｓｏｃｉｅｔｙ  ｏｆ  Ａｍｅｒｉｃａ  １３７，２５７３．
［３］Ｋｅｖｉｎ  Ｊ．Ｈａｗｏｒｔｈ，Ｋｅｎｎｅｔｈ  Ｂ．Ｂａｄｅｒ，Ｋｙｌｅ  Ｔ．Ｒｉｃｈ，Ｃｈｒｉｓｔｙ  Ｋ．Ｈｏｌｌａｎｄ，Ｍｅｍｂｅｒ  ａｎｄ  Ｔ．Ｄｏｕｇｌａｓ  Ｍａｓｔ（２０１７）．Ｑｕａｎｔｉｔａｔｉｖｅ  Ｆｒｅｑｕｅｎｃｙ−Ｄｏｍａｉｎ  Ｐａｓｓｉｖｅ  Ｃａｖｉｔａｔｉｏｎ  Ｉｍａｇｉｎｇ．ＩＥＥＥ  ｔｒａｎｓａｃｔｉｏｎｓ  ｏｎ  ｕｌｔｒａｓｏｎｉｃｓ，ｆｅｒｒｏｅｌｅｃｔｒｉｃｓ，ａｎｄ  ｆｒｅｑｕｅｎｃｙ  ｃｏｎｔｒｏｌ，ｖｏｌ．６４，ｎｏ．１．
［４］Ｂａｄｅｒ  Ｋｅｎｎｅｔｈ，Ｈａｗｏｒｔｈ  Ｋｅｖｉｎ，Ｄ．Ｍａｘｗｅｌｌ  Ａｄａｍ，Ｈｏｌｌａｎｄ  Ｃｈｒｉｓｔｙ（２０１７）．Ｐｏｓｔ  Ｈｏｃ  Ａｎａｌｙｓｉｓ  ｏｆ  Ｐａｓｓｉｖｅ  Ｃａｖｉｔａｔｉｏｎ  Ｉｍａｇｉｎｇ  ｆｏｒ  Ｃｌａｓｓｉｆｉｃａｔｉｏｎ  ｏｆ  Ｈｉｓｔｏｔｒｉｐｓｙ−Ｉｎｄｕｃｅｄ  Ｌｉｑｕｅｆａｃｔｉｏｎ  ｉｎ  Ｖｉｔｒｏ．ＩＥＥＥ  Ｔｒａｎｓａｃｔｉｏｎｓ  ｏｎ  Ｍｅｄｉｃａｌ  Ｉｍａｇｉｎｇ．
［５］Ｐ．Ｂｏｕｌｏｓ，Ｆ．Ｖａｒｒａｙ，Ａ．Ｐｏｉｚａｔ，Ｊ．Ｃ．Ｂｅｒａ，Ｃ．Ｃａｃｈａｒｄ（２０１５）．Ｐａｓｓｉｖｅ  ｃａｖｉｔａｔｉｏｎ  ｉｍａｇｉｎｇ  ｕｓｉｎｇ  ａｎ  ｏｐｅｎ  ｕｌｔｒａｓｏｎｉｃ  ｓｙｓｔｅｍ  ａｎｄ  ｔｉｍｅ  ｒｅｖｅｒｓａｌ  ｒｅｃｏｎｓｔｒｕｃｔｉｏｎ．２２ｅｍｅ  Ｃｏｎｇｒｅｓ  Ｆｒａｎｃａｉｓ  ｄｅ  Ｍｅｃａｎｉｑｕｅ  Ｌｙｏｎ．
【図１】
2021523776.tif 000007
【図２】
2021523776.tif 000008
【図３】
2021523776.tif 000009
【図４】
2021523776.tif 000010
【図５】
2021523776.tif 000011
【図６】
2021523776.tif 000012
【国際調査報告】
2021523776.tif 000013
2021523776.tif 000014
2021523776.tif 000015
2021523776.tif 000016
2021523776.tif 000017