CN102481137B - 用于剪切波分散振动测定的超声设备和方法 - Google Patents

用于剪切波分散振动测定的超声设备和方法 Download PDF

Info

Publication number
CN102481137B
CN102481137B CN201080029448.6A CN201080029448A CN102481137B CN 102481137 B CN102481137 B CN 102481137B CN 201080029448 A CN201080029448 A CN 201080029448A CN 102481137 B CN102481137 B CN 102481137B
Authority
CN
China
Prior art keywords
sampling
shearing wave
pulse
target position
trace pulse
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
CN201080029448.6A
Other languages
English (en)
Other versions
CN102481137A (zh
Inventor
H·谢
A·T·费尔南德斯
M·R·伯切
D·马克斯韦尔
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 Electronics 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 Electronics NV filed Critical Koninklijke Philips Electronics NV
Publication of CN102481137A publication Critical patent/CN102481137A/zh
Application granted granted Critical
Publication of CN102481137B publication Critical patent/CN102481137B/zh
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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0048Detecting, measuring or recording by applying mechanical forces or stimuli
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0048Detecting, measuring or recording by applying mechanical forces or stimuli
    • A61B5/0051Detecting, measuring or recording by applying mechanical forces or stimuli by applying vibrations
    • 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
    • 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

Landscapes

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

Abstract

执行剪切波分散振动测定(SDUV),从而使得在其推动脉冲(218)的单次请求之后,发出多个跟踪脉冲(222)以多于一次地采样在关联的单色剪切波(116)上的多个位置(120、148)中的每个,该采样以如下的采样方式:在分离的传送中扫描多个位置中的至少一个,并且利用多个跟踪脉冲中的一个脉冲并发地采样多个位置中的复数个。在补充的方面,通过考虑采样间延(156)来确定给定时刻的相位差,只要该确定依赖于不同时间进行的采样。

Description

用于剪切波分散振动测定的超声设备和方法
技术领域
本发明针对测量剪切波,更具体而言,针对通过剪切波分散振动测定(SDUV)来测量剪切波。
背景技术
活组织中的机械变化与病理学变化相关。因为在健康组织与病理组织之间,剪切弹性模量(硬度)和粘度可能显著变化。随着过去十年超声弹性成像的发展,很多临床研究已经显示:组织粘弹特性为医师提供了有用的信息以便于更好的癌症诊断和治疗评估。
剪切波分散振动测定(SDUV)为一种基于声辐射力的技术,该技术通过描述剪切波速度分散(即频率)的特征来测量组织剪切弹性和粘度。该技术的应用是无创测量肝硬度以区分肝纤维化和肝硬化的阶段。
用于医学成像的超声探查经常使用纵波。在体内组织中,超声以波形式传播。实际上,沿传播路径的所有粒子在适当位置前后振动,并且振动发生在传播方向上。振动产生压缩和稀疏。这些被模拟成正弦曲线的波峰和波谷。能量被传递给靶标,并且通过振荡的粒子移动返回。
与之相比,超声剪切波(或横波)特征为垂直于传播方向在适当位置的往复移动。振荡一方面产生波峰,另一方面产生波谷。
执行SDUV必须发出一系列聚焦的纵波推动脉冲。它们在焦点上建立了其传播方向与推动脉冲的传播方向垂直的剪切波。焦深已被选择,使得剪切波行进穿过感兴趣区域(ROI)。向ROI发出纵波跟踪脉冲以便在采样点评价剪切波的振幅。该测量被用于估计被采样位置处的剪切波的相位。为了采样另一位置,另一推动脉冲发往相同的推动焦点,接着向该位置发出跟踪脉冲。需要该第二个循环,因为两点间的相位差被用于确定弹性和粘度。
在两个循环中,在推动脉冲的发出与随后跟踪脉冲的发出之间的时间周期是相等的,这为测量给出了用于确定基于位置的相位差的重要手段。正是基于位置的相位差用作特定剪切波传播速度估计中的一因子,并且得到的若干波频上的速度值被用于推导弹性和粘度。
通常,如在Myao Clinic的两篇技术文章中描述的那样,在每个循环期间的SDUV采样来自单个位置,因为跟踪探头基本是不具有导引能力的活塞驱动的单元件换能器。该两篇技术文章为“Quantifying Elasticity andViscosity from Measurement of Shear Wave Speed Dispersion”,The Journal ofthe Acoustical Society of America,第115卷,第6期,第2781-2785页,2004,S.Chen、M.Fatemi和J.F.Greenleaf,以及“Error Estimates in Shear WaveSpeed and Tissue Material Properties in Shear Wave Dispersion UltrasoundVibrometry”,2007IEEE Ultrasonics Symposium,第664-667页,M.W.Urban,S.Chen和J.F.Greenleaf。
两篇文章都描述了实际的实验。
在两篇文章中都出现了使用两次测量之差作为相位差来计算波传播速度的公式,如公式(4)。Urban发表的出版物更详细示出了基础计算。
另一Mayo Clinic的出版物中,J.F.Greenleaf和S.Chen的名为“Ultrasound Vibrometry”、公开号2007/0038095的美国专利公开文本中还提出了在相同循环期间导引波束,以便:a)从两个或更多的位置上采样;并且b)进行两个或更多的相位测量。
对于所提出的方法,该专利公开文本中应用了完全相同的公式。
本发明的发明人已经注意到,上述现有技术似乎未引入实际的实验,未指定或建议与其提出的方法关联的任何时空采样方案,并且未公开或建议任何考虑到采样间延迟而对传播速度公式的改动。本发明人还注意到,未考虑当在相同循环期间从多个位置上采样时的这种延迟可能导致明显错误的结果。
更通常地,与使用SDUV检测辐射力引导的剪切波相关联的一个主要挑战在于相对低振幅的振动(轴向位移大约10μm)。采样测量在ROI中给定位置或当前位置上的剪切波的振幅。然而,系统电子噪声和诸如心脏运动或呼吸运动的患者运动将对剪切波位移估计贡献显著的噪声。
为了极大地减小严重的噪声影响,可以发射两个以上的推动/跟踪序列,以便探测在多个横向位置上的剪切位移的时空演变。可以将最小二乘拟合算法应用于多相位延迟估计以获得更加鲁棒的速度估计。
SDUV推动脉冲聚焦在组织中的相同激励位置上。来自多推动/跟踪序列的重复声暴露可能导致累积的局部组织加热。这在当连续激励间的时间延迟小于组织冷却时间时尤其明显。组织冷却时间将取决于被检查器官的灌注等。SDUV的热安全是研究者在该技术可以被应用到临床研究之前必须处理的关键问题之一。对于诊断超声,美国食品药品管理局(FDA)要求热指标(TI)小于6(组织加热小于6℃)。虽然来自单个SDUV推动的热增加大约为0.1℃到0.2℃,并且SDUV的暴露水平显著低于组织损害的阈值,但是SDUV的累积影响可能对组织和换能器加热带来风险。
因此,可能希望尽可能以提供鲁棒测量的方式实现,但是利用尽可能少的推动/跟踪序列,以避免不利的热影响。
发明内容
在一方面,一种新颖的剪切波分散振动测定(SDUV)方法包括,针对推动/跟踪序列,发出聚焦于感兴趣区域中的激励位置上的推动脉冲的单次请求,以建立剪切波;以及在针对所述推动/跟踪序列的所述推动脉冲的所述单次请求(instance)之后,发出跟踪脉冲以便多于一次地采样在所建立的剪切波的关联的单色剪切波上的多个靶标位置中的每个。其中,采样包括:(i)通过在所述多个跟踪脉冲的分离的传送中扫描所述多个靶标位置来获得采样,以及(ii)利用所述多个跟踪脉冲中的每个跟踪脉冲,通过并发地采样所述多个靶标位置中的复数个来获得采样,其中,所述采样提供用于确定被估计的剪切波传播速度的数据;所述方法还包括:通过考虑采样间延迟来确定给定时刻的相位差以用于确定所述被估计的剪切波传播速度,这依赖于在所述采样中不同时间获得的采样。
在补充的方面中,通过考虑采样间延迟来确定给定时刻的相位差。这在如果该确定依赖于不同时间进行的采样的情况下完成。
作为其他方面,该确定必须包括:使用两个不同时间的采样来计算该差,并且在所计算的差上加入基于采样间延迟的相位校正,该校正反映单色剪切波在时间上于进行所述两个采样之间发生的传播。
在替代的方面中,当从所述两个采样中的一个到另一个的扫描发生在与波传播相反的方向上时,所计算的差和该校正为符号相反的数字。
在另外的方面中,该校正与所述波的角频率成正比,并且与进行所述两个采样之间的延迟成正比。
在一个版本中,逐次传送地在相反方向上进行该采样,并且该考虑包括将来自两个相反方向的传送的相应测量进行组合以抵消该采样间延迟。
在再另外的方面中,在相同方向进行所述分离的传送。
在一些实施例中,逐个脉冲地顺次导引该跟踪脉冲。
在某些实施例中,该采样方式包括通过一个脉冲并发地采样多个位置。
在本发明的特别方面中,该脉冲被同时靶向到所述位置中的一个以上的位置。
在子方面中,跟踪脉冲被同时靶向所有位置。
在又另一方面中,响应于跟踪脉冲中的一个,并发地接收来自所述位置中对应的多个A线(A-line)。
在不同的方面中,当不存在推动脉冲振动时,发射参考脉冲,并且在计算该波的振幅中,将从该参考脉冲回波的数据与从跟踪脉冲回波的数据进行比较。
根据一些版本,一种用于执行剪切波分散振动测定(SDUV)的装置,所述装置包括:用于针对推动/跟踪序列,发出聚焦于感兴趣区域中的激励位置上的推动脉冲的单次请求,以建立剪切波的模块;以及用于在针对所述推动/跟踪序列的所述推动脉冲的所述单次请求之后,发出多个跟踪脉冲以便多于一次地采样在所建立的剪切波的关联的单色剪切波上的多个靶标位置中的每个的模块,其中,采样包括:(i)通过在所述多个跟踪脉冲的分离的传送中扫描所述多个靶标位置来获得采样,以及(ii)利用所述多个跟踪脉冲中的每个跟踪脉冲,通过并发地采样所述多个靶标位置中的复数个来获得采样,其中,所述采样提供用于确定被估计的剪切波传播速度的数据。所述装置还包括:用于通过考虑采样间延迟来确定给定时刻的相位差以用于确定所述被估计的剪切波传播速度的模块,确定所述被估计的剪切波传播速度依赖于在所述采样中不同时间获得的采样。
根据一些版本,一种制品包括其上编码有指令的机器可访问介质,该指令用于使进程执行上述方法。
同样地,在特别的实施例中,一种用于执行剪切波分散振动测定(SDUV)的计算机软件产品包括收录计算机程序的计算机可读介质,该计算机程序包括能够由处理器运行以执行上述方法的指令。
类似地,根据一个方面,一种被配置用于SDUV的超声设备,包括换能器装置,所述换能器装置被配置用于:针对推动/跟踪序列,发出聚焦于感兴趣区域中的激励位置上的推动脉冲的单次请求,以建立剪切波;以及在针对所述推动/跟踪序列的所述推动脉冲的所述单次请求之后,发出多个跟踪脉冲以便多于一次地采样在所建立的剪切波的关联的单色剪切波上的多个靶标位置中的每个,其中,采样包括:(i)通过在所述多个跟踪脉冲的分离的传送中扫描多个靶标位置来获得采样,以及(ii)利用所述多个跟踪脉冲中的每个跟踪脉冲,通过并发地采样所述多个靶标位置中的复数个来获得采样,其中,所述采样提供用于确定被估计的剪切波传播速度的数据;所述换能器装置还被配置用于:通过考虑采样间延迟来确定给定时刻的相位差以用于确定所述被估计的剪切波传播速度,这依赖于在所述采样中不同时间获得的采样。。
在又一其他实施例中,该设备被配置用于在没有用户干预的情况下自动基于图像深度和采样采集频率从每个采样单个A线接收模式切换到每个采样多个A线接收模式。
附图说明
参考下列附图,下文介绍了该新颖的SDUV推动/跟踪方案的细节。
图1为例示出根据本发明的连续两个SDUV采样的示意图;
图2为表示根据本发明的SDUV顺序跟踪的示例的时间线;
图3为根据本发明的对应于图2的时间线的流程图;
图4为根据本发明的两个可能的SDUV时空采样方案的概念图;
图5为示出了根据本发明的对于相应采样次序的可供选择的SDUV剪切波传播速度算法的流程图;
图6为表示SDUV并行跟踪的示例的时间线;
图7为表示根据本发明的SDUV的顺序跟踪并并行跟踪的组合的示例的时间线;以及
图8为根据本发明的在多线接收模式和单线接收模式之间自动切换类型的流程图和示意图组合。
具体实施方式
图1以说明性而非限制性示例的方式描绘了在两个连续SDUV采样之间的关系。示出了超声推动探头104,以及通过超声跟踪探头(未示出)分别在采样时间t1、t2发出的第一跟踪脉冲108和第二跟踪脉冲112。推动探头104发送一系列聚焦的推动脉冲(或“推动”)114以建立剪切波115,该剪切波115(或“单色剪切波”)的频率分量被表示为第一轨迹116。组织沿着y方向振动,并且剪切波沿着x方向传播。
可以分别将两个分离的探头用于推动和跟踪。然而这一类型的大体积配置对于临床应用是不实用的。
此外,在分离的单元(single-element)换能器的情况下,要求重复发射推动脉冲以便测量在不同横向位置处的剪切波相位延迟,从而导致重复的组织加热。
在另一方面,如果利用阵列换能器实施跟踪,则根据本文的提议缓和了对重复推动的需要。
此外,一个单阵列换能器能够同时实现推动和跟踪的功能。两用的单阵列换能器将使系统更加紧凑并且更容易控制。然而,该配置需要更多的硬件支持和软件支持,以便保证能够通过单个换能器产生声辐射力并且对得到的剪切波115进行成像。单个换能器可以为允许横向和垂向(elevation)推动和/或跟踪脉冲位置的单排阵列(1D换能器)或多排阵列(2D换能器)。
改进的2D阵列技术将允许以3D维发送跟踪脉冲;然而,下文的讨论涉及跟踪轴向-横向平面上的组织运动。在不失普遍性的情况下,相同的处置应用于轴向-垂向平面。
跟踪脉冲108、112被靶向到空间上分隔距离128、Δr的相应位置120、124。距离128通常在毫米的范围内,因为剪切波115随传播距离衰减。测量如在y轴136上表示的振幅132。假如采样频率满足奈奎斯特(Nyquist)阈值,则基于振幅132可以导出位置120处的相位140、Φ1。可以通过以足够频率在位置120上的多个采样传送(pass)来满足该奈奎斯特阈值。
如果采样时间t1、t2是同时的,则从跟踪脉冲112返回的数据(即回波数据)将表明与第一波形116上的采样点148相对应的相位142。
然而,如果对这两个位置120和124顺序发射跟踪脉冲,则采样时间t1、t2是非同时的。
实际上,在时间t2时刻,剪切波分量116已经向前传播,如由第二轨迹144表示的具有频率分量的波形。
在第二轨迹144上,对应的采样相位152、Φ2不同于在第一轨迹116上的相位140、Φ1,其中相位差为ΔΦ=Φ2–Φ1
如从图1中看到的那样,ΔΦ小于在提取同时跟踪脉冲的相位140、142之间的差别。
如从图1中进一步看到的那样,由基于采样间延迟的相位校正156、ωΔt使ΔΦ更小,“ω”表示单色剪切波116的角频率,“Δt”表示在相应采样时间t1、t2进行的两个采样之间的延迟。
另外,相位140、142之间的差别,与ω和Δr一起,定义了剪切波频率分量116的传播速度。
相位140、142之间的差别是通过将相位差ΔΦ加上基于采样间延迟的相位校正156ωΔt来计算得到的。
图2描绘了本文提出的SDUV顺序跟踪。在任何推动脉冲114振动感兴趣区域(ROI)之前,参考跟踪脉冲202、204、206、208发出,并被靶向在相应的采样位置x1、x2、x3、x4。位置x1、x2、x3和x4从用于推动的激励点向外放射性梯度排列(aligned radially),即在剪切波115的传播方向上。脉冲202、204、206、208每个都分别跟随有其各自的回波210、212、214、216。然后在位置x0(即激励点的位置)建立剪切波115。这通过以推动频率219(可能通常为大约100Hz)发射的推动脉冲218的初始系列(在图2中未示出)来完成。一旦波115已经建立,则可以发出推动脉冲218中的跟踪序列开始的一个(onet)。接着是逐个脉冲地顺次导引的跟踪脉冲222、224、226、228的传送220。发射跟踪脉冲222、224、226、228,每个跟随有其各自的回波230、232、234、236。以2KHz的相同位置跟踪频率238发出传送220。由此,脉冲218后跟随20个传送220。并且,2KHz的相同位置跟踪频率238允许对达到约1kHz(奈奎斯特阈值)的频率分量116进行相位测量。跟踪脉冲222、224、226、228中的每个都可以用于在奈奎斯特阈值的内在限制内同时提供对每个频率分量的测量。
利用常规B模式成像,即,将跟踪脉冲222、224、226、228在两个邻近的推动脉冲218之间的不同位置x1、x2、x3、x4上进行电子聚焦并且横向来回导引。因此顺序形成在不同位置x1、x2、x3、x4上的接收A线。应当限制在每个B模式帧内的传输线的数量,以使剪切波115的采样RPF(脉冲重复频率)足够高。如在图2中看到的,例如,对于每个位置x1、x2、x3、x4,相同位置跟踪频率238为2KHz(PRF2),同时传输A线频率242为8KHz(PRF3)。
有利地,仅仅需要单个推动/跟踪序列240来获得剪切波速度估计所需要的所有SDUV信息。如本文提出的,更快的数据采集是将SDUV技术从单点“虚拟活检”工具改进到可能的实时成像模态的关键要求。同时,当维持剪切波跟踪所观察的横向位置x1、x2、x3、x4、……的数量时,在组织和换能器中产生的热可以被降低。本发明的另一优点在于:其能够进一步减小在剪切波位移估计中的噪声,尤其避免大的时间-尺度噪声,例如在不同的推动/跟踪序列240期间组织经历的呼吸运动。
对于推动/跟踪序列的跟踪脉冲,本文定义了“该推动脉冲的单次请求”作为时间上紧在跟踪脉冲之前的推动脉冲。
图3为对应于图2中的时间线的流程图。焦点被设定在第一位置x1(步骤S304)。发射当前参考跟踪脉冲202(步骤S308)。跟着是其回波210(步骤S312)。如果存在更多待测位置(即x2、x3、x4)(步骤S316),则焦点被导引到下个位置(步骤S320),并且进程返回到参考跟踪脉冲发射步骤S308,其中将下个位置作为当前位置。另一方面,如果不存在其他待测位置(步骤S316),则连续发射推动脉冲218以建立剪切波115(步骤S324)。焦点然后被设定到第一位置x1(步骤S328)。发射当前跟踪参考脉冲222(步骤S330),并且返回其回波230(步骤S332)。如果存在更多的跟踪脉冲(即224、226、228)在当前传送220中发出(步骤S333),则焦点被设定到对应的下个位置(步骤S334),并且进程返回到跟踪脉冲发射步骤S330,其中将下个位置作为当前位置。否则,如果不存在下个待发射的跟踪脉冲,则询问是否将在位置x1、x2、x3、x4进行下次传送(步骤S336)。如果将进行下次传送,则将下次传送作为当前传送(S340),并且进程返回到下次传送220的开端,到步骤S328。然而,如果(在当前推动/跟踪序列240中)不存在待进行的下次传送220,则完成采样。在本文提出的技术中,单个推动/跟踪序列240足够提供剪切波速度的SDUV确定(在实际实施例中,组织弹性和粘度的计算)所需要的所有采样。
如本文所提出的,两个可能的SDUV时空采样方案402、404被描绘在图4中。
示出了对于统一方向采样方案402的四次传送220。每次都在三个位置x1、x2、x3上。每次都在图4中以左到右的方向表示。由411到422指示的数字1到12对应于采样发生的次序。在相同方向上进行分离的传送220。
同样地,示出了在三个位置x1、x2、x3上的四次传送220的交替方向采样方案404。然而,与统一方向采样方案402相比,交替方向采样方案404交替每次传送220的扫描方向。由431到442指示的数字1到12对应于采样发生的次序。如下文将更详细地讨论的那样,在相反方向上的采样是加入基于采样间延迟的相位校正156的替代选择。
在位置x2上的相同位置跟踪频率238在两种采样方案402、404中相同。对于交替方向采样方案404中的偏离中心位置x1、x3,采样在时间上不规则。然而,对于不规则采样保持奈奎斯特阈值。平均采样频率服从该阈值。由于平均采样频率对于两种方案402、404的所有位置x1、x2、x3都是相同的,因此任一方案可以避免混叠的方式实施。
图5说明了在不同时间t1、t2进行采样的相应次序的交替SDUV剪切波传播速度算法500。
针对在与剪切波115传播的方向相同的方向上的扫描,应用正向扫描公式510。公式510的推导如下:
cs(ω)=λf,其中cs为传播速度,λ为波长,并且f为波的频率。
并且,ω=2πf,并且Δr=λΔΦ/2π,其中ΔΦ表示对于给定时刻在同时进行的两个采样之间的相位差。
假设所有参数是正的。
代入值得出 c s ( ω ) = ωΔr Δφ .
如上文参考图1解释的那样,为处理在时间t1、t2进行的两个不同时间采样之间的采样间延迟,将基于采样间延迟的相位校正156ωΔt加入ΔΦ。该校正156与波115的角频率ω成正比,特别是与当前被测量速度的波分量116的角频率成正比,并且与进行该两个采样之间的延迟Δt成正比。
针对在与剪切波115传播的方向相反的方向上的扫描,应用反向扫描公式520。在反向扫描公式520中,基于采样间延迟的相位校正156ωΔt的符号与所计算的给定时刻(相位)差ΔΦ的符号相反。使该校正156为负的原因是基于类同于上文参考图1解释的推理。对于单色剪切波116,基于采样间延迟的相位校正156反映了时间上在分别于时间t1、t2进行的两个时间不同的采样之间发生的传播。
如果在一次传送220上的扫描与在另一(例如随后的)传送220上的扫描的方向相反,则可以将在相同的两个位置x1、x2上各自的ΔΦ测量进行平均来抵消采样间延迟。这假设了在两个位置x1、x2之间的采样间延迟在两次传送上是相等的。如果假设成立,公式530的分母表明:基于采样间延迟的相位校正156(如果确定了)将由单个扫描方向公式510、520的分母的平均进行抵消。由此,不需要计算该校正156,并且使用公式530不需要校正156。
图6描绘了如与图2中的顺序跟踪相反的SDUV并行跟踪,该描绘表示一个可能的实施例。在并行跟踪中,跟踪脉冲同时靶向多个位置中一个以上的位置,并且实际上靶向所有的位置。
为了监测由推动脉冲218建立的剪切波传播,在每个检测点或位置x1、x2、x3、x4上的相同位置跟踪频率238应当被设定为足够高以便满足与处于分析下的最高谐振频率相关的奈奎斯特限制。遵循上文所述并且如图6中所示的示例情形,通过2KHz的相同位置跟踪频率238来跟踪每个检测点,同时利用8KHz的传输A线频率242来顺序扫描四个检测点或位置x1、x2、x3、x4。如果期望更多的检测点,则将要求超声系统以甚至更高的RPF传输,该甚至更高的RPF在“顺序跟踪”模式下是难以实现或不可能实现的(受到最大成像深度以及从最大成像深度返回的超声传播时间物理上限制)。
为了克服该问题,本文提出了一种称为“并行跟踪”的第二方法,并且该第二方法在下文的示例性实施例中进行描绘。
在平行接收B模式中,宽的传输波束每个被微弱聚焦。在相应推动218之间,它们被重复发送到跨越不同位置x1、x2、x3、x4的相同区域。对于给定的单个检测传输,由多线波束形成器使用通过阵列立即接收的反向散射信号来平行-波束形成多个接收A线。换言之,利用每个跟踪脉冲并发地采样多个位置中的复数个。
发送参考跟踪传输602。对于相应的位置x1、x2、x3、x4形成多个平行接收参考A线604、606、608、610。在激励点x0的推动脉冲218之后,将若干跟踪传输脉冲612发到对应的位置x1、x2、x3、x4,每个跟踪传输脉冲跟随有相应的平行接收A线614、616、618、620。在平行推动/跟踪序列624期间以2kHz的并行跟踪频率622重复跟踪传输脉冲612及其相应的平行接收A线614、616、618、620。
如果期望,这里可以增加更多的接收A线,以便在不降低跟踪LPF的情况下增加空间采样精度,只要系统支持更宽的传输以及高次序多线接收波束形成器。与上文在图2和图3中描述的顺序方法类似,并行跟踪方法并不要求在剪切波激励原点重复发射聚焦超声。其使剪切波速度估计能够经由单个推动/跟踪序列来实施。
在并行跟踪方法中,因为覆盖不同横向位置x1、x2、x3、x4的检测脉冲在单个发射中发送,所以并发采样在不同横向位置的位移估计132。在该情况下,不需要做出任何补偿。应当使用以下公式简单估计剪切波速:
c s ( ω ) = ωΔr Δφ . (公式540)
当然,利用微弱焦点传输更宽的波束612将降低在焦点上的图像空间分辨率。在平面波传输的极端情况下,并行跟踪可以产生比顺序跟踪略微低的SNR的位移估计。如果SNR是非常关键的问题或者超声系统仅装配有低次序多线波束形成器,则“顺序跟踪”方法和“并行跟踪”方法可以组合成为“混合跟踪”方法。
图7为表示SDUV的顺序和并行跟踪的组合的示例的时间线。
除了一些位置之外,顺序和并行跟踪与上述并行跟踪类似,在此同时采样x1和x2(如通过跟踪传输脉冲702),此外相对于位置x3、x4顺序采样位置x1、x2。
在该情形中,2X多线波束形成器满足并行跟踪的需要。独立点的跟踪PRF为2KHz(PRF2),同时超声系统传输PRF为4KHz(PRF3)。当计算剪切波速时,如果分析点x1和x2,应当使用上面示出的公式540。与之相比,如果分析点x1和x3,应当使用适当的公式510、520、530。
图8为根据本发明的在多线接收模式和单线接收模式之间自动切换的类型的流程图和示意图组合。
超声探头804传输跟踪脉冲808、812到ROI816中最大图像深度817处的相应位置,这里为器官或血管。多线接收模式820中的传输如图6所示,或者单线接收模式824中的传输如图2所示。
监测进程828校验最大图像深度817(步骤S832)、传输A线频率242(步骤S836)以及当前模式836(步骤840)。如果最大图像深度817和传输A线频率242足够高,则穿过组织的声速可以具有以下限制:要求从单线接收模式824切换到多线接收模式820。同样地,可以依据发信号的节奏是否能被适应来做出从多线接收模式到单线接收模式的转变。如果指示切换(步骤S844),则在进程为了重校验在延迟之后返回到步骤S832、S836、S840之前,利用新的模式848交换当前模式836(步骤S852)。
如此执行剪切波分散振动测定(SDUV)以使在它们的推动脉冲的单次请求之后,发出多个跟踪脉冲以便多于一次地采样在关联的单色剪切波上的多个位置中的每个,该采样以如下的采样方式:在分离的传送中扫描多个位置中的至少一个,并且利用多个跟踪脉冲中的一个脉冲并发地采样多个位置中的复数个。在补充的方面,通过考虑采样间延迟来确定给定时刻的相位差,只要该确定依赖于不同时间进行的采样。
一种创新的SDUV推动/跟踪方案能够通过限制声辐射力激励218的数量来降低组织加热。
应当注意的是,上述实施例进行举例说明而不是限制本发明,并且本领域技术人员将能够在不脱离权利要求书的保护范围的情况下设计很多替代的实施例。例如,可以在各种实施例中采样比在图中示出的更多或更少的位置。并且,参考跟踪脉冲可以是已往的,并且可以由对于一个位置的重复A线之间的互相关来替代,以便确定波的振幅。在权利要求书中,任何置于括号内的附图标记不应解释为对权利要求保护范围的限制。动词“包括”及其词形变化的使用并不排除存在权利要求未记载的那些部件或步骤。部件之前的不定冠词“一个”或“一种”并不排除存在多个或多种这样的部件。本发明可以实现为包括若干分立部件的硬件,并且实现为具有计算机可读介质的适当编程的计算机。在相互不同的从属权利要求中记载某些手段并不代表不能使用这些手段的组合以获得优点。

Claims (15)

1.一种剪切波分散振动测定(SDUV)方法,包括:
针对推动/跟踪序列,发出聚焦于感兴趣区域中的激励位置上的推动脉冲的单次请求,以建立剪切波;以及
在针对所述推动/跟踪序列的所述推动脉冲(218)的所述单次请求之后,发出多个跟踪脉冲(222)以便多于一次地采样在所建立的剪切波的关联的单色剪切波(116)上的多个靶标位置(120)中的每个,其中,采样包括:(i)通过在所述多个跟踪脉冲的分离的传送中扫描所述多个靶标位置来获得采样,以及(ii)利用所述多个跟踪脉冲中的每个跟踪脉冲(612),通过并发地采样所述多个靶标位置中的复数个来获得采样,其中,所述采样提供用于确定被估计的剪切波传播速度的数据;所述方法还包括:
通过考虑采样间延迟来确定给定时刻的相位差以用于确定所述被估计的剪切波传播速度,这依赖于在所述采样中不同时间获得的采样。
2.根据权利要求1所述的方法,其中,对所述相位差的所述确定还包括:
使用两个不同时间的采样计算所述相位差;以及
将所计算的相位差加上基于采样间延迟的相位校正(156),所述基于采样间延迟的相位校正反映所述单色剪切波在时间上于获得所述两个采样之间发生的传播。
3.根据权利要求2所述的方法,其中,当从所述两个采样中的一个靶标位置到另一个靶标位置的扫描发生在与所述剪切波传播的方向相反的方向时,所计算的相位差与所述基于采样间延迟的相位校正是符号相反的数字(520)。
4.根据权利要求2所述的方法,其中,所述基于采样间延迟的相位校正(i)与所述剪切波的角频率成正比,并且(ii)与获得所述两个采样之间的延迟成正比。
5.根据权利要求1所述的方法,其中,所述采样方式逐次传送地以相反方向发生,并且所述考虑包括将来自所述传送中相反方向的两个的相应测量进行组合,以便抵消所述采样间延迟(530)。
6.根据权利要求1所述的方法,其中,逐个脉冲地顺次导引所述跟踪脉冲(108、112)。
7.根据权利要求6所述的方法,其中,所述采样方式利用所述多个跟踪脉冲中的单个跟踪脉冲(702)来并发地采样所述多个靶标位置中的复数个。
8.根据权利要求1所述的方法,其中,所述跟踪脉冲(702)被同时靶向到所述多个靶标位置中的一个以上的位置。
9.根据权利要求1所述的方法,还包括响应于所述多个跟踪脉冲中的一个,并发地从所述多个靶标位置中的对应多个位置接收多个接收A线(614)。
10.根据权利要求1所述的方法,还包括:
当不存在推动脉冲振动时,发射参考脉冲(202);以及,
在计算所述剪切波的振幅(132)中,将从所述参考脉冲回波的数据与从所述多个跟踪脉冲中的一跟踪脉冲回波的数据进行比较。
11.一种用于执行剪切波分散振动测定(SDUV)的装置,所述装置包括:
用于针对推动/跟踪序列,发出聚焦于感兴趣区域中的激励位置上的推动脉冲的单次请求,以建立剪切波的模块;以及
用于在针对所述推动/跟踪序列的所述推动脉冲的所述单次请求之后,发出多个跟踪脉冲以便多于一次地采样在所建立的剪切波的关联的单色剪切波上的多个靶标位置中的每个的模块,其中,采样包括:(i)通过在所述多个跟踪脉冲的分离的传送中扫描所述多个靶标位置(411-416)来获得采样,以及(ii)利用所述多个跟踪脉冲中的每个跟踪脉冲,通过并发地采样所述多个靶标位置中的复数个来获得采样,其中,所述采样提供用于确定被估计的剪切波传播速度的数据;所述装置还包括:
用于通过考虑采样间延迟来确定给定时刻的相位差以用于确定所述被估计的剪切波传播速度的模块,确定所述被估计的剪切波传播速度依赖于在所述采样中不同时间获得的采样。
12.一种被配置用于剪切波分散振动测定(SDUV)的超声设备,包括:
换能器装置(104),所述换能器装置(104)被配置用于:针对推动/跟踪序列,发出聚焦于感兴趣区域中的激励位置上的推动脉冲的单次请求,以建立剪切波;以及
在针对所述推动/跟踪序列的所述推动脉冲的所述单次请求之后,发出多个跟踪脉冲以便多于一次地采样在所建立的剪切波的关联的单色剪切波上的多个靶标位置中的每个,其中,采样包括:(i)通过在所述多个跟踪脉冲的分离的传送中扫描多个靶标位置来获得采样,以及(ii)利用所述多个跟踪脉冲中的每个跟踪脉冲,通过并发地采样所述多个靶标位置中的复数个来获得采样,其中,所述采样提供用于确定被估计的剪切波传播速度的数据;所述换能器装置还被配置用于:
通过考虑采样间延迟来确定给定时刻的相位差以用于确定所述被估计的剪切波传播速度,这依赖于在所述采样中不同时间获得的采样。
13.根据权利要求12所述的超声设备,被进一步配置用于:
使用在所述采样方式中不同时间进行的两个采样来计算相位差;以及
将所计算的相位差加上基于采样间延迟的相位校正,所述校正反映所述单色剪切波在时间上于获得所述两个采样之间发生的传播(156)。
14.根据权利要求12所述的超声设备,其中,所述设备被进一步配置用于逐次传送在相反方向(431-436)上进行采样,并且所述设备被进一步配置用于通过将来自所述传送中相反方向的两个的相应测量进行组合以抵消所述采样间延迟来考虑所述采样间延迟。
15.根据权利要求12所述的超声设备,被进一步配置用于,在没有用户干预的情况下,自动基于图像深度(817)以及传输A线频率(242)从每采样单个A线接收模式切换到每采样多个A线接收模式。
CN201080029448.6A 2009-06-30 2010-06-23 用于剪切波分散振动测定的超声设备和方法 Active CN102481137B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US22183109P 2009-06-30 2009-06-30
US61/221,831 2009-06-30
PCT/IB2010/052852 WO2011001333A1 (en) 2009-06-30 2010-06-23 Push/tracking sequences for shear wave dispersion vibrometry

Publications (2)

Publication Number Publication Date
CN102481137A CN102481137A (zh) 2012-05-30
CN102481137B true CN102481137B (zh) 2015-03-18

Family

ID=43029704

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201080029448.6A Active CN102481137B (zh) 2009-06-30 2010-06-23 用于剪切波分散振动测定的超声设备和方法

Country Status (7)

Country Link
US (2) US9125547B2 (zh)
EP (1) EP2448494B1 (zh)
JP (1) JP5877785B2 (zh)
CN (1) CN102481137B (zh)
BR (1) BRPI1010214B8 (zh)
RU (1) RU2559910C2 (zh)
WO (1) WO2011001333A1 (zh)

Families Citing this family (31)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101068039B1 (ko) * 2009-04-29 2011-09-28 알피니언메디칼시스템 주식회사 횡탄성파 생성 방법, 횡탄성파를 이용한 이미지 획득 방법 및 장치
WO2011007278A2 (en) * 2009-07-17 2011-01-20 Koninklijke Philips Electronics, N.V. Spatially-fine shear wave dispersion ultrasound vibrometry sampling
JP5886318B2 (ja) 2010-12-22 2016-03-16 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 質量中心を用いたせん断波速度の推定
KR101929198B1 (ko) * 2011-02-25 2018-12-14 메이오 파운데이션 포 메디칼 에쥬케이션 앤드 리써치 비집속식 초음파에 의한 초음파 바이브로메트리
WO2013073304A1 (ja) * 2011-11-17 2013-05-23 日立アロカメディカル株式会社 超音波診断装置及び超音波画像生成方法
CN103169493A (zh) * 2011-12-20 2013-06-26 通用电气公司 超声探针引导装置、方法及超声系统
CN103300890B (zh) 2012-03-16 2016-06-08 通用电气公司 用于测量组织机械特性的系统及方法
JP2014033943A (ja) * 2012-08-10 2014-02-24 Sony Corp 振動検出装置、振動検出方法、振動検出システムおよびプログラム
CN103800038B (zh) * 2012-11-12 2016-09-21 通用电气公司 改善的系统和装置以用于确定目标组织的机械特性
CN103908289B (zh) * 2012-12-31 2019-11-12 Ge医疗系统环球技术有限公司 消除剪切波中的背景噪声的方法和相应的超声成像系统
KR101874101B1 (ko) * 2013-08-12 2018-07-03 삼성전자주식회사 탄성 영상을 생성하는 방법 및 초음파 진단 장치
JP6288996B2 (ja) * 2013-09-11 2018-03-07 キヤノンメディカルシステムズ株式会社 超音波診断装置及び超音波イメージングプログラム
US10143442B2 (en) 2013-10-24 2018-12-04 Ge Medical Systems Global Technology, Llc Ultrasonic diagnosis apparatus
KR101654674B1 (ko) * 2013-11-28 2016-09-06 삼성전자주식회사 탄성 영상 제공 방법 및 이를 위한 초음파 장치
JP2015128554A (ja) * 2014-01-09 2015-07-16 日立アロカメディカル株式会社 超音波診断装置
WO2015151972A1 (ja) * 2014-04-02 2015-10-08 国立大学法人群馬大学 超音波映像システム
US10588604B2 (en) 2014-05-16 2020-03-17 Koninklijke Philips N.V. Autocorrelation guided cross-correlation in ultrasound shear wave elastography
CN110507360B (zh) * 2014-08-28 2022-06-03 深圳迈瑞生物医疗电子股份有限公司 剪切波成像方法及系统
CN106572838B (zh) * 2014-09-03 2019-09-06 深圳迈瑞生物医疗电子股份有限公司 弹性测量检测方法及系统
JP2017533031A (ja) * 2014-10-29 2017-11-09 メイヨ フォンデーシヨン フォー メディカル エジュケーション アンド リサーチ 超音波トランスデューサの連続振動による超音波エラストグラフィのための方法
CN104458081B (zh) * 2014-12-03 2016-07-20 中国矿业大学 一种基于超声表面波的应力测量主应力分离的方法及装置
RU2702090C2 (ru) * 2015-03-31 2019-10-04 Конинклейке Филипс Н.В. Калибровка ультразвукового, основанного на эластичности, отображения границы очага поражения
WO2016196631A1 (en) * 2015-06-01 2016-12-08 Duke University Methods, systems and computer program products for single track location shear wave elasticity imaging
US12023199B2 (en) * 2015-10-08 2024-07-02 Mayo Foundation For Medical Education And Research Systems and methods for ultrasound elastography with continuous transducer vibration
US11071524B2 (en) 2015-12-04 2021-07-27 Canon Medical Systems Corporation Analyzing apparatus
JP6987496B2 (ja) * 2015-12-04 2022-01-05 キヤノンメディカルシステムズ株式会社 解析装置
JP6672809B2 (ja) * 2016-01-13 2020-03-25 コニカミノルタ株式会社 超音波診断装置、及び超音波信号処理方法
JP6601320B2 (ja) * 2016-06-16 2019-11-06 コニカミノルタ株式会社 超音波診断装置、及び超音波診断装置の制御方法
US10646202B2 (en) 2017-01-23 2020-05-12 Siemens Medical Solutions Usa, Inc. Sheer speed imaging using coherence
CN112244813B (zh) * 2020-10-22 2022-08-26 无锡鸣石峻致医疗科技有限公司 一种低场核磁共振弹性测量方法及系统
KR102688130B1 (ko) * 2021-12-13 2024-07-25 알피니언메디칼시스템 주식회사 전단 파 파라미터 측정 방법 및 초음파 장치

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4489680A (en) 1984-01-23 1984-12-25 Borg-Warner Corporation Engine temperature control system
US5099848A (en) * 1990-11-02 1992-03-31 University Of Rochester Method and apparatus for breast imaging and tumor detection using modal vibration analysis
JPH06273397A (ja) * 1993-03-18 1994-09-30 Olympus Optical Co Ltd 超音波測定装置
US8041415B2 (en) * 2002-07-31 2011-10-18 Tsuyoshi Shiina Ultrasonic diagnosis system and strain distribution display method
FR2843290B1 (fr) * 2002-08-08 2005-06-24 Echosens Dispositif et procede pour la mesure de l'elasticite d'un organe humain ou animal
CN1809399B (zh) * 2003-04-17 2010-12-22 布赖汉姆妇女医院 剪切式治疗超声波
US7785259B2 (en) * 2003-10-03 2010-08-31 Mayo Foundation For Medical Education And Research Detection of motion in vibro-acoustography
US7753847B2 (en) 2003-10-03 2010-07-13 Mayo Foundation For Medical Education And Research Ultrasound vibrometry
FR2869521B1 (fr) * 2004-05-03 2007-02-02 Echosens Sa Dispositif pour la mesure de l'elasticite d'un organe humain ou animal au travers d'un conduit
RU2292530C1 (ru) * 2005-04-14 2007-01-27 Государственное образовательное учреждение высшего профессионального образования "Алтайский государственный технический университет им. И.И. Ползунова" (АлтГТУ) Способ измерения амплитуды колебаний
US8858441B2 (en) * 2005-05-12 2014-10-14 The Trustees Of Columbia University In The City Of New York System and method for electromechanical wave imaging of body structures
CA2685886C (en) * 2007-05-16 2016-02-23 Super Sonic Imagine Method and device for measuring a mean value of visco-elasticity of a region of interest

Also Published As

Publication number Publication date
RU2559910C2 (ru) 2015-08-20
US9125547B2 (en) 2015-09-08
EP2448494B1 (en) 2019-11-20
US10639001B2 (en) 2020-05-05
BRPI1010214B8 (pt) 2021-06-22
BRPI1010214A2 (pt) 2017-01-10
JP5877785B2 (ja) 2016-03-08
BRPI1010214B1 (pt) 2021-01-26
RU2012103016A (ru) 2013-08-10
JP2012531937A (ja) 2012-12-13
US20150335313A1 (en) 2015-11-26
EP2448494A1 (en) 2012-05-09
US20120123262A1 (en) 2012-05-17
CN102481137A (zh) 2012-05-30
WO2011001333A1 (en) 2011-01-06

Similar Documents

Publication Publication Date Title
CN102481137B (zh) 用于剪切波分散振动测定的超声设备和方法
CN102469980B (zh) 空间上精细的剪切波分散超声振动测定采样
KR101398948B1 (ko) 진폭―위상 변조된 초음파를 이용한 점탄성 측정
JP6129744B2 (ja) バックグランド動き効果に関する音響放射力効果の測定の調整
EP2654552B1 (en) Shear wave velocity estimation using center of mass
EP2504716B1 (en) Ultrasonic shear wave imaging with focused scanline beamforming
JP6353929B2 (ja) 超音波診断装置、及び弾性評価方法
CN109077754B (zh) 一种测量组织力学特性参数的方法及设备
WO2007138751A1 (ja) 超音波診断装置、医用画像処理装置及び医用画像処理プログラム
CN111885965A (zh) 用于三维中的剪切波成像的超声系统
JP5588924B2 (ja) 超音波診断装置
JPWO2019058753A1 (ja) 超音波診断装置および超音波診断装置の制御方法
KR102206496B1 (ko) 증가된 펄스 반복 인터벌을 갖는 초음파에 기반한 전단파 이미징
JP2017070317A (ja) 超音波診断装置及び脈波計測方法
JP4691185B2 (ja) 超音波画像表示装置
JP4581596B2 (ja) 超音波診断装置
Zhao et al. Robust shear wave speed measurement using comb-push ultrasound radiation force
JP5683210B2 (ja) 3次元pwスペクトル超音波診断装置および3次元pwスペクトルデータ生成方法
Dumont et al. Feasability of a ARFI/B-mode/Doppler system for real-time, freehand scanning of the cardiovascular system
WO2023278445A9 (en) System and method for non-invasive determination of pressure in a biological compartment
JP2008183118A (ja) 超音波診断装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant