JP2015506258A - 超音波画像における非定常多重反射を抑制するためのフィルタリングシステム及び方法 - Google Patents

超音波画像における非定常多重反射を抑制するためのフィルタリングシステム及び方法 Download PDF

Info

Publication number
JP2015506258A
JP2015506258A JP2014555825A JP2014555825A JP2015506258A JP 2015506258 A JP2015506258 A JP 2015506258A JP 2014555825 A JP2014555825 A JP 2014555825A JP 2014555825 A JP2014555825 A JP 2014555825A JP 2015506258 A JP2015506258 A JP 2015506258A
Authority
JP
Japan
Prior art keywords
filter
signal
ultrasound image
data
ultrasound
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2014555825A
Other languages
English (en)
Other versions
JP6161129B2 (ja
Inventor
クチェウィッツ ジョン
クチェウィッツ ジョン
ピー.カーラ フランチェスコ
ピー.カーラ フランチェスコ
ピー.ダーリントン グレゴリー
ピー.ダーリントン グレゴリー
ディー.ダンバー リー
ディー.ダンバー リー
ジェイ.カツノースキー ピーター
ジェイ.カツノースキー ピーター
リード ジャスティン
リード ジャスティン
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Mirabilis Medica Inc
Original Assignee
Mirabilis Medica Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Mirabilis Medica Inc filed Critical Mirabilis Medica Inc
Publication of JP2015506258A publication Critical patent/JP2015506258A/ja
Application granted granted Critical
Publication of JP6161129B2 publication Critical patent/JP6161129B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/50Image enhancement or restoration using two or more images, e.g. averaging or subtraction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10132Ultrasound image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20212Image combination
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20212Image combination
    • G06T2207/20221Image fusion; Image merging

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本技術は一般に、超音波画像における多重反射によるアーチファクトを抑制するためのフィルタリングシステム及び方法に向けられる。いくつかの実施形態では、フィルタされた超音波画像を得る方法は、アプリケータを用いて目標組織の第1の超音波画像をとることを含む。アプリケータの第1の位置に対する多重反射によるアーチファクトの超音波経路長が変化するようにアプリケータの少なくとも一部が動かされる。次いで、目標組織の第2の超音波画像がとられる。少なくとも1つのフィルタリング方法を用いて第1の超音波画像と第2の超音波画像とが合成される。フィルタリング方法は、合成された超音波画像における多重反射によるアーチファクトを減弱させる又は除去する。

Description

関連出願への相互参照
本願は、参照によりその全体が本明細書に組み込まれる2012年2月2日に出願された米国特許仮出願第61/594,244号に基づく優先権を主張するものである。
本技術は一般に、超音波画像における非定常多重反射を抑制するためのフィルタリングシステム及び方法に向けられる。
従来の医療用超音波イメージングシステムは、トランスデューサ要素のアレイから高周波音波の短いバーストを送波し、次に該音波が組織を通して伝搬してトランスデューサアレイの方にはね返されて戻った後で、該音波を受信する。通常は一定の1540m/sとみなされる、組織の中での音波の速さの知識により、音波がトランスデューサアレイから反射組織に直接移動し、且つトランスデューサアレイに直接戻ったと仮定して、超音波の送波と受波との間の時間遅延に基づいて組織の画像を形成することができる。しかしながら、音波が複数回反射される場合、時間遅延は超音波反射体の物理的場所にもはや対応せず、構造体がそれらの真の物理的場所以外の場所の超音波画像に現れることがある。これらの多重反射によるアーチファクトは、超音波画像内の組織特徴を不明瞭にすることがあり、超音波画像の品質に悪影響を及ぼすことがある。アーチファクトは、多重反射によるアーチファクトの見かけの位置が、画像化されている組織/物質に対して動いているという点で、「非定常」とすることができる。
例えば、流体を通じて患者に結合される別個のイメージングトランスデューサ及び治療トランスデューサを有する超音波治療装置では、非定常多重反射によるアーチファクトが生じる可能性が高く、この場合、イメージングトランスデューサの位置は患者に対して固定され、治療トランスデューサは或る体積の組織を治療するために治療トランスデューサの焦点を変化させるべく患者に対して動いている。このようなシステムでは、同じ装置内の静止した構成部品に超音波が複数回反射される結果として定常アーチファクトが存在することもある。このようなシステムでは、治療トランスデューサ又は他の表面からの二次反射が組織画像内にアーチファクトの多重反射像を形成して、組織の細部を不明瞭にする可能性がある。したがって、超音波画像から多重反射アーチファクトを除去する必要がある。
本技術の実施形態に従って構成される非定常多重反射によるアーチファクトをフィルタリングするためのシステムの略図である。 本技術の実施形態に従う、第1の位置で超音波画像を取得するイメージングトランスデューサの略図である。 本技術の実施形態に従う、図2Aの超音波取得からの多重反射によるアーチファクトの超音波画像の部分略図である。 本技術の実施形態に従う、第2の位置で超音波画像を取得するイメージングトランスデューサの略図である。 本技術の実施形態に従う、図2Cの超音波取得からの多重反射によるアーチファクトの超音波画像の部分略図である。 本技術の実施形態に従う、多重反射によるアーチファクトがフィルタで除去された後の図2B及び図2Dの超音波画像から合成された超音波画像の部分略図である。 本技術の実施形態に従う、超音波データを取得し非定常多重反射によるアーチファクトを除去するためにこれをフィルタリングする方法を示すブロック図である。 本技術の実施形態に従う、フィルタリング前の超音波画像の図である。 本技術の実施形態に従う、図4Aの超音波画像の取得中の超音波信号の速度分布の図である。 本技術の実施形態に従う、図4Bの速度分布に基づくイメージのフィルタリング後の図4Aの超音波画像の図である。 本技術の実施形態に従う、超音波データを取得し、非定常多重反射によるアーチファクトを除去するためにこれをフィルタリングする別の方法を示すブロック図である。 本技術の実施形態に従う、超音波データを取得し、非定常多重反射によるアーチファクトを除去するためにこれをフィルタリングする別の方法を示すブロック図である。 本技術の実施形態に従う、超音波データを取得し、非定常多重反射によるアーチファクトを除去するためにこれをフィルタリングする別の方法を示すブロック図である。
本技術は一般に、超音波画像における多重反射によるアーチファクトを抑制するためのフィルタリングシステム及び方法に向けられる。いくつかの実施形態では、フィルタされた超音波画像を得る方法は、アプリケータを用いて目標組織の第1の超音波画像をとることを含む。アプリケータの第1の位置に対する多重反射によるアーチファクトの超音波経路長が変化するようにアプリケータの少なくとも一部が動かされる。次いで、目標組織の第2の超音波画像がとられる。少なくとも1つのフィルタリング方法を用いて第1の超音波画像と第2の超音波画像が合成される。フィルタリング方法は、合成された超音波画像における多重反射によるアーチファクトを減弱させる又は除去する。
本技術のいくつかの実施形態の具体的な詳細について図1〜図7を参照して以下説明する。超音波イメージングにしばしば関連する周知の構造体及びシステムを説明する他の詳細は、本技術の種々の実施形態の説明を不必要に不明瞭にすることを避けるために以下の開示には記載されていない。図面に示される詳細、寸法、角度、及び他の特徴の多くは、本技術の特定の実施形態の単なる例証である。したがって、他の実施形態は、本技術の精神又は範囲から逸脱することなく他の詳細、寸法、角度、及び特徴を有することができる。当業者は、したがって、本技術が付加的な要素を伴う他の実施形態を有し得ること、又は本技術が図1〜図7を参照して以下に示され説明される特徴のいくつかを伴わない他の実施形態を有し得ることをそれゆえに理解するであろう。
図1は、本技術の実施形態に従って構成された非定常多重反射によるアーチファクトをフィルタリングするためのシステム100の略図である。システム100は、目標組織の超音波画像をとるように構成されたアプリケータ114を含むことができる。さらに詳細に後述するように、アプリケータ114は、患者境界面キャップ130を有するハウジング140の中に又はその近傍に位置決めされるイメージングトランスデューサ150及び高密度焦点式超音波(High Intensity Focused Ultrasound:HIFU)トランスデューサ112を含むことができる。アプリケータ114は、イメージングトランスデューサ150か又はHIFUトランスデューサ112のうちの少なくとも1つを目標組織に対して動かすように構成された移動機構110をさらに含むことができる。
HIFU源又はトランスデューサ112は、高レベルの超音波エネルギーを放出するように構成することができる。いくつかの実施形態では、システム100は、集合的に又は独立して作動される複数のトランスデューサ112を含む。HIFUトランスデューサ112は、HIFUコントローラ118からHIFU処置に関連する持続時間、周波数、振幅、及び/又は他の因子に関する命令を受信することができる。より詳細には、HIFUコントローラ118は、HIFUトランスデューサ112を駆動するために電気波形を生成することができるHIFU送信器122と通信することができる。
イメージングコントローラ102が、システム100の動作を制御することができる。例えば、イメージングコントローラ102は、HIFUコントローラ118と通信する(例えば、命令を送受信する)ことができる。イメージングコントローラ102は、HIFUトランスデューサ112の出力の態様を制御又は構成できるビームフォーマ104とさらに通信する及びこれに命令を送信することができる。例えば、ビームフォーマ104は、超音波ビームを操作するためにHIFUトランスデューサ112の波形遅延を設定することができる。ビームフォーマ104は、パルス駆動回路及び/又は波形駆動回路を含むことができるイメージング送信器106と通信することができる。例えば、イメージング送信器106は、イメージングトランスデューサ150に所望の波形を駆動することができる。イメージングトランスデューサ150は、波形命令を送信/受信スイッチ108からの電気信号として受信することができる。いくつかの実施形態では、送信/受信スイッチ108は、高振幅の伝送信号を遮断するためのステージを含むことができる。イメージングトランスデューサ150は、送信/受信スイッチ108からの電気信号を、HIFUトランスデューサ112によって用いられ得る音響信号に変換することができる。
前述のように、HIFUトランスデューサ112は、1つ以上の移動機構110によって移動又は調整することができる。例えば、いくつかの実施形態では、移動機構110は、HIFUトランスデューサ112の位置を角度θに沿って調整するためにステアリングウォブルを備える。さらなる実施形態では、移動機構110は、HIFUトランスデューサ112及びイメージングトランスデューサ150の一方又は両方のz軸に沿った高さを修正するために高さ調整部を備えることができる。さらなる実施形態では、HIFUトランスデューサ112、イメージングトランスデューサ150、又は他の構成部品は、他の機構によって若しくは他の方向又は角度にピボット運動する、平行移動する、調整する、又は他の方法で動かすことができる。移動機構110は、HIFUトランスデューサ112及び/又はイメージングトランスデューサ150を所望の位置に調整するために、モーションコントローラ116から移動命令及び/又は位置決め命令を受信することができる。図2A〜図7を参照してさらに詳細に後述するように、いくつかの実施形態では、モーションコントローラ116は、非定常多重反射によるアーチファクトを減弱させるためにフィルタリングプロセスの一部として用いられるべき繰り返しデータを得るべくHIFUトランスデューサ112及び/又はイメージングトランスデューサ150を一連のオフセット位置に動かすのに用いることができる。いくつかの実施形態では、モーションコントローラ116は、HIFUコントローラ118と通信することができ、又はHIFUコントローラ118と同じ装置を備えることができる。
いくつかの実施形態では、システム100は、HIFUトランスデューサ112と組織との間に配置することができるトランスデューサ流体120を含む。患者境界面キャップ130は、トランスデューサ流体120をハウジング140の中に保持する、目標組織の近傍の膜界面を備えることができる。いくつかの実施形態では、ハウジング140は、HIFUトランスデューサ112をさらに取り囲む、又は収容する。HIFUトランスデューサ112、イメージングトランスデューサ150、患者境界面キャップ130、ハウジング140、移動機構110、及び/又は他のコンポーネントが総称して「アプリケータ」114と呼ばれているが、さらなる実施形態では、「アプリケータ」114という用語は、超音波イメージング組立体全体、若しくはその個々の構成部品又はサブ構成部品を指すことがあることが理解されるであろう。例えば、いくつかの実施形態では、「アプリケータ」114という用語は、HIFUトランスデューサ112、又はイメージングトランスデューサ150のうちの1つを指すことがある。さらなる実施形態では、「アプリケータ」114という用語は、イメージングトランスデューサ150、HIFUトランスデューサ112、ハウジング140の反射面部分又は他の特徴、若しくはシステム100の他の部分を指すことがある。さらなる実施形態では、「アプリケータ」は、ユーザ(例えば、医者)又は患者とインターフェースするシステム100のあらゆる態様を含む。
音響信号をハウジング140内のHIFUトランスデューサ112又は他の表面に提供するとともに、イメージングトランスデューサ150は、HIFUトランスデューサ112から音響信号を受信し、該音響信号を、送信/受信スイッチ108に送られる及び下流の処理コンポーネントによって処理されるべき電気信号に変換することができる。例えば、送信/受信スイッチ108は、電気信号をイメージング受信器160に送ることができ、これは、トランスデューサ要素から電気信号を受信するための増幅器、フィルタ、及びデータ変換回路を含むことができる。イメージング受信器160は、受信した信号又はデータを無線周波数(RF)データプロセッサ162に送ることができる。より詳細には、前述のように、受信した信号(又は一続きの時間オフセット信号又は位置オフセット信号)を、フィルタリングのためにRFデータプロセッサ162に送ることができる。RFデータプロセッサ162は、RFデータからの多重反射及び/又は他のアーチファクトをフィルタすることができる。さらに詳細に後述するように、フィルタリングは、イメージをレンダリングする前に非定常多重反射信号が減衰されるあらゆるプロセスを含むことができる。RFデータプロセッサ162は、受信したデータを操作するためのRFチャネルデータ処理を行うことができるビームフォーマ164にデータを送ることができる。データは、次いで、多重反射及び/又は他のアーチファクトを除去又は低減させるべくデジタル信号をフィルタする及び/又は処理するために別のRFデータプロセッサ166に送ることができる。他の実施形態では、システム100は、ビームフォーマ164の前又は後に1つだけのRFデータプロセッサを含んでもよいが、さらなる実施形態では、システム100内のこれらの又は他の場所に2つ以上のRFデータプロセッサが存在してもよい。フィルタされたデータは、次いで、ビームにより生成されたデータから超音波画像を形成するためにイメージフォーマ168に送ることができる。イメージは、イメージデータからの多重反射及び/又は他のアーチファクトをフィルタするように構成されたイメージデータプロセッサ170に送信される。いくつかの実施形態では、イメージング受信器160、RFデータプロセッサ162、166、ビームフォーマ164、イメージフォーマ168、及びイメージデータプロセッサ170のうちの1つ以上が、システム制御全体の一部としてイメージングコントローラ102と通信する、命令を受信する、及び/又はフィードバックを提供する。最後に、フィルタされた画像がイメージディスプレイパネル172などの出力装置又はモニタ上に表示される。システム100は、異なる特徴を含んでもよい及び/又は異なる構成を有してもよいことが理解されるであろう。例えば、さらなる実施形態では、上記のシステム100の構成部品のうちの1つ以上を削除することができる。さらなる実施形態では、システム100の個々に説明された構成部品の機能性は、単一の装置又はサブ構成部品に組み合わされてもよい。
「フィルタリング」という用語は、本明細書で広義に用いられており、あらゆる適切なフィルタリング方法又は技術を含むことができる。信号処理経路内の種々の地点に種々のフィルタリング方法を個々に又は組み合わせて実装できる可能性がある。可能なフィルタは、有限インパルス応答(finite−impulse−response:FIR)フィルタ、無限インパルス(infinite−impulse−response:IIR)フィルタ、回帰フィルタ、及びパーセンタイルフィルタを含むがこれらに限定されない。フィルタリングは、ビームを形成する前に(例えば、データプロセッサ162で)、ビームを形成した後で(例えば、RFデータプロセッサ166で)、又はイメージデータプロセッサ170で画像を形成した後でチャネル毎に行うことができる可能性がある。フィルタの選択は、フィルタリング方法が実施される超音波システムのハードウェア能力及びソフトウェア能力によって決まることがある。いくつかの適切なフィルタに関するさらなる詳細が以下で提供される。
FIR/IIRフィルタ−FIRフィルタ及びIIRフィルタは、特に組織からの信号を抑制するために血流をドップラーイメージングするための超音波システムで用いることができる。いくつかの実施形態では、事前に計算されたフィルタ係数の固定の組と共にローパスフィルタ又は帯域消去フィルタを用いることができる。別の実施形態では、多重反射によるアーチファクトの速度(又は動きを示す他のメトリック)に基づいて、フィルタ係数を適応的に計算することができる。この実施形態では、ドップラー方法(又は動きを推定するための他の公知の超音波方法)が速度を推定するのに用いられることになり、この速度の推定から、フィルタ係数を、多重反射を最適に抑制するように適応的に調整することができる。別の実施形態では、パルスの送波/受波中の多重反射源の速度の事前知識に基づいて、フィルタ係数の複数の組を事前に計算する及び入れ替えることができる可能性がある。
回帰フィルタ−回帰フィルタは、例えば、特にカラードップラー及びパワードップラーに関するドップラー超音波処理に用いることができる。カラードップラー及びパワードップラーでは、超音波データは低次多項式にあてはめることができ、剰余は血流計算に用いることができる。多重反射の抑制のために、剰余は、多重反射信号によって支配されることになり、切り捨てることができ、あてはめは組織信号によって支配されるであろう。
平均化フィルタ−平均化フィルタは、複数のパルスからの信号の平均をとる。これらのフィルタは、FIRフィルタの概念的変形である。平均は、加重又は非加重とすることができる可能性がある。いくつかの実施形態では、各パルスは、最も強い信号が多重反射であるという動作上の仮定の下で、その振幅に逆比例して加重することができる可能性がある。
パーセンタイルフィルタ−パーセンタイルフィルタは、信号を振幅によってソートし、特定の閾値の、特定の閾値を上回る、又は特定の閾値を下回る信号を保持する。一実施形態では、ピクセルごとに最小振幅を有するパルスだけが保持される最小値フィルタが用いられる。
本明細書に記載のシステム100及び他のシステム/方法は、前述のフィルタ又は他の適切なフィルタのうちの1つ以上を使用してもよいことが理解されるであろう。
図2Aは、本技術の実施形態に係る第1の位置で超音波画像を取得するイメージングトランスデューサの略図である。図2Bは、本技術の実施形態に係る図2Aの超音波取得からの多重反射によるアーチファクトの超音波画像の部分略図である。図2Aと図2Bを一緒に参照すると、イメージングトランスデューサ150は、少なくとも部分的にハウジング140の中に深さzに位置決めされる。イメージングトランスデューサ150は、超音波を、経路Lに沿ってトランスデューサ流体120及び患者境界面キャップ130を通して患者組織200に送波する。組織200は、超音波画像エコー又は体積を生じる組織200における特徴である組織特徴202を含む。図2Bに示すように、超音波画像204は、医師などのユーザに表示されるように構成される。超音波画像204は、dによって表される等しい深さで示される多重反射信号経路長L+L+L+Lを有する多重反射された音波波形の結果であるアーチファクト206を含む。
図2C及び図2Dは類似の構成を示している。例えば、図2Cは、本技術の実施形態に従う、第2の位置で超音波画像を取得するイメージングトランスデューサの略図である。この配置では、しかしながら、イメージングトランスデューサ150は深さzに位置決めされる。図2Dは、本技術の実施形態に従う、図2Cの超音波取得からの多重反射によるアーチファクトの超音波画像の部分略図である。図2Dに示すように、超音波画像204は、dによって表される等しい深さに示される経路長L+L+L+Lを有する多重反射された音波波形の結果であるアーチファクト206を含む。
多重反射によるアーチファクト206は、トランスデューサ流体体積120内で音波波形の2回以上の反射が起き、反射した音波波形のすべて又は一部がイメージングトランスデューサ150によって受波されるときに生じる。多重反射によるアーチファクト条件は、セラミックトランスデューサ要素と組織境界面との間の物質内の従来のイメージングトランスデューサで生じるが、エネルギーは、普通は関心あるイメージング深さの前に十分なレベルに減弱する。イメージングトランスデューサ150は、多重反射された波形を組織から反射された波形と同時に受波する。これらの2つの受波された波形が、超音波システムによって処理され、多重反射によるアーチファクト206が多重反射経路長の半分に等しいトランスデューサからの距離((L+L+L+L)/2)に現れるように超音波画像204上に表示される。示される多重反射経路(L+L+L+L)は単純に2つの表面反射のものであるが、音波波形は、イメージングトランスデューサ150によって受波される前にトランスデューサ流体120内の多くの表面から反射される可能性がある。イメージデータが受信され、超音波画像に処理されるときに、イメージングトランスデューサ150の表面は、多重反射経路長(d又はd)の距離に等しいイメージングトランスデューサ150からの距離に多重反射によるアーチファクト206となって現れる。
図1を参照して前述したシステム100は、イメージングトランスデューサ150、又は反射面(HIFUトランスデューサ112など)のいずれかがシフトされて多重反射経路長(L+L+L+L)の変化をもたらすときの多重反射経路長のシフトを利用することができる。例えば、イメージングトランスデューサ150が位置zから位置zにシフトされるときに、結果的に多重反射によるアーチファクト206がそれぞれの深さd及びdに現れる。図2Eは、例えば、多重反射によるアーチファクトがフィルタで除去された後の図2B及び図2Dの超音波画像から合成された超音波画像の部分略図である。図2Eに示すように、それぞれz及びzで取り込まれた2つのデータセットを組み合わせることによって、多重反射によるアーチファクト206が存在しない第3のデータセットをもたらすことができる。本明細書で論じられる技術は2つのデータセットを処理することに言及できるが、それぞれのアルゴリズムを2つ以上のデータセットに同様に適用することができる。いくつかの実施形態では、トランスデューサの位置及びそれぞれのデータセットの数は4〜8の範囲内とすることができる。さらなる実施形態では、セラミックトランスデューサアレイのセラミック要素の異なる組から異なる操作角度で音波波形が送波される複合イメージング技術を用いるときに、多重反射経路長をシフトさせることができる。
本明細書に記載のいくつかの技術に関して、多重反射によるアーチファクト206は、該アーチファクト206が深さdのアーチファクト206と重ならない、深さdに存在することが望ましい場合があるが、アーチファクト206の「速度」を利用する技術などの他の技術は、d−dの分離距離を必要としない。これらの場合、受信されるエネルギーの位置の変化が小さく、エネルギーが時間及び空間において重なることがあるが、速度を依然として利用することができる。次いで、順次の又は一連のデータセット間の「移動(moving)」エネルギーをフィルタで除去するために、現在公知の信号処理技術及びアルゴリズムを適用することができる。
説明されたアルゴリズムは、複数の位置又は送波角度でデータを取り込み、単一のイメージに関する複数のデータセットを処理する、逐次プロセスを示すが、さらなる実施形態では、データは、表示されるイメージの数が取得したデータセットの数に等しい(例えばスライディングウィンドウ及び並列処理を用いる)ような様態で処理できることにも留意されたい。
超音波データを取得し、非定常多重反射によるアーチファクトを減弱させるべくデータをフィルタするために種々の方法を採用することができる。図3にブロック図として示される第1の方法300では、固定の又は可変の時間間隔だけ分離された一連の繰り返される超音波パルスが得られ、移動するHIFUトランスデューサなどの非定常多重反射によるアーチファクトに帰属する信号の部分を減弱させるために画像を形成する前にフィルタされる、パルスの集合体を画定する。ブロック302で、方法300は、時刻t(1)にイメージングトランスデューサを位置決めすること、及びカウンタkを2に設定することを含む。ブロック304で、時刻t(1)に超音波画像が取得される。ブロック306で、方法300は、時刻t(k)まで待つことを含む。ブロック308で後続の超音波画像kが取得される。このプロセスは、プロセスが最終時刻t(K)に到達するまで反復的超音波画像と共にブロック310及び312で繰り返される。前述のように、超音波パルス間の時間は、固定、可変、又はランダムとすることができる。
ブロック316で、カウンタnがイメージデータにおける最初のピクセルを示すポイントに設定される。ブロック318でイメージ1からイメージKまでのピクセルnでのピクセルデータベクトルが生成される。ブロック320で、フィルタ(図1を参照して上述したフィルタのいずれかなど)がベクトルにわたって適用される。フィルタは、組織又は体液よりも速く動いている多重反射によるアーチファクトなどの、所定の範囲を外れる速度で動いている超音波画像からのアーチファクトを除去するのに用いることができる。さらなる実施形態では、フィルタは、動いている組織及び流体が多重反射によるアーチファクトと共にフィルタで除去されるように設定することができ、こうしたフィルタは、例えばイメージの機能が静止した組織境界を伝えることである場合に有用であることがある。特定の実施形態では、フィルタは、回帰フィルタ、FIRフィルタ、又はIIRフィルタなどのローパスフィルタである。代替的に、その見かけの速度が事前に分かっている場合、多重反射によるアーチファクトをより選択的にフィルタするのに帯域消去フィルタを用いることができる可能性がある。
ブロック322で、フィルタリングの結果がイメージ1のピクセルnに格納される。このフィルタリングプロセスは、イメージファイルの終わりに到達するまでブロック324及び326で繰り返される。いくつかの実施形態では、フィルタ性能を改善するために、パルス間の時間間隔、集合体の長さ、及び/又はフィルタの次数及び係数を手動で又は適応的に変えることができる可能性がある。イメージがフィルタされた後で、ブロック328でイメージ1が表示される。他の実施形態では、方法300は、付加的なステップを含んでもよく、及び/又は上記のステップの1つ以上が削除されてもよいことが理解されるであろう。
図4Aは、本技術の実施形態に係るフィルタリング前の超音波画像の図である。図4Bは、本技術の実施形態に係る図4Aの超音波画像の取得中の超音波信号の速度分布の図である。図4Cは、図3を参照して上述した方法に係る図4Bの速度分布に基づくイメージのフィルタリング後の図4Aの超音波画像の図である。特定の実施形態では、16パルスの集合体と共に回帰フィルタを用いることができる。結果として得られるフィルタされた超音波画像は、高速多重反射によるアーチファクトの顕著な低減を示し、目標組織のよりはっきりとした画像を提供する。
図3を参照して説明される方法300の有効性は、組織の速度と多重反射の見かけの速度との相違に依存する。図5は、組織の速度と多重反射の速度との相違をもたらす又は増強するためにイメージングトランスデューサを平行移動することを含む方法500を示すブロック図である。イメージングトランスデューサが平行移動される際に、組織の見かけの位置が等しい量だけ平行移動することになり、より速く動く多重反射の見かけの位置がより大きい量だけ平行移動することになる。この方法500を用いると、定常アーチファクトと非定常アーチファクトとの両方を、減弱又は除去のためにこのようにより容易に識別することができる。
ブロック502で、例えば、方法500は、イメージングトランスデューサを位置z(1)に位置決めすること、及びカウンタkを2に設定することを含む。ブロック504で、位置z(1)での超音波画像が取得される。ブロック506で、方法500は、イメージングトランスデューサを位置z(k)に動かすことを含む。ブロック508で、位置z(k)での後続の超音波画像kが取得される。このプロセスは、プロセスが最終位置z(K)に到達するまで反復的超音波画像と共にブロック510及び512で繰り返される。いくつかの実施形態では、イメージングトランスデューサの位置はコンピュータの制御下にあり、一方、他の実施形態では、イメージングトランスデューサは手動制御の下にある。
ブロック514で、イメージ1とイメージkが位置合わせされ、ブロック516で、カウンタnがイメージデータにおける最初のピクセルを示すポイントに設定される。ブロック518で、イメージ1からイメージKまでのピクセルnでのピクセルデータベクトルが生成される。ブロック520で、フィルタ(図1を参照して上述したフィルタのいずれかなど)がベクトルにわたって適用される。フィルタは、2つ以上のイメージフレームにおける2つ以上の個別の及び別個の場所に入る超音波画像からのアーチファクトを除去するのに用いることができる。代替的に、フィルタは、組織又は体液よりも速く動いている多重反射によるアーチファクトなどの、所定の範囲を外れる速度で動いている超音波画像からのアーチファクトを除去するのに用いることができる。イメージングトランスデューサがデータを複数の位置で取得したために、多重反射の速度が誇張され(より容易にフィルタで除去され)得る。ブロック522で、フィルタリングの結果がイメージ1のピクセルnに格納される。このフィルタリングプロセスは、イメージファイルの終わりに到達するまでブロック524及び526で繰り返される。ブロック528で、フィルタされた画像が表示される。
図6は、超音波データを取得し、非定常多重反射によるアーチファクトを除去するためにフィルタリングする別の方法600を示すブロック図である。多重反射オブジェクトの見かけの位置を考慮に入れる方法600は、音の伝搬方向に依存するものとすることができる。方法600は、異なる方向に送波される一連のパルスを用いる。パルスは、次いで、全イメージング方向に共通でない特徴を抑制するためにフィルタされる。目標組織及び/又はHIFUシステムの部品の移動によって生じるアーチファクトを低減させるために、特定の焦点領域に向けて一連のパルスを行うことができる。
ブロック602で、例えば、方法600は、超音波波形を角度θ(1)で送波するようにビームフォーマをプログラミングすることを含む。ブロック604で波形が送波され、カウンタkが2に設定され、ブロック606で超音波画像1が取得される。反復的に、ブロック608〜616で、最終送波角度θ(K)が取得されるまで後続の超音波画像が反復的に取得される。
ブロック618で、イメージ1とイメージkが位置合わせされ、ブロック620で、カウンタnがイメージデータにおける最初のピクセルを示すポイントに設定される。ブロック622で、イメージ1からイメージKまでのピクセルnでのピクセルデータベクトルが生成される。ブロック624で、フィルタ(最小値フィルタ、加重平均フィルタ、又は他のフィルタなど)がベクトルにわたって適用される。これらの複数の操作された角度の空間的コンパウンディングを、多重反射オブジェクトを減弱させることに加えてスペックルを低減させるのに用いることができる。種々の操作されたパルスへの多重反射オブジェクトの反射は組織の分散よりも高い分散を有するであろう。この分散は、多重反射オブジェクトをセグメントにして除去するのに用いることができる。ブロック626で、フィルタリングの結果がイメージ1のピクセルnに格納される。このフィルタリングプロセスは、イメージファイルの終わりに到達するまでブロック628及び630で繰り返される。ブロック632で、フィルタされた画像が表示される。
図7は、超音波データを取得し、非定常多重反射によるアーチファクトを除去するためにフィルタリングするまた別の方法700を示すブロック図である。図7に示すように、多重反射をさらに抑制するために、一連のイメージにさらなるフィルタリングを適用できる可能性がある。単一のイメージ内で多重反射を抑制できる度合いは、多重反射によるアーチファクトを生成する構造体の速度と共に変化するであろう。抑制の度合いがイメージからイメージにかけて大いに変化する場合、さらなる多重反射の抑制のために一連のイメージをフィルタ又は「ステッチ」することができる可能性がある。可能なフィルタは、イメージ間に適用される平均値フィルタ、メジアンフィルタ、及び最小値フィルタを含むがこれらに限定されない。
ブロック702で、例えば、方法700は、イメージングトランスデューサを位置z(1)に位置決めすること、及びカウンタkを2に設定することを含む。ブロック704で、位置z(1)での超音波画像が取得される。ブロック706で、方法700は、イメージングトランスデューサを位置z(k)に動かすことを含む。ブロック708で、位置z(k)での後続の超音波画像kが取得される。このプロセスは、トランスデューサが最終位置z(K)に到達するまで反復的超音波画像と共にブロック710及び712で繰り返される。いくつかの実施形態では、イメージングトランスデューサの位置はコンピュータの制御下にあり、一方、他の実施形態では、イメージングトランスデューサは手動制御の下にある。
ブロック714で、イメージ1とイメージkが位置合わせされ、ブロック716で、カウンタnがイメージデータにおける最初のピクセルを示すポイントに設定される。ブロック718で、イメージ1からイメージKまでのピクセルnでのピクセルデータベクトルが生成される。ブロック720で、ピクセル輝度の分散がイメージにわたって比較され、フィルタリング技術が適用される。いくつかの実施形態では、ブロック722及び724で示されるように、分散が所定の閾値を上回る場合、最小の(すなわち最も暗い)ピクセル値がイメージ1上に保存される。分散が閾値内の場合、ブロック726でピクセルデータの平均がイメージ1上に保存される。このプロセスは、イメージファイルにおけるすべてのピクセルに関してブロック728及び730で繰り返される。1つの代替的な実施形態では、[p(1,n)、p(2,n)、…、p(K,n)]の最小値だけがイメージ1のピクセルnに格納される。さらに別の実施形態では、[p(1,n)、p(2,n)、…、p(K,n)]の平均値だけがイメージ1のピクセルnに格納される。さらに別の実施形態では、HIFUトランスデューサ100を動かすことによって多重反射経路長が変化する。ブロック732で、フィルタされた画像が表示される。
付記
以下の付記は本技術のいくつかの例証となる。
1.超音波画像を得る方法であって、
アプリケータを介して目標組織の第1の超音波画像をとることと、
前記アプリケータの少なくとも一部を動かすことと、
前記アプリケータを介して目標組織の第2の超音波画像をとることと、
前記第1の超音波画像と前記第2の超音波画像とを合成することと、
前記合成された超音波画像における多重反射によるアーチファクトを除去することと、を含む、方法。
2.前記アプリケータの少なくとも一部を動かすことが、前記アプリケータの反射面を動かすことを含む、付記1に記載の方法。
3.前記アプリケータの少なくとも一部を動かすことが、イメージングトランスデューサを動かすことを含む、付記1に記載の方法。
4.前記アプリケータの少なくとも一部を動かすことが、高密度焦点式超音波トランスデューサを動かすことを含む、付記1に記載の方法。
5.前記アプリケータの少なくとも一部を動かすことが、前記アプリケータの一部を平行移動すること又はその角度を調整することのうちの少なくとも1つを含む、付記1に記載の方法。
6.前記多重反射によるアーチファクトを除去することが、前記第1の超音波画像及び前記第2の超音波画像にローパスフィルタを用いることを含む、付記1に記載の方法。
7.前記多重反射によるアーチファクトを除去することが、回帰フィルタ、有限インパルス応答フィルタ、無限インパルス応答フィルタ、又は帯域消去フィルタのうちの少なくとも1つを用いることを含む、付記1に記載の方法。
8.前記多重反射によるアーチファクトを除去することが、平均値フィルタ、メジアンフィルタ、又は最小値フィルタのうちの少なくとも1つを用いることを含む、付記1に記載の方法。
9.超音波イメージングにおける非定常多重反射によるアーチファクトを抑制する方法であって、
イメージングトランスデューサによって、或る体積の組織に複数の超音波パルスを送波することであって、前記複数の超音波パルスの送波は或る時間間隔だけ分離される、ことと、
前記複数の超音波パルスの各送波に応答して、前記或る体積の組織から反射された信号を受信することであって、前記反射された信号は、前記或る体積の組織から発生したものではない多重反射信号を含む、ことと、
前記多重反射信号に帰属する前記反射された信号の部分を減弱させるために前記反射された信号をフィルタリングすることと、
前記フィルタされた信号に基づいて前記或る体積の組織の超音波画像を生成することと、
を含む、方法。
10.前記複数の超音波パルスを送波することが、固定の又は可変の時間間隔だけ分離されたパルスを送波することを含む、付記9に記載の方法。
11.前記反射された信号をフィルタリングすることが、ローパスフィルタを用いて前記信号をフィルタリングすることを含み、前記ローパスフィルタが、回帰フィルタ、有限インパルス応答フィルタ、無限インパルス応答フィルタ、及び帯域消去フィルタのうちの少なくとも1つである、付記9に記載の方法。
12.前記イメージングトランスデューサを前記複数の超音波パルスの送波間で固定量だけ平行移動することをさらに含む、付記9に記載の方法。
13.前記イメージングトランスデューサによって、前記或る体積の組織に第2の複数の超音波パルスを送波することと、
前記第2の複数の超音波パルスの各々に応答して、前記或る体積の組織から第2の反射された信号を受信することであって、前記第2の反射された信号は、前記或る体積の組織から発生したものではない多重反射信号を含む、ことと、
前記多重反射信号に帰属する前記第2の反射された信号の部分を減弱させるために前記第2の反射された信号をフィルタリングすることと、
前記第2のフィルタされた信号に基づいて前記或る体積の組織の第2の超音波画像を生成することと、
前記多重反射信号に帰属する前記イメージの部分をさらに減弱させるために前記第1の超音波画像と前記第2の超音波画像を比較することと、
をさらに含む、付記9に記載の方法。
14.前記第1の超音波画像と前記第2の超音波画像とを比較することが、平均値フィルタ、メジアンフィルタ、及び最小値フィルタのうちの少なくとも1つを用いて前記イメージをフィルタリングすることを含む、付記13に記載の方法。
15.超音波画像における非定常のアーチファクトを減弱させる方法であって、
第1の信号経路上に第1の信号を印加することによって前記アーチファクトを含む超音波画像データの第1の組を取り込むことと、
第2の信号経路上に第2の信号を印加することによって超音波画像データの第2の組を取り込むことと、
前記アーチファクトを除去するために前記データの第1の組と前記データの第2の組を組み合わせることと、
を含む、方法。
16.前記データの第1の組と前記データの第2の組とを組み合わせることが、ローパスフィルタを用いて前記データの第1の組及び前記データの第2の組をフィルタリングすることを含む、付記15に記載の方法。
17.前記超音波画像データの第1の組の取り込みが第1の時間間隔で行われ、前記超音波画像データの第2の組の取り込みが前記第1の時間間隔から離間した第2の時間間隔で行われる、付記15に記載の方法。
18.前記第2の信号経路上に第2の信号を印加することが、前記第1の信号経路とは異なる経路長又は経路ルートを有する第2の信号経路上に第2の信号を印加することを含む、付記15に記載の方法。
19.前記アーチファクトを除去するために前記データの第1の組と前記データの第2の組を組み合わせることが、前記データの第1の組及び前記データの第2の組に平均値フィルタ、メジアンフィルタ、又は最小値フィルタのうちの少なくとも1つを適用することを含む、付記15に記載の方法。
20.前記データの第1の組と前記データの第2の組との間の分散に少なくとも部分的に基づいて前記平均値フィルタ、前記メジアンフィルタ、又は前記最小値フィルタのうちの少なくとも1つを選択することをさらに含む、付記19に記載の方法。
21.コンピューティングシステムによって実行される場合に、前記コンピューティングシステムに、
第1の信号経路上に第1の信号を印加するようにイメージングトランスデューサに命令することと、
前記第1の信号に基づいてアーチファクトを含む超音波画像データの第1の組を受信することと、
第2の信号経路上に第2の信号を印加するように前記イメージングトランスデューサに命令することと、
前記第2の信号に基づいて超音波画像データの第2の組を受信することと、
前記アーチファクトを除去するために前記データの第1の組と前記データの第2の組とを組み合わせることと、
を含む動作を行わせるコンピュータで実行可能な命令を格納している物理的コンピュータ可読記憶媒体。
22.前記動作が、平均値フィルタ、メジアンフィルタ、又は最小値フィルタのうちの少なくとも1つを用いて前記超音波画像データの第1の組及び前記超音波画像データの第2の組をフィルタリングすることをさらに含む、付記21に記載の物理的コンピュータ可読記憶媒体。
23.前記動作が、ローパスフィルタを用いて前記データの第1の組及び前記データの第2の組をフィルタリングすることをさらに含む、付記21に記載の物理的コンピュータ可読記憶媒体。
24.前記第2の信号経路上に第2の信号を印加するようにイメージングトランスデューサに命令することが、前記第1の信号経路とは異なる経路長又は経路ルートを有する第2の信号経路上に前記第2の信号を印加するように前記イメージングトランスデューサに命令することを含む、付記21に記載の物理的コンピュータ可読記憶媒体。
25.前記動作が、前記イメージングトランスデューサを前記第1の信号の印加と前記第2の信号の印加との間で固定量だけ平行移動することをさらに含む、付記21に記載の物理的コンピュータ可読記憶媒体。
26.超音波イメージングシステムであって、
或る体積の組織に超音波パルスを送達するように構成されたイメージングトランスデューサと、
前記或る体積の組織から反射された信号を受信するように構成された受信器と、この場合、前記反射された信号は、前記或る体積の組織から発生したものではない多重反射信号を含み、
前記多重反射信号に帰属する前記反射された信号の部分を減弱させるために前記反射された信号をフィルタし、前記フィルタされた信号に基づいて前記或る体積の組織の超音波画像を生成するようにプログラムされたコントローラと、
を備える、システム。
27.前記イメージングトランスデューサが、固定の又は可変の時間間隔、分離された超音波パルスを送達するように構成される、付記26に記載のシステム。
28.前記反射された信号をフィルタするようにプログラムされたコントローラが、ローパスフィルタを用いて前記反射された信号をフィルタするようにプログラムされたコントローラを含み、前記ローパスフィルタが、回帰フィルタ、有限インパルス応答フィルタ、無限インパルス応答フィルタ、及び帯域消去フィルタのうちの少なくとも1つを含む、付記26に記載のシステム。
29.前記イメージングトランスデューサを動かし、それにより前記反射された信号の前記長さを修正するように構成された移動機構をさらに備える、付記26に記載のシステム。
30.前記反射された信号をフィルタするようにプログラムされたコントローラが、平均値フィルタ、メジアンフィルタ、又は最小値フィルタのうちの少なくとも1つを用いて前記反射された信号をフィルタするようにプログラムされたコントローラを含む、付記26に記載のシステム。
上記のことから、本技術の具体的な実施形態は、例示する目的で本明細書で説明されているが、本技術の精神及び範囲から逸脱することなく種々の修正がなされ得ることが理解されるであろう。例えば、本明細書に記載のシステム及び方法は、非定常多重反射によるアーチファクト以外のアーチファクトをフィルタするのに用いることができる。より詳細には、定常アーチファクトを除去する手段として、非定常多重反射によるアーチファクトの技術はまた、患者組織境界面に対するイメージングトランスデューサの距離を変化させることによって、他の方法では定常アーチファクトのみを有するであろう装置に用いることができる。さらに、特定の実施形態との関連で説明される新しい技術の或る態様は、他の実施形態において組み合わされてもよく、又は削除されてもよい。さらに、本技術の特定の実施形態に関連する利点がそれらの実施形態との関連で説明されているが、他の実施形態がこうした利点を呈してもよく、すべての実施形態が必ずしも本技術の範囲内に入るようにこうした利点を呈する必要はない。したがって、本開示及び関連する技術は、本明細書に明確に示されていない又は記載されていない他の実施形態を包含することができる。したがって、本開示は、付属の請求項以外によって限定されない。

Claims (30)

  1. 超音波画像を得る方法であって、
    アプリケータを介して目標組織の第1の超音波画像をとることと、
    前記アプリケータの少なくとも一部を動かすことと、
    前記アプリケータを介して前記目標組織の第2の超音波画像をとることと、
    前記第1の超音波画像と前記第2の超音波画像とを合成することと、
    前記合成された超音波画像における多重反射によるアーチファクトを除去することと、を含む、方法。
  2. 前記アプリケータの少なくとも一部を動かすことが、前記アプリケータの反射面を動かすことを含む、請求項1に記載の方法。
  3. 前記アプリケータの少なくとも一部を動かすことが、イメージングトランスデューサを動かすことを含む、請求項1に記載の方法。
  4. 前記アプリケータの少なくとも一部を動かすことが、高密度焦点式超音波トランスデューサを動かすことを含む、請求項1に記載の方法。
  5. 前記アプリケータの少なくとも一部を動かすことが、前記アプリケータの一部を平行移動すること又はその角度を調整することのうちの少なくとも1つを含む、請求項1に記載の方法。
  6. 前記多重反射によるアーチファクトを除去することが、前記第1の超音波画像及び前記第2の超音波画像にローパスフィルタを用いることを含む、請求項1に記載の方法。
  7. 前記多重反射によるアーチファクトを除去することが、回帰フィルタ、有限インパルス応答フィルタ、無限インパルス応答フィルタ、又は帯域消去フィルタのうちの少なくとも1つを用いることを含む、請求項1に記載の方法。
  8. 前記多重反射によるアーチファクトを除去することが、平均値フィルタ、メジアンフィルタ、又は最小値フィルタのうちの少なくとも1つを用いることを含む、請求項1に記載の方法。
  9. 超音波イメージングにおける非定常多重反射によるアーチファクトを抑制する方法であって、
    イメージングトランスデューサによって、或る体積の組織に複数の超音波パルスを送波することであって、前記複数の超音波パルスの送波は或る時間間隔だけ分離される、ことと、
    前記複数の超音波パルスの各送波に応答して、前記或る体積の組織から反射された信号を受信することであって、前記反射された信号は、前記或る体積の組織から発生したものではない多重反射信号を含む、ことと、
    前記多重反射信号に帰属する前記反射された信号の部分を減弱させるために前記反射された信号をフィルタリングすることと、
    前記フィルタされた信号に基づいて前記或る体積の組織の超音波画像を生成することと、
    を含む、方法。
  10. 前記複数の超音波パルスを送波することが、固定の又は可変の時間間隔だけ分離されたパルスを送波することを含む、請求項9に記載の方法。
  11. 前記反射された信号をフィルタリングすることが、ローパスフィルタを用いて前記信号をフィルタリングすることを含み、前記ローパスフィルタが、回帰フィルタ、有限インパルス応答フィルタ、無限インパルス応答フィルタ、及び帯域消去フィルタのうちの少なくとも1つである、請求項9に記載の方法。
  12. 前記イメージングトランスデューサを前記複数の超音波パルスの前記送波間で固定量だけ平行移動することをさらに含む、請求項9に記載の方法。
  13. 前記イメージングトランスデューサによって、前記或る体積の組織に第2の複数の超音波パルスを送波することと、
    前記第2の複数の超音波パルスの各々に応答して、前記或る体積の組織から第2の反射された信号を受信することであって、前記第2の反射された信号は、前記或る体積の組織から発生したものではない多重反射信号を含む、ことと、
    前記多重反射信号に帰属する前記第2の反射された信号の部分を減弱させるために前記第2の反射された信号をフィルタリングすることと、
    前記第2のフィルタされた信号に基づいて前記或る体積の組織の第2の超音波画像を生成することと、
    前記多重反射信号に帰属する前記イメージの部分をさらに減弱させるために前記第1の超音波画像と前記第2の超音波画像とを比較することと、
    をさらに含む、請求項9に記載の方法。
  14. 前記第1の超音波画像と前記第2の超音波画像とを比較することが、平均値フィルタ、メジアンフィルタ、及び最小値フィルタのうちの少なくとも1つを用いて前記イメージをフィルタリングすることを含む、請求項13に記載の方法。
  15. 超音波画像における非定常のアーチファクトを減弱させる方法であって、
    第1の信号経路上に第1の信号を印加することによって前記アーチファクトを含む超音波画像データの第1の組を取り込むことと、
    第2の信号経路上に第2の信号を印加することによって超音波画像データの第2の組を取り込むことと、
    前記アーチファクトを除去するために前記データの第1の組と前記データの第2の組を組み合わせることと、
    を含む、方法。
  16. 前記データの第1の組と前記データの第2の組とを組み合わせることが、ローパスフィルタを用いて前記データの第1の組及び前記データの第2の組をフィルタリングすることを含む、請求項15に記載の方法。
  17. 前記超音波画像データの第1の組の取り込みが第1の時間間隔で行われ、前記超音波画像データの第2の組の取り込みが前記第1の時間間隔から離間した第2の時間間隔で行われる、請求項15に記載の方法。
  18. 前記第2の信号経路上に第2の信号を印加することが、前記第1の信号経路とは異なる経路長又は経路ルートを有する第2の信号経路上に第2の信号を印加することを含む、請求項15に記載の方法。
  19. 前記アーチファクトを除去するために前記データの第1の組と前記データの第2の組を組み合わせることが、前記データの第1の組及び前記データの第2の組に、平均値フィルタ、メジアンフィルタ、又は最小値フィルタのうちの少なくとも1つを適用することを含む、請求項15に記載の方法。
  20. 前記データの第1の組と前記データの第2の組との間の分散に少なくとも部分的に基づいて前記平均値フィルタ、前記メジアンフィルタ、又は前記最小値フィルタのうちの少なくとも1つを選択することをさらに含む、請求項19に記載の方法。
  21. コンピューティングシステムによって実行される場合に、前記コンピューティングシステムに、
    第1の信号経路上に第1の信号を印加するようにイメージングトランスデューサに命令することと、
    前記第1の信号に基づいてアーチファクトを含む超音波画像データの第1の組を受信することと、
    第2の信号経路上に第2の信号を印加するように前記イメージングトランスデューサに命令することと、
    前記第2の信号に基づいて超音波画像データの第2の組を受信することと、
    前記アーチファクトを除去するために前記データの第1の組と前記データの第2の組とを組み合わせることと、
    を含む動作を行わせるコンピュータで実行可能な命令を格納している物理的コンピュータ可読記憶媒体。
  22. 前記動作が、平均値フィルタ、メジアンフィルタ、又は最小値フィルタのうちの少なくとも1つを用いて前記超音波画像データの第1の組及び前記超音波画像データの第2の組をフィルタリングすることをさらに含む、請求項21に記載の物理的コンピュータ可読記憶媒体。
  23. 前記動作が、ローパスフィルタを用いて前記データの第1の組及び前記データの第2の組をフィルタリングすることをさらに含む、請求項21に記載の物理的コンピュータ可読記憶媒体。
  24. 前記第2の信号経路上に前記第2の信号を印加するように前記イメージングトランスデューサに命令することが、前記第1の信号経路とは異なる経路長又は経路ルートを有する第2の信号経路上に前記第2の信号を印加するように前記イメージングトランスデューサに命令することを含む、請求項21に記載の物理的コンピュータ可読記憶媒体。
  25. 前記動作が、前記イメージングトランスデューサを前記第1の信号の印加と前記第2の信号の印加との間で固定量、平行移動することをさらに含む、請求項21に記載の物理的コンピュータ可読記憶媒体。
  26. 超音波イメージングシステムであって、
    或る体積の組織に超音波パルスを送達するように構成されたイメージングトランスデューサと、
    前記或る体積の組織から反射された信号を受信するように構成された受信器であって、前記反射された信号は、前記或る体積の組織から発生していない多重反射信号を含む、受信器と、
    前記多重反射信号に帰属する前記反射された信号の部分を減弱させるために前記反射された信号をフィルタし、前記フィルタされた信号に基づいて前記或る体積の組織の超音波画像を生成するようにプログラムされたコントローラと、
    を備える、システム。
  27. 前記イメージングトランスデューサが、固定の又は可変の時間間隔だけ分離された超音波パルスを送達するように構成される、請求項26に記載のシステム。
  28. 前記反射された信号をフィルタするようにプログラムされたコントローラが、ローパスフィルタを用いて前記反射された信号をフィルタするようにプログラムされたコントローラを含み、前記ローパスフィルタが、回帰フィルタ、有限インパルス応答フィルタ、無限インパルス応答フィルタ、及び帯域消去フィルタのうちの少なくとも1つを含む、請求項26に記載のシステム。
  29. 前記イメージングトランスデューサを動かし、それにより前記反射された信号の前記長さを修正するように構成された移動機構をさらに備える、請求項26に記載のシステム。
  30. 前記反射された信号をフィルタするようにプログラムされたコントローラが、平均値フィルタ、メジアンフィルタ、又は最小値フィルタのうちの少なくとも1つを用いて前記反射された信号をフィルタするようにプログラムされたコントローラを含む、請求項26に記載のシステム。
JP2014555825A 2012-02-02 2013-02-04 超音波画像における非定常多重反射を抑制するためのフィルタリングシステム Expired - Fee Related JP6161129B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201261594244P 2012-02-02 2012-02-02
US61/594,244 2012-02-02
PCT/US2013/024663 WO2013116855A1 (en) 2012-02-02 2013-02-04 Filtering systems and methods for suppression of non-stationary reverberation in ultrasound images

Publications (2)

Publication Number Publication Date
JP2015506258A true JP2015506258A (ja) 2015-03-02
JP6161129B2 JP6161129B2 (ja) 2017-07-12

Family

ID=47750805

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2014555825A Expired - Fee Related JP6161129B2 (ja) 2012-02-02 2013-02-04 超音波画像における非定常多重反射を抑制するためのフィルタリングシステム

Country Status (6)

Country Link
US (1) US10713758B2 (ja)
EP (1) EP2810248A1 (ja)
JP (1) JP6161129B2 (ja)
KR (1) KR20140121859A (ja)
CN (1) CN104160423A (ja)
WO (1) WO2013116855A1 (ja)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3382423A1 (en) * 2017-03-27 2018-10-03 Koninklijke Philips N.V. Methods and systems for filtering ultrasound image clutter
JP7333273B2 (ja) * 2017-05-11 2023-08-24 コーニンクレッカ フィリップス エヌ ヴェ 超音波診断画像における反響アーチファクト打ち消し
EP3622320A1 (en) * 2017-05-11 2020-03-18 Koninklijke Philips N.V. Methods and systems for controlling the generation of a compound ultrasound image
US11334974B2 (en) 2017-08-16 2022-05-17 Koninklijke Philips N.V. Systems, methods, and apparatuses for image artifact cancellation
WO2020011589A1 (en) * 2018-07-11 2020-01-16 Koninklijke Philips N.V. Ultrasound imaging system with pixel extrapolation image enhancement
CN113556979A (zh) * 2019-03-19 2021-10-26 奥林巴斯株式会社 超声波观测装置、超声波观测装置的工作方法以及超声波观测装置的工作程序
CN114727807A (zh) * 2019-11-21 2022-07-08 皇家飞利浦有限公司 减少超声图像中的混响伪影以及相关联的设备、系统和方法
CN111436971B (zh) * 2020-04-20 2022-06-10 汕头市超声仪器研究所股份有限公司 一种超声混响伪像抑制方法

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH03114451A (ja) * 1989-09-29 1991-05-15 Terumo Corp 超音波診断装置
JPH03114452A (ja) * 1989-09-29 1991-05-15 Terumo Corp 超音波診断装置
JPH04338462A (ja) * 1991-05-16 1992-11-25 Hitachi Ltd 超音波治療装置
JPH09103434A (ja) * 1995-03-31 1997-04-22 Toshiba Corp 超音波治療装置
JPH09182749A (ja) * 1995-12-28 1997-07-15 Aloka Co Ltd 超音波診断装置
JPH10286256A (ja) * 1997-04-15 1998-10-27 Fujitsu Ltd 超音波診断装置およびフィルタ
JP2006051356A (ja) * 2004-08-09 2006-02-23 General Electric Co <Ge> 空間合成画像生成のための距離依存重み付け
JP2006095151A (ja) * 2004-09-30 2006-04-13 Ge Medical Systems Global Technology Co Llc 超音波診断装置
JP2007236681A (ja) * 2006-03-09 2007-09-20 Toshiba Corp 医用画像診断装置及び医用画像診断装置の制御プログラム
JP2008284329A (ja) * 2007-05-21 2008-11-27 Kazuhiro Iinuma 超音波検査装置
WO2011001310A1 (en) * 2009-06-30 2011-01-06 Koninklijke Philips Electronics, N.V. Propagation-medium-modification-based reverberated-signal elimination

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5465723A (en) * 1994-06-02 1995-11-14 Vingmed Sound A/S Method and apparatus for ultrasound imaging
US5538004A (en) * 1995-02-28 1996-07-23 Hewlett-Packard Company Method and apparatus for tissue-centered scan conversion in an ultrasound imaging system
US6511426B1 (en) * 1998-06-02 2003-01-28 Acuson Corporation Medical diagnostic ultrasound system and method for versatile processing
US6210328B1 (en) * 1998-10-01 2001-04-03 Atl Ultrasound Ultrasonic diagnostic imaging system with variable spatial compounding
DE60113982D1 (de) 2000-01-31 2006-02-23 Tonni F Johansen Phasenabweichungskorrektur und pulsnachhallkorrektur bei der ultraschall-echographie
US6436041B1 (en) 2000-12-22 2002-08-20 Acuson Corporation Medical ultrasonic imaging method with improved ultrasonic contrast agent specificity
US6776759B2 (en) 2002-02-27 2004-08-17 Ge Medical Systems Global Technology Company, Llc Method and apparatus for high strain rate rejection filtering
US6905465B2 (en) 2002-04-05 2005-06-14 Angelsen Bjoern A. J. Corrections for pulse reverberations and phasefront aberrations in ultrasound imaging
US20050277835A1 (en) * 2003-05-30 2005-12-15 Angelsen Bjorn A Ultrasound imaging by nonlinear low frequency manipulation of high frequency scattering and propagation properties
US7311701B2 (en) * 2003-06-10 2007-12-25 Cierra, Inc. Methods and apparatus for non-invasively treating atrial fibrillation using high intensity focused ultrasound
US20050288587A1 (en) * 2004-06-25 2005-12-29 Yongrae Roh Drive machanism for mechanically scanned ultrasound transducers
EP1789136B1 (en) * 2004-09-16 2010-12-29 University of Washington Interference-free ultrasound imaging during hifu therapy, using software tools
US8699765B2 (en) 2006-11-13 2014-04-15 Siemens Medical Solutions Usa, Inc. Reducing jittering in medical diagnostic ultrasound imaging
US7780601B2 (en) * 2007-06-05 2010-08-24 Siemens Medical Solutions Usa, Inc. Adaptive clinical marker preservation in spatial compound ultrasound imaging
CN101305925B (zh) 2008-05-15 2010-11-03 上海理工大学 超声扫描探头
CN102136138B (zh) * 2011-03-28 2013-01-16 深圳市蓝韵实业有限公司 一种医学图像降噪与增强的处理方法

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH03114451A (ja) * 1989-09-29 1991-05-15 Terumo Corp 超音波診断装置
JPH03114452A (ja) * 1989-09-29 1991-05-15 Terumo Corp 超音波診断装置
JPH04338462A (ja) * 1991-05-16 1992-11-25 Hitachi Ltd 超音波治療装置
JPH09103434A (ja) * 1995-03-31 1997-04-22 Toshiba Corp 超音波治療装置
JPH09182749A (ja) * 1995-12-28 1997-07-15 Aloka Co Ltd 超音波診断装置
JPH10286256A (ja) * 1997-04-15 1998-10-27 Fujitsu Ltd 超音波診断装置およびフィルタ
JP2006051356A (ja) * 2004-08-09 2006-02-23 General Electric Co <Ge> 空間合成画像生成のための距離依存重み付け
JP2006095151A (ja) * 2004-09-30 2006-04-13 Ge Medical Systems Global Technology Co Llc 超音波診断装置
JP2007236681A (ja) * 2006-03-09 2007-09-20 Toshiba Corp 医用画像診断装置及び医用画像診断装置の制御プログラム
JP2008284329A (ja) * 2007-05-21 2008-11-27 Kazuhiro Iinuma 超音波検査装置
WO2011001310A1 (en) * 2009-06-30 2011-01-06 Koninklijke Philips Electronics, N.V. Propagation-medium-modification-based reverberated-signal elimination

Also Published As

Publication number Publication date
EP2810248A1 (en) 2014-12-10
JP6161129B2 (ja) 2017-07-12
WO2013116855A1 (en) 2013-08-08
CN104160423A (zh) 2014-11-19
US10713758B2 (en) 2020-07-14
KR20140121859A (ko) 2014-10-16
US20130204135A1 (en) 2013-08-08

Similar Documents

Publication Publication Date Title
JP6161129B2 (ja) 超音波画像における非定常多重反射を抑制するためのフィルタリングシステム
JP6188594B2 (ja) 超音波診断装置、画像処理装置及び画像処理方法
JP4931910B2 (ja) 超音波画像装置
JP5514911B2 (ja) 超音波撮像装置
JP2010183935A5 (ja)
JP2001245889A (ja) 超音波プローブおよび超音波診断装置
JP5398514B2 (ja) カラードプラ超音波診断装置
TW201922177A (zh) 具有高頻細節的超音波系統
US7666139B2 (en) Ultrasonic imaging apparatus and method
US9354300B2 (en) Ultrasound diagnostic apparatus and ultrasound image producing method
JP2007507271A (ja) 超音波イメージングにおけるスモールアンサンブル長によるクラッタフィルタリング
JP6169361B2 (ja) 超音波診断装置および輝度補正方法
JP5355327B2 (ja) 超音波診断装置
KR102263282B1 (ko) 이미지 이득을 조정하는 방법 및 시스템
JP7280713B2 (ja) 超音波診断装置
CN110433406B (zh) 超声聚焦定位图像优化方法
JP7578221B2 (ja) 高周波ディテールを有する超音波システム
JP2019025247A (ja) 超音波撮像装置、および、画像処理装置
KR101341092B1 (ko) 가변 인덕터를 이용한 초음파 진단 방법 및 장치
JP7015640B2 (ja) 超音波診断装置およびその制御方法
Suzuki et al. A study on ultrasound speckle reduction based on stochastic fluctuation of transmitted ultrasound beam
JP3663206B2 (ja) 超音波診断装置
JP2003339709A (ja) ドップラ信号処理装置および超音波診断装置
JP2002058670A (ja) 超音波診断装置
JP2023531979A (ja) 連続波ドプラ法のための超音波トランスデューサプローブに基づくアナログ-デジタル変換、並びに関連する装置、システム、及び方法

Legal Events

Date Code Title Description
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20150623

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20150804

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20151102

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20160204

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20160809

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20161109

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20170110

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20170509

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20170608

R150 Certificate of patent or registration of utility model

Ref document number: 6161129

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees