CN110710989A - 弹性成像方法、系统及计算机可读存储介质 - Google Patents

弹性成像方法、系统及计算机可读存储介质 Download PDF

Info

Publication number
CN110710989A
CN110710989A CN201910924688.5A CN201910924688A CN110710989A CN 110710989 A CN110710989 A CN 110710989A CN 201910924688 A CN201910924688 A CN 201910924688A CN 110710989 A CN110710989 A CN 110710989A
Authority
CN
China
Prior art keywords
propagation path
processor
controlling
main propagation
under test
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
CN201910924688.5A
Other languages
English (en)
Other versions
CN110710989B (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.)
Shenzhen Mindray Bio Medical Electronics Co Ltd
Original Assignee
Shenzhen Mindray Bio Medical Electronics 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 Shenzhen Mindray Bio Medical Electronics Co Ltd filed Critical Shenzhen Mindray Bio Medical Electronics Co Ltd
Priority to CN201910924688.5A priority Critical patent/CN110710989B/zh
Priority to CN202111470084.1A priority patent/CN114287968A/zh
Priority to CN202110659607.0A priority patent/CN113261991A/zh
Priority to PCT/CN2019/112507 priority patent/WO2021056645A1/zh
Priority to CN201980098269.9A priority patent/CN114072066B/zh
Publication of CN110710989A publication Critical patent/CN110710989A/zh
Application granted granted Critical
Publication of CN110710989B publication Critical patent/CN110710989B/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/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0833—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures
    • A61B8/085—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures for locating body or organic structures, e.g. tumours, calculi, blood vessels, nodules
    • 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/4444—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device related to the probe
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/46—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
    • A61B8/467—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient characterised by special input means
    • A61B8/469—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient characterised by special input means for selection of a region of interest
    • 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

Landscapes

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

Abstract

本申请实施例公开了一种弹性成像方法、系统及计算机可读存储介质。该弹性成像方法包括：获取受测组织的运动参数图像；确定所述运动参数图像中于所述受测组织内传播的剪切波的主传播路径；显示包含所述主传播路径的运动参数图像。本申请实施例通过获取运动参数图像中剪切波的主传播路径后显示包含主传播路径的运动参数图像，如此，可减少运动参数图像中因探头的振动会造成的各种余波、反射波等干扰信息的影响而造成医护人员无法较准确地解释运动参数图像所显示的意义的不足，亦有利于提高弹性测试结果的可理解性。

Description

弹性成像方法、系统及计算机可读存储介质
技术领域
本申请涉及医疗技术领域，尤其涉及一种弹性成像方法、系统及计算机可读存储介质。
背景技术
瞬时弹性成像利用探头振动产生剪切波传入受测组织，并发射超声波检测组织内部位移，从而计算并显示受测组织的弹性参数。在瞬时弹性技术中，除了给出弹性检测结果外，一般还会提供一个组织位移或应变的运动参数图像。然而，由于探头的振动会造成的各种余波、反射波等干扰信息，造成临床人员无法较准确地解释运动参数图像所显示的意义。
发明内容
本申请实施例提供了一种弹性成像方法、系统及计算机可读存储介质，可以提高弹性测试结果的可理解性。
一个实施例中，提供了一种弹性成像方法，应用于弹性成像系统，所述弹性成像系统包括探头、连接于所述探头的发射电路、连接于所述探头的接收电路、连接于所述接收电路的波束合成器、连接于所述波束合成器的处理器和显示所述处理器所传输的图像信息的显示屏，所述弹性成像方法包括：
控制所述探头在接收到所述发射电路的第一发射时序时向受测组织发射第一超声波，以跟踪在所述受测组织内传播的剪切波；
控制所述探头接收所述受测组织返回的第一超声回波，并将所述第一超声回波转换为电信号后传输至所述接收电路；
控制所述波束合成器对所述接收电路传输来的电信号进行波束合成，获得第一超声回波数据；
控制所述处理器基于所述第一超声回波数据获得所述受测组织的运动参数图像；
控制所述处理器确定所述运动参数图像中所述剪切波的主传播路径，获得主传播路径图，其中所述主传播路径图表征所述剪切波在所述受测组织中的主传播路径；
控制所述处理器根据所述第一超声回波数据或者所述运动参数图像或者所述主传播路径图确定所述受测组织的弹性信息；
控制所述处理器将所述主传播路径图和所述受测组织的弹性信息显示于所述显示屏内。
一个实施例中，提供了一种弹性成像方法，该方法包括：
向受测组织发射第一超声波，以跟踪在所述受测组织内传播的剪切波；
接收所述受测组织返回的第一超声回波，获得第一超声回波数据；
根据所述第一超声回波数据获得由所述剪切波在所述受测组织中的传播引起的所述受测组织在不同时刻及不同深度下的运动参数；
根据所述运动参数确定所述剪切波的主传播路径，获得主传播路径图，其中所述主传播路径图表征所述剪切波在所述受测组织中的主传播路径；
根据所述第一超声回波数据或者所述运动参数或者所述主传播路径确定所述受测组织的弹性信息；
显示所述主传播路径图和所述受测组织的弹性信息。
一个实施例中，提供了一种弹性成像方法，该方法包括：
向受测组织发射第一超声波，以跟踪在所述受测组织内传播的剪切波；
接收所述受测组织返回的第一超声回波，获得第一超声回波数据；
根据所述第一超声回波数据获得由所述剪切波在所述受测组织中的传播引起的所述受测组织在不同时刻及不同深度下的运动参数；
根据所述运动参数确定所述剪切波的主传播路径，获得主传播路径图，其中所述主传播路径图表征所述剪切波在所述受测组织中的主传播路径；
显示所述主传播路径图。
一个实施例中，提供了一种弹性成像方法，包括：
获取受测组织的运动参数图像；
根据所述运动参数图像确定在所述受测组织内传播的剪切波的主传播路径；
显示所述主传播路径。
一个实施例中，提供了一种弹性成像系统，包括：
探头，用于向受测组织发射第一超声波，以跟踪在所述受测组织内传播的剪切波，所述探头还用于接收所述受测组织返回的第一超声回波，获得第一超声回波数据；
处理器，连接于所述探头，所述处理器用于基于所述第一超声回波数据获得由所述剪切波在所述受测组织中的传播引起的所述受测组织在不同时刻及不同深度下的运动参数，根据所述运动参数确定所述剪切波的主传播路径，获得主传播路径图，并根据所述第一超声回波数据或者所述运动参数或者所述主传播路径图确定所述受测组织的弹性信息，其中所述主传播路径图表征所述剪切波在所述受测组织中的主传播路径；
显示屏，连接于所处理器，所述处理器用于将所述主传播路径图和所述受测组织的弹性信息显示于所述显示屏内。
一个实施例中，提供了一种计算机可读存储介质，所述计算机可读存储介质用于存储电子数据交换的计算机程序，其中，所述计算机程序使得计算机执行前述实施例中任一方法中所描述的部分或全部步骤。
本申请实施例的弹性成像方法、系统及计算机可读存储介质根据运动参数或运动参数图像获得剪切波的主传播路径，并显示该主传播路径，如此，可减少运动参数图像中因探头的振动会造成的各种余波、反射波等干扰信息的影响而造成医护人员无法较准确地解释运动参数图像所显示的意义的不足，亦有利于提高弹性测试结果的可理解性。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本申请一实施例中弹性成像系统的硬件结构示意图。
图2是本申请一实施例中弹性成像方法的步骤流程图。
图3是本申请一实施例中运动参数图像的示意图。
图4是本申请一实施例中探头的硬件结构框图。
图5是本申请一实施例中运动参数图像中多个带状区域的示意图。
图6是本申请一实施例中主传播路径图的二值化的示意图。
图7是本申请一实施例中主传播路径图的非二值化的示意图。
图8是本申请一实施例中主传播路径的示意图。
图9是本申请一实施例中弹性成像方法的步骤流程图。
图10是本申请一实施例中的弹性成像系统的框图示意图。
图11是本申请一实施例中运动参数图像的示意图。
图12是本申请又一实施例中运动参数图像的示意图。
图13是本申请再一实施例中运动参数图像的示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
本申请的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别不同对象，而不是用于描述特定顺序。此外，术语“包括”和“具有”以及它们任何变形，意图在于覆盖不排他的包含。例如包含了一系列步骤或单元的过程、方法、系统、产品或设备没有限定于已列出的步骤或单元，而是可选地还包括没有列出的步骤或单元，或可选地还包括对于这些过程、方法、产品或设备固有的其他步骤或单元。
在本文中提及“实施例”意味着，结合实施例描述的特定特征、结构或特性可以包含在本申请的至少一个实施例中。在说明书中的各个位置出现该短语并不一定均是指相同的实施例，也不是与其它实施例互斥的独立的或备选的实施例。本领域技术人员显式地和隐式地理解的是，本文所描述的实施例可以与其它实施例相结合。
请参阅图1，所示为本申请一实施例中弹性成像系统的硬件结构示意图。所述弹性成像系统10可包括探头100、连接探头100的发射电路102、连接探头 100的接收电路104、波束合成器106、处理器110及显示器112，其中，接收电路104、波束合成器106、处理器110及显示屏112可依次电性连接。本实施例中，弹性成像系统10可以获取受测组织的运动参数或运动参数图像，并根据该运动参数或运动参数图像获得剪切波在受测组织中的主传播路径，并可以在显示屏112内显示该主传播路径。由于主传播路径可以是可以准确代表不同深度下剪切波的传播位置的单一路径，在获取该主传播路径时可以排除各种余波、反射波等干扰信息，如此，医护人员直观地通过包含主传播路径的运动参数图像或者通过该主传播路径图来进行诊断，有利于提高弹性测试结果的可理解性。在一实施例中，波束合成器106及处理器110可由专门的电路或者可商业获得的芯片来实现。
请一并参阅图2，所示为本申请一实施例中弹性成像方法的步骤流程图。该弹性成像方法包括如下步骤：
步骤200，获取受测组织的运动参数图像。
本实施例中，发射电路102向探头100传输第一发射时序，以控制探头100 向受测组织发射第一超声波，其中，第一超声波用于跟踪于受测组织内传播的剪切波。在探头100向受测组织发射第一超声波后，在经一定延时，探头100 可接收从受测组织反射回来的带有检测对象的信息的第一超声回波。探头100 可将此超声回波转换为电信号。接收电路104接收探头100转换生成的电信号，获得第一超声回波数据，并将这些第一超声回波数据送入波束合成器106。波束合成器106对超声回波数据进行聚焦延时、加权和通道求和等波束合成处理，然后将波束处理后的超声回波数据送入处理器110，由处理器110根据第一超声回波数据获取受测组织的运动参数或运动参数图像，并用以在显示屏112上显示。
请一并参阅图3，所示为本申请一实施例中运动参数图像的示意图。该运动参数图像包括横向的时间属性及纵向的深度属性。本实施例中，在剪切波传入受测组织后，随着剪切波的传播，受测组织内部会发生振动，振动使受测组织的相应位置发生位移。通过持续一段时间向受测组织中发射第一超声波并接收其回波。
处理器110基于第一超声回波数据获得由剪切波在受测组织中的传播而引起的受测组织在不同时刻及不同深度下的运动参数，其中，这里的运动参数可以包括位移、速度或者应变。例如，处理器110将不同时刻所得第一超声回波数据进行比较分析(比如互相关算法)，可以计算出不同时刻下受测组织的位移，对来自不同深度的受测组织的第一回波数据分别进行位移计算，最终可以得到一个不同深度-不同时刻一一对应的位移矩阵。在位移矩阵中，每个数据各自代表某个深度的受测组织在某个时刻下的位移信息。在对上述位移矩阵沿着深度方向求梯度时，即可相应的得到应变矩阵。在应变矩阵中，每个数据各自代表某个深度的受测组织在某个时刻下的应变信息。在上述计算过程中，为了提高信噪比，还有可能加入一些时间方向或者深度方向的滤波操作。
处理器110可基于剪切波在不同时刻及不同深度下的运动参数确定受测组织的运动参数图像150。
请一并参阅图4，所示本申请一实施例中探头的硬件结构框图。探头100包括阵列型声头130、振动器132及位于阵列型声头130与振动器132之间的传感器134。在发射电路102向探头100传输第一发射时序之前，发射电路102可向探头100传输激励时序，以控制探头100的振动器132振动并在受测组织内产生剪切波。之后，探头100的阵列型声头130根据第一发射时序跟踪在受测组织内传播的剪切波。阵列式声头130包括预设数量的阵元，阵列式声头130的阵元的排列的方式为直线式排列或者扇形式排列等。传感器132用于感测探头100按压受测组织的力度。在一实施例中，探头100亦可不包括传感器134。
在一实施例中，医护人员可能需要针对受测组织其一目标位置范围进行检测，因此，医护人员需要在基础图像中选择该目标位置范围所对应的感兴趣区域，其中，基础图像包括B图像、C图像中的一种或多种。在获取受测组织的基础图像时，发射电路102向探头100传输第二发射时序，以控制探头100向受测组织发射第二超声波。在探头100向受测组织发射第二超声波后，在经一定延时，探头100可接收从受测组织反射回来的带有检测对象的信息的第二超声回波。探头100可将此超声回波转换为电信号。接收电路104接收探头100 转换生成的电信号，获得第二超声回波数据，并将这些第二超声回波数据送入波束合成器106。波束合成器106对超声回波数据进行聚焦延时、加权和通道求和等波束合成处理，然后将波束处理后的超声回波数据送入处理器110，由处理器110根据用户所需的成像模式的不同对信号进行不同的处理，以获得不同模式的组织图像数据，然后经对数压缩、动态范围调整、数字扫描变换等处理形成不同模式的超声组织图像，并用以在显示器112上显示，其中，不同模式的超声组织图像可包括M图像、B图像、C图像等，或者其他类型的二维超声组织图像或三维超声组织图像。在一实施例中，探头100发射的第一超声波与第二超声波可相同，即处理器110通过对探头100接收到的超声回波处理后可同时得到剪切波所对应的参数信息、生成瞬时弹性图及生成不同模式的超声组织图像；在一实施例中，探头100发射的第一超声波与第二超声波可不相同，即探头100可先后发射第一超声波及第二超声波、或者先后发射第二超声波及第一超声波、或者穿插式的发射第一超声波及第二超声波(如发射第一超声波后发射第二超声波，之后再发射第一超声波，如此穿插式的反复循环进行)，如此，处理器110可通过对探头100接收到对应第一超声波的第一超声回波处理后可得到剪切波所对应的参数信息、生成瞬时弹性图，并通过对探头100接收到对应第二超声波的第二超声回波处理后可生成不同模式的超声组织图像。
当显示屏112内显示基础图像时，医护人员可在基础图像中确定的感兴趣区域；处理器110可获取第一超声回波数据中与所述感兴趣区域相对应的目标位置范围，并基于目标位置范围内的第一超声回波数据确定受测组织的弹性信息，如受测组织的目标位置范围内的剪切波传播速度、剪切模量、杨氏模量等。
步骤202，确定所述运动参数图像中于所述受测组织内传播的剪切波的主传播路径，获得主传播路径图。
本实施例中，处理器110确定运动参数图像150中每一深度下位于预设范围的运动参数所对应的目标区域，并可确定目标区域在时间属性上的目标时刻范围，其中，目标时刻范围包括若干目标时刻。处理器110基于运动参数图像 150中相连续的目标区域确定带状区域。由于探头的振动会造成的各种余波、反射波等干扰信息的影响，处理器110在基于运动参数图像150中相连续的目标区域获取带状区域时可能得到一个或多个带状区域。
请一并参阅图5，所示为本申请一实施例中运动参数图像中多个带状区域的示意图。其中，运动参数图像150可包括第一带状区域S1、第二带状区域S2 及第三带状区域S3。本实施例中，运动参数图像150包括若干像素点。由于运动参数图像150中每一像素点的像素值与该像素点所对应的深度下的运动参数的大小相对应。例如，在运动参数图像150为非灰度图像(如运动参数图像150 为伪彩色图像)时，处理器110可对运动参数图像150进行灰度化处理。当运动参数所对应的数值较大时，在灰度化后的运动参数图像中，对应的深度下的像素点接近于白色(如像素点的像素值接近于255)；当运动参数所对应的数值较小时，在灰度化后的运动参数图像中，对应的深度下的像素点接近于黑色(如像素点的像素值接近于0)。在确定目标区域时，处理器110可确定最大极值范围或最小极值范围为该预设范围，其中，最大极值范围可为a至255，目标区域可为运动参数图像150(如图3)中的亮带区域；最小极值范围可为0至b，目标区域可为运动参数图像150(如图3)中的黑带区域。
在其他实施例中，目标区域亦可为运动参数图像150(如图3)中亮带区域与黑带区域之间的区域。例如，对于一设定深度下，当一像素点的像素值位于最小极值范围，且与该像素点间隔预设数量的其他像素点的像素值均位于最大极值范围，此时，处理器110可将该像素点及与该像素点间隔预设数量的其他像素点可作为该设定深度下所对应的目标区域；或者，当一像素点的像素值位于最大极值范围，且与该像素点间隔预设数量的其他像素点的像素值均位于最小极值范围，此时，处理器110可将该像素点及与该像素点间隔预设数量的其他像素点可作为该设定深度下所对应的目标区域。
本实施例中，在确定运动参数图像150中的第一带状区域S1、第二带状区域S2及第三带状区域S3时，若预设范围为最小极值范围，处理器110确定深度为V1时所对应的目标区域，包括目标区域AB及目标区域EF，其中，目标区域AB为运动参数图像150中在深度为V1下运动参数位于该最小极值范围内的像素点集合(如线段AB)，目标区域AB所对应的目标时刻范围为t1至t2目标时刻；目标区域EF为运动参数图像150中在深度为V1下运动参数位于该最小极值范围内的像素点集合(如线段EF)，目标区域EF所对应的目标时刻范围为t5至t6目标时刻；而运动参数图像150中在深度为V1时目标区域AB及目标区域EF之外的其他区域不满足最小极值范围。处理器110还可确定深度为 V2时所对应的目标区域，包括目标区域CD，其中，目标区域CD为运动参数图像150中在深度为V2下运动参数位于该最小极值范围内的像素点集合(如线段CD)，目标区域CD所对应的目标时刻范围为t3至t4目标时刻，而运动参数图像150中在深度为V2时目标区域CD之外的其他区域不满足最小极值范围。
处理器110基于运动参数图像150中相连续的目标区域确定一个或多个带状区域。由于剪切波在受测组织中是连续的，因此，不同深度下的目标区域是相连续的。如此，处理器110在运动参数图像150中相连续的目标区域所组成的带状区域，如处理器110确定出了运动参数图像150中的第一带状区域S1、第二带状区域S2及第三带状区域S3。
由于剪切波在受测组织内传播时，其对应的传播路径是单一的。若处理器 110确定运动参数图像150中存在多个带状区域时，表示运动参数图像150中存在干扰信息，因此，处理器110可确定多个带状区域中满足预设条件的目标带状区域为剪切波的主传播路径。
在一实施例中，由于剪切波是在探头100振动结束后产生，因此，处理器 110可获取探头100振动结束时所对应的参考时刻，并确定该一个或多个带状区域中目标时刻晚于参考时刻的目标区域所组成的带状区域为主传播路径。例如，若探头100振动结束时所对应的参考时刻为t0时，由于第三带状区域S3所对应的目标时刻均早于参考时刻t0，因此，处理器110确定第一带状区域S1及第二带状区域S2中目标时刻晚于参考时刻t0的目标区域。由于运动参数图像150 中还包括第一带状区域S1及第二带状区域S2，而运动参数图像150的主传播路径为一个带状区域。此时，处理器110可确定第一带状区域S1及第二带状区域 S2中具有最大长度或者最大面积的带状区域为剪切波的主传播路径，其中，每一带状区域包括第一斜边及第二斜边，带状区域的长度可表示为第一斜边或第二斜边的长度，或者是第一斜边与第二斜边中较长者；带状区域的面积可表示为第一斜边、第二斜边、第一斜边的第一端与第二斜边第一端在时间轴上的投影之差及第一斜边的第二端与第二斜边第二端在时间轴上的投影之差所围成的四边形的面积。由于第一带状区域S1中AC所在的斜边的长度大于第二带状区域S2中E所在的斜边的长度，因此，处理器110可确定第一带状区域为剪切波的主传播路径。
在一实施例中，处理器110可直接确定一个或多个带状区域中具有最大长度或者最大面积的带状区域为剪切波的主传播路径。例如，在第一带状区域S1、第二带状区域S2及第三带状区域S3中，第一带状区域S1具有最大的长度及最大面积，因此，处理器110可确定第一带状区域为剪切波的主传播路径。
在一实施例中，当处理器110确定该一个或多个带状区域中目标时刻晚于参考时刻的目标区域所组成的带状区域的数量为一个时，处理器110亦可无需再确定带状区域的长度或面积的属性信息。
在一实施例中，在确定受测组织的感兴趣区域时，感兴趣区域可能位于预设深度下，此时，处理器110可确定位于该预设深度下的目标带状区域为主传播路径。例如，当预设深度为深度V1时，处理器110可确定第一带状区域S1 中线段AB以下的目标区域满足条件、第二带状区域S2均满足条件、第三带状区域S3中线段EF以下的目标区域满足条件，处理器110可亦可得到三个带状区域。由于运动参数图像150中还包括第一带状区域S1及第二带状区域S2，而运动参数图像150的主传播路径为一个带状区域。此时，处理器110可确定第一带状区域S1中线段AB以下的目标区域、第二带状区域S2、第三带状区域 S3中线段EF以下的目标区域中具有最大长度或者最大面积的带状区域为剪切波的主传播路径，如处理器110可确定第一带状区域S1中线段AB以下的目标区域为剪切波的主传播路径。
在一实施例中，当确定位于该预设深度下的目标带状区域存在多个带状区域时，处理器110亦要结合探头100振动结束时所对应的参考时刻来确定目标带状区域。例如，由于第二带状区域S2所对应的目标时刻均早于参考时刻t0，因此，处理器110可确定第一带状区域S1中线段AB以下的目标区域、及第三带状区域S3中线段EF以下的目标区域满足条件，之后，处理器110可确定第一带状区域S1中线段AB以下的目标区域、及第三带状区域S3中线段EF以下目标区域中具有最大长度或者最大面积的带状区域为剪切波的主传播路径，即处理器110可确定第一带状区域S1中线段AB以下的目标区域为剪切波的主传播路径。
本文中，所说的主传播路径图是表征剪切波在受测组织中的实际主要传播路径的图，其基于由剪切波在受测组织中的传播引起的该受测组织在不同时刻及不同深度下的运动参数或者运动参数图而获得。
图11、图12和图13分别为本发明一些实施例中的由剪切波在受测组织中的传播引起的该受测组织在不同时刻及不同深度下的运动参数或者运动参数图。图11中，右侧所展示的向右上方延伸的带状区域即为剪切波的反射波，图12中虚线右侧的带状区域即为剪切波的余波，图13中虚线右侧的带状和长圆形区域即为剪切波的反射波和余波。可见，获得的运动参数或者运动参数图中，包含了剪切波的各种余波和反射波等干扰信息，会影响医生的查看和弹性参数的检测。
本文的实施例中，如前文以及后文所述，基于该运动参数或运动参数图获得表征剪切波在受测组织中的实际主要传播路径的主传播路径图，该主传播路径图中可以排除该运动参数或者运动参数图中的余波、反射波等干扰信息，从而能够更准确地反映剪切波在受测组织中的实际主要传播路径，便于医生查看。
步骤204，显示包含所述主传播路径的运动参数图像。
在处理器110获取到运动参数图像150中的剪切波的主传播路径时，处理器110通过显示屏112显示主传播路径160(示于图6)，其中，主传播路径图亦为包含主传播路径的运动参数图像。由于主传播路径图中包含主传播路径所对应的显示区域及主传播路径之外的其他区域，因此，主传播路径图中已排除了干扰信息的影响，提高弹性测试结果的可理解性，亦可使得医护人员可直观地根据主传播路径图进行诊断。
一个实施例中，处理器110还可以根据前述的第一超声回波数据或者前述的运动参数图像或者前述的主传播路径图确定受测组织的弹性信息。这里，受测组织的弹性信息可以是剪切波传播速度、受测组织的杨氏模量、受测组织的剪切模量等等参数。
请一并参阅图6，所示本申请一实施例中主传播路径图的二值化的示意图。在获取运动参数图像150中的主传播路径后，处理器110可对该运动参数图像 150进行二值化处理，以使得运动参数图像150中的主传播路径显示为第一颜色、主传播路径之外的区域显示为第二颜色，更好便于医护人员的辨识。例如，处理器110可将运动参数图像150中主传播路径所对应的目标带状区域设置为第一颜色(如白色)，并将主传播路径所对应的目标带状区域之外的区域设置为第二颜色(如黑色)。
在一实施例中，若干扰信息的影响较弱，当处理器110对运动参数图像150 进行二值化处理后，得到的主传播路径图中可区分出主传播路径的显示区域及主传播路径之外的其他区域。例如，当预设范围为最小极值范围时，处理器110 将运动参数图像150中像素值大于预设阈值的像素点设置为第一颜色，即将主传播路径之外的其他区域的颜色设置为黑色；处理器110还将运动参数图像中像素值不大于该预设阈值的像素点设置为第二颜色，即将主传播路径的显示区域的颜色设置为白色。如此，在处理器110对运动参数图像150进行二值化处理后可直接得到主传播路径160。
请一并参阅图7，所示为本申请一实施例中主传播路径图的非二值化的示意图。在一实施例中，由于剪切波在受测组织内传播时，随着传播的深度及传播的时间增加，受测组织的运动参数(如应变信息或位移信息)或能量逐渐减弱，为更清楚地显示这一变化过程，处理器110可将主传播路径160进行非二值化显示。处理器110可将主传播路径所对应的目标带状区域中每一深度处的像素点的像素值设置为具有对应运动参数的第三颜色，其中，不同的运动参数所对应的第三颜色不同；并可将主传播路径所对应的目标带状区域之外的区域设置为第四颜色。
请一并参阅图8，所示本申请又一实施例中主传播路径的示意图。当处理器 110确定第一带状区域S1为主传播路径时，处理器110亦可将主传播路径上某些特定的深度的点来做简易的显示。如处理器110可在显示屏112内只显示图8中的圆圈(不包括第一带状区域S1的两个斜边)，以进行简易的显示。
在一实施例中，在确定运动参数图像150中的主传播路径时，处理器110 可计算出受测组织的弹性参数，包括但不限于剪切波速度、杨氏模量、剪切模量等。如图5所示，处理器110可基于主传播路径拟合出图中虚线所示的直线，其中，虚线所示的直线的斜率可用于表示剪切波速度；处理器110可甘于剪切波速度来计算出受测组织的杨氏模量等。在一实施例中，处理器110也包括上述参数的多次测量统计结果等，如10次测量所得杨氏模量的中位数、四分位数、四分位数与中位数的比值等。
在一实施例中，处理器110可控制显示屏112显示基础图像及位于基础图像内的感兴趣区域、受测组织的杨氏模量及/或主传播路径图，以便医护人员进行诊断。
上述弹性成像方法通过获取运动参数图像中剪切波的主传播路径后显示包含主传播路径的运动参数图像，如此，可减少运动参数图像中因探头的振动会造成的各种余波、反射波等干扰信息的影响而造成医护人员无法较准确地解释运动参数图像所显示的意义的不足，亦有利于提高弹性测试结果的可理解性。
一个实施例中，一种弹性成像方法可以包括下列步骤：
向受测组织发射第一超声波，以跟踪在受测组织内传播的剪切波；
接收受测组织返回的该第一超声波的超声回波(本文称之为第一超声回波)，获得第一超声回波数据；
根据第一超声回波数据获得由该剪切波在该受测组织中的传播引起的该受测组织在不同时刻及不同深度下的运动参数，这里，该运动参数可以包括受测组织的位移、速度或者应变；
根据该运动参数确定该剪切波的主传播路径，并基于该主传播路径生成主传播路径图；
显示获得的主传播路径图。
本实施例中，还可以根据前述的第一超声回波数据或者该运动参数或者该主传播路径图确定该受测组织的弹性信息，并将该受测组织的弹性信息也显示给用户。该弹性信息可以是与主传播路径图同时显示，也可以是与主传播路径图非同时显示。
本实施例中，可以从这些运动参数中确定满足预设条件的运动参数(这里称之为“目标运动参数”)，并根据该目标运动参数对应的深度和时刻获得剪切波的主传播路径。这里，所说的预设条件可以是各种适合的条件，可以根据实际需要设置，例如，预设条件可以是大于预设阈值、在预设范围内，等等。
本实施例中，还可以根据获得的主传播路径确定受测组织的弹性信息的置信度。例如，可以根据该主传播路径确定该主传播路径所对应的一个或多个置信度参数，并根据该一个或多个置信度参数确定受测组织的弹性信息的置信度。这里，所说的置信度参数可以包括所述主传播路径的线性度参数、所述主传播路径直线拟合计算所述剪切波传播速度时的误差参数、所述主传播路径的长度参数、所述主传播路径的面积参数等等中的一个或多个。
本实施例中，还可以获取受测组织的基础图像，并将该受测组织的基础图像与前述的主传播路径图和受测组织的弹性信息同时显示。这里，所说的基础图像可以是受测组织的B图像、C图像或其他模式的图像中的一种或多种。该基础图像可以是由超声成像系统实时获得，即通过超声探头向受测组织发射超声波并接收超声回波，获得超声回波信号，并根据超声回波信号获得受测组织的基础图像；也可以从其他设备读取的预先采集获得并存储的受测组织的基础图像。
本实施例中，还可以显示该受测组织的弹性信息的置信度。
请参阅图9，所示为本申请又一实施例中弹性成像方法的步骤流程图。该弹性成像方法包括如下步骤：
步骤300，获取受测组织的运动参数图像。
本实施例中步骤300与上述实施例中步骤200类似，具体请参阅前述的步骤200。
步骤302，确定所述运动参数图像中于所述受测组织内传播的剪切波的主传播路径。
本实施例中步骤302与上述实施例中步骤202类似，具体请参阅前述的步骤202。
步骤304，显示包含所述主传播路径的运动参数图像。
本实施例中步骤304与上述实施例中步骤204类似，具体请参阅前述的步骤204。
步骤306，控制所述处理器确定所述受测组织的弹性信息的置信度。
在获取运动参数图像150中剪切波的主传播路径时，处理器110可根据主传播路径计算出受测组织的弹性信息，进而，主传播路径的可靠性影响到受测组织的弹性信息。因此，处理器150可基于主传播路径所对应的一个或多个置信度参数，并基于该一个或多个置信度参数确定受测组织的弹性信息的置信度，其中，置信度参数包括主传播路径的线性度参数、主传播路径直线拟合计算剪切波传播速度时的误差参数、主传播路径的长度参数、主传播路径的面积参数中的一个或多个。
例如，从主传播路径图上可以方便的判断出主传播路径路径的长短、能量的强弱、主传播路径是否线性。在临床肝纤维化检测时，主要通过判断剪切波在肝脏组织中平均传播速度的大小来计算出肝脏组织的杨氏模量值，从而反映肝纤维化的程度。一般来说，杨氏模量越大，说明肝组织越硬，肝纤维化程度越高。由于肝纤维化主要为弥漫性病变，所以剪切波传播时速度均匀，主传播路径表现为直线路径，且直线路径的斜率与剪切波的速度一一对应。在某些检查中，如果主传播路径太短(比如仅传播到50mm深度)，则说明剪切波能量弱或衰减较大，导致穿透力不足；如果主传播路径表现为曲线或者不太直的直线，则说明剪切波传播不均匀，那么所计算出来的斜率或者剪切波的传播速度或者组织的杨氏模量结果就可能不准确；如果主传播路径完全没有计算出来，那么说明本次检查的数据质量太差，难以获得有效的结果。
置信度可以是上述主传播路径的线性度参数(如确定主传播路径路径的线性程度)，当主传播路径的线性度越好，表示处理器110确定受测组织的弹性信息的可靠性越高；或者，置信度可以是处理器110对主传播路径进行直线拟合(如通过最小二乘法拟合)并计算拟合后直线的斜率时所对应的误差，当误差越小，表示处理器110确定受测组织的弹性信息的可靠性越高；或者，置信度可以是处理器110可确定主传播路径的长度参数，当主传播路径的长度参数越长，表示处理器110确定受测组织的弹性信息的可靠性越高；或者，置信度可以是处理器110可确定主传播路径的面积参数，面积越大，表示处理器110确定受测组织的弹性信息的可靠性越高。
在一实施例中，置信度可以是处理器110可基于上述的多个置信度参数的加权平衡所得的一个综合性评分参数，综合性评分参数越高，表示处理器110 确定受测组织的弹性信息的可靠性越高。
步骤308，控制输出与所述受测组织的弹性信息的置信度相对应的提示信息。
处理器110可将基于上述的一个或多个置信度参数确定受测组织的弹性信息的置信度显示于显示屏112内，每一置信度参数所对应的数值均具有对应的置信度，如置信度90％，如此，医护人员可根据显示的置信度来确定计算受测组织的弹性信息的可靠性。
在一实施例中，当处理器110获取得到了受测组织的基础图像后，处理器 110可将受测组织的基础图像、置信度及主传播路径图显示于显示屏112内。
在一实施例中，处理器110还可显示医护人员在基础图像上确定的感兴趣区域，并可显示根据计算得到的剪切波速度计算出受测组织的杨氏模量。
上述弹性成像方法通过主传播路径图确定受测组织的弹性信息的置信度，医护人员可方便地根据显示的置信度来确定受测组织的弹性信息的可靠性及可理解性。
请参阅图10，所示本申请又一实施例中的弹性成像系统80的框图示意图。如图10所示，所述弹性成像系统80可应用上述的各实施例。下面对本申请所提供的弹性成像系统80进行描述，所述弹性成像系统80可以包括处理器800、存储装置802、探头100、控制电路804及显示屏112，以及存储在所述存储装置802中并可向所述处理器800上运行的计算机程序(指令)，所述弹性成像系统 80还可以包括其他的硬件部分，例如通信装置、按键、键盘等，在此不再赘述。所述处理器800可通过信号线808与探头100、控制电路904存储装置802及显示屏112进行数据交换。
所述处理器800可以是中央处理单元(Central Processing Unit，CPU)，还可以是其他通用处理器、数字信号处理器(Digital Signal Processor，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现成可编程门阵列(Field-Programmable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等，所述处理器是所述弹性成像系统80的控制中心，利用各种接口和线路连接整个弹性成像系统80的各个部分。本实施例中，所述处理器800可用于实现前述图像处理模块110的全部功能，亦可集成有波束合成器106等元件的功能，具体元件的功能可参考前述实施例。例如，处理器800 可产生第一发射时序，以控制探头100产生第一超声波；处理器800可产生第二发射时序，以控制探头100产生第二超声波；处理器800可产生激励时序，并控制探头100振动后在受测组织内产生剪切波。
所述控制电路804可包括上述实施例中的发射电路102、接收电路104和/ 或波束合成器106等元件的功能，具体元件的功能可参考前述实施例。例如，控制电路804可产生第一发射时序，以控制探头100产生第一超声波；控制电路804可产生第二发射时序，以控制探头100产生第二超声波；控制电路804 可产生激励时序，并控制探头100振动后在受测组织内产生剪切波。
所述存储装置802可用于存储所述计算机程序和/或模块，所述处理器800 通过运行或执行存储在所述存储装置802内的计算机程序和/或模块，以及调用存储在存储装置802内的数据，实现上述弹性成像方法的各种功能。所述存储装置802可存储有超声回波数据，并可由处理器800进行根据超声回波数据确定剪切波的主传播路径，所述存储装置802可主要包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需的应用程序等。此外，存储装置802可以包括高速随机存取存储装置，还可以包括非易失性存储装置，例如硬盘、内存、插接式硬盘，智能存储卡(Smart Media Card,SMC)，安全数字(Secure Digital,SD)卡，闪存卡(Flash Card)、至少一个磁盘存储装置件、闪存器件、或其他易失性固态存储装置件。
所述显示屏112，可以显示用户界面(UI)、图形用户界面(GUI)，显示器112 可以包括液晶显示器(LCD)、薄膜晶体管LCD(TFT-LCD)、有机发光二极管(OLED)触摸显示器、柔性触摸显示器、三维(3D)触摸显示器等中的至少一种。
所述处理器800通过读取存储装置802中存储的可执行程序代码来运行与所述可执行程序代码对应的程序，以用于执行前面任一实施例中的弹性成像方法。
在上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详述的部分，可以参见其他实施例的相关描述。以上对本申请实施例进行了详细介绍，本文中应用了具体个例对本申请的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本申请的方法及其核心思想；同时，对于本领域的一般技术人员，依据本申请的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本申请的限制。

Claims (39)

1.一种弹性成像方法，应用于弹性成像系统，其特征在于，所述弹性成像系统包括探头、连接于所述探头的发射电路、连接于所述探头的接收电路、连接于所述接收电路的波束合成器、连接于所述波束合成器的处理器和显示所述处理器所传输的图像信息的显示屏，所述弹性成像方法包括：
控制所述探头在接收到所述发射电路的第一发射时序时向受测组织发射第一超声波，以跟踪在所述受测组织内传播的剪切波；
控制所述探头接收所述受测组织返回的第一超声回波，并将所述第一超声回波转换为电信号后传输至所述接收电路；
控制所述波束合成器对所述接收电路传输来的电信号进行波束合成，获得第一超声回波数据；
控制所述处理器基于所述第一超声回波数据获得所述受测组织的运动参数图像；
控制所述处理器确定所述运动参数图像中所述剪切波的主传播路径，获得主传播路径图，其中所述主传播路径图表征所述剪切波在所述受测组织中的主传播路径；
控制所述处理器根据所述第一超声回波数据或者所述运动参数图像或者所述主传播路径图确定所述受测组织的弹性信息；
控制所述处理器将所述主传播路径图和所述受测组织的弹性信息显示于所述显示屏内。
2.如权利要求1所述的弹性成像方法，其特征在于，所述控制所述处理器基于所述第一超声回波数据获取所述受测组织的运动参数图像，包括：
控制所述处理器基于所述第一超声回波数据获得所述剪切波在所述受测组织中的传播引起的所述受测组织在不同时刻及不同深度下的运动参数，其中，所述运动参数包括位移、速度或者应变；
控制所述处理器基于所述剪切波在不同时刻及不同深度下的运动参数确定所述受测组织的运动参数图像。
3.如权利要求1或者2所述的弹性成像方法，其特征在于，所述控制所述处理器确定所述运动参数图像中所述剪切波的主传播路径包括：
控制所述处理器确定每一深度下位于预设范围的运动参数所对应的目标区域；
控制所述处理器基于所述运动参数图像中相连续的目标区域确定一个或多个带状区域；
控制所述处理器确定所述一个或多个带状区域中满足预设条件的目标带状区域为所述主传播路径。
4.如权利要求1至3中任意一项所述的弹性成像方法，其特征在于，所述弹性成像方法还包括：
控制所述探头在接收到所述发射电路的激励时序时振动，以在所述受测组织内产生所述剪切波。
5.如权利要求3所述的弹性成像方法，其特征在于，每一目标区域对应目标时刻范围，每一目标时刻范围包括若干目标时刻，所述弹性成像方法还包括：
控制所述处理器获取所述探头振动结束时所对应的参考时刻；
其中，所述控制所述处理器确定所述一个或多个带状区域中满足预设条件的目标带状区域为所述主传播路径包括：
控制所述处理器确定所述一个或多个带状区域中目标时刻晚于所述参考时刻的目标区域；
控制所述处理器确定所述一个或多个带状区域中目标时刻晚于所述参考时刻的目标区域所组成的带状区域为所述主传播路径。
6.如权利要求3所述的弹性成像方法，其特征在于，所述控制所述处理器确定所述一个或多个带状区域中满足预设条件的目标带状区域为所述主传播路径包括：
控制所述处理器确定位于预设深度下的目标带状区域为所述主传播路径。
7.如权利要求3所述的弹性成像方法，其特征在于，所述控制所述处理器确定所述一个或多个带状区域中满足预设条件的目标带状区域为所述主传播路径包括：
控制所述处理器确定所述一个或多个带状区域中具有最大长度或者最大面积的目标带状区域为所述主传播路径。
8.如权利要求1所述的弹性成像方法，其特征在于，所述控制所述处理器确定所述运动参数图像中所述剪切波的主传播路径包括：
控制所述处理器对所述运动参数图像进行二值化处理，其中，二值化处理后的所述运动参数图像为所述主传播路径图。
9.如权利要求8所述的弹性成像方法，其特征在于，所述运动参数图像包括若干像素点，每一像素点的像素值与像素点所对应的深度下的运动参数的大小相对应，其中所述控制所述处理器对所述运动参数图像进行二值化处理包括：
控制所述处理器将所述运动参数图像中像素值大于预设阈值的像素点设置为第一颜色；
控制所述处理器将所述运动参数图像中像素值不大于所述预设阈值的像素点设置为第二颜色。
10.如权利要求8所述的弹性成像方法，其特征在于，所述控制所述处理器对所述运动参数图像进行二值化处理包括：
控制所述处理器将所述运动参数图像中所述主传播路径所对应的目标带状区域设置为第一颜色；
控制所述处理器将所述主传播路径所对应的目标带状区域之外的区域设置为第二颜色。
11.如权利要求3所述的弹性成像方法，其特征在于，所述控制所述处理器将所述主传播路径图显示于所述显示屏内，包括：
控制所述处理器将所述主传播路径所对应的目标带状区域中每一深度处的像素点的像素值设置为具有与运动参数对应的第三颜色，其中，不同的运动参数所对应的第三颜色不同；
控制所述处理器将所述主传播路径所对应的目标带状区域之外的区域设置为第四颜色。
12.如权利要求1至11中任意一项所述的弹性成像方法，其特征在于，所述弹性成像方法还包括：
控制所述处理器确定所述受测组织的弹性信息的置信度。
13.如权利要求12所述的弹性成像方法，其特征在于，所述弹性成像方法还包括：
控制所述处理器显示与所述受测组织的弹性信息的置信度相对应的提示信息。
14.如权利要求12或13所述的弹性成像方法，其特征在于，所述控制所述处理器确定所述受测组织的弹性信息的置信度，包括：
控制所述处理器确定所述主传播路径所对应的一个或多个置信度参数；
控制所述处理器基于所述一个或多个置信度参数确定所述受测组织的弹性信息的置信度。
15.如权利要求14所述的弹性成像方法，其特征在于，所述置信度参数包括所述主传播路径的线性度参数、所述主传播路径直线拟合计算所述剪切波传播速度时的误差参数、所述主传播路径的长度参数、所述主传播路径的面积参数中的一个或多个。
16.如权利要求13所述的弹性成像方法，其特征在于，所述控制所述处理器显示与所述受测组织的弹性信息的置信度相对应的提示信息，包括：
控制所述处理器获取所述受测组织的基础图像；
控制所述处理器将所述基础图像及与所述受测组织的弹性信息的置信度相对应的提示信息显示于所述显示屏内。
17.如权利要求16所述的弹性成像方法，其特征在于，所述控制所述处理器获取所述受测组织的基础图像，包括：
控制所述探头在接收到所述发射电路的第二发射时序时向所述受测组织发射第二超声波；
控制所述探头接收所述受测组织返回的第二超声回波，并将所述第二超声回波转换为电信号后传输至所述接收电路；
控制所述波束合成器对所述接收电路传输来的电信号进行波束合成，获得第二超声回波数据；
控制所述处理器基于所述第二超声回波数据生成所述基础图像。
18.如权利要求17所述的弹性成像方法，其特征在于，所述基础图像包括B图像、C图像中的一种或多种。
19.一种弹性成像方法，其特征在于，包括：
向受测组织发射第一超声波，以跟踪在所述受测组织内传播的剪切波；
接收所述受测组织返回的第一超声回波，获得第一超声回波数据；
根据所述第一超声回波数据获得由所述剪切波在所述受测组织中的传播引起的所述受测组织在不同时刻及不同深度下的运动参数；
根据所述运动参数确定所述剪切波的主传播路径，获得主传播路径图，其中所述主传播路径图表征所述剪切波在所述受测组织中的主传播路径；
根据所述第一超声回波数据或者所述运动参数或者所述主传播路径确定所述受测组织的弹性信息；
显示所述主传播路径图和所述受测组织的弹性信息。
20.如权利要求19所述的弹性成像方法，其特征在于，所述运动参数包括位移、速度或者应变。
21.如权利要求19或者20所述的弹性成像方法，其特征在于，所述根据所述运动参数确定所述剪切波的主传播路径包括：
从所述运动参数中确定满足预设条件的目标运动参数；
根据所述目标运动参数对应的深度和时刻获得所述主传播路径。
22.如权利要求21所述的弹性成像方法，其特征在于，所述预设条件为：大于预设阈值。
23.如权利要求21所述的弹性成像方法，其特征在于，所述预设条件为：在预设范围内。
24.如权利要求19至23中任意一项所述的弹性成像方法，其特征在于，还包括：
根据所述主传播路径确定所述受测组织的弹性信息的置信度。
25.如权利要求24所述的弹性成像方法，其特征在于，根据所述主传播路径确定所述受测组织的弹性信息的置信度包括：
根据所述主传播路径确定所述主传播路径所对应的一个或多个置信度参数；
根据所述一个或多个置信度参数确定所述受测组织的弹性信息的置信度。
26.如权利要求25所述的弹性成像方法，其特征在于，所述置信度参数包括所述主传播路径的线性度参数、所述主传播路径直线拟合计算所述剪切波传播速度时的误差参数、所述主传播路径的长度参数、所述主传播路径的面积参数中的一个或多个。
27.如权利要求19至26中任意一项所述的弹性成像方法，其特征在于，还包括：
获取所述受测组织的基础图像；
同时显示所述受测组织的基础图像、所述主传播路径图和所述受测组织的弹性信息。
28.如权利要求27所述的弹性成像方法，其特征在于，所述基础图像包括B图像、C图像中的一种或多种。
29.一种弹性成像方法，其特征在于，包括：
向受测组织发射第一超声波，以跟踪在所述受测组织内传播的剪切波；
接收所述受测组织返回的第一超声回波，获得第一超声回波数据；
根据所述第一超声回波数据获得由所述剪切波在所述受测组织中的传播引起的所述受测组织在不同时刻及不同深度下的运动参数；
根据所述运动参数确定所述剪切波的主传播路径，获得主传播路径图，其中所述主传播路径图表征所述剪切波在所述受测组织中的主传播路径；
显示所述主传播路径图。
30.一种弹性成像方法，其特征在于，所述弹性成像方法包括：
获取受测组织的运动参数图像；
根据所述运动参数图像确定在所述受测组织内传播的剪切波的主传播路径；
显示所述主传播路径。
31.如权利要求30所述的弹性成像方法，其特征在于，所述根据所述运动参数图像确定在所述受测组织内传播的剪切波的主传播路径包括：
获取所述运动参数图像中一个或多个带状区域；
确定所述一个或多个带状区域中满足预设条件的目标带状区域为所述主传播路径。
32.如权利要求31所述的弹性成像方法，其特征在于，所述运动参数图像包括深度及时间的属性信息，所述确定所述一个或多个带状区域中满足预设条件的目标带状区域为所述主传播路径包括：
选择所述一个或多个带状区域中时间的属性信息晚于预设时刻的目标带状区域为所述主传播路径；或者
选择所述一个或多个带状区域中深度的属性信息位于预设深度下的目标带状区域为所述主传播路径；或者
选择所述一个或多个带状区域中具有最大长度或者最大面积的目标带状区域为所述主传播路径。
33.如权利要求30至32中任意一项所述的弹性成像方法，其特征在于，所述主传播路径包括位于所述运动参数图像中的带状区域，所述显示所述主传播路径包括：
将所述运动参数图像中所述主传播路径所对应的带状区域设置为第一颜色；
将所述主传播路径所对应的带状区域之外的区域设置为第二颜色。
34.如权利要求30至32中任意一项所述的弹性成像方法，其特征在于，所述主传播路径包括位于所述运动参数图像中的带状区域，所述显示所述主传播路径包括：
将所述主传播路径所对应的带状区域中每一深度处的像素点的像素值设置为具有与运动参数对应的第三颜色，其中，不同的运动参数所对应的第三颜色不同；
将所述主传播路径所对应的目标带状区域之外的区域设置为第四颜色。
35.如权利要求30至34中任意一项所述的弹性成像方法，其特征在于，所述弹性成像方法还包括：
确定所述受测组织的弹性信息的置信度。
36.如权利要求35所述的弹性成像方法，其特征在于，所述确定所述受测组织的弹性信息的置信度包括：
确定所述主传播路径所对应的一个或多个置信度参数；
基于所述一个或多个置信度参数确定所述受测组织的弹性信息的置信度。
37.如权利要求36所述的弹性成像方法，其特征在于，所述置信度参数包括所述主传播路径的线性度参数、所述主传播路径直线拟合计算所述剪切波传播速度时的误差参数、所述主传播路径的长度参数、所述主传播路径的面积参数中的一个或多个。
38.一种弹性成像系统，其特征在于，所述弹性成像系统包括：
探头，用于向受测组织发射第一超声波，以跟踪在所述受测组织内传播的剪切波，所述探头还用于接收所述受测组织返回的第一超声回波，获得第一超声回波数据；
处理器，连接于所述探头，所述处理器用于基于所述第一超声回波数据获得由所述剪切波在所述受测组织中的传播引起的所述受测组织在不同时刻及不同深度下的运动参数，根据所述运动参数确定所述剪切波的主传播路径，获得主传播路径图，并根据所述第一超声回波数据或者所述运动参数或者所述主传播路径图确定所述受测组织的弹性信息，其中所述主传播路径图表征所述剪切波在所述受测组织中的主传播路径；
显示屏，连接于所处理器，所述处理器用于将所述主传播路径图和所述受测组织的弹性信息显示于所述显示屏内。
39.一种计算机可读存储介质，存储有计算机指令，其特征在于，所述计算机指令被处理器执行时实现如权利要求1至37中任意一项所述的弹性成像方法。
CN201910924688.5A 2019-09-27 2019-09-27 弹性成像方法、系统及计算机可读存储介质 Active CN110710989B (zh)

Priority Applications (5)

Application Number Priority Date Filing Date Title
CN201910924688.5A CN110710989B (zh) 2019-09-27 2019-09-27 弹性成像方法、系统及计算机可读存储介质
CN202111470084.1A CN114287968A (zh) 2019-09-27 2019-09-27 弹性成像方法、系统及计算机可读存储介质
CN202110659607.0A CN113261991A (zh) 2019-09-27 2019-09-27 弹性成像方法、系统及计算机可读存储介质
PCT/CN2019/112507 WO2021056645A1 (zh) 2019-09-27 2019-10-22 弹性成像方法、系统及计算机可读存储介质
CN201980098269.9A CN114072066B (zh) 2019-09-27 2019-10-22 弹性成像方法、系统及计算机可读存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910924688.5A CN110710989B (zh) 2019-09-27 2019-09-27 弹性成像方法、系统及计算机可读存储介质

Related Child Applications (2)

Application Number Title Priority Date Filing Date
CN202111470084.1A Division CN114287968A (zh) 2019-09-27 2019-09-27 弹性成像方法、系统及计算机可读存储介质
CN202110659607.0A Division CN113261991A (zh) 2019-09-27 2019-09-27 弹性成像方法、系统及计算机可读存储介质

Publications (2)

Publication Number Publication Date
CN110710989A true CN110710989A (zh) 2020-01-21
CN110710989B CN110710989B (zh) 2021-12-24

Family

ID=69212008

Family Applications (4)

Application Number Title Priority Date Filing Date
CN201910924688.5A Active CN110710989B (zh) 2019-09-27 2019-09-27 弹性成像方法、系统及计算机可读存储介质
CN202111470084.1A Pending CN114287968A (zh) 2019-09-27 2019-09-27 弹性成像方法、系统及计算机可读存储介质
CN202110659607.0A Pending CN113261991A (zh) 2019-09-27 2019-09-27 弹性成像方法、系统及计算机可读存储介质
CN201980098269.9A Active CN114072066B (zh) 2019-09-27 2019-10-22 弹性成像方法、系统及计算机可读存储介质

Family Applications After (3)

Application Number Title Priority Date Filing Date
CN202111470084.1A Pending CN114287968A (zh) 2019-09-27 2019-09-27 弹性成像方法、系统及计算机可读存储介质
CN202110659607.0A Pending CN113261991A (zh) 2019-09-27 2019-09-27 弹性成像方法、系统及计算机可读存储介质
CN201980098269.9A Active CN114072066B (zh) 2019-09-27 2019-10-22 弹性成像方法、系统及计算机可读存储介质

Country Status (2)

Country Link
CN (4) CN110710989B (zh)
WO (1) WO2021056645A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114072066A (zh) * 2019-09-27 2022-02-18 深圳迈瑞生物医疗电子股份有限公司 弹性成像方法、系统及计算机可读存储介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103300891A (zh) * 2013-05-29 2013-09-18 北京索瑞特医学技术有限公司 对组织的定量弹性信息和结构信息组合显示的方法和系统
US20160018364A1 (en) * 2014-07-16 2016-01-21 Duke University Methods, systems and computer program products for estimating shear wave speed using statistical inference
WO2017223312A1 (en) * 2016-06-22 2017-12-28 Duke University Ultrasound transducers for constructive shear wave interference and related methods and systems
WO2018058113A2 (en) * 2016-09-26 2018-03-29 Duke University Systems and methods for analysis of liquids by covered fluidic channels integrated onto sensor platforms
CN110292399A (zh) * 2018-05-04 2019-10-01 深圳迈瑞生物医疗电子股份有限公司 一种剪切波弹性测量的方法及系统
CN110353731A (zh) * 2019-06-12 2019-10-22 深圳迈瑞生物医疗电子股份有限公司 超声弹性成像方法及系统

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102283679B (zh) * 2011-08-04 2014-05-21 中国科学院深圳先进技术研究院 弹性测量的超声成像系统及测量生物组织弹性的方法
KR20140036650A (ko) * 2012-09-17 2014-03-26 삼성전자주식회사 1차원 초음파 프로브를 이용하여 조직의 탄성을 분석하는 방법, 장치 및 시스템.
US10743837B2 (en) * 2014-08-04 2020-08-18 Delphinus Medical Technologies, Inc. Ultrasound waveform tomography method and system
CN110507359B (zh) * 2014-08-28 2022-06-07 深圳迈瑞生物医疗电子股份有限公司 剪切波成像方法及系统
US11103216B2 (en) * 2014-12-24 2021-08-31 Supersonic Imagine Shear wave elastography method and apparatus for imaging an anisotropic medium
CN104605891B (zh) * 2014-12-31 2017-05-31 中国科学院苏州生物医学工程技术研究所 检测剪切波在生物组织中传播速度的方法、检测生物组织弹性的方法及生物组织弹性成像方法
CN108135568B (zh) * 2015-10-08 2021-11-26 梅约医学教育与研究基金会 用于利用持续换能器振动进行超声弹性成像的系统和方法
CN106618638B (zh) * 2016-11-04 2019-02-26 声泰特（成都）科技有限公司 一种定量剪切波弹性成像系统
KR20180054360A (ko) * 2016-11-15 2018-05-24 삼성메디슨 주식회사 초음파 진단 장치 및 초음파 진단 장치 제어 방법
CN114360727A (zh) * 2017-04-21 2022-04-15 深圳迈瑞生物医疗电子股份有限公司 超声弹性成像装置及弹性成像结果评价方法
CN108720869B (zh) * 2017-04-25 2022-05-27 深圳迈瑞生物医疗电子股份有限公司 超声弹性测量方法及装置
CN107510474B (zh) * 2017-09-21 2020-07-10 深圳开立生物医疗科技股份有限公司 剪切波弹性成像方法及系统
CN109919918A (zh) * 2019-02-21 2019-06-21 清华大学 弹性成像的控制方法及装置、计算机设备及可读存储介质
CN109893172B (zh) * 2019-02-22 2020-06-19 清华大学 基于弹性成像的力学参数的确定方法及装置、计算机设备
CN110710989B (zh) * 2019-09-27 2021-12-24 深圳迈瑞生物医疗电子股份有限公司 弹性成像方法、系统及计算机可读存储介质

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103300891A (zh) * 2013-05-29 2013-09-18 北京索瑞特医学技术有限公司 对组织的定量弹性信息和结构信息组合显示的方法和系统
US20160018364A1 (en) * 2014-07-16 2016-01-21 Duke University Methods, systems and computer program products for estimating shear wave speed using statistical inference
WO2017223312A1 (en) * 2016-06-22 2017-12-28 Duke University Ultrasound transducers for constructive shear wave interference and related methods and systems
WO2018058113A2 (en) * 2016-09-26 2018-03-29 Duke University Systems and methods for analysis of liquids by covered fluidic channels integrated onto sensor platforms
CN110292399A (zh) * 2018-05-04 2019-10-01 深圳迈瑞生物医疗电子股份有限公司 一种剪切波弹性测量的方法及系统
CN110353731A (zh) * 2019-06-12 2019-10-22 深圳迈瑞生物医疗电子股份有限公司 超声弹性成像方法及系统

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114072066A (zh) * 2019-09-27 2022-02-18 深圳迈瑞生物医疗电子股份有限公司 弹性成像方法、系统及计算机可读存储介质
CN114072066B (zh) * 2019-09-27 2024-06-07 深圳迈瑞生物医疗电子股份有限公司 弹性成像方法、系统及计算机可读存储介质

Also Published As

Publication number Publication date
CN110710989B (zh) 2021-12-24
WO2021056645A1 (zh) 2021-04-01
CN113261991A (zh) 2021-08-17
CN114072066A (zh) 2022-02-18
CN114287968A (zh) 2022-04-08
CN114072066B (zh) 2024-06-07

Similar Documents

Publication Publication Date Title
US11635514B2 (en) Imaging methods and apparatuses for performing shear wave elastography imaging
CN108733857B (zh) 超声弹性成像装置及弹性成像结果评价方法
US20210407084A1 (en) Analyzing apparatus and analyzing method
CN110353731B (zh) 超声弹性成像方法及系统
WO2018023336A1 (zh) 超声弹性测量显示方法及系统
US8475382B2 (en) Ultrasound diagnostic apparatus and method for tracing movement of tissue
CN102415902B (zh) 超声波诊断装置以及超声波图像处理装置
KR101922522B1 (ko) 전단파들을 사용하는 사운드 스피드 이미징
CN109303572B (zh) 超声诊断设备及控制超声诊断设备的方法
KR20130080640A (ko) 초음파 영상 제공 방법 및 초음파 영상 제공 장치
KR101776530B1 (ko) 초음파 화상 표시 장치 및 그 제어 프로그램
US20200337679A1 (en) Ultrasonic signal processing apparatus, ultrasonic diagnostic apparatus, ultrasonic signal processing method, and ultrasonic signal processing program
CN117224161A (zh) 一种超声成像方法以及超声成像系统
CN110710989B (zh) 弹性成像方法、系统及计算机可读存储介质
CN113967036B (zh) 超声成像设备及利用超声检测血管壁剪切指数的方法
CN112043307A (zh) 超声波诊断装置及其控制方法、计算机可读取的记录介质
JP5159326B2 (ja) 超音波診断装置
US10492762B2 (en) Ultrasound diagnostic device, ultrasound diagnostic method, and program
KR101107392B1 (ko) 가이드 정보를 제공하는 초음파 시스템 및 방법
WO2020019254A1 (zh) 一种剪切波成像方法及系统
US20240245383A1 (en) Shear wave elasticity imaging method and ultrasound imaging device
JP4392091B2 (ja) 超音波診断装置
CN116528771A (zh) 弹性成像方法和超声成像系统
CN116671960A (zh) 弹性成像方法、装置、电子设备及弹性成像系统
CN117159020A (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
GR01 Patent grant
GR01 Patent grant
EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20200121

Assignee: Shenzhen Mindray Animal Medical Technology Co.,Ltd.

Assignor: SHENZHEN MINDRAY BIO-MEDICAL ELECTRONICS Co.,Ltd.

Contract record no.: X2022440020009

Denomination of invention: Elastography method, system and computer readable storage medium

Granted publication date: 20211224

License type: Common License

Record date: 20220804