CN115414065A - 一种剪切波图像生成方法、超声设备及装置 - Google Patents

一种剪切波图像生成方法、超声设备及装置 Download PDF

Info

Publication number
CN115414065A
CN115414065A CN202211048860.3A CN202211048860A CN115414065A CN 115414065 A CN115414065 A CN 115414065A CN 202211048860 A CN202211048860 A CN 202211048860A CN 115414065 A CN115414065 A CN 115414065A
Authority
CN
China
Prior art keywords
particle
shear wave
vibration
shear
mass point
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.)
Pending
Application number
CN202211048860.3A
Other languages
English (en)
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.)
Qingdao Hisense Medical Equipment Co Ltd
Original Assignee
Qingdao Hisense Medical Equipment 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 Qingdao Hisense Medical Equipment Co Ltd filed Critical Qingdao Hisense Medical Equipment Co Ltd
Priority to CN202211048860.3A priority Critical patent/CN115414065A/zh
Priority to PCT/CN2022/116912 priority patent/WO2023071528A1/zh
Publication of CN115414065A publication Critical patent/CN115414065A/zh
Pending legal-status Critical Current

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/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4411—Device being modular
    • 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/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10132—Ultrasound image
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20092—Interactive image processing based on input by user
    • G06T2207/20104—Interactive definition of region of interest [ROI]

Landscapes

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

Abstract

本申请公开了一种剪切波图像生成方法、超声设备及装置。本申请利用声辐射力或低频振动在目标部位上进行激励产生剪切波，在剪切波传播路径上对剪切波进行检测，获取剪切波的检测数据；确定剪切波的各个质点在不同时刻分别对应的质点振动速度；针对任一个质点在不同时刻对应的质点振动速度，对质点在不同时刻对应的质点振动速度进行多个方向的滤波处理，得到每个质点振动速度在各个方向对应的参考质点振动速度，得到任一个方向上每个质点对应的剪切波速度；并根据预设的超声黑白图像以及各个方向上每个质点对应的剪切波速度，生成剪切波图像。由此，能够提高剪切波传播速度的方向性和准确性表征，提供一个更全面的剪切波信息。

Description

一种剪切波图像生成方法、超声设备及装置
技术领域
本申请涉及剪切波弹性成像技术领域，特别涉及一种剪切波图像生成方法、超声设备及装置。
背景技术
各向异性是一种方向依赖性的特性,存在于富含纤维的生物组织中。当骨骼肌、心肌、肌腱等生物软组织受到外力时，生物软组织的力学特性表现出一定的各向异性。在超声弹性成像系统中，用一种带方向性剪切波的传播速度的矢量场信息来描述当前生物软组织的力学特性。
但是，目前剪切波成像时一般会首先假设剪切波在各向同性的生物软组织中传播，然后计算剪切波的两个质点之间的剪切波速度，从而预估出生物软组织的力学特性。这种方法得到的是单一的剪切波速度信息，使得表征的剪切波传播速度信息的准确率较低。
发明内容
本申请的目的是提供一种剪切波图像生成方法、超声设备及装置，用以解决现有技术中表征的剪切波传播速度信息的准确率较低的问题。
第一方面，本申请提供一种剪切波图像生成方法，所述方法包括：
利用声辐射力或低频振动在对象的目标部位上进行激励产生剪切波，在剪切波传播路径上对所述剪切波进行检测，获取所述剪切波的检测数据；
基于所述检测数据，确定所述剪切波的各个质点在不同时刻分别对应的质点振动速度；
针对任一个质点在不同时刻对应的质点振动速度，对所述质点在不同时刻对应的质点振动速度进行多个方向的滤波处理，得到每个质点振动速度在各个方向对应的参考质点振动速度；
针对任一个方向，根据所述方向上各个质点的参考质点振动速度，确定所述方向上每个质点对应的剪切波速度；
根据预设的超声黑白图像，以及各个方向上每个质点对应的剪切波速度，生成剪切波图像。
在一种可能的实施方式中，所述基于所述检测数据，确定所述剪切波的各个质点在不同时刻分别对应的质点振动速度，包括：
针对任意一个质点，基于所述检测数据中所述质点对应的数据，利用二维自相关算法确定所述质点在每个时刻对应的预设时长内的轴向位移；
基于每个时刻对应的预设时长内的轴向位移以及所述预设时长，确定所述质点在不同时刻分别对应的质点振动速度。
在一种可能的实施方式中，所述对所述质点在不同时刻对应的质点振动速度进行多个不同方向的滤波处理，得到每个质点振动速度在各个方向对应的参考质点振动速度，包括：
对所述质点在不同时刻对应的质点振动速度分别进行傅里叶变换，得到每个质点振动速度对应的频域波数域数据；
针对任一个频域波数域数据，将所述频域波数域数据与所述多个方向对应的预设掩膜分别进行相乘；并将各个乘积均进行傅里叶逆变换，得到所述频域波数域数据在各个方向对应的参考质点振动速度。
在一种可能的实施方式中，所述预设掩膜为边缘采用截取函数进行截断之后的掩膜。
在一种可能的实施方式中，所述根据所述方向上各个质点的参考质点振动速度，确定所述方向上每个质点对应的剪切波速度，包括：
针对所述方向上的任一个质点分别执行以下操作：
基于所述质点的位置信息以及候选质点的位置信息，确定所述质点与所述候选质点之间的距离；所述候选质点为在所述方向上与所述质点相邻的一个或多个连续的质点；
基于所述质点与所述候选质点在所述方向上的参考质点振动速度，以及所述质点与所述候选质点之间的距离，确定所述质点对应的剪切波速度。
在一种可能的实施方式中，所述基于所述质点与所述候选质点在所述方向上的参考质点振动速度，以及所述质点与所述候选质点之间的距离，确定所述质点对应的剪切波速度，包括：
基于所述质点在所述方向上的参考质点振动速度的峰值点对应的第一时间，以及所述候选质点在所述方向上的参考质点振动速度的峰值点对应的第二时间，确定所述第一时间和所述第二时间之间的时间差值；
基于所述质点与所述候选质点之间的距离，以及所述第一时间和所述第二时间之间的时间差值，得到所述质点对应的剪切波速度。
在一种可能的实施方式中，所述利用声辐射力或低频振动在对象的目标部位上进行激励产生剪切波，在剪切波传播路径上对所述剪切波进行检测，获取所述剪切波的检测数据，包括：
利用声辐射力或低频振动在对象的目标部位上进行激励产生剪切波，并在剪切波传播路径上多次发射超声信号对所述剪切波进行检测，获取回波数据；
对所述回波数据进行正交解调处理和降频处理，得到所述剪切波的检测数据。
第二方面，本申请实施例提供一种超声设备，所述超声设备包括至少一个处理器、以及至少一个存储器；其中，所述存储器存储有程序代码，当程序代码被所述处理器执行时，使得所述处理器执行下列过程：
利用声辐射力或低频振动在对象的目标部位上进行激励产生剪切波，在剪切波传播路径上对所述剪切波进行检测，获取所述剪切波的检测数据；
基于所述检测数据，确定所述剪切波的各个质点在不同时刻分别对应的质点振动速度；
针对任一个质点在不同时刻对应的质点振动速度，对所述质点在不同时刻对应的质点振动速度进行多个方向的滤波处理，得到每个质点振动速度在各个方向对应的参考质点振动速度；
针对任一个方向，根据所述方向上各个质点的参考质点振动速度，确定所述方向上每个质点对应的剪切波速度；
根据预设的超声黑白图像，以及各个方向上每个质点对应的剪切波速度，生成剪切波图像。
在一种可能的实施方式中，在执行所述对所述质点在不同时刻对应的质点振动速度进行多个不同方向的滤波处理，得到每个质点振动速度在各个方向对应的参考质点振动速度时，所述处理器具体用于：
对所述质点在不同时刻对应的质点振动速度分别进行傅里叶变换，得到每个质点振动速度对应的频域波数域数据；
针对任一个频域波数域数据，将所述频域波数域数据与所述多个方向对应的预设掩膜分别进行相乘；并将各个乘积均进行傅里叶逆变换，得到所述频域波数域数据在各个方向对应的参考质点振动速度。
第三方面，本申请实施例提供一种剪切波图像生成装置，所述装置包括：
获取模块，用于利用声辐射力或低频振动在对象的目标部位上进行激励产生剪切波，在剪切波传播路径上对所述剪切波进行检测，获取所述剪切波的检测数据；
质点振动速度确定模块，用于基于所述检测数据，确定所述剪切波的各个质点在不同时刻分别对应的质点振动速度；
参考质点振动速度确定模块，用于针对任一个质点在不同时刻对应的质点振动速度，对所述质点在不同时刻对应的质点振动速度进行多个方向的滤波处理，得到每个质点振动速度在各个方向对应的参考质点振动速度；
剪切波速度确定模块，用于针对任一个方向，根据所述方向上各个质点的参考质点振动速度，确定所述方向上每个质点对应的剪切波速度；
图像生成模块，用于根据预设的超声黑白图像，以及各个方向上每个质点对应的剪切波速度，生成剪切波图像。
第四方面，本申请提供一种计算机可读存储介质，当所述计算机可读存储介质中的指令由超声设备执行时，使得所述超声设备能够执行如上述第一方面中任一项所述的剪切波图像生成方法。
第五方面，本申请提供一种计算机程序产品，包括计算机程序：所述计算机程序被处理器执行时实现如上述第一方面中任一项所述的剪切波图像生成方法。
本申请的实施例提供的技术方案至少带来以下有益效果：
本申请实施例利用声辐射力或低频振动在对象的目标部位上进行激励产生剪切波，在剪切波传播路径上对剪切波进行检测，获取剪切波的检测数据；基于检测数据，确定剪切波的各个质点在不同时刻分别对应的质点振动速度；针对任一个质点在不同时刻对应的质点振动速度，对质点在不同时刻对应的质点振动速度进行多个方向的滤波处理，得到每个质点振动速度在各个方向对应的参考质点振动速度；针对任一个方向，根据方向上各个质点的参考质点振动速度，确定方向上每个质点对应的剪切波速度；根据预设的超声黑白图像，以及各个方向上每个质点对应的剪切波速度，生成剪切波图像。
由此，通过将剪切波传播速度的不同方向的信息展现并叠加生成最终的剪切波图像，提高了剪切波传播速度的准确性，解决了单一维度(只有大小信息)的剪切波传播速度的准确率较低的问题，提供了一个更全面的剪切波信息，可以极大的提高临床医生的工作效率，对疾病诊断具有重要的临床意义。
本申请的其它特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本申请而了解。本申请的目的和其他优点可通过在所写的说明书、权利要求书、以及附图中所特别指出的结构来实现和获得。
附图说明
为了更清楚地说明本申请实施例的技术方案，下面将对本申请实施例中所需要使用的附图作简单地介绍，显而易见地，下面所介绍的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本申请实施例一种超声设备的结构示意图；
图2为本申请实施例一种超声设备的应用原理的示意图；
图3为本申请实施例提供的一种剪切波图像生成方法的应用场景图；
图4为本申请实施例提供的一种剪切波图像生成方法的流程示意图；
图5为本申请实施例提供的多个方向对应的预设掩膜的示意图；
图6为本申请实施例提供的多个方向对应的参考质点振动速度的示意图；
图7为本申请实施例提供的一种每个质点对应的剪切波速度的确定方法的流程示意图；
图8为本申请实施例提供的一种每个质点对应的剪切波速度的确定方法的流程示意图；
图9为本申请实施例提供的质点与候选质点的示意图；
图10为本申请实施例提供的质点与候选质点在不同时刻分别对应的参考质点振动速度的示意图；
图11为本申请实施例提供的一个方向上的剪切波图像的示意图；
图12为本申请实施例提供的另一个方向上的剪切波图像的示意图；
图13为本申请实施例提供的另一个方向上的剪切波图像的示意图；
图14为本申请实施例提供的另一个方向上的剪切波图像的示意图；
图15为本申请实施例提供的另一个方向上的剪切波图像的示意图；
图16为本申请实施例提供的另一个方向上的剪切波图像的示意图；
图17为本申请实施例提供的另一个方向上的剪切波图像的示意图；
图18为本申请实施例提供的另一个方向上的剪切波图像的示意图；
图19为本申请实施例提供的包含有各个方向的剪切波速度的剪切波图像；
图20为本申请实施例提供的一种剪切波图像生成方法的完整流程示意图；
图21为本申请实施例提供的一种剪切波图像生成装置的结构示意图。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述。其中，所描述的实施例是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本申请保护的范围。
并且，在本申请实施例的描述中，除非另有说明，“/”表示或的意思，例如，A/B可以表示A或B；文本中的“和/或”仅仅是一种描述关联对象的关联关系，表示可以存在三种关系，例如，A和/或B，可以表示：单独存在A，同时存在A和B，单独存在B这三种情况，另外，在本申请实施例的描述中，“多个”是指两个或多于两个。
以下，术语“第一”、“第二”仅用于描述目的，而不能理解为暗示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”、的特征可以明示或者隐含地包括一个或者更多个该特征，在本申请实施例的描述中，除非另有说明，“多个”的含义是两个或两个以上。
本申请实施例利用声辐射力或低频振动在对象的目标部位上进行激励产生剪切波，在剪切波传播路径上对剪切波进行检测，获取剪切波的检测数据；基于检测数据，确定剪切波的各个质点在不同时刻分别对应的质点振动速度；针对任一个质点在不同时刻对应的质点振动速度，对质点在不同时刻对应的质点振动速度进行多个方向的滤波处理，得到每个质点振动速度在各个方向对应的参考质点振动速度；针对任一个方向，根据方向上各个质点的参考质点振动速度，确定方向上每个质点对应的剪切波速度；根据预设的超声黑白图像，以及各个方向上每个质点对应的剪切波速度，生成剪切波图像。由此，通过将剪切波传播速度的不同方向的信息展现并叠加生成最终的剪切波图像，提高了剪切波传播速度的准确性，解决了单一维度(只有大小信息)的剪切波传播速度的准确率较低的问题，提供了一个更全面的剪切波信息，可以极大的提高临床医生的工作效率，对疾病诊断具有重要的临床意义。
在介绍完本申请实施例的主要发明思想之后，下面结合附图对本申请实施例提供的一种剪切波图像生成方法应用的超声设备进行介绍。如图1所示，为本申请实施例提供的一种超声设备的结构框图。
应该理解的是，图1所示超声设备仅是一个范例，并且超声设备可以具有比图1中所示的更多的或者更少的部件，可以组合两个或多个的部件，或者可以具有不同的部件配置。图中所示出的各种部件可以在包括一个或多个信号处理和/或专用集成电路在内的硬件、软件、或硬件和软件的组合中实现。
图1中示例性示出了超声设备的硬件配置框图。
如图1所示，超声设备例如可以包括：处理器110、存储器120、显示单元130和探头140；其中：
探头140，用于发射超声信号；
显示单元130，用于显示剪切波图像；
存储器120被配置为存储用于剪切波图像生成所需的数据，可包括软件程序，应用界面数据等；
处理器110，分别与探头140、显示单元130和存储器120相连接，被配置为执行：利用声辐射力或低频振动在对象的目标部位上进行激励产生剪切波，在剪切波传播路径上对剪切波进行检测，获取剪切波的检测数据；基于检测数据，确定剪切波的各个质点在不同时刻分别对应的质点振动速度；针对任一个质点在不同时刻对应的质点振动速度，对质点在不同时刻对应的质点振动速度进行多个方向的滤波处理，得到每个质点振动速度在各个方向对应的参考质点振动速度；针对任一个方向，根据方向上各个质点的参考质点振动速度，确定方向上每个质点对应的剪切波速度；根据预设的超声黑白图像，以及各个方向上每个质点对应的剪切波速度，生成剪切波图像。
通过显示单元130展示根据各个方向上每个质点对应的剪切波速度生成的各个方向的剪切波图像以及根据预设的超声黑白图像，以及各个方向上每个质点对应的剪切波速度，生成的最终的剪切波图像。
图2为根据本申请一个实施例的应用原理的示意图。其中，该部分可由图1所示超声设备的部分模块或功能组件实现，下面将仅针对主要的部件进行说明，而其它部件，如存储器、控制器、控制电路等，此处将不进行赘述。
如图2所示，应用环境中可以包括用户界面210、用于显示用户界面的显示单元220以及处理器230。
显示单元220可以包括显示面板221、背光组件222。其中，显示面板221被配置为对超声图像进行显示，背光组件222位于显示面板221背面，背光组件222可以包括多个背光分区(图中未示出)，各背光分区可以发光，以点亮显示面板221。
处理器230可以被配置为控制背光组件222中各背光分区的背光源亮度，以及控制探头发射超声信号，并接收超声回波信号并进行分析得到超声图像。
其中，处理器230可以对剪切波的回波信号进行处理，得到剪切波图像。
如图3所示，为本申请实施例提供的一种剪切波图像生成方法的应用场景图。图中包括：超声设备10、患者20、存储器30；
其中，超声设备10用于利用声辐射力或低频振动在患者20的目标部位上进行激励产生剪切波，在剪切波传播路径上对所述剪切波进行检测，获取所述剪切波的检测数据，并存储在存储器30中；
然后超声设备10利用存储在存储器30中的检测数据，确定所述剪切波的各个质点在不同时刻分别对应的质点振动速度；针对任一个质点在不同时刻对应的质点振动速度，对所述质点在不同时刻对应的质点振动速度进行多个方向的滤波处理，得到每个质点振动速度在各个方向对应的参考质点振动速度；针对任一个方向，根据所述方向上各个质点的参考质点振动速度，确定所述方向上每个质点对应的剪切波速度；根据预设的超声黑白图像，以及各个方向上每个质点对应的剪切波速度，生成剪切波图像，并进行显示。
当然，本申请实施例提供的方法并不限于图3所示的应用场景，还可以用于其它可能的应用场景，本申请实施例并不进行限制。对于图3所示的应用场景的各个设备所能实现的功能将在后续的方法实施例中一并进行描述，在此先不过多赘述。
为进一步说明本申请实施例提供的技术方案，下面结合附图以及具体实施方式对此进行详细的说明。虽然本申请实施例提供了如下述实施例或附图所示的方法操作步骤，但基于常规或者无需创造性的劳动在方法中可以包括更多或者更少的操作步骤。在逻辑上不存在必要因果关系的步骤中，这些步骤的执行顺序不限于本申请实施例提供的执行顺序。
参见图4，为本申请实施例提供的一种剪切波图像生成方法的流程示意图。如图4所示，该方法包括以下步骤：
在步骤401中，利用声辐射力或低频振动在对象的目标部位上进行激励产生剪切波，在剪切波传播路径上对剪切波进行检测，获取剪切波的检测数据。
在一种可能的实施方式中，利用声辐射力或低频振动在对象的目标部位上进行激励产生剪切波，在剪切波传播路径上对剪切波进行检测，获取剪切波的检测数据，可以执行为：利用声辐射力或低频振动在对象的目标部位上进行激励产生剪切波，并在剪切波传播路径上多次发射超声信号对剪切波进行检测，获取回波数据；对回波数据进行正交解调处理和降频处理，得到剪切波的检测数据。
其中，对回波数据进行正交解调处理和降频处理后会得到IQ数据，IQ数据即为剪切波的检测数据。
在步骤402中，基于检测数据，确定剪切波的各个质点在不同时刻分别对应的质点振动速度。
在一种可能的实施方式中，基于检测数据，确定剪切波的各个质点在不同时刻分别对应的质点振动速度，可以执行为：针对任意一个质点，基于检测数据中质点对应的数据，利用二维自相关算法确定质点在每个时刻对应的预设时长内的轴向位移；基于每个时刻对应的预设时长内的轴向位移以及预设时长，确定质点在不同时刻分别对应的质点振动速度。
具体可以实施为，首先在对回波数据进行正交解调处理和降频处理后会得到IQ数据，以及水平方向上的每一条剪切波检测线上的各个质点；然后针对任意一条剪切波检测线上的各个质点中的任意一个质点，在质点的预设轴向范围内，使用提前设置的垂直方向上的多个样点和时间方向上的多个样点，根据下列公式计算得到质点在每个时刻对应的预设时长内的轴向位移：
Figure BDA0003822952410000111
其中，
Figure BDA0003822952410000112
是质点在一个时刻对应的预设时长内的轴向位移；M是垂直方向上的样点个数；N是时间方向上的样点个数；m是垂直方向上的每个样点的坐标；n是时间方向上的每个样点的坐标；c是声音在目标部位的传播速度；π是圆周率；fc是回波数据信号的中心频率；I和Q分别是回波数据信号的同相和正交分量，是对回波数据信号进行正交解调处理和降频处理后得到的。
然后，对每个质点都使用上述公式计算后，会得到每个质点在每个时刻对应的预设时长内的轴向位移，然后基于每个时刻对应的预设时长内的轴向位移以及预设时长，根据下列公式就可以确定质点在该时刻对应的质点振动速度：
Figure BDA0003822952410000113
其中，
Figure BDA0003822952410000114
表示质点在该时刻对应的质点振动速度，
Figure BDA0003822952410000115
表示质点在该时刻对应的预设时长内的轴向位移，△t表示该时刻对应的预设时长。上述公式表示质点在一个时刻对应的预设时长内的轴向位移除以预设时长，就可以得到质点在该时刻对应的质点振动速度，也表示对质点在一个时刻对应的预设时长内的轴向位移进行微分，就可以得到质点在该时刻对应的质点振动速度。
在步骤403中，针对任一个质点在不同时刻对应的质点振动速度，对质点在不同时刻对应的质点振动速度进行多个方向的滤波处理，得到每个质点振动速度在各个方向对应的参考质点振动速度。
在一种可能的实施方式中，对质点在不同时刻对应的质点振动速度进行多个不同方向的滤波处理，得到每个质点振动速度在各个方向对应的参考质点振动速度，可以执行为：
对质点在不同时刻对应的质点振动速度分别进行傅里叶变换，得到每个质点振动速度对应的频域波数域数据；针对任一个频域波数域数据，将频域波数域数据与多个方向对应的预设掩膜分别进行相乘；并将乘积进行傅里叶逆变换，得到频域波数域数据在各个方向对应的参考质点振动速度。
在一种可能的实施方式中，为了防止出现吉布斯效应，本申请实施例中的预设掩膜为边缘采用截取函数进行截断之后的掩膜。
具体实施时，将每个质点在不同时刻对应的质点振动速度分别进行傅里叶变换，将时间空间域数据转换为频域波数域数据，得到每个质点振动速度对应的频域波数域数据，然后根据用户感兴趣的方向，将每个质点的所有频域波数域数据全部与用户感兴趣的方向对应的预设掩膜进行相乘，抑制其他方向传播的剪切波，然后将乘积进行傅里叶逆变换，就可以得到用户感兴趣的方向的每个质点在不同时刻对应的参考质点振动速度。
其中，本申请实施例中根据用户感兴趣的方向，将每个质点的所有频域波数域数据全部与用户感兴趣的方向对应的预设掩膜进行相乘时，可以使用边缘采用截取函数进行截断之后的掩膜作为预设掩膜，也可以使用边缘未采用截取函数进行截断的掩膜作为预设掩膜，本申请实施例对此不作限制，可以根据实际情况选择其中的一个进行使用。
示例性的，如图5所示，为本申请实施例提供的多个方向对应的预设掩膜的示意图，可以针对任一个频域波数域数据，将频域波数域数据分别与图5所示的0°、45°、90°、135°、180°、225°、270°、315°对应的预设掩膜分别进行相乘，再进行傅里叶逆变换，就可以得到如图6所示的0°，45°，90°，135°，180°，225°，270°、315°对应的参考质点振动速度。其中，将频域波数域数据与图5所示的0°对应的预设掩膜分别进行相乘，再进行傅里叶逆变换，得到图6所示的经过滤波处理后的0°对应的参考质点振动速度，其他方向和0°方向得到的参考质点振动速度的过程一样，在此不一一进行赘述。其中多个方向分别设置为0°，45°，90°，135°，180°，225°，270°、315°，也可以根据用户感兴趣的方向设置为其他的方向，本申请实施例对此不做限制。
在一种可能的实施方式中，为了消除孤立的质点，本申请实施例中在对质点在不同时刻对应的质点振动速度进行多个不同方向的滤波处理之前，还可以对任一个质点在不同时刻对应的质点振动速度进行低通滤波和中值滤波。
由此，可以排除剪切波之间的干涉，以便在下面的步骤404中能够更加精确的计算出每个方向的剪切波的相邻质点之间的时延。
在步骤404中，针对任一个方向，根据方向上各个质点的参考质点振动速度，确定方向上每个质点对应的剪切波速度。
在一种可能的实施方式中，根据方向上各个质点的参考质点振动速度，确定方向上每个质点对应的剪切波速度，可以针对方向上的任一个质点分别执行图7所示的步骤：
在步骤701中，基于质点的位置信息以及候选质点的位置信息，确定质点与候选质点之间的距离；候选质点为在方向上与质点相邻的一个或多个连续的质点。
在步骤702中，基于质点与候选质点在方向上的参考质点振动速度，以及质点与候选质点之间的距离，确定质点对应的剪切波速度。
在一种可能的实施方式中，基于质点与候选质点在方向上的参考质点振动速度，以及质点与候选质点之间的距离，确定质点对应的剪切波速度，可以执行为如图8所示的步骤：
在步骤801中，基于质点在方向上的参考质点振动速度的峰值点对应的第一时间，以及候选质点在方向上的参考质点振动速度的峰值点对应的第二时间，确定第一时间和第二时间之间的时间差值；
在步骤802中，基于质点与候选质点之间的距离，以及第一时间和第二时间之间的时间差值，得到质点对应的剪切波速度。
示例性的，如图9所示，计算90°方向上质点A对应的剪切波速度，若候选质点为在90°方向上与质点A相邻的一个质点，则首先在90°方向上选取与质点A相邻的一个质点，即质点B，然后基于质点A的位置信息以及质点B的位置信息，确定质点A与质点B之间的距离d；接着基于如图10所示的质点A与质点B在90°方向上的参考质点振动速度，确定基于质点A在90°方向上的参考质点振动速度的峰值点对应的第一时间tA，以及质点B在90°方向上的参考质点振动速度的峰值点对应的第二时间tB，利用互相关的原理对在90°方向上的质点A计算时延，即确定第一时间和第二时间之间的时间差值tB-tA，最后基于质点A与质点B之间的距离d，以及第一时间和第二时间之间的时间差值tB-tA，根据下列公式得到质点A对应的剪切波速度：
Figure BDA0003822952410000141
示例性的，计算90°方向上质点A对应的剪切波速度，若候选质点为在90°方向上与质点A相邻的多个连续的质点，则在90°方向上选取与质点A相邻的多个连续的质点，如质点B、质点C、质点D，然后分别确定质点A与质点B、质点C、质点D之间的距离d1、d2、d3，同时基于质点A、质点B、质点C、质点D在90°方向上的参考质点振动速度，分别确定质点A与质点B、质点C、质点D之间的时间差值t1、t2、t3，然后使用上述公式确定质点A对应的多个剪切波速度为：
Figure BDA0003822952410000142
最后使用V1、V2、V3的均值作为质点A最终对应的剪切波速度。
由此，可以计算得到带有方向信息的质点A对应的剪切波速度大小，使得剪切波速度具有大小和方向两个特征信息。
在步骤405中，根据预设的超声黑白图像，以及各个方向上每个质点对应的剪切波速度，生成剪切波图像。
具体实施时，在计算得到各个方向上每个质点对应的剪切波速度后，可以将各个方向上每个质点对应的剪切波速度按照矢量信息(方向信息与大小信息)，和预设的超声黑白图像进行叠加之后显示，分别得到图11所示的0°方向上的剪切波图像，图12所示的45°方向上的剪切波图像，图13所示的90°方向上的剪切波图像，图14所示的135°方向上的剪切波图像，图15所示的180°方向上的剪切波图像，图16所示的225°方向上的剪切波图像，图17所示的270°方向上的剪切波图像，图18所示的315°方向上的剪切波图像，最后可以将图11-图18所示的每一个方向对应的剪切波图像按照每个方向对应的权重进行融合，得到如图19所示的一张包含有各个方向的剪切波速度的剪切波图像。然后可以将图11-图19所示的图像全部存储在超声设备中，就可以根据用户的需求显示剪切波图像，例如用户需要查看90°方向上的剪切波图像，就可以选取图13所示的剪切波图像进行显示；例如用户需要查看各个方向上的剪切波图像，则可以选取图19所示的剪切波图像进行显示。
在图11-图19中，每一个箭头的长度代表了每个质点对应的剪切波速度大小，箭头越长，质点对应的剪切波速度越大；每一个箭头的方向代表了每个质点对应的剪切波速度的方向。
在另一种示例中，也可以使用箭头的长度和颜色共同代表每个质点对应的剪切波速度大小，箭头越长，颜色的RGB值越高，则表示质点对应的剪切波速度越大；使用箭头的方向代表了每个质点对应的剪切波速度的方向。
为了便于理解，下面结合图20对本申请实施例提供的一种剪切波图像生成方法的完整流程进行说明。如图20所示，该方法包括以下步骤：
在步骤2001中，利用声辐射力或低频振动在对象的目标部位上进行激励产生剪切波，在剪切波传播路径上对剪切波进行检测，获取剪切波的检测数据。
在步骤2002中，针对任意一个质点，基于检测数据中质点对应的数据，利用二维自相关算法确定质点在每个时刻对应的预设时长内的轴向位移；基于每个时刻对应的预设时长内的轴向位移以及预设时长，确定质点在不同时刻分别对应的质点振动速度。
在步骤2003中，针对任一个质点在不同时刻对应的质点振动速度，对质点在不同时刻对应的质点振动速度分别进行傅里叶变换，得到每个质点振动速度对应的频域波数域数据。
在步骤2004中，针对任一个频域波数域数据，将频域波数域数据与多个方向对应的预设掩膜分别进行相乘；并将乘积进行傅里叶逆变换，得到频域波数域数据在各个方向对应的参考质点振动速度。
在步骤2005中，针对任一个方向上的任一个质点，基于质点的位置信息以及候选质点的位置信息，确定质点与候选质点之间的距离；基于质点在方向上的参考质点振动速度的峰值点对应的第一时间，以及候选质点在方向上的参考质点振动速度的峰值点对应的第二时间，确定第一时间和第二时间之间的时间差值。
在步骤2006中，基于质点与候选质点之间的距离，以及第一时间和第二时间之间的时间差值，得到质点对应的剪切波速度。
在步骤2007中，根据预设的超声黑白图像，以及各个方向上每个质点对应的剪切波速度，生成剪切波图像。
基于前文的描述，本申请实施例利用声辐射力或低频振动在对象的目标部位上进行激励产生剪切波，在剪切波传播路径上对剪切波进行检测，获取剪切波的检测数据；基于检测数据，确定剪切波的各个质点在不同时刻分别对应的质点振动速度；针对任一个质点在不同时刻对应的质点振动速度，对质点在不同时刻对应的质点振动速度进行多个方向的滤波处理，得到每个质点振动速度在各个方向对应的参考质点振动速度；针对任一个方向，根据方向上各个质点的参考质点振动速度，确定方向上每个质点对应的剪切波速度；根据预设的超声黑白图像，以及各个方向上每个质点对应的剪切波速度，生成剪切波图像。
由此，通过将剪切波传播速度的不同方向的信息展现并叠加生成最终的剪切波图像，提高了剪切波传播速度的准确性，解决了单一维度(只有大小信息)的剪切波传播速度的准确率较低的问题，提供了一个更全面的剪切波信息，可以极大的提高临床医生的工作效率，对疾病诊断具有重要的临床意义。
根据相同的发明构思，本申请实施例提供了一种剪切波图像生成装置。如图21所示，为本申请实施例提供的一种剪切波图像生成装置的示意图，包括：
获取模块2101，用于利用声辐射力或低频振动在对象的目标部位上进行激励产生剪切波，在剪切波传播路径上对所述剪切波进行检测，获取所述剪切波的检测数据；
质点振动速度确定模块2102，用于基于所述检测数据，确定所述剪切波的各个质点在不同时刻分别对应的质点振动速度；
参考质点振动速度确定模块2103，用于针对任一个质点在不同时刻对应的质点振动速度，对所述质点在不同时刻对应的质点振动速度进行多个方向的滤波处理，得到每个质点振动速度在各个方向对应的参考质点振动速度；
剪切波速度确定模块2104，用于针对任一个方向，根据所述方向上各个质点的参考质点振动速度，确定所述方向上每个质点对应的剪切波速度；
图像生成模块2105，用于根据预设的超声黑白图像，以及各个方向上每个质点对应的剪切波速度，生成剪切波图像。
在一种可能的实施方式中，在执行所述基于所述检测数据，确定所述剪切波的各个质点在不同时刻分别对应的质点振动速度时，质点振动速度确定模块2102，具体用于：
针对任意一个质点，基于所述检测数据中所述质点对应的数据，利用二维自相关算法确定所述质点在每个时刻对应的预设时长内的轴向位移；
基于每个时刻对应的预设时长内的轴向位移以及所述预设时长，确定所述质点在不同时刻分别对应的质点振动速度。
在一种可能的实施方式中，在执行所述对所述质点在不同时刻对应的质点振动速度进行多个不同方向的滤波处理，得到每个质点振动速度在各个方向对应的参考质点振动速度时，参考质点振动速度确定模块2103，具体用于：
对所述质点在不同时刻对应的质点振动速度分别进行傅里叶变换，得到每个质点振动速度对应的频域波数域数据；
针对任一个频域波数域数据，将所述频域波数域数据与所述多个方向对应的预设掩膜分别进行相乘；并将各个乘积均进行傅里叶逆变换，得到所述频域波数域数据在各个方向对应的参考质点振动速度。
在一种可能的实施方式中，所述预设掩膜为边缘采用截取函数进行截断之后的掩膜。
在一种可能的实施方式中，在执行所述根据所述方向上各个质点的参考质点振动速度，确定所述方向上每个质点对应的剪切波速度时，剪切波速度确定模块2104，具体用于：
针对所述方向上的任一个质点分别执行以下操作：
基于所述质点的位置信息以及候选质点的位置信息，确定所述质点与所述候选质点之间的距离；所述候选质点为在所述方向上与所述质点相邻的一个或多个连续的质点；
基于所述质点与所述候选质点在所述方向上的参考质点振动速度，以及所述质点与所述候选质点之间的距离，确定所述质点对应的剪切波速度。
在一种可能的实施方式中，在执行所述基于所述质点与所述候选质点在所述方向上的参考质点振动速度，以及所述质点与所述候选质点之间的距离，确定所述质点对应的剪切波速度时，剪切波速度确定模块2104，具体用于：
基于所述质点在所述方向上的参考质点振动速度的峰值点对应的第一时间，以及所述候选质点在所述方向上的参考质点振动速度的峰值点对应的第二时间，确定所述第一时间和所述第二时间之间的时间差值；
基于所述质点与所述候选质点之间的距离，以及所述第一时间和所述第二时间之间的时间差值，得到所述质点对应的剪切波速度。
在一种可能的实施方式中，在执行所述利用声辐射力或低频振动在对象的目标部位上进行激励产生剪切波，在剪切波传播路径上对所述剪切波进行检测，获取所述剪切波的检测数据时，获取模块2101，具体用于：
利用声辐射力或低频振动在对象的目标部位上进行激励产生剪切波，并在剪切波传播路径上多次发射超声信号对所述剪切波进行检测，获取回波数据；
对所述回波数据进行正交解调处理和降频处理，得到所述剪切波的检测数据。
在示例性实施例中，还提供了一种包括指令的计算机可读存储介质，例如包括指令的存储器，上述指令可由处理器执行以完成上述剪切波图像生成方法。可选地，存储介质可以是非临时性计算机可读存储介质，例如，所述非临时性计算机可读存储介质可以是ROM、随机存取存储器(RAM)、CD-ROM、磁带、软盘和光数据存储设备等。
在示例性实施例中，还提供一种计算机程序产品，包括计算机程序，所述计算机程序被处理器执行时实现如本申请提供的剪切波图像生成方法的任一方法。
在示例性实施例中，本申请提供的一种剪切波图像生成方法的各个方面还可以实现为一种程序产品的形式，其包括程序代码，当程序产品在计算机设备上运行时，程序代码用于使计算机设备执行本说明书上述描述的根据本申请各种示例性实施方式的剪切波图像生成方法中的步骤。
程序产品可以采用一个或多个可读介质的任意组合。可读介质可以是可读信号介质或者可读存储介质。可读存储介质例如可以是但不限于电、磁、光、电磁、红外线、或半导体的系统、装置或器件，或者任意以上的组合。可读存储介质的更具体的例子(非穷举的列表)包括：具有一个或多个导线的电连接、便携式盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。
本申请的实施方式的用于剪切波图像生成方法的程序产品可以采用便携式紧凑盘只读存储器(CD-ROM)并包括程序代码，并可以在电子设备上运行。然而，本申请的程序产品不限于此，在本文件中，可读存储介质可以是任何包含或存储程序的有形介质，该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。
可读信号介质可以包括在基带中或者作为载波一部分传播的数据信号，其中承载了可读程序代码。这种传播的数据信号可以采用多种形式，包括但不限于电磁信号、光信号或上述的任意合适的组合。可读信号介质还可以是可读存储介质以外的任何可读介质，该可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。
可读介质上包含的程序代码可以用任何适当的介质传输，包括但不限于无线、有线、光缆、RF等等，或者上述的任意合适的组合。可以以一种或多种程序设计语言的任意组合来编写用于执行本申请操作的程序代码，程序设计语言包括面向对象的程序设计语言诸如Java、C++等，还包括常规的过程式程序设计语言诸如“如”语言或类似的程序设计语言。程序代码可以完全地在用户电子设备上执行、部分地在用户设备上执行、作为一个独立的软件包执行、部分在用户电子设备上部分在远程电子设备上执行、或者完全在远程电子设备或服务端上执行。在涉及远程电子设备的情形中，远程电子设备可以通过任意种类的网络包括局域网(LAN)或广域网(WAN)连接到用户电子设备，或者，可以连接到外部电子设备(例如利用因特网服务提供商来通过因特网连接)。
应当注意，尽管在上文详细描述中提及了装置的若干单元或子单元，但是这种划分仅仅是示例性的并非强制性的。实际上，根据本申请的实施方式，上文描述的两个或更多单元的特征和功能可以在一个单元中具体化。反之，上文描述的一个单元的特征和功能可以进一步划分为由多个单元来具体化。
此外，尽管在附图中以特定顺序描述了本申请方法的操作，但是，这并非要求或者暗示必须按照该特定顺序来执行这些操作，或是必须执行全部所示的操作才能实现期望的结果。附加地或备选地，可以省略某些步骤，将多个步骤合并为一个步骤执行，和/或将一个步骤分解为多个步骤执行。
本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
显然，本领域的技术人员可以对本申请进行各种改动和变型而不脱离本申请的精神和范围。这样，倘若本申请的这些修改和变型属于本申请权利要求及其等同技术的范围之内，则本申请也意图包含这些改动和变型在内。

Claims (10)

1.一种剪切波图像生成方法，其特征在于，所述方法包括：
利用声辐射力或低频振动在对象的目标部位上进行激励产生剪切波，在剪切波传播路径上对所述剪切波进行检测，获取所述剪切波的检测数据；
基于所述检测数据，确定所述剪切波的各个质点在不同时刻分别对应的质点振动速度；
针对任一个质点在不同时刻对应的质点振动速度，对所述质点在不同时刻对应的质点振动速度进行多个方向的滤波处理，得到每个质点振动速度在各个方向对应的参考质点振动速度；
针对任一个方向，根据所述方向上各个质点的参考质点振动速度，确定所述方向上每个质点对应的剪切波速度；
根据预设的超声黑白图像，以及各个方向上每个质点对应的剪切波速度，生成剪切波图像。
2.根据权利要求1所述的方法，其特征在于，所述基于所述检测数据，确定所述剪切波的各个质点在不同时刻分别对应的质点振动速度，包括：
针对任意一个质点，基于所述检测数据中所述质点对应的数据，利用二维自相关算法确定所述质点在每个时刻对应的预设时长内的轴向位移；
基于每个时刻对应的预设时长内的轴向位移以及所述预设时长，确定所述质点在不同时刻分别对应的质点振动速度。
3.根据权利要求1所述的方法，其特征在于，所述对所述质点在不同时刻对应的质点振动速度进行多个不同方向的滤波处理，得到每个质点振动速度在各个方向对应的参考质点振动速度，包括：
对所述质点在不同时刻对应的质点振动速度分别进行傅里叶变换，得到每个质点振动速度对应的频域波数域数据；
针对任一个频域波数域数据，将所述频域波数域数据与所述多个方向对应的预设掩膜分别进行相乘；并将各个乘积均进行傅里叶逆变换，得到所述频域波数域数据在各个方向对应的参考质点振动速度。
4.根据权利要求3所述的方法，其特征在于，所述预设掩膜为边缘采用截取函数进行截断之后的掩膜。
5.根据权利要求1所述的方法，其特征在于，所述根据所述方向上各个质点的参考质点振动速度，确定所述方向上每个质点对应的剪切波速度，包括：
针对所述方向上的任一个质点分别执行以下操作：
基于所述质点的位置信息以及候选质点的位置信息，确定所述质点与所述候选质点之间的距离；所述候选质点为在所述方向上与所述质点相邻的一个或多个连续的质点；
基于所述质点与所述候选质点在所述方向上的参考质点振动速度，以及所述质点与所述候选质点之间的距离，确定所述质点对应的剪切波速度。
6.根据权利要求5所述的方法，其特征在于，所述基于所述质点与所述候选质点在所述方向上的参考质点振动速度，以及所述质点与所述候选质点之间的距离，确定所述质点对应的剪切波速度，包括：
基于所述质点在所述方向上的参考质点振动速度的峰值点对应的第一时间，以及所述候选质点在所述方向上的参考质点振动速度的峰值点对应的第二时间，确定所述第一时间和所述第二时间之间的时间差值；
基于所述质点与所述候选质点之间的距离，以及所述第一时间和所述第二时间之间的时间差值，得到所述质点对应的剪切波速度。
7.根据权利要求1所述的方法，其特征在于，所述利用声辐射力或低频振动在对象的目标部位上进行激励产生剪切波，在剪切波传播路径上对所述剪切波进行检测，获取所述剪切波的检测数据，包括：
利用声辐射力或低频振动在对象的目标部位上进行激励产生剪切波，并在剪切波传播路径上多次发射超声信号对所述剪切波进行检测，获取回波数据；
对所述回波数据进行正交解调处理和降频处理，得到所述剪切波的检测数据。
8.一种超声设备，其特征在于，所述超声设备包括至少一个处理器、以及至少一个存储器；其中，所述存储器存储有程序代码，当程序代码被所述处理器执行时，使得所述处理器执行下列过程：
利用声辐射力或低频振动在对象的目标部位上进行激励产生剪切波，在剪切波传播路径上对所述剪切波进行检测，获取所述剪切波的检测数据；
基于所述检测数据，确定所述剪切波的各个质点在不同时刻分别对应的质点振动速度；
针对任一个质点在不同时刻对应的质点振动速度，对所述质点在不同时刻对应的质点振动速度进行多个方向的滤波处理，得到每个质点振动速度在各个方向对应的参考质点振动速度；
针对任一个方向，根据所述方向上各个质点的参考质点振动速度，确定所述方向上每个质点对应的剪切波速度；
根据预设的超声黑白图像，以及各个方向上每个质点对应的剪切波速度，生成剪切波图像。
9.根据权利要求8所述的超声设备，其特征在于，在执行所述对所述质点在不同时刻对应的质点振动速度进行多个不同方向的滤波处理，得到每个质点振动速度在各个方向对应的参考质点振动速度时，所述处理器具体用于：
对所述质点在不同时刻对应的质点振动速度分别进行傅里叶变换，得到每个质点振动速度对应的频域波数域数据；
针对任一个频域波数域数据，将所述频域波数域数据与所述多个方向对应的预设掩膜分别进行相乘；并将各个乘积均进行傅里叶逆变换，得到所述频域波数域数据在各个方向对应的参考质点振动速度。
10.一种剪切波图像生成装置，其特征在于，所述装置包括：
获取模块，用于利用声辐射力或低频振动在对象的目标部位上进行激励产生剪切波，在剪切波传播路径上对所述剪切波进行检测，获取所述剪切波的检测数据；
质点振动速度确定模块，用于基于所述检测数据，确定所述剪切波的各个质点在不同时刻分别对应的质点振动速度；
参考质点振动速度确定模块，用于针对任一个质点在不同时刻对应的质点振动速度，对所述质点在不同时刻对应的质点振动速度进行多个方向的滤波处理，得到每个质点振动速度在各个方向对应的参考质点振动速度；
剪切波速度确定模块，用于针对任一个方向，根据所述方向上各个质点的参考质点振动速度，确定所述方向上每个质点对应的剪切波速度；
图像生成模块，用于根据预设的超声黑白图像，以及各个方向上每个质点对应的剪切波速度，生成剪切波图像。
CN202211048860.3A 2021-10-27 2022-08-30 一种剪切波图像生成方法、超声设备及装置 Pending CN115414065A (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202211048860.3A CN115414065A (zh) 2022-08-30 2022-08-30 一种剪切波图像生成方法、超声设备及装置
PCT/CN2022/116912 WO2023071528A1 (zh) 2021-10-27 2022-09-02 剪切波传播速度的确定方法及超声设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211048860.3A CN115414065A (zh) 2022-08-30 2022-08-30 一种剪切波图像生成方法、超声设备及装置

Publications (1)

Publication Number Publication Date
CN115414065A true CN115414065A (zh) 2022-12-02

Family

ID=84199544

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211048860.3A Pending CN115414065A (zh) 2021-10-27 2022-08-30 一种剪切波图像生成方法、超声设备及装置

Country Status (1)

Country Link
CN (1) CN115414065A (zh)

Similar Documents

Publication Publication Date Title
US8864671B2 (en) Methods and systems for color flow imaging
CN113040878B (zh) 超声穿刺针的位置信息处理方法、超声设备及存储介质
JP2010075586A (ja) 超音波診断装置
EP3429476B1 (en) Shear wave group velocity estimation using spatiotemporal peaks and amplitude thresholding
CN113077394A (zh) 图像处理方法、电子设备和存储介质
US11517289B2 (en) Method and apparatus for acquiring motion information
JP2007202617A (ja) 超音波診断装置
EP3656310B1 (en) Method and device for quantifying medium viscoelasticity
US20090030319A1 (en) Doppler signal spectrum calculating method and apparatus
CN104127207A (zh) 基于超声成像检测小血流信号灵敏度的方法及系统
Kranemann et al. Real-time magnetomotive ultrasound imaging using a recursive estimator
US20130237828A1 (en) Method and apparatus for obtaining movement velocity and direction of tissue
KR20130037112A (ko) 진단영상 생성장치, 의료영상시스템 및 빔포밍 수행방법
JP2014036735A (ja) 超音波診断装置
JP2010125203A (ja) 超音波診断装置
Song et al. Combining autocorrelation signals with delay multiply and sum beamforming algorithm for ultrasound imaging
CN115414065A (zh) 一种剪切波图像生成方法、超声设备及装置
KR101570194B1 (ko) 조직의 이동 속도 및 방향 획득 방법 및 장치
CN103778620A (zh) 医学图像的处理方法和装置
JP2013174473A (ja) 水中探知装置およびプログラム
JP4918344B2 (ja) 画像生成方法および超音波診断装置
Hermawati et al. Ultrasound Image Formation from Doppler Transducer
CN113545806A (zh) 前列腺弹性成像方法和超声弹性成像系统
Rossi et al. Toward the Real Time Implementation of the 2-D Frequency-Domain Vector Doppler Method
WO2023071528A1 (zh) 剪切波传播速度的确定方法及超声设备

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