JP2015532848A - 画像組織(tissue)動作推定 - Google Patents
画像組織(tissue)動作推定 Download PDFInfo
- Publication number
- JP2015532848A JP2015532848A JP2015532513A JP2015532513A JP2015532848A JP 2015532848 A JP2015532848 A JP 2015532848A JP 2015532513 A JP2015532513 A JP 2015532513A JP 2015532513 A JP2015532513 A JP 2015532513A JP 2015532848 A JP2015532848 A JP 2015532848A
- Authority
- JP
- Japan
- Prior art keywords
- motion
- phase
- matrix
- pair
- estimate
- 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.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/485—Diagnostic techniques involving measuring strain or elastic properties
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8979—Combined Doppler and pulse-echo imaging systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52023—Details of receivers
- G01S7/52036—Details of receivers using analysis of echo signal for target characterisation
- G01S7/52042—Details of receivers using analysis of echo signal for target characterisation determining elastic properties of the propagation medium or of the reflective target
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8979—Combined Doppler and pulse-echo imaging systems
- G01S15/8986—Combined Doppler and pulse-echo imaging systems with measures taken for suppressing velocity ambiguities, i.e. anti-aliasing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radiology & Medical Imaging (AREA)
- Public Health (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Heart & Thoracic Surgery (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biophysics (AREA)
- Acoustics & Sound (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
動きプロセッサー(118)は、予め処理されたエコーの連続するフレームのペア間の動きを反復して推定し、動き推定器(306)は、少なくとも1つの反復に基づいて、推定された動きを生成する。方法は、少なくとも1つの反復を介して、予め処理されたエコーの連続するフレームのペア間の組織動きを反復して推定することを含む。【選択図】図1
Description
以下は、一般に画像に関し、より具体的には、画像組織動作推定に関し、ここでは特に、超音波画像について記述される。
超音波(US)画像のような画像は、検査対象の内部の特性について有用な情報を提供してきた。組織動作推定は、エラストグラフィー(組織弾性画像法:elastography)のような画像に関連して使用されており、歪み画像(strain imaging)、歪み率画像(strain rate imaging)および弾力画像(elasticity imaging)とも呼ばれる。そのような画像では、軟組織の硬直(stiffness)または歪み(strain)の画像が、健康か正常な軟組織の周囲に比べてより硬直した不健康な組織(腫瘍、壊死組織など)を検出し、分類するために使用される。これについては、機械的圧縮又は振動が関心領域に適用され、歪み(strain)が周囲の組織よりも少ないので、不健康な組織は周囲の組織よりも圧縮が少ない。
軸の動き(motion)(超音波ビームの方向に沿った動き)で、ブロック一致、位相シフトおよび周波数推定は、動きを推定するために使用される。スペックル(speckle)トラッキングとして知られているタイムシフト技術が、ブロック一致に使用される。スペックルトラッキングで使用されるデータは、包絡線検出されたデータ(位相情報なし)又は無線周波数(RF)データである。RFデータは包絡線検出の前に現れれるデータである。これについては、データ(カーネル)の1次元(1D)又は2次元(2D)の組が、最初のフレームから選ばれる。このカーネルは、予め定義されたサーチ範囲内で次のフレームからの同じサイズのデータセットと比較される。ベストな一致が見つけられた場合、置換(displacement)が、2つの領域の位置の差から推定される。
異なる一致アルゴリズムは、ブロック一致と使用される。一致アルゴリズムの例は、絶対差の合計、2乗差の合計及び最大相互相関を含む。望ましくないことに、これらの技術は、高いサンプリング周波数(例えば、8倍の搬送周波数)を要求し、それはナイキスト要求の4倍である。更に、絶対差の合計のアプローチ、及び2乗差の合計のアプローチは、最小誤り識別の影響を受け易く、特に、相互相関がRFデータ上で行なわれる場合に、最大相互相関アプローチは、誤りピーク識別の影響を受け易い。
エラストグラフィーは、また、推定のサブサンプルの精密性を要求する。これについて、多項式近似が、最大相互相関関数又は絶対差の合計関数にフィットされる。全ての行列が乗算される場合、相互相関の1つの実例は最も速く、その結果(それも行列である)は、一時的なマトリックスと比較される。相互相関に内挿(intepolation)を加えることは、異なる遅れの結果がメモリに維持されることを意味し、それはRAMのニーズを増加させる。メモリアクセスは、しばしな現代のシステムのボトルネックである。高精度置換推定のためのブロック一致を使用することは、大量のデータの転送を必要とする。
位相シフト推定は、RFデータ信号又は複素数値(complex-valued)ベースバンド信号で行なわれ、同相(in-phase)及び直交位相信号としても知られている。位相シフト置換推定の基礎は、測定された信号が狭周波数帯であり、x(t) = a(t) exp(-jωt)、角周波数ω=2πfとして記述されることができるという仮定である。超音波信号の周波数fは、しばしば組織における送信されたパルスの搬送周波数である。時間tはカレントラインの放射のスタートから測定される。古典的推定手続きは、フレーム1およびフレーム2からの同じ深さからサンプルをとることである。2つのフレーム間の違いだけが動きによると仮定される。
さらに、フレーム2からの信号x2(t)が、フレーム1からの信号x1(t)の時間シフトバージョン、すなわち、x2(t) = x1(t-tm) であると仮定される。tmがポジティブな場合には、動きは、変換器(transducer)から離れている。狭帯域モデルを x1(t) = a(t)と仮定する。x1(t)とx2(t)との間の遅れゼロ相互相関は、
である。包絡線a(t)は、典型的にはゆっくり変化する関数であり、時間シフトは、複素相関関数R12(0):φm = <(R12(0)) = -jωtm の位相から得られる。深さと時間の関係は、t = 2z/ cによって与えられ、ここでtは時間、zは深さ、cは音の速度である。
位相シフトφmに対応する深さuの置換は、u = (-c/4πf)φm である。周波数は一定で選ばれるが、高精度な推定については、所与の深さでの信号の瞬時平均(mean)周波数を反映しなければならない。実際、R12(0)の推定は、フレーム内の同じ場所に位置した、推定カーネルの内部のフレーム2からのすべてのサンプルに推定カーネルの内部のフレーム1からのすべてのサンプルを乗算することにより行なわれる。次に、最終的な積は共に加えられる。アークタンジェントは位相シフトを決定するために使用される。不運にも、位相シフト方法は、エイリアシングの影響を受け易い。位相はπ周辺をラップする(wrap)。したがって、ロバスト位相アンラッピングアルゴリズムが必要とされる。
この方法の別の問題は、推定の変化(variance)が動きにつれて増加するということである。すなわち、小さい動きは確実に検出されるが、大きな動きは確実に検出されない。2つのフレームからのIQ信号の位相が引かれる場合、変化(variants)は、R12のゼロ位相サーチ及ぶ位相分離を含んでいる。これらの変化は、エイリアシングの問題を克服し、変化を減少しようとする。変換器表面からスタートするそれらの推定は、次に、蓄積された位相シフトの軌跡を追う。しかしながら、間違いが1つの深さで起こった場合、そのエラーはフレームの終了まで伝搬する。
組み合わされた自己相関アプローチは、ブロック一致および位相推定に基づいて推定する。ブロック一致および位相推定の双方について、推定が必要な場合、相互相関関数は画素毎に推定される。第1のフレームからの相関ウィンドウは、第2のフレームとは異なる置換での相関ウィンドウに対してマッチされる。すべての置換 (k, l) について、標準化された相互相関関数が推定される。R12(k,l)の最大値は維持される。
であり、ここで、fsはサンプリング周波数、kは深さでの例である。このアプローチは、R12が、正しい期間の倍数内のサンプルで最大を持つことが期待されるので、エイリアシングを克服するために使用されることができる。しかしながら、R12の正しいピークの検出は、広帯域のパルスを必要とし、信頼できるφの推定は、狭帯域パルスを必要とし、誤ったピークが選択された場合、その結果は、画像に現われ、このアプローチは、”包絡線圧縮”の利点を利用することができない。
歪み(strain)は、動きの一次微分であるので、エラストグラフィーは、特に、動きの推定の変化に敏感である。カーネル(相関ウィンドウ)内の動きは一様ではない。推定手続きは、相関ウィンドウの中心の動きを推定することを意図する。ウィンドウの端に強い散乱がある場合、それは動作推定を支配する。周期的な変化が、動きに推定おいて存在し、周期的な変化は、歪み画像における”ゼブラ”アーティファクトとして知られているアーティファクトを引き起こす。ゼブラアーティファクトの影響を低減するために、文献は、推定に先立ってRF信号の包絡線を圧縮することを示唆する。文献から知られている圧縮関数は、以下の通りである。xc(t) = sgn(s(t))×log10(1+cs×abs(x(t)) ここで、sgn(x)は、値の符号(sign)を返し、csは圧縮のレベルを制御する定数である。
本願のアスペクトは、上記及びその他の問題に取り組む。1つのアスペクトでは、動きプロセッサー(118)は、予め処理されたエコーの連続するフレームのペア間の動きを反復して推定する動き推定器(306)を含み、動き推定器(306)は、少なくとも1つの反復に基づいた、推定された動きを生成する。
他のアスペクトでは、方法は、少なくとも1つの反復を通して、予め処理されたエコーの連続するフレームのペア間の組織の動きを反復的に推定することを含む。
他のアスペクトでは、コンピュータ可読記憶媒体は、コンピュータ実行可能な命令でエンコードされ、プロセッサによって実行された場合に、プロセッサに、予め処理されたエコーの連続する圧縮されたフレームのペアと、遅れ−ゼロ相互相関及び平均(mean)周波数に基づいて、最初の動き推定を決定させ、初期動き推定に基づいて非ゼロ相互相関遅れを決定させ、要素ごとに、非ゼロ相互相関遅れでの予め処理されたエコーの圧縮された連続するフレームのペアにおける要素を乗算させ、積行列を生成させ、前記積行列を2D重みづけ関数で畳み込ませ、端での同相推定における非連続を取り除かせ、デグリッチされた畳み込み行列を生成させ、デグリッチされた畳み込み行列の位相シフトを決定させ、決定された位相シフトに基づいて、デグリッチされた畳み込み行列の位相をアンラップ(unwrap)させ、アンラップされた位相を動き推定に変換させる。
当業者は、添付された記述を読み理解することにより、本願の他のアスペクトを認識することができる。
本願は、例示によって図示されており、添付図面の図に限定されるものではなく、同様の参照は、同様の要素を示している。
図1は、超音波画像システムのような画像システム102の例を概略的に示す。図示された画像システム102は、コンソール106とインターフェイスされた変換器アレイ104を含む。変換器アレイ104は、超音波信号を視界に送信し、かつ視界における組織に対する送信超音波信号の相互作用に応じて生成されたエコー信号を受け取るように構成された多くの変換器要素(見えない)を含む。変換器アレイ104は、線形、曲げられ、及び/又は、形作られ、完全に密集され、希薄化され、及び/又は、それらのコンビネーション、1次元(1D)又は2次元(2D)などであっても良い。
コンソール106は、超音波信号を送信するために変換器要素を動作させる送信回路108を含む。1つの実例では、これは、異なる機械的振動を、エラストグラフィー応用のために、少なくとも2つのフレームの視界領域における組織に適用するために変換器要素を動作することを含んでいる。受信回路110は、送信された超音波信号に応じて生成されたエコーの組を受信する。このエコーは、一般に、視界の走査範囲で発せられた超音波信号と組織(例えば器官、腫瘍など)との間の相互作用の結果である。エラストグラフィーについては、エコーは、データのフレーム間の組織の動き又は置換情報を含む。ビームフォーマ112は、受信エコーを処理し、例えば、エコーに時間遅れと重みづけを適用し、そのエコー結果を加算して、RF信号を生成又はラインを走査する。フレームは、視界の2D画像をともに形成する走査ラインの組から成る。
コンソール106は、さらに、RF信号を予め処理するプリ・プロセッサ114を含む。プリ・プロセッサ114の制限しない例が、図2に示される。この例において、プリ・プロセッサ114は、例えば、バンドパス・スライディングフィルタを介して、RF信号をフィルタするフィルタ202を含む。一般に、そのようなフィルタは信号対雑音比を最大にし、RF信号の特性と一致するように、深さの関数としてその特性を変更する。変換器204は、ヒルバート(Hilbert)又は他の変換を適用し、実際の入力から複素信号を作り出す。周波数シフタ206は、例えば、複素信号を所定の混合周波数と混合することにより、ベースバンド(ゼロヘルツ、又は0.0 Hz)に向けて複素信号のスペクトルをシフトする。デシメータ208は、情報の損失をすることなく、データレートを維持しながら、この信号(例えば2以上のファクタにより)をデシメートする。プリ・プロセッサ114は、複素ベースバンド信号を出力する。
図1に戻り、動きプロセッサー118は、複素ベースバンド信号の隣接したフレームのペアを処理する。そのような処理は、データのフレーム間の組織の動きを推定することを含む。下記に詳細に述べるように、1つの実例では、処理は、相互相関を使用する位相シフト推定を通した動き推定の反復改善を含む。このようなアプローチは、低減されたデータセットで高精度を実現しながら、滑らかな動きマップ(すなわち、組織の動き推定の減少した変化)、低減されたノイズアーティファクト、及び/又は、エラストグラフィーを有する増大したコントラスト(すなわち、構造の検出性)、及び/又は他の画像に導く。動き推定器306は、コンソール106の一部として示されているが、動き推定器306は、例えば、コンピューター及び/又は他のコンピューティング装置の一部として、コンソール106から、択一的に遠隔にあっても良いことが認識される。
レンダリングエンジン120は、処理されたフレームに基づいて、エラストグラフィー又は他の画像を少なくとも生成するように構成される。エラストグラフィー画像は、ディスプレイ122及び/又は他のディスプレイを介して視覚的に示され、他の装置に送られて格納され、及び/又は使用される。ユーザ・インターフェース(UI)124は、1つ以上の入力装置(例えば、マウス、キーボードなど)を含み、システム102とのユーザ相互作用を可能にする。コントローラー126は、画像システム102の種々のコンポーネントを制御する。例えば、そのような制御は、Aモード、Bモード、C-プレーン、及び/又は他のデータ取得モードのために変換器アレイ104の変換器要素の個々又はグループを動作し又は駆動し、送信され及び/又は受信された信号を、ステアリング及び/又はフォーカスすることなどを含む。
コンソール106は、物理メモリーおよび他の非一時的媒体のようなコンピューター可読記憶媒体上でエンコードされ、又は埋め込まれた1つ以上のコンピューター可読命令を実行する1台以上のプロセッサー(例えば、中央演算ユニット(CPU)、グラフィック処理ユニット(GPU)など)を含んでも良い。付加的又は択一的に、命令は信号、搬送波、および他の一時的又は非一時的コンピュータ可読媒体に伝送されることができる。1つの実例では、命令の実行は、1台以上のプロセッサーに関連して、画像システム102のビームフォーマ112、プリ・プロセッサ114、動き推定器306、レンダリングエンジン120及び/又は画像システム102の他のコンポーネントの1つ以上を実行する。
簡潔に上に議論したように、動きプロセッサー118は、複素ベースバンド信号の隣接したフレームのペアを処理し、反復アプローチを使用してフレーム間の組織の動きを推定する。図3は、動きプロセッサー118の例を概略的に示す。
包絡線コンプレッサー302は、入力として、複素ベースバンド信号xi(m,n)及びxj(m,n)の2つの隣接するフレームを受信し、ここで、i及びjは2つの隣接するフレーム(例えば、フレーム1及びフレーム2、フレーム2及びフレーム3、など)のフレームインデックスであり、mは走査ラインに沿ったサンプルインデックスであり、nは、走査ラインインデックスである。xi(m,n)及びxj(m,n)は、行列として処理され、ベクトル及び/又は個々の要素として、mは行列の行のインデックス、nは行列の列のインデックスである。包絡線コンプレッサー302は、位相を維持し、xi(m,n)及びxj(m,n)を生成する一方、包絡線を圧縮し、ここで、cは、圧縮されたことを示す、それは”ゼブラ”アーティファクトを減少することを容易にする。
限定されない例として、入力信号は複素サンプルから作られ、x = xr + jxi であり、ここで、xrは実(同相)成分であり、xiは虚(直交)成分である。圧縮された信号は、xc = xcr + jxci として、
として表わされる。1つの実例では、定数1(1)は、ゼロによる除算を防ぐために付け加えられる。CPU及び/又はGPUの内蔵命令は、平方根を計算するために使用されることができる。
周波数推定器304は、入力として、圧縮されたxic(m,n)を受信し、平均周波数f(m)を、例えば、各深さに基づいて、推定する。超音波の減衰は、周波数依存であり、したがって、深さの関数として、平均(中心)周波数は変化する。結果の中心周波数は、送信パルスの中心周波数、変換器のインパルス応答、周波数依存の減衰及びスライディングフィルタの特性の関数である。スキャナは、アプリケーションに依存して、異なる送信パルス及びスライディングフィルタを使用する。更に、周波数依存の減衰は、走査されている組織の関数である。
一般に、新しいフレームが入力されるごとに、中間周波数f(m)が走査された組織及びシステムセットアップと一致するように推測される。入力は、行列を示す離散複素信号xic(m,n)であり、ここで、インデックスmは軸方向(深さ)であり、インデックスnは横方向である。周波数の推定は、遅れ1の複素自己相関関数
の位相から行なわれ、ここで、Rii(m,n;1,0)は画素/位置(m,n)での複素自己相関であり、遅れ1は軸方向(mに沿う)、遅れ0は、横方向(nに沿う)である。
平均ウィンドウは、V行(サンプル)及びW列(ライン)からなり、上付き文字*は、複素共役を示す。自己相関関数は、複素であり、Rii(m,n;1,0) = Riir(m,n;1,0) + jRiii(m,n;1,0)である。 画素(m、n)の瞬時周波数は、
は、最初に平均化され、次に、f(m)がn次(例えば、2nd, 3rdなど)として得られ、多項式が最小2乗エラーフィッティングを介して得られる。
図4は、位相アンラップアプローチの例を示す。図4では、位相アンラップは、位相アンラッパ 400によって実現される。位相アンラッパ 400は、入力として、入力マトリクスΦiを受信し、マトリックスΦ0を出力する。この例において、処理はカラム毎に行なわれる。しかしながら、他のアプローチ(例えば、行、要素、マトリックス全体などによって)も、ここでは考えられる。微分器402は、入力位相φi(m,n)を微分する。ここで、Δφ(m,n) = φi(m,n) - φi(m-1,n)である。コンパレーター404は、エイリアシングを決定するために、結果を所定のしきい値と比較する。1つの実例では、しきい値は±βπであり、ここで、βは1以下の定数(β≦1)である。コンパレーター404の出力は、エイリアシングを示す。エイリアシングが表われるごとに、差の符号に依存して、±2πの値が位相に加えられなければならない。積分器406は、調節値を積分する。加算器408は、積分値及び入力位相を加算する。
図3に戻り、動き推定器306は、位相シフト推定を使用し、平均周波数推定f(m)を介して位相を動きに変換して、xic(m,n)及びxjc(m,n)に基づいて、動きul(m,n)を推定する(ここで、lは、0から1,2,3,5,17,100などのようなまでの反復インデックスである。)。より詳しく以下に記述されるように、1つの実例では、最初(l=0)の動作推定u0(m,n)は、遅れ(0,0)、0個のサンプル及び0ラインオフセットで、反復(l = 1, 2, 3など)毎に、相互相関を使用し、動き推定器306は、遅れkl(m,n)で2つのフレーム間の位相シフトを見つけ、この位相シフトを使用して、新しい動きul(m,n)を推定し、より改善され又は正確な動き推定ul(m,n)を生成する。相互相関推定は、複素数値であり、相互相関関数の角度(位相)は動きに一致する。
表面フィルタ308は、表面を動き推定ul(m,n)にフィットし、変化を減らし、アウトライアーを取り除き、その結果、雑音動きが動き推定を推定し、滑らかにする。1つの実例では、表面が最小二乗2D多項式フィッティングを使用してフィルタされ、ここで、多項式係数はP = Gul、Gは3次の表面に対しての最小二乗フィットを解くことから導出された行列である。多項式の他の順序(例えばn次)は、正確及び計算上の要求の間のトレードオフに依存して、使用される。ここで、ulは、動きの推定を備えたマトリックスである。その表面は、ufit(m,n) = P(0) + P(1)n + P(2)m + P(3)n2 + P(4)mn + P(5)m2 + P(6)n3 + P(7)n2m + P(8)nm2 + P(9)m3 として表わされる。表面のフィルタ210は、さらに異なる遅れ間の滑らかな端を識別し、それは、下記に議論されるように、位相アンラップを容易にする。
量子化器310は、フィットされた動き推定を量子化し、2つのフレーム間の相互相関がul(m,n)よりも正確な推定を与えるk(m,n)の遅れを見つける。例として、フィットされた表面が量子化され、次の反復で相互相関を計算する場所での遅れを与える。
決定コンポーネント212は、入力として、現在の遅れ及び新しい遅れを受信し、動き推定は反復を通して改善されたか、或いは最終的な動き推定ul(m,n)が種々の停止基準に基づいて出力されたかを決定する。このような基準は、これに限定されるものではないが、所定数の反復(例えば、l = 1 又は1つの反復及び動き推定)、現在及び新しい遅れの間の予め決定された差異、予め決定されたタイムインターバル、オンデマンドベースのユーザ入力及び/又は他の停止基準を含んでも良い。
上記概略的に述べたように、動き推定器306は動きul(m,n)を推定する。図5は、反復l=0における最初の動きを推定するサンプルコンポーネントを示し、図6は、l-1, ... Sにおける続く移動推定を推定するサンプルコンポーネントを示す。
乗算器502は、2つの入力行列xic(m,n)及びxjc(m,n)を要素毎に乗算し、Pij(m,n)を生成する。コンボルバー 504は、平均ウィンドウH(v,w)でPij(m,n)を畳み込み(convolve)する。平均ウインドウ H(v,w)が2Dの長方形関数である場合、畳み込みは、サンプル当たり、4つの(4)加算及び2つの(2)減算に減らされる-移動平均の再帰的実行。位相シフト決定器506は、画素(m、n)で、
としての最初の動き推定に変換する。
図6では、計算機602は、本例において見えていないが、乗算器502及びコンボルバー 504を含む。乗算器502は、要素毎に、遅れ(k,m)に基づいて、2つの入力行列xic(m,n)及びxjc(m,n)を乗算する。このために、k(m,n)の最小値及び最大値が決定される。乗算器502は、miin(k(m,n))からmax(k(m,n))まで繰り返す。全ての遅れχについて、全ての要素xicが、それぞれの遅れχについて、xjcの要素の複素共役によって乗算される。コンボルバー 504は、フィルタH(w,v)で出力行列Pを畳み込みし、Rijにおける位置(m,n)での要素では、k(m,n) = χ が更新される。
として、各相互相関の位相を決定する。グリッチ除去器604は、その要素がk(m,n)である遅れKを有する行列の推移を検出し、取り除くことにより、グリッチを削除する。Φtmp2= unwrap_with_lags(Φtmp1, K)k(m,n)において遷移が発生する境界では、エイリアシングによらず、アンラップアルゴリズムによって捕らえられないステップがあることに留意すべきである。位相アンラッパ 508は、ここに記述された、アンラップアプローチ及び/又は他のアプローチを使用して、位相をアンラップする。遅れの推移が間違った画素にある場合、遅れは同じであるが、位相ラップが生ずる。そのため、ここに議論された、アンラップアルゴリズムが適用される。Φ1 = unwrap(Φtmp2)位相から距離への変換器510は、
として、動きを決定する。
図7は、グリッチ削除アプローチを有する適切な位相アンラップの例を示す。図7は、702で示されるような遅れの包含を備えた図4と同様である。遅れを有する位相アンラップは、同じラインに沿って得られるデータで動作する。ラインインデックスは、簡潔さのために省略される。入力は、推定された位相Φi(m)及び推定された位相での遅れk(m)である。示されるように、遅れと位相の両方が微分される。遅れは整数であり、必要な情報は変化の大きさである。そのような変化が検出される場合、このステップはその推移での位相ジャンプと等しい。メインパスは、位相差をしきい値と比較することにより、位相ラップを検出する。その結果は、位相がラップする位置及び位相ラップの大きさを示す一連のパルスである。これら一連のインパルスは積分され、階段状の位相補償関数を生成し、入力位相に加えられる。
図8は、エラストグラフィー画像に関連する動きを推定する方法を示す。
以下の動作の順番は、例示的な目的のためであり、限定されないいないことに注意すべきである。
そのため、1つ以上の動作は、異なる順番で行なわれることができ、限定されるものではないが、同時に行なわれることも含まれる。更に、1つ以上の動作は省略することができ、及び/又は、1つ以上の他の動作を加えることができる。
一般に、最初の動き推定(l=0)は、ゼロ(0)へセットされた全ての遅れk(m,n)を有する、その後の動き推定(l>0)の特定のケースである。最初の及び後の動き推定のためのデータは、同じ又は異なるモジュールによって処理されることができる。
800では、k(m,n), Rij(m,n),φ(m,n),ul(m,n)及びl がゼロ(0)に初期化される。
802では、ここで述べられたように、f(m)が決定される。
804では、ここで述べられたように、Rij(m,n;kl(m,n))が決定される。
806では、ここで述べられたように、φ(m,n)が決定される。
810では、ここで述べられたように、φ(m,n)がアンラップされる。
812では、ここで述べられたように、ul(m,n)が決定される。
814では、表面ufit(m,n)がul(m,n)にフィットされる。
816では、ここで述べられたように、新しい遅れkl+1(m,n)が決定される。
818では、停止基準がチェックされる。この例では、停止基準は、ここで述べられたように、遅れkl(m,n)と新しい遅れkl+1(m,n)との間の絶対差に基づく。
820では、止まる基準が満たされない場合、ここで述べられたように、k(m,n)は新しい遅れkl+1(m,n)にセットされ、動作804-816が繰り返される。
822では、停止基準が満たされる場合、ul(m,n)が最終動き推定として出力される。最終動き推定は、エラストグラフィー及び/又は他の画像を生成するために使用されることができる。
上記のものはコンピューター可読命令により実行されることができ、コンピューター可読命令は、コンピューター可読記録媒体上でエンコードされ、埋め込まれ、コンピュータープロセッサーによって実行された場合に、プロセッサーに述べられた動作を実行させる。さらに、又は択一的に、コンピューター可読命令の少なくとも1つは、信号、搬送波又は他の一時的媒体によって実行される。
本願は、様々の実施形態に関して記述されている。修正と代用が、本願を読むことで他者に生ずる。本願発明は、そのような全ての修正及び代用が含まれ、添付された請求項及びそれらの均等物の範囲内の限りのものを含むように構成されることが意図されている。
Claims (20)
- 予め処理されたエコーの連続するフレームのペア間の動きを反復的に推定する動き推定器を具備する動きプロセッサー(118)。
- 予め処理されたエコーの連続するフレームの各ペアの包絡線を圧縮する包絡線コンプレッサー(302)と、
前記予め処理されたエコーの連続するフレームの平均周波数を推定する周波数推定器(304)とをさらに具備し、前記動き推定器は、予め処理されたエコーの連続するフレームのペア、遅れゼロ相互相関及び平均周波数に基づき最初の動き推定を生成する、
請求項1記載の動きプロセッサ。 - 前記予め処理されたエコーの連続するフレームの各ペアは行列であり、
前記動き推定器は、
前記予め処理されたエコーの連続するフレームのペアにおける要素を要素毎に乗算し、第1の積行列を生成する乗算器(502)と、
前記第1の積行列を所定の変換で畳み込むコンボルバー(502)と、
前記畳み込み行列の位相シフトを決定する位相シフト決定器(506)と、
前記決定された位相シフトに基づいて、前記畳み込み行列の位相をアンラップするアンラッパ(508)と、
前記アンラップされた位相を前記最初の動き推定に変換する位相変換器(510)と
を具備する、請求項2記載の動きプロセッサ。 - 前記変換は、ヒルバート変換又は2D矩形関数である、請求項3記載の動きプロセッサ。
- 最初の動き推定に、n次の表面を最初の動き推定にフィットする表面フィルタ(308)を具備する請求項3又は請求項4のいずれかに記載の動きプロセッサ。
- 前記表面フィルタは、n次の表面を2D多項式フィッティングでフィットする請求項5記載の動きプロセッサ。
- 少なくとも1つの表面は、最初の動き推定を滑らかにし、又は誤り変化及び雑音を除去する請求項5又は請求項6いずれかに記載の動きプロセッサ。
- n次表面にフィットされた最初の動き推定を量子化し、非ゼロ相互相関遅れを決定する量子化器(310)をさらに具備する請求項5乃至請求項7のいずれかに記載の動きプロセッサ。
- 前記動き推定器は、予め処理されたエコーの連続するフレームのペア及び決定された非ゼロ相互相関遅れに基づいて、予め処理されたエコーの連続するフレームのペア間の動きを推定する、請求項8記載の動きプロセッサ。
- 前記非ゼロ相互相関遅れは行列であり、前記動き推定器は、
乗算器と、前記乗算器は、要素ごとに、前記予め処理されたエコーの連続するフレームのペアにおける要素及び前記非ゼロ相互相関遅れを乗算し、第2の積行列を生成し、
コンボルバーと、前記コンボルバーは、前記所定の変換で第2の積を生成し、
境界で位相推定のステップを除去するグリッチ除去器(604)と、
位相シフト決定器と、前記位相シフト推定器は、デグリッッチされた畳み込み行列の位相シフトを決定し、
アンラッパと、前記アンラッパは、前記決定された位相シフトに基づいてデグリッチされた畳み込み行列の位相をアンラップし、
位相変換器と、前記位相変換器は、アンラップされた位相を後の位相推定に変換し、
とを具備する計算機(602)
を有する請求項9記載の動きプロセッサ。 - 現在の動きが、停止基準に基づく反復を介して改良されるか否かを決定する決定コンポーネント(312)をさらに具備する請求項2乃至請求項10のいずれかに記載の動きプロセッサ。
- 前記動きプロセッサーは、超音波画像システムの一部である請求項2乃至請求項11のいずれかに記載の動きプロセッサ。
- 予め処理されたエコーの連続するフレームのペア間の組織動きを少なくとも1つの反復を通して、反復的に推定することを含む方法。
- 予め処理されたエコーの連続するフレームの各ペアの包絡線を圧縮し、
予め処理されたエコーの連続するフレームのペアの平均周波数を推定し、
予め処理されたエコーの連続するフレームのペア、遅れゼロ相互相関及び平均周波数に基づいて、最初の動き推定を推定すること
をさらに具備する請求項13記載の方法。 - 予め処理されたエコーの連続するフレームの各ペアは行列であり、
要素ごとに、前記予め処理されたエコーの連続するフレームのペアにおける要素を乗算し、第1の積行列を生成し、
所定の変換で第1の積行列を畳み込み、
前記畳み込み行列の位相シフトを決定し、
前記決定された位相シフトに基づいて、前記畳み込み行列の位相をアンラップし、
前記アンラップされた位相を前記最初の動作推定に変換すること
を含む請求項14記載の方法。 - 表面を最初の動き推定にフィッティングすることをさらに具備する請求項15記載の方法。
- 前記表面にフィットされた前記最初の動き推定を量子化し、非ゼロ相互相関遅れを決定することをさらに具備する請求項16記載の方法。
- 予め処理されたエコーの連続するフレームのペア及び前記決定された非ゼロ相互相関遅れに基づいて、後の動き推定を決定することをさらに具備する請求項17記載の方法。
- 前記非ゼロ相互相関遅れは行列であり、前記動き推定器は、
要素ごとに、前記予め処理されたエコーの連続するフレームのペアにおける要素及び前記非ゼロ相互相関遅れを乗算し、第2の積行列を生成し、
前記所定の変換で第2の積を畳み込み、
境界で位相推定のステップを除去し、デグリッチされた畳み込み行列を生成し、、
前記デグリッチされた畳み込み行列の位相シフトを決定し、
前記決定された位相シフトに基づいて前記デグリッチされた畳み込み行列の位相をアンラップし、
アンラップされた位相を第2の動き推定に変換することを具備する請求項18記載の動きプロセッサ。 - コード化されたコンピューター実行可能な命令を有するコンピューター可読記録媒体であって、前記命令は、プロセッサーによって実行された場合に、前記プロセッサーに、
予め処理されたエコーの連続する圧縮されたフレームのペア、遅れゼロ相互相関及び平均周波数に基づいて、最初の動き推定を決定させ、
前記最初の動き推定に基づいて非ゼロ相互相関遅れを決定させ、
前記予め処理されたエコーの圧縮された連続するフレームのペアにおける要素を非ゼロ相互相関遅れで要素毎に乗算させ、積行列を生成させ、
2D重み関数で前記積行列を畳み込みさせ、
境界での位相推定における非連続を除去させ、デグリッチされた畳み込み行列を生成させ、
前記デグリッチされた畳み込み行列の位相シフトを決定させ、
前記決定された位相シフトに基づき、前記デグリッチされた畳み込み行列の位相をアンラップさせ、
前記アンラップされた位相を前記動き推定に変換させる、
コンピューター可読記録媒体。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/IB2012/001854 WO2014045074A1 (en) | 2012-09-20 | 2012-09-20 | Imaging tissue motion estimation |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2015532848A true JP2015532848A (ja) | 2015-11-16 |
Family
ID=47073479
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2015532513A Pending JP2015532848A (ja) | 2012-09-20 | 2012-09-20 | 画像組織(tissue)動作推定 |
Country Status (4)
Country | Link |
---|---|
US (1) | US9964634B2 (ja) |
EP (1) | EP2898346B1 (ja) |
JP (1) | JP2015532848A (ja) |
WO (1) | WO2014045074A1 (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11832364B2 (en) | 2018-09-28 | 2023-11-28 | Apple Inc. | Systems and methods for wavelength locking in optical sensing systems |
US11835836B1 (en) | 2019-09-09 | 2023-12-05 | Apple Inc. | Mach-Zehnder interferometer device for wavelength locking |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9772219B2 (en) | 2012-09-11 | 2017-09-26 | S.P.M. Instrument Ab | Apparatus for monitoring the condition of a machine |
US10302741B2 (en) * | 2015-04-02 | 2019-05-28 | Texas Instruments Incorporated | Method and apparatus for live-object detection |
WO2018060820A1 (en) | 2016-09-29 | 2018-04-05 | Koninklijke Philips N.V. | Ultrasonic shear wave imaging with background motion compensation |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3462584B2 (ja) | 1994-02-14 | 2003-11-05 | フクダ電子株式会社 | 超音波診断装置 |
US6508768B1 (en) | 2000-11-22 | 2003-01-21 | University Of Kansas Medical Center | Ultrasonic elasticity imaging |
CN100475154C (zh) | 2002-07-31 | 2009-04-08 | 株式会社日立医药 | 超声诊断系统和应变分布显示方法 |
US7682311B2 (en) * | 2005-09-22 | 2010-03-23 | Siemens Medical Solutions Usa, Inc. | Phase unwrapped velocity display for ultrasound medical imaging |
-
2012
- 2012-09-20 US US14/429,639 patent/US9964634B2/en active Active
- 2012-09-20 EP EP12777945.2A patent/EP2898346B1/en active Active
- 2012-09-20 WO PCT/IB2012/001854 patent/WO2014045074A1/en active Application Filing
- 2012-09-20 JP JP2015532513A patent/JP2015532848A/ja active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11832364B2 (en) | 2018-09-28 | 2023-11-28 | Apple Inc. | Systems and methods for wavelength locking in optical sensing systems |
US11835836B1 (en) | 2019-09-09 | 2023-12-05 | Apple Inc. | Mach-Zehnder interferometer device for wavelength locking |
Also Published As
Publication number | Publication date |
---|---|
WO2014045074A1 (en) | 2014-03-27 |
EP2898346A1 (en) | 2015-07-29 |
US9964634B2 (en) | 2018-05-08 |
EP2898346B1 (en) | 2022-11-23 |
US20150234040A1 (en) | 2015-08-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR101140482B1 (ko) | 빔형성 시간 지연 수정 방법, 초음파 시스템 및 시간 지연수정 시스템 | |
US6760486B1 (en) | Flash artifact suppression in two-dimensional ultrasound imaging | |
US20130258805A1 (en) | Methods and systems for producing compounded ultrasound images | |
CN108135568B (zh) | 用于利用持续换能器振动进行超声弹性成像的系统和方法 | |
KR102041433B1 (ko) | 초음파 컬러 흐름에서의 스파클 아티팩트 검출 | |
JP2003010178A (ja) | 超音波診断装置 | |
KR101935514B1 (ko) | 탄성 이미징에서의 주파수 컴파운딩 | |
WO2012066568A1 (en) | An improved ultrasound imaging method/technique for speckle reduction/suppression in an improved ultra sound imaging system | |
KR20110053621A (ko) | 클러터 필터링을 적응적으로 수행하는 초음파 시스템 및 방법 | |
JP2015532848A (ja) | 画像組織(tissue)動作推定 | |
EP2830508A1 (en) | Methods and apparatus for ultrasound imaging | |
US20230012638A1 (en) | Adaptive clutter filtering in acoustic radiation force-based ultrasound imaging | |
US6775400B1 (en) | Medical diagnostic ultrasonic imaging method and apparatus for suppressing electronic noise | |
US10342509B2 (en) | Beamforming device, ultrasonic imaging device, and beamforming method allowing simple spatial smoothing operation | |
US10993701B2 (en) | Ultrasonic imaging device | |
CN110536647B (zh) | 用于在存在运动的情况下处理和复合超声图像的系统、方法和计算机可读介质 | |
US20170307725A1 (en) | Method for reducing sidelobe in ultrasound images | |
JP2007029703A (ja) | 弾性映像信号の非相関度を減少させて超音波映像を形成する方法 | |
US20220330920A1 (en) | Ultrasonic diagnostic apparatus and medical image processing apparatus | |
CN112842382B (zh) | 用于对信道数据进行流处理以应用非线性波束形成的方法和系统 | |
JP5634819B2 (ja) | 超音波診断装置 | |
KR100830152B1 (ko) | 탄성영상신호의 비상관도를 감소시켜 초음파 영상을형성하는 방법 | |
Srinivasan et al. | Theoretical derivation of SNR, CNR and spatial resolution for a local adaptive strain estimator for elastography | |
KR101610877B1 (ko) | 공간 일관성 기초 초음파 신호 처리 모듈 및 그에 의한 초음파 신호 처리 방법 |