JP6259953B2 - せん断波を用いる測定に関する超音波システム及び方法 - Google Patents

せん断波を用いる測定に関する超音波システム及び方法 Download PDF

Info

Publication number
JP6259953B2
JP6259953B2 JP2017522080A JP2017522080A JP6259953B2 JP 6259953 B2 JP6259953 B2 JP 6259953B2 JP 2017522080 A JP2017522080 A JP 2017522080A JP 2017522080 A JP2017522080 A JP 2017522080A JP 6259953 B2 JP6259953 B2 JP 6259953B2
Authority
JP
Japan
Prior art keywords
parameter
interest
region
characteristic
shear wave
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.)
Expired - Fee Related
Application number
JP2017522080A
Other languages
English (en)
Other versions
JP2017535327A (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.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips NV
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 Koninklijke Philips NV filed Critical Koninklijke Philips NV
Publication of JP2017535327A publication Critical patent/JP2017535327A/ja
Application granted granted Critical
Publication of JP6259953B2 publication Critical patent/JP6259953B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0833Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures
    • A61B8/085Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures for locating body or organic structures, e.g. tumours, calculi, blood vessels, nodules
    • 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
    • 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/5215Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5223Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
    • 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
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8909Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration
    • G01S15/8915Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array
    • G01S15/8927Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array using simultaneously or sequentially two or more subarrays or subapertures
    • 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/52019Details of transmitters
    • G01S7/5202Details of transmitters for pulse systems
    • G01S7/52022Details of transmitters for pulse systems using a sequence of pulses, at least one pulse manipulating the transmissivity or reflexivity of the medium
    • 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/52079Constructional features
    • G01S7/5208Constructional features with integration of processing functions inside probe or scanhead
    • 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/52085Details related to the ultrasound signal acquisition, e.g. scan sequences
    • G01S7/52095Details related to the ultrasound signal acquisition, e.g. scan sequences using multiline receive beamforming
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment

Landscapes

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

Description

本発明は、せん断波を用いた測定、特に、対象の関心領域の特性をせん断波を用いて測定する超音波システム及び測定方法に関する。
生体組織の機械的変化は病理学的変化と相関する。例えば、組織の粘性、組織の剛性(弾性としても知られる)、及び組織の減衰係数(組織の縦波の減衰、及びせん断波の減衰としても知られる組織の横波の減衰を含む)は、臨床診療にとって重要な物理的パラメータである。組織の機械的特性を遠隔から調べるための様々な手段が開発されている。これは、超音波ビームの放射力を利用して、患者といった対象の体内の組織の領域に遠隔から力(音響放射力。これは、プッシュパルス又はプッシュビームとも呼ばれる)を加えて組織を変形させる。変形のパターンに従う縦波ベースの超音波イメージングを使用して直接変形を追跡することにより、変形のポイント(焦点と呼ばれる)で局所的に、又は変形される領域(即ち、焦点)から横方向に伝搬するせん断波をせん断波速度撮像を介して追跡することにより、焦点の隣接領域で、弾性特性が測定されることができるよう、音響放射力が適用されることができる。
医用イメージングの目的で、超音波による探査は、しばしば縦波を利用する。縦波は、伝搬方向における前後の運動により特徴付けられる。縦波に基づかれる従来の組織減衰測定では、組織縦波減衰(組織縦減衰とも呼ばれる)が、超音波エコー信号に基づき推定される。このエコー信号は、組織の後方散乱及び減衰の両方により同時に影響を受ける。斯かる従来の組織縦波減衰測定は、後方散乱の影響と組織減衰の影響とを分離することが困難であるため、精度が制限されるという欠点を持つ。
対照的に、超音波せん断(又は横波)波は、伝搬方向に垂直な前後運動により特徴付けられる。今日、多くの市販の超音波スキャナは、組織せん断弾性を測定するための超音波せん断波エラストグラフィ製品を提供する。せん断波エラストグラフィによる組織せん断粘性とせん断波減衰の推定はまだ実用化されていない。それらは、科学的に活発な研究課題のままである。なぜなら、臨床応用可能性が、特定の用途で浮上しているためである。図1は、従来技術によるせん断波を用いた測定を示す。プッシュパルス(プッシュビームとしても知られる)110は、焦点130に向かって送信され、せん断波150を生成する。これは、プッシュパルスの伝搬方向(即ち、長手方向z)に垂直な横方向xといった方向において焦点から伝搬する。1つ又は複数の追跡パルス(追跡ビームとしても知られる)が送信され、超音波エコー信号が、横方向に沿って間隔を置いて配置された複数のサンプリング位置140、142、144のそれぞれにおいて、せん断波の位相又は伝搬時間を推定するため、複数の追跡ライン(「Aライン」と呼ばれる)120、122、124に沿って受信される。複数のサンプリング位置におけるせん断波の推定された位相又は伝搬時間が、せん断波の速度値を得るために更に使用される。せん断波の得られた速度値は、その処理が放射力インパルス撮像として知られる超音波画像を生成するために、及び/又は組織粘性若しくは弾性といった機械的特性を得るために使用され得る。斯かる従来のせん断波エラストグラフィ技術は、組織せん断弾性、せん断粘度及びせん断波減衰を提供することができる。
US2012/0089019A1及びUS2011/0263978A1号は共に、従来のせん断波エラストグラフィ技術、即ち、せん断波の伝搬に対するその影響に基づき、組織の機械的特性を推定することに向けられている。
従って、対象における関心領域の特性(例えば、組織機械的特性)を測定するための改良された超音波システム及び方法を提供することが有利である。
本発明の第1の側面の実施形態によれば、超音波撮像システムが提案される。この超音波システムは、それぞれが相互に異なる深度値を持つ上記関心領域における複数の焦点の各々に、せん断波を発生させるためのプッシュパルスを順次送信し、上記複数の焦点のそれぞれに隣接する超音波エコー信号を受信するよう構成される超音波プローブと、上記受信された超音波エコー信号に基づき、上記生成されたせん断波が上記焦点で持つ特性を示す第1のパラメータを、複数の焦点の各々対して得るよう構成されるせん断波検出器と、上記得られた第1のパラメータの関数として、上記関心領域の特性を示す第2のパラメータを推定するよう構成される特性推定器とを有する。
こうして、関心領域内の異なる深度で複数のせん断波が励起され、複数のせん断波の強度に基づき、関心領域の特性が推定される。本発明の発明者らは、関心領域の特性は、プッシュパルスが関心領域を通って伝搬するとき、各プッシュパルスの伝搬に影響を及ぼし、次に、関心領域における焦点においてプッシュパルスによって励起されたそれぞれのせん断波に影響を及ぼし、及び従って焦点における励起されたせん断波の特性(強度、位相など)は、関心領域の特性(組織粘性、組織剛性、組織減衰など)に依存することを認識した。例えば、組織の長手方向の減衰係数が高いほど、プッシュパルスが組織を通って伝搬するときプッシュパルスがより多く減衰する態様で、組織は、プッシュパルスの伝搬を減衰させ、及び従って、焦点におけるプッシュパルスの強度、ひいては励起されたせん断波の強度が低下する。斯かる認識に基づき、本発明者らは、関心領域における異なる焦点深度で複数のせん断波を励起させ、長手方向に沿ったプッシュパルスの伝搬に対する関心領域の特性の影響に基づき、複数の焦点で生成されるせん断波の特性の関数として、関心領域の特性を推定することを提案する。言い換えると、関心領域の特性は、複数の焦点で生成されたせん断波の特性間の関係に基づき推定される。例えば、ある焦点に関して得られた第1のパラメータは、生成されたせん断波が焦点で持つ特性(強度、位相など)を示すことができる。これとは対照的に、従来のせん断波エラストグラフィ技術では、横波/横方向に沿ったせん断波の伝搬に関する関心領域の特性の影響に基づき、関心領域の特性が推定される。従って、従来のせん断波エラストグラフィ技術では、横方向に沿って伝搬する前に焦点でせん断波が持つ任意の特性が推定又は得られるのではなく、せん断波の伝搬速度又は遅延といったせん断波の伝搬特性が、推定される。
せん断波は、超音波信号の送信/受信方向(即ち、長手方向)に垂直な方向に伝搬するので、せん断波の測定は、長手方向に沿った超音波信号の後方散乱によって影響を受けない。従って、上記の組織特性推定は、従来の縦波ベースの組織減衰測定の前述の欠点を持たない。
従来のせん断波エラストグラフィ技法と比較して、関心領域の追加の特性は、複数のせん断波の測定されたパラメータを組み合わせることにより推定されることができる。いくつかの実施形態では、上記第2のパラメータは、例えば組織長手方向減衰、組織長手方向粘性等の上記関心領域長手方向の特性を示すことができる。良く知られるように、組織特性は、すべての方向で同じであってもよいが、異なる方向において異なることもできる。例えば、長手方向における組織の減衰、組織の粘性又は組織の弾性といった組織の機械的特性は、横方向の特性とは異なっていてもよい。関心領域の長手方向特性は、例えば長手方向の減衰、長手方向の粘度といった長手方向に沿った関心領域の特性(即ち、長手方向の波の伝搬)として知られる。同様に、せん断特性とも呼ばれる横方向特性は、例えば横方向減衰(せん断減衰とも呼ばれる)、横方向粘度など、横方向に沿った関心領域の特性として知られる。
焦点の数は、任意の適切な値に設定されることができる。一般に、焦点の数が多いほど、関心領域の推定された特性がより正確になる。超音波プローブは、1次元又は2次元のトランスデューサアレイを含むことができる。関心領域は、例えば、肝臓組織、腎臓組織、前立腺組織、乳房組織といった組織の領域とすることができる。対象は、患者といったヒト、動物、ファントム等とすることができる。焦点で生成されたせん断波の特性は、時間領域又は周波数領域における特性とすることができ、これは、せん断波の変位、速度、及び/又はスペクトルを含むがこれに限定されるものではない。
一実施形態によれば、関心領域の特性は、関心領域の減衰、剛性及び粘度のいずれか1つを含む。関心領域の特性は、組織減衰、組織剛性及び組織粘性のいずれか1つを含むことができる。組織の減衰は、長手方向の減衰又はせん断減衰であり得る。関心領域の特性は、関心領域の長手方向特性を含むことができる。例えば、関心領域の特性は、長手方向の減衰であり得る。
特別なせん断波励起及び追跡スキームを開発することにより、縦波減衰といった追加パラメータの推定を可能にする統合アプローチを提案する。提案されたアプローチは、縦波及びせん断波の両方に基づき多パラメータ組織定量化を提供するワンボタンプッシュツールとして構想される。臨床研究は、縦波の音響減衰係数が軟組織における脂肪定量化(例えば脂肪肝定量)に有用であることを長い間指摘してきた。最近開発されたせん断波エラストグラフィは、組織特徴化(例えば、肝線維症病期分類)及び癌検出(例えば、肝臓、乳房、前立腺及び甲状腺などの癌検出)のための診断ツールとして登場した。
一実施形態によれば、上記焦点に関する第1のパラメータは、焦点に隣接する追跡スポットにおける上記生成されたせん断波の変位から得られる。言い換えると、各第1のパラメータは、生成されたせん断波がそれぞれの焦点に隣接する追跡スポットで示す変位から得られる。時間領域又は周波数領域のいずれかにおいて、エコー信号からのせん断波変位を得るのに、様々な既知の手法が使用されることができる。様々な実施形態において、第1のパラメータは、追跡スポットにおける生成されたせん断波のピーク変位、平均変位などを示す。追跡スポットが、それぞれの焦点に隣接し、追跡スポットでの変位は、焦点における変位に比例すると考えられるので、得られた第1のパラメータは、生成されたせん断波のそれぞれの焦点位置における例えばピーク変位、平均変位等の変位を示すために用いられることができる。
一実施形態によれば、第1のパラメータは、焦点における生成されたせん断波のピーク変位を示すため、追跡スポットにおける生成されたせん断波のピーク変位から得られる。
一実施形態によれば、上記特性推定器は、上記得られた第1のパラメータを用いて指数曲線を適合させることにより、上記関心領域の組織長手方向減衰特性を推定するよう構成される。
一実施形態によれば、上記超音波システムは、インターフェースを更に有し、上記インターフェースが、上記複数の焦点を示すための入力を受信するよう構成され、上記特性推定器は更に、上記示された複数の焦点のそれぞれの焦点深度値に基づき、上記第2のパラメータを推定するよう構成される。一例では、ユーザ入力は関心領域とすることができ、複数の焦点は、所定の規則に基づき自動的に設定されることができる。例えば、複数の焦点は、関心領域にわたって等間隔に配置される所定数のスポットとして設定されることができる。別の例では、入力は、複数の焦点とすることができる。インターフェースは、ユーザインターフェース、又は病院情報システム、患者情報管理システムといった必要とされる情報を提供することができるシステムとのインターフェースとすることができる。
一実施形態によれば、超音波システムは、複数の焦点の各々に関して、プッシュパルスのビーム形成設定によってもたらされるプッシュパルスの強度変動を推定するよう構成される強度変動推定器を更に有する。
本発明の発明者らは、関心領域の特定の特性の推定は、プッシュパルスの初期音響強度を利用することができるが、プッシュパルスの初期音響強度は、プッシュパルスのビームプロファイルとともに変化することを認識した。関心領域の特性に関する推定を更に向上させるため、本発明者らは、斯かる変動を推定し、その変動を補償することを提案する。
いくつかの実施形態では、上記超音波プローブが、上記推定された変動に基づき上記プッシュパルスの強度を調整するよう更に構成される。いくつかの他の実施形態では、上記特性推定器が、上記推定された変動に基づき上記第2のパラメータを推定するよう更に構成される。
一実施形態によれば、上記せん断波検出器が、上記複数の焦点の各々に関して、上記焦点の深度における上記関心領域の組織せん断を示す第3のパラメータを推定し、上記推定された上記第3のパラメータに基づき上記第1のパラメータを調整するよう更に構成され、上記特性推定器は、上記調整された第1のパラメータの関数として、上記第2のパラメータを得るよう更に構成される。例えば、関心領域のせん断特性は、関心領域のせん断剛性、せん断粘度、又はせん断減衰とすることができる。例えば、各焦点の深度におけるせん断特性は、その焦点で生成された個別のせん断波の速度に基づき推定されることができる。
本発明者らは、異なる深度における局部せん断特性が、異なる場合があり、従って関心領域の特性の推定に影響を及ぼす場合があることを認識し、そこで、関心領域の特性の推定の前に、第1のパラメータを局所的なせん断波特性を用いて調整することを提案する。
一実施形態によれば、超音波プローブは、複数の焦点の各々に関して、この焦点に隣接する追跡スポットに向かって追跡パルスを送信し、追跡スポットから超音波エコー信号を受信するよう構成され、上記せん断波検出器は、上記複数の焦点のそれぞれに関して、上記焦点に隣接する追跡スポットから受信される上記超音波エコー信号に基づき、上記第1のパラメータを得るよう構成される。
本発明の第2の側面の実施形態によれば、せん断波を用いて対象の関心領域の特性を測定する方法が提案される。この方法は、それぞれが相互に異なる深度値を持つ上記関心領域における複数の焦点の各々に、せん断波を発生させるためのプッシュパルスを順次送信し、上記複数の焦点のそれぞれに隣接する超音波エコー信号を受信するステップと、上記複数の焦点の各々に関して、上記受信された超音波エコー信号に基づき、上記生成されたせん断波が上記焦点において持つ特性を示す第1のパラメータを得るステップと、上記得られた第1のパラメータの関数として、上記関心領域の特性を示す第2のパラメータを推定するステップとを有する。
第3の本発明の側面の実施形態によれば、実行されるときコンピュータに、複数のせん断波に基づき、対象における関心領域の特性を推定する方法を実行させるコンピュータプログラムが提供され、複数のせん断波の各々が、異なる深度の焦点にプッシュパルスを送信することにより生成され、上記方法が、複数の焦点の各々に関して、上記複数の焦点のそれぞれに隣接して受信される超音波エコー信号に基づき、上記焦点において上記せん断波が持つ特性を示す第1のパラメータを得るステップと、上記得られた第1のパラメータの関数として、上記関心領域の特性を示す第2のパラメータを推定するステップとを有する。
従来技術における従来のせん断波ベースの組織の粘性又は弾性測定を示す図である。 本発明のいくつかの実施形態に従って構成された超音波撮像システムをブロック図形式で示す図である。 本発明のいくつかの実施形態による、関心領域における複数の焦点で生成された複数のせん断波を示す図である。 本発明のいくつかの実施形態による、関心領域における複数の焦点を示す図である。 本発明の一実施形態による、焦点深度で集束されたプッシュパルスのビームプロファイルを示す図である。
本発明の他の目的及び利点が、添付の図面と共になされる説明を参照して、より明らかとなり、これを用いて容易に理解される。
本発明が、以下、実施形態及び図面を参照し、これらを組み合わせてより詳細に記載及び説明される。
図面における同じ参照符号は、類似する又は対応する特徴及び/又は機能を示す。
本発明は、特定の実施形態及び特定の図面を参照して説明されることになるが、本発明はそれらに限定されるものではなく、添付された請求項によってのみ定まるものである。記載された図面は、概略的なものに過ぎず非限定的なものである。図面において幾つかの要素の大きさが誇張されている場合があり、説明目的のため実際のスケール通りに描かれていない場合がある。
まず図2を参照すると、本発明のいくつかの実施形態に従って構成された超音波システムがブロック図形式で示される。
超音波プローブ10は、超音波信号を送信及び受信するトランスデューサ要素の1次元又は2次元アレイ12を持つ。トランスデューサアレイ12は、ビームを送信し、身体内の単一平面上で反射エコー信号を受信することにより、2次元(2D)平面を走査することができ、体のボリュメトリック(3D)領域の異なる方向及び/又は平面においてビームを送信することにより、ボリュメトリック領域を走査するのに使われることもできる。アレイ要素は、プローブ内に配置されたマイクロビーム形成器38に結合される。これは、この要素による送信を制御し、この要素のグループ又はサブアレイから受信されるエコー信号を部分的にビーム形成された信号へと処理する。部分的にビーム形成された信号は、送信/受信(T/R)スイッチ14により、プローブから超音波システム内のマルチライン受信ビーム形成器20に結合される。ビーム形成器による送信及び受信の調整は、ビーム形成器コントローラ16により制御され、これは、マルチライン受信ビーム形成器及び送信コントローラ18に結合される。送信コントローラは、制御信号をマイクロビーム形成器に提供する。ビーム形成器コントローラは、ユーザ制御パネル40のユーザ操作に基づき生成される信号に応答して、超音波システム及びそのプローブの動作を制御する。マルチライン受信ビーム形成器20は、単一の送信−受信間隔中にエコー信号の1つ又は複数の空間的に分離した追跡ライン(受信ライン又はAラインとしても知られる)を生成する。
本発明の一実施形態によれば、超音波プローブは、関心領域内の複数の焦点の各々に、せん断波を生成するためのプッシュパルスを順次送信し、上記複数の焦点の各々に隣接する超音波エコー信号を受信するよう構成される。複数の焦点の各々は、互いに異なる焦点深度値を持つ。例えば、超音波エコー信号は、各焦点に隣接するAラインに沿って受信される。
受信されたエコー信号は、更なる処理のため、せん断波検出器に結合される。せん断波検出器は、受信した超音波エコー信号に基づき、複数の焦点のそれぞれについて、個別のせん断波の特性の第1のパラメータを得るよう構成される。せん断波検出器は、信号プロセッサ22と、Aラインメモリ24と、Aライン動き推定器(例えば、相互相関器)26と、強度検出器28とを有する。エコー信号は、信号プロセッサ22により、フィルタリング、ノイズ除去等の処理が施され、Aラインメモリ24に格納される。連続するAラインサンプリングのエコー信号が、Aライン動き推定器26によって追跡され、各サンプリング位置について組織変位のサンプルの時間シーケンスが生成される。強度検出器28は、Aラインに沿ったせん断波変位の検出に基づき、せん断波の強度を決定する。特性推定器32は、得られた第1のパラメータの関数に基づき、関心領域の特性を示す第2パラメータを推定するように結合される。せん断波検出器及び特性推定器32は、図3を参照してより詳細に説明される。
いくつかの実施形態では、得られた特性情報は、画像プロセッサ34に結合されることができる。これは、画像ディスプレイ36に表示するため、好ましくは組織の解剖学的超音波画像と組み合わせて、特性情報を処理する。
いくつかの実施形態では、超音波システムは、プッシュパルスのビームプロファイルによってもたらされるプッシュパルスの初期音響強度の変化を、各プッシュパルスについて推定するよう構成される強度変動推定器40を更に有することができる。一例では、超音波プローブ10は、推定された変動に基づきプッシュパルスの強度を調整するよう更に構成される。この場合、強度変動推定器40は、ビーム形成器コントローラ16の一部とすることができる。別の例では、特性推定器32は、推定された変動に基づき第2のパラメータを推定するよう更に構成される。ビームプロファイルによりもたらされる斯かる強度変動及び対応する補償が、図5を参照してより詳細に説明される。
図3は、本発明の一実施形態による関心領域内の複数の焦点で生成された複数のせん断波を示す。3つのせん断波が図3に示されるが、異なる数のせん断波が、様々な実施形態に基づき生成されてもよい。
図3を参照すると、複数のプッシュパルス310、312、314の各々が、皮膚表面下の異なる深度z1、z2、z3で集束される。各焦点320、322、324においてもたらされる組織変形は、長手方向zに直交する横方向(例えば図示される方向x)に沿って伝搬するせん断波330、332、334を生じさせる。
プッシュパルスは、集束された高MIプッシュパルス(例えば、FDA診断限度内にあるように1.9以下のMI)とすることができる。これは、焦点において組織を下方に変位させ、せん断波の発生を引き起こす。典型的には、超音波信号、例えばプッシュパルスは、組織内で約1560メートル/秒の速度で移動し、せん断波は、組織内で約1〜5メートル/秒の速度で移動し、組織内のせん断波周期は、2〜10ミリ秒のオーダーである。いくつかの実施形態では、複数のプッシュパルスが、所定の時間間隔で逐次送信されることができる。例えば、所定の時間間隔は、せん断波の少なくとも1周期より長くすることができる。
図4は、本発明のいくつかの実施形態による関心領域内の複数の焦点を示す。図4において、方向zは、超音波プローブ410の長手方向、即ち深度方向としても知られる超音波信号の伝搬方向を表す。方向x及びyは、超音波プローブ410の長手方向に垂直な方向、即ち、方位角及び仰角方向をそれぞれ表す。2次元トランスデューサアレイを含む超音波プローブ410が、3次元撮像野420を持つことが示されるが、1次元トランスデューサアレイを持つ超音波プローブが同様に使用され、関心領域は、2次元又は3次元であり得る点を理解されたい。
図4を参照すると、関心領域430内の複数の空間位置440が焦点として使用される。3つの焦点が図示されるが、本発明はこれに限定されず、任意の数の焦点が規定されることができる。焦点440は、異なる深度に、言い換えると、長手方向zに沿って異なる位置に配置される。一実施形態において、複数の焦点は、長手方向に沿って互いに同じ所定の距離だけ離間されることができる。所定の距離は、例えば、1センチメートルのオーダーとすることができる。複数の焦点440は、様々な実施形態において、他の2つの方向x及びyに沿って同じ位置又は異なる位置を持つことができる。
複数の焦点は、ユーザインターフェースを介して手動で選択されることができる。例えば、関心領域を含む超音波画像がディスプレイに提示され、次いでユーザは、任意の適切なユーザ入力手段を介して、表示された超音波画像における対応する空間位置を入力することにより、焦点を選択することができる。この入力手段は、ジョイスティック、マウス、物理的なボタン、タッチスクリーンなどを含むが、これらに限定されるものではない。追加的又は代替的に、焦点が関心領域において自動的に生成されることができる。例えば、所定の数の焦点が生成されることができ、及び/又は同じ所定の距離だけ離れた複数の焦点が生成されることができる。所定の数及び/又は所定の距離は、ユーザによって選択される、又は予め規定されることができる。
関心領域は、表示された超音波画像におけるユーザインターフェースを介して手動で選択されることもでき、及び/又は例えば超音波画像を処理することにより自動的若しくは半自動的に規定されることもできる。
再び図3を参照すると、追跡スポット350、352、354における組織変位を追跡するため、焦点に隣接する追跡ライン340、342、344に沿って超音波エコー信号が受信される。本発明の発明者は、例えば焦点に隣接する追跡スポットでのプッシュパルス反射又は残響のため、焦点でのせん断波追跡がそれほど正確ではないことを認識し、従って、焦点に隣接する追跡スポットで、せん断波変位を追跡することを提案する。更に、追跡スポットは、焦点から離れすぎてはならない。なぜなら、せん断波変位が伝搬方向に沿って急速に減少するからである。各焦点と最も近い追跡スポットとの間の距離は、任意の適切な値とすることができ、異なる種類の関心領域に対して異なることができる。例えば、追跡スポットと焦点との間の距離は、0.1cmのオーダーとすることができる。例えば、追跡スポットと焦点との間の距離は、プッシュフォーカスから離れたプッシュビーム幅の半分に設定されることができる。
複数の焦点の各々について、せん断波の特性を示す第1のパラメータが、それぞれの追跡ラインに沿って受信されたエコー信号から得られる。焦点320を例にとると、焦点320に対して得られる第1のパラメータは、追跡ライン340に沿って受信されたエコー信号から得られる。いくつかの実施形態では、焦点について得られた第1のパラメータは、その焦点における個別のせん断波の特性値を示すことができる。本発明の一実施形態によれば、第1のパラメータは、追跡スポット350、352、354におけるせん断波変位のピーク値として得られる。他の例では、第1のパラメータは、せん断波の平均変位などとして得られることができる。
次に、関心領域の特性が、複数の焦点におけるせん断波特性値の関数として推定される。いくつかの実施形態では、各焦点におけるせん断波特性値は、プッシュパルスの強度、関心領域の特性、及び焦点の深度値の関数としてモデル化され、斯かるモデルに基づき、関心領域の未知の特性が得られることができる。得られる未知の特性は、関心領域内で同じであると仮定される点に留意されたい。一例として組織の減衰を考慮して、一実施形態による推定が以下に説明される。
本発明の一実施形態によれば、個別の追跡スポットで測定されたk番目のせん断波のピーク又は最大変位Mkは、プッシュパルスの強度Ikに指数式を掛けたものに比例するものとしてモデル化されることができる。指数式における指数は、それぞれの焦点の深度dkと負の組織減衰係数αとの積であり、即ち
Figure 0006259953
となる。
例えば、推定された組織減衰係数αは、数学的に
Figure 0006259953
として得られることができる。
例えば、組織減衰係数の推定は、得られた第1のパラメータMkを用いて指数曲線をフィットさせることにより実行されることができる。いくつかの実施形態では、プッシュパルス310、312、314の強度Ikが実質的に同じであるよう設定されることができる。いくつかの他の実施形態では、プッシュパルス310、312、314の強度が異なるよう設定されることができる。
本発明の発明者は、関心領域の特定の特性の推定が、プッシュパルスの初期音響強度を利用することができるが、プッシュパルスの初期音響強度は、プッシュパルスのビームプロファイルと共に変化することを認識した。例えば、本発明者らは、深度値dを持つ焦点におけるピーク又は最大せん断波変位Mdが、プッシュパルスの初期音響強度Iに指数式を乗じたものに比例し、指数式における指数が、負の組織減衰係数αと深度dとの積であり、即ち
Figure 0006259953
であることを理論的に得た。
本発明者らは、関心領域の特性に関する推定を更に向上させるため、斯かる変動を推定し、この変動を補償することを提案する。
図5は、本発明の一実施形態による焦点深度で集束されたプッシュパルスのビームプロファイルを示す。図5を参照すると、皮膚表面520における超音波プローブ510が、ビームプロファイル530を持つプッシュパルスを体500を通して関心領域内の所与の焦点に送信する。焦点は、焦点深度値dを持つ陰影領域550によって示される。このプッシュパルスは、焦点における組織を下方向に変位させ、変位された組織から外側に発するせん断波面540を生じさせる。特定のビームプロファイルを生成するために、異なるトランスデューサ要素が使用されることができ、異なるエネルギー重み、位相シフトが、使用されるトランスデューサ要素に適用されることができる。これは、プッシュパルスの強度が同じに設定されるとしても、プッシュパルスの初期音響強度Iが異なることができることを生じさせる。初期音響強度の変動は、所与のビームプロファイルに対して決定されることができる。
本発明の一実施形態によれば、強度推定器が、斯かる強度変化を推定するよう構成される。例えば、ビームプロファイルによって提供される場合、強度推定器は、所定のモデル化に基づきオンライン計算により、又は事前に格納されたルックアップテーブルから強度変動を取得することにより、このビームプロファイルの強度変動を推定するよう構成することができる。
本発明の実施形態によれば、超音波プローブは、推定された強度変動に基づきプッシュパルスの強度を調整するよう更に構成される。例えば、所定の初期音響強度値を得るため、所定の初期音響強度値と、ビームプロファイルによりもたらされる推定された強度変動とに基づき、あるビームプロファイルを持つプッシュパルスの強度が決定され、超音波プローブは、この決定された強度でプッシュパルスを送信する。代替的に、特性推定器は、推定された変動に基づき第2のパラメータを推定するよう更に構成される。例えば、ビームプロファイルを持つプッシュパルスに関して、補償された初期音響強度値は、プッシュパルスの強度と、ビームプロファイルによりもたらされる推定された強度変動とに基づき決定されることができ、特性推定器は、補償された初期音響強度値を用いて関心領域の特性を推定する。
一実施形態によれば、せん断波検出器は、複数の焦点のそれぞれに対して、焦点の深度における関心領域のせん断特性を示す第3パラメータを推定し、推定された第3のパラメータに基づき、第1のパラメータを調整するよう更に構成され、特性推定器は更に、調整された第1のパラメータの関数として上記第2のパラメータを得るよう構成される。
例えば、関心領域のせん断特性は、関心領域のせん断剛性又は粘度とすることができる。例えば、各焦点の深度におけるせん断特性は、その焦点で生成された個別のせん断波の速度に基づき推定されることができる。
焦点の深度における関心領域のせん断特性は、従来のせん断波エラストグラフィ技術又は将来開発される他の適切な手段によって推定されることができる。従来のせん断波エラストグラフィ技術を使用する場合、生成された各せん断波に対して、1つ以上の追跡スポットが、個別の焦点から離れて横方向xに沿って離間され、超音波エコー信号は、この1つ以上の追跡スポットに沿って受信され、各せん断波の速度が、受信された超音波エコー信号から得られ、組織せん断剛性及び粘性といったせん断特性が、個別のせん断波の得られた速度から推定される。
本書に記載される技術プロセスは、様々な手段により実現されることができる。例えば、これらの技術は、ハードウェア、ソフトウェア、又はこれらの組み合わせで実現されることができる。ハードウェアによる実現の場合、処理ユニットは、1つ又は複数の特定用途向け集積回路(ASIC)、デジタル信号プロセッサ(DSP)、デジタル信号処理デバイス(DSPD)、プログラマブルロジックデバイス(PLD)、フィールドプログラマブルゲートアレイ(FPGA)、プロセッサ、コントローラ、マイクロコントローラ、マイクロプロセッサ、本書に記載される機能を実行するよう設計された他の電子ユニット、又はこれらの組み合わせにおいて実現されることができる。ソフトウェアの場合、実現は、本書に記載される機能を実行するモジュール(例えば、プロシージャ、関数等)を介して行われることができる。ソフトウェアコードは、メモリユニットに記憶され、プロセッサにより実行されることができる。
更に、請求項に記載される主題の側面は、この主題の様々な側面を実現するようコンピュータ又はコンピューティングコンポーネントを制御するためのソフトウェア、ファームウェア、ハードウェア、又はこれらの任意の組み合わせを生成するための標準プログラミング及び/又はエンジニアリング技術を用いて、方法、装置、システム又は製品として実現されることができる。本書で使用される用語「製品」は、任意のコンピュータ可読デバイス、担体又は媒体からアクセス可能なコンピュータプログラムを包含するものとして意図される。例えば、コンピュータ可読媒体は、磁気記憶デバイス(例えば、ハードディスク、フロッピーディスク、磁気ストリップ...)、光ディスク(例えば、コンパクトディスク(CD)、デジタル多用途ディスク(DVD)...)、スマートカード、及びフラッシュメモリデバイス(例えば、カード、スティック、キードライブ...)を含むことができるが、これらに限定されるものではない。もちろん、当業者であれば、本書に記載された範囲又は趣旨から逸脱することなく、多くの変更がこの構成に加えられることができることを理解されるであろう。
本願において用いられる用語「ビーム形成器」、「コントローラ」、「プロセッサ」、「相互相関器」、「検出器」、「推定器」、例えば「せん断波検出器」、「特性推定器」及び「撮像エンコーダ」は、ハードウェア、ハードウェアとソフトウェアとの組み合わせ、ソフトウェア、又は実行中のソフトウェアのいずれかであるプロセッサ又はコンピュータ関連エンティティを指すものとして意図される。例えば、コンポーネントは、プロセッサで実行中のプロセス、プロセッサ、オブジェクト、実行可能ファイル、実行スレッド、プログラム、及び/又はコンピュータとすることができるが、これに限定されるものではない。例として、サーバー上で実行中のアプリケーション及びサーバーの両方が、コンポーネントになることができる。1つ又は複数のコンポーネントが、プロセス及び/又は実行スレッド内に存在することができ、コンポーネントは、1つのコンピュータ上にローカライズされ、及び/又は2つ以上のコンピュータ間に分散されることができる。
上述されたことは、1つ又は複数の実施形態の例を含む。もちろん、前述の実施形態を説明する目的で構成要素又は方法論の考えられるすべての組み合わせを説明することは不可能であるが、当業者であれば、様々な実施形態の多くのさらなる組み合わせ及び置換が可能であることを理解されるであろう。従って、記載された実施形態は、添付の特許請求の範囲の趣旨及び範囲に含まれる斯かる変更、修正及び変形をすべて包含するものとして意図される。更に、用語「含む(include)」は、詳細な説明又は特許請求の範囲のいずれかで使用される限り、斯かる用語は、用語「有する(comprising)」が、請求項において伝統的な用語として使用されるとき解釈されるのと同様の態様で包括的であるものとして意図される。

Claims (15)

  1. せん断波を用いて対象の関心領域の特性を測定する超音波システムであって、
    それぞれが相互に異なる深度値を持つ前記関心領域における複数の焦点の各々に、せん断波を発生させるためのプッシュパルスを順次送信し、前記複数の焦点のそれぞれに隣接する超音波エコー信号を受信するよう構成される超音波プローブと、
    前記受信された超音波エコー信号に基づき、前記生成されたせん断波が前記焦点で持つ特性を示す第1のパラメータを、複数の焦点の各々に関して得るよう構成されるせん断波検出器と、
    前記複数の焦点における前記得られた第1のパラメータの関数として、前記関心領域の特性を示す第2のパラメータを推定するよう構成される特性推定器とを有する、超音波システム。
  2. 前記関心領域の特性が、前記関心領域の長手方向特性を含み、前記関心領域の長手方向特性は、前記プッシュパルスの伝搬方向に沿った関心領域の特性である、請求項1に記載の超音波システム。
  3. 前記焦点に関して得られた第1のパラメータが、前記生成されたせん断波が焦点で持つ強度を示す、請求項1に記載の超音波システム。
  4. 前記焦点に関する第1のパラメータが、前記焦点に隣接する追跡スポットにおける前記生成されたせん断波の変位から得られる、請求項1乃至3の任意の一項に記載の超音波システム。
  5. 前記第1のパラメータが、前記追跡スポットにおける前記生成されたせん断波のピーク変位から得られる、請求項4に記載の超音波システム。
  6. 前記特性推定器が、前記得られた第1のパラメータを用いて指数曲線を適合させることにより、前記関心領域の組織長手方向減衰特性を推定するよう構成される、請求項5に記載の超音波システム。
  7. インターフェースを更に有し、
    前記インターフェースが、前記複数の焦点を示すための入力を受信するよう構成され、
    前記特性推定器は更に、前記示された複数の焦点のそれぞれの焦点深度値に基づき、前記第2のパラメータを推定するよう構成される、請求項1乃至3の任意の一項に記載の超音波システム。
  8. 前記特性推定器が、前記得られた第1のパラメータ及びプッシュパルスの強度の関数として、前記第2のパラメータを推定するよう構成される、請求項1乃至3の任意の一項に記載の超音波システム。
  9. 前記プッシュパルスのビームプロファイルによりもたらされる前記プッシュパルスの初期音響強度の変動を、各プッシュパルスに関して推定するよう構成される強度変動推定器を更に有する、請求項8に記載の超音波システム。
  10. 前記超音波プローブが、各プッシュパルスに関して、前記推定された変動に基づき前記プッシュパルスの強度を調整するよう更に構成される、請求項9に記載の超音波システム。
  11. 前記特性推定器が、前記推定された変動に基づき前記第2のパラメータを推定するよう更に構成される、請求項9に記載の超音波システム。
  12. 前記せん断波検出器が、前記複数の焦点の各々に関して、前記焦点の深度における前記関心領域のせん断特性を示す第3のパラメータを推定し、前記推定された前記第3のパラメータに基づき前記第1のパラメータを調整するよう更に構成され、
    前記特性推定器は、前記複数の焦点における前記調整された第1のパラメータの関数として、前記第2のパラメータを得るよう更に構成される、請求項1乃至3の任意の一項に記載の超音波システム。
  13. せん断波を用いて対象の関心領域の特性を測定する方法において、
    それぞれが相互に異なる深度値を持つ前記関心領域における複数の焦点の各々に、せん断波を発生させるためのプッシュパルスを順次送信し、前記複数の焦点のそれぞれに隣接する超音波エコー信号を受信するステップと、
    前記複数の焦点の各々に関して、前記受信された超音波エコー信号に基づき、前記生成されたせん断波が前記焦点において持つ特性を示す第1のパラメータを得るステップと、
    前記複数の焦点における前記得られた第1のパラメータの関数として、前記関心領域の特性を示す第2のパラメータを推定するステップとを有する、方法。
  14. 前記第2のパラメータが、前記関心領域の長手方向特性を示し、前記関心領域の長手方向特性は、前記プッシュパルスの伝搬方向に沿った前記関心領域の特性である、請求項13に記載の方法。
  15. 実行されるときコンピュータに、複数のせん断波に基づき、対象における関心領域の特性を推定する方法を実行させるコンピュータプログラムであって、複数のせん断波の各々が、異なる深度の焦点にプッシュパルスを送信することにより生成され、前記方法が、
    複数の焦点の各々に関して、前記複数の焦点のそれぞれに隣接して受信される超音波エコー信号に基づき、前記焦点において前記せん断波が持つ特性を示す第1のパラメータを得るステップと、
    前記複数の焦点における前記得られた第1のパラメータの関数として、前記関心領域の特性を示す第2のパラメータを推定するステップとを有する、コンピュータプログラム。
JP2017522080A 2015-08-03 2016-07-29 せん断波を用いる測定に関する超音波システム及び方法 Expired - Fee Related JP6259953B2 (ja)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
CN2015085935 2015-08-03
CNPCT/CN2015/085935 2015-08-03
EP15186920.3 2015-09-25
EP15186920 2015-09-25
PCT/EP2016/068250 WO2017021341A1 (en) 2015-08-03 2016-07-29 Ultrasound system and method for measurement using shear wave

Publications (2)

Publication Number Publication Date
JP2017535327A JP2017535327A (ja) 2017-11-30
JP6259953B2 true JP6259953B2 (ja) 2018-01-10

Family

ID=56787422

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2017522080A Expired - Fee Related JP6259953B2 (ja) 2015-08-03 2016-07-29 せん断波を用いる測定に関する超音波システム及び方法

Country Status (6)

Country Link
US (1) US10729404B2 (ja)
EP (1) EP3190976B1 (ja)
JP (1) JP6259953B2 (ja)
CN (1) CN107072639B (ja)
RU (1) RU2641067C1 (ja)
WO (1) WO2017021341A1 (ja)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11523774B2 (en) * 2017-04-06 2022-12-13 Siemens Medical Solutions Usa, Inc. Tissue property estimation with ultrasound medical imaging
US11553901B2 (en) 2017-04-06 2023-01-17 Siemens Medical Solutions Usa, Inc. Liver disease activity estimation with ultrasound medical imaging
US20180310918A1 (en) * 2017-04-27 2018-11-01 Siemens Medical Solutions Usa, Inc. Variable focus for shear wave imaging
US11726195B2 (en) * 2017-09-19 2023-08-15 Koninklijke Philips N.V. Ultrasound systems and methods for measuring acoustic attenuation coefficients
CN111885965A (zh) * 2018-03-21 2020-11-03 皇家飞利浦有限公司 用于三维中的剪切波成像的超声系统
WO2020037674A1 (zh) * 2018-08-24 2020-02-27 深圳迈瑞生物医疗电子股份有限公司 一种超声成像装置和方法、超声弹性检测装置和方法
CN110742644B (zh) * 2019-09-29 2022-07-08 深圳大学 弹性成像系统、弹性成像方法及存储介质
KR102617857B1 (ko) * 2020-03-12 2023-12-22 지멘스 메디컬 솔루션즈 유에스에이, 인크. 초음파 의료 이미징을 사용한 간 질환 활성도 추정

Family Cites Families (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2232547C2 (ru) * 2002-03-29 2004-07-20 Общество с ограниченной ответственностью "АММ - 2000" Способ и устройство для получения ультразвуковых изображений структур и сосудов головного мозга
WO2008141220A1 (en) * 2007-05-09 2008-11-20 University Of Rochester Shear modulus estimation by application of spatially modulated impulse acoustic radiation force approximation
US9364194B2 (en) * 2008-09-18 2016-06-14 General Electric Company Systems and methods for detecting regions of altered stiffness
US8602994B2 (en) 2009-03-09 2013-12-10 Mayo Foundation For Medical Education And Research Method for ultrasound vibrometry using orthogonal basis functions
US8343050B2 (en) 2009-05-04 2013-01-01 Siemens Medical Solutions Usa, Inc. Feedback in medical ultrasound imaging for high intensity focused ultrasound
US8734352B2 (en) 2009-07-17 2014-05-27 Koninklijke Philips N.V. Spatially-fine shear wave dispersion ultrasound vibrometry sampling
US8500639B2 (en) 2009-09-11 2013-08-06 Mr Holdings (Hk) Limited Systems and methods for shear wave field formation
CN102667522B (zh) * 2009-11-25 2014-10-08 皇家飞利浦电子股份有限公司 采用聚焦扫描线波束形成的超声剪切波成像
US8753277B2 (en) * 2009-12-10 2014-06-17 The University Of Rochester Methods and systems for spatially modulated ultrasound radiation force imaging
US9986973B2 (en) * 2010-04-23 2018-06-05 Mayo Foundation For Medical Education And Research Method for shear wave ultrasound vibrometry with interleaved push and detection pulses
US8961418B2 (en) 2010-10-06 2015-02-24 Siemens Medical Solutions Usa, Inc. Solving for shear wave information in medical ultrasound imaging
US8758248B2 (en) 2010-11-30 2014-06-24 General Electric Company Systems and methods for acoustic radiation force imaging with enhanced performance
EP2649944A4 (en) * 2010-12-08 2017-02-22 Hitachi, Ltd. Ultrasound diagnosis apparatus
JP5868419B2 (ja) * 2010-12-13 2016-02-24 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 超音波材料特性測定と画像化のための超音波音響放射力励起
US8494791B2 (en) * 2010-12-16 2013-07-23 General Electric Company Methods and systems for improved correlation of shear displacement waveforms
CN103347450B (zh) * 2011-02-04 2015-07-08 株式会社日立医疗器械 超声波诊断装置及方法
JP6067590B2 (ja) 2011-02-25 2017-01-25 メイヨ フォンデーシヨン フォー メディカル エジュケーション アンド リサーチ 非合焦超音波による超音波振動法
US8734350B2 (en) * 2011-03-04 2014-05-27 Mayo Foundation For Medical Education And Research System and method for correcting errors in shear wave measurements arising from ultrasound beam geometry
US9237878B2 (en) 2011-04-22 2016-01-19 Mayo Foundation For Medical Education And Research Generation and assessment of shear waves in elasticity imaging
US9239373B2 (en) 2011-11-16 2016-01-19 Siemens Medical Solutions Usa, Inc. Adaptive image optimization in induced wave ultrasound imaging
US9244169B2 (en) 2012-06-25 2016-01-26 Siemens Medical Solutions Usa, Inc. Measuring acoustic absorption or attenuation of ultrasound
WO2014055410A1 (en) 2012-10-01 2014-04-10 Mayo Foundation For Medical Education And Research Shear wave attenuation from k-space analysis system
EP2833792B1 (en) 2012-10-07 2024-06-12 Mayo Foundation For Medical Education And Research Method for shear wave elastography by transmitting ultrasound with subgroups of ultrasound transducer elements
US10743814B2 (en) 2013-03-15 2020-08-18 Siemens Medical Solutions Usa, Inc. Fat fraction estimation using ultrasound with shear wave propagation
FR3017041B1 (fr) * 2014-01-31 2016-03-04 Centre Nat Rech Scient Procede et dispositif ultrasonore de caracterisation des milieux mous anisotropes, et ensemble de sonde ultrasonore pour un tel dispositif de caracterisation

Also Published As

Publication number Publication date
US20180228463A1 (en) 2018-08-16
WO2017021341A1 (en) 2017-02-09
EP3190976A1 (en) 2017-07-19
JP2017535327A (ja) 2017-11-30
US10729404B2 (en) 2020-08-04
RU2641067C1 (ru) 2018-01-15
EP3190976B1 (en) 2018-03-14
CN107072639A (zh) 2017-08-18
CN107072639B (zh) 2019-03-08

Similar Documents

Publication Publication Date Title
JP6259953B2 (ja) せん断波を用いる測定に関する超音波システム及び方法
US11446006B2 (en) Adjusting measurements of the effects of acoustic radiation force for background motion effects
Deng et al. Ultrasonic shear wave elasticity imaging sequencing and data processing using a verasonics research scanner
JP6053860B2 (ja) 医療用超音波イメージングシステムにおいて剪断波情報を求めるためのコンピュータプログラム
JP6230827B2 (ja) 超音波の音響吸収量の測定方法、超音波の音響吸収量を測定するシステム、および、超音波の音響減衰量を測定するためのコンピュータで読み出し可能な非変換性記憶媒体
EP2881041B1 (en) Apparatus and method for ultrasonic diagnosis
JP5829882B2 (ja) せん断変位波形の相関を向上させた方法
JP6140953B2 (ja) 医用超音波剪断波イメージングにおける分類前処理
US10448924B2 (en) Ultrasonic acoustic radiation force excitation for ultrasonic material property measurement and imaging
CN102469980B (zh) 空间上精细的剪切波分散超声振动测定采样
JP2010526626A (ja) 関心領域の粘弾性の平均値を測定するための方法および装置
KR20150037689A (ko) 의료용 초음파 영상의 전단파 검출
CN106419961A (zh) 声学辐射力成像中的自适应运动估计
KR20190049566A (ko) 초음파 의료 이미징에서 전단 속도로부터의 조직 점탄성 추정
KR20220036339A (ko) 매질의 초음파 특성화 방법 및 시스템
CN109310401A (zh) 使用反射剪切波来监测热消融中的损伤生长
CN111885965A (zh) 用于三维中的剪切波成像的超声系统
KR20220036338A (ko) 매질의 초음파 특성화 방법 및 시스템
KR20220036863A (ko) 매질의 초음파 특성화 방법 및 시스템
CN111970972A (zh) 以超声粒子速度估计器映射的声波场
JP2020536668A (ja) 患者適応型のせん断波生成による超音波せん断波イメージング
KR102206496B1 (ko) 증가된 펄스 반복 인터벌을 갖는 초음파에 기반한 전단파 이미징
WO2018060820A1 (en) Ultrasonic shear wave imaging with background motion compensation
KR20220036337A (ko) 매질의 초음파 특성화 방법 및 시스템
JP5491671B2 (ja) 関心領域の粘弾性の平均値を測定するための方法

Legal Events

Date Code Title Description
A975 Report on accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A971005

Effective date: 20170828

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20170831

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

R150 Certificate of patent or registration of utility model

Ref document number: 6259953

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

LAPS Cancellation because of no payment of annual fees