CN111449681A - 一种剪切波成像方法、装置、设备及可读存储介质 - Google Patents

一种剪切波成像方法、装置、设备及可读存储介质 Download PDF

Info

Publication number
CN111449681A
CN111449681A CN202010268506.6A CN202010268506A CN111449681A CN 111449681 A CN111449681 A CN 111449681A CN 202010268506 A CN202010268506 A CN 202010268506A CN 111449681 A CN111449681 A CN 111449681A
Authority
CN
China
Prior art keywords
shear wave
data
acceleration
scanning line
wave velocity
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
CN202010268506.6A
Other languages
English (en)
Other versions
CN111449681B (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 Medical Corp
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 Medical Corp filed Critical Sonoscape Medical Corp
Priority to CN202010268506.6A priority Critical patent/CN111449681B/zh
Publication of CN111449681A publication Critical patent/CN111449681A/zh
Application granted granted Critical
Publication of CN111449681B publication Critical patent/CN111449681B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/48—Diagnostic techniques
    • A61B8/485—Diagnostic techniques involving measuring strain or elastic properties
    • 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
    • A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30—Assessment of water resources

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)
  • Veterinary Medicine (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)
  • Computer Vision & Pattern Recognition (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本发明公开了一种剪切波成像方法、装置、设备及可读存储介质，该方法包括以下步骤：获取剪切波速度数据；对剪切波速度数据进行一阶导处理，得到剪切波加速度数据；利用剪切波加速度数据生成剪切波加速度图像。在本方法中，在获得剪切波速度数据之后，对剪切波速度数据进行一阶导处理，得到剪切波加速度数据，然后基于剪切波加速度数据便可生成剪切波加速度图像。即基于剪切波加速度图像可确定出软硬不一的区域分界，也就是说，可以得到组织弹性变化边界。

Description

一种剪切波成像方法、装置、设备及可读存储介质
技术领域
本发明涉及超声成像技术领域，特别是涉及一种剪切波成像方法、装置、设备及可读存储介质。
背景技术
基于声辐射力的剪切波弹性成像技术是一种评估组织弹性的超声弹性成像技术，其被广泛应用于人体组织病变的分析和诊断中。基于声辐射力的剪切波弹性成像技术的原理为：由探头向人体的软组织发射高能量的超声波后，在声辐射力和组织的剪切应力的作用下，特定区域内的软组织会产生向四周传播的振动，从而产生剪切波，而由于人体软组织的弹性与剪切波的速度之间存在着密切的关系，进而可以基于剪切波的速度确定是否存在软组织的病变。
由于人体组织错综复杂，包括大量的软硬度不一的区域。剪切波在传播的过程中遇到软硬度不一的区域时，剪切波会发生反射、折射等现象，从而导致剪切波弹性图像在分界处得到错误的结果，无法准确地定位出病变区域的边界。
综上所述，如何有效地提高剪切波成像质量，是目前本领域技术人员急需解决的技术问题。
发明内容
本发明的目的是提供一种剪切波成像方法、装置、设备及可读存储介质，通过计算剪切波的加速度，并利用剪切波加速度精准定位软组织软硬度变化的区域边缘，使得剪切波弹性图像在分界处得到准确的结果，进而准确地定位出病变区域的边界。
为解决上述技术问题，本发明提供如下技术方案：
一方面提供了一种剪切波成像方法，包括：
获取剪切波速度数据；
对所述剪切波速度数据进行一阶导处理，得到剪切波加速度数据；
利用所述剪切波加速度数据生成剪切波加速度图像。
优选地，还包括：
对所述剪切波加速度数据和剪切波速度数据进行融合成像处理，得到剪切波弹性图像。
优选地，所述获取剪切波速度数据，包括：
向目标组织发射超声波，并采集超声回波数据；所述超声回波数据中包括剪切波传输信息；
对所述超声回波数据进行处理，得到所述剪切波速度数据。
优选地，在对所述超声回波数据进行处理之前，还包括：
利用方向滤波器滤除干扰数据，得到无干扰的超声回波数据。
优选地，对所述超声回波数据进行处理，得到所述剪切波速度数据，包括：
利用互相关算法或时间飞跃算法计算出剪切波经过每个目标扫描线的前后扫描线对应的所述时间间隔；
计算每个目标扫描线的前后扫描线的距离与所述时间间隔的比值，得到剪切波速度；
将所述剪切波速度和所述时间间隔作为所述剪切波速度数据。
优选地，对所述剪切波速度数据进行一阶导处理，得到剪切波加速度数据，包括：
计算出剪切波经过每个所述目标扫描线的前后扫描线对应的速度差值；
计算所述速度差值与所述时间间隔的比值，得到所述目标扫描线的所述加速度。
优选地，利用所述剪切波加速度数据生成剪切波加速度图像，包括：
利用灰度映射或彩色映射，对所述剪切波加速度数据进行线性映射，得到所述剪切波加速度图像。
另一方面，提供了一种剪切波成像装置，包括：
剪切波速度获取模块，用于获取剪切波速度数据；
剪切波加速度获取模块，用于对所述剪切波速度数据进行一阶导处理，得到剪切波加速度数据；
弹性分界图像获取模块，用于利用所述剪切波加速度数据生成剪切波加速度图像。
另一方面，提供了一种剪切波成像设备，包括：
发射器，用于发射超声波；
采集器，用于采集超声回波数据；
存储器，用于存储计算机程序；
处理器，用于执行所述计算机程序时实现如上述剪切波成像方法的步骤。
另一方面，提供了一种可读存储介质，所述可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现上述剪切波成像方法的步骤。
应用本发明实施例所提供的方法，获取剪切波速度数据；对剪切波速度数据进行一阶导处理，得到剪切波加速度数据；利用剪切波加速度数据生成剪切波加速度图像。
剪切波在传播过程中当遇到软硬度不一的区域时，剪切波会发生反射、折射等现象，使得剪切波速度发生变化，进一步造成剪切波速度与组织软硬度不匹配，从而导致了剪切波弹性图像在分界处得到错误的结果。为了使得剪切波弹性图像在分界处得到正确的结果，即需要得知边界的准确位置。由于剪切波速度在分界处会发生变化，也就是说，找出剪切波速度的加速度，即可确定出分界的准确位置，而速度的一阶导即为加速度。基于此，在本方法中，在获得剪切波速度数据之后，对剪切波速度数据进行一阶导处理，得到剪切波加速度数据，然后基于剪切波加速度数据便可生成剪切波加速度图像。即基于剪切波加速度图像可确定出软硬不一的区域分界，也就是说，可以得到组织弹性变化边界。
相应地，本发明实施例还提供了与上述剪切波成像方法相对应的剪切波成像装置、设备和可读存储介质，具有上述技术效果，在此不再赘述。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例中一种剪切波成像方法的实施流程图；
图2为本发明实施例中一种剪切波生成原理示意图；
图3为本发明实施例中与图2相对应的各个扫描线上的剪切波速度示意图；
图4为本发明实施例中一种剪切波速度图像；
图5本发明实施例一种剪切波加速度图像；
图6为本发明实施例中一种剪切波成像装置的结构示意图；
图7为本发明实施例中一种剪切波成像设备的结构示意图；
图8为本发明实施例中一种剪切波成像设备的结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案，下面结合附图和具体实施方式对本发明作进一步的详细说明。显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
请参考图1，图1为本发明实施例中一种剪切波成像方法的流程图，该方法包括以下步骤：
S101、获取剪切波速度数据。
在本实施例中可直接从可读存储介质中读取预先存储的剪切波速度数据。也可通过实时发射超声波，并采集超声回波数据的方式获取剪切波速度数据，可参见图4。
具体的，实时获取剪切波速度图像的过程，包括：
步骤一、向目标组织发射超声波，并采集超声回波数据；超声回波数据中包括剪切波传输信息；
步骤二、对超声回波数据进行处理，得到剪切波速度数据。
为便于描述，下面将上述两个步骤结合起来进行说明。
其中，目标组织即为需要进行剪切波弹性成像的软组织。利用超声长脉冲产生的声辐射力激励目标组织，以产生振动，进一步在感兴趣区域内产生剪切波。可利用发射跟踪脉冲序列采集感兴趣区域内包括剪切波传播的超声回波数据。
其中，超声波发射与剪切波检测可以为同一探头，也可以为不同探头。
为便于理解，下面以图2为例，对剪切波产生的过程进行详细说明。在图2中的超声探头5发射聚焦长脉冲1，聚焦于某一深度2。通过动量原理，即声辐射力与组织相互作用，该声辐射力被组织大量吸收，使组织内部局部发生振动而产生剪切形变，从而在感兴趣区域内产生剪切波3。
得到超声回波数据之后，对超声回波数据进行处理，便可得到剪切波速度数据。具体的，可对超声回波速度进行波速合成，得到射频数据(RF数据)或调制数据(IQ数据)。基于波束合成后的数据类型，可利用基于RF数据的互相关或基于IQ数据的相移等算法得到组织运动数据。
优选地，为了避免损失原始信息，缩小数据量，得到相位信息，可优选采用IQ数据生成剪切波速度数据。其中，相位信息可用于得到更加准确的剪切波速度数据。具体的实现过程，包括：
步骤2.1、对超声回波数据进行波速合成，得到调制数据；
步骤2.2、对调制数据进行相移处理，得到组织运动数据；
步骤2.3、利用时间飞跃算法或互相关算法对组织运动数据进行处理，得到剪切波速度数据。
得到调整数据之后，可对调整数据进行相移处理，得到组织运动数据。然后利用时间飞跃算法或互相关算对组织运动数据进行处理，即可得到剪切波速度数据。其中，互相关算法，即利用相邻扫描线处的数据来做处理，以得到时间间隔；而对于时间飞跃法，则是利用所采集到的回波数据进行处理得到时间间隔。在实际应用中可选择互相关算法或时间飞跃法来获得剪切波速度，本实施例对于用于剪切波速度的算法并不做限定。
优选地，为了避免组织运动数据中的干扰数据影响剪切波速度数据的质量。还可对组织运动数据中的干扰进行滤波。具体的，可利用方向滤波器滤除超声回波数据中的干扰数据，得到无干扰的超声回波数据。如此，对无干扰的超声回波数据进行处理，可得到质量更高的剪切波速度数据。其中，干扰数据包括噪声数据。由于剪切波在组织传播的过程中，会遇到软硬度不一的界面，因此会造成剪切波的反射。为了解决该问题，使用方向滤波器，即将指定传播方向的剪切波速度给滤出来。
其中，关于步骤2.3，可具体包括：
步骤2.3.1、利用互相关算法或时间飞跃算法计算出剪切波经过每个目标扫描线的前后扫描线对应的时间间隔；
步骤2.3.2、计算每个目标扫描线的前后扫描线的距离与时间间隔的比值，得到剪切波速度。
步骤2.3.3、将剪切波速度和时间间隔作为剪切波速度数据。
也就是说，剪切波速度数据包括每一个扫描线对应的剪切波速度，以及该扫描线的前后扫描线的时间间隔。
S102、对剪切波速度数据进行一阶导处理，得到剪切波加速度数据。
对速度求一阶导，即可获得加速度。具体的，即基于剪切波速度数据，计算出剪切波在分界处发生折射、反射等变化，从而导致剪切波速度发生变化的加速度。
根据各个采集的剪切波速度以及时间间隔，即可计算出相应采集点的加速度。具体的，加速度的计算过程，包括：
步骤一、计算出剪切波经过每个目标扫描线的前后扫描线对应的速度差值；
步骤二、计算速度差值与时间间隔的比值，得到目标扫描线的加速度。
加速度即为速度变化与时间间隔的比值。在本实施例中，可直接从剪切波速度数据中获得各个采集点的剪切波速度、时间间隔。也可采用自相关算法或时间飞跃法，对剪切波速度进行处理，得到相应时间间隔以及速度差值，进而进行计算出加速度。
举例说明：请参考图3，图3为本发明实施例中与图2相对应的各个扫描线上的剪切波速度示意图。其中SwA为图2中扫描线A处的剪切波速度，SwB为图2中扫描线B处的剪切波速度，SwC为图2中扫描线C处的剪切波速度，SwD为图2中扫描线D处的剪切波速度，SwE为图2中扫描线E处的剪切波速度。利用图2中扫描线A处的组织运动数据和图2中扫描线C处的组织运动数据，利用互相关或时间飞跃法等算法处理，即可求出剪切波在经过扫描线A和扫描线C之间的时间间隔TAC。利用SwA和SwC的剪切波速度，结合下式可计算出在SwB处相应的剪切波加速度：
Figure BDA0002442257270000071
以此类推，即可求出在其它采集点的剪切波加速度，
S103、利用剪切波加速度数据生成剪切波加速度图像。
得到若干个未知的加速度之后，便可得到剪切波加速度图像。具体的，可利用灰度映射或彩色映射，对剪切波加速度数据进行线性映射，得到剪切波加速度图像。
在软硬区域不一的情况下，剪切波发生映射、反射等情况时，速度会发生变化，即加速度不为0。即，当加速度不为0的情况下，即可确定组织的软硬度发生了变化，即可定位分界。因此剪切波加速度图像可呈现出清晰分界，具体可参见图5。
应用本发明实施例所提供的方法，获取剪切波速度数据；对剪切波速度数据进行一阶导处理，得到剪切波加速度数据；利用剪切波加速度数据生成剪切波加速度图像。
剪切波在传播过程中当遇到软硬度不一的区域时，剪切波会发生反射、折射等现象，使得剪切波速度发生变化，进一步造成剪切波速度与组织软硬度不匹配，从而导致了剪切波弹性图像在分界处得到错误的结果。为了使得剪切波弹性图像在分界处得到正确的结果，即需要得知边界的准确位置。由于剪切波速度在分界处会发生较大变化，也就是说，找出剪切波速度的加速度，即可确定出分界的准确位置，而速度的一阶导即为加速度。基于此，在本方法中，在获得剪切波速度数据之后，对剪切波速度数据进行一阶导处理，得到剪切波加速度数据，然后基于剪切波加速度数据便可生成剪切波加速度图像。即基于剪切波加速度图像可确定出软硬不一的区域分界，也就是说，可以得到组织弹性变化边界。
优选地，为了更好的呈现组织的弹性情况，还可在上述实施例一的基础上，得到既有清晰边缘有能够展示各个区域的弹性情况的图像。具体的，对剪切波加速度数据和剪切波速度数据进行融合成像处理，得到剪切波弹性图像。即，在进行成像处理时，可利用灰度映射或彩色映射，对剪切波速度数据和剪切波加速度数据进行线性映射，得到剪切波弹性图像。
当然，也可基于剪切波速度数据生成剪切波速度图像，然后叠加剪切波速度图像和剪切波加速度图像，得到剪切波弹性图像。剪切波速度图像即一个关于剪切波速度的信息矩阵，剪切波加速度图像即一个关于剪切波加速度的信息矩阵。在得到感兴趣区域内的剪切波速度和剪切波加速度的信息矩阵时，可利用线性映射的方法，将这些信息矩阵以灰度映射、伪彩或彩色算法将其显示出来。
相应于上面的方法实施例，本发明实施例还提供了一种剪切波成像装置，下文描述的剪切波成像装置与上文描述的剪切波成像方法可相互对应参照。
参见图6所示，该装置包括以下模块：
剪切波速度获取模块101，用于获取剪切波速度数据；
剪切波加速度获取模块102，用于对剪切波速度数据进行一阶导处理，得到剪切波加速度数据；
弹性分界图像获取模块103，用于利用剪切波加速度数据生成剪切波加速度图像。
应用本发明实施例所提供的装置，获取剪切波速度数据；对剪切波速度数据进行一阶导处理，得到剪切波加速度数据；利用剪切波加速度数据生成剪切波加速度图像。
剪切波在传播过程中当遇到软硬度不一的区域时，剪切波会发生反射、折射等现象，使得剪切波速度发生变化，进一步造成剪切波速度与组织软硬度不匹配，从而导致了剪切波弹性图像在分界处得到错误的结果。为了使得剪切波弹性图像在分界处得到正确的结果，即需要得知边界的准确位置。由于剪切波速度在分界处会发生较大变化，也就是说，找出剪切波速度的加速度，即可确定出分界的准确位置，而速度的一阶导即为加速度。基于此，在本装置中，在获得剪切波速度数据之后，对剪切波速度数据进行一阶导处理，得到剪切波加速度数据，然后基于剪切波加速度数据便可生成剪切波加速度图像。即基于剪切波加速度图像可确定出软硬不一的区域分界，也就是说，可以得到组织弹性变化边界。
在本发明的一种具体实施方式中，还包括：
融合成像模块，用于对剪切波加速度数据和剪切波速度数据进行融合成像处理，得到剪切波弹性图像。
在本发明的一种具体实施方式中，剪切波速度获取模块，具体用于向目标组织发射超声波，并采集超声回波数据；超声回波数据中包括剪切波传输信息；对超声回波数据进行处理，得到剪切波速度数据。
在本发明的一种具体实施方式中，剪切波速度获取模块，在对超声回波数据进行处理之前，还可利用方向滤波器滤除干扰数据，得到无干扰的超声回波数据。
在本发明的一种具体实施方式中，剪切波速度获取模块，具体用于利用互相关算法或时间飞跃算法计算出剪切波经过每个目标扫描线的前后扫描线对应的时间间隔；计算每个目标扫描线的前后扫描线的距离与时间间隔的比值，得到剪切波速度；将剪切波速度和时间间隔作为剪切波速度数据。
在本发明的一种具体实施方式中，剪切波加速度获取模块，具体用于计算出剪切波经过每个目标扫描线的前后扫描线对应的速度差值；计算速度差值与时间间隔的比值，得到目标扫描线的加速度。
在本发明的一种具体实施方式中，弹性分界图像获取模块，具体用于利用灰度映射或彩色映射，对剪切波加速度数据进行线性映射，得到剪切波加速度图像。
相应于上面的方法实施例，本发明实施例还提供了一种剪切波成像设备，下文描述的一种剪切波成像设备与上文描述的一种剪切波成像方法可相互对应参照。
参见图7所示，该剪切波成像设备包括：
发射器310，用于发射超声波；
采集器320，用于采集超声回波数据；
存储器332，用于存储计算机程序；
处理器322，用于执行计算机程序时实现上述方法实施例的剪切波成像方法的步骤。需要说明的是，310和320可以是同一个实体。
具体的，请参考图8，为本实施例提供的一种剪切波成像设备的具体结构示意图，该剪切波成像设备可因配置或性能不同而产生比较大的差异，可以包括一个或一个以上的发射器310，一个或一个以上的采集器320，一个或一个以上处理器(central processingunits，CPU)322(例如，一个或一个以上处理器)和存储器332。其中，存储器332可以是短暂存储或持久存储。更进一步地，中央处理器322可以设置为与存储器332通信，在剪切波成像设备301上执行存储器332中的一系列指令操作。
剪切波成像设备301还可以包括一个或一个以上电源326，一个或一个以上有线或无线网络接口350，一个或一个以上输入输出接口358，和/或，一个或一个以上操作系统341。例如，Windows ServerTM，Mac OS XTM，UnixTM，LinuxTM，FreeBSDTM等。
上文所描述的剪切波成像方法中的步骤可以由剪切波成像设备的结构实现。
相应于上面的方法实施例，本发明实施例还提供了一种可读存储介质，下文描述的一种可读存储介质与上文描述的一种剪切波成像方法可相互对应参照。
一种可读存储介质，可读存储介质上存储有计算机程序，计算机程序被处理器执行时实现上述方法实施例的剪切波成像方法的步骤。
该可读存储介质具体可以为U盘、移动硬盘、只读存储器(Read-Only Memory，ROM)、随机存取存储器(Random Access Memory，RAM)、磁碟或者光盘等各种可存储程序代码的可读存储介质。
专业人员还可以进一步意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、计算机软件或者二者的结合来实现，为了清楚地说明硬件和软件的可互换性，在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。

Claims (10)

1.一种剪切波成像方法，其特征在于，包括：
获取剪切波速度数据；
对所述剪切波速度数据进行一阶导处理，得到剪切波加速度数据；
利用所述剪切波加速度数据生成剪切波加速度图像。
2.根据权利要求1所述的剪切波成像方法，其特征在于，还包括：
对所述剪切波加速度数据和剪切波速度数据进行融合成像处理，得到剪切波弹性图像。
3.根据权利要求1所述的剪切波成像方法，其特征在于，所述获取剪切波速度数据，包括：
向目标组织发射超声波，并采集超声回波数据；所述超声回波数据中包括剪切波传输信息；
对所述超声回波数据进行处理，得到所述剪切波速度数据。
4.根据权利要求3所述的剪切波成像方法，其特征在于，在对所述超声回波数据进行处理之前，还包括：
利用方向滤波器滤除干扰数据，得到无干扰的超声回波数据。
5.根据权利要求4所述的剪切波成像方法，其特征在于，对所述超声回波数据进行处理，得到所述剪切波速度数据，包括：
利用互相关算法或时间飞跃算法计算出剪切波经过每个目标扫描线的前后扫描线对应的所述时间间隔；
计算每个目标扫描线的前后扫描线的距离与所述时间间隔的比值，得到剪切波速度；
将所述剪切波速度和所述时间间隔作为所述剪切波速度数据。
6.根据权利要求5所述的剪切波成像方法，其特征在于，对所述剪切波速度数据进行一阶导处理，得到剪切波加速度数据，包括：
计算出剪切波经过每个所述目标扫描线的前后扫描线对应的速度差值；
计算所述速度差值与所述时间间隔的比值，得到所述目标扫描线的所述加速度。
7.根据权利要求1所述的剪切波成像方法，其特征在于，利用所述剪切波加速度数据生成剪切波加速度图像，包括：
利用灰度映射或彩色映射，对所述剪切波加速度数据进行线性映射，得到所述剪切波加速度图像。
8.一种剪切波成像装置，其特征在于，包括：
剪切波速度获取模块，用于获取剪切波速度数据；
剪切波加速度获取模块，用于对所述剪切波速度数据进行一阶导处理，得到剪切波加速度数据；
弹性分界图像获取模块，用于利用所述剪切波加速度数据生成剪切波加速度图像。
9.一种剪切波成像设备，其特征在于，包括：
发射器，用于发射超声波；
采集器，用于采集超声回波数据；
存储器，用于存储计算机程序；
处理器，用于执行所述计算机程序时实现如权利要求1至7任一项所述剪切波成像方法的步骤。
10.一种可读存储介质，其特征在于，所述可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现如权利要求1至7任一项所述剪切波成像方法的步骤。
CN202010268506.6A 2020-04-08 2020-04-08 一种剪切波成像方法、装置、设备及可读存储介质 Active CN111449681B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010268506.6A CN111449681B (zh) 2020-04-08 2020-04-08 一种剪切波成像方法、装置、设备及可读存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010268506.6A CN111449681B (zh) 2020-04-08 2020-04-08 一种剪切波成像方法、装置、设备及可读存储介质

Publications (2)

Publication Number Publication Date
CN111449681A true CN111449681A (zh) 2020-07-28
CN111449681B CN111449681B (zh) 2023-09-08

Family

ID=71671650

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010268506.6A Active CN111449681B (zh) 2020-04-08 2020-04-08 一种剪切波成像方法、装置、设备及可读存储介质

Country Status (1)

Country Link
CN (1) CN111449681B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113076577A (zh) * 2021-03-25 2021-07-06 西南交通大学 高铁路基剪切波速计算方法、装置、设备及可读存储介质

Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102667672A (zh) * 2009-07-07 2012-09-12 韩铮 一种加速度动作识别系统及其方法
CN102754097A (zh) * 2009-12-11 2012-10-24 诺基亚公司 用于呈现内容的第一人世界视图的方法和装置
CN102892358A (zh) * 2010-04-05 2013-01-23 日立阿洛卡医疗株式会社 用于超声成像的方法和装置
US20130046175A1 (en) * 2011-06-16 2013-02-21 Chikayoshi Sumi Imaging method, displacement measurement method and apparatus
CN103098108A (zh) * 2010-11-25 2013-05-08 松下电器产业株式会社 通信设备
US20130261453A1 (en) * 2012-03-30 2013-10-03 Tadashi Tamura Methods and apparatus for ultrasound imaging
CN106618635A (zh) * 2017-01-12 2017-05-10 清华大学 剪切波弹性成像方法和装置
CN107049360A (zh) * 2017-01-26 2017-08-18 清华大学 剪切波弹性成像方法和装置
CN107510474A (zh) * 2017-09-21 2017-12-26 深圳开立生物医疗科技股份有限公司 剪切波弹性成像方法及系统
CN108663265A (zh) * 2018-07-12 2018-10-16 大连理工大学 测量土体剪切波速的组合式扭剪试验装置、系统及方法
CN108733857A (zh) * 2017-04-21 2018-11-02 深圳迈瑞生物医疗电子股份有限公司 超声弹性成像装置及弹性成像结果评价方法
JP2018188887A (ja) * 2017-05-09 2018-11-29 国立大学法人 東京大学 液状化強度比の推定方法
CN109223503A (zh) * 2018-10-09 2019-01-18 深圳市智城华业科技有限公司 测量心肺复苏按压深度的方法、装置、设备和存储介质
CN109596078A (zh) * 2019-01-28 2019-04-09 吉林大学 多信息融合路面谱实时测试系统和测试方法
CN110840488A (zh) * 2019-11-27 2020-02-28 深圳开立生物医疗科技股份有限公司 一种基于剪切波的成像方法、系统及装置

Patent Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102667672A (zh) * 2009-07-07 2012-09-12 韩铮 一种加速度动作识别系统及其方法
CN102754097A (zh) * 2009-12-11 2012-10-24 诺基亚公司 用于呈现内容的第一人世界视图的方法和装置
CN102892358A (zh) * 2010-04-05 2013-01-23 日立阿洛卡医疗株式会社 用于超声成像的方法和装置
CN103098108A (zh) * 2010-11-25 2013-05-08 松下电器产业株式会社 通信设备
US20130046175A1 (en) * 2011-06-16 2013-02-21 Chikayoshi Sumi Imaging method, displacement measurement method and apparatus
US20130261453A1 (en) * 2012-03-30 2013-10-03 Tadashi Tamura Methods and apparatus for ultrasound imaging
CN106618635A (zh) * 2017-01-12 2017-05-10 清华大学 剪切波弹性成像方法和装置
CN107049360A (zh) * 2017-01-26 2017-08-18 清华大学 剪切波弹性成像方法和装置
CN108733857A (zh) * 2017-04-21 2018-11-02 深圳迈瑞生物医疗电子股份有限公司 超声弹性成像装置及弹性成像结果评价方法
JP2018188887A (ja) * 2017-05-09 2018-11-29 国立大学法人 東京大学 液状化強度比の推定方法
CN107510474A (zh) * 2017-09-21 2017-12-26 深圳开立生物医疗科技股份有限公司 剪切波弹性成像方法及系统
CN108663265A (zh) * 2018-07-12 2018-10-16 大连理工大学 测量土体剪切波速的组合式扭剪试验装置、系统及方法
CN109223503A (zh) * 2018-10-09 2019-01-18 深圳市智城华业科技有限公司 测量心肺复苏按压深度的方法、装置、设备和存储介质
CN109596078A (zh) * 2019-01-28 2019-04-09 吉林大学 多信息融合路面谱实时测试系统和测试方法
CN110840488A (zh) * 2019-11-27 2020-02-28 深圳开立生物医疗科技股份有限公司 一种基于剪切波的成像方法、系统及装置

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
YOSHIKI YAMAKOSHI等: "Novel imaging method of continuous shear wave by ultrasound color flow imaging" *
丁玉琴: "场地非线性地震反应分析方法及其应用研究" *
何琼: "超高速超声成像的研究进展", 《中国医学影像技术》 *
傅磊等: "龙门山及其相邻区域的剪切波高频衰减特性", 《世界地震译丛》 *
冯乃章等: "二维Savitzky-Golay数字差分器剪切波弹性成像方法", 《深圳大学学报(理工版)》 *
刘荣桂等: "糖尿病肾病患者肾实质剪切波速度与肾脏血流动力学改变相关性研究", 《中国超声医学工程学会第十二届全国腹部超声医学学术大会》 *
展义臻等: "环氧树脂型相变微胶囊的制备及其在海藻纤维中的应用" *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113076577A (zh) * 2021-03-25 2021-07-06 西南交通大学 高铁路基剪切波速计算方法、装置、设备及可读存储介质

Also Published As

Publication number Publication date
CN111449681B (zh) 2023-09-08

Similar Documents

Publication Publication Date Title
CN107510474B (zh) 剪切波弹性成像方法及系统
CN102641137B (zh) 使用幅度-相位调制超声波的粘弹性测量
EP2356941B1 (en) Ultrasonographic system
US6689063B1 (en) Method and apparatus for acquiring images by recursive ultrasound images
CN111819467B (zh) 波传播和散射参数的估计方法和仪器
US6193665B1 (en) Doppler angle unfolding in ultrasound color flow and Doppler
CN103239258B (zh) 采用超声波的同轴切变波表征
JP2008188424A (ja) 超音波撮像スキャナ向けのガイド画像を形成するための方法及び装置
CN110974296B (zh) 一种获取剪切波弹性图像的方法、装置及超声成像系统
US9883851B2 (en) System and method for shear wave generation with steered ultrasound push beams
CN106419961A (zh) 声学辐射力成像中的自适应运动估计
US11719613B2 (en) Method and device for quantifying viscoelasticity of a medium
CN105530870A (zh) 一种超声成像方法和系统
CN106580372A (zh) 一种超声彩色血流成像的脉冲重复频率调整方法及装置
EP1504284B1 (en) Ultrasound imaging system with a high lateral resolution
JP2000197636A (ja) 複数のアルゴリズムを使用する超音波カラ―・フロ―画像
CN111449681B (zh) 一种剪切波成像方法、装置、设备及可读存储介质
JP4369427B2 (ja) ドプラ速度検出装置及びそれを用いた超音波診断装置
CN111388012B (zh) 用于检测组织硬度的方法、设备及系统
CN111789632B (zh) 超声波诊断装置、信号处理装置以及存储介质
Marion et al. Toward a real-time simulation of ultrasound image sequences based on a 3-D set of moving scatterers
CN116256433A (zh) 用于处理波束成形数据的方法和系统
CN108784744A (zh) 用于剪切波成像的可变焦点
CN114601496A (zh) 一种基于线性阵列的三维超声剪切波弹性成像方法
US5048528A (en) Alias suppression in pulsed doppler systems

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