[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN118000784A - 剪切波弹性成像测量显示方法及系统 - Google Patents

剪切波弹性成像测量显示方法及系统 Download PDF

Info

Publication number
CN118000784A
CN118000784A CN202410260898.XA CN202410260898A CN118000784A CN 118000784 A CN118000784 A CN 118000784A CN 202410260898 A CN202410260898 A CN 202410260898A CN 118000784 A CN118000784 A CN 118000784A
Authority
CN
China
Prior art keywords
elastic
graphic
interest
group
abscissa
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
CN202410260898.XA
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.)
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 CN202410260898.XA priority Critical patent/CN118000784A/zh
Publication of CN118000784A publication Critical patent/CN118000784A/zh
Pending legal-status Critical Current

Links

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/46—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
    • A61B8/461—Displaying means of special interest
    • A61B8/463—Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
    • 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
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/13—Tomography
    • A61B8/14—Echo-tomography
    • 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/461—Displaying means of special interest
    • 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
    • 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
    • A61B8/5223—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/14—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object using acoustic emission techniques
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/20—Drawing from basic elements, e.g. lines or circles
    • G06T11/206—Drawing of charts or graphs
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2200/00—Indexing scheme for image data processing or generation, in general
    • G06T2200/24—Indexing scheme for image data processing or generation, in general involving graphical user interfaces [GUIs]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2210/00—Indexing scheme for image generation or computer graphics
    • G06T2210/41—Medical

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Surgery (AREA)
  • Public Health (AREA)
  • Biophysics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Veterinary Medicine (AREA)
  • Animal Behavior & Ethology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Acoustics & Sound (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Immunology (AREA)
  • Theoretical Computer Science (AREA)
  • Physiology (AREA)
  • Vascular Medicine (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本发明公开了一种剪切波弹性成像测量显示方法及系统，该系统包括：探头用于在目标组织内部产生剪切波使之行经感兴趣区域，对所述剪切波行经的所述感兴趣区域发送追踪脉冲，接收电路及波束合成模块用于接收所述追踪脉冲的回波信息；弹性相关结果计算模块用于根据所述追踪脉冲的回波信息进行剪切波参数的计算，获得所述感兴趣区域对应的弹性测量结果，记录所述弹性测量结果与统计变量之间的对应关系，产生图形标记，所述图形标记的属性由所述弹性测量结果确定。实施该方法和系统，可以为用户提供更加直观的、可参考性更强的弹性测量统计数据。

Description

剪切波弹性成像测量显示方法及系统
技术领域
本发明涉及超声弹性成像的技术领域，特别是涉及利用剪切波进行超声弹性成像的技术、以及其测量显示方法和系统。
背景技术
超声弹性成像是近年来临床研究关心的热点之一，其主要反映组织的弹性或软硬程度，在组织癌症病变的辅助检测、良恶性判别、预后恢复评价等方面得到越来越多应用。
超声弹性成像主要通过对感兴趣区域内的弹性相关参数进行成像，从而反映组织的软硬程度。近二十年来，已经出现了许多种不同的弹性成像方法，比如基于探头按压组织造成应变的准静态弹性成像，基于声辐射力产生剪切波的弹性成像或弹性测量，基于外部振动产生剪切波的瞬时弹性成像等。其中，剪切波弹性成像或测量技术，由于可以提供定量弹性指标，受到医生的广泛关注和欢迎。
但是，目前的剪切波弹性成像技术，由于容易受到病人呼吸、心跳、血管搏动、探头移动等干扰的影响，测量结果重复性和稳定性相对欠佳，给医生诊断带来一定的困扰。所以，有待对剪切波弹性成像技术的测量显示进行改进。
发明内容
基于前述提到的问题，本发明提供一种新的有关超声弹性测量显示方法，其可以为医生提供弹性测量统计分析数据作为诊断参考，弹性测量结果的显示更为直观的、可参考性更强。
本发明第一方面提供了一种弹性超声测量显示方法，其中，包括：
在目标组织内部产生剪切波使之行经感兴趣区域；
对所述剪切波行经的所述感兴趣区域发送追踪脉冲，接收所述追踪脉冲的回波信息；
根据所述追踪脉冲的回波信息进行剪切波参数的计算，获得所述感兴趣区域对应的弹性测量结果；
记录所述弹性测量结果与统计变量之间的对应关系；
产生图形标记，所述图形标记的属性由所述弹性测量结果确定；
在显示器上显示所述图形标记与所述统计变量之间的关系，形成与所述感兴趣区域关联的弹性统计提示图。
本发明第二方面提供了一种弹性超声测量显示方法，其中，包括：
在目标组织内部产生第一剪切波；
对所述第一剪切波行经的感兴趣区域发送至少一次追踪脉冲，接收所述追踪脉冲的回波信息，获得第一组回波信息；
根据所述第一组回波信息计算剪切波参数，获得第一组弹性测量结果；
产生第一组图形标记，所述第一组图形标记的属性由所述第一组弹性测量结果确定；
在弹性统计提示图窗内显示所述第一组图形标记。
本发明第三方面提供了一种弹性超声测量显示系统，其包括：
探头，用于在目标组织内部产生剪切波使之行经感兴趣区域，对所述剪切波行经的所述感兴趣区域发送追踪脉冲，
接收电路及波束合成模块，用于接收所述追踪脉冲的回波信息；
弹性相关结果计算模块，用于根据所述追踪脉冲的回波信息进行剪切波参数的计算，获得所述感兴趣区域对应的弹性测量结果，记录所述弹性测量结果与统计变量之间的对应关系，产生图形标记，所述图形标记的属性由所述弹性测量结果确定；
显示器，用于显示所述图形标记与所述统计变量之间的关系，形成与所述感兴趣区域关联的弹性统计提示图。
本发明第三方面提供了一种弹性超声测量显示系统，其包括：
探头，用于在目标组织内部产生第一剪切波，对所述第一剪切波行经的感兴趣区域发送至少一次追踪脉冲；
接收电路及波束合成模块，用于接收所述追踪脉冲的回波信息，获得第一组回波信息；
弹性相关结果计算模块，用于根据所述第一组回波信息计算剪切波参数，获得第一组弹性测量结果，产生第一组图形标记，所述第一组图形标记的属性由所述第一组弹性测量结果确定；
显示器，用于在弹性统计提示图窗内显示所述第一组图形标记。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明提供的一种超声成像设备的结构示意图。
图2为本发明提供的一种超声弹性测量显示方法的第一实施例流程示意图；
图3为图2方法执行完成后的一个实施例的显示效果图；
图4为剪切波传播时追踪脉冲的发射接收示意图；
图5为本发明提供的一种超声弹性测量显示方法的一个实施例流程示意图；
图6为在图5所示实施例之后，提供的另一种实施例的方法流程示意图；
图7提供了针对一个感兴趣区域时的一种显示效果示意图；
图8提供了针对二个感兴趣区域时的一种显示效果示意图；
图9提供了一种带悬浮窗的显示效果示意图；
图10提供了针对多个感兴趣区域时的一种显示效果示意图；
图11提供了针对多个感兴趣区域时的一种显示效果示意图；
图12提供了针对多个感兴趣区域时在超声图像上叠加图形标记的显示效果示意图；
图13提供了针对多个感兴趣区域时独立的弹性统计提示图窗的一种显示效果示意图；
图14提供了针对多个感兴趣区域时独立的弹性统计提示图窗的另一种显示效果示意图；
图15提供了针对一个感兴趣区域时的一种显示效果示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明基于剪切波超声成像技术，主要通过在组织内部产生剪切波的传播并检测其传播参数并成像，从而反映组织间的硬度差异。此种方法具有较佳的稳定性和可重复操作性。本发明实施例提供一种测量、统计和显示组织硬度的方法和系统，当用户采集和测量组织弹性结果时，使用一个特殊的统计提示图来显示弹性测量结果的统计分析结果，并展现随时间变量和/或位置变化的弹性测量结果的分布统计情况。时间变量表示任意一种因时间变化或随时间变化而引起变化的参数，例如执行剪切波的测量次数，剪切波参数的采集时间，帧数和系统时间等等；而位置变量表示图像像素点或图像区域的坐标位置。本发明实施例的方法，在用户采集弹性结果的过程中，可直观的表现每次测量结果的稳定性和重复性，不仅可帮助用户方便的剔除异常的测量结果，还可以自动对多次测量结果进行进一步统计计算，从而提升最终测量结果的准确性。
图1提供了本发明实施例中的一种超声成像系统的结构示意图，其包括：探头1、发射电路2、发射/接收选择开关3、接收电路4、波束合成模块5、信号处理模块6、组织图像处理模块71、成像波发射控制模块13和显示器8。
成像波发射控制模块13产生成像波发射时序，控制发射电路2将经过延迟聚焦的具有一定幅度和极性的发射脉冲通过发射/接收选择开关3发送到探头1。探头1受发射脉冲的激励，向目标组织(例如，人体或者动物体内的器官、组织、血管等等，图中未示出)发射超声波束，经一定延时后接收从目标组织反射回来的带有检测对象的信息的超声回波，并将此超声回波重新转换为电信号。接收电路4接收探头1转换生成的电信号，获得超声回波信号，并将这些超声回波信号送入波束合成模块5。波束合成模块5对超声回波信号进行聚焦延时、加权和通道求和等处理，然后将超声回波信号送入信号处理模块6进行相关的信号处理。经过信号处理模块6处理的超声回波信号送入组织图像处理模71，由组织图像处理模71根据用户所需成像模式的不同，对信号进行不同的处理，获得不同模式的超声组织图像数据，然后经对数压缩、动态范围调整、数字扫描变换等处理形成不同模式的超声组织图像，并用以在显示器8上显示，不同模式的超声组织图像包括B图像，C图像，D图像等等，或者其他类型的二维超声组织图像或三维超声组织图像。
在组织弹性超声成像过程中，图1中的系统结构还包括：推动脉冲发射控制模块12和弹性相关结果计算模块72。推动脉冲发射控制模块12产生推动波发射时序，控制发射电路2将经过延迟聚焦的具有一定幅度和极性的发射脉冲通过发射/接收选择开关3发送到探头1，用以在目标组织内部产生剪切波。成像波发射控制模块13产生的成像波发射时序，控制发射电路2将发射脉冲通过发射/接收选择开关3发送到探头1，由探头1受发射脉冲的激励，向目标组织中剪切波行经的感兴趣区域发射追踪脉冲，经一定延时后接收从目标组织内感兴趣区域反射回来的追踪脉冲的超声回波信号，并将此超声回波重新转换为电信号。接收电路4接收超声回波信号，并将这些超声回波信号送入波束合成模块5，然后依次经过信号处理模块6后送入弹性相关结果计算模块72，计算感兴趣区域对应的剪切波参数，获得感兴趣区域对应的弹性测量结果，还可以根据弹性测量进行统计记录分析后计算弹性统计结果，并送入到显示器8上与超声组织图像同时显示。
此外，超声成像系统还包括：操作控制模块9，用于接收来自操作用户输入的调节信号，该调节信号包括对推动波发射序列、成像波发射时序等成像参数进行的调整，或者还可以包括对组织图像处理模块的图像、感兴趣区域或弹性相关结果计算模块的计算结果进行的调整。操作控制模块9可以为人机交互接口，例如键盘、滚轮、带触摸功能的显示屏、鼠标、有关手势控制信号的收发模块等等。
成像波发射控制模块13和推动脉冲发射控制模块12可以集成在一个或多个处理器上，由同一个处理模块分时序或按照用户输入命令分布产生成像波发射时序和推动波发射序列。当然，波束合成模块5、信号处理模块6、组织图像处理模块71、弹性相关结果计算模块72、成像波发射控制模块13和推动脉冲发射控制模块12可以分布集成在不同的处理器或电路模块上以实现不同的功能，或者也可以集成在一个或多个处理器上，并通过软件编程来实现多个功能模块的划分。
图1中给出的系统结构基于声辐射力弹性成像原理，此外基于图1中虚线框22所示的振动装置以及虚线箭头，其实提供了另一种组织弹性成像系统结构，在本实施例中，推动脉冲发射控制模块12产生推动波发射时序，提供给振动装置22用以产生机械振动力，并刺激目标组织在其内部产生剪切波，然后参见前述过程，通过成像波发射控制模块13产生成像波发射时序，控制发射电路2使探头1受发射脉冲的激励，向目标组织中剪切波行经的感兴趣区域发射追踪脉冲，接收追踪脉冲的回波信息，依次通过接收电路4、波束合成模块5、信号处理模块6、弹性相关结果计算模块72获得弹性统计结果，并送入到显示器8上与超声组织图像叠加显示。
可见，无论是声辐射力弹性成像系统，还是瞬时弹性成像系统均可以适用本发明图2等相关的弹性超声测量显示方法。
图2提供了本发明实施例中的一种弹性超声测量显示方法的流程示意图。
在步骤S100中，利用探头在目标组织内部产生剪切波，并使之行经感兴趣区域。
在上述步骤S100之前，还可以基于图1所示的系统结构，利用探头向目标组织发射超波束，用以获取目标组织的超声组织图像，其具体过程可参见图1的相关说明。本实施例中的超声组织图像可以为不同模式的超声组织图像，如B图像，C图像，D图像等等，或者其他类型的二维超声组织图像或三维超声组织图像。如图3所示，在显示器8的显示界面81上，超声组织图像显示区811内显示超声组织图像。超声组织图像显示区811可以是显示界面81上的任意一个区域。利用探头向目标组织发送的超声波束可以是平面波、聚焦波、发散波等任何一种超声波束，或两者的组合。
超声组织图像可以包括利用超声成像系统获取的实时超声组织图像，利用超声成像系统获取的一时间段内的超声组织图像视频数据，或者利用超声成像系统获取的某一帧超声组织图像等等。当然，图1所示的系统结构中还包括用于存储这些缓存图像数据和分析处理中间数据的存储设备，例如，硬盘、缓存模块等等。
基于获得的超声组织图像，组织图像处理模块识别在超声组织图像中的感兴趣区域；或者基于系统默认获得感兴趣区域。本文中提到的感兴趣区域可以包括：超声组织图像中的一个像素点、多个离散分布的像素点、多个连续分布的像素点、表征整个探头扫描区域的图像区域和多个离散分布的图像区域等等中的其中之一。如图3所示，感兴趣区域813为一个小的像素区域或图像区域，然后在根据步骤S600获得的弹性统计提示图812显示在超声组织图像显示区811的下方，其中有多组图形标记814构成，表示连续对感兴趣区域813进行多次弹性测量后的统计结果。从图3中可以很直观的获得同一个感兴趣区域的弹性测量结果的整体分析结果，而很容易直观剔除图中虚线框所包含的两次测量误差过大的测量结果，为用户提供更加精确的诊断结果。
感兴趣区域可以基于系统默认，根据操作控制模块的用户输入，和/或基于图像自动分割处理方法来获得。例如，基于系统默认，可以将表征整个探头扫描区域的图像区域或者表征部分探头扫描区域的多个连续分布像素点构成的图像区域作为感兴趣区域。
又例如，根据用户利用操作控制模块在获取的超声组织图像中输入的选择指令，来确定感兴趣区域。这里的操作控制模块可以为人机交互设备，例如鼠标、键盘，或者滚轮等等，或者当显示器为触摸屏时，操作控制模块可以为手势检测模块，用于检测用户在触摸显示屏上给出的选择指令。本实施例中的感兴趣区域可以为超声组织图像中任意被选择的一个像素点、多个离散分布的像素点、多个连续分布的像素点构成的图像区域或者还可以是表征整个探头扫描区域的图像区域。
还比如，感兴趣区域的识别还可以基于系统自动运行的图像自动分割处理方法。例如，在获得的超声心脏图像中进行图像自动分割处理后，识别出诸如动脉根血管壁、心室壁等等组织器官的图像区域或部分图像区域。因此，利用图像自动分割算法获得的组织器官的图像区域或部分图像区域亦可为步骤S100中的感兴趣区域。
更进一步地，感兴趣区域的识别还可以是半自动过程。例如，基于系统自动运行的图像自动分割处理方法获得目标组织器官区域(例如心脏壁、肝部、胃壁、血管壁等等组织结构)，然后接收用户利用操作控制模块在目标组织器官区域中输入的选择指令，用以确定前述感兴趣区域。还例如，基于系统默认，将表征整个探头扫描区域的图像区域或者表征部分探头扫描区域的多个连续分布像素点构成的图像区域作为目标区域，然后接收用户利用操作控制模块在目标区域中输入的选择指令，用以确定前述感兴趣区域。
基于上述识别获得的感兴趣区域，可以用于定位利用剪切波进行组织弹性测量的目标位置，或者还可以用于确定弹性测量结果对应的测量位置。例如，当利用剪切波对表征整个探头扫描区域的图像区域和部分图像区域(包括多个连续分布的像素点构成的区域)进行组织弹性测量时，则感兴趣可以基于二次定位识别来获得，例如在表征整个探头扫描区域的图像区域和部分图像区域内，接收来自用户输入的选择指令，获得相应的感兴趣区域。在步骤S300中，计算剪切波参数可以是基于表征整个探头扫描区域的图像区域和部分图像区域来进行整体计算，而在获取感兴趣区域对应的弹性测量结果时，可以是从整体计算的结果中提取二次定位识别获得的感兴趣区域对应的剪切波参数来计算弹性测量结果。
因此，上述感兴趣区域的识别可以手动识别，也可以采用自动识别，或者是半自动识别。此外，感兴趣区域的个数可以为一个，也可以为多个(本文中的多个表示两个以上，包含两个，本文同)。在本发明的其中一个实施例中，还可以基于多个感兴趣区域来同时获得多个感兴趣区域对应的弹性统计提示图。在本发明的其中一个实施例中，感兴趣区域的个数可以由用户选择确定或者由系统默认确定。当感兴趣区域的个数为多个时，用户可以同时获取多个感兴趣区域对应的组织弹性参数进行对比观测，显示可用于对比观察多个感兴趣区域的弹性统计提示图。
此外，在本发明的其中一个实施例中，感兴趣区域的个数、方位和大小可以由用户来自行选择，例如，用户通过操作控制模块9输入关于感兴趣区域的个数、方位和大小的调整命令，组织图像处理模块获取该调整命令，更新在显示器上显示的感兴趣区域。此外，还可以依据该调整命令，清空已经经过步骤S600显示的图形标记，即弹性统计提示图中的图形标记，视为感兴趣区域进行了调整，之前的数据可以清除，重新输出显示图形标记。
在步骤S200中，利用探头对剪切波行经的感兴趣区域发送追踪脉冲，接收电路及波束合成模块接收追踪脉冲的回波信息。剪切波的产生方式可以参考前文中针对图1的详细说明。而产生剪切波的位置，可以依据系统默认设定，或者依据感兴趣区域的位置来确定。
在本步骤中，可采用各种方法于组织内部产生剪切波，如通过组织外部的外力振动于组织内部产生剪切波，还可通过向组织内部发射声辐射力脉冲(ARFI，acousticradiation force impulse)等方式于组织内部产生剪切波。其中，所述声辐射力脉冲可聚焦，也可以无聚焦，因此，追踪脉冲可以为平面波、聚焦波、发散波等任何一种超声波束，或两者的组合。
可以理解的是，由于发射声辐射力脉冲产生的剪切波本身幅度较小，又由于剪切波会随着传播迅速衰减，因此可通过发射一系列声辐射力脉冲以提高剪切波的强度，或者拓宽剪切波的传播范围，或者改变剪切波波形特征已提高检测灵敏度等，从而避免由于剪切波的衰减影响成像。
更进一步地，在其中一个实施例中，可向同一个测量位置连续发射多次聚焦脉冲，以提高所产生剪切波的强度。还可以，向改变连续发射的聚焦脉冲的纵向(指聚焦发射的方向)及横向(指与聚焦发射垂直的方向)位置，以拓宽剪切波的传播范围，并使得剪切波沿着某特定方向传播。或者，可同时在不同横向位置发射脉冲，以使先后不同时间到达的两个剪切波波形叠加起来，方便检测。
例如，在其中一个实施例中，步骤S200进一步可以包括以下步骤：
步骤S210，预估剪切波于目标组织内的传播速度。
剪切波产生后，开始在组织中传播，随着组织弹性特性的不同，其传播快慢是不同的。为了预估追击剪切波，需要根据目标组织估计一个平均速度并估计可能的速度变化范围cl至ch，这个平均速度及范围可以参考既有的学术测量数据、或者测量经验等，由系统根据情况预先指定。比如假定目标组织内剪切波的平均传播速度约为2m/s，可能的变化范围为1～4m/s等，或者假定其平均传播速度为1m/s，可能的变化范围为0.5～2m/s等。
步骤S220，根据剪切波于目标组织内的传播速度，预估剪切波在各个时刻于目标组织内的位置，获取剪切波在各个时刻的所述剪切波预估位置。
在剪切波产生后的不同时刻tk，假设剪切波起始传播时刻为t0，则可以估测其距波源的位置距离满足下式关系：
假设相邻两次检测时刻的时间间隔为Δt，即满足：Δt＝tk-tk-1，则可以估测两次检测时刻间剪切波的传播范围Δdl～Δdh，满足：
Δdl＝clΔt
Δdh＝chΔt
于是，相邻两次检测时刻之间剪切波检测位置的移动距离小于或等于Δdl，避免剪切波传播太慢时预估定位超前，同时每次检测时剪切波的检测宽度大于或等于Δdh-Δdl＝Δt(ch-cl)，以保证每次预估定位能包含该时刻剪切波所有可能位置。
步骤S230，分别于各个时刻向相应的剪切波预估位置发送追踪脉冲，并接收各个追踪脉冲的回波信息。
如图4所示，从剪切波传播起始开始，系统以间隔Δt发送追踪脉冲以持续检测，每次检测均保持一定的检测横向波束宽度，即同时收回一定宽度的回波信息，回波信息包括上述一定宽度中各个横向位置的信息，横向位置的间隔不能太大以保证一定的横向分辨力。同时，相邻检测之间波束中心保持小于Δdl的移动距离，或者如果Δt较小导致Δdl太小，则也等效于每隔n个检测时刻之间波束中心距离保持小于nΔdl的移动距离。当然，系统可以从任意时刻开始检测或者从距离剪切波波源一定距离开始检测，需先根据平均传播速度估计下在当前位置处或当前时刻下剪切波的可能的位置，等剪切波传到该位置并离开后才开始改变每次检测的中心位置。
由于每次检测必须保持一定的横向波束宽度，而横向线间距不能太大以保证一定的横向分辨力，则可能要求系统具备超宽波束合成能力，即同时收回多个横向位置的回波信息的能力。波束数目比如1～1024，系统根据需要调节，比如4波束、16波束、32波束、64波束、96波束、128波束等。波束越宽，则意味着发射声场的聚焦越弱，声场能量横向分布越均匀越不集中，这也同时会带来波束中每个检测位置的信噪比的降低。为了改善检测质量，对同一个中心位置，可以连续发射多次，每次发射波束的角度不同，然后将不同角度的回波信号进行合成以增加信噪比。角度数目及偏转角度的大小都由系统根据实际需要来调节，比如采用3个角度，偏转-5°、0°、5°等。
在步骤S300中，弹性相关结果计算模块根据追踪脉冲的回波信息进行剪切波参数的计算，获得感兴趣区域对应的弹性测量结果。根据所述追踪脉冲的回波信息可以计算出多种剪切波参数，比如剪切波传播距离、剪切波传播速度、杨氏模量、剪切模量等等各种相关指标参数。而对于感兴趣区域对应的弹性测量结果，可以是剪切波参数计算的直接结果，也可以是感兴趣区域内相关指标参数的统计量，例如感兴趣区域内相关指标参数的最大值、最小值、平均值、标准差、平均方差、等级评分等。
在本步骤中，可通过对各个时刻的追踪脉冲的回波信息进行整合，从而获取剪切波在传播过程中目标组织各个位置处一小段时间内的回波信息，且剪切波刚好在这一小段时间内通过该对应位置。
在本发明其中一个实施例中，上述步骤S300进一步包括以下步骤：
步骤S301，获取参考信息；可以理解的是，所述参考信息可根据需要自行选用。如选用相应位置某一时刻的追踪脉冲的回波信息作为参考信息。也可在剪切波传播前发送参考脉冲，并将所述参考脉冲的回波信息作为参考回波信息。所述参考需要用于与对所述剪切波进行追击的追踪脉冲做互相关比较。
步骤S302，将目标区域内各个位置的不同时刻的追踪脉冲的回波信息与该位置对应的参考信息做互相关比较，获取该位置处不同时刻的质点位移数据。进一步地，可形成该位置处的位移-时间曲线，在这段时间内，剪切波会经历接近、到达并离开该位置的全过程，对应曲线出现波峰。由于预估追击检测，每个横向位置都能得到相应的一小段位移-时间曲线，只是曲线对应的时刻不同，相邻位置对应的时刻可能会有一部分重叠。位移-时间曲线上峰值所在位置即对应着剪切波到达该位置的时刻。
剪切波的传播速度有多种计算方法可用，比如，对同一深度上两个不同横向位置对应的位移-时间曲线，进行互相关比较，可以得到两横向位置之间对应的时间差异，该时间差对应着这两个横向位置之间的剪切波传播时间。横向位置间的距离与传播时间之比即这两个横行位置间的传播速度。
比如，对某个位置，取出剪切波到达该位置时刻附近两个时刻分别对应的各个横向位置的位移数据，形成两个时刻的位移-横向位置曲线，对两曲线进行互相关比较可以得到两个时刻之间的横向位置差异，该位置差对应了这两个时刻之间的剪切波的传播距离。传播距离与两时刻时间差之比即为该位置附近的传播速度。
比如可以直接利用波的传播方程推导出近似计算公式如下：
式中，c表示传播速度，uz可以视为纵向位移数据，也可采用纵向速度数据进行计算，x代表横向坐标，z代表纵向坐标。还可以将上述公式变换到频域进行计算。
在一定条件下，剪切波的传播速度与组织硬度有近似固定的关系：
E＝3ρc2
式中，ρ表示组织密度，E表示组织的杨氏模量值。一定条件下，杨氏模量越大，意味着组织硬度越大。
此外，从各位置的剪切波的传播速度值，还可以进一步计算得到剪切模量、某段固定时间内的传播距离、目标区域内的平均传播速度、杨氏模量等。
本实施例中的弹性测量结果可以采用杨氏模量、剪切模量、剪切波速度和剪切波传播距离等多种度量方式之一，在本实施例的系统中提供用户自由选择弹性测量结果的度量方式，例如，用户通过操作控制模块9输入关于弹性测量结果的度量方式的改变指令，弹性相关结果计算模块获取该改变指令，更新在显示器上显示的图形标记，更新的内容包括图形标记的属性，或图形标记所在的弹性统计提示图坐标系等等。
上述步骤S100至步骤S300的过程中可包括利用探头进行多次剪切波发射的测量过程。在本发明的其中一个实施例中，利用探头在目标组织内部产生多次剪切波，使之分别行经感兴趣区域；利用探头对多次剪切波行经的感兴趣区域，分别发送多次追踪脉冲，接收电路及波束合成模块接收多次追踪脉冲分别对应的回波信息，获得多组回波信息；弹性相关结果计算模块根据所述多组回波信息进行多次剪切波参数的计算，分别获得多次测量所述感兴趣区域对应的多组弹性测量结果，记录多组弹性测量结果与时间变量的对应关系，可以产生多组图形标记。每组弹性测量结果决定一组图形标记的属性，在显示器上显示多组图像标记与时间变量的关系形成弹性统计提示图，用以展示有关感兴趣区域的多次剪切波测量的统计分析结果，如图3所示。
上述步骤S300中的根据追踪脉冲的回波信息进行剪切波参数的计算过程，可以是基于系统默认的包含感兴趣区域的初始范围来进行，例如，在本发明的其中一个实施例中，基于系统默认或用户的首次输入，将表征整个探头扫描区域的图像区域或者表征部分探头扫描区域的多个连续分布像素点构成的图像区域作为目标区域，该目标区域包含感兴趣区域，然后依次执行上述步骤S100至步骤S300中进行剪切波参数计算的过程，如下所示：
首先，利用探头在目标组织内部产生剪切波使之行经包含感兴趣区域的目标区域；
其次，利用探头对所述剪切波行经的所述目标区域发送追踪脉冲，接收电路及波束合成模块接收所述追踪脉冲的回波信息；
然后，弹性相关结果计算模块根据所述追踪脉冲的回波信息进行剪切波参数的计算，获得所述目标区域对应的弹性测量结果；
最后，弹性相关结果计算模块从所述目标区域对应的弹性测量结果中提取获得感兴趣区域对应的弹性测量结果。
当然上述目标区域不限于只是表征整个探头扫描区域的图像区域或者表征部分探头扫描区域的多个连续分布像素点构成的图像区域，还可以是用户在超声组织图像上选择的任意一个图像区域，或者还可以是系统运行图像自动分割处理方法获取的表征某一或多个生物组织结构的图像区域。
当然在从目标区域对应的弹性测量结果中提取获得所述感兴趣区域对应的弹性测量结果之前还包括，基于用户的选择或者系统的自动识别，来从上述超声组织图像中识别感兴趣区域。
当需要显示弹性测量结果的感兴趣区域的个数为多个时，可以采用上述方式，将多个感兴趣区域所在的图像区域作为上述目标区域进行剪切波的测量过程，获得目标区域对应的弹性测量结果，再从中提取获得多个感兴趣区域分别对应的弹性测量结果。或者，还可以对多个感兴趣区域分别进行剪切波的测量过程，依次获得多个感兴趣区域分别对应的弹性测量结果。
在步骤S400中，弹性相关结果计算模块记录所述弹性测量结果与统计变量之间的对应关系。所述统计变量至少包括位置变量和时间变量中的其中之一。
时间变量表示任意一种因时间变化或随时间变化而引起变化的参数，例如执行剪切波的测量次数，剪切波参数的采集时间，帧数和系统时间等等。
位置变量表示在超声组织图像中图像像素点或图像区域的坐标位置。该位置变量可以是固定某一轴坐标，而另一个或两个轴坐标发生变化；还可以是随超声波束或追踪波束的传播方向发生的坐标位置变化；或者还可以是随剪切波的传播方向发生的坐标位置变化，当然，还可以是随超声组织图像中任意选定的方向上发生的坐标位置变化。
当同时针对多个感兴趣区域时，可以记录多个感兴趣区域分别对应的弹性测量结果与位置变量之间的对应关系，用以根据多个感兴趣区域分别对应的弹性测量结果产生多个图形标记，在弹性统计提示图中依据位置变量的空间位置关系，同步显示多个图形标记，用以对比展示多个感兴趣区域的弹性测量结果。本文中提到的空间位置关系可以包括多个感兴趣区域在超声组织图像中的坐标位置空间关系。
当针对一个感兴趣区域时，可以记录该感兴趣区域对应的弹性测量结果与时间变量之间的对应关系，用以产生随时间变量变化的多组图形标记，在弹性统计提示图中依据时间变量的递变顺序，依次显示相应的图形标记，可以逐一展示多次测量同一位置的测量结果，便于医生整体评估，剔除误差较大的测量结果。
当同时针对多个感兴趣区域组合上述两种情况时，还可以记录多个感兴趣区域分别对应的弹性测量结果与位置变量和时间变量之间的对应关系，用以产生随时间变量变化的多组图形标记，每组图形标记中包括：根据多个感兴趣区域在同一时间变量值时分别对应的弹性测量结果产生的多个图形标记。在弹性统计提示图中对比展示，同时测量组织弹性的多个感兴趣区域，便于用户能够连续观测一定区域范围内某些位置的弹性参数的对比变化情况，更清楚了解整个区域范围内弹性结果的分布。在本发明的其中一个实施例中，在存储设备中记录一段时间内采集获得的所有弹性测量结果与统计变量之间的对应关系。在本发明的其中一个实施例中，在存储设备中记录一段时间内采集获得的所有剪切波参数的计算结果与统计变量之间的对应关系，便于弹性相关结果计算模块从目标区域对应的弹性测量结果中提取获得所述感兴趣区域对应的弹性测量结果，形成所记录的感兴趣区域对应的弹性测量结果与统计变量之间的对应关系。
在步骤S500中，弹性相关结果计算模块产生图形标记，该图形标记的属性由上述弹性测量结果确定。基于时间变量的变化依次产生多组图形标记，而基于多个感兴趣区域的坐标位置变化，分别产生多个图形标记，两者可以相互组合呈现，即每一组图形标记中可以包括多个图形标记，具体可参照前文所述。一组图形标记的属性由一个时间变量值对应的弹性测量结果确定；而一个图形标记的属性由一个感兴趣区域对应的弹性测量结果确定；同理，一组图形标记中的一个图形标记的属性，由一个感兴趣区域在一个时间变量值对应的弹性测量结果确定。图形标记的形状可以表现为矩形条、立体柱、圆点、线段、扇形和气泡等等规则曲线构成的形状中的其中之一，还可以表现为不规则曲线构成的形状，例如图14所示，而弹性统计提示图则是利用上述图形标记展现弹性测量结果与统计变量之间的对应关系的统计图表，通常统计图表可以基于三维直角坐标系、二维直角坐标系、或者极坐标系，还可以基于自定义的坐标系来绘制图形标记。
图形标记的属性包括形状大小、在弹性统计提示图中的坐标和色彩等等。本文的形状大小包括面积、或外形尺寸等。在弹性统计提示图中的坐标包括：单个图形标记在弹性统计提示图中的横纵坐标或极坐标。色彩表示用于渲染图形标记使之相互区分的颜色、透明度、对比度等等属性。弹性测量结果的值的大小可以通过图形标记的属性来展示，例如，不同的弹性测量结果的值对应不同的形状大小、在在弹性统计提示图中的坐标和/或色彩。
在本发明的其中一个实施例中，一个二维坐标平面上的彩色矩形棒(bar)来提示当前测量结果，每获得一次组织弹性测量，即在上述坐标平面上当前时刻为之显示一个bar，bar的纵向高度(即在直角坐标系中的纵坐标)表示弹性结果的大小。与图形标记的形状大小、在弹性统计提示图中的坐标和色彩等等关联的数据，决定了在弹性统计提示图中的图形标记的形状大小、在弹性统计提示图中的坐标和色彩属性中其中之一的值，因此，弹性测量结果的值的大小确定了这一数据的大小。
本实施例的系统中提供多种模型的图形标记及相应的弹性统计提示图供用户自主选择，例如，在本发明的其中一个实施例中，接收来自用户的有关图形标记形状的选择指令，根据图形标记形状的选择指令产生相应形状的图形标记，生成相应的弹性统计提示图。当然，在本发明的其中一个实施例中，系统可以默认图形标记为矩形条或立体柱，弹性统计提示图为柱状或条形统计图。在本发明的其他实施例中可以应用统计学中的多种统计图表来展现弹性测量结果与统计变量之间的对应关系。
参见图5所示，具体展示了相邻两组图形标记产生过程的一种实施例。
在步骤510中，利用探头在目标组织内部产生第一剪切波；
在步骤520中，利用探头对所述第一剪切波行经的感兴趣区域发送至少一次追踪脉冲，接收电路与波束合成模块接收所述追踪脉冲的回波信息，获得第一组回波信息；
在步骤530中，弹性相关结果计算模块根据所述第一组回波信息计算剪切波参数，获得第一组弹性测量结果；
在步骤540中，弹性相关结果计算模块产生第一组图形标记，所述第一组图形标记的属性由所述第一组弹性测量结果确定；
在步骤550中，在显示器上，在弹性统计提示图窗内显示所述第一组图形标记。
上述实施例中提供了产生一组图形标记的过程，一组图形标记对应一次剪切波的测量过程。此外，在进行第二次剪切波的测量时，继续下述过程：
如图6所示，在步骤S550之后还包括：
在步骤610中，利用探头在所述目标组织内部产生第二剪切波；
在步骤620中，利用探头对所述第二剪切波行经的所述感兴趣区域发送至少一次追踪脉冲，接收电路及波束合成模块接收所述追踪脉冲的回波信息，获得第二组回波信息；
在步骤630中，弹性相关结果计算模块根据所述第二组回波信息计算剪切波参数，获得第二组弹性测量结果；
在步骤640中，弹性相关结果计算模块产生第二组图形标记，所述第二组图形标记的属性由所述第二组弹性测量结果确定；
在步骤650中，在显示器上，在所述弹性统计提示图窗内显示所述第二组图形标记。本实施例中，第二组图形标记对应第二次剪切波测量时获得的弹性测量结果。
上述各个实施例中的感兴趣区域的获取方式和定义可参见前文相关说明。另外，第一组图形标记和/或第二组图形标记是基于前后两次剪切波测量过程所获得的展示结果，分别与剪切波的测量次数对应，可参见前文中有关图形标记的相关说明，在此不再累述。当然，在图5和图6中的过程中也可包含有关超声组织图像的产生和显示过程，具体也可参见前文的相关说明。此外，上述提到的弹性统计提示图窗，包括在显示界面上划出的用于显示弹性统计提示图的区域，或者在显示界面上浮动的用于显示弹性统计提示图的窗口区域。
在步骤S600中，在显示器上显示上述图形标记与统计变量之间的关系，形成与感兴趣区域关联的弹性统计提示图，在显示器上用于显示部分或全部弹性统计提示图的图像区域都可以称之为弹性统计提示图窗。
在本实施例中的系统中可根据所需要的效果来呈现不同形式的弹性统计提示图。
例如，如图7所示，在显示器8的显示界面上，超声组织图像显示区704内显示超声组织图像703。超声组织图像显示区704可以是显示界面上的任意一个区域。解剖学组织结构701位于超声组织图像上703，这里的超声组织图像可以是二维图像，也可以是三维图像。感兴趣区域702为一个包含多个像素点的区域，其对应的弹性测量结果可以取该区域内的统计量，例如平均值、均方差等。本实施例中弹性统计提示图采用柱状图统计图表，显示在弹性统计提示图窗705上，图7中弹性统计提示图窗705位于超声组织图像显示区704之外。前几次获得的多组图形标记715已标记为黑色的矩形条，并已展示，而下一次获得的一组图形标记716(图7中表示为虚线框)将要在相应的时间变量值的对应处显示。图7中依据时间变量的递变顺序，依次排列显示多组图形标记，一组图形标记的属性由一个时间变量值对应的弹性测量结果确定。当然，也可以缓存多个时间变量值分别对应的弹性测量结果，并一次性产生并显示相应的多组图形标记。这里的时间变量可以取超声组织图像的帧数、系统时间、剪切测量次数等等与时间相关的参数。在图6所示的实施例中，在弹性统计提示图窗705内，获得的第二组图形标记显示在位于第一组图形标记的临近位置，从而实现沿时间变量的递变顺序依次显示第一组图形标记和第二组图形标记。当然，沿时间变量的递变顺序依次显示第一组图形标记和第二组图形标记时，还可以是按照时间变量的递变顺序用第二组图形标记替换第一组图形标记的显示，从而在同一位置实现图形标记的按时间的依次显示。
更进一步地，如图8所示，当存在两个以上的感兴趣区域702a和702b时，则可以设置相应个数的弹性统计提示图窗705，并用于显示与感兴趣区域702a和702b相对应的弹性统计提示图。为了便于用户识别与之对应的感兴趣区域，则可以在弹性统计提示图上设置标识符717a和717b，具有标识符717a的弹性统计提示图对应于感兴趣区域702a，根据感兴趣区域702a的弹性测量结果获得；具有标识符717b的弹性统计提示图对应于感兴趣区域702b，根据感兴趣区域702b的弹性测量结果获得，当然感兴趣区域702a和702b采用与标识符717a和717b相同属性的标识。当然，除了采用标识符，还可以通过色彩来区分，例如，与感兴趣区域相对应的弹性统计提示图或弹性统计提示图窗采用与感兴趣区域相同的区分色彩区分。当然，图8中区分为两个弹性统计提示图窗，其实，可以直接在显示器的显示界面上划分一个独立的区域作为弹性统计提示图窗，用于显示与多个感兴趣区域相对应的所有弹性统计提示图或显示被选择的弹性统计提示图。此外，弹性统计提示图窗除了是独立于超声组织图像显示区域之外的显示区外，还可是悬浮窗。例如如图9所示，当移动光标的实时位置位于或靠近感兴趣区域702b时显示悬浮窗式的弹性统计提示图窗705，并在弹性统计提示图窗705内显示与感兴趣区域702b对应的弹性统计提示图。当然，如图9所示，可以同时在界面上展现悬浮窗式的弹性统计提示图窗705和独立的区域式的弹性统计提示图窗。
上述图7至图9中的各个弹性统计提示图均以时间变量作为横坐标，弹性测量结果的值作为纵坐标构建直角坐标系，形成弹性统计提示图坐标系，然后将各组图形标记显示在此坐标系中构成弹性统计提示图。相应地，当弹性测量结果的度量方式发生改变时，则弹性统计提示图坐标系的纵坐标也将相应改变。弹性统计提示图坐标系除了二维坐标系，还可以采用三维坐标系来展现。具体可参见下述说明。
又例如，在显示器上依据位置变量的空间位置关系，同步显示多个图形标记，一个图形标记的属性由一个感兴趣区域对应的弹性测量结果确定。当然，还可以使图形标记的显示随光标的移动而移动显示。比如，当移动光标的实时位置位于或靠近所述感兴趣区域时，显示当前相应感兴趣区域对应的图形标记(或一个图像标记)。
如图10所示，在显示界面上有两个感兴趣区域702a和702b，在独立于超声组织图像显示区域之外的弹性统计提示图窗705内显示与感兴趣区域702a和702b相对应的弹性统计提示图。本实施例中弹性统计提示图采用可变的矩形条Sa和Sb，当然，还可以采用立体矩形柱，带色彩的矩形条或立体柱等替代。为了区分矩形条与感兴趣区域702a和702b的对应关系，可采用将矩形条Sa和Sb与感兴趣区域702a和702b标记成对应相同的色彩的方式，或者将矩形条Sa和Sb与感兴趣区域702a和702b对应标记相同的文本提示信息，或者还可以依据位置变量的空间位置关系来展示感兴趣区域702a和702b对应的矩形条Sa和Sb，如图10中，在超声组织图像中感兴趣区域702a在702b的左边，则在弹性统计提示图窗705内，矩形条Sa在Sb的左边。或者还如图11所示，在超声组织图像中感兴趣区域702a在702b的下边，则在弹性统计提示图窗705内，矩形条Sa在Sb的下边。图10中依据时间变量获得的多组图形标记在原始显示位置处替换显示，例如基于图6所示的方法流程和图10所示的展示效果，第一组图形标记显示在弹性统计提示图窗705内呈现为具有填充颜色的矩形条Sa和Sb，而第二组图形标记在第一组图形标记的显示位置处替换第一组图形显示标记，也就是说，在弹性统计提示图窗705内显示第二组图形标记时，矩形条Sa和Sb将呈现在虚线框所示的位置。当然，还可以将第一组图形标记和第二组图形标记都显示在同一个弹性统计提示图坐标系中，比如，分区域罗列显示每个感兴趣区域对应的多组图形标记，参见图8和图9所示。或者还可以如图11所示，构建一个以时间变量为横坐标轴、以弹性测量结果C为纵坐标的弹性统计提示图坐标系，然后按照时间变量将获得多组图形标记依次排列显示，而在每一组图像标记显示时，同步显示与多个感兴趣区域位置相关的多个图形标记，参见图11中将与感兴趣区域702a和702b对应的矩形条Sa和Sb叠加显示，更进一步的，在叠加时还可按照感兴趣区域的空间位置关系进行叠加。当弹性测量结果C的度量方式经用户选择发生改变时，则弹性统计提示图坐标系的纵坐标单位也发生改变。
当然，弹性统计提示图窗也可以采用悬浮窗，图11所示的弹性统计提示图窗705可以采用图9所示的悬浮窗。或者，图10中的矩形条Sa和Sb也可以采用悬浮窗式的结构，如图12所示，将超声组织图像作为弹性统计提示图坐标系，然后当移动光标的实时位置位于或靠近感兴趣时，则显示与该感兴趣区域相应的矩形条。在另一个实施例中，如图12所示，还可以直接在超声组织图像上感兴趣区域的图像显示位置处显示矩形条Sa和Sb，这里的矩形条，也可以是外形采用立体柱的图形标记。
还例如，组合上述两种情况，可以在显示器上即依据时间变量的递变顺序显示图形标记，也同时展现位置变量的空间位置关系。比如，当移动光标的实时位置位于或靠近所述感兴趣区域时，依据时间变量的递变顺序、依次显示当前感兴趣区域对应的多组图形标记。还可以是，在显示器上同时依据时间变量的递变顺序依次显示多组图形标记，和依据位置变量的空间位置关系，同步显示每组图形标记中的多个图形标记，一组图形标记中的一个图形标记的属性，由一个感兴趣区域在一个时间变量值对应的弹性测量结果确定。具体可参加图10、图11、图9、图8所示，均可以同时显示多个感兴趣区域的弹性测量结果与时间变量和位置变量之间的对应关系。
从上述实施例中，可见，弹性统计提示图除了位于超声组织图像显示区之外，还可以采用悬浮窗，或者如图12所示直接将图像标记代表的弹性测量结果直接显示在超声组织图像上，例如，在感兴趣区域的图像显示位置处显示与该感兴趣区域对应的图形标记，当存在多个感兴趣区域时，则在多个感兴趣区域的图像显示位置处分别对应显示多个图形标记。此外，除上述方法之外，本发明的其中一个实施例中还提供一种显示方式，如图13所示。在显示器8的显示界面上，超声组织图像显示区704内显示超声组织图像。超声组织图像显示区704可以是显示界面上的任意一个区域。解剖学组织结构701位于超声组织图像上，在解剖学组织结构701内设置有多个感兴趣区域702C1、702C2、702C3、702C4、702C5。在独立于超声组织图像显示区704之外的弹性统计提示图窗705内，以位置变量为坐标轴设置弹性统计提示图坐标系，其横纵坐标为超声组织图像上的像素点的坐标，坐标系所展示的超声组织图像区域的最大范围包含多个感兴趣区域702C1、702C2、702C3、702C4、702C5。然后按照多个感兴趣区域702C1、702C2、702C3、702C4、702C5的空间位置关系，在坐标系中标记与各个感兴趣区域的图像位置关联对应的图形标记(SC1、SC2、SC3、SC4、SC5)，从而通过间接展示的方式，体现了所要观察的超声组织图像区域的弹性组织分布。识别光标F1的移动位置，当光标F1移动到超声组织图像中相应的感性兴趣区域(如702C1)，则在弹性统计提示图窗705内突出标记与感兴趣区域702C1相对应的图形标记SC1(例如图13中相应突出标记为黑色)。以位置变量(例如图13中的多个感兴趣区域在超声组织图像中的坐标位置关系)建立弹性统计提示图坐标系，然后将多个图形标记分别显示在弹性统计提示图坐标系中与相关感兴趣区域图像显示位置对应的位置处，因此实现了依据位置变量的空间位置关系，将感兴趣区域在超声组织图像中的位置和弹性测量结果映射在弹性统计提示图中，同步显示与多个感兴趣区域关联对应的多个图形标记。此外，当存在连续的多次测量时，图13中的图形标记(图中显示的圆圈)，将会随时间变量的递变顺序，同步更新替换为最新测量的弹性测量结果。例如按照图13的显示方式显示第一组图形标记中与多个感兴趣区域对应的多个图形标记，在获得第二组图形标记时，则将弹性统计提示图窗705中的数据同时更新第二组图形标记中与多个感兴趣区域对应的多个图形标记。
上述图13所示实施例的变形方案如图14所示，弹性统计提示图坐标系的展现形式可以不限制为矩形或者严格的坐标系表达方式，还可以是圆形或者环形展示结构。并且，在弹性统计提示图坐标系中显示的图形标记可以不限于是圆点、圆形块、立体柱、矩形条等等结构方式，还可以采用图14所示的不规则曲线围成的区域，当然，图14的分割区也可以采用规则曲线来包围。这些不规则曲线或规则曲线包围形成的区域在空间位置关系上与多个感兴趣区域在超声组织图像中的坐标位置空间关系一一对应，同样的，这些不规则曲线或规则曲线包围形成的区域表示为图形标记SC1’、SC2’、SC3’、SC4’、SC5’,分别对应于多个感兴趣区域702C1、702C2、702C3、702C4、702C5。多个感兴趣区域702C1、702C2、702C3、702C4、702C5对应的弹性测量结果的值决定了不规则曲线或规则曲线包围形成的区域(图形标记SC1’、SC2’、SC3’、SC4’、SC5’)的形状大小和/或色彩等属性参数，以便用户可以更加直观的了解到包含多个感兴趣区域的超声组织图像区域组织弹性分布。当然，这些多个感兴趣区域702C1、702C2、702C3、702C4、702C5，可以是用户自主选择的，或者还可以是，根据系统自动在一定区域内离散分割获得的多个区域，或者多个像素点。更进一步地，在其中一些实施例中，这些不规则曲线或规则曲线包围形成的区域SC1’、SC2’、SC3’、SC4’、SC5’的边界坐标，与多个感兴趣区域702C1、702C2、702C3、702C4、702C5在超声组织图像中所覆盖的区域边界坐标一一对应。
此外，在显示的超声组织图像上识别光标F1的移动位置，当光标的移动位置位于或靠近感兴趣区域时，在弹性统计提示图内突出标记该感兴趣区域对应的图形标记，例如图13中，采用填充为黑色来突出标记感兴趣区域702C1对应的图形标记SC1。在图14中，采用突出标记不规则曲线围成的区域SC1’的边界。
在其中一些实施例中，利用探头向目标组织发射超声波束，获取所述目标组织的超声组织图像，显示超声组织图像并标记所述感兴趣区域。还可以，显示超声组织图像并标记所述剪切波的产生位置。用以向用户展示对应的感兴趣区域的位置，和/或剪切波的产生位置。
在其中一些实施例中，在显示器的显示界面上设置一定的区域用于显示弹性统计提示图，而弹性统计提示图可以位于超声组织图像显示区之外，或者，还可以，当光标的移动位置位于或靠近所述感兴趣区域时，再在光标的移动位置附近显示所述弹性统计提示图，实现弹性统计提示图与光标随动的效果。
在其中一些实施例中，上述步骤S600还包括：在所述弹性统计提示图中突出渲染当前显示的图形标记。例如图5所示的实施例中，还包括：突出渲染第一组图形标记，或者显示关联指示第一组图形标记的突出标记。突出渲染的方式可以采用加深着色、微放大、填充色彩、提高亮度等渲染处理方式。或者还可以如图15所示，相邻两次获得的第一组图形标记719a和第二组图形标记719b，增加突出标记720(图示中表示为三角形)，当显示到第二组图形标记719b时，则在第二组图形标记719b的下方显示三角形突出标记720，以示突出渲染了当前显示的图形标记。
此外，更进一步地，本实施例的系统中还提供用户可以自主的进行数据的回放和二次统计功能。例如，在本发明的其中一个实施例中，在步骤S600之后还包括：
首先，获取用户在所述弹性统计提示图上批量选择的图形标记；然后，提取所述批量选择的图形标记分别对应的弹性测量结果，统计分析后获得统计评估值并显示。这里的统计评估值可以通过对所述批量选择的图形标记分别对应的弹性测量结果进行求平均、取极大值、取极小值和求标准差等统计分析操作中的任意一种统计分析操作或两种以上统计分析操作而获得的，因此，统计评估值统计评估值至少包括平均值、极大值、极小值和标准差等中的任意一个。如图6所示的实施例中，也可以采用这种批量选择进行二次回放统计的功能，批量选择的图形标记包含第一组图形标记和第二组图形标记的部分或全部。当然，批量选择的图形标记不限于是基于时间变量获得的多组图形标记的部分或全部，还是基于位置变量获得的多个图形标记的部分或全部。也就是说，统计评估值可以是基于一时间段内的弹性测量结果的二次统计分析结果，也可以是基于一定超声组织图像区域范围内的弹性测量结果的二次统计分析结果。
此外，在显示器上还提供增加有关弹性测量统计量的文本提示信息。这里的统计量可以包括上述经过二次回放统计功能获得的统计评估值，或者包括上述步骤S300中计算获得的感兴趣区域内相关指标参数的统计量。例如，当在所述弹性统计提示图中突出渲染当前显示的图形标记时，文本提示信息标识对应显示当前显示的图形标记对应的弹性测量结果，例如，当前显示的图形标记对应的感兴趣区域内相关指标参数的统计量。
在操作二次统计功能时，在所述弹性统计提示图中突出渲染用户在所述弹性统计提示图上批量选择的图形标记。当一段时间的采集完成后，系统可根据用户选定的其中一部分结果进行自动的统计计算，并显示最终的计算结果。被选中的bar需要特殊标记出来，方便用户识别，比如令其颜色发生改变等。
针对弹性统计提示图，本实施例的系统中提供用户多种功能操作，例如可以改变弹性统计提示图的坐标，图形标记的属性表征等等。在本发明的其中一个实施例中，在所述弹性统计提示图中获取有关用户变更纵坐标量程的更新指令；基于所述更新指令，更新所述弹性统计提示图的显示。便于用户根据自己的需要来改变弹性统计提示图的展现形式。在本发明的其中一个实施例中，在所述弹性统计提示图中获取有关用户变更图形标记形状属性的更新指令；基于所述更新指令，更新所述弹性统计提示图中图像标记的形状。
此外，在上述弹性统计提示图中还可以叠加彩色编码信息，其中，利用彩色编码等级表征所述弹性测量结果的值的变化。正如前文所述，图形标记也可以通过更改色彩属性来显示弹性测量结果的值的变化。
在本发明的其中一个实施例中，所述显示所述弹性统计提示图的步骤还包括：
计算所述批量选择的图形标记分别对应的弹性测量结果的平均值，或获取关于感兴趣区域的弹性参考值，本实施例中弹性参考值表示的是基于某一解剖学组织结构本身的特性而基于经验获得的正常值；
在所述弹性统计提示图中叠加显示用于表征所述平均值或弹性参考值的曲线(如图15中所示的折线721)。这样可以将实时测量的弹性结果与正常值结果或平均值结果进行对比判断，为用户提供更多的参考意见。
在本发明的其中一个实施例中，在所述弹性统计提示图中，突出渲染高于或低于所述曲线的用于表征弹性测量结果的图形。还比如，如图6所示的实施例中，识别第一组弹性测量结果或第二组弹性测量结果与标准值的大小关系，标准值包括上述弹性参考值、或关于用户批量选择的图形标记分别对应的弹性测量结果的统计评估值；当所述第一组弹性测量结果或第二组弹性测量结果高于或低于所述标准值时，突出标记所述第一组图形标记或第二组图形标记。本文中提到的突出标记或者突出渲染，旨在区分其他未标记的图形标记，因此渲染方式除了前文中提到的方式外，还可以采用其他渲染突出方式。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory，ROM)或随机存储记忆体(Random AccessMemory，RAM)等。上述各个实施例中的方法步骤并不局限于上述各个实施例所陈述的方式，在不违反前后文逻辑的情况下，各个步骤的执行顺序并不严格的按照附图中箭头所指示的方向进行，同样不影响本发明方案的实现。

Claims (39)

1.一种剪切波弹性测量显示方法，其特征在于，包括：
在目标组织内部产生剪切波使之行经感兴趣区域；
对所述剪切波行经的所述感兴趣区域发送追踪脉冲，接收所述追踪脉冲的回波信息；
根据所述追踪脉冲的回波信息进行剪切波参数的计算，获得所述感兴趣区域对应的弹性测量结果，所述弹性测量结果与时间变量之间具有对应关系，所述弹性测量结果至少包括第一组弹性测量结果和第二组弹性测量结果，所述第一组弹性测量结果先于所述第二组弹性测量结果获得；
生成与所述弹性测量结果对应的多组图形标记，所述多组图形标记和所述时间变量之间具有对应关系，其中，一组图形标记的属性由一个时间变量值对应的弹性测量结果确定，所述多组图形标记至少包括与所述第一组弹性测量结果对应的第一组图形标记和与所述第二组弹性测量结果对应的第二组图形标记；
依据时间变量的递变顺序，在显示界面中依次显示所述第一组图形标记和所述第二组图形标记，以形成弹性统计提示图，其中，所述弹性统计提示图以弹性测量结果为纵坐标轴，以时间变量为横坐标轴；
获取所述弹性测量结果的度量方式的改变指令，所述度量方式包括杨氏模量、剪切模量、剪切波速度和剪切波传播距离中的其中之一；
根据所述改变指令更新所述第一组图形标记和/或所述第二组图形标记的属性，或更新所述弹性统计提示图中所述纵坐标轴的纵坐标单位。
2.如权利要求1所述的方法，其特征在于，所述依据时间变量的递变顺序，在显示界面中依次显示所述第一组图形标记和所述第二组图形标记，包括：
依据时间变量的递变顺序，在显示界面中依次排列显示所述第一组图形标记和所述第二组图形标记。
3.如权利要求2所述的方法，其特征在于，所述依据时间变量的递变顺序，在显示界面中依次排列显示所述第一组图形标记和所述第二组图形标记，包括：
确定所述第一组弹性测量结果对应的时间变量值在所述横坐标轴上的第一横坐标，在所述横坐标轴上与第一横坐标对应的位置显示第一组图形标记，确定所述第二组弹性测量结果对应的时间变量值在所述横坐标轴上的第二横坐标，在所述横坐标轴上与第二横坐标对应的位置显示第二组图形标记。
4.如权利要求1所述的方法，其特征在于，所述依据时间变量的递变顺序，在显示界面中依次显示所述第一组图形标记和所述第二组图形标记，包括：
在显示界面中显示第一组图形标记，而后将所述第一组图形标记替换为所述第二组图形标记进行显示。
5.如权利要求1所述的方法，其特征在于，所述感兴趣区域为多个时；
所述第一组弹性测量结果包括同一时间变量值对应的多个第一弹性结果，一个第一弹性结果对应一个感兴趣区域，所述第一组图形标记包括多个第一图形标记，一个第一图形标记对应一个第一弹性结果，且一个第一图形标记的属性由对应的第一弹性结果确定；和/或
所述第二组弹性测量结果包括同一时间变量值对应的多个第二弹性结果，一个第二弹性结果对应一个感兴趣区域，所述第二组图形标记包括多个第二图形标记，一个第二图形标记对应一个第二弹性结果，且一个第二图形标记的属性由对应的第二弹性结果确定。
6.如权利要求5所述的方法，其特征在于，所述依据时间变量的递变顺序，在显示界面中依次显示所述第一组图形标记和所述第二组图形标记，包括：
确定所述第一组弹性测量结果对应的时间变量值在所述横坐标轴上的第一横坐标，在所述横坐标轴上与第一横坐标对应的位置显示多个所述第一图形标记，确定所述第二组弹性测量结果对应的时间变量值在所述横坐标轴上的第二横坐标，在所述横坐标轴上与第二横坐标对应的位置显示多个所述第二图形标记。
7.如权利要求6所述的方法，其特征在于，所述弹性统计提示图包括至少两个所述横坐标轴，一个所述横坐标轴与一个所述感兴趣区域以及一个第一图形标记对应；
所述确定所述第一组弹性测量结果对应的时间变量值在所述横坐标轴上的第一横坐标，在所述横坐标轴上与第一横坐标对应的位置显示多个所述第一图形标记，包括：
对于至少两个所述横坐标轴中的其中一个所述横坐标轴，确定所述第一组弹性测量结果对应的时间变量值在该横坐标轴的第一横坐标，在该横坐标轴上与第一横坐标对应的位置显示对应的所述第一图形标记；
确定所述第二组弹性测量结果对应的时间变量值在所述横坐标轴上的第二横坐标，在所述横坐标轴上与第二横坐标对应的位置显示多个所述第二图形标记，包括：
对于至少两个所述横坐标轴中的另一个所述横坐标轴，确定所述第二组弹性测量结果对应的时间变量值在该横坐标轴的第二横坐标，在该横坐标轴上与第二横坐标对应的位置显示对应的所述第二图形标记。
8.如权利要求6所述的方法，其特征在于，在显示所述第一组图形标记中的多个第一图形标记时，根据多个所述感兴趣区域之间的空间位置关系确定多个所述第一图形标记之间的相对位置；和/或
在显示所述第二组图形标记中的多个第二图形标记时，根据多个所述感兴趣区域之间的空间位置关系确定多个所述第二图形标记之间的相对位置。
9.如权利要求5所述的方法，其特征在于，所述依据时间变量的递变顺序，在显示界面中依次显示所述第一组图形标记和所述第二组图形标记，包括：
在显示界面中同时显示多个所述第一图形标记，而后将多个所述第一图形标记替换为多个所述第二图形标记进行显示，其中，将各所述感兴趣区域对应的第一图形标记替换为该感兴趣区域对应的第二图形标记进行显示。
10.如权利要求1至9中任一项所述的方法，其特征在于，所述方法还包括：
获取所述弹性统计提示图内批量选择的图形标记，所述批量选择的图形标记包含所述第一组图形标记和第二组图形标记的部分或全部；
提取所述批量选择的图形标记分别对应的弹性测量结果，统计分析后获得统计评估值并显示，其中，所述统计评估值至少包括平均值、极大值、极小值和标准差中的任意一个。
11.如权利要求10所述的方法，其特征在于，还包括：显示有关所述统计评估值的文本提示信息。
12.如权利要求10所述的方法，其特征在于，在所述弹性统计提示图中突出渲染所述弹性统计提示图上批量选择的图形标记。
13.如权利要求10所述的方法，其特征在于，还包括；
识别所述第一组弹性测量结果与标准值的大小关系，所述标准值包括弹性参考值或所述统计评估值，当所述第一组弹性测量结果高于或低于所述标准值时，突出标记所述第一组图形标记；和/或
识别所述第二组弹性测量结果与标准值的大小关系，所述标准值包括弹性参考值或所述统计评估值，当所述第二组弹性测量结果高于或低于所述标准值时，突出标记所述第二组图形标记。
14.如权利要求10所述的方法，其特征在于，还包括：
计算所述批量选择的图形标记分别对应的弹性测量结果的平均值，或获取关于感兴趣区域的弹性参考值；
在所述弹性统计提示图中叠加显示用于表征所述平均值或弹性参考值的曲线。
15.如权利要求1至14中任一项所述的方法，其特征在于，还包括：
利用探头向目标组织发射超声波束用以获取超声组织图像；
显示所述超声组织图像；
将所述弹性统计提示图显示在超声组织图像显示区域之外；或者，在所述超声组织图像显示区域中识别光标的位置，当光标的移动位置位于或靠近所述感兴趣区域时显示所述弹性统计提示图；或者，在所述超声组织图像显示区域中识别光标的位置，当光标的移动位置位于或靠近所述感兴趣区域时突出标记该感兴趣区域关联对应的图形标记；或者，将所述第一组图形标记和/或所述第二组图形标记显示在所述超声组织图像中感兴趣区域的图像显示位置处。
16.如权利要求1所述的方法，其特征在于，所述第一组图形标记和/或第二组图形标记的属性包括：图形标记的形状、图形标记在所述弹性统计提示图坐标系中的坐标和图形标记的色彩中的至少一个。
17.如权利要求1所述的方法，其特征在于，所述在目标组织内部产生剪切波使之行经感兴趣区域，对所述剪切波行经的所述感兴趣区域发送追踪脉冲，接收所述追踪脉冲的回波信息，根据所述追踪脉冲的回波信息进行剪切波参数的计算，获得所述感兴趣区域对应的弹性测量结果，包括：
在目标组织内部产生第一剪切波；
对所述第一剪切波行经的感兴趣区域发送至少一次追踪脉冲，接收所述追踪脉冲的回波信息，获得第一组回波信息；
根据所述第一组回波信息计算剪切波参数，获得第一组弹性测量结果；
在所述目标组织内部产生第二剪切波；
对所述第二剪切波行经的所述感兴趣区域发送至少一次追踪脉冲，接收所述追踪脉冲的回波信息，获得第二组回波信息；
根据所述第二组回波信息计算剪切波参数，获得第二组弹性测量结果。
18.如权利要求1至17任一项所述的方法，其特征在于，所述第一组弹性测量结果和/或第二组弹性测量结果采用杨氏模量、剪切模量、剪切波速度和剪切波传播距离中之一来表述。
19.一种剪切波弹性测量显示方法，其特征在于，包括：
在目标组织内部产生剪切波使之行经多个感兴趣区域，多个所述感兴趣区域之间具有表征空间位置关系的位置变量；
对所述剪切波行经的多个所述感兴趣区域发送追踪脉冲，接收所述追踪脉冲的回波信息；
根据所述追踪脉冲的回波信息进行剪切波参数的计算，获得多个感兴趣区域对应的多个弹性测量结果；
生成与多个所述弹性测量结果对应的多个图形标记，其中，一个图形标记的属性由一个感兴趣区域对应的弹性测量结果确定；
根据所述位置变量在显示界面中同步显示多个所述图形标记，形成弹性统计提示图，其中，所述弹性统计提示图中多个图形标记之间的空间位置关系与多个感兴趣区域之间的坐标位置空间关系一一对应。
20.如权利要求19所述的方法，其特征在于，所述在目标组织内部产生剪切波使之行经多个感兴趣区域之前，还包括：利用探头向目标组织发射超声波束，获取所述目标组织的超声组织图像；
所述弹性统计提示图以位置变量为横纵坐标轴设置弹性统计提示图坐标系；
所述根据所述位置变量在显示界面中同步显示多个所述图形标记，形成弹性统计提示图，包括：
对于任意一个图形标记，根据该图形标记对应的感兴趣区域在所述超声组织图像中的位置，确定该图形标记在所述弹性统计提示图坐标系的横纵坐标，在所述横纵坐标对应的位置显示该图形标记。
21.如权利要求19或20所述的方法，其特征在于，所述位置变量随剪切波的传播方向发生坐标位置变化；
或者，所述在目标组织内部产生剪切波使之行经多个感兴趣区域之前，还包括：利用探头向目标组织发射超声波束，获取所述目标组织的超声组织图像；所述位置变量随超声波束或追踪波束的传播方向发生坐标位置变化。
22.一种剪切波弹性测量结果的显示方法，其特征在于，包括：
在显示界面上显示目标组织的超声组织图像；
在所述超声组织图像上确定所述目标组织的感兴趣区域；
获取所述感兴趣区域的弹性测量结果，其中，所述弹性测量结果与时间变量之间具有对应关系，所述弹性测量结果至少包括第一组弹性测量结果和第二组弹性测量结果，所述第一组弹性测量结果先于所述第二组弹性测量结果获得；
生成所述弹性测量结果对应的多组图形标记，所述多组图形标记和所述时间变量之间具有对应关系，其中，一组图形标记的属性由一个时间变量值对应的弹性测量结果确定，所述多组图形标记至少包括与所述第一组弹性测量结果对应的第一组图形标记和与所述第二组弹性测量结果对应的第二组图形标记；
依据时间变量的递变顺序，在所述显示界面上依次显示所述第一组图形标记和所述第二组图形标记，以形成弹性统计提示图，其中，所述弹性统计提示图以弹性测量结果为纵坐标轴，以时间变量为横坐标轴；
获取所述弹性测量结果的度量方式的改变指令，所述度量方式包括杨氏模量、剪切模量、剪切波速度和剪切波传播距离中的其中之一；
根据所述改变指令更新所述第一组图形标记和/或所述第二组图形标记的属性，或更新所述弹性统计提示图中所述纵坐标轴的纵坐标单位。
23.如权利要求22所述的方法，其特征在于，所述依据时间变量的递变顺序，在所述显示界面中依次显示所述第一组图形标记和所述第二组图形标记，包括：
依据时间变量的递变顺序，在所述显示界面中依次排列显示所述第一组图形标记和所述第二组图形标记。
24.如权利要求23所述的方法，其特征在于，所述依据时间变量的递变顺序，在所述显示界面中依次排列显示所述第一组图形标记和所述第二组图形标记，包括：
确定所述第一组弹性测量结果对应的时间变量值在所述横坐标轴上的第一横坐标，在所述横坐标轴上与第一横坐标对应的位置显示第一组图形标记，而后确定所述第二组弹性测量结果对应的时间变量值在所述横坐标轴上的第二横坐标，在所述横坐标轴上与第二横坐标对应的位置显示第二组图形标记。
25.如权利要求22所述的方法，其特征在于，所述依据时间变量的递变顺序，在显示界面中依次显示所述第一组图形标记和所述第二组图形标记，包括：
在所述显示界面中显示第一组图形标记，而后将所述第一组图形标记替换为所述第二组图形标记进行显示。
26.如权利要求22所述的方法，其特征在于，所述感兴趣区域为多个时:
所述第一组弹性测量结果包括同一时间变量值对应的多个第一弹性结果，一个第一弹性结果对应一个感兴趣区域，所述第一组图形标记包括多个第一图形标记，一个第一图形标记对应一个第一弹性结果，且一个第一图形标记的属性由对应的第一弹性结果确定；和/或
所述第二组弹性测量结果包括同一时间变量值对应的多个第二弹性结果，一个第二弹性结果对应一个感兴趣区域，所述第二组图形标记包括多个第二图形标记，一个第二图形标记对应一个第二弹性结果，且一个第二图形标记的属性由对应的第二弹性结果确定。
27.如权利要求26所述的方法，其特征在于，所述依据时间变量的递变顺序，在所述显示界面中依次显示所述第一组图形标记和所述第二组图形标记，包括：
确定所述第一组弹性测量结果对应的时间变量值在所述横坐标轴上的第一横坐标，在所述横坐标轴上与第一横坐标对应的位置显示多个所述第一图形标记，而后确定所述第二组弹性测量结果对应的时间变量值在所述横坐标轴上的第二横坐标，在所述横坐标轴上与第二横坐标对应的位置显示多个所述第二图形标记。
28.如权利要求27所述的方法，其特征在于，所述弹性统计提示图包括至少两个所述横坐标轴，一个所述横坐标轴与一个所述感兴趣区域以及一个第一图形标记对应；
所述确定所述第一组弹性测量结果对应的时间变量值在所述横坐标轴上的第一横坐标，在所述横坐标轴上与第一横坐标对应的位置显示多个所述第一图形标记，包括：
对于至少两个所述横坐标轴中的其中一个所述横坐标轴，确定所述第一组弹性测量结果对应的时间变量值在该横坐标轴的第一横坐标，在该横坐标轴上与第一横坐标对应的位置显示对应的所述第一图形标记；
确定所述第二组弹性测量结果对应的时间变量值在所述横坐标轴上的第二横坐标，在所述横坐标轴上与第二横坐标对应的位置显示多个所述第二图形标记，包括：
对于至少两个所述横坐标轴中的另一个所述横坐标轴，确定所述第二组弹性测量结果对应的时间变量值在该横坐标轴的第二横坐标，在该横坐标轴上与第二横坐标对应的位置显示对应的所述第二图形标记。
29.如权利要求26所述的方法，其特征在于，在显示所述第一组图形标记中的多个第一图形标记时，根据多个所述感兴趣区域之间的空间位置关系确定多个所述第一图形标记之间的相对位置；和/或
在显示所述第二组图形标记中的多个第二图形标记时，根据多个所述感兴趣区域之间的空间位置关系确定多个所述第二图形标记之间的相对位置。
30.如权利要求26所述的方法，其特征在于，所述依据时间变量的递变顺序，在所述显示界面中依次显示所述第一组图形标记和所述第二组图形标记，包括：
在显示界面中同时显示多个所述第一图形标记，而后将多个所述第一图形标记替换为多个所述第二图形标记进行显示，其中，将各所述感兴趣区域对应的第一图形标记替换为该感兴趣区域对应的第二图形标记进行显示。
31.如权利要求22至30任一项所述的方法，其特征在于，所述方法还包括：显示所述感兴趣区域的弹性测量结果；其中，所述感兴趣区域的弹性测量结果包括剪切波参数计算的直接结果或者所述感兴趣区域内相关指标参数的统计量。
32.如权利要求31所述的方法，其特征在于，所述感兴趣区域内相关指标参数的统计量包括感兴趣区域内相关指标参数的最大值、最小值、平均值、标准差和平均方差中的至少一个；所述相关指标参数包括剪切波传播距离、剪切波传播速度、杨氏模量和剪切模量中的至少一个。
33.如权利要求22至32任一项所述的方法，其特征在于，所述方法还包括：
获取所述弹性统计提示图内批量选择的图形标记，所述批量选择的图形标记包含所述第一组图形标记和第二组图形标记的部分或全部；
提取所述批量选择的图形标记分别对应的弹性测量结果，统计分析后获得统计评估值并显示，其中，所述统计评估值至少包括平均值、极大值、极小值和标准差中的任意一个。
34.如权利要求22所述的方法，其特征在于，还包括；
识别所述第一组弹性测量结果与标准值的大小关系，所述标准值包括弹性参考值或所述统计评估值，当所述第一组弹性测量结果高于或低于所述标准值时，突出标记所述第一组图形标记；和/或
识别所述第二组弹性测量结果与标准值的大小关系，所述标准值包括弹性参考值或所述统计评估值，当所述第二组弹性测量结果高于或低于所述标准值时，突出标记所述第二组图形标记。
35.如权利要求22至30中任一项所述的方法，其特征在于，还包括：
将所述弹性统计提示图显示在超声组织图像显示区域之外；
或者，在所述超声组织图像显示区域中识别光标的位置，当光标的移动位置位于或靠近所述感兴趣区域时显示所述弹性统计提示图；
或者，在所述超声组织图像显示区域中识别光标的位置，当光标的移动位置位于或靠近所述感兴趣区域时突出标记该感兴趣区域关联对应的图形标记；
或者，将所述第一组图形标记和/或所述第二组图形标记显示在所述超声组织图像中感兴趣区域的图像显示位置处。
36.如权利要求35所述的方法，其特征在于，所述将所述弹性统计提示图显示在超声组织图像显示区域之外，包括：
在显示界面上划分出独立于超声组织图像显示区域外的区域以显示所述弹性统计提示图；或者
在显示界面上显示浮动窗口区域，所述浮动窗口区域用于显示所述弹性统计提示图。
37.如权利要求22至36任一项所述的方法，其特征在于，所述第一组图形标记和/或第二组图形标记的属性包括：图形标记的形状、图形标记在所述弹性统计提示图坐标系中的坐标和图形标记的色彩中至少一个。
38.如权利要求22至37任一项所述的方法，其特征在于，所述第一组弹性测量结果和/或第二组弹性测量结果采用杨氏模量、剪切模量、剪切波速度和剪切波传播距离中之一来表述。
39.一种弹性超声测量显示系统，其特征在于，包括：
探头，用于在目标组织内部产生剪切波使之行经感兴趣区域，对所述剪切波行经的所述感兴趣区域发送追踪脉冲，接收所述追踪脉冲的回波信息；
处理器，用于执行如权利要求1-38中任一项所述的方法。
CN202410260898.XA 2016-08-01 2016-08-01 剪切波弹性成像测量显示方法及系统 Pending CN118000784A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202410260898.XA CN118000784A (zh) 2016-08-01 2016-08-01 剪切波弹性成像测量显示方法及系统

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
CN202410260898.XA CN118000784A (zh) 2016-08-01 2016-08-01 剪切波弹性成像测量显示方法及系统
PCT/CN2016/092706 WO2018023336A1 (zh) 2016-08-01 2016-08-01 超声弹性测量显示方法及系统
CN201680087715.2A CN109475343B (zh) 2016-08-01 2016-08-01 剪切波弹性成像测量显示方法及系统

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
CN201680087715.2A Division CN109475343B (zh) 2016-08-01 2016-08-01 剪切波弹性成像测量显示方法及系统

Publications (1)

Publication Number Publication Date
CN118000784A true CN118000784A (zh) 2024-05-10

Family

ID=61072576

Family Applications (2)

Application Number Title Priority Date Filing Date
CN202410260898.XA Pending CN118000784A (zh) 2016-08-01 2016-08-01 剪切波弹性成像测量显示方法及系统
CN201680087715.2A Active CN109475343B (zh) 2016-08-01 2016-08-01 剪切波弹性成像测量显示方法及系统

Family Applications After (1)

Application Number Title Priority Date Filing Date
CN201680087715.2A Active CN109475343B (zh) 2016-08-01 2016-08-01 剪切波弹性成像测量显示方法及系统

Country Status (3)

Country Link
US (1) US11331076B2 (zh)
CN (2) CN118000784A (zh)
WO (1) WO2018023336A1 (zh)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014144948A1 (en) * 2013-03-15 2014-09-18 Stuart Micheal D Visible audiovisual annotation of infrared images using a separate wireless mobile device
JP6515615B2 (ja) * 2015-03-19 2019-05-22 セイコーエプソン株式会社 制御装置、電子機器の制御方法及びプログラム
WO2018178379A1 (en) * 2017-03-31 2018-10-04 Koninklijke Philips N.V. System and method for ultrasound shear wave elastography using external mechanical vibrations
JP6996035B2 (ja) * 2017-11-02 2022-01-17 富士フイルムヘルスケア株式会社 超音波診断装置、および、生体組織の物性評価方法
WO2019218141A1 (zh) 2018-05-15 2019-11-21 深圳迈瑞生物医疗电子股份有限公司 一种剪切波弹性测量方法及剪切波弹性成像系统
US20230181154A1 (en) * 2019-01-22 2023-06-15 UltraVision Corporation Method of detection of microcalcifications by ultrasound
CN110044828B (zh) * 2019-04-26 2022-12-16 华南理工大学 一种声光融合的珊瑚活性探测系统及方法
CN111202519B (zh) * 2020-01-17 2023-04-14 首都医科大学宣武医院 一种在体血栓软硬度检测的方法及其系统
WO2021212494A1 (zh) * 2020-04-24 2021-10-28 深圳迈瑞生物医疗电子股份有限公司 粘弹性测量方法和超声测量系统
WO2021217475A1 (zh) * 2020-04-29 2021-11-04 深圳迈瑞生物医疗电子股份有限公司 弹性成像方法、系统和存储介质
CN111671464A (zh) * 2020-06-19 2020-09-18 深圳开立生物医疗科技股份有限公司 一种图像标记方法、装置、设备及可读存储介质
CN114376602B (zh) * 2021-11-17 2024-09-17 深圳迈瑞生物医疗电子股份有限公司 超声设备获取生理参数的方法和超声设备
CN115644921B (zh) * 2022-10-12 2024-08-16 逸超医疗科技(北京)有限公司 一种自动弹性测量方法
CN115876888B (zh) * 2023-03-02 2023-10-24 华中科技大学 一种适用于脉冲强磁场下的超声回波测量装置

Family Cites Families (32)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3697292B2 (ja) * 1995-05-17 2005-09-21 株式会社東芝 超音波診断装置
US20020173721A1 (en) * 1999-08-20 2002-11-21 Novasonics, Inc. User interface for handheld imaging devices
US6371912B1 (en) * 2000-04-05 2002-04-16 Duke University Method and apparatus for the identification and characterization of regions of altered stiffness
US20030028401A1 (en) * 2001-07-17 2003-02-06 Leon Kaufman Customizable lung report generator
JP4233808B2 (ja) * 2002-06-04 2009-03-04 株式会社日立メディコ 超音波診断装置
US6835177B2 (en) * 2002-11-06 2004-12-28 Sonosite, Inc. Ultrasonic blood vessel measurement apparatus and method
US6994673B2 (en) * 2003-01-16 2006-02-07 Ge Ultrasound Israel, Ltd Method and apparatus for quantitative myocardial assessment
JP3932485B2 (ja) * 2003-05-30 2007-06-20 株式会社日立メディコ 超音波診断装置
JP2005218713A (ja) * 2004-02-06 2005-08-18 Aloka Co Ltd 超音波診断装置
EP1744671B1 (en) * 2004-04-26 2013-10-23 David F. Yankelevitz Medical imaging system for accurate measurement evaluation of changes in a target lesion
WO2006013916A1 (ja) * 2004-08-05 2006-02-09 Hitachi Medical Corporation 弾性像表示方法及び超音波診断装置
US7930014B2 (en) * 2005-01-11 2011-04-19 Volcano Corporation Vascular image co-registration
JP4787570B2 (ja) * 2005-08-17 2011-10-05 日立アロカメディカル株式会社 超音波診断装置
US20090143676A1 (en) * 2005-10-19 2009-06-04 Takeshi Matsumura Ultrasonograph for Creating Elastic Image
JP5028416B2 (ja) * 2006-05-25 2012-09-19 株式会社日立メディコ 超音波診断装置
JP5156206B2 (ja) * 2006-08-28 2013-03-06 株式会社日立メディコ 超音波診断装置
WO2009154133A1 (ja) * 2008-06-16 2009-12-23 株式会社 日立メディコ 超音波診断装置、超音波画像表示方法及び超音波診断プログラム
US8647276B2 (en) * 2009-09-04 2014-02-11 Hitachi Medical Corporation Ultrasonic diagnostic device
CN102667522B (zh) * 2009-11-25 2014-10-08 皇家飞利浦电子股份有限公司 采用聚焦扫描线波束形成的超声剪切波成像
JP4709937B2 (ja) * 2010-10-01 2011-06-29 株式会社東芝 超音波診断装置及び画像処理装置
US20120253194A1 (en) * 2011-03-30 2012-10-04 Tadashi Tamura Methods and apparatus for ultrasound imaging
RU2633915C2 (ru) * 2012-09-27 2017-10-19 Конинклейке Филипс Н.В. Автоматизированная последовательность операций в двухплоскостном режиме pw для ультразвуковой оценки стеноза
WO2014061370A1 (ja) * 2012-10-18 2014-04-24 日立アロカメディカル株式会社 超音波診断装置及び画像表示方法
CN103845074B (zh) * 2012-11-28 2017-12-12 深圳迈瑞生物医疗电子股份有限公司 一种超声弹性成像系统和方法
JP6292836B2 (ja) * 2012-12-28 2018-03-14 キヤノン株式会社 被検体情報取得装置、表示方法、プログラム、処理装置
JP2015013099A (ja) * 2013-06-04 2015-01-22 セイコーエプソン株式会社 超音波測定装置および超音波測定方法
WO2014207605A1 (en) * 2013-06-26 2014-12-31 Koninklijke Philips N.V. System and method for mapping ultrasound shear wave elastography measurements
JP5730979B2 (ja) * 2013-11-08 2015-06-10 日立アロカメディカル株式会社 超音波診断装置、及び弾性評価方法
EP3091907A1 (en) * 2014-01-02 2016-11-16 Koninklijke Philips N.V. Ultrasound navigation/tissue characterization combination
KR102185726B1 (ko) * 2014-01-28 2020-12-02 삼성메디슨 주식회사 관심 영역에 대응하는 초음파 영상 표시 방법 및 이를 위한 초음파 장치
JP6325270B2 (ja) * 2014-02-07 2018-05-16 キヤノンメディカルシステムズ株式会社 超音波診断装置及び医用画像処理プログラム
US10813624B2 (en) * 2015-10-30 2020-10-27 Carestream Health, Inc. Ultrasound display method

Also Published As

Publication number Publication date
US11331076B2 (en) 2022-05-17
CN109475343A (zh) 2019-03-15
US20190254629A1 (en) 2019-08-22
CN109475343B (zh) 2024-04-02
WO2018023336A1 (zh) 2018-02-08

Similar Documents

Publication Publication Date Title
CN109475343B (zh) 剪切波弹性成像测量显示方法及系统
CN108784735B (zh) 用于显示采集质量水平的超声成像系统和方法
US11715202B2 (en) Analyzing apparatus and analyzing method
CN108733857B (zh) 超声弹性成像装置及弹性成像结果评价方法
KR101654674B1 (ko) 탄성 영상 제공 방법 및 이를 위한 초음파 장치
CN100556360C (zh) 超声波探头轨迹显示装置及方法和超声波诊断装置及方法
KR102188176B1 (ko) 침 시술 가이드 기능을 가진 초음파 영상 기기
CN110072466B (zh) 产前超声成像
US11622743B2 (en) Rib blockage delineation in anatomically intelligent echocardiography
KR101579740B1 (ko) 초음파 진단장치, 그에 따른 초음파 진단 방법 및 그에 따른 컴퓨터 판독 가능한 저장매체
JP7285826B2 (ja) 肺超音波検査におけるｂラインの検知、提示及び報告
JP4879263B2 (ja) 超音波診断装置及び超音波診断方法
CN105025807B (zh) 用来对静止流体池进行识别和定界的扫描
JP2023160986A (ja) 超音波診断装置及び解析装置
CN111281425B (zh) 用于显示目标对象质量水平的超声成像系统和方法
CN114072060A (zh) 一种超声成像方法以及超声成像系统
JP2013099386A (ja) 超音波診断装置および医用画像処理装置
CN113616238A (zh) 用于确定声学接触的超声成像系统和方法
CN114072066A (zh) 弹性成像方法、系统及计算机可读存储介质
JP7299100B2 (ja) 超音波診断装置及び超音波画像処理方法
CN118695814A (zh) 胎儿心率监测
KR101882325B1 (ko) 스펙트럴 도플러 영상을 생성하는 방법 및 그 초음파 도플러 영상 장치
CN113040822A (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