JP6263380B2 - 剪断波中の背景ノイズを除去するための方法及び関連超音波イメージング・システム - Google Patents

剪断波中の背景ノイズを除去するための方法及び関連超音波イメージング・システム Download PDF

Info

Publication number
JP6263380B2
JP6263380B2 JP2013266224A JP2013266224A JP6263380B2 JP 6263380 B2 JP6263380 B2 JP 6263380B2 JP 2013266224 A JP2013266224 A JP 2013266224A JP 2013266224 A JP2013266224 A JP 2013266224A JP 6263380 B2 JP6263380 B2 JP 6263380B2
Authority
JP
Japan
Prior art keywords
group
positions
displacement information
imaging system
ultrasound imaging
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.)
Active
Application number
JP2013266224A
Other languages
English (en)
Other versions
JP2014128661A (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.)
GE Medical Systems Global Technology Co LLC
Original Assignee
GE Medical Systems Global Technology Co LLC
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 GE Medical Systems Global Technology Co LLC filed Critical GE Medical Systems Global Technology Co LLC
Publication of JP2014128661A publication Critical patent/JP2014128661A/ja
Application granted granted Critical
Publication of JP6263380B2 publication Critical patent/JP6263380B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5269Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
    • A61B8/5276Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts due to motion
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/13Tomography
    • A61B8/14Echo-tomography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/485Diagnostic techniques involving measuring strain or elastic properties
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52023Details of receivers
    • G01S7/52036Details of receivers using analysis of echo signal for target characterisation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52023Details of receivers
    • G01S7/52036Details of receivers using analysis of echo signal for target characterisation
    • G01S7/52042Details of receivers using analysis of echo signal for target characterisation determining elastic properties of the propagation medium or of the reflective target
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52077Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging with means for elimination of unwanted signals, e.g. noise or interference

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Surgery (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Veterinary Medicine (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Description

本発明は、一般的に云えば、超音波検出の分野に関し、より詳しく云えば、剪断波中の背景ノイズを除去するための方法及び装置、並びに該装置を有する超音波イメージング・システムに関するものである。
超音波検出は現在の医療分野において一般的な検出手段である。音波が不均質媒体中を伝播するとき、音響放射力が生成される。音響放射力は組織(tissue)の生体力学的弾性を特徴付けることができ、これは材料化学及び医学的診断の分野での研究において関心を引き起こしている。
長い期間にわたって、診断用イメージングの目的の1つは、組織の精密な特徴付けである。臨床医学者が望んでいることは、身体器官の診断領域を求めることであり、また診断システムを用いて、組織の特徴を画像中で識別することである。理想的には、臨床医学者が期待していることは、イメージング・システムにより、病理学的状態が悪性であるか又は良性であるかを識別することである。上記の目的を満たすために様々な努力が為されているが、診断用イメージングは組織の組成に関する手掛かりを臨床医学者に対して依然として示すことができる。当該分野に関する技術の1つは、身体内の組織の弾性又は硬さの程度を測定するのに使用するための弾性イメージングである。例えば、非常に硬い乳房腫瘍又はしこりは悪性の可能性があり、他方、比較的軟らかく且つ可撓性のしこりは良性である可能性がある。しこりの硬さの程度としこりの悪性度又は良性度との間の関連性により、弾性イメージングは、診断を行い且つ関連した治療の問題解決手法を考え出すための有用な別の証拠を、臨床医学者に対して提供する。
1つの問題解決手法は、剪断波を測定することである。身体中の或る1つの点における圧縮が解放されるとき、下側層の組織が下向きに圧縮され、次いで圧縮力の解放の際に上向きに跳ね返る。その際、圧縮力を受けている組織は周囲の組織に連続的に連結されているので、力ベクトルの側方向に位置する未圧縮の組織もまた、圧縮された組織の上下の動きに応答するようになる。側方向におけるこの凹み(dimple)効果はまた、剪断波と呼ばれており、これは、下向きの圧縮力に対する周囲組織によって生じる応答として作用する。更に、組織を下向きに押すのに必要な力は、超音波パルスからの放射圧力によって生成することができること、また更に超音波は、剪断波によって誘起された組織の動きを受けて、認識し及び測定するために使用することができることが判明している。剪断波速度が局部組織の機械的性質によって決定される。剪断波は軟組織を或る特定の速度で通過し、また硬組織をより高い速度で通過する。身体中の或る一点で剪断波の速度を測定することによって、組織の剪断弾性率、ヤング率及び動的剪断粘度のような、関連組織の性質に関する情報を得ることができる。側方向に伝播される剪断波は、容易に検出できるように、低速で、通常は毎秒数メートル又はそれよりも低い速度で伝播されるが、剪断波は数センチメートル以下の距離で急速に減衰する。各測定のために同じ「プッシュ・パルス(push puls) 」を繰り返すことができるので、剪断波技術は、超音波によって組織の性質を客観的に定量化するのに役立つ。更に、剪断波速度はプッシュ・パルス強度から独立しており、これによって、測定がユーザーによって左右されないようにする。
音響放射力イメージング(ARFI)は、放射力を小さな組織ボリュームへ集めるために集束超音波を用いている。次いで、通常の超音波を使用することにより、組織の剪断波変位を追跡し、超音波に基づく関連方法により組織の弾性を計算する。AFRI法による画像は、組織の機械的インパルス応答情報を提供する。Bモード(白黒モード)の場合、AFRI法は、特に肺野及び胸部に適用されたとき、コントラスト比を改善することができる。組織が病理学的変化を受けているとき、弾性は著しく変化することがある。多くの臨床研究により、組織の弾性についての情報が癌の診断及び治療に不可欠であることが実証されている。ARFI画像は、Bモード画像を補完するものであって、組織の機械的特性を特徴つけるための有用な方法である。
それにも拘わらず、ARFIはまた弱点もある。その主な弱点は、剪断波の振幅が非常に小さく、10μm程度に過ぎず、それによりAFRIがシステムの電子ノイズ及び患者の体動、特に心臓の動き又は呼吸運動、などのようなノイズの影響を受け易いことであり、このようなノイズは剪断波変位についての推定に影響を与える。従来技術における、上記の問題に対する解決策は、音波ビーム送出電圧又はその持続時間を増大することである。それにも拘わらず、該解決策は、米国食品医薬品局(FDA)によって課せられる制約を超える程度までも、音響出力を著しく増大させている。
本発明の目的は、上記の問題に対処することである。
前に述べた目的のために、本発明は、剪断波中の背景ノイズを除去するための方法を提供する。本方法は、プッシュ・パルス・ベクトルに沿ってプッシュ・パルスを送信する段階と、前記プッシュ・パルス・ベクトルから遠く離れた第1群の複数の位置で第1群の焦点追跡(focus track) パルスを送信する段階と、前記第1群の焦点追跡パルスに応答して第1群の焦点エコー信号を受信する段階と、前記第1群の焦点エコー信号を処理して、時間につれて変化する第1群の組織変位情報を決定する段階と、を有する。
本発明は更に、超音波イメージング・システムを提供し、該超音波イメージング・システムは、プッシュ・パルス・ベクトルに沿ってプッシュ・パルスを送信し、前記プッシュ・パルス・ベクトルから遠く離れた第1群の複数の位置で第1群の焦点追跡パルスを送信し、前記第1群の焦点追跡パルスに応答して第1群の焦点エコー信号を受信することのできる超音波プローブを有し、また該超音波イメージング・システムは、前記第1群の焦点エコー信号を処理して、時間につれて変化する第1群の組織変位情報を決定する。
本発明は、ARFI超音波システムに適用することができる。本発明によれば、剪断波中に含まれる背景ノイズを除くことができ、またプッシュ・パルスの持続時間及び電圧を低減することができる。
本発明のこれらの及び他の目的及び特徴は、添付の図面を参照して以下の詳しい説明を読めばより良く理解されよう。図面では、全体を通じて同様な参照符号は同じ又は同様な要素を表す。
図1Aは、本発明の原理を示す概略図である。 図1Bは、本発明の原理を示す概略図である。 図2は、本発明の一実施形態に従った、剪断波中の背景ノイズを除去するための方法を示す概略図である。 図3は、組織変位情報に従って剪断波速度値を計算する原理を例示する概略図である。
本書で用いられる選ばれた用語が定義され、それらには、それらの用語の範囲内に入り且つ発明の実現に使用することのできる構成部品の様々な例や形態が含まれる。それらの例は制限ではない。これらの定義は単数形及び複数形の両方の用語を含むことができる。
「一実施形態」、「様々な実施形態」、「一例」、「様々な例」などと云う場合、そのような実施形態(1つ又は複数)又は例(1つ又は複数)が、特別な特徴、構造、属性、性質又は制限を含み得ることを表す。とは云え、様々な実施形態又は様々な例の全てが、必ずしも特別な特徴、構造、属性、性質、要素又は制限を包含しているものではない。また、語句「一実施形態では」を繰り返し使用するが、これは、必ずしも同じ実施形態を表しているものではない。
説明を簡潔にするために、以下の記載では、当業者に周知の幾つかの技術的特徴は省略する。
図1A及び図1Bは、本発明の原理を図式的に示す。図1A及び図1Bにおいて、符号「10」は超音波プローブを表し、符号「11」は患者の身体又は模擬体を表し、符号「Push」はプッシュ・パルスを表し、「Track1」及び「Track2」は一群の追跡パルスを表し、及び「Track1’」及び「Track2’」は別の一群の追跡パルスを表す。
超音波プローブを用いることにより、プッシュ・パルス・ベクトルの方向に沿って組織について超音波押出しが実行されるように、組織に対して1つ以上のプッシュ・パルスを送信する。
従来のARFIイメージングでは、プッシュ・パルスはカラー・ドップラー・パルスと同様である。プッシュ・パルスは、組織中に剪断波を生成するように充分長くすべきである。プッシュ・パルスの後に一連の追跡パルスが続き、これらの追跡パルスは通常のBパルスと同様である。従来技術では、追跡パルスは常にプッシュ・パルスに隣接する。
本発明によれば、追跡パルス対が2つの異なる位置で供給される。一方の位置は、図1Aに「Track1」及び「Track2」によって例示されているような位置である。この場合、追跡パルス「Track1」及び「Track2」は、プッシュ・パルス「Push」から遠く離れた位置に配置されて、組織からのノイズを捕獲するようにする。組織中の剪断波は急速に減衰し、従って剪断波は「Track1」及び「Track2」の位置では完全に減衰している、すなわち、その位置ではノイズ以外には何ら剪断波はない。
他方の位置は、従来技術における位置と同じであり、図1Bに「Track1’」及び「Track2’」で例示されているような、ノイズを含む剪断波を収集する位置である。
最後に、第1の位置における信号が第2の位置における信号から減算されて、ノイズを含まない剪断波が求められる。このようにして、組織中のノイズの影響を除去することができる。
図2は、本発明の一実施形態に従った、剪断波中の背景ノイズを除去するための方法200を図式的に示す。
本発明の一実施形態によれば、方法200は段階210〜240を含むことができる。
段階210において、プッシュ・パルス・ベクトルに沿ってプッシュ・パルスを送信する。本発明の好ましい一実施形態によれば、先ずプッシュ・パルスの比較的長いビームを1つのフレームで送信する。パルスによる励起時に、水平に伝播される剪断波が焦点組織へ生じる。
段階220において、プッシュ・パルス・ベクトルから遠く離れた第1群の複数の位置で第1群の焦点追跡パルスを送信する。第1群の焦点追跡パルスは、図1Aの「Track1」及び「Track2」に対応することができる。本発明の好ましい一実施形態によれば、複数の追跡パルス「Track1」及び「Track2」は、2つの音波ビーム「Track1」及び「Track2」の位置において交互に送信することができ、この場合、「Track1」及び「Track2」を送信するための時間切換えが既定で既知であり、また「Track1」及び「Track2」の間の物理的距離が既定で既知である。
段階230において、第1群の焦点追跡パルスに応答して第1群の焦点エコー信号を受信する。例えば、各追跡パルス・ベクトルは時間的にインターリーブする態様で繰り返しサンプリングすることができ、従って、各追跡パルス・ベクトル位置が剪断波の結果としてシフトしたとき、そのシフトは、好ましくは、該ベクトルから得られた相次いで照会されたエコー・データを相関させることによって検出することができる。剪断波がプッシュ・パルス・ベクトルから離れる側方向に沿って動くとき、追跡パルスの位置決めもまた、剪断波の伝播を追跡するように側方向に沿ってシフトさせることができる。
段階240において、第1群の焦点エコー信号を処理して、時間につれて変化する第1群の組織変位情報を決定する。第1群の焦点追跡パルスがプッシュ・パルスから遠く離れて位置決めされているので、剪断波は完全に減衰している。第1群の組織変位情報は、組織中の背景ノイズによって引き起こされた変位を含むが、剪断波の伝播によって引き起こされた変位を含まない。
本発明の別の実施形態によれば、方法200はまた、段階250〜280を含むことができる。
段階250において、第2群の焦点追跡パルスがプッシュ・パルス・ベクトルに隣接した第2群の複数の位置で送信される。第2群の焦点追跡パルスは、図1Bの「Track1’」及び「Track2’」に対応することができる。本発明の好ましい一実施形態によれば、複数の追跡パルス「Track1’」及び「Track2’」は、2つの音波ビーム「Track1’」及び「Track2’」の位置において交互に送信することができ、この場合、「Track1’」及び「Track2’」を送信するための時間切換えが既定で既知であり、また「Track1’」及び「Track2’」の間の物理的距離が既定で既知である。
段階260において、第2群の焦点追跡パルスに応答して第2群の焦点エコー信号を受信する。
段階270において、第2群の焦点エコー信号を処理して、時間につれて変化する第2群の組織変位情報を決定する。第2群の組織変位情報は、剪断波の伝播によって引き起こされた変位と組織中の背景ノイズによって引き起こされた変位との両方を含む。背景ノイズは、例えば、患者の心臓の動き又は呼吸運動に由来する。
段階280において、第1群の組織変位情報を第2群の組織変位情報から減算して、減算後の一群の組織変位情報を求める。減算後の一群の組織変位情報は、剪断波によって引き起こされた組織変位のみを含み、これによって剪断波信号中の背景ノイズが除去されている。
本発明の一実施形態によれば、第1群の複数の位置の中の位置の数は、第2群の複数の位置の中の位置の数に等しい。
本発明の一実施形態によれば、図1A及び図1Bに例示されているように、第1群の複数の位置の中の位置の数及び第2群の複数の位置の中の位置の数はそれぞれ2である。
本発明の一実施形態によれば、第1群の複数の位置の中の位置の数は第2群の複数の位置の中の位置の数に等しく、その各々の数は3つ以上である。この実施形態では、各群の位置の中の相異なる位置に対応する組織変位情報は、先ず、各群の組織変位情報の中に含まれる組織変位情報の数が2になるように、平均化することができる。例えば、第1群の位置は、追跡パルス「Track1」、「Track2」及び「Track3」に対応する3つの位置より成り、また第2群の位置は、追跡パルス「Track1’」、「Track2’」及び「Track3’」に対応する3つの位置より成る。次いで、「Track1」及び「Track2」に対応する組織変位情報を平均化して、「Track12」を得ることができ、また「Track2」及び「Track3」に対応する組織変位情報を平均化して、「Track23」を得ることができる。このようにして、ノイズのみを含む剪断波及びノイズの両方を含む一群の組織変位情報「Track12」及び「Track23」を得ることができる。その後、「Track1’」及び「Track2’」に対応する組織変位情報を平均化して、「Track1’2’」を得ることができ、また「Track2’」及び「Track3’」に対応する組織変位情報を平均化して、「Track2’3’」を得ることができる。このようにして、剪断波及びノイズの両方を含む一群の組織変位情報「Track1’2’」及び「Track2’3’」を得ることができる。最後に、前者の一群の組織変位情報「Track12」及び「Track23」を後者の一群の組織変位情報「Track1’2’」及び「Track2’3’」から減算し、これによって、剪断波のみを含む一群の組織変位情報が得られる。ここで、関連した平均化方法は上記方法に限定されず、当該分野で知られている又は将来開発される任意の平均化方法とすることができることに留意されたい。
本発明の一実施形態によれば、第1群の複数の位置の中の位置の数は、第2群の複数の位置の中の位置の数に等しくない。この実施形態では、相対的に数の多い方の位置に関する一群の組織変位情報を先ず平均化して、これにより2つの群の組織変位情報の数を互いに等しくすることができる。
本発明の一実施形態によれば、減算後の一群の組織変位情報に従って剪断波速度値を計算することができる。例えば、繰り返しサンプリングされた追跡パルス・ベクトルからのデータは、剪断波が追跡パルス・ベクトルの各点におけるピーク値変位をトリガする時間を見い出すように処理することができ、これは、好ましくは、曲線当てはめ又は補間を介して、相次ぐ変位測定結果を相間させることによって実現される。隣接するサンプリングされたベクトルの中の様々な点においてピーク値変位が生じる時間を分析することにより、特定のベクトル位置に対応する剪断波の速度の測定結果が得られよう。
図3について説明すると、図3は、組織変位情報に従って剪断波速度値を計算する原理を例示する概略図である。例えば、「Track1」の位置において時間につれて変化する組織の変位を、後処理アルゴリズムにより算定することができ、この変位は剪断波の伝播により生じる。同様にして、「Track2」の位置において時間につれて変化する組織の変位を算定することができる。「Track1」の位置における変位の最大値d1(max) と「Track2」の位置における変位の最大値d2(max) との間の時間差Tは、「Track1」の位置から「Track2」の位置への剪断波の伝播のために必要とされる期間と見なすことができ、また更に、関心のある組織中の剪断波の伝播速度を算定することができる。具体的に述べると、
剪断波速度=(d1(max) −d2(max) )/T
である。
本発明の一実施形態によれば、ヤング率を剪断波速度値に従って計算することができる。別の方法で述べると、速度変化は組織中の硬さ又は弾性の異なる程度を示す。剪断波の伝播速度と組織のヤング率との間に下記のような関係があるので、ヤング率を導き出すことができる。
剪断波速度=[E/(2(1+ν)ρ)]1/2
ここで、Eはヤング率であり、ρは組織の密度であり、またνはポアソン比である。
本発明の一実施形態によれば、背景ノイズは、患者の心臓の動き又は呼吸運動により生じる。
剪断波は非常に急速に減衰するので、一般に、単一のプッシュ・パルス・ベクトルを用いることによって画像フィールド全体の剪断波データを収集することは不可能である。従って、組織の別の領域における剪断波速度測定結果を収集するように、組織の別の位置において上記のプロセスを繰り返すことが必要である。該プロセスは、期待される画像フィールド全体における剪断波の収集が完了するまで繰り返される。
本発明はまた、下記の超音波プローブを有する超音波イメージング・システムを提供する。
超音波プローブは、プッシュ・パルス・ベクトルに沿ってプッシュ・パルスを送信することができる。本発明の好ましい一実施形態によれば、先ずプッシュ・パルスの比較的長いビームを1つのフレームで送信する。パルスによる励起時に、水平に伝播される剪断波が焦点組織へ生じる。
超音波プローブはまた、プッシュ・パルス・ベクトルから遠く離れた第1群の複数の位置で第1群の焦点追跡パルスを送信することができる。第1群の焦点追跡パルスは、図1Aの「Track1」及び「Track2」に対応することができる。本発明の好ましい一実施形態によれば、複数の追跡パルス「Track1」及び「Track2」は、2つの音波ビーム「Track1」及び「Track2」の位置において交互に送信することができ、この場合、「Track1」及び「Track2」を送信するための時間切換えが既定で既知であり、また「Track1」及び「Track2」の間の物理的距離が既定で既知である。
超音波プローブはまた、第1群の焦点追跡パルスに応答して第1群の焦点エコー信号を受信することができる。例えば、各追跡パルス・ベクトルは時間的にインターリーブする態様で繰り返しサンプリングすることができ、従って、各追跡パルス・ベクトル位置が剪断波の結果としてシフトしたとき、そのシフトは、好ましくは、該ベクトルから得られた相次いで照会されたエコー・データを相関させることによって検出することができる。剪断波がプッシュ・パルス・ベクトルから離れる側方向に沿って動くとき、追跡パルスの位置決めもまた、剪断波の伝播を追跡するように側方向に沿ってシフトさせることができる。
本発明の超音波イメージング・システムは、第1群の焦点エコー信号を処理して、時間につれて変化する第1群の組織変位情報を決定することができる。第1群の焦点追跡パルスがプッシュ・パルスから遠く離れて位置決めされているので、剪断波は減衰している。第1群の組織変位情報は、組織中の背景ノイズによって引き起こされた変位を含むが、剪断波の伝播によって引き起こされた変位を含まない。
超音波プローブはまた、第2群の焦点追跡パルスがプッシュ・パルス・ベクトルに隣接した第2群の複数の位置で第2群の焦点追跡パルスを送信することができる。第2群の焦点追跡パルスは、図1Bの「Track1’」及び「Track2’」に対応することができる。本発明の好ましい一実施形態によれば、複数の追跡パルス「Track1’」及び「Track2’」は、2つの音波ビーム「Track1’」及び「Track2’」の位置において交互に送信することができ、この場合、「Track1’」及び「Track2’」を送信するための時間切換えが既定で既知であり、また「Track1’」及び「Track2’」の間の物理的距離が既定で既知である。
超音波プローブはまた、第2群の焦点追跡パルスに応答して第2群の焦点エコー信号を受信することができる。
本発明の超音波イメージング・システムはまた、第2群の焦点エコー信号を処理して、時間につれて変化する第2群の組織変位情報を決定することができる。第2群の組織変位情報は、剪断波の伝播によって引き起こされた変位と組織中の背景ノイズによって引き起こされた変位との両方を含む。背景ノイズは、例えば、患者の心臓の動き又は呼吸運動に由来する。
本発明の超音波イメージング・システムはまた、第1群の組織変位情報を第2群の組織変位情報から減算して、減算後の一群の組織変位情報を求める。減算後の一群の組織変位情報は、剪断波によって引き起こされた組織変位のみを含み、これによって剪断波信号中の背景ノイズを除去する。
本発明の一実施形態によれば、第1群の複数の位置の中の位置の数は、第2群の複数の位置の中の位置の数に等しい。
本発明の一実施形態によれば、図1A及び図1Bに例示されているように、第1群の複数の位置の中の位置の数及び第2群の複数の位置の中の位置の数は、それぞれ2である。
本発明の一実施形態によれば、第1群の複数の位置の中の位置の数及び第2群の複数の位置の中の位置の数は、それぞれ3つ以上である。この実施形態では、各群の位置の中の相異なる位置に対応する組織変位情報は、先ず、各群の組織変位情報の中に含まれる組織変位情報の数が2になるように、平均化することができる。例えば、第1群の位置は、追跡パルス「Track1」、「Track2」及び「Track3」に対応する3つの位置より成り、また第2群の位置は、追跡パルス「Track1’」、「Track2’」及び「Track3’」に対応する3つの位置より成る。次いで、「Track1」及び「Track2」に対応する組織変位情報を平均化して、「Track12」を得ることができ、また「Track2」及び「Track3」に対応する組織変位情報を平均化して、「Track23」を得ることができる。このようにして、ノイズのみを含む一群の組織変位情報「Track12」及び「Track23」を得ることができる。その後、「Track1’」及び「Track2’」に対応する組織変位情報を平均化して、「Track1’2’」を得ることができ、また「Track2’」及び「Track3’」に対応する組織変位情報を平均化して、「Track2’3’」を得ることができる。。このようにして、剪断波及びノイズの両方を含む一群の組織変位情報「Track1’2’」及び「Track2’3’」を得ることができる。最後に、前者の一群の組織変位情報「Track12」及び「Track23」を後者の一群の組織変位情報「Track1’2’」及び「Track2’3’」から減算し、これによって、剪断波のみを含む一群の組織変位情報が得られる。ここで、関連した平均化方法は上記方法に限定されず、当該分野で知られている又は将来開発される任意の平均化方法とすることができることに留意されたい。
本発明の一実施形態によれば、第1群の複数の位置の中の位置の数は、第2群の複数の位置の中の位置の数に等しくない。この実施形態では、相対的に数の多い方の位置に関する一群の組織変位情報を先ず平均化して、これにより2つの群の組織変位情報の数を互いに等しくすることができる。
本発明の一実施形態によれば、超音波イメージング・システムはまた、減算後の一群の組織変位情報に従って剪断波速度値を計算することができる。例えば、繰り返しサンプリングされた追跡パルス・ベクトルからのデータは、剪断波が追跡パルス・ベクトルの各点におけるピーク値変位をトリガする時間を見い出すように処理することができ、これは、好ましくは、曲線当てはめ又は補間を介して、相次ぐ変位測定結果を相間させることによって実現される。隣接するサンプリングされたベクトルの中の様々な点においてピーク値変位が生じる時間を分析することにより、特定のベクトル位置に対応する剪断波の速度の測定結果が得られよう。
本発明の一実施形態によれば、超音波イメージング・システムはまた、ヤング率を剪断波速度値に従って計算することができる。別の方法で述べると、速度変化は組織中の硬さ又は弾性の異なる程度を示す。剪断波の伝播速度と組織のヤング率との間に下記のような関係があるので、ヤング率を導き出すことができる。
剪断波速度=[E/(2(1+ν)ρ)]1/2
ここで、Eはヤング率であり、ρは組織の密度であり、またνはポアソン比である。
本発明の一実施形態によれば、背景ノイズは、例えば、患者の心臓の動き又は呼吸運動に由来する。
本発明の技術的解決策はARFI超音波システムに適用することができる。本発明によれば、剪断波中に含まれている背景ノイズを除去することができる。従って、本発明は、音波ビームの持続時間又は送出電圧を増大する必要がないと云う点で従来技術とは異なる。
本発明によれば、プッシュ・パルスの持続時間を短くすることができ、これによって音波出力を著しく低減することができる。従って、本発明は、FDAによって規定された制約内に音波出力を抑えるばかりでなく、患者から見ても一層信頼性がある。
本発明は、後処理に採用されていた複雑なフィルタよりも良好な効果を持つ。
上述した様々な実施形態は模範的なものであり且つ制限するものではないことに留意されたい。また更に、当業者には、特許請求の範囲から逸脱せずに、より多くの様々な代替実施形態を考案することができよう。本書で用いられる用語「有する」は、「特許請求の範囲」及び明細書に列挙したもの以外の他の構成要素又は処理段階を排除することを意味していない。更に、本書に含まれている構成要素の数は、1つ又は複数とすることができる。
10 超音波プローブ
11 患者身体又は模擬体
200 方法

Claims (16)

  1. 剪断波中の背景ノイズを除去するように超音波イメージング・システムを動作させるための方法であって、
    プッシュ・パルス・ベクトルに沿ってプッシュ・パルスを送信する段階と、
    前記背景ノイズを捕獲するために、前記プッシュ・パルス・ベクトルから遠く離れて同じ方向において異なる距離を有する第1群の複数の位置で第1群の焦点追跡パルスを送信する段階であって、前記焦点追跡パルスは前記同じ方向において前記異なる距離に沿って交互に順番に送信される段階と、
    前記第1群の焦点追跡パルスに応答して第1群のインターリーブされた焦点エコー信号を受信する段階と、
    前記第1群の焦点エコー信号を処理して、時間につれて変化する第1群の組織変位情報を決定する段階と、
    前記プッシュ・パルス・ベクトルに隣接した第2群の複数の位置で第2群の焦点追跡パルスを交互に順番に送信する段階と、
    前記第2群の焦点追跡パルスに応答して第2群の焦点エコー信号を受信する段階と、
    前記第2群の焦点エコー信号を処理して、時間につれて変化する第2群の組織変位情報を決定する段階と、
    前記第1群の組織変位情報を前記第2群の組織変位情報から減算して、前記背景ノイズを除外した減算後の一群の組織変位情報を求める段階と、
    を有する方法。
  2. 前記第1群の複数の位置の中の位置の数は、前記第2群の複数の位置の中の位置の数に等しい、請求項記載の方法。
  3. 前記第1群の複数の位置の中の位置の数及び前記第2群の複数の位置の中の位置の数はそれぞれ2である、請求項記載の方法。
  4. 前記第1群の複数の位置の中の位置の数及び前記第2群の複数の位置の中の位置の数は、それぞれ3つ以上である、請求項記載の方法。
  5. 前記第1群の複数の位置の中の位置の数は、前記第2群の複数の位置の中の位置の数に等しくない、請求項記載の方法。
  6. 更に、前記減算後の一群の組織変位情報に従って剪断波速度値を計算する段階を有する請求項記載の方法。
  7. 更に、前記剪断波速度値に従ってヤング率を計算する段階を有する請求項記載の方法。
  8. 前記背景ノイズが患者の心臓の動き又は呼吸運動に由来する、請求項1記載の方法。
  9. プッシュ・パルス・ベクトルに沿ってプッシュ・パルスを送信し、前記背景ノイズを捕獲するために、前記プッシュ・パルス・ベクトルから遠く離れて同じ方向において異なる距離を有する第1群の複数の位置で、前記同じ方向において前記異なる距離に沿って交互に順番に送信されるように第1群の焦点追跡パルスを送信し、前記第1群の焦点追跡パルスに応答して第1群のインターリーブされた焦点エコー信号を受信することのできる超音波プローブを有し、
    前記第1群の焦点エコー信号を処理して、時間につれて変化する第1群の組織変位情報を決定することができる超音波イメージング・システムであって、
    前記超音波プローブが更に、前記プッシュ・パルス・ベクトルに隣接した第2群の複数の位置で第2群の焦点追跡パルスを交互に順番に送信し、前記第2群の焦点追跡パルスに応答して第2群の焦点エコー信号を受信することができ、
    前記超音波イメージング・システムが更に、前記第2群の焦点エコー信号を処理して、時間につれて変化する第2群の組織変位情報を決定し、次いで前記第1群の組織変位情報を前記第2群の組織変位情報から減算して、前記背景ノイズを除外した減算後の一群の組織変位情報を求めることができる超音波イメージング・システム。
  10. 前記第1群の複数の位置の中の位置の数は、前記第2群の複数の位置の中の位置の数に等しい、請求項記載の超音波イメージング・システム。
  11. 前記第1群の複数の位置の中の位置の数及び前記第2群の複数の位置の中の位置の数はそれぞれ2である、請求項1記載の超音波イメージング・システム。
  12. 前記第1群の複数の位置の中の位置の数及び前記第2群の複数の位置の中の位置の数は、それぞれ3つ以上である、請求項1記載の超音波イメージング・システム。
  13. 前記第1群の複数の位置の中の位置の数は、前記第2群の複数の位置の中の位置の数に等しくない、請求項記載の超音波イメージング・システム。
  14. 前記超音波イメージング・システムは更に、前記減算後の一群の組織変位情報に従って剪断波速度値を計算することができる、請求項記載の超音波イメージング・システム。
  15. 前記超音波イメージング・システムは更に、前記剪断波速度値に従ってヤング率を計算することができる、請求項14記載の超音波イメージング・システム。
  16. 前記背景ノイズが患者の心臓の動き又は呼吸運動に由来する、請求項記載の超音波イメージング・システム。
JP2013266224A 2012-12-31 2013-12-25 剪断波中の背景ノイズを除去するための方法及び関連超音波イメージング・システム Active JP6263380B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201210593023.9 2012-12-31
CN201210593023.9A CN103908289B (zh) 2012-12-31 2012-12-31 消除剪切波中的背景噪声的方法和相应的超声成像系统

Publications (2)

Publication Number Publication Date
JP2014128661A JP2014128661A (ja) 2014-07-10
JP6263380B2 true JP6263380B2 (ja) 2018-01-17

Family

ID=51017974

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2013266224A Active JP6263380B2 (ja) 2012-12-31 2013-12-25 剪断波中の背景ノイズを除去するための方法及び関連超音波イメージング・システム

Country Status (3)

Country Link
US (1) US10143446B2 (ja)
JP (1) JP6263380B2 (ja)
CN (1) CN103908289B (ja)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6288996B2 (ja) * 2013-09-11 2018-03-07 キヤノンメディカルシステムズ株式会社 超音波診断装置及び超音波イメージングプログラム
CN106572838B (zh) * 2014-09-03 2019-09-06 深圳迈瑞生物医疗电子股份有限公司 弹性测量检测方法及系统
CN104605890B (zh) * 2014-12-18 2017-04-05 深圳开立生物医疗科技股份有限公司 一种剪切波峰值波形校正的方法、装置、系统及其应用
CN110149286A (zh) * 2015-05-11 2019-08-20 蒋华丽 加性噪声环境下数字通信的提高信噪比的信号处理方法
EP3419516B1 (en) * 2016-02-23 2020-12-16 Mayo Foundation for Medical Education and Research Ultrasound blood flow imaging
US10675007B2 (en) * 2016-04-19 2020-06-09 Siemens Medical Solutions Usa, Inc. Frequency compounding in elasticity imaging
KR102608323B1 (ko) * 2016-04-21 2023-11-30 듀크 유니버시티 보강적인 전단파 초음파 이미징을 위한 방법, 시스템 및 컴퓨터 프로그램 제품
US11138723B2 (en) * 2016-07-22 2021-10-05 Canon Medical Systems Corporation Analyzing apparatus and analyzing method
CN108158610B (zh) * 2018-01-16 2024-02-06 苏州国科昂卓医疗科技有限公司 一种弹性成像方法、装置、设备及超声成像探头
WO2019210479A1 (zh) * 2018-05-03 2019-11-07 深圳迈瑞生物医疗电子股份有限公司 剪切波弹性成像方法及超声成像设备
CN110292399B (zh) * 2018-05-04 2022-03-08 深圳迈瑞生物医疗电子股份有限公司 一种剪切波弹性测量的方法及系统

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101431942B (zh) * 2006-03-20 2013-08-28 松下电器产业株式会社 超声波诊断装置
JP5437820B2 (ja) * 2008-02-21 2014-03-12 株式会社日立メディコ 超音波診断装置、超音波画像処理方法
US9364194B2 (en) * 2008-09-18 2016-06-14 General Electric Company Systems and methods for detecting regions of altered stiffness
US20100286520A1 (en) * 2009-05-11 2010-11-11 General Electric Company Ultrasound system and method to determine mechanical properties of a target region
US9125547B2 (en) * 2009-06-30 2015-09-08 Koninklijke Philips N.V. Push/tracking sequences for shear wave dispersion vibrometry
US8727995B2 (en) * 2010-09-09 2014-05-20 Siemens Medical Solutions Usa, Inc. Reduction of motion artifacts in ultrasound imaging with a flexible ultrasound transducer
RU2603051C2 (ru) * 2010-12-13 2016-11-20 Конинклейке Филипс Электроникс Н.В. Коррекция результатов измерений воздействия силы акустического излучения с учетом эффектов фонового движения

Also Published As

Publication number Publication date
CN103908289B (zh) 2019-11-12
JP2014128661A (ja) 2014-07-10
US20140187951A1 (en) 2014-07-03
US10143446B2 (en) 2018-12-04
CN103908289A (zh) 2014-07-09

Similar Documents

Publication Publication Date Title
JP6263380B2 (ja) 剪断波中の背景ノイズを除去するための方法及び関連超音波イメージング・システム
Feigin et al. A deep learning framework for single-sided sound speed inversion in medical ultrasound
EP2654552B1 (en) Shear wave velocity estimation using center of mass
JP4890554B2 (ja) 超音波診断装置
Song et al. Quantitative assessment of left ventricular diastolic stiffness using cardiac shear wave elastography: A pilot study
CN106999162B (zh) 超声波诊断装置以及弹性评价方法
US10338203B2 (en) Classification preprocessing in medical ultrasound shear wave imaging
JP5509437B2 (ja) 超音波診断装置
US8864671B2 (en) Methods and systems for color flow imaging
US20120108968A1 (en) Tissue Density Quantification Using Shear Wave Information in Medical Ultrasound Scanning
US20150080730A1 (en) Ultrasonic diagnostic apparatus, medical image processing apparatus, and medical image processing method
WO2006082966A1 (ja) 超音波診断装置
JPWO2014162966A1 (ja) 超音波診断装置、及び弾性評価方法
US20150094579A1 (en) Shear Wave Detection in Medical Ultrasound Imaging
CN103190932B (zh) 一种冠状动脉血管壁应力和应变的估算方法
JP2007222533A (ja) 超音波診断装置及び超音波画像処理方法
WO2019087741A1 (ja) 超音波診断装置、および、生体組織の物性評価方法
JP2008161674A (ja) 超音波診断装置
JP2009039277A (ja) 超音波診断装置
JP5844175B2 (ja) 超音波診断装置および超音波画像生成方法
SE531787C2 (sv) Global och lokal detektion av blodkärlselasticitet
EP2853918B1 (en) Shear wave detection in medical ultrasound imaging
Correia et al. Non-invasive myocardial shear wave elastography device for clinical applications in cardiology
Joseph et al. Non-invasive estimation of arterial compliance
WO2007072720A1 (ja) 医用画像診断装置および生体組織の同定方法

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20161219

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20170905

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20170913

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20171027

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: 20171205

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20171218

R150 Certificate of patent or registration of utility model

Ref document number: 6263380

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250