CN103462643A - 一种剪切波速度测量方法、装置和系统 - Google Patents

一种剪切波速度测量方法、装置和系统 Download PDF

Info

Publication number
CN103462643A
CN103462643A CN2013104547253A CN201310454725A CN103462643A CN 103462643 A CN103462643 A CN 103462643A CN 2013104547253 A CN2013104547253 A CN 2013104547253A CN 201310454725 A CN201310454725 A CN 201310454725A CN 103462643 A CN103462643 A CN 103462643A
Authority
CN
China
Prior art keywords
wave
soft tissue
velocity
scan line
line position
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2013104547253A
Other languages
English (en)
Other versions
CN103462643B (zh
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.)
Sonoscape Medical Corp
Original Assignee
Sonoscape Co Ltd
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 Sonoscape Co Ltd filed Critical Sonoscape Co Ltd
Priority to CN201310454725.3A priority Critical patent/CN103462643B/zh
Publication of CN103462643A publication Critical patent/CN103462643A/zh
Application granted granted Critical
Publication of CN103462643B publication Critical patent/CN103462643B/zh
Active 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/48Diagnostic techniques
    • A61B8/485Diagnostic techniques involving measuring strain or elastic properties

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本发明公开了一种剪切波速度测量方法、装置和系统,该方法包括:获取待检测的软组织区域针对多次发射的超声脉冲检测波束反射回来的多个回波信息,其中,待检测的软组织区域内存在剪切波的传播;根据所述多个回波信息,分别确定所述剪切波在所述软组织区域中预设的任意两条扫描线之间的波速;依据每个所述波速对应的两条扫描线的距离,确定每个所述波速的权重,其中,所述波速所对应的两条扫描线的距离越大,所述波速的权重越大;根据每个所述波速的权重,对所有所述波速进行加权平均,得到所述剪切波在所述软组织区域内的传播速度。该方法可以提高软组织内剪切波传播速度测量的准确度和可靠性。

Description

一种剪切波速度测量方法、装置和系统
技术领域
本发明涉及软组织病变分析技术领域,更具体的说是涉及一种剪切波速度测量方法、装置和系统。
背景技术
基于声辐射力的剪切波弹性成像技术是一种评估组织弹性的超声弹性成像技术,其被广泛应用于人体组织病变分析诊断。基于声辐射力的剪切波弹性成像技术的原理为:由探头向生物体的软组织发射超声脉冲后,在剪切应力的作用下,特定区域内的软组织会产生振动,从而产生剪切波,而由于生物体软组织的弹性与剪切波的速度有一定的关联关系,因此可以通过检测剪切波的速度分析软组织的弹性,进而确定是否存在软组织的病变。
然而由于剪切波在软组织中传播的复杂度,给剪切波速度的测量造成了很大影响,使得测得剪切波在软组织中的传播速度存在较大误差,进而影响了正常的组织病变分析。
发明内容
有鉴于此,本发明提供一种剪切波速度测量方法、装置和系统,以提高软组织内剪切波传播速度测量的准确度和可靠性。
为实现上述目的,本发明提供如下技术方案:一种剪切波速度测量方法,包括:
获取待检测的软组织区域针对多次发射的超声脉冲检测波束反射回来的多个回波信息,其中,所述待检测的软组织区域内存在剪切波的传播;
根据所述多个回波信息,分别确定所述剪切波在所述软组织区域中预设的任意两条扫描线之间的波速;
依据每个所述波速对应的两条扫描线的距离,确定每个所述波速的权重,其中,所述波速所对应的两条扫描线的距离越大,所述波速的权重越大;
根据每个所述波速的权重,对所有所述波速进行加权平均,得到所述剪切波在所述软组织区域内的传播速度。
优选地,所述依据每个所述波速对应的两条扫描线的距离,确定每个所述波速的权重,包括:
对于任意一个所述波速,计算所述波速对应的两条扫描线的距离与第一预设值的乘积,将所述乘积作为所述波速的权重,其中,所述第一预设值为大于零的任意值。
优选地,所述根据所述多个回波信息,分别确定所述剪切波在所述软组织区域中预设的任意两条扫描线之间的波速,包括:
根据所述多个回波信息,并结合波束合成技术,确定所述软组织区域中,预设的每个扫描线位置针对所述多次发射的超声脉冲检测波束反射的多个第一回波信号;
根据所述每个扫描线位置所对应的多个第一回波信号,分别确定所述剪切波在所述软组织区域中预设的任意两条扫描线之间的波速。
优选地,所述根据所述多个回波信息,分别确定所述剪切波在所述软组织区域中预设的任意两条扫描线之间的波速,包括:
根据所述多个回波信息,确定所述剪切波到达所述软组织区域内预设的任意两条扫描线的第一时间差;
对于任意两条扫描线,依据所述两条扫描线的第一距离以及所述剪切波到达所述两条扫描线的第一时间差,计算所述剪切波在所述两条扫描线之间的波速。
优选地,所述根据所述多个回波信息,确定所述剪切波到达所述软组织区域内预设的任意两条扫描线的第一时间差,包括:
根据所述多个回波信息,并结合波束合成技术,确定所述软组织区域中,预设的每个扫描线位置处针对所述多次发射的超声脉冲检测波束反射回来的多个第一回波信号;
对于任意一个扫描线位置,根据预先存储的所述扫描线位置的静态组织回波信号以及所述多个第一回波信号,分别确定所述扫描线位置处的软组织在多个不同时刻的位移量,其中,所述在多个不同时刻位移量分别为所述软组织在反射回所述多个第一回波信号的多个时刻所对应的位移量;
根据每个所述扫描线位置处的软组织在多个不同时刻的位移量,确定每个扫描线位置处对应的所述剪切波的波形;
对于任意两个扫描线位置,根据所述两个扫描线位置处对应的所述剪切波的波形,确定所述剪切波到达所述两个扫描线位置的第一时间差。
优选地,所述根据所述两个扫描线位置处对应的所述剪切波的波形,确定所述剪切波到达所述两个扫描线位置的第一时间差,包括:
对所述两个扫描线位置处对应的所述剪切波的波形进行互相关运算,则互相关最大值所对应的时间延迟量为所述第一时间差;
或者,分别将每个扫描线位置处对应的所述剪切波的波形中,位移最大值所对应的时刻作为所述剪切波到达所述扫描线位置的时刻,根据剪切波到达所述两个扫描线位置的时间,确定出所述剪切波到达所述两个扫描线位置的第一时间差。
优选地,所述获取待检测的软组织区域针对发射的超声脉冲检测波束反射回来的回波信息之前,还包括:
控制第一超声阵列探头向所述软组织区域对应的聚焦点位置发射超声脉冲激励波束,以在所述软组织区域内形成剪切波;
控制第二超声阵列探头按照预设周期,向所述软组织区域发射多次超声脉冲检测波束,并接收每次反射回来的所述回波信息。
另一方面,本发明还提供了一种剪切波速度测量装置,包括:
信息获取单元,用于获取待检测的软组织区域针对多次发射的超声脉冲检测波束反射回来的多个回波信息,其中,所述待检测的软组织区域内存在剪切波的传播;
分段波速确定单元,用于根据所述多个回波信息,分别确定所述剪切波在所述软组织区域中预设的任意两条扫描线之间的波速;
权重设置单元,用于依据每个所述波速对应的两条扫描线的距离,确定每个所述波速的权重,其中,所述波速所对应的两条扫描线的距离越大,所述波速的权重越大;
剪切波速度确定单元,用于根据每个所述波速的权重,对所有所述波速进行加权平均,得到所述剪切波在所述软组织区域内的传播速度。
优选地,所述权重设置单元,包括:
第一权重设置单元,用于对于任意一个波速,计算所述波速对应的两条扫描线的距离与第一预设值的乘积,将所述乘积作为所述波速的权重,其中,所述第一预设值为大于零的任意值。
优选地,所述分段波速确定单元,包括:
时间差确定单元,用于根据所述多个回波信息,确定所述剪切波到达所述软组织区域内预设的任意两条扫描线的第一时间差;
分段波速确定子单元,用于对于任意两条扫描线,依据所述两条扫描线的第一距离以及所述剪切波到达所述两条扫描线的第一时间差,计算所述剪切波在所述两条扫描线之间的波速。
优选地,所述时间差确定单元,包括:
信号合成单元,用于根据所述多个回波信息,并结合波束合成技术,确定所述软组织区域中,预设的每个扫描线位置针对所述多次发射的超声脉冲检测波束反射回来的多个第一回波信号;
位移量确定单元,用于对于任意一个扫描线位置,根据预先存储的所述扫描线位置的静态组织回波信号以及所述多个第一回波信号,分别确定所述扫描线位置处的软组织在多个不同时刻的位移量,其中,所述在多个不同时刻位移量分别为所述软组织在反射回所述多个第一回波信号的多个时刻所对应的位移量;
波形确定单元,用于根据每个所述扫描线位置处的软组织在多个不同时刻的位移量,确定每个扫描线位置处对应的所述剪切波的波形;
时间差确定子单元,用于对于任意两个扫描线位置,根据所述两个扫描线位置处对应的所述剪切波的波形,确定所述剪切波到达所述两个扫描线位置的第一时间差。
优选地,所述时间差确定子单元,包括:
第一时间差确定子单元,用于对所述两个扫描线位置处对应的所述剪切波的波形进行互相关运算,则互相关最大值所对应的时间延迟量为所述第一时间差;
或者,第二时间差确定子单元,用于分别将每个扫描线位置处对应的所述剪切波的波形中,位移最大值所对应的时刻作为所述剪切波到达所述扫描线位置的时刻,根据剪切波到达所述两个扫描线位置的时间,确定出所述剪切波到达所述两个扫描线位置的第一时间差。
另一方面,本发明还提供了一种剪切波速度测量系统,包括:用于向所述待检测的软组织区域对应的聚焦点位置发射超声脉冲激励波束的第一超声阵列探头;
用于按照预设周期,向所述软组织区域发射多次超声脉冲检测波束,并接收每次反射回来的多个回波信息的第二超声阵列探头;
以及,与所述第一超声阵列探头和第二超声阵列探头连接,并控制所述第一超声阵列探头和第二超声阵列探头的控制节点;
所述控制节点内置有如上所述的剪切波速度测量装置。
经由上述的技术方案可知,与现有技术相比,本发明根据软组织区域内反射回来的回波信息,确定出剪切波在任意两条扫描线间的波速后,根据波速所对应的两条扫描线间的距离,设定波速所对应的权重。而且在设定波速的权重时,考虑到了该波速所对应的扫描线间的距离越小,该波速的偏差就越大,进而将该波速的权重设置的相对较低,从而降低了偏差较大的波速对确定剪切波在软组织区域内的传播速度所造成的影响,进而提高了确定出的剪切波在软组织区域内的传播速度的准确度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。
图1示出了本发明一种剪切波速度测量方法一个实施例的流程示意图;
图2示出了利用超声激励脉冲在软组织区域内产生剪切波,以及向该软组织区域内发射超声脉冲检测波束的原理示意图;
图3示出了本发明一种剪切波速度测量方法另一个实施例的流程示意图;
图4示出了根据一个扫描线位置处的软组织在多个不同时刻的位移量,确定出的剪切波的波形示意图;
图5示出了本发明一种剪切波速度测量装置一个实施例的结构示意图;
图6示出了本发明一种剪切波速度测量装置另一个实施例的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例公开了一种剪切波速度测量方法,以提高对在软组织中传播的剪切波的速度测量的精准度和可靠性。
参见图1,示出了本发明一种剪切波速度测量方法一个实施例的流程示意图,本实施例可以应用于基于声辐射力的剪切波弹性成像技术中,对于剪切波速度的测量;同时,对于通过其他方式在组织内产生剪切波以进行组织病变分析的情况,本发明实施例的剪切波速度测量方法也同样适用。本实施例的方法可以包括:
步骤101:获取待检测的软组织区域针对多次发射的超声脉冲检测波束反射回来的多个回波信息。
其中,待检测的软组织区域内存在剪切波的传播。在本发明中所指的软组织可以是生物体内的肌肉,脂肪、血液以及实质性脏器等。其中,如肝、脾、肾等器官为实质性脏器。
为了能够对该软组织区域进行病变检测,需要使得该软组织区域内存在剪切波传播,以通过测量剪切波在该软组织区域内的传播速度,来判断该软组织区域内是否存在病变。其中,使得软组织区域内存在剪切波传播可以采用现有的任意方式,如,可以通过在组织表面外加低频机械振动,使软组织产生同频振动,机械能量以剪切波的形式向软组织内传播。又如,通过向软组织区域中指定的焦点发射高能量的超声脉冲激励波束,在声辐射力的作用下该焦点处的组织也会发生位移,焦点附近的软组织由于剪切应力的作用也会产生振动,从而在该待检测的软组织区域内产生剪切波的传播。
为了分析出剪切波在软组织区域内的传播速度,需要在该软组织区域产生剪切波传播的情况下,向该软组织区域内发射多次超声脉冲波束,以便接收到该软组织反射回来的回波信息。其中,在软组织内存在剪切波传播的情况下,为了检测剪切波而发射的超声脉冲波束可以称为超声脉冲检测波束。系统接收针对每次发射的超声脉冲检测波束反射回来的回波信息,得到多个回波信息。其中,每次的回波信息是以回波信号的形式接收,因此,回波信息也可以认为是回波信号。例如,向该软组织区域内连续发射了10次超声脉冲检测波束,相应的,会接收到10次回波信号,从而得到不同时刻接收到10个回波信息。
步骤102:根据该多个回波信息,分别确定该剪切波在该软组织区域中预设的任意两条扫描线之间的波速。
为了便于进行剪切波传播速度的分析,在该待检测的软组织区域上可以预先标定多个扫描线,扫描线的位置和数量可以根据需要设定。扫描线也可以称为检测线,本发明实施例中软组织区域内预设的扫描线与现有技术中标定扫描线的方式相同。具体的,每条扫描线均垂直于剪切波的传播方向。
获取到软组织区域内反射回来的回波信息后,可以根据回波信息进行剪切波的波速的分析,具体分析方式可以采用现有技术中任意剪切波的波速计算方式。为了提高计算波速精度,本发明实施例中,需要分别确定出剪切波在该软组织区域中任意两条扫描线之间的波速。例如,当软组织区域内存在3条扫描线:扫描线1、扫描线2和扫描线3,则需要确定剪切波在扫描线1与扫描线2间的波速V1,剪切波在扫描线1与扫描线3之间的波速V2,以及剪切波在扫描线2和扫描线3之间的波速V3,得到对应3个不同子区域内3个波速。
步骤103:依据每个波速对应的两条扫描线的距离,确定每个波速的权重。
其中,波速所对应的两条扫描线的距离越大,该波速的权重越大。当然,该每个波速的权重均为大于零的值。
步骤104:根据每个波速的权重,对所有波速进行加权平均,得到该剪切波在该软组织区域内的传播速度。
在本发明实施例中,在最终确定剪切波在软组织中的传播速度时,并不是直接利用对步骤102中确定出的多个波速求平均值,而是需要根据每个波速所对应的两条扫描线的距离,来确定每个波速的权重,且如果某两条扫描线之间距离相对其他扫描线之间的距离较大,则剪切波在这两条扫描线间的波速所对应权重也相对较大。
例如,仍以软组织区域内设置有扫描线1、扫描线2和扫描线3为例,并假设扫描线1与扫描线3之间的距离最大,扫描线2和扫描线3之间的距离最小,则剪切波在扫描线1和扫描线3之间的波速V1所对应的权重W1最大,而剪切波在扫描线2和扫描线3之间的波速V3所对应的权重W3最小,剪切波在扫描线1和扫描线2之间的波速V2所对应的权重W2则居中。相应的,最终确定出的剪切波在该软组织区域内的传播速度V为:
V = V 1 * W 1 + V 2 * W 2 + V 3 * W 3 W 1 + W 2 + W 3
可选的,对于步骤102中确定出的任意一个波速,计算该波速对应的两条扫描线的距离与第一预设值的乘积,将该乘积作为所述波速的权重,其中,所述第一预设值为大于零的任意值。如,当该第一预设值为一时,该波速对应的扫描线的距离即为该波速对应的权重。
由于软组织区域中传播的剪切波的幅度非常小,一般为微米μm级别的,相应的组织位移幅度也相当弱小。而发射的用于检测剪切波速度的超声脉冲检测波束的波长一般都是100μm以上,因此,对于如此弱小的组织位移幅度,很容易出现受到噪声以及其他干扰等影响,导致测量出的剪切波的速度存在偏差。然而可以理解的是,在外界干扰一定的情况下,在计算剪切波的波速的过程中,两条扫描线之间的距离越近,计算出的剪切波在这两条扫描线之间的波速的偏差也就越大。因此,根据扫描线之间的距离大小,来设定计算出的剪切波在扫描线间的波速所对应的权重,以降低由于扫描线间的波速的偏差较大,对确定剪切波在软组织中的传播速度所造成的影响。
剪切波在两条扫描线之间的波速,仅仅是剪切波在该软组织区域内的一部分子区域的波速,并不能代表该剪切波在整个软组织区域的传播速度。为了描述方便,将剪切波在任意两条扫描线间的波速称为扫描线间的波速,为了确定出该剪切波在该软组织区域内的传播速度,需要利用各个扫描线间的波速的权重,对所有确定出的扫描线间的波速进行加权平均,加权平均的结果即为该剪切波在软组织区域内的传播速度。
在本发明实施例中,根据软组织区域内反射回来的回波信息,确定出剪切波在任意两条扫描线间的波速后,根据波速所对应的两条扫描线间的距离,设定波速所对应的权重。而且在设定波速的权重时,考虑到了该波速所对应的扫描线间的距离越小,确定出的该波速的偏差就越大,进而将该波速的权重设置的相对较低,从而降低了偏差较大的波速对确定剪切波在软组织区域内的传播速度所造成的影响,进而提高了确定出的剪切波在软组织区域内的传播速度的准确度。
可以理解的是,在本发明实施例中,只需获取到待检测的软组织反射回来的回波信息,便可以对软组织内产生的剪切波的传播速度进行分析。对于何时以及如何在软组织区域内产生剪切波,以及如何由软组织区域反射回针对超声脉冲检测波束的回波信息,则可以采用现有的任意方式。当然,为了能够更好的对软组织区域进行检测,本发明的方法可以是基于声辐射力的剪切波弹性成像技术进行软组织区域的病变检测,相应的,可以采用超声激励脉冲来激励软组织区域内产生剪切波。
可选的,在本发明任意实施例中,在获取软组织区域针对多次发射的超声脉冲检测波束反射回来的回波信息之前,还可以包括:控制第一超声阵列探头向待检测的软组织区域对应的聚焦点位置发射超声脉冲激励波束,以在软组织区域内形成剪切波。其中,聚焦点为在该待检测的软组织区域中指定的某一点或者是临近该待检测的软组织区域中的某一点,通过超声脉冲激励波束对该聚焦点进行照射,可以在该软组织区域内产生剪切波。
进一步的,在软组织区域内产生剪切波之后,还可以控制第二超声阵列探头向该软组织区域发射超声脉冲检测波束,并接收反射回来的回波信息,从而使得系统获取接收到的回波信息,并进行后续分析。其中,该第二超声阵列探头和该第一超声阵列探头为同一个超声阵列探头,也可以为不同的超声阵列探头。
为了便于理解利用超声脉冲激励波束在该软组织区域内产生剪切波,并接收回波信息的过程,可以参见图2,其示出了利用超声激励脉冲在软组织区域内产生剪切波,以及通过超声脉冲检测波束对该软组织区域内的剪切波进行测量的原理示意图。
位于图2中间的方框为感兴趣区域(ROI,region of interest)1,所谓感兴趣区域为组织图像中选择的一个图像区域,该图像区域即为需要分析的区域,在该图2中该感兴趣区域可以理解为本发明中待检测的软组织区域的一个平面图像。在该感兴趣区域1右侧用圆圈标出的部分2为超声脉冲激励波束照射的部分,也就是该感兴趣区域对应的聚焦点2。
在该感兴趣区域的上方为超声阵列探头3,在该图2中以发射超声脉冲激励波束和超声脉冲检测波束的阵列探头为同一个为例进行介绍。通过该超声阵列探头3向该聚焦点2发射超声脉冲激励波束31,如图中指向该聚焦点2的带箭头的较细实体直线所示。其中,发射的该超声脉冲激励波束一般为高能量、长持续时间的超声波束。由于超声脉冲激励波束的照射,聚焦点处的组织在声辐射力的作用下会发生位移,该聚焦点附近的组织也会产生振动,从而在该感兴趣区域1的软组织区域内产生剪切波,如感兴趣区域1内的波形所示。其中,剪切波是传播方向与介质质点振动方向垂直的波,在该图2中,感兴趣区域内的软组织的振动方向为上下振动,该剪切波则为从聚焦点向左沿水平方向传播。
为了测量该剪切波在该感兴趣区域的软组织中的传播速度,该超声阵列探头3按照预设的周期向该感兴趣区域发射超声脉冲检测波束32,如图2中带箭头的较粗的实体直线所示。由图2可以看出超声脉冲检测波束的聚焦焦点在该感兴趣区域1的正下方,使得整个感兴趣区域能够完全被该超声脉冲检测波束覆盖,使得不同位置处的超声脉冲的声强均匀一致,且该感兴趣区域1内有足够强的超声脉冲检测波束经过。
其中,该超声阵列探头本次发射超声脉冲检测波束距离下一次发射超声脉冲检测波束的时间为脉冲重复时间。每隔一个脉冲重复时间,该超声阵列探头均会发送一次超声脉冲检测波束,发射该超声脉冲检测波束的总次数可以根据实际需要设定。每次向该感兴趣区域发射了超声脉冲检测波束后,该超声阵列探头都会接收由该感兴趣区域反射回来的回波信息,进而系统可以获取该回波信息,来分析感兴趣区域内剪切波的传播速度。
可选的,在实际应用中,当获取到软组织区域反射回来的回波信息后,为了能够最终确定出任意两条扫描线之间的波速,可以根据获取到的多个回波信息并结合波束合成技术,确定该软组织区域中,预设的每个扫描线位置针对多次发射的超声脉冲检测波束反射的多个回波信号。为了便于与整个待检测的软组织区域内反射回来的回波信息相区分,将根据这回波信息,并利用波束合成技术,合成的每个扫描线位置对应的回波信号称为第一回波信号。也就是说,对于每次发射超声检测脉冲波束,获取到该软组织区域针对本次超声脉冲检测波束反射回的回波信息后,均需要利用波束合成技术,确定出每个扫描线位置处反射回的第一回波信号,其中,不同扫描线位置所对应的第一回波信号不同。当然,根据回波信息并结合波束合成技术,合成每个扫描线位置处反射回的第一回波信号的过程与现有的合成扫描线位置处或者某个位置处的回波信号的过程相似,在此不再赘述。
在此之后,可以根据每个扫描线位置所对应的多个第一回波信号,分别确定剪切波在该软组织区域中预设的任意两条扫描线之间的波速。
其中,软组织区域中预设的扫描线并不是实际存在于生物体的软组织内,而是为了剪切波波速分析的需要,人为标定的线条,以根据合成的该扫描线位置处反射的第一回波信号来分析剪切波速度。为了便于理解,可以参看图2所示,在该感兴趣区域上标示有4个扫描线位置,如图中标示有A、B、C、D的四点,在该4点位置处上分别设置有4条扫描线,该扫描线垂直于该剪切波的传播方向,如图在感兴趣区域内分别位于A、B、C、D位置处相互平行虚线所示。。因此,对于每次接收到的回波信息,均需要分别合成位于位置A、B、C、D处扫描线位置处的软组织反射回的第一回波信号。为了便于描述,将位于位置A、B、C、D处扫描线位置处分别称为扫描线A、扫描线B、扫描线C和扫描线D,对于同一次接收到的回波信息,合成的该扫描线A位置处对应的第一回波信号,该扫描线B位置处对应的第一回波信号,扫描线C位置处对应的第一回波信号,以及扫描线D位置处对应的第一回波信号可能均不相同。
为了便于理解合成的各个扫描线位置处反射回去的第一回波信号,在图2中,在每个扫描线位置处还对应了画出了代表通过合成技术,确定出的各个扫描线位置处返回的回波信号,如图中虚线所示将感兴趣区域内的4条扫描线延伸至超声阵列探头处,以表示确定出的各个扫描线位置处反射回的回波信号。假设超声脉冲检测波束的脉冲重复时间为100微秒,连续对该软组织区域发射了20次超声脉冲检测波束,那么相应的会接收到20次回波信息,每次接收回波信息的时刻也相差100微秒。依据结合波束合成技术,分别对于任意一次接收到回波信息,确定出各个扫描线所对应的第一回波信号,则每个扫描线也均对应了20个不同时刻反射的第一回波信号。
可选的,为了计算剪切波在任意两条扫描线之间的波速,可以采用先根据获取到的多个回波信息,确定剪切波到达该软组织区域内预设的任意两条扫描线的第一时间差。当然,也可以先根据获取到的该多个回波信息,并结合波束合成技术,确定该软组织区域中,预设的每个扫描线位置针对多次发射的超声脉冲检测波束反射回来的多个第一回波信号。对于任意两条扫描线位置,如假设为第一扫描线位置和第二扫描线位置,则可以根据第一扫描线位置处针对每次发射的超声波反射回来的第一回波信息,以及第二扫描线位置处针对每次发射的超声波反射回来的第一回波信息,确定出剪切波到达该第一扫描线位置和第二扫描线位置的时间差,即第一时间差。
当剪切波到达任意两条扫描线的第一时间差确定后,对于任意两条扫描线,依据这两条扫描线的第一距离以及该第一时间差,便可以计算出剪切波在该两条扫描线之间的波速。
例如,仍以图2所示为例,在软组织区域中设定有4条扫描线:扫描线A、扫描线B、扫描线C和扫描线D,假设需要确定剪切波在扫描线A和扫描线B之间的波速时,则可以首先确定出剪切波到达扫描线A和扫描线B的第一时间差,然后确定出扫描线A与扫描线B之间的距离,将该距离除以该第一时间差即可得到该剪切波在这两条扫描线之间的波速。
其中,确定剪切波到达两条扫描线的第一时间差可以是先分别确定剪切波到达每条扫描线的时间,然后再计算该第一时间差。当然,剪切波到达扫描线的时间可以采用现有任意方式,在此不加以限制。
参见图3,示出了一种计算剪切波速度测量方法另一个实施例的流程示意图,本实施例的方法可以应用于基于声辐射力的剪切波弹性成像技术中,对于剪切波速度的测量;同时,对于通过其他方式在组织内产生剪切波以进行组织病变分析的情况,本发明实施例的剪切波速度测量方法也同样适用。本实施例中提供了计算剪切波到达两条扫描线的时间差的一种优选方式。本实施例的方法可以包括:
步骤301:获取待检测的软组织区域针对多次发射的超声脉冲检测波束反射回来的多个回波信息。
其中,该待检测的软组织区域内存在剪切波传播。
该步骤可以参见以上实施例中相应步骤的描述,在此不再赘述。
步骤302:根据多个回波信息,并结合波束合成技术,确定软组织区域中,预设的每个扫描线位置针对所述多次发射的超声脉冲检测波束反射回来的多个第一回波信号。
该步骤可以参照前面相关描述,在此不再赘述。
步骤303:对于任意一条扫描线,根据预先存储的扫描线位置的静态组织回波信号以及该多个第一回波信号,分别确定该扫描线位置处的软组织在多个不同时刻的位移量。
其中,该扫描线位置处的软组织在多个不同时刻位移量分别为该软组织在反射回这多个第一回波信号的多个时刻所对应的位移量。也就是说,对于任意一个扫描线位置的软组织而言,每次发射超声脉冲检测波束,该扫描线位置处偏离静态组织的位移量均不相同,需要分别确定每次发射的超声波脉冲检测波束,引起该扫描线位置处的软组织移动的位移变化量。
其中,静态组织回波信号是指在向该软组织区域发射超声脉冲激励波束之前,向该软组织区域发射超声脉冲检测波束,并根据接收到的回波信息,确定出的各个扫描线位置处所对应的回波信号。简而言之,该静态组织回波信号就是在软组织区域中不存在剪切波传播的情况下,反映各个扫描线位置处软组织状态的回波信号。根据扫描线位置处的软组织在剪切波传播前的静态组织回波信号,以及存在剪切波传播的状态下,该扫描线位置处的软组织多次反射的第一回波信号可以分别确定出每次发射超声脉冲检测波束后该扫描线位置处的软组织偏离静态位置的位移量,得到多个时刻的多个位移量。
如,其中一种计算一个扫描线位置处软组织的位移量方式可以是在该扫描线位置处的扫描线上选取多个点,例如在该图2中扫描线A的位置上选取若干个点,然后根据该扫描线上选取的各个点反射回来的回波信号,来确定出该扫描线位置处软组织偏离静态组织的位移量。具体的公式如下:
[ v j , r j ] = max ( Σ t = 1 m S ‾ ( t ) × S j ( t + τ ) ) (公式一)
其中,t=1,2….m,m表示在该扫描线位置处选取的总点数,也就是软组织区域上标定的该扫描线处选取的总点数。该是该扫描线上第t个点处的静态组织回波的射频信号,Sj(t)代表在软组织区域内存在剪切波的情况下,第j次发射超声脉冲检测波束后,该扫描线上第t个点所对应的第一回波信号,j=1,2·····N,N为发射超声脉冲检测波束的总次数。τ是第t次
Figure BDA0000389913640000143
与Sj互相关运算的偏移量,τ的取值为1,2·····m。vj
Figure BDA0000389913640000144
与Sj(t)的互相关的最大值。当
Figure BDA0000389913640000145
与Sj的自相关值达到最大值Vj时,所得到的rj值即为组织偏离静态位置的位移量。这样,对于根据每次发射的超声脉冲检测波束,得到的该扫描线位置处的第一回波信号,均分别与该扫描线位置处的静态组织回波信号分别进行m次互相关。也就是说,需要分别针对该扫描线上选取的各个点对应的第一回波信号Sj(t)与扫描线上相应点的静态组织回波信号
Figure BDA0000389913640000146
进行互相关,从而对该扫描线上m个位置点处分别进行第一回波信号与静态组织回波信号的互相关,并找出每次互相关的延迟量τ。比较各个点Sj(t)与S_(t)的互相关值,并确定Sj(t)与
Figure BDA0000389913640000147
的互相最大值,将该Sj(t)与
Figure BDA0000389913640000148
的互相关最大值时所对应的延迟量。τ赋值给rj,则此时得到的rj即为在第j次发射超声脉冲检测波束时,该扫描线位置处的软组织偏离静态位置的位移量。
又如,另一种计算一个扫描线位置处的软组织的位移量的过程如下:
I=I1×Qj-Q1×Ij               (公式二)
Q=I1×Ij-Q1×Qj               (公式三)
S j = c × f 4 × π × arctan ( I Q ) (公式四)
其中,I1和Q1分别是由该扫描线位置处的静态组织回波信号计算出的同相信号和正交信号,Qj和Ij是对该扫描线位置的第j次发射超声检测脉冲后,得到对应的第一回波信号解调得到的同相信号和正交信号。f是发射超声检测脉冲的中心频率,c是超声检测脉冲在软组织中的传播速度,Sj即为扫描线位置处对应第j次发射超声检测脉冲波束时,该软组织偏离静态位置的位移量。
步骤304:根据每个扫描线位置处的软组织在多个不同时刻的位移量,确定每个扫描线位置处对应的该剪切波的波形。
在确定出多个时刻发射超声脉冲检测波束之后,该扫描线位置处的软组织在这多个时刻所对应的位移量,便可以根据该扫描线位置的软组织在多个时刻的位移量,构造成该扫描线位置处所对应的剪切波的波形。由确定该扫描线位置处对应的剪切波的波形的过程可知,该扫描线位置处对应的剪切波的波形也可以理解为是在剪切波传播中,该扫描线位置处的组织的位移波形。
如,以图2为例,仍假设超声脉冲检测波束的脉冲重复时间为100微秒,连续对该软组织区域发射了20次超声脉冲检测波束,那么相应的会接收到20次回波信息,每次接收回波信息的时刻也近似相差100微秒。依据结合波束合成技术,分别对于任意一次接收到回波信息,确定出各个扫描线所对应的第一回波信号,则每个扫描线也均对应了20个不同时刻反射的第一回波信号。假设在第0微秒第一次发射超声脉冲检测波束,然后接收该次反射的回波信号;之后会在第100微秒时,第二次发射超声脉冲检测波束,并接收该次反射回来的回波信息,依次类推。
以确定扫描线A位置处对应的剪切波的波形为例进行介绍,根据该扫描线A位置处的静态组织回波信号,以及根据第一次反射回来的回波信息确定出的该扫描线A位置处对应的第一回波信号,可以确定出在该0微秒时,该扫描线A位置处的软组织的位移量。然后依次确定后续第100微秒时,该扫描线A位置处的软组织的为位移量,依次类推,得到该扫描位置处在20个不同时刻的位移量。根据这20个不同时刻的20个位移量,可以得到该扫描线位置处的软组织对应的剪切波的波形,如图4所示。在图4中横坐标代表发射超声脉冲检测波束的各个时刻,而纵坐标对应各个时刻的位移量。其中,图4中任意两点之间的横坐标之差为一个脉冲重复时间即100微秒。如,表示第一次确定的位移量的时刻t1,表示第二次确定出的位移量的时刻t2,t2与t1的时间差即为一个脉冲重复时间。
当然,图4仅仅是一种确定该扫描线位置处对应的剪切波的一种示意图,在实际应用中,还可以对该扫描线位置处在多个不同时刻的位移量进行插值,然后确定出该扫描线位置处对应的剪切波波形,具体过程和现有的对多个时刻点的离散数据进行插值,来确定波形的方式相似,在此不再赘述。
步骤305:对于任意两个扫描线位置,根据这两个扫描线位置处对应的剪切波的波形,确定剪切波到达该两个扫描线位置的第一时间差。
确定出每个扫描线位置处对应的剪切波的波形后,确定剪切波到达任意两个扫描线位置的第一时间差也可以有多种方式。
如,其中一种方式为:分别将该每个扫描线位置处对应的剪切波的波形中,位移最大值所对应的时刻作为剪切波到达该扫描线位置的时刻。例如,对于图4所示的波形中,位移最大值所对应的时刻为时刻t6,则该时刻t6即为剪切波到达该扫描线位置A的时刻。这样,对于任意两个扫描线位置,则可以将剪切波到达这两个扫描线位置的时刻相减,便可以得到剪切波到达这两个扫描线位置的第一时间差。
又如,其中另一种计算该第一时间差的方式可以为:对于任意两个扫描线位置,对这两个扫描线位置处对应的剪切波的波形进行互相关运算,则互相关最大值所对应的时间延迟量为该第一时间差。其中,互相关运算如下所示:
[ v , r ij ] = max ( Σ t = 1 n S i ( t ) × S j ( t + τ ) ) (公式五)
其中,Si(t)为第i条扫描线对应的剪切波的波形,Sj(t)为第j条扫描线对应的剪切波的波形。
t为剪切波的波形中对应的时刻中选取的时刻点,n为选取的时刻点的总个数,其中t取值为1,2····n,n为大于一的自然数,n的具体值可以根据需要设定,t=1表示选取出的第一个时刻点,t=2表示选取出的第二个时刻点,依次类推。τ为互相关偏移量,τ=1,2….n。相应的,τ=1,表示偏移第一个时刻点的时长。可选的,t可以为第t次发射超声脉冲检测波束的时刻,相应的,n为发射脉检测波束的总次数,如,t=n,则第n次发射超声脉冲检测波束的时刻。
由于可以有n个取值,因此,Si(t)与Sj(t)之间需要作n次互相关。v表示这n次互相关运算的最大值。rij该第一时间差。找出这n次互相关运算的最大值v,确定得到该互相关运算最大值v时,所对应的偏移量τ,将该偏移量rij赋值给rij,得到剪切波到达这两条扫描线的第一时间差。
步骤306:对于任意两条扫描线,依据两条扫描线的第一距离以及该剪切波到达该两条扫描线的第一时间差,计算剪切波在该两条扫描线之间的波速。
步骤307:依据每个波速对应的两条扫描线的距离,确定每个波速的权重。
步骤308:根据每个波速的权重,对所有波速进行加权平均,得到该剪切波在该软组织区域内的传播速度。
在本实施例中,在确定扫描线位置处对应的剪切波的波形时,任意相邻的两个位移量所对应的时间差均为一个脉冲重复时间。而剪切波到达任意两条扫描线位置处的时间差,实际上也是根据多个时刻的位移量这样一个离散序列来确定的,因此,剪切波到达任意两条扫描线位置处的时间差也一定是脉冲重复时间的倍数。换言之,剪切波到达任意两条扫描线的时间差是以脉冲重复时间为最小计量单元,该时间差一定是若干个脉冲重复时间。这样,当假设剪切波到达任意两个扫描线位置的时间差的误差量都在同一个范围内时,则两条扫描线之间的距离越远,则确定出的剪切波在这两条扫描线之间的波速也就更加趋近真实值。
例如,假设剪切波到达任意两条扫描线的时间差的误差都是一个脉冲重复时间,如果两个扫描线的距离比较远,如图2中的扫描线A和扫描线C,剪切波需要10个脉冲重复时间从能够从扫描线A传播到扫描线C,那么计算出的剪切波在这两条扫描线之间传播的波速的误差就只有10%。如果两个扫描线距离比较近的话,如图中C线和B线,剪切波只需要5个脉冲重复时间就能从扫描线B到达扫描线C,那么在剪切波到达时间差的误差为一个脉冲重复时间的前提下,那么计算出剪切波在扫描线B和扫描线C间传播的波速的误差就达到20%。
可见,在一定范围内,在同样误差量的前提下,如果两条扫描线之间的距离越大,计算出的剪切波在这两条扫描线之间的波速会更加精确。因此,在对所有扫描线间波速进行加权平均时,将距离较大的两条扫描线间的波速所对应的权重增大,有利于提高软组织区域内剪切波的传播速度的测量精度。
本发明还提供了一种剪切波速度测量装置,用于实现本发明上述实施例提供的剪切波速度测量方法。以下举例说明。
参见图5,示出了本发明一种剪切波速度测量装置一个实施例的结构示意图,本实施例装置可以包括:信息获取单元501、分段波速确定单元502、权重设置单元503和剪切波速度确定单元504。
其中,信息获取单元501,用于获取待检测的软组织区域针对多次发射的超声脉冲检测波束反射回来的多个回波信息,其中,待检测的软组织区域内存在剪切波的传播。
为了能够对该软组织区域进行病变检测,需要使得该软组织区域内存在剪切波传播,以通过测量剪切波在该软组织区域内的传播速度,来判断该软组织区域内是否存在病变。其中,使得软组织区域内存在剪切波传播可以采用现有的任意方式,如,可以通过在组织表面外加低频机械振动,使软组织产生同频振动,机械能量以剪切波的形式向软组织内传播。又如,通过向软组织区域中指定的焦点发射高能量的超声脉冲激励波束,在声辐射力的作用下该焦点处的组织也会发生位移,焦点附件的软组织由于剪切应力的作用,也会产生振动,从而在该待检测的软组织区域内产生剪切波的传播。
分段波速确定单元502,用于根据所述多个回波信息,分别确定所述剪切波在所述软组织区域中预设的任意两条扫描线之间的波速。
权重设置单元503,用于依据每个所述波速对应的两条扫描线的距离,确定每个所述波速的权重,其中,所述波速所对应的两条扫描线的距离越大,所述波速的权重越大。
剪切波速度确定单元504,用于根据每个所述波速的权重,对所有所述波速进行加权平均,得到所述剪切波在所述软组织区域内的传播速度。
其中,该权重设置单元503设置权重的方式有多种,对应其中一种方式,该权重设置单元,可以包括:
第一权重设置单元,用于对于任意一个波速,计算所述波速对应的两条扫描线的距离与第一预设值的乘积,将所述乘积作为所述波速的权重,其中,所述第一预设值为大于零的任意值。
可选地,该剪切波速度测量装置还可以包括:第一控制单元,用于控制第一超声阵列探头向所述待检测的软组织区域对应的聚焦点位置发射超声脉冲激励波束,以在所述软组织区域内形成剪切波;
第二控制单元,用于控制第二超声阵列探头按照预设周期,向所述软组织区域发射多次超声脉冲检测波束,并接收每次反射回来的所述回波信息。
参见图6示出了本发明一种剪切波速度测量装置另一个实施例的结构示意图,本实施例的装置与上一装置的实施例不同之处在于:
在本实施例中,该分段波速确定单元502,可以包括:
时间差确定单元5021,用于根据所述多个回波信息,确定所述剪切波到达所述软组织区域内预设的任意两条扫描线的第一时间差;
分段波速确定子单元5022,用于对于任意两条扫描线,依据所述两条扫描线的第一距离以及所述剪切波到达所述两条扫描线的第一时间差,计算所述剪切波在所述两条扫描线之间的波速。
可选地,该时间差确定单元5021,可以包括:
信号合成单元,用于根据所述多个回波信息,并结合波束合成技术,确定所述软组织区域中,预设的每个扫描线位置针对所述多次发射的超声脉冲检测波束反射回来的多个第一回波信号;
位移量确定单元,用于对于任意一个扫描线位置,根据预先存储的所述扫描线位置的静态组织回波信号以及所述多个第一回波信号,分别确定所述扫描线位置处的软组织在多个不同时刻的位移量,其中,所述在多个不同时刻位移量分别为所述软组织在反射回所述多个第一回波信号的多个时刻所对应的位移量;
波形确定单元,用于根据每个所述扫描线位置处的软组织在多个不同时刻的位移量,确定每个扫描线位置处对应的所述剪切波的波形;
时间差确定子单元,用于对于任意两个扫描线位置,根据所述两个扫描线位置处对应的所述剪切波的波形,确定所述剪切波到达所述两个扫描线位置的第一时间差。
进一步地,该时间差确定子单元,可以包括:
第一时间差确定子单元,用于对所述两个扫描线位置处对应的所述剪切波的波形进行互相关运算,则互相关最大值所对应的时间延迟量为所述第一时间差;
或者,第二时间差确定子单元,用于分别将每个扫描线位置处对应的所述剪切波的波形中,位移最大值所对应的时刻作为所述剪切波到达所述扫描线位置的时刻,根据剪切波到达所述两个扫描线位置的时间,确定出所述剪切波到达所述两个扫描线位置的第一时间差。
另外,本发明还提供了一种剪切波速度测量系统,该系统包括:
用于向所述待检测的软组织区域对应的聚焦点位置发射超声脉冲激励波束的第一超声阵列探头;
用于按照预设周期,向所述软组织区域发射多次超声脉冲检测波束,并接收每次反射回来的多个回波信息的第二超声阵列探头;
以及,与所述第一超声阵列探头和第二超声阵列探头连接,并控制所述第一超声阵列探头和第二超声阵列探头的控制节点。
其中,该控制节点内置有如上任意一个实施例所描述的剪切波速度测量装置。
需要说明的是,该第一超声波阵列探头和第二超声波阵列探头可以是同一个,也可以分别是两个阵列探头。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (13)

1.一种剪切波速度测量方法,其特征在于,包括:
获取待检测的软组织区域针对多次发射的超声脉冲检测波束反射回来的多个回波信息,其中,所述待检测的软组织区域内存在剪切波的传播;
根据所述多个回波信息,分别确定所述剪切波在所述软组织区域中预设的任意两条扫描线之间的波速;
依据每个所述波速对应的两条扫描线的距离,确定每个所述波速的权重,其中,所述波速所对应的两条扫描线的距离越大,所述波速的权重越大;
根据每个所述波速的权重,对所有所述波速进行加权平均,得到所述剪切波在所述软组织区域内的传播速度。
2.根据权利要求1所述的方法,其特征在于,所述依据每个所述波速对应的两条扫描线的距离,确定每个所述波速的权重,包括:
对于任意一个所述波速,计算所述波速对应的两条扫描线的距离与第一预设值的乘积,将所述乘积作为所述波速的权重,其中,所述第一预设值为大于零的任意值。
3.根据权利要求1所述的方法,其特征在于,所述根据所述多个回波信息,分别确定所述剪切波在所述软组织区域中预设的任意两条扫描线之间的波速,包括:
根据所述多个回波信息,并结合波束合成技术,确定所述软组织区域中,预设的每个扫描线位置针对所述多次发射的超声脉冲检测波束反射的多个第一回波信号;
根据所述每个扫描线位置所对应的多个第一回波信号,分别确定所述剪切波在所述软组织区域中预设的任意两条扫描线之间的波速。
4.根据权利要求1至3任一项所述的方法,其特征在于,所述根据所述多个回波信息,分别确定所述剪切波在所述软组织区域中预设的任意两条扫描线之间的波速,包括:
根据所述多个回波信息,确定所述剪切波到达所述软组织区域内预设的任意两条扫描线的第一时间差;
对于任意两条扫描线,依据所述两条扫描线的第一距离以及所述剪切波到达所述两条扫描线的第一时间差,计算所述剪切波在所述两条扫描线之间的波速。
5.根据权利要求4所述的方法,其特征在于,所述根据所述多个回波信息,确定所述剪切波到达所述软组织区域内预设的任意两条扫描线的第一时间差,包括:
根据所述多个回波信息,并结合波束合成技术,确定所述软组织区域中,预设的每个扫描线位置处针对所述多次发射的超声脉冲检测波束反射回来的多个第一回波信号;
对于任意一个扫描线位置,根据预先存储的所述扫描线位置的静态组织回波信号以及所述多个第一回波信号,分别确定所述扫描线位置处的软组织在多个不同时刻的位移量,其中,所述在多个不同时刻位移量分别为所述软组织在反射回所述多个第一回波信号的多个时刻所对应的位移量;
根据每个所述扫描线位置处的软组织在多个不同时刻的位移量,确定每个扫描线位置处对应的所述剪切波的波形;
对于任意两个扫描线位置,根据所述两个扫描线位置处对应的所述剪切波的波形,确定所述剪切波到达所述两个扫描线位置的第一时间差。
6.根据权利要求5所述的方法,其特征在于,所述根据所述两个扫描线位置处对应的所述剪切波的波形,确定所述剪切波到达所述两个扫描线位置的第一时间差,包括:
对所述两个扫描线位置处对应的所述剪切波的波形进行互相关运算,则互相关最大值所对应的时间延迟量为所述第一时间差;
或者,分别将每个扫描线位置处对应的所述剪切波的波形中,位移最大值所对应的时刻作为所述剪切波到达所述扫描线位置的时刻,根据剪切波到达所述两个扫描线位置的时间,确定出所述剪切波到达所述两个扫描线位置的第一时间差。
7.根据权利要求1所述的方法,其特征在于,所述获取待检测的软组织区域针对发射的超声脉冲检测波束反射回来的回波信息之前,还包括:
控制第一超声阵列探头向所述软组织区域对应的聚焦点位置发射超声脉冲激励波束,以在所述软组织区域内形成剪切波;
控制第二超声阵列探头按照预设周期,向所述软组织区域发射多次超声脉冲检测波束,并接收每次反射回来的所述回波信息。
8.一种剪切波速度测量装置,其特征在于,包括:
信息获取单元,用于获取待检测的软组织区域针对多次发射的超声脉冲检测波束反射回来的多个回波信息,其中,所述待检测的软组织区域内存在剪切波的传播;
分段波速确定单元,用于根据所述多个回波信息,分别确定所述剪切波在所述软组织区域中预设的任意两条扫描线之间的波速;
权重设置单元,用于依据每个所述波速对应的两条扫描线的距离,确定每个所述波速的权重,其中,所述波速所对应的两条扫描线的距离越大,所述波速的权重越大;
剪切波速度确定单元,用于根据每个所述波速的权重,对所有所述波速进行加权平均,得到所述剪切波在所述软组织区域内的传播速度。
9.根据权利要求8所述的装置,其特征在于,所述权重设置单元,包括:
第一权重设置单元,用于对于任意一个波速,计算所述波速对应的两条扫描线的距离与第一预设值的乘积,将所述乘积作为所述波速的权重,其中,所述第一预设值为大于零的任意值。
10.根据权利要求8或9所述的装置,其特征在于,所述分段波速确定单元,包括:
时间差确定单元,用于根据所述多个回波信息,确定所述剪切波到达所述软组织区域内预设的任意两条扫描线的第一时间差;
分段波速确定子单元,用于对于任意两条扫描线,依据所述两条扫描线的第一距离以及所述剪切波到达所述两条扫描线的第一时间差,计算所述剪切波在所述两条扫描线之间的波速。
11.根据权利要求10所述的装置,其特征在于,所述时间差确定单元,包括:
信号合成单元,用于根据所述多个回波信息,并结合波束合成技术,确定所述软组织区域中,预设的每个扫描线位置针对所述多次发射的超声脉冲检测波束反射回来的多个第一回波信号;
位移量确定单元,用于对于任意一个扫描线位置,根据预先存储的所述扫描线位置的静态组织回波信号以及所述多个第一回波信号,分别确定所述扫描线位置处的软组织在多个不同时刻的位移量,其中,所述在多个不同时刻位移量分别为所述软组织在反射回所述多个第一回波信号的多个时刻所对应的位移量;
波形确定单元,用于根据每个所述扫描线位置处的软组织在多个不同时刻的位移量,确定每个扫描线位置处对应的所述剪切波的波形;
时间差确定子单元,用于对于任意两个扫描线位置,根据所述两个扫描线位置处对应的所述剪切波的波形,确定所述剪切波到达所述两个扫描线位置的第一时间差。
12.根据权利要求11所述的装置,其特征在于,所述时间差确定子单元,包括:
第一时间差确定子单元,用于对所述两个扫描线位置处对应的所述剪切波的波形进行互相关运算,则互相关最大值所对应的时间延迟量为所述第一时间差;
或者,第二时间差确定子单元,用于分别将每个扫描线位置处对应的所述剪切波的波形中,位移最大值所对应的时刻作为所述剪切波到达所述扫描线位置的时刻,根据剪切波到达所述两个扫描线位置的时间,确定出所述剪切波到达所述两个扫描线位置的第一时间差。
13.一种剪切波速度测量系统,其特征在于,包括:用于向所述待检测的软组织区域对应的聚焦点位置发射超声脉冲激励波束的第一超声阵列探头;
用于按照预设周期,向所述软组织区域发射多次超声脉冲检测波束,并接收每次反射回来的多个回波信息的第二超声阵列探头;
以及,与所述第一超声阵列探头和第二超声阵列探头连接,并控制所述第一超声阵列探头和第二超声阵列探头的控制节点;
所述控制节点内置有如权利要求8至12任一项所述的剪切波速度测量装置。
CN201310454725.3A 2013-09-29 2013-09-29 一种剪切波速度测量方法、装置和系统 Active CN103462643B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310454725.3A CN103462643B (zh) 2013-09-29 2013-09-29 一种剪切波速度测量方法、装置和系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310454725.3A CN103462643B (zh) 2013-09-29 2013-09-29 一种剪切波速度测量方法、装置和系统

Publications (2)

Publication Number Publication Date
CN103462643A true CN103462643A (zh) 2013-12-25
CN103462643B CN103462643B (zh) 2015-02-11

Family

ID=49787860

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310454725.3A Active CN103462643B (zh) 2013-09-29 2013-09-29 一种剪切波速度测量方法、装置和系统

Country Status (1)

Country Link
CN (1) CN103462643B (zh)

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104605890A (zh) * 2014-12-18 2015-05-13 深圳开立生物医疗科技股份有限公司 一种剪切波峰值波形校正的方法、装置、系统及其应用
WO2016029402A1 (zh) * 2014-08-28 2016-03-03 深圳迈瑞生物医疗电子股份有限公司 剪切波成像方法及系统
CN105877783A (zh) * 2015-02-15 2016-08-24 深圳开立生物医疗科技股份有限公司 二维剪切波弹性成像方法和装置
CN107510474A (zh) * 2017-09-21 2017-12-26 深圳开立生物医疗科技股份有限公司 剪切波弹性成像方法及系统
CN107833218A (zh) * 2017-11-22 2018-03-23 深圳中科乐普医疗技术有限公司 一种生物组织剪切波弹性图像质量评估显示方法
CN107970043A (zh) * 2017-12-28 2018-05-01 深圳开立生物医疗科技股份有限公司 一种剪切波的检测方法及装置
CN108065964A (zh) * 2018-01-16 2018-05-25 中国科学院苏州生物医学工程技术研究所 一种超声成像方法、装置、设备及超声成像探头
CN108338808A (zh) * 2017-01-23 2018-07-31 美国西门子医疗解决公司 使用相干的剪切速度成像
CN108852416A (zh) * 2018-05-17 2018-11-23 深圳开立生物医疗科技股份有限公司 一种剪切波传播速度的确定方法及装置
CN109100801A (zh) * 2018-07-04 2018-12-28 中国科学院地质与地球物理研究所 一种高压介质地震波传播模拟方法
CN110327077A (zh) * 2019-07-09 2019-10-15 深圳开立生物医疗科技股份有限公司 一种血流显示方法、装置及超声设备和存储介质
CN110432926A (zh) * 2014-09-03 2019-11-12 深圳迈瑞生物医疗电子股份有限公司 弹性测量检测方法及系统
WO2020019254A1 (zh) * 2018-07-26 2020-01-30 深圳迈瑞生物医疗电子股份有限公司 一种剪切波成像方法及系统
CN111671462A (zh) * 2020-06-18 2020-09-18 京东方科技集团股份有限公司 超声成像数据采集的扫描间隔计算方法及其设备、装置
CN111700639A (zh) * 2019-03-18 2020-09-25 三星麦迪森株式会社 超声诊断设备及其控制方法
CN111728642A (zh) * 2017-06-06 2020-10-02 深圳迈瑞生物医疗电子股份有限公司 一种在超声扫描中成像的方法、装置及系统
US11540809B2 (en) * 2019-07-12 2023-01-03 Konica Minolta, Inc. Ultrasonic diagnostic apparatus and method for propagation speed analysis of shear wave and elastic modulus measurement of a tissue

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100016718A1 (en) * 2008-07-16 2010-01-21 Siemens Medical Solutions Usa, Inc. Shear Wave Imaging
CN102551801A (zh) * 2010-10-06 2012-07-11 美国西门子医疗解决公司 对医学超声成像中剪切波信息的求解
WO2012105152A1 (ja) * 2011-02-04 2012-08-09 株式会社日立メディコ 超音波診断装置,および方法
CN102667522A (zh) * 2009-11-25 2012-09-12 皇家飞利浦电子股份有限公司 采用聚焦扫描线波束形成的超声剪切波成像
CN102892358A (zh) * 2010-04-05 2013-01-23 日立阿洛卡医疗株式会社 用于超声成像的方法和装置
CN103269639A (zh) * 2010-12-22 2013-08-28 皇家飞利浦电子股份有限公司 利用质心估计剪切波速度

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100016718A1 (en) * 2008-07-16 2010-01-21 Siemens Medical Solutions Usa, Inc. Shear Wave Imaging
CN102667522A (zh) * 2009-11-25 2012-09-12 皇家飞利浦电子股份有限公司 采用聚焦扫描线波束形成的超声剪切波成像
CN102892358A (zh) * 2010-04-05 2013-01-23 日立阿洛卡医疗株式会社 用于超声成像的方法和装置
CN102551801A (zh) * 2010-10-06 2012-07-11 美国西门子医疗解决公司 对医学超声成像中剪切波信息的求解
CN103269639A (zh) * 2010-12-22 2013-08-28 皇家飞利浦电子股份有限公司 利用质心估计剪切波速度
WO2012105152A1 (ja) * 2011-02-04 2012-08-09 株式会社日立メディコ 超音波診断装置,および方法

Cited By (27)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016029402A1 (zh) * 2014-08-28 2016-03-03 深圳迈瑞生物医疗电子股份有限公司 剪切波成像方法及系统
CN110432926A (zh) * 2014-09-03 2019-11-12 深圳迈瑞生物医疗电子股份有限公司 弹性测量检测方法及系统
CN104605890A (zh) * 2014-12-18 2015-05-13 深圳开立生物医疗科技股份有限公司 一种剪切波峰值波形校正的方法、装置、系统及其应用
CN105877783A (zh) * 2015-02-15 2016-08-24 深圳开立生物医疗科技股份有限公司 二维剪切波弹性成像方法和装置
CN105877783B (zh) * 2015-02-15 2019-04-16 深圳开立生物医疗科技股份有限公司 二维剪切波弹性成像方法和装置
CN108338808A (zh) * 2017-01-23 2018-07-31 美国西门子医疗解决公司 使用相干的剪切速度成像
CN108338808B (zh) * 2017-01-23 2024-05-24 美国西门子医疗解决公司 使用相干的剪切速度成像
CN111728642A (zh) * 2017-06-06 2020-10-02 深圳迈瑞生物医疗电子股份有限公司 一种在超声扫描中成像的方法、装置及系统
US11717270B2 (en) 2017-06-06 2023-08-08 Shenzhen Mindray Bio-Medical Electronics Co., Ltd. Method, apparatus and system for imaging in ultrasonic scanning
CN107510474B (zh) * 2017-09-21 2020-07-10 深圳开立生物医疗科技股份有限公司 剪切波弹性成像方法及系统
CN107510474A (zh) * 2017-09-21 2017-12-26 深圳开立生物医疗科技股份有限公司 剪切波弹性成像方法及系统
CN107833218B (zh) * 2017-11-22 2021-01-26 深圳中科乐普医疗技术有限公司 一种生物组织剪切波弹性图像质量评估显示方法
CN107833218A (zh) * 2017-11-22 2018-03-23 深圳中科乐普医疗技术有限公司 一种生物组织剪切波弹性图像质量评估显示方法
CN107970043A (zh) * 2017-12-28 2018-05-01 深圳开立生物医疗科技股份有限公司 一种剪切波的检测方法及装置
CN108065964B (zh) * 2018-01-16 2021-04-20 中国科学院苏州生物医学工程技术研究所 一种超声成像方法、装置、设备及超声成像探头
CN108065964A (zh) * 2018-01-16 2018-05-25 中国科学院苏州生物医学工程技术研究所 一种超声成像方法、装置、设备及超声成像探头
US11832992B2 (en) 2018-05-17 2023-12-05 Sonoscape Medical Corp. Shear wave propagation speed determination method and device
WO2019218669A1 (zh) * 2018-05-17 2019-11-21 深圳开立生物医疗科技股份有限公司 一种剪切波传播速度的确定方法及装置
CN108852416B (zh) * 2018-05-17 2021-03-30 深圳开立生物医疗科技股份有限公司 一种剪切波传播速度的确定方法及装置
CN108852416A (zh) * 2018-05-17 2018-11-23 深圳开立生物医疗科技股份有限公司 一种剪切波传播速度的确定方法及装置
CN109100801A (zh) * 2018-07-04 2018-12-28 中国科学院地质与地球物理研究所 一种高压介质地震波传播模拟方法
CN112584768A (zh) * 2018-07-26 2021-03-30 深圳迈瑞生物医疗电子股份有限公司 一种剪切波成像方法及系统
WO2020019254A1 (zh) * 2018-07-26 2020-01-30 深圳迈瑞生物医疗电子股份有限公司 一种剪切波成像方法及系统
CN111700639A (zh) * 2019-03-18 2020-09-25 三星麦迪森株式会社 超声诊断设备及其控制方法
CN110327077A (zh) * 2019-07-09 2019-10-15 深圳开立生物医疗科技股份有限公司 一种血流显示方法、装置及超声设备和存储介质
US11540809B2 (en) * 2019-07-12 2023-01-03 Konica Minolta, Inc. Ultrasonic diagnostic apparatus and method for propagation speed analysis of shear wave and elastic modulus measurement of a tissue
CN111671462A (zh) * 2020-06-18 2020-09-18 京东方科技集团股份有限公司 超声成像数据采集的扫描间隔计算方法及其设备、装置

Also Published As

Publication number Publication date
CN103462643B (zh) 2015-02-11

Similar Documents

Publication Publication Date Title
CN103462643B (zh) 一种剪切波速度测量方法、装置和系统
US11446006B2 (en) Adjusting measurements of the effects of acoustic radiation force for background motion effects
CN102809610B (zh) 一种基于改进的动态深度聚焦的相控阵超声检测方法
US10448924B2 (en) Ultrasonic acoustic radiation force excitation for ultrasonic material property measurement and imaging
CN103505243B (zh) 测量超声波的声吸收或衰减
CN103492855B (zh) 使用非聚焦超声的超声测振
CN102551791B (zh) 一种超声成像方法和装置
CN101852773B (zh) 超声波成像装置
CN107049361B (zh) 使用剪切波的声音速度成像
EP3153884A1 (en) Detection apparatus, fish finder, and radar
CN104605890B (zh) 一种剪切波峰值波形校正的方法、装置、系统及其应用
WO1998000719A2 (en) Apparatus and method for determining movements and velocities of moving objects
CN105877783B (zh) 二维剪切波弹性成像方法和装置
CN103800038A (zh) 改善的系统和装置以用于确定目标组织的机械特性
CN103149274B (zh) 一种混凝土缺陷检测方法
CN107003403A (zh) 声成像的方法和设备
CN1886094B (zh) 多普勒速度检测装置以及使用该装置的超声波诊断装置
CN110507359B (zh) 剪切波成像方法及系统
CN110244305A (zh) 一种水下目标信号散射的仿真方法
CN106443646B (zh) 一种超声波测距系统、回波处理方法及装置
US20040166478A1 (en) All ship speeds correlation SONAR simulator
Liang et al. An Ice Thickness Measurement Method based on Up-looking Sonar
Jensen et al. Improved accuracy in the estimation of blood velocity vectors using matched filtering
Laokulrat et al. Compensation of a transmitter's radiation pattern for synthetic transmit aperture acoustic imaging
JP6165024B2 (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
C56 Change in the name or address of the patentee

Owner name: SHENGZHEN SONOSCAPE MEDICAL CORP.

Free format text: FORMER NAME: SONOSCAPE INC.

CP03 Change of name, title or address

Address after: 518000 Guangdong city of Shenzhen province Nanshan District Yuquanlu Road Yizhe building 4, 5, 8, 9, 10 floor

Patentee after: Biologic medical Science and Technology Co., Ltd. is opened in Shenzhen

Address before: 518051 Guangdong city of Shenzhen province Nanshan District Yuquanlu Road Yizhe building 4 floor

Patentee before: Sonoscape Inc.