JP5661329B2 - 医用超音波イメージングにおける高密度焦点式超音波に対するフィードバックを行うシステムおよびコンピュータプログラム - Google Patents

Description

本発明は、超音波イメージングに関する。殊に超音波イメージングは、高密度焦点式超音波（ＨＩＦＵ high intensity focused ultrasound）におけるフィードバックに使用される。

連邦政府の後援による研究または開発
米国政府は、本願発明において使用料の支払いを要しない実施権を有しており、また特殊な場合にＤＡＲＰＡによる契約条件の第HR 0011-08-3-0004号の条項によって規定された適当な条件に基づいて第３者に特許の使用許可を与えることを特許権者に要求する権利を有するものである。

ＨＩＦＵは、癌、腫瘍、病変または他の好ましくない組織構造を治療するのに使用される。超音波エネルギは、このような組織を十分に熱して、好ましくない組織を壊死させる。この超音波エネルギは、健康な組織の損傷を回避するために集束される。超音波による治療により、侵襲性の処置、例えば、手術または高周波切除処置などを回避することができる。

超音波イメージングは、ＨＩＦＵ療法をガイドするのに使用される。このイメージングにより、上記のような所望の組織への治療パルスの集束が支援される。例えば、同一のアレイがＨＩＦＵのイメージングおよび送信に使用されるため、このＨＩＦＵは、所望の組織に集束されるのである。しかしながらＨＩＦＵでは、イメージングに使用されるのとは異なるアレイが使用されることがある。

また上記のような治療の最中に組織の温度的および生物学的な変化を監視するためにさまざまな試みも行われてきている。例えば、超音波エネルギは、熱膨張係数（例えば、スペックルトラッキングによる組織膨張の測定）、組織における音の速度、またはスチフネスの変化（例えば、ストレインイメージング）を測定するのに使用される。しかしながらこれらの診断ベースの超音波組織の特徴付けは、十分な信号対雑音分解能を有しないか、または臨床場面において活用可能ではない。

本発明の課題は、従来技術の欠点を解消する、高密度焦点式超音波に対してフィードバックを提供するための方法、コンピュータ読み出し可能媒体および高密度焦点式超音波に対するフィードバックを行うシステムを提供することである。

上記の方法についての課題は、本願発明の請求項１により、高密度焦点式超音波に対するフィードバックを行う方法において、高密度焦点式超音波トランスデューサからの励起を患者の組織に送信するステップと、この励起によって生じた上記の変位を検出するステップと、この変位の関数として、上記の励起のビームプロフィールの画像を生成するステップとを有することを特徴とする、高密度焦点式超音波に対するフィードバックを行う方法によって解決される。

またコンピュータ読み出し可能媒体についての課題は、本願発明の請求項１２により、コンピュータ読み出し可能記憶媒体において、このコンピュータ読み出し可能記憶媒体には、高密度焦点式超音波に対するフィードバックを行うためにプログラムされたプロセッサによって実行される命令を表すデータが記憶されており、この記憶媒体には、組織の変位の関数として、高密度焦点式超音波治療に対するビーム位置を求めるための命令およびこのビーム位置の画像を生成するための命令が含まれていることを特徴とする、コンピュータ読み出し可能記憶媒体によって解決される。

さらに上記のシステムについての課題は、本願発明の請求項１９により、高密度焦点式超音波に対するフィードバックを行うシステムにおいて、このシステムは、高密度焦点式超音波治療波形を生成する高密度焦点式超音波トランスデューサと、上記の高密度焦点式超音波トランスデューサとは別個に移動可能でありかつ上記の高密度焦点式超音波トランスデューサからの送信に応答して信号を受信するイメージングトランスデューサと、この受信信号の関数として、空間位置を表すデータを出力する受信ビームフォーマと、この出力データの関数として組織の変位を推定しかつこの組織の変位の関数としてビームプロフィールの画像を生成するプロセッサと、当該のビームプロフィールの画像を表示するためのディスプレイとを有することを特徴とするシステムによって解決される。

高密度焦点式超音波に対するフィードバックを行う方法の１実施形態の流れ図である。 高密度焦点式超音波に対するフィードバックを行うシステムの１実施形態のブロック図である。 組織領域および検出したビームプロフィールを示す医用画像である。

最初に述べておきたいのは、以下に説明する好適な実施形態には、高密度焦点式超音波用にフィードバックを提供するための方法、命令およびシステムが含まれていることである。ＨＩＦＵトランスデューサからのビームの位置は、超音波イメージングを使用して決定される。この超音波イメージングにより、ＨＩＦＵトランスデューサから送信されたビームによって生じた組織の変位が検出される。この変位またはこの変位から導出された情報は、上記のビームの位置の検出およびイメージングに使用可能である。ＨＩＦＵおよびイメージングに対して別個のトランスデューサを使用可能である。ユーザは、超音波イメージングからの変位フィードバックによってＨＩＦＵ送信の狙いを定める。ここでは接続されたトランスデューサまたは共通のトランスデューサを使用することができる。

第１の態様によれば、高密度焦点式超音波に対してフィールドバックを行う方法が提供される。励起が高密度焦点式超音波トランスデューサから患者の組織に送信される。この励起によって生じた組織の変位が検出される。上記の励起のビームプロフィールの画像が、上記の変位の関数として生成される。

第２の態様では、コンピュータ読み出し可能記憶媒体が、高密度焦点式超音波に対してフィードバックを行うためにプログラムされたプロセッサによって実行される命令を表すデータを記憶している。この記憶媒体には、組織の変位の関数として、高密度焦点式超音波に対するビーム位置を決定するための命令と、このビーム位置の画像を生成するための命令とが含まれている。

第３の態様では高密度焦点式超音波に対するフィールドバックを行うシステムが提供される。高密度焦点式超音波トランスデューサは、高密度焦点式治療波形を生成するために使用可能である。イメージングトランスデューサは、上記の高密度式超音波トランスデューサとは別個に移動可能である。このイメージングトランスデューサは、上記の高密度式超音波トランスデューサからの送信に応答した信号を受信するために使用可能である。この受信信号の関数として空間位置を表すデータを出力するため、受信ビームフォーマを使用可能である。この出力データの関数として、組織の変位を推定するため、またこの組織の変位の関数としてビームプロフィールの画像を生成するためにプロセッサを使用可能である。このビームプロフィールの画像を表示するためにディスプレイを使用可能である。

本発明は、特許請求の範囲によって規定されるものであり、この明細書に記載されているものは、特許請求の範囲を限定するものとみなすべきではない。以下、好適な実施形態に関連して本発明の別の態様および利点を説明する。これらは後に独立してまたは組み合わせで権利主張され得る。

コンポーネントおよび図は、必ずしも縮尺通りではなくむしろ本発明の原理を説明するために強調されている。さらに複数の図において、類似の参照番号は、見方が変わっても全体を通じて相応する部分を示すものとする。

変位からの上記のビームのエミュレーションまたは治療ビームの検出は、位相収差および／または減衰を説明することができる。適用量分布は、弾性、ひずみ、剪断などの変位イメージングを使用することによってより良好に推定することができる。変位のイメージングの結果は、上記の治療ビームの分布を補正するのに使用可能である。上記の治療トランスデューサの出力および／または焦点位置は、上記のフィードバックに基づいて調整可能である。ＨＩＦＵトランスデューサは、上記のフィードバックに基づいて照準などの調整を行うことができる。

図１には、高密度焦点式超音波に対してフィールドバックを行う方法が示されている。この方法は、図２のシステムまたは別のシステムによって実現される。付加的なステップ、別のステップを設けるか、またはステップを少なくすることが可能である。例えば、ステップ３８〜４６が行われない実施形態もある。またこれらのステップは、説明する順序または図示の順序で実行されるが、別の順序で実行することも可能である。例えば、ビーム位置は、ステップ３６で画像を生成する前にステップ３８で補正される。別の例として、ステップ３２〜２６は、ステップ４２のインタリーブの一部としてステップ４０と共に繰り返される。

ステップ３０では、励起が高密度焦点式超音波トランスデューサから送信される。この励起は、音響エネルギである。この音響エネルギは集束されて、３次元のビームプロフィールが得られる。上記の励起は、フェーズドアレイおよび／またはメカニカルフォーカスを使用して集束される。この励起は、１つの次元において、例えば高さの次元において集束しないことも可能である。励起は、患者の組織に送信される。励起は、治療に対する位置、例えば腫瘍に集束される。しかしながら音響的な作用、位置合わせの不良または別の要因に起因して、上記のビームは、正しくオーバラップしないかまたカバーしないことがあり、または治療に対する組織の位置で交わらないことさえもある。

上記の励起は、治療励起として生成される。択一的には励起より、この治療励起をエミュレートする。上記の高密度焦点式超音波治療波形がエミュレートされる。一般的には上記の治療励起と同じ焦点、振幅、周波数または別の特性が、上記のエミュレーションに使用される。このエミュレーションは、実質的に治療効果を生じさせないようにするために使用される。例えば、振幅、持続時間またはその双方が、治療波形と比べて低減される。「実質的に」治療効果を生じさせなくすることによって、領域、例えば治療の領域を生成することができる。収差または焦点歪みに起因して閾値レベル以上に１つの点が加熱されることはあるが、上記の治療の領域全体では、上記のエミュレーションによる治療効果が生じることはない。治療効果が生じないようにすることにより、組織が変質されるかまたはキャビティが形成される点まで加熱されることを回避することができる。例えば、生物学的な作用には、約４１〜４５℃の組織温度における温熱療法、４５℃以上の温度におけるプロテイン変質、および５０℃以上の温度における組織壊死が含まれ得る。組織スチフネスは、４５℃以下の温度であっても生じることがある。４５℃以上の温度では、粘度が増大しおよび／またはスチフネスが生じることがある。５０℃以上の温度では、組織はスチフネスが高くなり、および／または減衰が大きくなる。

ステップ３２では上記の組織の変位が検出される。上記の励起によってこの組織の変位が発生する。この変位は、縦波によって生じ得る。変位は、択一的または付加的に剪断波によって生じ得る。

上記の変位は、超音波スキャニングによって検出される。この変位の位置、大きさ、タイミングおよび／または別の特性が検出される。例えば、上記のエミュレーションによって生じた変位の閾値量に関連した位置が検出される。

上記の変位を検出するため、変位中の組織に超音波エネルギが送信され、このエネルギの反射が受信される。この送信および受信は複数回行われ、変位による変化が求められる。任意の送信および受信シーケンスを使用することができる。上記の変位は、空間位置毎の違いから検出される。例えば、速度、分散、強度パターンのシフト（例えば、スペックルトラッキング）または別の情報が、上記の受信したデータから、変位として検出されるのである。

任意の公知の変位イメージングまたは将来開発される変位イメージングを使用することができる。例えば、診断用超音波に対して調整されるレベル以下の強度および持続時間を有する診断パルスが送信される。例えば、１〜５サイクル持続時間のパルスが、７２０ｍＷ／ｃｍ²以下の強度で使用される。他の強度のパルス、例えば１０００ｍＷ／ｃｍ²以下のパルスを使用することも可能である。この超音波送信は、治療する組織を含む領域に集束される。この送信は１つまたは複数のラインをカバーすることができる。例えば、ビーム幅の広い送信パルスを使用して、２次元または３次元的な分布を有する２つまたはそれ以上の受信スキャンラインに沿って受信する。択一的には、単一の受信ビームが、送信に応答して形成される。領域は順次にスキャンすることができる。１つまたは複数の測定が、受信スキャンライン毎に行われる。

２つまたはそれ以上、例えば２〜１０個のパルスが、測定毎にまたは組み合わせ式の測定に対して同じ位置に送信される。択一的には測定毎に単一のパルスを送信することができる。治療強度および休止後時間が既知である場合、単一のパルスを使用し、またプレエミュレーションまたは励起測定と比較して位置の変化を求めることができる。

上記の治療用超音波をエミュレートする励起が休止した後、組織は、弛緩位置に移動する。比較的小さい複数の診断イメージングパルスからのエコーが受信される。これらのエコーは１つまたは複数の画像を生成するために使用されて、上記の治療に関連する励起によって生じた変位を有する位置が識別される。

これらのエコーは、Ｂモードまたはドップラ検出を使用して検出される。Ｂモードデータを使用して、複数のパルスから得られた上記のデータが相関付けされる。この相関は、１次元、２次元、または３次元である。例えば、１つのスキャンラインに沿いかつ上記のトランスデューサに向かう方向またはこれから離れる方向の相関が使用される。任意の公知の相関または将来開発される相関を使用することができ、例えば、相互相関、パターンマッチング、または絶対差の最小和を使用することができる。組織構造および／またはスペックルが相関される。ドップラ検出を使用する場合、クラッタフィルタにより、移動する組織に関連した情報を通過させる。上記の組織の速度は、複数のエコーから導き出される。この速度は、上記のトランスデューサに向かう方向またはこれから離れる方向の変位を求めるのに使用される。択一的には、異なる位置における速度間の相対関係または差分によって、ひずみまたは変位を示すことも可能である。

上記の変位の量は、ＨＩＦＵに関連した上記の励起によって影響を受ける領域を表す。特定の変位に関連付けられる時間により、減衰曲線を推定することができる。時間の関数として上記の変位を測定することにより、励起の休止以来のひずみの減衰を測定することができる。変位の単独で測定するかまたは上記の減衰の任意の変位特性を測定することができる。

別の実施形態では、ストレインイメージングまたは弾性イメージングを使用する。上記の組織の変位は、時間の関数として求められる。この変位は、Ｂモード超音波データのような組織データから測定することができる。相関、相互相関、絶対差の最小和または別の類似の量を使用して、複数のスキャンの間における上記の変位が求められる。これらの変位は、１次元、２次元または３次元に沿って求められる。１実施形態では、米国特許第５１０７８３７号、第５２９３８７０号、第５１７８１４７号、第６５０８７６８号または第６５５８３２４号明細書に記載されている任意の１つまたは複数の方法またはシステムを使用して、ストレイン情報としての画像またはデータの弾力フレームを生成する。ここで上記の明細書の開示内容は、参照によって本明細書に組み込まれるものとする。ひずみを測定する別の方法も使用可能である。上記の変位は、組織速度および／または加速度を求めることによって測定することができる。

上記の変位は、任意サイズの領域にわたって検出される。１実施形態では、この変位は、治療しようとする組織が含まれていそうな関心対象領域において検出され、これは、例えばＢモードイメージングに対する完全なスキャン領域の約１／３〜１／２である。図３には１実施形態が示されており、ここでは左側の画像により、関心対象領域ボックスが変位イメージングに対して示されている。イメージング領域全体にわたって変位を検出するなど、より大きな関心対象領域、より小さな関心対象領域を使用することもでき、または関心対象領域を使用しないことも可能である。関心対象領域を狭くすることにより、ＨＩＦＵに関連する励起波形を送信する繰り返し回数を少なくして変位を検出することができる。形成される受信ビームの数およびサンプル密度に依存して、繰り返しなし、１回の繰り返し、またはそれ以上の繰り返しを使用可能である。Ｂモードサンプル位置毎に変位をサンプリングするなどの完全サンプリングを使用することができる。Ｂモードスキャングリッドよりも多くの変位をサンプリングすること、または少ない変位をサンプリングすること（例えば、疎な）が可能である。

１実施形態では、縦波に加えてまたはこれとは択一的に剪断波が検出される。上記の励起によって１つのビームが形成され、このビームによって複数の空間位置に１つの剪断波が生成される。上記のビームが十分に強い個所で、１つの剪断波が生成される。この剪断波は、音響波放射方向に沿って縦波よりもゆっくりと組織を通って伝搬する。この剪断波は、加わる応力の方向に対して垂直な方向を含めてさまざまな方向に伝搬する。この剪断波の変位は、上記の励起ビームに相応する位置において比較的大きい。

超音波データを取得する。この超音波データのうちの少なくともいくつかは剪断波に応答するデータである。関心対象領域を監視してこの剪断波を検出する。この関心対象領域のサイズは、例えば横方向に６ｍｍ，軸方向に１０ｍｍなど任意である。この検出領域は、超音波によって監視される。例えば、Ｂモードスキャンを実行して、上記の剪断波によって生じる組織の変位が検出される。ドップラモード、カラーフローモード、または別の超音波モードを使用して上記の剪断波を監視することができる。

この監視は、任意個のスキャンラインに対して実行される。例えば、各送信に応答して４つの受信ビームが形成される。上記の励起を送信して剪断波を生成した後、単一のスキャンラインに沿ってＢモード送信が繰り返し行われ、また隣接する４つのスキャンラインに沿った受信が行われる。任意の回数の繰り返し、例えば、約１２０回の繰り返しを使用することができる。上記の超音波データのうちのいくつか、例えば、上記の繰り返しの開始時または終了時における超音波データは上記の剪断波に応答するものではないことがある。別の実施形態では、単一の受信ビームだけまたは別の個数の受信ビームが各送信に応答して形成される。

上記の剪断波は、スキャンラインを通って伝搬するため、上記のＢモード強度は変化することがある。上記の監視されるスキャンラインに対し、上記の剪断波によって発生する組織運動の時間プロフィールを表す一連のデータが得られる。例えば、（例えば、複数のスキャンラインに沿った）複数の空間位置から得られるデータは、時間の関数として相関付けされる。任意の弾性検出を使用することができる。深さまたは空間位置毎に複数の深さまたは空間位置にわたって相関付けが行われる（例えば、プロフィールの計算される点である中心深さを有する６４個の深さの核）。空間における２次元または３次元の変位を使用可能である。上記のスキャンラインまたはビームとは異なる方向に沿った１次元変位を使用可能である。

所定の時点における最大の相関または十分な相関を有する空間オフセットにより、変位の量が示される。所定の位置に対して異なる時点に複数の変位が求められる。所定の位置に対する時間についてのこのプロフィールにより、上記の剪断波の検出が示される。このプロフィールは、ノイズのないまたは単一の変化の事例について調べられる。時間ローパスフィリングを伴うまたはこれを伴わない上記のプロフィールにおけるピークは、上記の剪断波面の通過を表す。最大変位が選択されるが、平均または別の変位統計を使用可能である。所定の位置における最大剪断が検出される。択一的には平均または別の剪断が検出される。

より大きな領域を監視するためには、監視送信ビームに応答して付加的な複数の受信ビームを形成する。択一的には別の剪断波を生成し、この剪断波生成点から異なる距離で送信ビームおよび受信ビームを準備する。上記の６ｍｍ×１０ｍｍの例では、３６個の受信スキャンラインを準備することができる。上記の処理を、１送信ビーム当たりに４つの受信ビームで、異なる横方向間隔に対して９回繰り返す。受信ビーム位置毎に、超音波データによって表される動き情報の時間プロフィールが得られる。同一の剪断波を監視するために異なるスキャンラインに沿って送信を行うことは、高い時間方向解像度を得るため、上記の時間プロフィールの形成中には回避されるが、インタリーブされたまたはシフトするスキャニング位置を設けることが可能である。

上で述べたことは１つの深さに対するものである。上記の関心対象領域の軸方向の拡がり全体をカバーする１つのゲートを得るため、上記のサンプリングを準備することができる。別の１実施形態では、受信ビーム毎に複数の深さにおいてサンプルが取得される。軸方向の深さおよび横方向位置毎に別の時間プロフィールを得る。例えば５ｍｍに対して約２００個または１０ｍｍに対して４００個などの任意個の深さを使用することができる。

上記の関心対象領域における異なる位置を表す超音波データが得られる。この超音波データは、スキャニングによってリアルタイムに得られるか、またはメモリから得られる。上記の動き情報は、位置毎に異なる時点における上記の応答を表す。別のスキャニング、別の監視または別の技術を使用し、超音波データを得て、剪断の大きさを推定することができる。

上記の剪断波の検出と、治療ＨＩＦＵの実際の適用と交互に行われる個所において、上記の組織は時間と共に変化することがある。上記の組織のスチフネスは、上記の治療個所において増大することがある。スチフネスのこの変化により、上記の剪断の検出した大きさが変わることがある。剪断速度および／または弾性率または剪断の別の複合的な表現を使用して、上記の組織における変化の影響を最小化または回避することができる。別の方法で弾性率が利用できないかまたは求めるのが難しいところでは、剪断速度が好適である。吸収係数は、治療する組織に依存して推定することができる（例えば、０．５，０．６，０．８または別の値）。上記の剪断速度および／または弾性率は、圧力および吸収係数に基づき、少なくとも部分的に求めることができる。この加えられる圧力は、上記の送信される励起から、また減衰を考慮することによって既知になる。

上記の組織の異なる空間位置に対し、剪断速度が検出される。位置毎に時間の関数として上記の変位が求められる。上記の剪断速度は、剪断波の生成から、異なる位置においてこの剪断波が検出されるまでの時間を求めることによって得られる。この時間および上記の位置までの距離によって上記の速度が求められる。この距離は、スキャンライン間隔（すなわち、上記の剪断波を生成するための送信ビーム位置およびこの剪断波を検出するための受信ビーム位置）からわかる。上記の時間は、剪断波の生成と、検出との間の相対的な時間からわかる。

上記のプロフィールにおけるピークを検出するために別の手法を使用することできる。例えば、回帰が適用される。上記の剪断波は線形であるため、自動化されたアウトライヤ検出を伴うロバスト線形回帰により、剪断波速度を示すことができる。上記の関心対象領域におけるすべてのサンプル点に対する超音波データは、距離に対し、時間の関数としてプロットされるか、または時間および距離によってプロットされる。上記のプロットまたはデータに対して線形回帰が適用され、これらのデータに対する直線当てはめが得られる。この線に対する傾斜が、剪断波速度を表す。

別のアプローチを使用することも可能である。例えば、異なる時点から得られるデータを相関付けて、上記の剪断波によって生じる組織のシフトを検出する。別の例では、上記の時間プロフィールから特徴が抽出される。主成分分解（Principle component decomposition）を使用することができる。異なる時間プロフィール間の相関付けが実行される。異なる時間プロフィールに対し、異なる距離に関連した遅延により、上記の速度が得られる。択一的にはウェーブレット解析を行うことができる。ウェーブレット変換を上記の時間プロフィールに適応して上記の剪断波に相応するピークを識別する。各空間位置へのこのピークの移動時間から速度値が識別される。

全関心対象領域から得られるすべてのピーク移動時間を一緒に、例えば線形回帰において使用することができる。特徴抽出または回帰に対して１つの深さに対するデータを使用するなど上記のデータのサブセットだけを使用することが可能である。剪断速度は、上記の関心対象領域における位置毎に計算することができる。択一的には関心対象領域内の剪断波速度分散の空間表現を得ることができる。

他のアプローチにおいて、任意の弾性率または剪断値を推定することができる。組織弾性率値は、その位置における硬さまたはスチフネスを表す。例えば、組織の剪断弾性率が推定される。択一的な実施形態では、ヤング弾性率が推定される。別の実施形態では、質的または量的な別の剪断値が推定される。

上記の剪断弾性率は、ｇ＝ρｖ_s ²によって得られ、ここでρは密度、ｖは推定した剪断速度である。１実施形態では、上記の組織弾性率または剪断情報、例えば剪断弾性率は、ひずみまたは変位および上記の弾性率または剪断情報の関数として求められる。例えば、サンプル位置毎の剪断弾性率が、拡散方程式を繰り返して解くことによって求められる。０．５のポワソン比を仮定するかまたは既知のポワソン比を使用することにより、異なる時間におけるストレインフィールドの関数として、または異なる位置および剪断弾性率に対する異なる応力の下で、異なる位置における剪断弾性率が繰り返して計算される。

ステップ３４では、ＨＩＦＵに対するビーム位置が、上記の組織情報の変位を使用して求められる。剪断弾性率、剪断、剪断速度または変位が十分に大きい位置が識別される。上記の変位が比較的大きい位置が、閾値を適用することによって識別される。閾値に対し、またはあらかじめプログラミングされるかまたは所定のデータセットに適合された別の値に対して相対的に行われる。上記のビームの予想される位置から離れた空間位置におけるデータに基づく閾値など、上記の閾値は正規化することができる。別の例では、関心対象領域にわたる平均変位が求められる。この平均よりも大きな極大変位に関連した位置がビーム位置を示す。択一的には、あらかじめプログラミングされた閾値または別の閾値が適用される。別の実施形態または付加的な実施形態では、閾値は適用されないか、またはノイズ閾値が使用される。変位の範囲は、ディスプレイ値にマッピングされて、変位の少ないまたは変位のない領域が、ダイナミックレンジの一方の端になり、また最大の変位がこのダイナミックレンジの他方の端部になるようにされる。線形のマッピングまたは非線形のマッピングを使用可能である。

上記の変位データは、上記の閾値を適用する前に空間的にフィルタリングすることも、またはフィルタリングしないことも可能である。これらの変位は、上記の閾値を適用した後、ローパスフィルタリングすることができる。

ステップ３６では、上記の変位の関数として画像データが生成される。例えば、この画像は、上記の励起のビームプロフィールを表す。このビームプロフィールは、十分な変位を有する位置に相応する。この画像は、上記のビームの空間分布を含めて、このビームの位置を示す。上記の画像は、１つまたはそれ以上の次元においてビームプロフィールの空間的な拡がりを表す。図３には上記のビームの右側における画像が示されている。比較的明るい位置または高輝度の位置は、大きい変位に相応する。この画像には、ローパスフィルタリングなどのフィルタリングを行うことができる。

上記の変位データは、ディスプレイフォーマットであるか、またはディスプレイフォーマットにスキャンコンバートすることができる。変位データはカラーまたはグレイスケールデータとすることが可能であるが、グレイスケールまたはカラースケールとマッピングする前のデータとすることも可能である。この情報は、ディスプレイ値に線形または非線形にマッピングすることができる。

この画像は、異なる位置に対する剪断または弾性率（例えば剪断弾性率）などの変位情報を表す。関心対象領域または視野においてすべてのグリッドポイントに対して値を決めた場合、上記のディスプレイのピクセルは、その領域に対するビームを表す。このディスプレイグリッドは、スキャングリッドおよび／または変位が計算されるグリッドとは異なるものとすることも可能である。カラー、明るさ、輝度、色相またはべつの特性が、変位の関数として変調される。

上記の画像は別のデータを含むことも可能である。例えば、Ｂモードまたは同じ領域における組織、液体または造影剤を表す別のデータが含まれる。上記の変位データは、上記の別のデータのオーバレイまたは別のデータとの組み合わせに使用される。上記の別のデータによって支援されるのは、治療する組織を基準にしてビームの位置をユーザが決定することである。

択一的な１実施形態において上記の画像は、ビームプロフィールモデルおよび変位の関数として生成される。上記の変位は、ノイズを含んでいることがある。上記の治療ビームのモデルが準備される。このモデルは、経験的な情報または理論的な計算を使用して生成することができる。上記のビームプロフィールは、このビームプロフィールの予想分布を求める問題の逆問題として作成される。上記の空間および時間的な変位は、このビームプロフィールモデルに基づいて定式化することができる。上記の測定した変位およびモデルが得られれば、上記のモデルに対し、上記の測定の最小二乗当てはめまたは最善の当てはめの表現が求められる。この当てはめ表現は、上記の当てはめに基づいてサイズおよび形が決められるグラフィックオーバレイなどの画像を生成するのに使用される。この画像によって表されるビームプロフィールは、上記の測定した変位にモデルを当てはめることによって再構成される。

ステップ３８では、上記のビームプロフィールが音響伝搬の関数として補正される。上記の測定した変位は、加わる応力の、深さに依存する減衰について補正されることも補正されないこともある。圧力は組織を通して伝搬するため、この圧力は減衰する。圧力源から離れた位置（上記の発生源を基準にした深さ）では、上記の減衰により、生じる動きまたは変位は小さい。上記の変位は、減衰を表すように調整され、これにより、異なる深さにおいてさらに正規化された変位または応力が得られる。

上記の補正は、応力の発生源（例えば、トランスデューサ）の点または領域から離れた距離の関数として、線形である。組織モデルまたは異なるタイプの組織に基づいて、非線形の補正を使用することも可能である。経験的なデータに基づいてまたは伝搬モデルに基づいて、上記の線形の関数または非線形の関数が想定される。音響の及ぼす力に対し、組織における音の減衰が、距離および周波数の関数として補正される。別の実施形態では、減衰および／または周波数に対して補正は行われない。

ステップ４０では、高密度焦点式超音波治療波形が送信される。高電圧の波形が、上記の高密度焦点式超音波トランスデューサに加えられ、このトランスデューサにより、音響ドメインにおいてＨＩＦＵ治療波形が生成される。このＨＩＦＵパルスは、フェーズドアレイまたは機械式焦点調節を使用して集束され、焦点位置またはビーム位置の組織に高密度音響エネルギが供給される。上記の治療超音波パルスは、任意の所望の周波数において複数のサイクルを有する。１実施形態において上記の治療パルスは、１つの超音波周波数で、例えば５００ｋＨｚ〜２０ＭＨｚの周波数で１秒の数分の１または数秒にわたって持続する。例えば、１平方センチメールあたり１００ワット以上、１平方センチメールあたり５００ワット以上、１平方センチメールあたり１０００〜２０００ワット、または１平方センチメールあたり約１０００ワットの任意のピーク強度を供給することができる。任意の強度、周波数および／またはサイクル数を有する任意の公知の治療波形または将来開発される治療波形を使用可能である。上記の波形は、連続的かまたは間欠的である。

上記の治療超音波パルスは、所望の組織位置において熱を生成することによって組織を治療する。上記の強度はまた組織に応力を発生させる。上記のパルスは、負および正の音響圧力により、トランスデューサに向かってまたはこれから離れるように上記の組織を押す。十分に長い治療パルスに対して、上記の組織に実質的に一定のひずみが形成される。ひずみεは、組織のスチフネスＥと、粘度ηと、ＨＩＦＵ放射圧からの応力との関数である。上記の治療パルス中の安定状態応力は、組織における音速ｃに対する平均ＨＩＦＵ強度Ｉの比に正比例する。

上記のＨＩＦＵ波形はまた、検出可能なバイオメカニカルな変化を発生させる。上記の治療音響エネルギの熱効果は、熱膨張に起因する体積の変化、音速（ｃ）の変化、組織のスチフネス（Ｅ）の変化、および／または組織における液体の粘度（η）の変化を発生させることがある。上記の治療音響エネルギはまたは、放射圧、流動および／またはキャビテーションのような機械的な作用を引き起こすことがある。生物学的な作用には、約４１〜４５℃の組織温度における温熱療法、４５℃以上の温度におけるプロテイン変質、および５０℃以上の温度における組織壊死が含まれ得る。組織スチフネスは、４５℃以下の温度であっても影響を受けることがある。４５℃以上の温度では、粘度および／またはスチフネスが増大することがある。５０℃以上の温度では、組織はスチフネスが高くなり、および／または減衰が大きくなる。

上記の組織に治療を施す前、ステップ３６で生成した画像からＨＩＦＵビームの位置を求める。ユーザは、トランスデューサ、焦点位置を再配置するかまたは別の設定を変更して、治療する組織にわたってビームを位置決めすると同時にＨＩＦＵの影響を受ける健康な組織を最小化する。別の実施形態では、自動化された位置決めを使用することができ、ここでは、治療する組織がコンピュータ支援診断を使用して求められる。上記のビームを所望のように位置決めした後、ステップ４０においてＨＩＦＵを送信する。

このＨＩＦＵは、連続的または散発的なものとすることができる。任意の治療レジメンを使用可能である。進行中の治療の間に、または治療の異なる部分の間に、ステップ３６のイメージングを行うことができる。ステップ４０の治療波形には、ステップ３６のイメージングが挟み込まれる。ステップ３６のイメージングは、ステップ３０の送信と、ステップ３２の検出と、ステップ３４の決定とを使用して行われる。上記のＨＩＦＵ治療は、上記のビーム位置が求められる間、中止される。択一的にはステップ４０のＨＩＦＵ送信は、ステップ３２において変位を検出するための励起として使用される。上記のＨＩＦＵは変位を検出するために中止することができ、または上記のＨＩＦＵ波形に起因する変位が使用される。別の択一的な実施形態では、上記のＨＩＦＵは１つの周波数またはコーディングで実行され、複数の励起の送信および相応の受信が別の１つの周波数またはコーディングで行われ、これにより、前処理組織位置情報が与えられれば、同時に処理を行うことができる。このような交互実行方式により、ユーザまたはシステムは、進行中の状態をベースにしてＨＩＦＵビームを位置決めすることができる。患者またはトランスデューサの位置がずれた場合、ビームを変化させて適切な組織が治療されるようにすることが可能である。組織における音速が治療に起因して変化した場合、上記のビームを変化させて、適切な組織が治療されるようにすることができる。

ステップ４４では、高密度焦点式超音波治療による組織の変化が求められる。例えば、変位の変化が求められる。ＨＩＦＵは、組織な弾性を低下させるないしは硬くすることがある。アブレーション、コラーゲン変性、凝固または別の作用により、上記の剪断速度または別の特性が変化し得る。応力が同じであるか、または既知であるが応力が異なるのに応じて、変位の量を求めることができる。この変位の量または大きさは、測定可能である。１つの領域に対する変位における変化の中央値または平均などの任意の測定値を使用することができる。剪断、ひずみ、弾性、弾性率、速度または別の組織特性における変化を測定することができる。

上記の変化は、ユーザまたはシステムによって求められる。例えば、数量が決定される。別の例では、１つまたは複数の画像に基づいてユーザが上記の変化を検出する。

この変化は、適用量のフィードバック制御に使用することができる。ステップ４２では、上記の変化に基づき、ステップ４０におけるＨＩＦＵの適用を変更または中止する。高密度焦点式超音波治療の適用量は、上記の変化に依存して変更される。健康な組織へのダメージを最小化するため、十分に治療が行われた個所ではＨＩＦＵ強度および持続時間を低減することができる。組織における上記の変化により、治療が十分であることが示される。効果のない治療を回避するため、予想される変化よりも小さいことによってわかる不十分な治療が行われた個所ではＨＩＦＵ強度または持続時間を増大させることができる。

図２には、高密度焦点式超音波に対するフィードバックを行うシステム１０の１実施形態が示されている。システム１０により、図１の方法または別の方法が実現される。システム１０には、送信ビームフォーマ１２と、トランスデューサ１４と、受信ビームフォーマ１６と、画像プロセッサ１８と、ディスプレイ２０と、記憶装置２２と、ＨＩＦＵトランスデューサ２４とが含まれている。付加的なコンポーネント、別のコンポーネントを設けるかまたはコンポーネントの数を減らすことが可能である。例えば、ユーザ入力部を設けて、ビームプロフィールが手動で変更されるかまたは支援によって変更されるようにする（例えば、周波数、焦点深度、スキャンライン角度、アパーチャ、焦点位置および／またはアポディゼーションの選択）。システム１０は、医用診断超音波イメージングシステムである。択一的な実施形態において、システム１０は、パーソナルコンピュータ、ワークステーション、ＰＡＣＳステーションまたは同じ位置にあるまたはネットワークにわたって分散されているリアルタイムまたは取得後イメージング用の別の装置である。

送信ビームフォーマ１２は、超音波送信器、記憶装置、パルス発生器、アナログ回路、ディジタル回路またはこれらの組み合わせである。送信ビームフォーマ１２は、異なるまたは相対的な振幅、遅延量および／または位相調整を有する複数のチャネルに対して波形を形成する。生成された波に応答してトランスデューサ１４から音波が送信されると、１つまたは複数のビームが形成される。一連の送信ビームが形成されて２次元または３次元領域がスキャンされる。セクタ、Ｖｅｃｔｏｒ(R)、リニア、または別のスキャンフォーマットを使用することができる。同じ領域が複数回スキャンされる。フローイメージングまたはドップライメージングに対して、またストレインイメージングに対して一連のスキャンが使用される。ドップライメージングでは、上記の一連のスキャンには、隣のスキャンラインをスキャンする前に得られる同じスキャンラインに沿った複数のビームを含むことができる。ストレインイメージングに対して、スキャンまたはフレームインタリービングを使用可能である（すなわち、再スキャンする前に全体領域をスキャンする）。択一的な実施形態では、送信ビームフォーマ１２により、より高速なスキャニングに対して平面波または八字波が生成される。

ＨＩＦＵの生成に対し、および／またはＨＩＦＵトランスデューサ２４による波形のエミュレーションおよびイメージングトランスデューサ１４による変位検出に対し、同じ送信ビームフォーマ１２が示されている。択一的な実施形態では、変位イメージングおよび治療に対して異なる送信ビームフォーマ１２が設けられている。例えば、別個の治療システムが使用される。トランスデューサ１４および送信ビームフォーマ１２は、上記の別個の治療システムを操作するため、変位のイメージングに使用される。別の択一的な実施形態では同じトランスデューサ１４が、変位を検出するためおよび治療を適用するために両方に使用される。１つまたは複数のエレメントが治療の送信および診断の送信の両方に使用される。

高密度焦点式超音波トランスデューサ２４は、高密度焦点式超音波治療波形を生成する。ＨＩＦＵトランスデューサ２４は、電気的な波形から音響エネルギを生成するためのアレイである。１次元または多次元のアレイを用いることができる。択一的にはメカニカルフォーカスを有する単一のエレメントが使用される。１つのアレイに対し、相対的な遅延により、上記の音響エネルギが集束される。１つの送信イベントは、音響エネルギを送信することに相応し、この音響エネルギの送信は、異なるエレメントによって実質的に同じ時間において所定のように遅延させることによって行われる。上記の送信イベントにより、組織を治療するための超音波エネルギの１つのパルスが得られる。択一的にはメカニカルフォーカスが上記のアレイに対して設けられる。任意の公知の治療トランスデューサまたは将来開発される超音波トランスデューサ２４を使用することができる。

１実施形態においてＨＩＦＵトランスデューサ２４は、イメージングトランスデューサ１４とは別である。イメージングトランスデューサ１４は、ＨＩＦＵトランスデューサ２４とは別個に移動可能である。イメージングは、治療位置を決定するのに使用される。イメージングトランスデューサ１４は、ＨＩＦＵトランスデューサ２４からの１つのまたは複数の送信に応答するエコー信号を受信する。例えば、治療波形のエミュレーションに応答する信号が受信される。択一的または付加的に２つのトランスデューサ１４，２４には、磁気位置センサなどの空間レジストレーションシステムが含まれる。トランスデューサ１４，２４は互いに接続されていないが、接続して同じケーシング内に配置することなども可能である。

トランスデューサ１４は、圧電素子または容量性メンブレン素子の１次元、１．２５次元、１．５次元、１．７５次元または２次元のアレイである。トランスデューサ１４には、音響エネルギと電気エネルギとの間を変換するための複数の素子が含まれている。受信信号は、トランスデューサ１４の素子に衝突する超音波エネルギ（エコー）に応答して生成される。これらの素子は、送信ビームフォーマ１２および受信ビームフォーマ１６のチャネルに接続される。

受信ビームフォーマ１６には、増幅器、遅延装置および／または位相回転子および１つまたは複数の加算器を有する複数のチャネルが含まれている。各チャネルは、１つまたは複数のトランスデューサ素子に接続される。受信ビームフォーマ１６は、相対的な遅延、位相および／またはアポディゼーションを適用して、各送信に応じて１つまたは複数の受信ビームを形成する。受信ビームフォーマ１６は、上記の受信信号を使用して空間位置を表すデータを出力する。異なる素子からの信号の相対的な遅延および／または位相および和により、ビーム形成が行われる。択一的な実施形態では、受信ビームフォーマ１６は、フーリエ変換または別の変換を用いてサンプルを生成するプロセッサである。

受信ビームフォーマ１６は、フィルタを含むことができ、送信周波数帯域に対して２次高調波または別の周波数帯域における情報を分離するフィルタなどを含むことができる。このような情報には、所望の組織情報、造影剤情報および／またはフロー情報が含まれていることが多い。別の実施形態では受信ビームフォーマ１６には記憶装置またはバッファおよびフィルタまたは加算器が含まれる。２つまたはそれ以上の受信ビームが結合され、これによって２次高調波、３次基本波または他の帯域などの所望の周波数帯域において情報が分離される。

受信ビームフォーマ１６は、空間位置を表すビームを加算したデータを出力する。単一の位置、ラインに沿った位置、領域における位置または体積おける位置に対するデータが出力される。動的な集束を行うことも可能である。これらのデータの目的は異なっていてもよい。例えば、Ｂモードまたは組織データに対し、変位に対するのと異なるスキャンが行われる。択一的には上記のＢモードデータは、変位を求めるためにも使用される。

イメージプロセッサ１８は、Ｂモード検出器、ドップラ検出器、パルス波ドップラ検出器、相関プロセッサ、フーリエ変換プロセッサ、特定用途向け集積回路、汎用プロセッサ、制御プロセッサ、画像プロセッサ、ＦＰＧＡ、ＤＳＰ、アナログ回路、ディジタル回路またはこれらの組み合わせ、あるいはビーム形成された超音波サンプルから情報を検出および処理して表示する現在公知の別の装置または将来開発される装置である。

１実施形態では、プロセッサ１８には１つまたは複数の検出器および別個のプロセッサが含まれている。この別個のプロセッサは、制御プロセッサ、汎用プロセッサ、ＤＳＰ、特定用途向け集積回路、ＦＰＧＡ、ネットワーク、サーバ、プロセッサ群、データパス、それらの組み合わせ、または変位を決定しかつ組織特性を計算するための現在公知の別の装置または将来開発される装置である。例えば、上記の別個のプロセッサにより、図１に示したステップのうちの１つまたは複数の任意の組み合わせが実行される。

プロセッサ１８により、ビーム形成されたサンプルに依存して組織の変位が推定される。任意のタイプの変位を推定することができる。例えば、プロセッサ１８により、剪断波変位情報が検出される。別の実施例では、プロセッサ１８は、弾性率または剪断速度推定が使用される。プロセッサ１８により、受信ビームフォーマ１６からの出力データに依存して情報が求められる。

プロセッサ１８は、上記の組織特性からマッピングされた表示値または画像をディスプレイ２０に出力する。例えば、最大変位、剪断速度、剪断弾性率または別の値が位置毎に求められる。上記の値の大きさにより、色、色相、輝度および／または別の表示特性が変調される。組織変位として表されたビームプロフィールの画像は、これらの変調された表示値から生成される。この画像は、単独で、またはオーバレイされて、または別の画像（例えば、Ｂモード画像）と組み合わされて表示することができる。

図３にはビームプロフィールの画像が示されている。明るい領域は、大きな変位に相当する。上記の変位値は、上記のＨＩＦＵ治療波形のエミュレーションまたは治療波形そのものに起因する変位をベースとしている。このエミュレーション波形または実際の治療波形は、ＨＩＦＵトランスデューサ２４から送信されるため、このビーム位置は治療ビームに相応する。上記の画像は、患者における相対的な位置および／またはビームプロフィールの空間的な拡がりを表す。

１実施形態においてプロセッサ１８は制御プロセッサである。プロセッサ１８により、ＨＩＦＵ治療波形の使用が制御される。上記の組織変位において検出した変化に基づき、プロセッサ１８により、予定された終了の前に治療を止めるか否かが決定される。

変位を求めるため、複数のスキャンまたは測定によって得られたデータが取得されて記憶される。データは、記憶装置２２または別の記憶装置に記憶される。処理の１つまたは複数のステージから得られたデータ、例えば、高周波データ、チャネルデータ、ビームサムデータ、検出データ、ストレインデータ、剪断データ、弾性率データ、剪断弾性率データおよび／または計算値などが記憶される。

プロセッサ１８は、記憶装置２２または別の記憶装置に記憶されている命令にしたがって動作する。プロセッサ１８は、高密度焦点式超音波に対するフィールドバックを行うようにプログラムされている。記憶装置２２は、コンピュータ読み出し可能記憶媒体である。本明細書において説明した処理、方法および/または技術を実現するための命令は、コンピュータ読み出し可能記憶媒体または記憶装置、例えば、キャッシュ、バッファ、ＲＡＭ、リムーバル媒体、ハードディスクドライブまたは別のコンピュータ読み出し可能記憶媒体において提供される。コンピュータ読み出し可能記憶媒体には、種々のタイプの揮発性および不揮発性記憶媒体が含まれる。本明細書において説明したまたは図面に示した機能、ステップまたはタスクは、コンピュータ読み出し可能記憶媒体に記憶されている１つまたは複数の命令のセットに応答して実行される。上記の機能、ステップまたはタスクは特定のタイプの命令セット、記憶媒体、プロセッサまたは処理ストラテジなどには依存しておらず、またソフトウェア、ハードウェア、集積回路、ファームウェア、マイクロコードその他などにより、単独でまたは任意の組み合わせで実行可能である。同様に処理ストラテジは、マルチプロセシング、マルチタスク処理、並列処理などを含んでいてもよい。１つの実施形態において、上記の命令は、リムーバブル媒体装置に記憶されており、ローカルシステムまたはリモートシステムによって読み出される。別の実施形態において上記の命令は、遠隔地に記憶されており、またコンピュータネットワークまたは電話回線を介して転送される。さらに別の実施形態において上記の命令は所定のコンピュータ、ＣＰＵ、ＧＰＵまたはシステム内に記憶されている。

ディスプレイ２０は、２次元画像または３次元表現を表示するＣＲＴ，ＬＣＤ，プロジェクタ、プラズマディスプレイまたは別のディスプレイである。ディスプレイ２０により、上記のＨＩＦＵビームプロフィールを表す１つまたは複数の画像が表示される。上記のＨＩＦＵに関連した変位の空間分布は、上記の画像に示される。上では本発明を複数の実施形態に基づいて説明してきたが、本発明の範囲から逸脱することなく多くの変更および変形を加えることが可能であると理解されたい。したがってここで意図したのは、上記における詳細な説明が制限ではなくむしろ説明上のものであるとみなされるべきであり、また本発明の精神および範囲を定めるのは、同等のものすべてを含む添付の特許請求の範囲であると理解すべきことである。

１０ システム、 １２ 送信ビームフォーマ、 １４ トランスデューサ、 １６ 受信ビームフォーマ、 １８ 画像プロセッサ、 ２０ ディスプレイ、 ２２ 記憶装置、 ２４ ＨＩＦＵトランスデューサ、 ３０〜４６ ステップ

Claims (19)

1. 高密度焦点式超音波に対するフィードバックを行うシステムにおいて、
   該システムは、治療波形およびエミュレーション用の波形を生成する高密度焦点式超音波トランスデューサと、
   前記の高密度焦点式超音波トランスデューサとは別個に移動可能でありかつ反射信号を受信するイメージングトランスデューサとを有しており、前記反射信号は、前記高密度焦点式超音波トランスデューサからのエミュレーション用の波形を患者の組織へ送信した際の、当該患者組織からの反射信号であり、
   前記システムはさらに
   当該の受信信号に応じて、空間位置を表すデータを出力する受信ビームフォーマと、
   当該の出力データに応じて組織の変位を検出しかつ当該の組織の変位に応じてビームプロフィールの画像を生成するプロセッサと、
   当該のビームプロフィールの画像を表示するためのディスプレイとを有することを特徴とする、
   高密度焦点式超音波に対するフィードバックを行うシステム。
2. 前記のビームプロフィールは、前記のエミュレーション用であり、
   前記の画像は、患者における相対的な位置および当該のビームプロフィールの空間的な拡がりを表す、
   請求項１に記載のシステム。
3. 前記のプロセッサにより、前記の組織の変位から検出された変化に応じて高密度焦点式超音波治療波形が制御される、
   請求項１に記載のシステム。
4. 前記のエミュレーション用の波形は、前記の治療波形とは振幅または持続時間が異なるか、または振幅および持続時間の両方が異なる、
   請求項１に記載のシステム。
5. 前記エミュレーション用波形により、治療の熱およびキャビテーションが実質的に生じないようにし、かつ、前記の治療波形により、熱またはキャビテーションが生じる、
   請求項２に記載のシステム。
6. 前記の超音波の送信および受信を複数回行って、当該の複数回の受信によって得られる前記の反射から前記の変位を検出する、
   請求項１に記載のシステム。
7. 前記の組織の変位を検出する際に、前記のエミュレーション用の波形によって生じた剪断波を検出する、
   請求項１に記載のシステム。
8. 前記の剪断波を検出する際に、剪断速度または剪断弾性率を検出する、
   請求項７に記載のシステム。
9. 前記のエミュレーション用の波形のビームプロフィールの画像を生成する際に、前記の組織内における当該のエミュレーション用の波形のビームプロフィールの空間的な拡がりおよび位置を示す多次元画像を生成する、
   請求項１に記載のシステム。
10. さらに、非線形音響伝搬に応じて、前記のビームプロフィールを補正する、
    請求項１に記載のシステム。
11. さらに、
    前記のエミュレーション用の波形の送信および検出を繰り返しにより、前記の治療波形をインタリーブする、
    請求項１に記載のシステム。
12. 前記の変化に応じて、当該の高密度焦点式超音波治療の適用量を変更する、
    請求項１に記載のシステム。
13. 前記のビームプロフィールの画像を生成する際に、前記の変位およびビームプロフィールモデルに応じてビームプロフィールを再構成する、
    請求項１に記載のシステム。
14. コンピュータプログラムにおいて、
    該コンピュータプログラムは、高密度焦点式超音波に対するフィードバックを行うためにプログラムされたプロセッサによって実行される命令を表しており、
    当該命令には、組織の変位に応じて、高密度焦点式超音波治療に対するビーム位置を求めるための命令および当該のビーム位置の画像を生成するための命令が含まれており、
    前記の組織の変位に応じてビーム位置を求めることには、当該のビーム位置に相応するビームによって生じた剪断波を検出することが含まれている
    ことを特徴とする、
    コンピュータプログラム。
15. コンピュータプログラムにおいて、
    該コンピュータプログラムは、高密度焦点式超音波に対するフィードバックを行うためにプログラムされたプロセッサによって実行される命令を表しており、
    当該命令には、組織の変位に応じて、高密度焦点式超音波治療に対するビーム位置を求めるための命令および当該のビーム位置の画像を生成するための命令が含まれており、
    前記のビーム位置を求めることには、前記の組織の異なる空間位置に対して剪断速度または弾性率を求めることが含まれており、
    前記のビーム位置は、別の空間位置よりも剪断速度または弾性率が大きな空間位置に相応する、
    コンピュータプログラム。
16. コンピュータプログラムにおいて、
    該コンピュータプログラムは、高密度焦点式超音波に対するフィードバックを行うためにプログラムされたプロセッサによって実行される命令を表しており、
    当該命令には、組織の変位に応じて、高密度焦点式超音波治療に対するビーム位置を求めるための命令および当該のビーム位置の画像を生成するための命令が含まれており、
    前記のビーム位置の画像を生成することには、前記の組織内でビームプロフィールの空間的な拡がりを示す多次元画像を生成することが含まれている、
    コンピュータプログラム。
17. コンピュータプログラムにおいて、
    該コンピュータプログラムは、高密度焦点式超音波に対するフィードバックを行うためにプログラムされたプロセッサによって実行される命令を表しており、
    当該命令には、組織の変位に応じて、高密度焦点式超音波治療に対するビーム位置を求めるための命令および当該のビーム位置の画像を生成するための命令が含まれており、
    さらに前記の組織の変位に応じてビームプロフィールの画像を生成する際に、非線形音響伝搬に応じて、前記のビームプロフィールを補正する命令が含まれている、
    コンピュータプログラム。
18. 前記のビームの位置を求めることには、高密度焦点式超音波治療波形のエミュレーションを送信することが含まれており、
    当該のエミュレーションは、前記の治療波形とは振幅、持続時間、または振幅および持続時間の双方が異なり、
    前記のエミュレーションにより、治療の熱およびキャビテーションが実質的に生じないようにし、また超音波スキャニングによって前記の変位が求められ、ここで変位の位置が前記のエミュレーションによって生じる、
    請求項１４から１７までのいずれか１項に記載のコンピュータプログラム。
19. さらに、
    前記の変位に応じて、高密度焦点式超音波治療による前記の組織における変化を求めるための命令と、
    当該の変化の関数として、当該の高密度焦点式超音波治療の適用量を変更する命令とを有する、
    請求項１４から１８までのいずれか１項に記載のコンピュータプログラム。