CN107049360B - 剪切波弹性成像方法和装置 - Google Patents
剪切波弹性成像方法和装置 Download PDFInfo
- Publication number
- CN107049360B CN107049360B CN201710061588.5A CN201710061588A CN107049360B CN 107049360 B CN107049360 B CN 107049360B CN 201710061588 A CN201710061588 A CN 201710061588A CN 107049360 B CN107049360 B CN 107049360B
- Authority
- CN
- China
- Prior art keywords
- shear wave
- interest
- region
- detected
- shear
- 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
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Radiology & Medical Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Biomedical Technology (AREA)
- Physics & Mathematics (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
技术领域
本发明涉及超声波医疗技术领域,更具体地涉及一种剪切波弹性成像方法和装置。
背景技术
超声成像由于具有实时、廉价、非侵入性和非电离辐射等优点而广泛地用于临床诊断。超声弹性成像,特别是基于声辐射力的剪切波成像,在组织弹性定性和定量的测量中发挥巨大的作用,例如,血管剪切波弹性成像对管壁的弹性定量起到至关重要的作用。然而由于一些组织本身的一些特性,导致相关技术中的剪切波弹性成像不能准确地对其进行弹性定性和定量。以血管为例来说,当血管短轴切面由于组织结构原因或是病理原因导致管壁分布不是圆形时,相关技术中的剪切波弹性成像将不能准确地得到血管壁的弹性模量,而不能准确地对血管壁的病理定性。
综上,需要提供一种剪切波弹性成像方法和装置,以至少部分地解决上述问题。
发明内容
考虑到上述问题而提出了本发明一种剪切波弹性成像方法和装置,可以针对一些由于组织结构原因或者病理原因导致其本身形状结构变形的待检测对象,准确地得到待检测对象的弹性模量,提高待检测对象的病理检测精度。
根据本发明一方面,提供了一种剪切波弹性成像方法,所述方法包括:获取待检测对象的感兴趣区域的超声图像数据;基于所述超声图像数据得到所述感兴趣区域的剪切波的二维位移矢量;对所述二维位移矢量的处理得到所述感兴趣区域的剪切波弹性成像数据,以得到所述感兴趣区域的弹性模量。
其中,所述方法还包括:对所述剪切波弹性成像数据进行可视化处理,以用于显示所述感兴趣区域的弹性模量。
其中,所述基于所述超声图像数据得到所述感兴趣区域的剪切波的二维位移矢量包括:采用块运动匹配的方法对所述超声图像数据进行处理,以得到所述感兴趣区域的剪切波传输的轴向位移;根据所述感兴趣区域的剪切波传输的轴向位移得到所述感兴趣区域的剪切波的二维位移矢量。
其中,所述根据所述感兴趣区域的剪切波传输的轴向位移得到所述感兴趣区域的剪切波的二维位移矢量包括:基于所述超声图像数据,利用空间向量复合计算不同角度的轴向位移,以得到所述感兴趣区域的剪切波的垂直位移和水平位移;所述感兴趣区域的剪切波的垂直位移和水平位移形成所述感兴趣区域的剪切波的二维位移矢量。
其中,利用空间向量符合计算不同角度的轴向位移采用如下公式:
其中,uver表示垂直位移,uhor表示水平位移,uax1至uaxn表示不同角度的轴向位移,θ1至θn表示不同角度,AT为A的转置矩阵。
其中,所述对所述二维位移矢量的处理得到所述感兴趣区域的剪切波弹性成像数据包括:根据所述二维位移矢量得到所述剪切波的传播方向;基于所述剪切波的传播方向进行坐标变换,以将所述剪切波从直角坐标系转换为极坐标系;在所述极坐标系下,对所述剪切波的传播方向进行方向滤波,以得到沿待检测对象的环向传播的剪切波;通过对所述沿待检测对象的环向传播的剪切波得到所述感兴趣区域的剪切波弹性成像数据,以得到所述感兴趣区域的弹性模量。
其中,所述通过对所述沿待检测对象的环向传播的剪切波得到所述感兴趣区域的剪切波弹性成像数据包括:在极坐标系下,计算沿待检测对象的环向传播的剪切波的传播速度;以及获取沿待检测对象的环向传播的剪切波的频散曲线;其中,所述传播速度和所述频散曲线形成所述感兴趣区域的剪切波弹性成像数据,用于表征感兴趣区域的弹性模量。
其中,所述获取沿待检测对象的环向传播的剪切波的频散曲线包括:在极坐标系中,根据沿待检测对象的环向传播的剪切波的波前在不同时间到达位置得到的位移时间曲线;对所述位移时间曲线进行二维傅里叶变换,得到傅里叶变换结果,其中,所述傅里叶变换结果包括多个频率和每个频率对应的多个波数值;通过每个频率下最大的波数值,得到待检测对象的剪切波的频散曲线。
根据本发明另一方面,提供了一种剪切波弹性成像装置,所述装置包括:
数据获取单元,用于获取待检测对象的感兴趣区域的超声图像数据;二维位移矢量计算单元,用于所述超声图像数据得到所述感兴趣区域的剪切波的二维位移矢量;剪切波弹性成像单元,用于对所述二维位移矢量的处理得到所述感兴趣区域的剪切波弹性成像数据,以得到所述感兴趣区域的弹性模量。
其中,所述装置还包括:可视化处理单元,用于对所述剪切波弹性成像数据进行可视化处理,以用于显示所述感兴趣区域的弹性模量。
本发明实施例提供的剪切波弹性成像方法和装置,利用二维位移矢量估计,得到剪切波沿待检测对象(例如,血管壁)的径向的运动,为待检测对象的弹性的定量研究提供基础,最大限度地提高待检测对象的弹性模量的准确度,提高对待检测对象的病理检测精度。不仅适用于形状规则的待检测对象,对于受周围组织或疾病(例如斑块)影响导致待检测对象(例如血管)形状不规则的情况,也可以用本发明的方法和装置进行基于剪切波的弹性定量,适用性较高,适用范围更广。
附图说明
通过结合附图对本发明实施例进行更详细的描述,本发明的上述以及其它目的、特征和优势将变得更加明显。附图用来提供对本发明实施例的进一步理解,并且构成说明书的一部分,与本发明实施例一起用于解释本发明,并不构成对本发明的限制。在附图中,相同的参考标号通常代表相同部件或步骤。
图1示出根据本发明一个实施例的剪切波弹性成像方法的流程示意图;
图2示出根据本发明一个实施例的空间向量复合计算的原理示意图;
图3示出根据本发明另一个实施例的剪切波弹性成像方法的流程示意图;
图4示出根据本发明一个实施例的剪切波弹性成像装置的结构框图。
具体实施方式
为了使得本发明的目的、技术方案和优点更为明显,下面将参照附图详细描述根据本发明的示例实施例。显然,所描述的实施例仅仅是本发明的一部分实施例,而不是本发明的全部实施例,应理解,本发明不受这里描述的示例实施例的限制。基于本发明中描述的本发明实施例,本领域技术人员在没有付出创造性劳动的情况下所得到的所有其它实施例都应落入本发明的保护范围之内。
基于声辐射力的剪切波弹性成像技术是一种评估组织硬度的超声弹性成像技术。其基本原理为:由探头向生物体的软组织发射高能量的超声波,在声辐射力和组织的剪切应力的作用下,特定区域内的软组织会产生向四周传播的振动,从而产生剪切波,由于生物体软组织的硬度和剪切波速度存在着关联关系,因此可以通过检测剪切波的速度以分析生物体软组织的硬度。
相关技术的剪切波弹性成像的可以包括两个关键点:1、考虑偏振方向,关心偏振方向的运动,即剪切波的振动方向;2、方向滤波器设置时仅考虑某一个感兴趣的方向。然而由于一些组织本身的一些特性,导致相关技术中的剪切波弹性成像不能准确地对其进行弹性定性和定量。以血管为例来说,当血管短轴切面由于组织结构原因或是病理原因导致管壁分布不是圆形时,相关技术中的剪切波弹性成像只对其长轴切面进行处理,并不能准确地得到血管壁的弹性模量,因此,很难准确地对血管壁的病理定性。
为了解决上文所述的问题,本发明实施例提出一种剪切波弹性成像方法和装置,利用二维位移矢量估计,得到剪切波沿待检测对象(例如,血管壁)的径向的运动,为一些产生形变的待检测对象的弹性的定量研究提供基础,最大限度地提高待检测对象的弹性模量的准确度,提高对待检测对象的病理检测精度。
下面将结合附图对本发明提供的剪切波弹性成像方法和装置进行详细描述,以使本领域技术人员能够清楚、准确地理解本发明的技术方案。
图1示出根据本发明一个实施例的剪切波弹性成像方法的流程示意图。
本发明实施例提供的一种剪切波弹性成像方法包括如下步骤:
步骤110,获取待检测对象的感兴趣区域的超声图像数据。
在本步骤中,可以利用声辐射力激励待检测对象的感兴趣区域,对感兴趣区域进行超高速超声成像,以得到待检测对象的感兴趣区域的超声图像数据。示例性地,超声图像数据可以包括RF(Radio Frequency,射频)数据、包络数据、B模数据中的至少一种。
以血管为例,血管内弹性成像可以利用气囊、血压变化或者外部挤压来激励血管,估计血管的运动即位移(一般为纵向),得到血管的应变分布,从而表征血管的弹性。在本发明实施例中,可以利用声辐射力激励血管,得到血管的超声图像数据。示例性地,在本步骤中获得的超声图像数据为直角坐标系的超声图像数据。
步骤120,基于超声图像数据得到感兴趣区域的剪切波的二维位移矢量。
在本步骤中,可以通过如下方式得到感兴趣区域的剪切波的二维位移矢量:采用块运动匹配的方法对超声图像数据进行处理,以得到感兴趣区域的剪切波传输的轴向位移,然后根据感兴趣区域的剪切波传输的轴向位移得到感兴趣区域的剪切波的二维位移矢量。
其中,块运动匹配的方法例如可以包括互相关、自相关、光流法等。采用块运动匹配方法得到轴向位移例如可以实现为:对超声图像数据中的任一数据点进行加窗处理,得到对应该数据点的数据窗,当连续两帧的超声图像数据中数据窗的相关性(评价相关性的参数诸多,例如,可以是自相关系数,也可以是互相关系数等)最大时,确定该数据点的在极坐标系中的横向位移和轴向位移。在本发明实施例中,相关性系数可以通过自相关、互相关、绝对误差和幅度差平方和以及光流法等方式计算确定,可以理解,相关性系数的计算不限于这几种实现方式,还可以通过其他块匹配的运动估计算法实现,在此不再一一赘述。
而根据轴向位移得到感兴趣区域的剪切波的二维位移矢量可以实现为:利用空间向量复合计算不同角度(或者不同角度不同时刻)的轴向位移,以得到感兴趣区域的剪切波的垂直位移和水平位移,由感兴趣区域的剪切波的垂直位移和水平位移形成感兴趣区域的剪切波的二维位移矢量。如图2所示,为空间向量复合计算的原理示意图。在图2中,示意了θ1至θ2两个角度方向,假设已通过块运动匹配计算得到两个角度方向的横向位移ulat1和ulat2,以及轴向位移uax1和uax2,因此可以通过如下公式进行空间向量复合计算,以得到垂直位移uver和水平位移uhor:
uver=uax cosθ+ulat sinθ 式(1)
uhor=uax(-sinθ)+ulat cosθ 式(2)
uver=uax1cosθ1+ulat1sinθ1 式(3)
uhor=uax1(-sinθ1)+ulat1cosθ1 式(4)
uver=uax2cosθ2+ulat2sinθ2 式(5)
uhor=uax2(-sinθ2)+ulat2cosθ2 式(6)
在空间向量复合计算以进行剪切波弹性成像的条件下,超声回波信号会沿着偏转角方向传播,因此通过块运动匹配方法可以得到其传播的轴向位移和侧向位移,基于轴向位移和侧向位移,可以通过公式(1)和公式(2)计算得到剪切波的垂直位移和水平位移。但是由于超声成像过程中轴向采样率很高而侧向采样率远低于轴向采样率,因此,超声成像的侧向分辨率比较低,使得侧向位移估计质量远低于轴向位移,基于此,通过两个角度θ1至θ2得到不同的轴向位移和横向位移估计结果,根据公式(3)至公式(6)消去横向位移得到公式(7)和公式(8),即仅通过轴向位移复合得出垂直位移和水平位移。
同理,多个角度空间向量复合时,可以通过提高平均次数提高位移估计的质量,因此,在本发明实施例中,采用如下公式计算复合得到垂直位移和水平位移:
其中,uver表示垂直位移,uhor表示水平位移,uax1至uaxn表示不同角度的轴向位移,θ1至θn表示不同角度,AT为A的转置矩阵。
根据本发明另一实施例,当待检测对象(例如,血管)为规则或近似规则的圆形时,可以认为剪切波的运动是沿着半径方向传播,可以通过手动确定圆心或者通过机器学习方法确定圆心,并由轴向位移和横向位移得到径向位移,进而基于径向位移得到剪切波弹性成像数据,得到待检测对象的感兴趣区域的弹性模量。
步骤130,对二维位移矢量的处理得到感兴趣区域的剪切波弹性成像数据,以得到感兴趣区域的弹性模量。
可以通过如下步骤实现本步骤130:根据二维位移矢量得到剪切波的传播方向,基于剪切波的传播方向进行坐标变换,以将剪切波从直角坐标系转换为极坐标系,在极坐标系下,对剪切波的传播方向进行方向滤波,以得到沿待检测对象的环向传播的剪切波,通过对沿待检测对象的环向传播的剪切波得到感兴趣区域的剪切波弹性成像数据,以得到感兴趣区域的弹性模量。
基于剪切波的传播方向进行坐标变换,以将剪切波从直角坐标系转换为极坐标系例如可以实现为:通过空间插值法实现直角坐标系至极坐标系的坐标转换。作为一个示例,可以将待检测对象的几何中心确定为极坐标原点,根据预设采样率分别对极坐标系的角度方向和半径方向进行数据采样,以获得采样数据点,用直角坐标系的超声图像数据对采用数据点进行空间插值,以得到极坐标系的超声图像数据。示例性地,在极坐标系中,横坐标为角度,纵坐标为半径。关于确定待检测对象的几何中心,例如可以通过手动选择的方式确定,也可以通过对待检测对象的图像进行图像自动分割以实现对其几何中心的选择。
仍以血管为例,将血管的几何中心(例如,圆心)作为极坐标原点,从原点发射线状采样,其中,沿着射线方向为半径方向,垂直于射线方向为角度方向。在本发明实施例中,根据预设采样率对血管的半径方向和角度方向均进行采样,以获得采用数据点,再利用直角坐标系的超声图像数据对采用数据点进行空间插值,实现直角坐标系到极坐标系的坐标转换,从而得到极坐标系的超声图像数据。预设采样率例如可以是极坐标系的任一个方向采集的数据点的数量,其可以根据实际需求设定,在本发明实施例中不对此进行限定。空间插值法可以例如采用线性插值、多项式拟合插值等任一种实现,本发明不对此进行限制。
以及在极坐标系下,对剪切波的传播方向进行方向滤波,以得到沿待检测对象的环向传播的剪切波可以实现为:基于垂直位移和水平位移建立位移时间曲线,示例性地,例如可以在待检测对象的感兴趣区域的传播路径上,任意一个时刻对应波前到达的空间位置,将连续的若干个时刻整合在一起,一个维度代表时间,另一个维度代表空间位置,由于波前传播是连续的,因此可以在二维空间内确定波前传播的情况,即建立位移时间曲线。然后对位移时间曲线进行二维傅里叶变换,以得到傅里叶变换结果,其中,傅里叶变换结果的不同角度范围表示剪切波的不同方向的运动信息。最后,对感兴趣的角度进行加窗处理,后经过逆傅里叶变换,以得到沿待检测对象的感兴趣方向传播的剪切波,即实现对剪切波的传播方向的方向滤波,得到感兴趣方向传播的剪切波。
进一步地,通过对沿待检测对象的环向传播的剪切波得到感兴趣区域的剪切波弹性成像数据可以实现为:在极坐标系下,计算沿待检测对象的环向传播的剪切波的传播速度,以及获取沿待检测对象的环向传播的剪切波的频散曲线,传播速度和频散曲线形成感兴趣区域的剪切波弹性成像数据,用于表征感兴趣区域的弹性模量。
具体地,可以基于任一数据点的位移,得到剪切波波前在不同时间到达位置得到的位移时间曲线,对位移时间曲线进行二维傅里叶变换,以得到傅里叶变换结果,傅里叶变换结果中包括频率值和任一频率值对应的所有波数值,针对感兴趣的角度,基于感兴趣的角度中每个频率值对应的最大波数值建立剪切波频散曲线。示例性地,选择感兴趣的角度方向,例如,当感兴趣的角度方向向右时,对应的傅里叶变换中的信息位于第一象限,其中,傅里叶变换中的信息分别包括频率信息和波数信息(频率值和波数值),取任一个频率值下最大的波数值,得到待检测对象的剪切波频散曲线,即不同频率的剪切波相速度曲线。基于剪切波的频散曲线评价待检测对象的剪切模量,准确率更高,而对于像血管这样的待检测对象,基于极坐标系下对感兴趣的角度的选择,可以同时处理长轴切面和短轴切面,弥补相关技术的弊端,得到更精准的待检测对象的弹性定量。
而在本发明实施例中,将频散曲线和传播速度相结合共同对待检测对象进行定性和定量。示例性地,在极坐标系中,根据剪切波波前在不同时间到达位置对剪切波波速进行测量,从而得到剪切波的传播速度。
综上,结合剪切波的传播速度及剪切波的频散曲线评价待检测对象的剪切模量,两者可以相辅相成,进一步提高待检测对象的弹性定量的精准度。
应用本发明实施例,利用二维位移矢量估计,得到剪切波沿待检测对象(例如,血管壁)的径向的运动,为待检测对象的弹性的定量研究提供基础,最大限度地提高待检测对象的弹性模量的准确度,提高对待检测对象的病理检测精度。
图3示出根据本发明另一个实施例的剪切波弹性成像方法的流程示意图。
如图3所示,一种剪切波弹性成像方法可以包括如下步骤:
步骤110,获取待检测对象的感兴趣区域的超声图像数据。
步骤120,基于超声图像数据得到感兴趣区域的剪切波的二维位移矢量。
步骤130,对二维位移矢量的处理得到感兴趣区域的剪切波弹性成像数据,以得到感兴趣区域的弹性模量。
其中,步骤110至步骤130的具体实现方式及技术效果可以参考图1所示实施例,为了简洁,在此不再对步骤110至步骤130进行详细叙述。
步骤140,对剪切波弹性成像数据进行可视化处理,以用于显示感兴趣区域的弹性模量。
在本步骤中,由于方向滤波是在极坐标下进行,得到方向滤波后的剪切波传播后,可以对剪切波传播的数据(例如,可以包括散射子震动位移、速度、加速度等)进行坐标变换,即由极坐标转换成直角坐标,转换原理及过程与步骤130中涉及到的坐标转换的原理相同,在此不再赘述,通过坐标转换实现剪切波传播的可视化,即实现剪切波弹性成像的可视化处理。
应用本发明实施例提供的剪切波弹性成像方法,采用二维位移估计的方法来得到剪切波激励后各个时刻的二维位移矢量,并通过坐标变换来解决方向滤波的问题,且通过坐标转换完成剪切波传播的可视化,进而得到剪切波传播速度、频散曲线等,为待检测对象,尤其是血管壁的弹性的定量,特别是各向异性信息的提取打下坚实基础。
并且,本发明的剪切波弹性成像方法不仅适用于形状规则的待检测对象,对于受周围组织或疾病(例如斑块)影响导致待检测对象(例如血管)形状不规则的情况,也可以用本发明的方法进行基于剪切波的弹性定量,适用性较高,适用范围更广。
另外,通过待检测对象(例如血管)横截面的剪切波成像,进一步可以得到管壁的剪切波速度、频散曲线,结合剪切波传播速度和频散曲线对待检测对象的弹性进行定量测量,可以提高测量精度和准确度。
本发明实施例还提供了一种应用上述剪切波成像方法的剪切波成像装置,下面将结合图4对剪切波弹性成像装置进行详述。
图4示出根据本发明一个实施例的剪切波弹性成像装置的结构框图。
如图4所示,剪切波弹性成像装置400可以包括数据获取单元410、二维位移矢量计算单元420和剪切波弹性成像单元430。
数据获取单元410可以用于获取待检测对象的感兴趣区域的超声图像数据。
二维位移矢量计算单元420可以用于超声图像数据得到感兴趣区域的剪切波的二维位移矢量。
剪切波弹性成像单元430可以用于对二维位移矢量的处理得到感兴趣区域的剪切波弹性成像数据,以得到感兴趣区域的弹性模量。
根据本发明另一实施例,剪切波弹性成像装置400还可以包括可视化处理单元(图中未示出),其可以用于对剪切波弹性成像数据进行可视化处理,以用于显示感兴趣区域的弹性模量。
需要说明的是,剪切波弹性成像装置400中各个模块的功能及配置结构可以参考前述剪切波弹性成像方法的对应步骤的实现过程,在此不再赘述。
应用本发明实施例,利用二维位移矢量估计,得到剪切波沿待检测对象(例如,血管壁)的径向的运动,为待检测对象的弹性的定量研究提供基础,最大限度地提高待检测对象的弹性模量的准确度,提高对待检测对象的病理检测精度,且不仅适用于形状规则的待检测对象,对于受周围组织或疾病(例如斑块)影响导致待检测对象(例如血管)形状不规则的情况,也可以用本发明的方法进行基于剪切波的弹性定量,适用性较高,适用范围更广。
尽管本文已经参考附图描述了示例实施例,应理解上述示例实施例仅仅是示例性的,并且不意图将本发明的范围限制于此。本领域普通技术人员可以在其中进行各种改变和修改,而不偏离本发明的范围和精神。所有这些改变和修改意在被包括在所附权利要求所要求的本发明的范围之内。
本领域普通技术人员可以意识到,结合本文中所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本发明的范围。
在本申请所提供的几个实施例中,应该理解到,所揭露的设备和方法,可以通过其它的方式实现。例如,以上所描述的设备实施例仅仅是示意性的,例如,单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个设备,或一些特征可以忽略,或不执行。
在此处所提供的说明书中,说明了大量具体细节。然而,能够理解,本发明的实施例可以在没有这些具体细节的情况下实践。在一些实例中,并未详细示出公知的方法、结构和技术,以便不模糊对本说明书的理解。
类似地,应当理解,为了精简本发明并帮助理解各个发明方面中的一个或多个,在对本发明的示例性实施例的描述中,本发明的各个特征有时被一起分组到单个实施例、图、或者对其的描述中。然而,并不应将该本发明的方法解释成反映如下意图:即所要求保护的本发明要求比在每个权利要求中所明确记载的特征更多的特征。更确切地说,如相应的权利要求书所反映的那样,其发明点在于可以用少于某个公开的单个实施例的所有特征的特征来解决相应的技术问题。因此,遵循具体实施方式的权利要求书由此明确地并入该具体实施方式,其中每个权利要求本身都作为本发明的单独实施例。
本领域的技术人员可以理解,除了特征之间相互排斥之外,可以采用任何组合对本说明书(包括伴随的权利要求、摘要和附图)中公开的所有特征以及如此公开的任何方法或者设备的所有过程或单元进行组合。除非另外明确陈述,本说明书(包括伴随的权利要求、摘要和附图)中公开的每个特征可以由提供相同、等同或相似目的的替代特征来代替。
此外,本领域的技术人员能够理解,尽管在此所述的一些实施例包括其它实施例中所包括的某些特征而不是其它特征,但是不同实施例的特征的组合意味着处于本发明的范围之内并且形成不同的实施例。例如,在权利要求书中,所要求保护的实施例的任意之一都可以以任意的组合方式来使用。
本发明的各个部件实施例可以以硬件实现,或者以在一个或者多个处理器上运行的软件模块实现,或者以它们的组合实现。本领域的技术人员应当理解,可以在实践中使用微处理器或者数字信号处理器(DSP)来实现根据本发明实施例的导线悬挂点定位装置中的一些模块的一些或者全部功能。本发明还可以实现为用于执行这里所描述的方法的一部分或者全部的装置程序(例如,计算机程序和计算机程序产品)。这样的实现本发明的程序可以存储在计算机可读介质上,或者可以具有一个或者多个信号的形式。这样的信号可以从因特网网站上下载得到,或者在载体信号上提供,或者以任何其他形式提供。
应该注意的是上述实施例对本发明进行说明而不是对本发明进行限制,并且本领域技术人员在不脱离所附权利要求的范围的情况下可设计出替换实施例。在权利要求中,不应将位于括号之间的任何参考符号构造成对权利要求的限制。单词“包含”不排除存在未列在权利要求中的组件或步骤。位于组件之前的单词“一”或“一个”不排除存在多个这样的组件。本发明可以借助于包括有若干不同组件的硬件以及借助于适当编程的计算机来实现。在列举了若干装置的单元权利要求中,这些装置中的若干个可以是通过同一个硬件项来具体体现。单词第一、第二、以及第三等的使用不表示任何顺序。可将这些单词解释为名称。
以上所述,仅为本发明的具体实施方式或对具体实施方式的说明,本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。本发明的保护范围应以权利要求的保护范围为准。
Claims (9)
1.一种剪切波弹性成像方法,其特征在于,所述方法包括:
获取待检测对象的感兴趣区域的超声图像数据;
基于所述超声图像数据得到所述感兴趣区域的剪切波的二维位移矢量;
对所述二维位移矢量的处理得到所述感兴趣区域的剪切波弹性成像数据,以得到所述感兴趣区域的弹性模量;
其中,所述对所述二维位移矢量的处理得到所述感兴趣区域的剪切波弹性成像数据包括:
根据所述二维位移矢量得到所述剪切波的传播方向;
基于所述剪切波的传播方向进行坐标变换,以将所述剪切波从直角坐标系转换为极坐标系;
在所述极坐标系下,对所述剪切波的传播方向进行方向滤波,以得到沿待检测对象的环向传播的剪切波;
通过所述沿待检测对象的环向传播的剪切波得到所述感兴趣区域的剪切波弹性成像数据,以得到所述感兴趣区域的弹性模量。
2.根据权利要求1所述的方法,其特征在于,所述方法还包括:
对所述剪切波弹性成像数据进行可视化处理,以用于显示所述感兴趣区域的弹性模量。
3.根据权利要求1所述的方法,其特征在于,所述基于所述超声图像数据得到所述感兴趣区域的剪切波的二维位移矢量包括:
采用块运动匹配的方法对所述超声图像数据进行处理,以得到所述感兴趣区域的剪切波传输的轴向位移;
根据所述感兴趣区域的剪切波传输的轴向位移得到所述感兴趣区域的剪切波的二维位移矢量。
4.根据权利要求3所述的方法,其特征在于,所述根据所述感兴趣区域的剪切波传输的轴向位移得到所述感兴趣区域的剪切波的二维位移矢量包括:
基于所述超声图像数据,利用空间向量复合计算不同角度的轴向位移,以得到所述感兴趣区域的剪切波的垂直位移和水平位移;
所述感兴趣区域的剪切波的垂直位移和水平位移形成所述感兴趣区域的剪切波的二维位移矢量。
6.根据权利要求1所述的方法,其特征在于,所述通过对所述沿待检测对象的环向传播的剪切波得到所述感兴趣区域的剪切波弹性成像数据包括:
在极坐标系下,计算沿待检测对象的环向传播的剪切波的传播速度;
以及获取沿待检测对象的环向传播的剪切波的频散曲线;
其中,所述传播速度和所述频散曲线形成所述感兴趣区域的剪切波弹性成像数据,用于表征感兴趣区域的弹性模量。
7.根据权利要求6所述的方法,其特征在于,所述获取沿待检测对象的环向传播的剪切波的频散曲线包括:
在极坐标系中,根据沿待检测对象的环向传播的剪切波的波前在不同时间到达位置得到的位移时间曲线;
对所述位移时间曲线进行二维傅里叶变换,得到傅里叶变换结果,其中,所述傅里叶变换结果包括多个频率和每个频率对应的多个波数值;
通过每个频率下最大的波数值,得到待检测对象的剪切波的频散曲线。
8.一种剪切波弹性成像装置,其特征在于,所述装置包括:
数据获取单元,用于获取待检测对象的感兴趣区域的超声图像数据;
二维位移矢量计算单元,用于所述超声图像数据得到所述感兴趣区域的剪切波的二维位移矢量;
剪切波弹性成像单元,用于根据所述二维位移矢量得到所述剪切波的传播方向;基于所述剪切波的传播方向进行坐标变换,以将所述剪切波从直角坐标系转换为极坐标系;在所述极坐标系下,对所述剪切波的传播方向进行方向滤波,以得到沿待检测对象的环向传播的剪切波;通过所述沿待检测对象的环向传播的剪切波得到所述感兴趣区域的剪切波弹性成像数据,以得到所述感兴趣区域的弹性模量。
9.根据权利要求8所述的装置,其特征在于,所述装置还包括:
可视化处理单元,用于对所述剪切波弹性成像数据进行可视化处理,以用于显示所述感兴趣区域的弹性模量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710061588.5A CN107049360B (zh) | 2017-01-26 | 2017-01-26 | 剪切波弹性成像方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710061588.5A CN107049360B (zh) | 2017-01-26 | 2017-01-26 | 剪切波弹性成像方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107049360A CN107049360A (zh) | 2017-08-18 |
CN107049360B true CN107049360B (zh) | 2020-12-08 |
Family
ID=59598700
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710061588.5A Active CN107049360B (zh) | 2017-01-26 | 2017-01-26 | 剪切波弹性成像方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107049360B (zh) |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109171814B (zh) * | 2018-07-26 | 2020-08-28 | 清华大学 | 血管超声弹性成像的方法及装置 |
CN109259801B (zh) * | 2018-09-12 | 2021-05-14 | 深圳开立生物医疗科技股份有限公司 | 一种剪切波弹性成像方法及装置 |
WO2020077598A1 (zh) * | 2018-10-18 | 2020-04-23 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声弹性检测方法及其系统 |
CN113229846B (zh) * | 2018-12-24 | 2023-10-13 | 深圳迈瑞生物医疗电子股份有限公司 | 超声成像方法与设备 |
CN111399051B (zh) * | 2019-01-02 | 2021-10-08 | 无锡海斯凯尔医学技术有限公司 | 数据处理方法、装置及介质 |
WO2020144077A1 (en) * | 2019-01-09 | 2020-07-16 | Koninklijke Philips N.V. | Ultrasound system and method for shear wave characterization of anisotropic tissue |
CN109919918A (zh) * | 2019-02-21 | 2019-06-21 | 清华大学 | 弹性成像的控制方法及装置、计算机设备及可读存储介质 |
CN109893172B (zh) * | 2019-02-22 | 2020-06-19 | 清华大学 | 基于弹性成像的力学参数的确定方法及装置、计算机设备 |
JP7236312B2 (ja) * | 2019-04-04 | 2023-03-09 | 富士フイルムヘルスケア株式会社 | 超音波診断装置、信号処理装置、および、プログラム |
CN112438751A (zh) * | 2019-09-04 | 2021-03-05 | 通用电气精准医疗有限责任公司 | 用于剪切波弹性成像的方法和系统以及存储有对应程序的介质 |
CN110742644B (zh) * | 2019-09-29 | 2022-07-08 | 深圳大学 | 弹性成像系统、弹性成像方法及存储介质 |
CN110742645B (zh) * | 2019-09-29 | 2022-09-27 | 深圳大学 | 多模成像系统、多模成像方法及存储介质 |
CN110811689B (zh) * | 2019-10-31 | 2020-11-27 | 汕头市超声仪器研究所股份有限公司 | 一种剪切波运动速度的一阶估算方法 |
CN111275706B (zh) * | 2020-03-04 | 2023-06-02 | 中山大学附属第一医院 | 一种基于剪切波弹性成像超声组学深度分析方法及系统 |
CN111449681B (zh) * | 2020-04-08 | 2023-09-08 | 深圳开立生物医疗科技股份有限公司 | 一种剪切波成像方法、装置、设备及可读存储介质 |
WO2021212494A1 (zh) * | 2020-04-24 | 2021-10-28 | 深圳迈瑞生物医疗电子股份有限公司 | 粘弹性测量方法和超声测量系统 |
CN113081039B (zh) * | 2021-04-06 | 2024-02-13 | 无锡海斯凯尔医学技术有限公司 | 弹性成像方法、装置、电子设备及存储介质 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003530941A (ja) * | 2000-04-26 | 2003-10-21 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 剪断波パラメター評価のための超音波方法およびシステム |
JP5411699B2 (ja) * | 2007-09-06 | 2014-02-12 | 株式会社日立メディコ | 超音波撮像装置 |
US8197410B2 (en) * | 2007-09-07 | 2012-06-12 | Kabushiki Kaisha Toshiba | Ultrasonic diagnostic apparatus, ultrasonic image processing apparatus and ultrasonic image processing method |
CN102423264B (zh) * | 2011-09-01 | 2014-05-21 | 中国科学院深圳先进技术研究院 | 基于图像的生物组织弹性的测量方法及装置 |
US20150094580A1 (en) * | 2012-04-13 | 2015-04-02 | Hitachi Aloka Medical, Ltd. | Ultrasonic diagnostic device and locus display method |
CN104055540B (zh) * | 2014-05-09 | 2016-02-24 | 西安交通大学 | 一种无创高精度血管壁弹性成像的方法 |
-
2017
- 2017-01-26 CN CN201710061588.5A patent/CN107049360B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN107049360A (zh) | 2017-08-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107049360B (zh) | 剪切波弹性成像方法和装置 | |
Deng et al. | Ultrasonic shear wave elasticity imaging sequencing and data processing using a verasonics research scanner | |
JP4258015B2 (ja) | 超音波診断システム、歪み分布表示方法及び弾性係数分布表示方法 | |
CN106618635B (zh) | 剪切波弹性成像方法和装置 | |
EP2889004B1 (en) | Motion correction in three-dimensional elasticity ultrasound imaging | |
KR102025328B1 (ko) | 평면파 합성을 이용한 초음파 벡터 도플러 영상의 생성 장치 및 방법 | |
EP3082616B1 (en) | Computed ultrasound tomography in echo mode (cute) for imaging speed of sound using pulse-echo sonography | |
JP2013000349A (ja) | イメージング方法及び変位計測方法及び装置、並びに、超音波画像診断装置 | |
JP4221555B2 (ja) | 超音波診断システム、歪み分布表示方法及び弾性係数分布表示方法 | |
JP6564059B2 (ja) | 組織形態及び弾性情報の処理方法、並びに弾性検出装置 | |
US20150272551A1 (en) | Module for Processing Ultrasonic Signal Based on Spatial Coherence and Method for Processing Ultrasonic Signal | |
CN114176639A (zh) | 用于介质的超声表征的方法和系统 | |
JP4389091B2 (ja) | 超音波診断システム、歪み分布表示方法及び弾性係数分布表示方法 | |
Yuan et al. | Analytical phase-tracking-based strain estimation for ultrasound elasticity | |
US20190076126A1 (en) | Shear wave group velocity estimation using spatiotemporal peaks and amplitude thresholding | |
Engel et al. | A new method for shear wave speed estimation in shear wave elastography | |
CN114176640A (zh) | 用于介质的超声表征的方法和系统 | |
CN114176626A (zh) | 用于介质的超声表征的方法和系统 | |
Lu et al. | A real time displacement estimation algorithm for ultrasound elastography | |
Fredriksen et al. | Investigations of spectral resolution and angle dependency in a 2-D tracking Doppler method | |
CN107530051B (zh) | 使用向量速度超声(us)对血管内压变化的无创估计 | |
CN114176625A (zh) | 用于介质的超声表征的方法和系统 | |
KR101117544B1 (ko) | 탄성영상 형성을 위한 초음파 진단 시스템 및 방법 | |
Jiang et al. | Modified phase zero method for ultrasound freehand strain imaging | |
Zhao et al. | Liver elasticity imaging using external vibration multi-directional ultrasound shearwave elastography (EVMUSE) |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |