(57)【要約】  本開示は、標的組織の弾性を決定する超音波システム及び方法を説明する。システムは、剛性が増加した領域を含みうる組織に向けて送信された超音波パルスに応答するエコーを取得する超音波トランスデューサを含むことができる。システムはまた、組織内にプッシュパルスを送信し、これにより、剛性が増加した領域内にせん断波を生成するようにトランスデューサを制御するビームフォーマを含むことができる。ビームフォーマは、プッシュパルスに隣接して追跡パルスを放出するようにトランスデューサを制御することができる。システムはさらに、せん断波の変位振幅を決定し、その振幅に基づいて、組織の定性的組織弾性マップを生成するプロセッサを含むことができる。プロセッサは、定性的マップを同じ組織の定量的マップと組み合わせ、その組み合わせに基づいて、剛性が増加した領域の境界を決定することができる。
2021522004.tif 000002
【特許請求の範囲】
【請求項１】
せん断波撮像用の超音波撮像システムであって、
標的組織に向けて送信された超音波パルスに応答するエコーを取得する超音波トランスデューサと、
ビームフォーマと、
前記ビームフォーマと通信するプロセッサと、
を含み、
前記ビームフォーマは、プッシュパルスに応答して追跡パルスを、前記超音波トランスデューサから送信するように構成され、前記プッシュパルスは、前記標的組織内にせん断波を生成し、前記追跡パルスは、１つ又は複数の場所において、前記せん断波と交差するように空間的に計画され、
前記ビームフォーマは、前記追跡パルスが前記せん断波と交差した場所からのエコー信号を、前記超音波トランスデューサから受信するように構成され、
前記プロセッサは、受信した前記エコー信号から生成される追跡エコーデータを記憶し、前記追跡エコーデータに応答して、前記標的組織を通って伝搬する前記せん断波の変位振幅を決定し、決定された前記変位振幅に基づいて、前記標的組織の定性的組織弾性マップを生成する、
超音波撮像システム。
【請求項２】
前記プロセッサは、決定された前記変位振幅を参照変位振幅と比較することによって、前記定性的組織弾性マップを生成する、請求項１に記載の超音波撮像システム。
【請求項３】
決定された前記変位振幅は、前記標的組織内の２つ以上の横方向に間隔を置いた点において決定されていて、前記参照変位振幅は、参照組織内の２つ以上の横方向に間隔を置いた点において決定されるか、又は、シミュレートされたモデルから数値的に決定される、請求項２に記載の超音波撮像システム。
【請求項４】
前記参照組織は、前記標的組織のファントムモデルを含む、請求項３に記載の超音波撮像システム。
【請求項５】
前記参照組織は、前記標的組織に対応する組織タイプで、剛性が増加した領域を欠く患者サンプルを含む、請求項３に記載の超音波撮像システム。
【請求項６】
前記プロセッサはさらに、
前記標的組織内の前記横方向に間隔を置いた点の連続する対における、決定された前記変位振幅を比較することによって、変位振幅非相関化値を決定し、
前記参照組織内の前記横方向に間隔を置いた点の連続する対における、決定された前記変位振幅を比較することによって、参照変位振幅非相関化値を決定し、
前記変位振幅非相関化値を前記参照変位振幅非相関化値と比較し、
前記比較に基づいて、前記定性的組織弾性マップを生成する、
請求項３に記載の超音波撮像システム。
【請求項７】
前記超音波トランスデューサは、焼灼デバイスに結合され、前記焼灼デバイスは、剛性が増加した領域又は前記剛性が増加した領域を含むより大きな領域を焼灼する、請求項１に記載の超音波撮像システム。
【請求項８】
前記超音波トランスデューサ、前記ビームフォーマ及び前記プロセッサは、前記焼灼デバイスと同時に動作する、請求項７に記載の超音波撮像システム。
【請求項９】
前記定性的組織弾性マップを表示するユーザーインターフェースをさらに含む、請求項１に記載の超音波撮像システム。
【請求項１０】
前記参照変位振幅は、参照マップから導出される、請求項２に記載の超音波撮像システム。
【請求項１１】
複数の参照マップを記憶するメモリをさらに含む、請求項１０に記載の超音波撮像システム。
【請求項１２】
前記標的組織は、熱的病変からなる剛性が増加した領域を含む、請求項１に記載の超音波撮像システム。
【請求項１３】
標的組織に向けて送信された超音波パルスに応答する超音波エコーを取得するステップと、
前記標的組織内にせん断波を生成するように、前記標的組織内へとプッシュパルスを送信するステップと、
１つ又は複数の場所で、前記せん断波と交差するように空間的に計画された追跡パルスを送信するステップと、
前記追跡パルスが前記せん断波と交差した場所からエコー信号を受信するステップと、
受信した前記エコー信号から生成される追跡エコーデータを記憶するステップと、
前記追跡エコーデータに基づいて、前記標的組織を通って伝播する前記せん断波の変位振幅を決定するステップと、
決定された前記変位振幅に基づいて、前記標的組織の定性的組織弾性マップを生成するステップと、
を含む、せん断波撮像方法。
【請求項１４】
決定された前記変位振幅に基づいて、前記標的組織の前記定性的組織弾性マップを生成するステップは、決定された前記変位振幅を参照変位振幅と比較するステップを含む、請求項１３に記載の方法。
【請求項１５】
決定された前記変位振幅は、前記標的組織内の２つ以上の横方向に間隔を置いた点において決定されていて、前記参照変位振幅は、参照組織内の２つ以上の横方向に間隔を置いた点において決定されるか、又は、シミュレートされたモデルから数値的に決定される、請求項１４に記載の方法。
【請求項１６】
前記標的組織内の前記横方向に間隔を置いた点の連続する対における、決定された前記変位振幅を比較することによって、変位振幅非相関化値を決定するステップと、
前記参照組織内の前記横方向に間隔を置いた点の連続する対における、決定された前記変位振幅を比較することによって、参照変位振幅非相関化値を決定するステップと、
前記変位振幅非相関化値を前記参照変位振幅非相関化値と比較するステップと、
前記比較に基づいて、前記定性的組織弾性マップを生成するステップと、
をさらに含む、請求項１５に記載の方法。
【請求項１７】
前記定性的組織弾性マップを、ユーザーインターフェース上に表示するステップをさらに含む、請求項１３に記載の方法。
【請求項１８】
前記参照変位振幅は、参照マップから導出される、請求項１４に記載の方法。
【請求項１９】
前記標的組織は、熱的病変からなる剛性が増加した領域を含む、請求項１３に記載の方法。
【請求項２０】
実行されると、超音波撮像システムのプロセッサに、請求項１３から１９のいずれか一項に記載の方法を行わせる、実行可能命令を含む、非一時的コンピュータ可読媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
[001]  本開示は、組織弾性を決定するための超音波システム及び方法に関する。特定の実装形態には、せん断波の振幅減衰及び非相関化に基づく組織弾性の決定が含まれる。
【背景技術】
【０００２】
[002]  高周波焼灼（ＲＦＡ）は、現在、世界で第２位の癌による死亡原因となっている肝臓癌の治療処置の最も広く使用されている方式である。ＲＦＡは、侵襲が最小限であり、腫瘍部位に挿入された焼灼電極、針、又はタイン（くし歯）を使用して、腫瘍を凝固壊死のポイントまで加熱することを伴う。ＲＦＡの標的となる腫瘍の明確な境界線の描写は、正確に癌組織を標的化するために重要である。現在のＲＦＡ治療プロトコルは、多くの場合、楕円焼灼ボリューム予測法を実装している。しかし、このような方法は、あまりにも単純である。結果として、実際の治療ボリュームは、予測されたボリュームから大幅に逸脱する場合があり、健康な組織の、標的から外れた焼灼及び／又は腫瘍組織の不完全な焼灼につながる。
【０００３】
[003]  超音波撮像は、ＲＦＡ処置中の誘導に一般的に使用される。特に、超音波せん断波エラストグラフィ撮像（ＳＷＩ）が、組織弾性を測定することによって焼灼の範囲を推定するために使用されてきた。超音波ＳＷＩは、トランスデューサから組織に「プッシュパルス」（音響放射力として知られる現象）を送信し、それによって内部を横方向に伝播するせん断波を生成することによって、肝臓組織を含むさまざまな組織の局所的な剛性レベルを決定することができる。次に、トランスデューサによって放出された追跡パルスを使用して、せん断波が伝播するときのその速度を測定することができる。これは、多くの場合、組織の剛性に比例する。例えば、同じプッシュパルスを使用して各組織タイプ内にせん断波を生成すると仮定すると、軟組織内のせん断波速度は、通常、硬組織内のせん断波速度よりも遅い。焼灼された壊死組織は通常、未処置の組織よりもはるかに硬いため、焼灼された組織の境界は、理論的にはせん断波速度に基づいて決定することができる。
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかし、現在のＳＷＩモダリティ及び関連する組織再構成技術では、ＲＦＡによって生じた熱的病変の剛性が高いことや、治療部位に剛性の焼灼電極が存在することにより、このような境界を確実に判断することができない。これらの要因は、ＳＷＩシステムによって生成されるせん断波の振幅を減少し、さまざまな組織特徴から反射する外部信号の中でこのような波を検出することがさらに困難になり、その結果、焼灼ゾーン内の信号対ノイズ比（ＳＮＲ）が低くなる。したがって、ＳＷＩを介して調べられる焼灼ゾーンの弾性推定値及び定量的弾性マップは、広く不正確であり、信頼性がない。焼灼治療及び組織マッピングの精度を高めるために、改良された組織弾性測定及び焼灼モニタリング技術が必要である。
【課題を解決するための手段】
【０００５】
[004]  本開示は、せん断波超音波撮像を介して標的組織の弾性を決定するシステム及び方法を説明する。標的組織は、剛性が増加した領域を含み得、局在的で、サイズが様々であり得る。いくつかの例では、領域は、焼灼処置によって生じた熱的病変を含む。せん断波撮像は、既存のシステムにおいて、組織弾性の情報を決定するために利用されてきたが、このようなシステムは、熱焼灼によって生じた硬い病変など、剛性の高い領域で生成される低振幅のせん断波によって妨げられることがよくある。硬組織の位置及び境界を正確に決定するために、本明細書のシステムは、伝播するせん断波の変位振幅を決定し、この振幅に基づいて、硬い領域の位置を正確かつ精密に示す定性的組織弾性マップを導出するように独自に構成されている。硬い領域の境界を明確に定め、計算効率を向上させるために、本明細書のシステムはまた、同じ領域の定量的組織マップからノイズを除去する。特定の例には、定性的組織弾性マップを定量的組織弾性マップと組み合わせ、結果として得られる組み合わせに基づいて、硬い領域の正確な境界を定めるプロセッサが含まれる。さまざまな実施形態では、２つのマップタイプを組み合わせることには、最初に、定性的組織弾性マップから組織弾性の等高線プロットを生成し、次に、等高線プロットを定量的マップにオーバーレイすることが含まれる。剛性が増加した領域に最適に適合する等高線が決定され、選択される。
【０００６】
[005]  本開示の原理によれば、超音波撮像システムは、標的組織に向けて送信された超音波パルスに応答するエコーを取得する超音波トランスデューサを含む。システムはまた、プッシュパルスに応答して超音波トランスデューサから追跡パルスを送信するように構成されるビームフォーマを含み、プッシュパルスは、標的組織内にせん断波を生成し、追跡パルスは、１つ又は複数の場所において、せん断波と交差するように空間的に計画されている。ビームフォーマはまた、追跡パルスがせん断波と交差した場所からのエコー信号を超音波トランスデューサから受信する。システムはまた、ビームフォーマと通信し、受信したエコー信号から生成された追跡エコーデータを記憶し、追跡エコーデータに応答して、標的組織を通って伝搬するせん断波の変位振幅を決定し、決定された変位振幅に基づいて、標的組織の定性的組織弾性マップを生成するプロセッサを含む。
【０００７】
[006]  いくつかの例では、プロセッサは、決定された変位振幅を参照変位振幅と比較することによって、定性的組織弾性マップを生成する。いくつかの実施形態では、決定された変位振幅は、標的組織内の２つ以上の横方向に間隔を置いた点において決定されていて、参照変位振幅は、参照組織内の２つ以上の横方向に間隔を置いた点において決定されるか、又は、シミュレートされたモデルから数値的に決定される。いくつかの実施形態では、参照組織は、標的組織のファントムモデル、又は、標的組織に対応する組織タイプで、剛性が増加した領域を欠く患者サンプルを含む。
【０００８】
[007]  いくつかの例では、プロセッサはさらに、標的組織内の横方向に間隔を置いた点の連続する対における、決定された変位振幅を比較することによって、変位振幅非相関化値を決定し、参照組織内の横方向に間隔を置いた点の連続する対における、決定された変位振幅を比較することによって、参照変位振幅非相関化値を決定し、変位振幅非相関化値を参照変位振幅非相関化値と比較し、比較に基づいて、定性的組織弾性マップを生成する。いくつかの実施形態では、超音波トランスデューサは、焼灼デバイスに結合され、焼灼デバイスは、剛性が増加した領域又は剛性が増加した領域を含むより大きな領域を焼灼する。いくつかの例では、超音波トランスデューサ、ビームフォーマ、及びプロセッサは、焼灼デバイスと同時に動作する。例示的なシステムはさらに、定性的組織弾性マップを表示するユーザーインターフェースを含む。いくつかの実施形態では、参照変位振幅は、参照マップから導出される。いくつかの例では、システムはまた、複数の参照マップを記憶するメモリを含む。いくつかの実施形態では、標的組織は、熱的病変からなる剛性が増加した領域を含む。
【０００９】
[008]  本開示の原理によれば、せん断波撮像用の超音波撮像システムは、標的組織に向けて送信された超音波パルスに応答するエコーを取得する超音波トランスデューサを含む。システムはまた、プッシュパルスに応答して超音波トランスデューサから追跡パルスを送信するように構成されたビームフォーマを含み、ここで、プッシュパルスは、標的組織内にせん断波を生成し、追跡パルスは、１つ又は複数の場所で、せん断波と交差するように空間的に計画される。ビームフォーマは、追跡パルスがせん断波と交差した場所からのエコー信号を超音波トランスデューサから受信する。システムはまた、ビームフォーマと通信するプロセッサを含む。プロセッサは、受信したエコー信号に基づいて、標的組織の定性的組織弾性マップを生成し、受信したエコー信号に基づいて、標的組織の定量的組織弾性マップを生成し、定性的組織弾性マップを定量的組織弾性マップと組み合わせることによって、標的組織内の剛性が増加した領域の境界を決定する。
【００１０】
[009]  いくつかの例では、プロセッサはさらに、定性的組織弾性マップから均一剛性の領域の境界を定める等高線プロットを導出する。いくつかの実施形態では、定性的組織弾性マップを定量的組織弾性マップと組み合わせることは、等高線プロットを定量的組織弾性マップにオーバーレイすることを含む。いくつかの実施形態では、プロセッサは、剛性が増加した領域の周りに最適適合の等高線をさらに決定することによって、剛性が増加した領域の境界を決定する。いくつかの実装形態では、最適適合の等高線は、ユーザーによって選択可能な基準のセットに基づいて、最適適合の２つ以上の候補等高線から選択される。いくつかの例では、基準のセットは、最適適合の等高線が組織の最大領域を画定する等高線を含むように、過剰封入バイアスを含む。いくつかの実施形態では、プロセッサはさらに、剛性が増加した領域の決定された境界に基づいてハイブリッド組織弾性マップを生成する。いくつかの実施形態では、プロセッサは、決定された境界の外側のハイブリッドマップの少なくとも１つの領域をマスキングすることによってハイブリッドマップを生成する。例示的なシステムはさらに、ハイブリッドマップを表示するユーザーインターフェースを含む。いくつかの実装形態では、定量的組織弾性マップは、組織弾性勾配マップを含む。いくつかの実施形態では、プロセッサは、等高線プロット内の各等高線に沿って勾配値を合計し、各等高線の勾配値の合計を各等高線の長さで割り、平均勾配が最大の等高線を選択することによって、最適適合の等高線を決定する。いくつかの例では、超音波トランスデューサは、焼灼デバイスに結合される。焼灼デバイスは、剛性が増加した領域又は剛性が増加した領域を含むより大きな領域を焼灼する。いくつかの実施形態では、剛性が増加した領域は、熱的病変を含む。
【００１１】
[010]本開示の原理によれば、せん断波撮像方法は、標的組織に向けて送信された超音波パルスに応答する超音波エコーを取得するステップと、プッシュパルスを標的組織に送信して、標的組織内にせん断波を生成するステップと、１つ又は複数の場所でせん断波と交差するように空間的に計画された追跡パルスを送信するステップと、追跡パルスがせん断波と交差した場所からエコー信号を受信するステップと、受信したエコー信号から生成された追跡エコーデータを記憶するステップと、追跡エコーデータに基づいて、標的組織を通って伝播するせん断波の変位振幅を決定するステップと、決定された変位振幅に基づいて、標的組織の定性的組織弾性マップを生成するステップとを含む。
【００１２】
[011]  いくつかの例では、決定された変位振幅に基づいて、標的組織の定性的組織弾性マップを生成するステップは、決定された変位振幅を参照変位振幅と比較するステップを含む。いくつかの実施形態では、決定された変位振幅は、標的組織内の２つ以上の横方向に間隔を置いた点において決定されていて、参照変位振幅は、参照組織内の２つ以上の横方向に間隔を置いた点において決定されるか、又はシミュレートされたモデルから数値的に決定される。いくつかの例では、方法はさらに、標的組織内の横方向に間隔を置いた点の連続する対において決定された変位振幅を比較することによって、変位振幅非相関化値を決定するステップと、参照組織内の横方向に間隔を置いた点の連続する対において決定された変位振幅を比較することによって、参照変位振幅非相関化値を決定するステップと、変位振幅非相関化値を参照変位振幅非相関化値と比較するステップと、比較に基づいて、定性的組織弾性マップを生成するステップとを含む。
【００１３】
[012]  いくつかの実施形態では、方法はさらに、定性的組織弾性マップをユーザーインターフェース上に表示するステップを含む。いくつかの例では、参照変位振幅は、参照マップから導出される。いくつかの実施形態では、標的組織は、熱的病変からなる剛性が増加した領域を含む。
【００１４】
[013]  本開示の原理によれば、せん断波撮像方法は、標的組織に向けて送信された超音波パルスに応答する超音波エコーを取得するステップと、プッシュパルスを標的組織に送信して、標的組織内にせん断波を生成するステップと、１つ又は複数の場所でせん断波と交差するように空間的に計画された追跡パルスを送信するステップと、追跡パルスがせん断波と交差した場所からエコー信号を受信するステップと、受信したエコー信号から生成された追跡エコーデータを記憶するステップと、受信したエコー信号に基づいて、標的組織の定性的組織弾性マップを生成するステップと、受信したエコー信号に基づいて、標的組織の定量的組織弾性マップを生成するステップと、定性的組織弾性マップを定量的組織弾性マップと組み合わせることによって、標的組織内の剛性が増加した領域の境界を決定するステップとを含む。
【００１５】
[014]  いくつかの例では、方法はさらに、定性的組織弾性マップから均一剛性の領域の境界を定める等高線プロットを導出するステップを含む。いくつかの実施形態では、定性的組織弾性マップを定量的組織弾性マップと組み合わせるステップは、等高線プロットを定量的組織弾性マップにオーバーレイするステップを含む。いくつかの例では、剛性が増加した領域の境界を決定するステップは、剛性が増加した領域の周りの最適適合の等高線を決定するステップを含む。いくつかの実施形態では、決定された境界の外側の少なくとも１つの領域をマスキングすることによって、剛性が増加した領域の決定された境界に基づいてハイブリッド組織弾性マップが生成される。いくつかの例では、剛性が増加した領域は、熱的病変を含む。
【００１６】
[015]  本明細書に説明される方法のいずれか又はそのステップは、実行されると、医用撮像システムのプロセッサに、本明細書に具体化される方法又はステップを行わせることができる実行可能命令を含む非一時的コンピュータ可読媒体に具体化することができる。
【図面の簡単な説明】
【００１７】
【図１】[016]図１は、本開示の原理に従って構成された超音波撮像システムのブロック図である。
【図２】[017]図２は、本開示の原理に従って組織焼灼システムと同時に動作する超音波せん断波撮像システムの写真である。
【図３】[018]図３は、図２に示されるシステムによって生成された焼灼ゾーン及び隣接する関心領域の超音波画像及び対応する組織弾性マップである。
【図４Ａ】[019]図４Ａは、本開示の原理に従って第１の時点において構成された参照せん断波振幅変位マップである。
【図４Ｂ】[020]図４Ｂは、本開示の原理に従って第２の時点において構成された参照せん断波振幅変位マップである。
【図４Ｃ】[021]図４Ｃは、本開示の原理に従って構成された参照せん断波振幅変位グラフである。
【図４Ｄ】[022]図４Ｄは、本開示の原理に従って構成された参照せん断波非相関化グラフである。
【図５Ａ】[023]図５Ａは、本開示の原理に従って第１の時点において構成されたせん断波振幅変位マップである。
【図５Ｂ】[024]図５Ｂは、本開示の原理に従って第２の時点において構成されたせん断波振幅変位マップである。
【図５Ｃ】[025]図５Ｃは、本開示の原理に従って構成されたせん断波振幅変位グラフである。
【図５Ｄ】[026]図５Ｄは、本開示の原理に従って構成されたせん断波非相関化グラフである。
【図６Ａ】[027]図６Ａは、本開示の原理に従って構成されたせん断波変位振幅再構成マップである。
【図６Ｂ】[028]図６Ｂは、本開示の原理に従って構成されたせん断波非相関化再構成マップである。
【図７】[029]図７は、本開示の原理に従って実行される方法を示すフローチャートである。
【図８Ａ】[030]図８Ａは、本開示の原理に従って構成されたせん断波非相関化再構成マップである。
【図８Ｂ】[031]図８Ｂは、図８Ａのせん断波非相関化再構成マップを使用して生成された等高線マップである。
【図８Ｃ】[032]図８Ｃは、図８Ｂの等高線マップでオーバーレイされた定量的弾性マップである。
【図９Ａ】[033]図９Ａは、本開示の原理による、病変の境界を定める候補等高線に適用される最適適合メトリックのグラフである。
【図９Ｂ】[034]図９Ｂは、図９Ａのグラフにおいて特定された最適な等高線でオーバーレイされたせん断波非相関化再構成マップである。
【図９Ｃ】[035]図９Ｃは、図９Ｂに示される最適な等高線でオーバーレイされた定量的弾性マップである。
【図１０Ａ】[036]図１０Ａは、本開示の原理による、病変の境界を定める候補等高線に適用される最適適合メトリックのグラフである。
【図１０Ｂ】[037]図１０Ｂは、候補の最適な等高線でオーバーレイされたせん断波非相関化再構成マップである。
【図１０Ｃ】[038]図１０Ｃは、図１０Ｂに示される候補の最適な等高線でオーバーレイされた定量的弾性マップである。
【図１０Ｄ】[039]図１０Ｄは、候補の最適な等高線でオーバーレイされた組織弾性勾配マップである。
【図１１】[040]図１１は、本開示の原理に従って実行される方法を示すフローチャートである。
【発明を実施するための形態】
【００１８】
[041]  特定の実施形態の以下の説明は、本質的に単なる例示であり、本発明又はその用途又は使用を制限することを決して意図するものではない。本システム及び方法の実施形態の以下の詳細な説明では、本明細書の一部を形成し、また、説明されるシステム及び方法を実施することができる特定の実施形態を例示として示す添付図面を参照する。これらの実施形態は、当業者が現在開示されているシステム及び方法を実践することを可能にするのに十分詳細に説明され、また、他の実施形態を利用することができ、構造的及び論理的変更を本システムの精神及び範囲から逸脱することなく行うことができることを理解されたい。さらに、明確にするために、いくつかの特徴の詳細な説明は、それらが当業者に明らかである場合には、本システムの説明を曖昧にしないように議論されない。したがって、以下の詳細な説明は、限定的な意味で解釈されるべきではなく、本システムの範囲は、添付の特許請求の範囲によってのみ定義される。
【００１９】
[042]  本技術はまた、本実施形態による方法、装置（システム）及び／又はコンピュータプログラム製品のブロック図及び／又はフローチャート図を参照して以下に説明される。ブロック図及び／又はフローチャート図のブロック、ならびにブロック図及び／又はフローチャート図のブロックの組み合わせは、コンピュータ実行可能命令によって実装され得ることが理解される。これらのコンピュータ実行可能命令は、汎用コンピュータ、特殊目的コンピュータ、及び／又は他のプログラム可能なデータ処理装置のプロセッサ、コントローラ、又は制御ユニットに提供されて、コンピュータ及び／又は他のプログラム可能なデータ処理装置のプロセッサを介して実行するこのような命令が、ブロック図及び／又はフローチャートブロック（複数可）で指定された機能／行為を実施するための手段を作成するようにマシンを作成する。
【００２０】
[043]  本明細書において、焼灼処置の精度を向上させるやり方で組織焼灼の正確でリアルタイムのモニタリング及び表示を提供する超音波ベースのＳＷＩシステムが提供される。既存の焼灼モニタリングシステムは、ＳＷＩを利用して、組織の領域を通過するせん断波の伝播速度を測定することにより、組織弾性を推定することができる。ピーク到達時間（time-to-peak）再構成とも呼ばれるこのような飛行時間型アプローチは、通常、超音波トランスデューサから放出される横方向に間隔を置いた追跡ビームを介して領域全体で検出された波の遅延を測定することによってせん断波速度を推定する。ＲＦＡに関連する複数の要因により、このような手法の精度が低下する。例えば、飛行時間型アプローチは、ＲＦＡによって生成された熱的病変を通過するせん断波の振幅が小さいことによる低信号対ノイズ比によって妨げられる。熱的病変は、焼灼後の組織の乾燥状態によって、周囲の健康又は癌性の組織よりも６〜８倍硬い場合がある。せん断波は、このような硬い組織をすばやく通過し、波の振幅を減少させる。焼灼処置が続くにつれて、熱的病変のボリュームが増加し、低振幅のせん断波によって生成される信頼性の低い定量的データの割合も増加する。この結果、残っている腫瘍組織の量が最小である処置の終わり近くで、焼灼モニタリングの精度が最低になる可能性がある。したがって、継続焼灼の精度は、処置の最も重要な時点、すなわち、最小量の癌性組織が残っているときに、低い点に達する可能性があり、それによって、健康な組織の、標的から外れた焼灼及び癌性組織の不完全な焼灼のリスクを増幅する。ＳＷエラストグラフィは、自由場のせん断波伝播を想定しており、ＲＦＡで使用される固体焼灼電極又はタイン（くし歯）の存在は、この想定を無効にし、電極がせん断波撮像の視野内に直接置かれていなくても、伝播するせん断波の振幅をさらに減少させる可能性がある。これらの要因の累積的な影響は、一貫して低いＳＮＲであり、中程度の解像度又は精度でさえ組織弾性マップの解明を妨げる。これらの問題を解決するために、本明細書に開示されるシステムは、硬い焼灼針に近接している場合でも高感度で病変境界を識別し、焼灼ゾーン全体にわたって包括的でリアルタイムの焼灼モニタリング及び表示を提供することによって超音波ＳＷＩを向上させる。例示的なシステムは、治療部位内の異なる空間位置においてせん断波振幅変位及び相関値を検出することによって、定性的な組織弾性マップを生成するように具体的に構成することができる。システムはまた、病変境界の描写の精度を高めるやり方で、定量的組織弾性マップを定性的組織再構成マップと組み合わせるように構成される。本明細書に説明される実施例のほとんどは、焼灼によって生じた組織病変の位置及び境界の決定に関連しているが、当業者は、開示されたシステムを利用して、病変が存在しているかいないかに関わらず、硬組織の領域や、プッシュパルスの受信に応答して低振幅のせん断波を生成する組織の領域を含む多くの組織タイプを調べることができることを理解すべきである。
【００２１】
[044]  図１は、焼灼ゾーンを横切ってせん断波エラストグラフィを実行し、正確な組織弾性データを生成することによって焼灼ゾーンの境界を決定し、生成されたデータに基づいて、焼灼ゾーンの定性的弾性マップを生成する例示的な超音波システム１００を示す。システム１００は、超音波信号を送受信する超音波センサーアレイ１１４を収容する超音波プローブ１１２を含む超音波取得ユニット１１０を含む。アレイ１１４は、腫瘍又は硬組織封入体などの１つ又は複数の組織異常１２０を含む標的領域１１８に、高振幅プッシュパルス１１６を放出する。追加の又は代替の実施形態では、プッシュパルス１１６は、アレイ１１４以外のアレイによって生成される。例えば、いくつかの例では、１つのアレイが音響放射力（ＡＲＦ）を加えるために使用され、異なるアレイが撮像のために使用される。さらに他の例では、組織は、例えば、外部の機械力を加えるように構成された外部の機械的アクチュエータを使用して、機械的に刺激される。治療段階に応じて、凝固領域又は熱的病変とも呼ばれることがある少なくとも１つの焼灼ゾーン１２２が存在する。標的領域１１８は、ヒトの肝臓、膵臓、腎臓、肺、心臓、又は脳を含むがこれらに限定されない臓器、又は組織の領域、例えば、筋肉組織を含む。アレイ１１４はまた、複数の追跡パルス又はビーム１２４を標的ゾーン１１８に送信して、プッシュパルス１１６によって生成されたせん断波１１９の伝播を検出する。追跡パルス１２４は、プッシュパルス１１６に隣接して送信することができ、いくつかの例では、プッシュパルスに対して横方向に間隔を空けられている。いくつかの実施形態では、追跡パルス１２４は、例えば、線形プローブが追跡パルスを放出するために利用される場合、プッシュパルス１１６に平行である。他の例では、追跡パルス１２４は、プッシュパルス１１６と平行に送信されなくてもよい。例えば、湾曲したプローブが、追跡パルスを半径方向に、それらの間に角距離を置いて送信することができる。このようなパルスは、デカルト空間では平行ではない場合があるが、極座標フレーム又は円筒座標フレームでは同じ方向に送信される。アレイ１１４は、送信／受信（Ｔ／Ｒ）スイッチ１３０を介して送信ビームフォーマ１２６及びマルチライン受信ビームフォーマ１２８に結合されている。ビームフォーマ１２６、１２８による送信及び受信の調整は、ビームフォーマコントローラ１３２によって制御することができる。動作中、送信ビームフォーマ１２６は、アレイ１１４を制御して、１つ又は複数の、例えば、一連のプッシュパルス１１６を、焼灼ゾーン１２２（焼灼処置がまだ始まっていない場合、予想焼灼ゾーンを含んでもよい）に隣接する標的領域１１８に送信することができる。マルチライン受信ビームフォーマ１２８は、アレイ１１４によって受信されるエコー信号１３４であって、信号プロセッサ１３６によってフィルタリング、ノイズ低減などによって処理されるエコー信号１３４の空間的に異なる受信ライン（Ａライン）を生成する。いくつかの実施形態では、取得ユニット１１０の構成要素は、超音波エコー１３４から複数の連続する超音波画像フレーム１３８を生成する。
【００２２】
[045]  システム１００はまた、焼灼ゾーン１２２を通って伝播するせん断波１１９の変位振幅を決定する定性的処理モジュール１４０などの１つ又は複数のプロセッサを含む。実施形態では、変位振幅は、焼灼ゾーン１２２内の２つ以上の横方向に間隔を置いた点において検出することができ、これにより、プッシュパルスから離れた組織全体の変位振幅の減衰を決定する。次に、定性的処理モジュール１４０は、振幅減衰を、選択された参照マップ１４２から導出された参照振幅減衰と比較する。図１に示されるように、参照マップ１４２は、定性的処理モジュール１４０と通信可能に結合されたメモリ構成要素１４４から選択的に抽出する。メモリ構成要素１４４は、ライブラリ又はデータベース１４６に複数の参照マップを記憶する。決定された振幅減衰と参照振幅減衰との比較に基づいて、定性的処理モジュール１４０は、プッシュパルス１１６によって生成されるせん断波１１９の振幅が小さいにもかかわらず、焼灼ゾーン１２２内及びその近くの下にある組織の弾性を高精度で再構成する。
【００２３】
[046]  さまざまな実施形態では、システム１００はまた、定性的処理モジュール１４０と結合された表示プロセッサ１４８を、表示プロセッサの出力を表示するユーザーインターフェース１５０とともに含む。表示プロセッサ１４８は、画像フレーム１３８及び定性的せん断波再構成マップ１５４から超音波画像１５２を生成する。以下に説明するように、定性的せん断波再構成マップ１５４は、せん断波変位振幅再構成マップ又はせん断波非相関化再構成マップを含み、これらは両方とも、組織焼灼ゾーン内及びその近くの組織弾性の定性的表現を具体化する。ユーザーインターフェース１５０は、超音波スキャン及び／又は焼灼処置が実行されているときに、画像１５２及び定性的再構成マップ１５４をリアルタイムで表示することができ、また、このような処置の前、最中、又は後にいつでもユーザー入力１５６を受け取ることができる。いくつかの例では、ユーザーインターフェース１５０に表示される超音波画像及び／又はマップは、超音波ＳＷＩスキャン中、また、いくつかの実施形態では、焼灼処置中に、データ取得ユニット１１０によって受信及び処理される取得フレームごとに更新される。
【００２４】
[047]  いくつかの実装形態では、システム１００はまた、定性的せん断波再構成マップ１５４及び定量的弾性マップ１６０を受信し、２つのマップに基づいて、焼灼によって生成された熱的病変または硬組織の領域を含む病変の決定された境界を精緻化する境界モジュール１５８を含む。簡単にするために、本明細書では病変境界を参照するが、追加の組織境界、例えば、臓器境界又は一部の癌性病変などの剛性が増加した局所領域もまた、開示されるシステム及び方法を介して識別可能であることが当業者に理解されるべきである。病変の境界を精緻化することにより、境界モジュール１５８は、以前は病変の境界をぼやけさせていた信号ノイズを低減又は排除し、これにより、システム１００の精度及び精密さを高めることができる。境界モジュール１５８は、病変の実際の物理的境界を正確に追跡する絶対閾値を決定するように具体的に構成される。定量的弾性マップ１６０は、ＳＷＩ中にデータ取得ユニット１１０によって生成された超音波画像フレーム１３８から生成される。特定の実施形態では、定量的処理モジュール１６２が、画像フレーム１３８から定量的弾性マップ１６０を生成する。定量的剛性マップは、メモリ構成要素１６２に記憶される。さらに示されるように、ユーザーインターフェース１５０と結合される境界表示プロセッサ１６４が含まれる。表示プロセッサ１６４は、単独で又は境界モジュール１５８及び／又はユーザーインターフェース１５０とともに構成されて、病変１６８及び病変境界の外側の少なくとも１つのマスクされた領域１７０を含むハイブリッドマップ１６６を生成する。
【００２５】
[048]  図１に示されるシステム１００の構成は、異なっていてもよい。例えば、システム１００は、携帯型であっても固定式であってもよい。さまざまな携帯機器、例えば、ラップトップ、タブレット、スマートフォンなどを使用して、システム１００の１つ又は複数の機能を実装することができる。このようなデバイスを組み込んだ例では、超音波センサーアレイ１１４は、例えば、ＵＳＢインターフェースを介して接続可能である。システム１００は、ＲＦＡ、マイクロ波焼灼、又は高密度焦点式超音波焼灼を含むがこれらに限定されない焼灼療法を実行する焼灼デバイス１７２と同時に動作することができる。いくつかの例では、焼灼デバイス１７２は、焼灼ゾーン１２２を構成する組織に挿入される針、タイン（くし歯）又はプローブを含み、デバイスは、いくつかの実施形態では、局所的な熱を向ける。システム１００の１つ又は複数の構成要素は、焼灼デバイス１７２と通信可能に、動作可能に、及び／又は物理的に結合し、これにより、デバイス及びシステムが、他方の動作に応答できる。例えば、いくつかの実施形態では、焼灼デバイス１７２は、定性的処理モジュール１４０及び／又は境界モジュール１５８によって生成された組織弾性マップに基づいて、焼灼ゾーンのサイズ及び／又は位置を調整する。したがって、システム１００は、焼灼処置中又はその前に、焼灼デバイス１７２及び／又はその操作者に病変の位置及び／又は境界に関する情報を提供することによって、焼灼精度を向上させ、及び／又は焼灼時間を短縮する。いくつかの例では、システム１００はまた、生検針の存在により、生検標的の有効剛性がその間に増加することが多い画像誘導生検処置と連携して動作する。システム１００はまた、疾患、例えば脂肪肝疾患の結果として起こりうる病変など、強いせん断波減衰を引き起こす組織における組織弾性を決定するために使用される。このような組織は、焼灼された領域又は針を含まない可能性がある。図１に示される１つ又は複数の構成要素は、いくつかの例では、組み合わせることができる。例えば、いくつかの実施形態では、定性的処理モジュール１４０、定量的処理モジュール１６２、及び／又は境界モジュール１５８を組み合わせて、１つの処理モジュールとすることもできる。定量的表示プロセッサ１６４及び定性的表示プロセッサ１４８もまた、いくつかの例では、統合されてもよい。
【００２６】
[049]  参照マップ１４２は、さまざまなソースから導出する。いくつかの例では、参照マップ１４６は、特定の組織タイプを模倣するようにデザインされ、また、組織内の病変の存在を模倣するために剛性が増加した１つ又は複数の領域を含む人工的に作成された組織ファントムから導出される。例えば、参照マップ１４２は、１つ又は複数の腫瘍及び／又は焼灼された領域を含む肝臓ファントム、膀胱ファントム、肺ファントムなどについて作成される。さらに又はあるいは、参照マップ１４２は、実際の患者組織から導出される。参照マップ１４６を作成するために使用される組織は、健康な組織、例えば、病変、腫瘍、又は他の異常のない弾性的に均質な組織である。このような組織から得られたＳＮＲ比は、異種組織よりも高いので、ベースライン参照として組織の信頼性が高まる。しかし、いくつかの例では、少なくとも１つの異常がある組織が参照として利用される。例えば、いくつかの実施形態では、参照マップ１４２は、より早い時点で超音波ＳＷＩを受けた患者から導出された癌性組織のマップである。異なる時点で得られた同じ組織の振幅減衰値を比較することによって、異常組織の拡大又は縮小、したがって治療の進捗状況を決定できる。進捗状況は、治療全体で、又は個別の治療後に追跡できる。例えば、参照マップは、焼灼処置の前に作成され、メモリ１４４に書き込まれ、次いで、比較目的で、処置中又は処置の直後に１回又は複数回メモリから読み出される。参照マップ１４２は、新しいデータが蓄積されるにつれて時間とともに変更されても、又は新しい超音波スキャンが実行されるたびに単に置き換えられてもよい。さらに又はあるいは、参照マップ１４２は、数値モデル、例えば、シミュレートされ、数値計算されたモデルから導出される。モデル特性を変更するために、定量的ＳＷＩ測定値を使用してもしなくてもよい。いくつかの例では、参照組織情報は、組織サンプルや数値モデルから導出されるのではなく、代わりに、特定の媒体の先験的知識から導出される。いくつかの実装形態では、参照マップ１４２は、特定の組織タイプ、例えば、肝臓組織について決定された振幅減衰を、同じ組織タイプから導出された参照マップから決定された振幅減衰値と比較できるように、さまざまな組織について作成される。参照マップ１４２は、システム１００を使用してＳＷＩが実行されるたびに作成して、メモリ１４４に記憶され、これにより、参照マップ１４６のライブラリを時間の経過とともに追補できる。包括的なライブラリは、複数の患者や組織タイプに固有の参照マップを含み得る。
【００２７】
[050]  定性的処理モジュール１４０は、１つ又は複数の要因に基づいて特定の参照マップ１４２を選択的に抽出する。例えば、定性的処理モジュールは、現在検査されている及び／又は焼灼の標的とされている組織タイプに対応する参照マップを選択する。処理モジュールはまた、現在検査されている特定の患者に対応する参照マップを選択する。同じ患者から導出された参照マップを、時間の経過とともに記憶し、処理モジュールは、特定の時点で記憶された参照マップを選択する。いくつかの実施形態では、ユーザーが、特定の参照マップを手動で選択する。さらに又はあるいは、定性的処理モジュールは、ユーザー入力なしで、特定の参照マップを自動的に選択する。このような例によれば、定性的処理モジュールは、参照マップ選択基準を適用する。例えば、定性的処理モジュールは、組織タイプよりも患者のアイデンティティを優先し、これにより、モジュールは、現在検査されている特定の患者に対応する参照マップを最初に要求する。その特定の患者に利用可能な参照マップがない場合、モジュールは、患者のアイデンティティに関係なく、検査されている特定の組織タイプに対応する参照マップの要求に進む。患者の年齢及び／又は健康状態などの追加の基準もまた、処理モジュールによって適用されて、超音波スキャン及び／又は組織焼灼処置中にメモリから特定のマップを選出する。
【００２８】
[051]  図２は、焼灼システムと同時に動作する超音波システムの写真である。超音波トランスデューサ２０２及び焼灼針２０４は、肝臓２０６内の関心領域に向けられている。超音波トランスデューサ２０２は、焼灼針２０４によって標的とされる焼灼部位に隣接する肝臓２０６の領域にプッシュパルスを放出する。示されているように、撮像及び焼灼が同じ組織領域で同時に実行できるように、撮像構成要素及び焼灼構成要素を物理的に近接して結合する。図２に示される組織は、肝臓組織であるが、当業者は、本明細書に説明されるシステム及び方法が肝臓に限定されず、広範囲の組織タイプに適用できることを理解すべきである。
【００２９】
[052]  図３は、図２に示される超音波ＳＷＩシステムによって得られた組織領域の超音波画像である。関心のせん断波領域３０２が、焼灼ゾーン境界３０４及び焼灼針標的部位３０６と共に示されている。せん断波ＲＯＩの弾性マップ３０８が画像にオーバーレイされている。せん断波ＲＯＩ３０２は、各プッシュパルスが組織内へと（この画像では上から下に）放出された後に、せん断波が伝播する組織の領域を包含し、焼灼針標的部位３０６に隣接して配置される焼灼ゾーン境界３０４と重なる。伝播するせん断波の変位振幅を検出するために、横方向に間隔を置いた追跡ビームをせん断波ＲＯＩ３０２内へと放出する。示されているように、プッシュビームは、予想される病変領域のすぐ外側の組織内へと伝達する。したがって、システムは、隣接する、例えば、焼灼された領域とは異なる弾性特性を有する領域の組織弾性特性を決定する。
【００３０】
[053]  さまざまな実施形態では、図３に示されるゾーンのうちの１つ又は複数は、単独で又は焼灼システムと連携して、超音波ＳＷＩシステムによって自動的に定義される。いくつかの例では、システムを操作するユーザーが、せん断波ＲＯＩ３０２、焼灼ゾーン境界３０４、及び／又は焼灼針３０６の標的部位を修正、例えば、再配置及び／又はサイズ変更できる。このような変更は、超音波画像を表示するユーザーインターフェースで入力する。いくつかの実施形態では、あるゾーンの調整が別のゾーンの自動調整を促すように、２つ以上のゾーン間の空間的関係を事前定義する。例えば、焼灼ゾーン境界３０４の調整は、せん断波ＲＯＩ３０２の調整を引き起こすことができ、逆もまた同様である。さらに又はあるいは、例えば、画面上で針を表すグラフィックを動かすことによって、ユーザーインターフェース、例えば、ユーザーインターフェース１５０で焼灼針標的部位３０６を指定する。せん断波ＲＯＩ３０２及び焼灼ゾーン境界３０４は、針標的部位３０６の移動に応答して自動的に動いて、焼灼ゾーンと撮像領域との一貫した空間配置を維持する。超音波トランスデューサによって撮像されるカバレッジエリアは、約１から約１０ｃｍ、約２から約８ｃｍ、約３から約６ｃｍ、約４から約５ｃｍ、又は約２から約３ｃｍの幅の範囲で変化する。さまざまな実施形態では、本明細書のシステムによって生成された組織マップは、カバレッジエリア全体に及ぶ。せん断波がせん断波ＲＯＩ３０２を通って伝播するにつれて、波の振幅は減少する可能性がある。焼灼ゾーン境界３０４内の焼灼された硬い領域の存在など、基礎となる組織特性の変動は、波の振幅に追加の変化を引き起こし、波の減衰パターンを変化させる可能性がある。熱焼灼によって引き起こされる剛性の増加など、組織の基礎となる材料特性の変化は、せん断波振幅の減衰の変化を引き起こす。本明細書のシステムは、せん断波変位振幅を参照変位振幅と比較することによって、このような変動を検出する。例えば、標的領域内の横方向に間隔を置いた点での現在のせん断波振幅と参照せん断波振幅とのピーク間差を決定することによって、減衰変動を定性的処理モジュール１４０によって検出できる。
【００３１】
[054]  図４Ａ及び図４Ｂは、ほぼ一貫した弾性によって特徴付けられる、実際の又は人工の病変のない媒体を含む均質組織媒体内のＲＯＩを通るせん断波伝播の参照振幅変位マップを示す。変位マップは、図１に示されるデータ取得及び表示構成要素と協調して動作する定性的処理モジュール１４０によって生成される。参照マップの生成後、これらはメモリ構成要素１４４に記憶され、後で超音波ＳＷＩ処置の開始に応答して、定性的処理モジュールによって抽出される。図４Ａは、プッシュパルスの送信から９ミリ秒後のせん断波伝播のスナップショットを提供し、図４Ｂは、プッシュパルスの送信から１２ミリ秒のせん断波伝播のスナップショットを提供する。ｙ軸に組織の深さを示し、ｘ軸に組織内の横方向距離を示す。各マップ上に、３つの横方向に間隔を置いた点４０２、４０４、４０６が示されている。各点の座標は、例えば、処理モジュールによって自動的に指定される。いくつかの例では、座標は、例えば、ユーザー入力の受信を介して手動で指定されるが、いくつかの例では、ＲＯＩのみがユーザーによって指定される。各図に示されている点の特定の配置及び数は、説明のみを目的としており、限定的なものと見なされるべきではない。点は、再構成マップに最終的に具体化される弾性情報を決定することに従って指定される。図４Ｃは、経時的に測定された、点４０２、４０４、４０６のそれぞれで検出されたせん断波変位プロファイルのグラフであり、図４Ｄは、横方向に間隔を置いた点の連続する対間の対応する相互相関値のグラフである。図４Ｃにおいて明らかなように、プッシュパルスの送信から約３ミリ秒後に、プッシュパルスの部位の最も近くに位置付けられた第１の点４０２において、最大絶対せん断波変位が検出された。
【００３２】
[055]  図５Ａ及び図５Ｂは、少なくとも１つの病変を含んでいる可能性があり、媒体の残りの部分と比較して異なる弾性値によって特徴付けられる不均質媒体内のＲＯＩを通るせん断波伝播の振幅変位マップを示す。振幅変位マップは、図１の構成要素を使用する超音波ＳＷＩ及び焼灼処置中に、定性的処理モジュール１４０によって生成される。本明細書で説明されるように、変位振幅マップを参照振幅変位マップと比較して、マップ間のせん断波減衰の変動を特定する。図５Ａは、プッシュパルスの送信から９ミリ秒後のせん断波伝搬のスナップショットを提供し、図５Ｂは、プッシュパルス送信から１２ミリ秒後のせん断波伝搬のスナップショットを提供する。各マップ上に、３つの横方向に間隔を置いた点５０２、５０４、５０６が示されている。図４Ｂに示される点のように、図５Ｂに示される点も、自動で、又は、いくつかの例では手動で指定する。上記のように、図５Ｂに示される点は、説明のみを目的としており、限定的なものと見なされるべきではない。図５Ｃは、経時的に測定された、点５０２、５０４、５０６のそれぞれで検出されたせん断波変位プロファイルのグラフであり、図５Ｄは、横方向に間隔を置いた点の連続する対の間の対応する相互相関値のグラフである。
【００３３】
[056]  図４Ｃ及び図５Ｃに示される変位プロファイルから、各点５０２、５０４、５０６における最大振幅変位（「ピーク変位」）を決定するか、又は各点５０２、５０４、５０６における最大値と最小値との差（「ピーク間変位」）を決定して、このような値を、各参照点４０２、４０４、４０６において決定された対応するピーク変位又はピーク間変位と比較することによって、定性的なせん断波振幅変位再構成マップを計算する。例えば、ｘ０，ｙ０にある単一ピクセルの変位振幅値は、ｐ（ｘ０，ｙ０）＝Ｄ（ｘ０，ｙ０）／Ｄｒｅｆ（ｘ０，ｙ０）として与えられる。ここで、Ｄは、変位ピーク振幅又はピーク間振幅を表し、ｒｅｆは、選択した参照マップから抽出された対応する振幅値を示す。せん断波変位振幅は、受信した超音波ビームごとに検出される。いくつかの実施形態では、受信ビームの数は、約４０から約４８の範囲の値をとる。各受信ビームに沿って、変位振幅は、垂直軸に沿った高解像度で、例えば、約３００から約６００ピクセルのそれぞれで決定される。垂直軸の点数は、マップを表示する前に、後のスキャン変換ステップにおいて減少する。定性的マップ１５４は、いくつかの例では、データ取得ユニット１１０で受信した空間的に異なるＡラインのそれぞれについてせん断波振幅のピーク又はピーク間分析を実行することによって得られる。
【００３４】
[057]  いくつかの例では、本明細書のシステムはまた、横方向に間隔を置いた２つの点でのせん断波振幅変位プロファイル間の相関を具体化するせん断波非相関化マップを生成する。横方向に間隔を置いた点間の相関値を、参照マップ内の同等の横方向に間隔を置いた点から導出された相関値と比較することによって、せん断波相関の変動、したがって、組織弾性を特定する。いくつかの実施形態では、非相関化マップは、波形及び周波数の変化にも敏感であることから、せん断波振幅変位マップよりも好ましい場合がある。図４Ｄの線４０８は、第１の点４０２と第２の点４０４との間の相互相関を表し、線４１０は、第２の点４０６及び第３の点４０８との間の経時的な相互相関を表す。同様に、図５Ｄにおいて、線５０８は、図５Ａ及び図５Ｂに示される第１の点５０２と第２の点５０４との間の相互相関を表し、線５１０は、第２の点５０６と第３の点５０８との間の経時的な相互相関を表す。図４Ｄに示されるように、相互相関値は、均質媒体内の各横方向の対の間では、比較的高い。対照的に、図５の線５１０では、相互相関値ははるかに低く、不均質媒体中の点５０４と点５０６との間に病変が存在する可能性が高いことを示している。振幅変位及び相互相関値を対応する参照値と比較することによって、本明細書のシステムは、図６Ａ及び図６Ｂに示される組織弾性の定性的表現を再構成する。これは、せん断波の飛行時間に基づく定量的マップとは異なり、高いせん断波変位ＳＮＲを必要としない。さまざまな実施形態において、参照マップに対するものも含む、定性的処理モジュール１４０によって行われたすべての振幅変位決定は、最大変位量に対して正規化され、したがって、変位振幅の差に基づく相対的な組織弾性変化の正確な描写が提供される。
【００３５】
[058][059]  特定の態様では、参照マップは、評価されている媒体に対応する参照相関値と比較された、ピクセルについて測定された相関値を使用して生成する。例えば、隣接する空間位置でのせん断波の相関係数は、空間剛性分布のマッピングに使用される。このようなマップは、封入体がない場合に最大相関が観察され、また、封入体の検出は相関の喪失（したがって、「非相関化」の呼称の所以である）に基づいているため、非相関化マップと呼ぶ。せん断変位振幅に基づく「従来の」マップは、波の振幅の変化にのみ敏感であるが、２つの空間点の変位プロファイル間の相関係数の変化（非相関化）も波形及び周波数の変化に敏感であり、したがって、材料特性の変化に対する高い感度を実現できる。この高い感度は、低変位振幅および低ＳＮＲの場合の病変の検出及びモニタリングに不可欠である。また、参照マップが実験的に得られ、ピクセルについて測定された相関値が参照相関値で除算されて、均質媒体と比較した組織剛性の違いによる相関の喪失が得られる。（ｘ０，ｙ０）にあるピクセルの非相関化値は、次のように与えられる。
【数１】
2021522004.tif 000003
ここで、＊は、２つの信号間の相互相関演算子を示し、Ｄ（ｘ，ｙ）は、２つの連続する点について測定された、座標（ｘ，ｙ）における時間の関数としての変位信号を示し、下付きのｒｅｆは、参照測定を示す。計算された非相関化マップは、次に、視覚化を良くするためにデシベルスケールで表示される。
【００３６】
[060][061][062][063]  また、比較は、変更された相互相関関数＊を用いて実行できる。
【数２】
2021522004.tif 000004
これは、積分範囲を、適切に選択された最大時間差Ｔに制限し、したがって、信号の潜在的にノイズの多い部分が非常に早い時間／遅い時間に除外される。
【００３７】
[064]  あるいは、制限された積分範囲は、非対称である。つまり、ＴｍｉｎからＴｍａｘの範囲であって、積分時間の最小値及び最大値の値は適切に選択される。
【００３８】
[065][066][067]  あるいは、比較は、式（１）の相互相関演算子を適用する前に、適切に選択された重み付け関数ｗを関数Ｄ、Ｄｒｅｆに掛けることによって、信号Ｄ及び／又はＤｒｅｆの重み付けされたバージョンを用いて実行する。ピークから離れた低信号レベルからのノイズを低減するための適切な重み付け関数の１つは、ガウス関数である。
【数３】
2021522004.tif 000005
スケールパラメータとして任意のａ及び定数ｂは、重み付け関数のピークが期待信号Ｄ（ｘ）のピークと一致するように選択される。
【００３９】
[068]  図６Ａは、図４Ｃ及び図５Ｃから抽出されたデータを比較することによって生成されたせん断波変位振幅再構成マップである。マップは、参照に関連して、振幅変動が高い領域６０２ａを含み、振幅変動が比較的低い領域６０４ａに囲まれた、焼灼によって誘発された病変の存在の可能性を示している。図６Ｂは、図４Ｄ及び図５Ｄから抽出されたデータを比較することによって生成されたせん断波非相関化再構成マップである。焼灼によって誘発された病変の存在及びおおよその位置は、図６Ａに示される振幅変動の高い領域及び低い領域に密接に位置合わせされる、非相関化が比較的低い領域６０４ｂに囲まれた非相関化が高い局所領域６０２ｂによって確証される。図６Ａのせん断波振幅変位再構成マップ及び／又は図６Ｂのせん断波非相関化再構成マップを生成するとき、本明細書の１つ又は複数の構成要素、例えば、定性的処理モジュール１４０及び／又は表示プロセッサ１４８は、横方向に間隔を置いた受信ビームのそれぞれの間に配置された振幅変位値／非相関化値を補間でき、これにより、カバレッジエリア全体の滑らかで空間的に連続なマップが生成される。
【００４０】
[069]  定性的再構成マップにおいて特定された１つ又は複数の病変の外側境界を明確に解明するために、本明細書のシステム構成要素、例えば、境界モジュール１５８は、例えば、図６Ａ及び図６Ｂのマップなどの定性的マップを、同じ領域の少なくとも１つの定量的弾性マップと組み合わせる。そうすることで、定量的弾性マップからノイズを取り除くことができる。したがって、本明細書のシステムは、定量的弾性マップに通常見られるノイズの多い低ＳＮＲデータを省き、代わりに明確な境界によって描写される定義された病変位置を具体化する向上された弾性データを生成する。このようにして、本明細書のシステムは、ノイズを除去し、病変同定の精度を向上させ、処理時間を短縮することによって、超音波ＳＷＩ技術を向上させる。
【００４１】
[070]  図７は、本開示の原理に従って実行されるせん断波撮像方法のフロー図である。例示的な方法７００は、本明細書に説明されるシステム及び／又は装置によって、任意の順序で利用するステップを示している。方法７００は、システム１００などの超音波撮像システム、又は例えば、コーニンクレッカ・フィリップス社（「フィリップス」）によるＬＵＭＩＦＹなどの移動システムを含む他のシステムによって実行できる。追加の例示的なシステムには、同じくフィリップスによって製造されているＳＰＡＲＱ及び／又はＥＰＩＱが含まれる。
【００４２】
[071]  図示の実施形態では、方法７００は、「標的組織に向けて送信された超音波パルスに応答する超音波エコーを取得する」ことによって、ステップ７０２で開始する。
【００４３】
[072]  ステップ７０４において、方法は、「プッシュパルスを標的組織に送信して、標的組織内にせん断波を生成する」ことを含む。
【００４４】
[073]  ステップ７０６において、方法は、「１つ又は複数の場所で、せん断波と交差するように空間的に計画された追跡パルスを送信する」ことを含む。
【００４５】
[074]  ステップ７０８において、方法は、「追跡パルスがせん断波と交差した場所からエコー信号を受信する」ことを含む。
【００４６】
[075]  ステップ７１０において、方法は、「受信したエコー信号から生成された追跡エコーデータを記憶する」ことを含む。
【００４７】
[076]  ステップ７１２において、方法は、「追跡エコーデータに基づいて、標的組織を通って伝播するせん断波の変位振幅を決定する」ことを含む。
【００４８】
[077]  ステップ７１４において、方法は、「決定された変位振幅に基づいて、標的組織の定性的組織弾性マップを生成する」ことを含む。
【００４９】
[078]  図８Ａは、せん断波非相関化再構成マップの例である。図６Ｂに示される非相関化マップのように、図８Ａのマップは、定性的であり、非相関化が高い１つの局所領域８０２ａを含む。境界モジュール１５８は、同様の非相関化の個別の領域を非相関化再構成マップから識別しグループ化する。境界モジュール１５８は、各個別の領域間の境界を定め、これらを描写する等高線を生成し、これにより、図８Ｂに示されるように、非相関化等高線プロットが生成される。
【００５０】
[079][080][081][082]  実際の病変境界に関して等高線の精度を高めるために、境界モジュール１５８はまた、図８Ｃのオーバーレイプロットに示されるように、等高線プロットを定量的組織弾性マップ上にオーバーレイできる。定量的弾性マップは、例えば飛行時間型アプローチに従って、文献に説明されているさまざまなやり方で生成される。図８Ｃに示されているように、非相関化マップから生成された等高線は、定量的マップにおいて明らかである明確な剛性レベルの領域と正確に位置合わせしていない場合がある。等高線が下にある組織内に存在する任意の焦点病変を正確に包むように、等高線プロットと定量的弾性マップとの間の適合を評価及び向上させるために、境界モジュール１５８は２つのメトリックを生成する。最初のメトリックである局限性（focality）（ｆ）は、等高線の重心に関する所与の等高線内の平均剛性分布を具体化する。２番目のメトリックであるペナルティ（ｐ）は、所与の等高線の外側の組織剛性の平均分布を具体化する。オーバーレイプロット上の点（ｘ，ｚ）が等高線内にある場合、その点でのヤング率は、境界モジュール１５８によって等高線の局限性に寄与する。その点が等高線の外側にある場合、その点でのヤング率は、その等高線のペナルティに寄与する。この関係は、数学的に表すことができる。
（ｘ，ｙ）∈Ｃｉの場合、
【数４】
2021522004.tif 000006
であり、
そうではない場合、
【数５】
2021522004.tif 000007
であり、Ｅ（ｘ，ｚ）は、等高線によって画定される領域の重心から距離ｒにある特定のピクセル位置（ｘ，ｚ）における剛性値を示し、Ｃｉは、非相関化ｄＢレベル「ｉ」において等高線によって囲まれた領域を示す。境界モジュール１５８は、ｃｏｓｔ＝ｐ−ｆと定義される費用関数を操作する。各等高線と定量的マップに含まれる局限性病変との適合を向上させ、これにより、いくつかの例において、病変を包含する最適な等高線を生成するために、境界モジュール１５８はさらに、１つ又は複数の候補等高線の局限性値を最大化し、ペナルティ値を最小化することによって、費用関数を最小化する。
【００５１】
[083]  図９Ａは、図９Ｂの定性的非相関化マップから生成された、さまざまな等高線候補値（デシベル（ｄＢ）単位）についてのペナルティ９０２、局限性９０４、及び費用関数９０６のグラフである。示されるように、費用関数９０６は最初に、約２９ｄＢで最小点に到達するまで減少する。次に、費用関数９０６は、ペナルティ９０２項が上昇するにつれて増加し、異なる断面積を包含する各等高線候補の外側にある硬組織を表すピクセルの数が増加していることを示す。２９ｄＢにおける等高線は、図９Ｂの非相関化マップ上に等高線がオーバーレイされたときに、等高線９０８と実際の病変との緊密な適合によって示されるように、図９Ｂに示される病変９１０の最適な等高線を表す。最適な等高線９０８はまた、境界モジュール１５８によって、定量的剛性マップ上にオーバーレイされる。図９Ｃに示されるように、これにより、病変を表さない高剛性の領域を明らかにできる。このような領域は、代わりに、実際の病変に寄与しない信号ノイズ又は高剛性の構造を表す場合がある。ノイズの存在下で実際の病変を明確に描写するために、境界モジュール１５８はさらに、領域９１２などの最適な等高線の外側にある定量的マップ内の領域をマスクし、このようなマスクされた領域がマスクされているハイブリッドマップをディスプレイ上に生成する。図１を参照すると、図示のハイブリッドマップ１６６は、病変１６８及びマスクされた領域１７０を含む。図示の特定の実施形態では、病変１６８のみがマップ１６６に示されている。追加の例では、組織弾性変動の１つ又は複数の領域がマップ上に維持される。例えば、境界モジュール１５８は、単独で、又は境界表示プロセッサ１６４及び／又はユーザーインターフェース１５０と組み合わせて、異なるマスキング閾値を適用して、最終的なハイブリッドマップ１６６に表される詳細を増減する。
【００５２】
[084]  いくつかの例では、境界モジュール１５８は、ペナルティ関数及び局限性関数に最大値演算子を適用して費用関数を定義し、次に費用関数を最小化する等高線を特定する。このような実装形態によれば、費用関数は、図１０Ａに示されるように、ペナルティ関数と局限性関数との交点において最小化される。ペナルティ関数１００２は、約３２．８ｄＢにおいて、局限性関数１００４と交差する。この等高線閾値は、図１０Ｂの定性的非相関化マップ上に表示される第１の最適等高線１００６によって示される。あるいは、最適等高線は、図１０Ｄの勾配マップに示されている定性的非相関化マップの勾配に、少なくとも部分的に基づいて決定される。勾配に基づいて最適等高線を特定するために、さまざまな等高線を勾配マップに描画し、各線について、勾配値を合計し、合計を等高線の長さで割って、各線の平均勾配を算出する。最大平均勾配を有する等高線を選択して、図１０Ｂ及び図１０Ｃの線１００８によって示されるように、非相関化マップ及び／又は定量的剛性マップ上にオーバーレイする。したがって、境界は、最も顕著な勾配の変化の境界を定める。さまざまな例において、境界モジュール１５８は、線１００６及び１００８などの２つ以上の候補等高線を評価し、病変の境界を定めるために１つを選択する。各線について信頼水準を計算し、信頼水準が最も高い線が選択される。さらに又はあるいは、ユーザーによる選択のために、２つ以上の候補等高線を表示する。ユーザーは、候補等高線の視覚的評価に基づいて、及び／又は事前設定された基準を適用することによって、単一の等高線を選択する。いくつかの例では、基準には、過剰封入の組織についてプログラムされた選好が含まれてよく、これにより、標的組織を焼灼しない可能性を減少させる。対照的に、基準には、過小封入の組織についてプログラムされた選好が含まれてもよく、これにより、健康な組織を焼灼する可能性が最小限に抑えられる。
【００５３】
[085]  いくつかの実施形態では、意味のある定量的剛性マップを得るにはＳＮＲが低すぎる場合がある。このような場合、境界モジュール１５８は、焼灼プロセスの熱的モデルから導出されたヒートマップ又は弾性マップ上のペナルティ関数及び局限性関数を決定する。このような実装形態によれば、定性的非相関化マップは、定量的熱データによって具体化できる。
【００５４】
[086]  図１１は、本開示の原理に従って実行されるせん断波撮像方法のフロー図である。例示的な方法１１００は、本明細書に説明されるシステム及び／又は装置によって、任意の順序で利用できるステップを示している。方法１１００は、システム１００などの超音波撮像システム、又は例えば、コーニンクレッカ・フィリップス社（「フィリップス」）によるＬＵＭＩＦＹなどの移動システムを含む他のシステムによって実行できる。追加の例示的なシステムには、同じくフィリップスによって製造されているＳＰＡＲＱ及び／又はＥＰＩＱが含まれる。
【００５５】
[087]  図示の実施形態では、方法１１００は、「標的組織に向けて送信される超音波パルスに応答する超音波エコーを取得する」ことによって、ステップ１１０２で開始する。
【００５６】
[088]  ステップ１１０４において、方法は、「プッシュパルスを標的組織に送信して、標的組織内にせん断波を生成する」ことを含む。
【００５７】
[089]  ステップ１１０６において、方法は、「１つ又は複数の場所で、せん断波と交差するように空間的に計画された追跡パルスを送信する」ことを含む。
【００５８】
[090]  ステップ１１０８において、方法は、「追跡パルスがせん断波と交差した場所からエコー信号を受信する」ことを含む。
【００５９】
[091]  ステップ１１１０において、方法は、「受信したエコー信号から生成された追跡エコーデータを記憶する」ことを含む。
【００６０】
[092]  ステップ１１１２で、方法は、「受信したエコー信号に基づいて、標的組織の定性的組織弾性マップを生成する」ことを含む。
【００６１】
[093]  ステップ１１１４において、方法は、「受信したエコー信号に基づいて、標的組織の定量的組織弾性マップを生成する」ことを含む。
【００６２】
[094]  ステップ１１１６において、方法は、「定性的組織弾性マップを定量的組織弾性マップと組み合わせることによって、標的組織内の剛性が増加した領域の境界を決定する」ことを含む。
【００６３】
[095]  構成要素、システム、及び／又は方法が、コンピュータベースのシステム又はプログラマブルロジックなどのプログラム可能デバイスを使用して実装されるさまざまな実施形態では、上記のシステム及び方法は、「Ｃ」、「Ｃ＋＋」、「ＦＯＲＴＲＡＮ」、「Ｐａｓｃａｌ」、「ＶＨＤＬ」などのさまざまな既知の又は後に開発されたプログラミング言語のいずれかを使用して実装することができることが理解されるべきである。したがって、コンピュータなどのデバイスに上記のシステム及び／又は方法を実装するように指示することができる情報を含むことができる、磁気コンピュータディスク、光ディスク、電子メモリなどのようなさまざまな記憶媒体を準備することができる。適切なデバイスが記憶媒体に含まれる情報及びプログラムにアクセスできるようになると、記憶媒体は、情報及びプログラムをデバイスに提供することができ、したがって、デバイスが本明細書に説明されるシステムの機能及び／又は方法を実行できるようにする。例えば、ソースファイル、オブジェクトファイル、実行可能ファイルなどの適切な資料を含むコンピュータディスクが、コンピュータに提供された場合、コンピュータはその情報を受け取り、それ自体を適切に構成し、さまざまな機能を実装するように上記の図及びフローチャートに概説されるさまざまなシステムの機能及び方法を実行できる。すなわち、コンピュータは、上記のシステム及び／又は方法の異なる要素に関連する情報のさまざまな部分をディスクから受信し、個々のシステム及び／又は方法を実装し、前述の個々のシステムの機能及び／又は方法を調整できる。
【００６４】
[096]  本開示を考慮して、本明細書に説明されるさまざまな方法及びデバイスは、ハードウェア、ソフトウェア、及びファームウェアで実装することができることに留意されたい。さらに、さまざまな方法及びパラメータは、例としてのみ含まれており、限定的な意味では含まれていない。本開示を考慮して、当業者は、本発明の範囲内に留まりながら、それらの技術及びこれらの技術に影響を与えるために必要な機器を決定する際に本発明の教示を実装することができる。本明細書に説明される１つ又は複数のプロセッサの機能は、より少ない数の又は単一の処理ユニット（例えば、ＣＰＵ）に組み込まれてもよく、また、特定用途向け集積回路（ＡＳＩＣ）又は本明細書に説明される機能を実行するための実行可能命令に応答するようにプログラムされた汎用処理回路を使用して実装することができる。
【００６５】
[097]  本システムは、特に超音波撮像システムを参照して説明されたが、本システムは、１つ又は複数の画像が体系的に得られる他の医用撮像システムに拡張できることも想定されている。したがって、本システムは、腎臓、精巣、乳房、卵巣、子宮、甲状腺、肝臓、肺、筋骨格、脾臓、心臓、動脈及び血管系に関連するがこれらに限定されない画像情報だけでなく、超音波誘導介入に関連する他の撮像アプリケーションを取得及び／又は記録するために使用できる。さらに、本システムはまた、本システムの特徴及び利点を提供できるように、従来の撮像システムと共に使用できる１つ又は複数のプログラムを含む。本開示の特定の追加の利点及び特徴は、本開示を研究する際に当業者に明らかであるか、又は本開示の新規のシステム及び方法を使用する人によって経験され得る。本システム及び方法の別の利点は、従来の医用画像システムを容易にアップグレードして、本システム、デバイス、及び方法の特徴及び利点を組み込む点である。
【００６６】
[098]  当然ながら、本明細書に説明される例、実施形態又はプロセスのいずれか１つは、１つ又は複数の他の例、実施形態及び／又はプロセスと組み合わせることができ、又は本システム、デバイス、及び方法に従って、別個のデバイス又はデバイス部分の間で分離及び／又は実行できることが理解されたい。
【００６７】
[099]  最後に、上記の議論は、本システムを単に例示することを意図しており、添付の特許請求の範囲を特定の実施形態又は実施形態のグループに限定するものとして解釈されるべきではない。したがって、本システムは、例示的な実施形態を参照して特に詳細に説明されたが、以下の特許請求の範囲に記載されている本システムのより広く意図された精神及び範囲から逸脱することなく、当業者によって多数の修正及び代替実施形態が考案されうることも理解されたい。したがって、明細書及び図面は、例示的に見なされるべきであり、添付の特許請求の範囲を限定することを意図するものではない。
【図１】
2021522004.tif 000008
【図２】
2021522004.tif 000009
【図３】
2021522004.tif 000010
【図４Ａ】
2021522004.tif 000011
【図４Ｂ】
2021522004.tif 000012
【図４Ｃ】
2021522004.tif 000013
【図４Ｄ】
2021522004.tif 000014
【図５Ａ】
2021522004.tif 000015
【図５Ｂ】
2021522004.tif 000016
【図５Ｃ】
2021522004.tif 000017
【図５Ｄ】
2021522004.tif 000018
【図６Ａ】
2021522004.tif 000019
【図６Ｂ】
2021522004.tif 000020
【図７】
2021522004.tif 000021
【図８Ａ】
2021522004.tif 000022
【図８Ｂ】
2021522004.tif 000023
【図８Ｃ】
2021522004.tif 000024
【図９Ａ】
2021522004.tif 000025
【図９Ｂ】
2021522004.tif 000026
【図９Ｃ】
2021522004.tif 000027
【図１０Ａ】
2021522004.tif 000028
【図１０Ｂ】
2021522004.tif 000029
【図１０Ｃ】
2021522004.tif 000030
【図１０Ｄ】
2021522004.tif 000031
【図１１】
2021522004.tif 000032
【国際調査報告】
2021522004.tif 000033
2021522004.tif 000034
2021522004.tif 000035
2021522004.tif 000036
2021522004.tif 000037