CN114903519A - 超声波成像方法、超声成像设备及图像显示装置 - Google Patents

超声波成像方法、超声成像设备及图像显示装置 Download PDF

Info

Publication number
CN114903519A
CN114903519A CN202210409389.XA CN202210409389A CN114903519A CN 114903519 A CN114903519 A CN 114903519A CN 202210409389 A CN202210409389 A CN 202210409389A CN 114903519 A CN114903519 A CN 114903519A
Authority
CN
China
Prior art keywords
microbubble
signal
target
data
frequency signals
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
CN202210409389.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 CN202210409389.XA priority Critical patent/CN114903519A/zh
Publication of CN114903519A publication Critical patent/CN114903519A/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/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
    • 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/0891—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of blood vessels
    • 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/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4483—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer
    • A61B8/4488—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer the transducer being a phased array
    • 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

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Veterinary Medicine (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Vascular Medicine (AREA)
  • Gynecology & Obstetrics (AREA)
  • Image Processing (AREA)

Abstract

本申请实施例公开了一种超声波成像方法、超声成像设备及图像显示装置，通过采用两种处理效果不同的时空滤波处理方法，来对获取到的组织结构信号进行处理，以过滤空间射频信号中的静态组织射频信号，提取微泡射频信号，得到第一微泡信号和第二微泡信号两组微泡信号，并进行融合处理，可以使得成像的图像兼具两种时空滤波处理的优点，能够提高对微血管成像的精确性，提升微血管的成像效果。

Description

超声波成像方法、超声成像设备及图像显示装置
技术领域
本申请实施例涉及超声成像领域，尤其涉及一种超声波成像方法、超声成像设备及图像显示装置。
背景技术
医用超声成像诊断设备利用超声波在人体中的传播，能够得到人体组织和器官结构的超声图像信息。由于超声诊断具有安全、适应面广、直观、可反复检查、对软组织鉴别力强、灵活性强及价廉等优点，因此已成为当代医学图像诊断中的首选技术，并已在现代诊断技术中占有极为重要的地位。其中，利用超声成像诊断设备对微血管、血管微循环进行检查、观察，对于各类疾病的早期诊断及治疗具有重要的价值，但由于超声波在远场面临衍射极限限制，导致临床常规超声造影显示微血管结构细节的能力有限。
为了解决微血管结构细节的显示问题，目前的做法是借鉴光学超分辨成像中的荧光显微定位技术原理，通过引入超声造影剂，定位与跟踪超声造影剂孤立微泡而获得空间分辨率为数十微米的图像，为了能对超声造影剂的微泡进行精确的识别、定位，需要对超声信号进行数据处理，然而，现有的数据处理方法得到的微泡定位结果并不理想，成像效果不理想。
发明内容
以下是对本文详细描述的主题的概述。本概述并非是为了限制权利要求的保护范围。
本申请实施例提供了一种超声波成像方法、超声成像设备及图像显示装置，能够提高对微血管成像的精确性，提升微血管的成像效果。
第一方面，本申请实施例提供了一种超声波成像方法，包括：
向目标组织发射第一超声波，接收所述目标组织返回的所述第一超声波的回波，以获得第一超声回波信号；
根据所述第一超声回波信号，得到所述目标组织的组织结构信号，所述组织结构信号包括多帧空间射频信号；
对所述组织结构信号进行奇异值分解滤波处理，以提取各帧所述空间射频信号中的微泡射频信号，得到第一微泡信号；
对所述组织结构信号进行时域差分滤波处理，以提取各帧所述空间射频信号中的微泡射频信号，得到第二微泡信号；
通过对所述第一微泡信号和所述第二微泡信号进行融合处理，得到目标微血管图像数据；
根据所述目标微血管图像数据显示微血管图像。
第二方面，本申请实施例提供了一种超声波成像方法，包括：
向目标组织发射第一超声波，接收所述目标组织返回的所述第一超声波的回波，以获得第一超声回波信号；
根据所述第一超声回波信号，得到所述目标组织的组织结构信号，所述组织结构信号包括多帧空间射频信号；
对所述组织结构信号进行第一微泡信号处理，以提取各帧所述空间射频信号中的微泡射频信号，得到第一微泡信号；
对所述组织结构信号进行第二微泡信号处理，以提取各帧所述空间射频信号中的微泡射频信号，得到第二微泡信号，其中，所述第一微泡信号处理和所述第二微泡信号处理的处理效果不同，所述第一微泡信号处理和所述第二微泡信号处理均为时空滤波处理，所述时空滤波处理为根据多帧所述空间射频信号之间数据的比较来过滤所述空间射频信号中的静态组织射频信号，以在所述空间射频信号中提取所述微泡射频信号；
通过对所述第一微泡信号和所述第二微泡信号进行融合处理，得到目标微血管图像数据；
根据所述目标微血管图像数据显示微血管图像。
第三方面，本申请实施例提供了一种超声波成像方法，包括：
向目标组织发射第一超声波，接收所述目标组织返回的所述第一超声波的回波，以获得第一超声回波信号；
根据所述第一超声回波信号，得到所述目标组织的组织结构信号，所述组织结构信号包括多帧空间射频信号；
对所述组织结构信号进行多种时空滤波处理，提取微泡射频信号，得到多组微泡信号，其中，所述时空滤波处理为根据多帧所述空间射频信号之间数据的比较来过滤所述空间射频信号中的静态组织射频信号，在所述空间射频信号中提取所述微泡射频信号，进而得到所述微泡信号；
对各组所述微泡信号进行融合处理，得到目标微血管图像数据；
根据所述目标微血管图像数据显示微血管图像。
第四方面，本申请实施例提供了一种超声波成像方法，包括：
获取组织结构信号，所述组织结构信号包括多帧空间射频信号；
对所述组织结构信号进行第一微泡信号处理，以提取各帧所述空间射频信号中的微泡射频信号，得到第一微泡信号；
对所述组织结构信号进行第二微泡信号处理，以提取各帧所述空间射频信号中的微泡射频信号，得到第二微泡信号，其中，所述第一微泡信号处理和所述第二微泡信号处理的处理效果不同，所述第一微泡信号处理和所述第二微泡信号处理均为时空滤波处理，所述时空滤波处理为根据多帧所述空间射频信号之间数据的比较来过滤所述空间射频信号中的静态组织射频信号，在所述空间射频信号中提取所述微泡射频信号；
通过对所述第一微泡信号和所述第二微泡信号进行融合处理，得到目标微血管图像数据；
根据所述目标微血管图像数据显示微血管图像。
第五方面，本申请实施例提供了一种超声波成像方法，包括：
获取组织结构信号，所述组织结构信号包括多帧空间射频信号；
对所述组织结构信号进行奇异值分解滤波处理，以提取各帧所述空间射频信号中的微泡射频信号，得到第一微泡信号；
对所述组织结构信号进行时域差分滤波处理，以提取各帧所述空间射频信号中的微泡射频信号，得到第二微泡信号；
通过对所述第一微泡信号和所述第二微泡信号进行融合处理，得到目标微血管图像数据；
根据所述目标微血管图像数据显示微血管图像。
第六方面，本申请实施例提供了一种超声波成像方法，包括：
获取组织结构信号，所述组织结构信号包括多帧空间射频信号；
对所述组织结构信号进行多种时空滤波处理，提取微泡射频信号，得到多组微泡信号，其中，所述时空滤波处理为根据多帧所述空间射频信号之间数据的比较来过滤所述空间射频信号中的静态组织射频信号，在所述空间射频信号中提取所述微泡射频信号，进而得到所述微泡信号；
对各组所述微泡信号进行融合处理，得到目标微血管图像数据；
根据所述目标微血管图像数据显示微血管图像。
第七方面，本申请实施例提供了一种超声成像装置，包括：
超声探头；
发射/接收电路，所述发射/接收电路用于控制所述超声探头向目标组织发射超声波并接收超声回波，获得超声回波信号；
处理器，所述处理器用于处理所述超声回波信号，获得所述目标组织的组织结构信号以及通过对所述组织结构信号进行微泡信号处理，获得微泡信号并根据所述微泡信号得到目标微血管图像数据；
显示器，所述显示器用于显示所述微血管图像数据；
所述处理器还用于执行上述第一方面、第二方面或第三方面实施例的超声波成像方法。
第八方面，本申请实施例提供了一种图像显示装置，包括：
处理器，所述处理器用于获取并处理组织结构信号，其中所述组织结构信号包括多帧空间射频信号；
显示器，所述显示器用于显示所述微血管图像数据；
所述处理器还用于执行上述第四方面、第五方面或第六方面实施例的超声波成像方法。
本申请实施例提供的超声波成像方法、超声成像设备及图像显示装置，通过采用两种处理效果不同的时空滤波处理方法，来对获取到的组织结构信号进行处理，以过滤空间射频信号中的静态组织射频信号，提取微泡射频信号，得到第一微泡信号和第二微泡信号两组微泡信号，并进行融合处理，可以兼具两种时空滤波处理的优点，能够提高对微血管成像的精确性，提升微血管的成像效果。
附图说明
附图用来提供对本申请技术方案的进一步理解，并且构成说明书的一部分，与本申请的实施例一起用于解释本申请的技术方案，并不构成对本申请技术方案的限制。
图1是本申请一个实施例提供的超声成像装置的结构示意图；
图2是本申请一个实施例提供的超声波成像方法的流程图；
图3是本申请一个实施例提供的奇异值分解滤波处理的流程图；
图4是图2中的步骤S250所包括的方法步骤的流程图；
图5是图4中的步骤S430所包括的方法步骤的流程图；
图6是图5中的步骤S520所包括的方法步骤的流程图；
图7是另一实施例提供的图2中的步骤S250所包括的方法步骤的流程图；
图8是又一实施例提供的图2中的步骤S250所包括的方法步骤的流程图；
图9是本申请另一个实施例提供的超声波成像方法的流程图；
图10是本申请又一个实施例提供的超声波成像方法的流程图；
图11是本申请一个实施例提供的图像显示装置的结构示意图；
图12是本申请一个实施例提供的超声波成像方法的流程图；
图13是本申请另一个实施例提供的超声波成像方法的流程图；
图14是本申请又一个实施例提供的超声波成像方法的流程图；
图15是本申请一个实施例提供的超声波成像效果与SVD滤波成像效果、DI处理成像效果的对比示意图。
具体实施方式
下面结合说明书附图和具体的实施例对本申请进行进一步的说明。所描述的实施例不应视为对本申请的限制，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
在以下的描述中，涉及到“一些实施例”，其描述了所有可能实施例的子集，但是可以理解,“一些实施例”可以是所有可能实施例的相同子集或不同子集，并且可以在不冲突的情况下相互结合。
除非另有定义，本文所使用的所有的技术和科学术语与属于本申请的技术领域的技术人员通常理解的含义相同。本文中所使用的术语只是为了描述本申请实施例的目的，不是旨在限制本申请。
医用超声成像诊断设备利用超声波在人体中的传播，能够得到人体组织和器官结构的超声图像信息。由于超声诊断具有安全、适应面广、直观、可反复检查、对软组织鉴别力强、灵活性强及价廉等优点，因此已成为当代医学图像诊断中的首选技术，并已在现代诊断技术中占有极为重要的地位。其中，利用超声成像诊断设备对微血管、血管微循环进行检查、观察，对于各类疾病的早期诊断及治疗具有重要的价值，但由于超声波在远场面临衍射极限限制，导致临床常规超声造影显示微血管结构细节的能力有限。
为了解决微血管结构细节的显示问题，目前的做法借鉴光学超分辨成像中的荧光显微定位技术原理，通过引入超声造影剂，定位与跟踪超声造影剂孤立微泡而获得空间分辨率为数十微米的图像，为了能对超声造影剂的微泡进行精确的识别、定位，需要对超声信号进行数据处理，然而，现有的数据处理方法得到的微泡定位结果并不理想，成像效果不理想。
基于此，本发明提出了一种超声波成像方法、超声成像设备及图像显示装置图像处理装置，能够提高对微血管成像的精确性，提升微血管的成像效果。
如图1所示为一种超声成像装置的结构示意图。其中，超声成像装置100包括超声探头110、发射/接收电路120、处理器130以及显示器140；发射/接收电路120用于控制超声探头110向目标组织发射超声波并接收超声回波，获得超声回波信号；处理器130用于处理超声回波信号，获得目标组织的组织结构信号以及通过对组织结构信号进行微泡提取处理，获得微泡信号并根据微泡信号得到目标微血管图像数据；显示器140用于显示微血管图像数据。
超声探头110包括由阵列式排布的多个阵元组成的换能器(图中未示出)，多个阵元排列成一排构成线阵，或排布成二维矩阵构成面阵，多个阵元也可以构成凸阵列。阵元用于根据激励电信号发射超声波束，或将接收的超声波束变换为电信号。因此每个阵元可用于实现电脉冲信号和超声波束的相互转换，从而实现向人体组织的目标组织(例如心脏、肺部、子宫等)发射超声波、也可用于接收经组织反射回的超声波的回波。参与超声波发射的阵元可以同时被电信号激励，从而同时发射超声波；或者参与超声波发射的阵元也可以被具有一定时间间隔的若干电信号激励，从而持续发射具有一定时间间隔的超声波。
发射电路可以根据处理器130的控制产生发射序列，发射序列用于控制多个阵元中的部分或者全部向目标组织发射超声波，发射序列参数包括发射用的阵元位置、阵元数量和超声波束发射参数(例如幅度、频率、发射次数、发射间隔、发射角度、波型、聚焦位置等)。某些情况下，发射电路可以对发射的波束进行相位延迟，使不同的发射阵元按照不同的时间发射超声波，以便各发射超声波束能够在预定的感兴趣区域聚焦。
接收电路可以从超声探头110接收超声回波的电信号，并对超声回波的电信号进行处理。接收电路可以包括一个或多个放大器、模数转换器(ADC)等。放大器用于在适当增益补偿之后放大所接收到的超声回波的电信号，模数转换器用于对模拟回波信号按预定的时间间隔进行采样，从而转换成数字化的信号，数字化后的回波信号依然保留有幅度信息、频率信息和相位信息。
处理器130的图像处理模块用于对数字化后的回波信号数据进行处理，例如获得目标组织的组织结构信号以及通过对组织结构信号进行微泡信号处理，获得微泡信号并根据微泡信号得到目标微血管图像数据。图像处理模块可以将微血管图像数据输出至人机交互装置的显示器140进行显示。其中人机交互装置用于进行人机交互，即接收用户的输入和输出可视化信息；其接收用户的输入可采用键盘、操作按钮、鼠标、轨迹球等，也可以采用与显示器集成在一起的触控屏；其输出可视化信息采用显示器140。
应理解，图1所示的超声成像装置100所包括的部件只是示意性的，其可以包括更多或更少的部件。本发明对此不限定。图1所示的超声成像装置100中的处理器130还用于执行下文图2至图10中的超声波成像方法或者用于执行下文图12至图14中的超声波成像方法。
参照图2，图2所示为本申请一个实施例的超声波成像方法的流程图。超声波成像方法可以包括如下步骤：
步骤S210：向目标组织发射第一超声波，接收目标组织返回的第一超声波的回波，以获得第一超声回波信号。
可以理解的是，向目标组织发射超声波之前，需要先向目标组织注射超声造影剂，然后通过例如图1中的发射/接收电路120将有一定幅度和极性的发射脉冲发送到超声探头110，以激励超声探头110向目标组织发射超声波。经一定延时后，发射/接收电路120接收超声回波，从而获得超声回波信号。
步骤S220：根据第一超声回波信号，得到目标组织的组织结构信号，组织结构信号包括多帧空间射频信号。
需要说明的是，在本实施例中，目标组织的组织结构信号包括多帧空间射频信号，因此本申请实施例的后续步骤是对多帧空间射频信号进行处理。可以理解的是，发射/接收电路120所获取到的第一超声回波信号也应包括多组信号，从而可以根据多组信号得到多帧空间射频信号。
步骤S230：对组织结构信号进行奇异值分解滤波处理，以提取各帧空间射频信号中的微泡射频信号，得到第一微泡信号。
采用奇异值分解(Singular Value Decomposition，SVD)滤波处理方法对组织结构信号进行处理，也即对多帧空间射频信号进行批处理，可以根据图像矩阵的奇异值提取各帧空间射频信号中的微泡射频信号，使得后续成像的图像的超微血管结构更完整。
步骤S240：对组织结构信号进行时域差分滤波处理，以提取各帧空间射频信号中的微泡射频信号，得到第二微泡信号。
采用时域差分滤波处理方法对组织结构信号进行处理，利用相邻两帧空间射频信号的差异来提取各帧空间射频信号中的微泡射频信号，可以使得成像的图像中上级微血管具有较优的连续性。
步骤S250：通过对第一微泡信号和第二微泡信号进行融合处理，得到目标微血管图像数据。
步骤S260：根据目标微血管图像数据显示微血管图像。
在本实施例中，通过采用奇异值分解滤波处理和时域差分滤波处理两种处理效果不同的时空滤波处理方法，来对组织结构信号进行处理，以过滤空间射频信号中的静态组织射频信号，提取微泡射频信号，得到第一微泡信号和第二微泡信号两组微泡信号，然后进行融合处理，可以使得成像的图像兼具两种时空滤波处理的优点，能够提高对微血管成像的精确性，提升微血管的成像效果。
在本申请的一个实施例中，步骤S240中所提及的时域差分滤波处理，包括以下步骤：
根据两帧预设间距的空间射频信号中的其中一帧空间射频信号与另一帧作差，过滤静态组织射频信号，得到微泡射频信号。其中，作差包括幅度差、相位差或频率差中的至少一种。
可以理解的是，两帧空间射频信号数据之间的信息差异可以包括幅度差异、相位差异或频率差异，将一帧空间射频信号与另一帧作差，计算组织结构信号中两帧空间射频信号的差分数据，然后可以根据差分数据对两帧空间射频信号中的其中一帧空间射频信号进行区域划分，以确定空间射频信号的动态区域和静态区域，其中，动态区域表征微泡射频信号，静态区域表征静态组织射频信号；最后通过过滤静态区域的射频信号，得到微泡射频信号。帧间差分(Differential Imaging,DI)处理是一种简单的时域差分处理，利用相邻两个时刻组织相对静止、微泡存在运动的状态差异，达到消除组织信号且保留运动微泡信号的目的。DI处理方法除了能够获取运动微泡的信号，还可检测到存在***以及溶解等行为的微泡信号。
参照图3，在本申请一个实施例中，步骤S230中所提及的奇异值分解滤波处理，包括以下步骤：
步骤S310：计算多帧空间射频信号的奇异值，得到奇异值集合。
步骤S320：根据预设统计参数确定奇异值集合中的目标奇异值。
步骤S330：根据各个目标奇异值提取各个空间射频信号中的目标特征信号,得到微泡射频信号。
需要说明的是，奇异值分解(Singular Value Decomposition，SVD)滤波，是一种多帧图像数据批处理的方法，图像矩阵的奇异值及其特征空间反映了图像中的不同成分和特征。从物理角度来看，SVD相当于把一个矩阵分解为了若干个子信号空间的线性组合，其中奇异值就是各个子信号空间对应的能量大小。在未作微泡提取的数据中，主要有三个成分：组织信号、微泡信号和噪声信号。其中，组织信号的能量最强，对应较大奇异值；微泡信号次之，对应中间大小奇异值；噪声信号最弱，对应较小奇异值。因此，利用预设统计参数确定微泡信号对应的奇异值，然后通过这些奇异值和其对应的特征向量来得到微泡射频信号。
在本申请的一个实施例中，图3中的步骤S320，包括以下步骤：
通过将奇异值集合中各个奇异值与预设的阈值范围进行对比，将落入阈值范围的奇异值确定为目标奇异值。
在本实施例中，假设空间射频信号的帧数为P，P帧空间射频信号通过SVD滤波后得到的奇异值分别记为λ1,λ2,...,λP，P个奇异值降序排列。预设统计参数包括信号阈值a和噪声阈值b，其中信号阈值a大于噪声阈值b。当奇异值大于或者等于信号阈值a时，被认为是组织信号对应的奇异值，应放弃；当奇异值小于或者等于噪声阈值b时，被认为是噪声信号对应的奇异值，应放弃；当奇异值大于噪声阈值b且小于信号阈值a时，被认为是微泡信号对应的奇异值，应保留；也即预设的阈值范围为(信号阈值a，噪声阈值b)。此时根据预设的阈值范围所确定的需要保留的目标奇异值，记为
Figure BDA0003603526020000071
目标奇异值对应的左奇异矢量、右奇异矢量分别记为
Figure BDA0003603526020000072
其中R为目标奇异值的数量，R小于P。然后根据各个目标奇异值提取的目标特征信号分别为
Figure BDA0003603526020000073
Figure BDA0003603526020000074
最后所得到的微泡射频信号为
Figure BDA0003603526020000075
其中，H为转置符号。
在本申请的另一个实施例中，图3中的步骤S320，包括以下步骤：
将奇异值集合中的奇异值依据大小进行排序，并获取排序序列中预设数量的中间部分奇异值作为目标奇异值。
可以理解的是，区别于上述实施例各个奇异值与预设的阈值范围进行对比来确定目标奇异值的方法，本实施例先将各个奇异值依据大小进行排序，再从中间部分选取预设数量的奇异值作为目标奇异值，能够避免预设的阈值范围选取不合适所带来的精确性不足的问题。可以理解的是，以上两种确定为目标奇异值的方式的侧重点不同，本领域技术人员可以在实际应用中按照实际情况任意选取其中一种。
参照图4，在本申请一个实施例中，图2中的步骤S250，包括以下步骤：
步骤S410：对第一微泡信号中的多帧微泡射频信号进行微泡定位处理，得到多帧第一微泡定位数据。
步骤S420：对第二微泡信号中的多帧微泡射频信号进行微泡定位处理，得到多帧第二微泡定位数据。
步骤S430：对多帧第一微泡定位数据和多帧第二微泡定位数据进行融合处理，得到目标微血管图像数据。
其中，微泡定位处理为根据微泡射频信号在空间上的分布情况，确定微泡在空间上的定位，以得到对应微泡射频信号的微泡定位数据。
可以理解的是，在本实施例中，先分别对第一微泡信号和对第二微泡信号进行微泡定位处理，再进行融合处理，可以使得所得到的目标微血管图像数据兼具两种时空滤波处理的优点，能够提高对微血管成像的精确性，提升微血管的成像效果。其中，微泡定位处理可以使用现有技术中已有的处理方法，本申请对此不再进行详细介绍。
参照图5，在本申请的一个实施例中，图4中的步骤S430，包括以下步骤：
步骤S510：根据各帧第一微泡定位数据及其对应帧序的各帧第二微泡定位数据，确定微血管等级分布数据。
具体地，步骤S510包括：根据各帧第一微泡定位数据及其对应帧序的各帧第二微泡定位数据，确定微血管等级，进而得到多帧微血管等级分布数据。
步骤S520：对微血管等级分布数据进行图像处理，得到目标微血管图像数据。
具体地，参照图6，图5中的步骤S520包括以下步骤：
步骤S610：分别确定各帧微血管等级分布数据对应的显示参数，得到多帧空间显示分布数据。
步骤S620：根据多帧空间显示分布数据得到目标微血管图像数据。
可以理解的是，各帧空间显示分布数据包括有各帧微血管等级分布数据对应的显示参数。因此，可以对多帧空间显示分布数据进行多帧追踪累加，从而得到目标微血管图像数据。
具体地，在一实施例中，第一微泡定位数据和第二微泡定位数据为二值化矩阵数据；
当第一微泡定位数据中的目标元素的值为第一数据值，且第二微泡定位数据中与目标元素对应的元素的值为第一数据值，微血管等级为一级微血管；
当第一微泡定位数据中的目标元素的值为第二数据值，且第二微泡定位数据中与目标元素对应的元素的值为第一数据值，微血管等级为二级微血管；
当第一微泡定位数据中的目标元素的值为第一数据值，且第二微泡定位数据中与目标元素对应的元素的值为第二数据值，微血管等级为三级微血管；
当第一微泡定位数据中的目标元素的值为第二数据值，且第二微泡定位数据中与目标元素对应的元素的值为第二数据值，微血管等级为无微血管；
其中，第一数据值表征空间射频信号中对应目标元素的空间位置存在微泡射频信号，第二数据值表征空间射频信号中对应目标元素的空间位置没有微泡射频信号，二级微血管表征的血管直径大于三级微血管表征的血管直径且小于一级微血管表征的血管直径。
在本实施例中，各帧的第一微泡定位数据的对应的二值化矩阵数据记为Location_1，第一微泡定位数据中的目标元素记为Location_1(h,w)，第一数据值示例性取值为1，即Location_1(h,w)＝1时空间射频信号中对应目标元素的空间位置存在微泡射频信号；第一数据值示例性取值为0，即Location_1(h,w)＝0时空间射频信号中对应目标元素的空间位置没有微泡射频信号；同理，各帧的第二微泡定位数据的对应的二值化矩阵数据记为Location_2，第二微泡定位数据中的目标元素记为Location_2(h,w)，第一数据值示例性取值为1，第一数据值示例性取值为0；其中，h表示目标元素的图像深度，w表示目标元素的图像宽度。
对两组微泡定位数据的定位结果进行融合，融合处理后包括以下四种情况：
情况1：Location_1(h,w)＝1且Location_2(h,w)＝1，对应微血管等级为最粗的一级微血管，则对应的显示参数Location(h,w)＝1*Enhance_1；
情况2：Location_1(h,w)＝1且Location_2(h,w)＝1，对应微血管等级为SVD滤波连续性较差的二级微血管，则对应的显示参数Location(h,w)＝1*Enhance_2；
情况3：Location_1(h,w)＝1且Location_2(h,w)＝0，对应微血管等级为DI缺失的三级微血管，则对应的显示参数Location(h,w)＝1*Enhance_3；
情况4：Location_1(h,w)＝0且Location_2(h,w)＝0，对应无微泡区域，微血管等级为无微血管，则对应的显示参数Location(h,w)＝0；
其中，Enhance_1、Enhance_2、Enhance_3可以表示为显示亮度参数，可以通过优化Enhance_1、Enhance_2、Enhance_3的取值，从而达到不同的成像效果。
可以理解的是，对融合处理得到多帧空间显示分布数据进行多帧追踪累加，从而获得超分辨图像结果。图15是本申请一个实施例提供的超声波成像效果与SVD滤波成像效果、DI处理成像效果的对比示意图，可以看到，本实施例提供的超声波成像方法，兼具了SVD滤波处理和DI处理两种处理方法的优点，提升了微血管成像的层次感。
在本申请的另一个实施例中，图4中的步骤S430，包括以下步骤：
将多帧第一微泡定位数据和对应帧序的多帧第二微泡定位数据在空间上分别进行加权求和，得到多帧目标微泡定位数据，根据目标微泡定位数据确定微血管分布情况，得到目标微血管图像数据。
区别于上述图5中的实施例通过确定微血管等级得到目标微血管图像数据的方式，本实施例通过将多帧第一微泡定位数据和对应帧序的多帧第二微泡定位数据在空间上分别进行加权求和，从而得到多帧目标微泡定位数据，可以使得目标微泡定位数据在空间分布上更加连续平滑，得到的微血管分布情况和目标微血管图像数据的成像效果更好。
在本申请的又一个实施例中，图4中的步骤S430，包括以下步骤：
根据多帧第一微泡定位数据确定第一微血管图像数据，根据多帧第二微泡定位数据确定第二微血管图像数据，根据第一微血管图像数据和第二微血管图像数据得到目标微血管图像数据。
可以理解的是，上述图5先确定微血管等级的实施例以及上述对第一微泡定位数据和第二微泡定位数据在空间上进行加权求和的实施例，均是先对两组微泡定位数据进行融合处理，再形成微血管图像数据。区别于以上两种实施例，本实施例先根据第一微泡定位数据和第二微泡定位数据分别形成两组微血管图像数据，再将两组微血管图像数据进行融合，得到目标微血管图像数据。其中，第一微血管图像数据和第二微血管图像数据，一般均是微血管的亮度显示数据，将两组微血管图像数据进行叠加融合，可以使得成像的图像兼具两种时空滤波处理的优点，能够提高对微血管成像的精确性，提升微血管的成像效果。
参照图7，在本申请一个实施例中，图2中的步骤S250，包括以下步骤：
步骤S710：对第一微泡信号和第二微泡信号进行融合处理，得到目标微泡信号，目标微泡信号包括多帧融合微泡射频信号。
步骤S720：对目标微泡信号中的各帧融合微泡射频信号进行微泡定位处理，得到多帧目标微泡定位数据。
步骤S730：根据多帧目标微泡定位数据确定微血管分布情况，得到目标微血管图像数据。
可以理解的是，上述图4的实施例，先对第一微泡信号、第二微泡信号中的多帧微泡射频信号分别进行微泡定位处理，再将得到的两组微泡定位数据进行融合处理，得到目标微血管图像数据。本实施例则是先对第一微泡信号、第二微泡信号两组微泡信号进行融合处理，再进行微泡定位处理，进而得到目标微血管图像数据。需要说明的是，两种不同顺序的处理方式，均能得到较好的成像效果，本领域技术人员可以根据实际情况进行选用。
参照图8，在本申请一个实施例中，图2中的步骤S250，包括以下步骤：
步骤S810：根据第一微泡信号确定第一微血管图像数据，根据第二微泡信号确定第二微血管图像数据。
步骤S820：根据第一微血管图像数据和第二微血管图像数据得到目标微血管图像数据。
可以理解的是，区别于上述图4的实施例和上述图7的实施例提供的两种融合处理方式，本实施例直接根据两组微泡信号分别生成微血管图像数据，然后对微血管图像数据进行融合，例如将两组微血管图像数据在空间上进行叠加，或者进行加权求和，得到目标微血管图像数据。
另外，参照图9，图9所示为本申请另一个实施例的超声波成像方法的流程图。超声波成像方法可以包括如下步骤：
步骤S910：向目标组织发射第一超声波，接收目标组织返回的第一超声波的回波，以获得第一超声回波信号。
步骤S920：根据第一超声回波信号，得到目标组织的组织结构信号，组织结构信号包括多帧空间射频信号。
步骤S930：对组织结构信号进行第一微泡信号处理，以提取各帧空间射频信号中的微泡射频信号，得到第一微泡信号。
步骤S940：对组织结构信号进行第二微泡信号处理，以提取各帧空间射频信号中的微泡射频信号，得到第二微泡信号，其中，第一微泡信号处理和第二微泡信号处理的处理效果不同，第一微泡信号处理和第二微泡信号处理均为时空滤波处理，时空滤波处理为根据多帧空间射频信号之间数据的比较来过滤空间射频信号中的静态组织射频信号，以在空间射频信号中提取微泡射频信号。
步骤S950：通过对第一微泡信号和第二微泡信号进行融合处理，得到目标微血管图像数据。
步骤S960：根据目标微血管图像数据显示微血管图像。
在本实施例中，通过采用第一微泡信号处理和第二微泡信号处理两种处理效果不同的时空滤波处理方法，来对组织结构信号进行处理，以过滤空间射频信号中的静态组织射频信号，提取微泡射频信号，得到第一微泡信号和第二微泡信号两组微泡信号，然后进行融合处理，可以使得成像的图像兼具两种时空滤波处理的优点，能够提高对微血管成像的精确性，提升微血管的成像效果。
在一些实施例中，第一微泡信号处理为时域差分滤波处理或者奇异值分解滤波处理；第二微泡信号处理为时域差分滤波处理或者奇异值分解滤波处理。
可以理解的是，当第一微泡信号处理为时域差分滤波处理且第二微泡信号处理为奇异值分解滤波处理，以及当第一微泡信号处理为奇异值分解滤波处理且第二微泡信号处理为时域差分滤波处理，第一微泡信号处理和第二微泡信号处理为不同种类的时空滤波处理；当第一微泡信号处理和第二微泡信号处理均为时域差分滤波处理，以及当第一微泡信号处理和第二微泡信号处理均为奇异值分解滤波处理，第一微泡信号处理和第二微泡信号处理需要选用不同的滤波参数，第一微泡信号处理和第二微泡信号处理为不同滤波参数的同种类时空滤波处理。
需要说明的是，本实施例中对于时域差分滤波处理的具体描述，可以参照图2所示的实施例对步骤S240的具体描述；本实施例中对于奇异值分解滤波处理的具体描述，可以参照图2所示的实施例对步骤S230的具体描述，例如参照图3所示的实施例的具体描述；本实施例中对于步骤S950的具体描述，可以参照图2所示的实施例对步骤S250的具体描述，例如参照图4至图8的相关具体描述。
另外，参照图10，图10所示为本申请又一个实施例的超声波成像方法的流程图。超声波成像方法可以包括如下步骤：
步骤S1010：向目标组织发射第一超声波，接收目标组织返回的第一超声波的回波，以获得第一超声回波信号。
步骤S1020：根据第一超声回波信号，得到目标组织的组织结构信号，组织结构信号包括多帧空间射频信号。
步骤S1030：对组织结构信号进行多种时空滤波处理，提取微泡射频信号，得到多组微泡信号，其中，时空滤波处理为根据多帧空间射频信号之间数据的比较来过滤空间射频信号中的静态组织射频信号，在空间射频信号中提取微泡射频信号，进而得到微泡信号。
步骤S1040：对各组微泡信号进行融合处理，得到目标微血管图像数据。
步骤S1050：根据目标微血管图像数据显示微血管图像。
可以理解的是，图9所示的实施例，对组织结构信号分别进行第一微泡信号处理和第二微泡信号处理两种处理效果不同的时空滤波处理。而本实施例则是可以进行多种时空滤波处理，例如对组织结构信号进行三种微泡信号处理，分别为第一微泡信号处理、第二微泡信号处理和第三微泡信号处理，然后得到第一微泡信号、第二微泡信号和第三微泡信号，最后对三组微泡信号进行融合处理，得到目标微血管图像数据，并显示微血管图像，可以使得成像的图像兼具三种时空滤波处理的优点，能够提高对微血管成像的精确性，提升微血管的成像效果。对组织结构信号进行四种微泡信号处理的情况可以以此类推，此处不再详述。
本实施例中所采用的多个时空滤波处理，可以是包括不同滤波参数的同种类时空滤波处理，也可以是包括不同种类的时空滤波处理，例如包括时域差分滤波处理或者奇异值分解滤波处理。
另外，如图11所示为一种图像显示装置的结构示意图。其中，图像显示装置1100包括处理器1110以及显示器1120；处理器1110用于获取并处理组织结构信号，其中组织结构信号包括多帧空间射频信号；显示器1120用于显示微血管图像数据。
应理解，图11所示的图像显示装置1100所包括的部件只是示意性的，其可以包括更多或更少的部件。本发明对此不限定。图11所示的图像显示装置1100中的处理器1110还用于执行下文图12至图14中的超声波成像方法。
参照图12，图12所示为本申请一个实施例的超声波成像方法的流程图。超声波成像方法可以包括如下步骤：
步骤S1210：获取组织结构信号，组织结构信号包括多帧空间射频信号。
步骤S1220：对组织结构信号进行第一微泡信号处理，以提取各帧空间射频信号中的微泡射频信号，得到第一微泡信号。
步骤S1230：对组织结构信号进行第二微泡信号处理，以提取各帧空间射频信号中的微泡射频信号，得到第二微泡信号，其中，第一微泡信号处理和第二微泡信号处理的处理效果不同，第一微泡信号处理和第二微泡信号处理均为时空滤波处理，时空滤波处理为根据多帧空间射频信号之间数据的比较来过滤空间射频信号中的静态组织射频信号，在空间射频信号中提取微泡射频信号。
步骤S1240：通过对第一微泡信号和第二微泡信号进行融合处理，得到目标微血管图像数据。
步骤S1250：根据目标微血管图像数据显示微血管图像。
可以理解的是，本实施例提供的超声波成像方法，相较于图9所示的实施例提供的超声波成像方法，没有发射超声波和接收超声波的回波步骤，只有对获取到的组织结构信号进行处理，并得到目标微血管图像数据，显示微血管图像等信号处理步骤，因此一般应用于诸如图11所示的图像显示装置1100。
参照图13，图13所示为本申请另一个实施例的超声波成像方法的流程图。超声波成像方法可以包括如下步骤：
步骤S1310：获取组织结构信号，组织结构信号包括多帧空间射频信号。
步骤S1320：对组织结构信号进行奇异值分解滤波处理，以提取各帧空间射频信号中的微泡射频信号，得到第一微泡信号。
步骤S1330：对组织结构信号进行时域差分滤波处理，以提取各帧空间射频信号中的微泡射频信号，得到第二微泡信号。
步骤S1340：通过对第一微泡信号和第二微泡信号进行融合处理，得到目标微血管图像数据。
步骤S1350：根据目标微血管图像数据显示微血管图像。
可以理解的是，本实施例提供的超声波成像方法，相较于图2所示的实施例提供的超声波成像方法，没有发射超声波和接收超声波的回波步骤，只有对获取到的组织结构信号进行处理，并得到目标微血管图像数据，显示微血管图像等信号处理步骤，因此一般应用于诸如图11所示的图像显示装置1100。
参照图14，图14所示为本申请又一个实施例的超声波成像方法的流程图。超声波成像方法可以包括如下步骤：
步骤S1410：获取组织结构信号，组织结构信号包括多帧空间射频信号；
步骤S1420：对组织结构信号进行多种时空滤波处理，提取微泡射频信号，得到多组微泡信号，其中，时空滤波处理为根据多帧空间射频信号之间数据的比较来过滤空间射频信号中的静态组织射频信号，在空间射频信号中提取微泡射频信号，进而得到微泡信号；
步骤S1430：对各组微泡信号进行融合处理，得到目标微血管图像数据；
步骤S1440：根据目标微血管图像数据显示微血管图像。
可以理解的是，本实施例提供的超声波成像方法，相较于图10所示的实施例提供的超声波成像方法，没有发射超声波和接收超声波的回波步骤，只有对获取到的组织结构信号进行处理，并得到目标微血管图像数据，显示微血管图像等信号处理步骤，因此一般应用于诸如图11所示的图像显示装置1100。
本申请实施例提供了一种计算机存储介质，其上存储有计算机程序，应用于超声成像装置或者图像显示装置，该计算机程序被处理器执行时实现如上述实施例的超声波成像方法。
本申请实施例提供了一种计算机程序产品或计算机程序，该计算机程序产品或计算机程序包括计算机指令，该计算机指令存储在计算机可读存介质中。计算机设备的处理器从计算机可读存储介质读取该计算机指令，处理器执行该计算机指令，使得该计算机设备执行如上述实施例的超声波成像方法。
在本申请所提供的几个实施例中，应该理解到，所揭露的***，装置和方法，可以通过其他的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个***，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外，在本申请各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机装置(可以是个人计算机，服务器，或者网络装置等)执行本申请各个实施例方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(Read-OnlyMemory，简称ROM)、随机存取存储器(RandomAccessMemory，简称RAM)、磁碟或者光盘等各种可以存储程序代码的介质。
还应了解，本申请实施例提供的各种实施方式可以任意进行组合，以实现不同的技术效果。
以上是对本申请的较佳实施进行了具体说明，但本申请并不局限于上述实施方式，熟悉本领域的技术人员在不违背本申请精神的共享条件下还可作出种种等同的变形或替换，这些等同的变形或替换均包括在本申请权利要求所限定的范围内。

Claims (21)

1.一种超声波成像方法，其特征在于，所述方法包括：
向目标组织发射第一超声波，接收所述目标组织返回的所述第一超声波的回波，以获得第一超声回波信号；
根据所述第一超声回波信号，得到所述目标组织的组织结构信号，所述组织结构信号包括多帧空间射频信号；
对所述组织结构信号进行奇异值分解滤波处理，以提取各帧所述空间射频信号中的微泡射频信号，得到第一微泡信号；
对所述组织结构信号进行时域差分滤波处理，以提取各帧所述空间射频信号中的微泡射频信号，得到第二微泡信号；
通过对所述第一微泡信号和所述第二微泡信号进行融合处理，得到目标微血管图像数据；
根据所述目标微血管图像数据显示微血管图像。
2.一种超声波成像方法，其特征在于，所述方法包括：
向目标组织发射第一超声波，接收所述目标组织返回的所述第一超声波的回波，以获得第一超声回波信号；
根据所述第一超声回波信号，得到所述目标组织的组织结构信号，所述组织结构信号包括多帧空间射频信号；
对所述组织结构信号进行第一微泡信号处理，以提取各帧所述空间射频信号中的微泡射频信号，得到第一微泡信号；
对所述组织结构信号进行第二微泡信号处理，以提取各帧所述空间射频信号中的微泡射频信号，得到第二微泡信号，其中，所述第一微泡信号处理和所述第二微泡信号处理的处理效果不同，所述第一微泡信号处理和所述第二微泡信号处理均为时空滤波处理，所述时空滤波处理为根据多帧所述空间射频信号之间数据的比较来过滤所述空间射频信号中的静态组织射频信号，以在所述空间射频信号中提取所述微泡射频信号；
通过对所述第一微泡信号和所述第二微泡信号进行融合处理，得到目标微血管图像数据；
根据所述目标微血管图像数据显示微血管图像。
3.根据权利要求2所述的一种超声波成像方法，其特征在于，所述第一微泡信号处理和第二微泡信号处理为不同滤波参数的同种类时空滤波处理，或者为不同种类的时空滤波处理。
4.根据权利要求2所述的一种超声波成像方法，其特征在于，所述第一微泡信号处理为时域差分滤波处理或者奇异值分解滤波处理；所述第二微泡信号处理为时域差分滤波处理或者奇异值分解滤波处理。
5.根据权利要求1或4所述的一种超声波成像方法，其特征在于，所述时域差分滤波处理包括：
根据两帧预设间距的所述空间射频信号中的其中一帧空间射频信号与另一帧作差，过滤静态组织射频信号，得到所述微泡射频信号。
6.根据权利要求5所述的一种超声波成像方法，其特征在于，所述作差包括幅度差、相位差或频率差中的至少一种。
7.根据权利要求1或4所述的一种超声波成像方法，其特征在于，所述奇异值分解滤波处理包括：
计算多帧所述空间射频信号的奇异值，得到奇异值集合；
根据预设统计参数确定所述奇异值集合中的目标奇异值；
根据各个所述目标奇异值提取各个所述空间射频信号中的目标特征信号,得到所述微泡射频信号。
8.根据权利要求7所述的一种超声波成像方法，其特征在于，所述根据预设统计参数确定所述奇异值集合中的目标奇异值，包括以下之一：
通过将所述奇异值集合中各个所述奇异值与预设的阈值范围进行对比，将落入所述阈值范围的奇异值确定为所述目标奇异值；
或者，将所述奇异值集合中的奇异值依据大小进行排序，并获取排序序列中预设数量的中间部分奇异值作为所述目标奇异值。
9.根据权利要求1或2所述的一种超声波成像方法，其特征在于，所述通过对所述第一微泡信号和所述第二微泡信号进行融合处理，得到目标微血管图像数据，包括：
对所述第一微泡信号中的多帧微泡射频信号进行微泡定位处理，得到多帧第一微泡定位数据；
对所述第二微泡信号中的多帧微泡射频信号进行微泡定位处理，得到多帧第二微泡定位数据；
对所述多帧第一微泡定位数据和所述多帧第二微泡定位数据进行融合处理，得到所述目标微血管图像数据；
其中，所述微泡定位处理为根据微泡射频信号在空间上的分布情况，确定微泡在空间上的定位，以得到对应所述微泡射频信号的微泡定位数据。
10.根据权利要求9所述的一种超声波成像方法，其特征在于，所述对所述多帧第一微泡定位数据和所述多帧第二微泡定位数据进行融合处理，得到所述目标微血管图像数据，包括以下至少之一：
根据各帧第一微泡定位数据及其对应帧序的各帧第二微泡定位数据，确定微血管等级分布数据，对所述微血管等级分布数据进行图像处理，得到所述目标微血管图像数据；
或者，将所述多帧第一微泡定位数据和对应帧序的所述多帧第二微泡定位数据在空间上分别进行加权求和，得到多帧目标微泡定位数据，根据所述目标微泡定位数据确定微血管分布情况，得到所述目标微血管图像数据；
或者，根据多帧所述第一微泡定位数据确定第一微血管图像数据，根据多帧所述第二微泡定位数据确定第二微血管图像数据，根据所述第一微血管图像数据和所述第二微血管图像数据得到所述目标微血管图像数据。
11.根据权利要求10所述的一种超声波成像方法，其特征在于，所述根据各帧第一微泡定位数据及其对应帧序的各帧第二微泡定位数据，确定微血管等级分布数据，包括：
根据各帧第一微泡定位数据及其对应帧序的各帧第二微泡定位数据，确定微血管等级，进而得到多帧微血管等级分布数据；
所述对所述微血管等级分布数据进行图像处理，得到目标微血管图像数据，包括：
分别确定各帧微血管等级分布数据对应的显示参数，得到多帧空间显示分布数据；
根据所述多帧空间显示分布数据得到所述目标微血管图像数据。
12.根据权利要求11所述的一种超声波成像方法，其特征在于，所述第一微泡定位数据和所述第二微泡定位数据为二值化矩阵数据；
当所述第一微泡定位数据中的目标元素的值为第一数据值，且所述第二微泡定位数据中与所述目标元素对应的元素的值为第一数据值，所述微血管等级为一级微血管；
当所述第一微泡定位数据中的目标元素的值为第二数据值，且所述第二微泡定位数据中与所述目标元素对应的元素的值为第一数据值，所述微血管等级为二级微血管；
当所述第一微泡定位数据中的目标元素的值为第一数据值，且所述第二微泡定位数据中与所述目标元素对应的元素的值为第二数据值，所述微血管等级为三级微血管；
当所述第一微泡定位数据中的目标元素的值为第二数据值，且所述第二微泡定位数据中与所述目标元素对应的元素的值为第二数据值，所述微血管等级为无微血管；
其中，所述第一数据值表征所述空间射频信号中对应所述目标元素的空间位置存在所述微泡射频信号，所述第二数据值表征所述空间射频信号中对应所述目标元素的空间位置没有所述微泡射频信号，所述二级微血管表征的血管直径大于所述三级微血管表征的血管直径且小于所述一级微血管表征的血管直径。
13.根据权利要求1或2所述的一种超声波成像方法，其特征在于，所述通过对所述第一微泡信号和所述第二微泡信号进行融合处理，得到目标微血管图像数据，包括以下之一：
对所述第一微泡信号和所述第二微泡信号进行融合处理，得到目标微泡信号，所述目标微泡信号包括多帧融合微泡射频信号，对所述目标微泡信号中的各帧融合微泡射频信号进行微泡定位处理，得到多帧目标微泡定位数据，根据多帧所述目标微泡定位数据确定微血管分布情况，得到所述目标微血管图像数据；
或者，根据所述第一微泡信号确定第一微血管图像数据，根据所述第二微泡信号确定第二微血管图像数据，根据所述第一微血管图像数据和所述第二微血管图像数据得到所述目标微血管图像数据。
14.一种超声波成像方法，其特征在于，所述方法包括：
向目标组织发射第一超声波，接收所述目标组织返回的所述第一超声波的回波，以获得第一超声回波信号；
根据所述第一超声回波信号，得到所述目标组织的组织结构信号，所述组织结构信号包括多帧空间射频信号；
对所述组织结构信号进行多种时空滤波处理，提取微泡射频信号，得到多组微泡信号，其中，所述时空滤波处理为根据多帧所述空间射频信号之间数据的比较来过滤所述空间射频信号中的静态组织射频信号，在所述空间射频信号中提取所述微泡射频信号，进而得到所述微泡信号；
对各组所述微泡信号进行融合处理，得到目标微血管图像数据；
根据所述目标微血管图像数据显示微血管图像。
15.根据权利要求14所述的一种超声波成像方法，其特征在于，所述多个时空滤波处理，包括不同滤波参数的同种类时空滤波处理，或者包括不同种类的时空滤波处理。
16.根据权利要求14所述的一种超声波成像方法，其特征在于，所述时空滤波处理为时域差分滤波处理或者奇异值分解滤波处理。
17.一种超声波成像方法，其特征在于，所述方法包括：
获取组织结构信号，所述组织结构信号包括多帧空间射频信号；
对所述组织结构信号进行第一微泡信号处理，以提取各帧所述空间射频信号中的微泡射频信号，得到第一微泡信号；
对所述组织结构信号进行第二微泡信号处理，以提取各帧所述空间射频信号中的微泡射频信号，得到第二微泡信号，其中，所述第一微泡信号处理和所述第二微泡信号处理的处理效果不同，所述第一微泡信号处理和所述第二微泡信号处理均为时空滤波处理，所述时空滤波处理为根据多帧所述空间射频信号之间数据的比较来过滤所述空间射频信号中的静态组织射频信号，在所述空间射频信号中提取所述微泡射频信号；
通过对所述第一微泡信号和所述第二微泡信号进行融合处理，得到目标微血管图像数据；
根据所述目标微血管图像数据显示微血管图像。
18.一种超声波成像方法，其特征在于，所述方法包括：
获取组织结构信号，所述组织结构信号包括多帧空间射频信号；
对所述组织结构信号进行奇异值分解滤波处理，以提取各帧所述空间射频信号中的微泡射频信号，得到第一微泡信号；
对所述组织结构信号进行时域差分滤波处理，以提取各帧所述空间射频信号中的微泡射频信号，得到第二微泡信号；
通过对所述第一微泡信号和所述第二微泡信号进行融合处理，得到目标微血管图像数据；
根据所述目标微血管图像数据显示微血管图像。
19.一种超声波成像方法，其特征在于，所述方法包括：
获取组织结构信号，所述组织结构信号包括多帧空间射频信号；
对所述组织结构信号进行多种时空滤波处理，提取微泡射频信号，得到多组微泡信号，其中，所述时空滤波处理为根据多帧所述空间射频信号之间数据的比较来过滤所述空间射频信号中的静态组织射频信号，在所述空间射频信号中提取所述微泡射频信号，进而得到所述微泡信号；
对各组所述微泡信号进行融合处理，得到目标微血管图像数据；
根据所述目标微血管图像数据显示微血管图像。
20.一种超声成像装置，其特征在于，包括：
超声探头；
发射/接收电路，所述发射/接收电路用于控制所述超声探头向目标组织发射超声波并接收超声回波，获得超声回波信号；
处理器，所述处理器用于处理所述超声回波信号，获得所述目标组织的组织结构信号以及通过对所述组织结构信号进行微泡信号处理，获得微泡信号并根据所述微泡信号得到目标微血管图像数据；
显示器，所述显示器用于显示所述微血管图像数据；
所述处理器还用于执行上述权利要求1至16中任意一项所述的超声波成像方法。
21.一种图像显示装置，其特征在于，包括：
处理器，所述处理器用于获取并处理组织结构信号，其中所述组织结构信号包括多帧空间射频信号；
显示器，所述显示器用于显示所述微血管图像数据；
所述处理器还用于执行上述权利要求17至19所述的超声波成像方法。
CN202210409389.XA 2022-04-19 2022-04-19 超声波成像方法、超声成像设备及图像显示装置 Pending CN114903519A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210409389.XA CN114903519A (zh) 2022-04-19 2022-04-19 超声波成像方法、超声成像设备及图像显示装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210409389.XA CN114903519A (zh) 2022-04-19 2022-04-19 超声波成像方法、超声成像设备及图像显示装置

Publications (1)

Publication Number Publication Date
CN114903519A true CN114903519A (zh) 2022-08-16

Family

ID=82765556

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210409389.XA Pending CN114903519A (zh) 2022-04-19 2022-04-19 超声波成像方法、超声成像设备及图像显示装置

Country Status (1)

Country Link
CN (1) CN114903519A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115530875A (zh) * 2022-10-26 2022-12-30 杭州永锦科技有限公司 超声波成像方法、装置、设备及可读存储介质

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115530875A (zh) * 2022-10-26 2022-12-30 杭州永锦科技有限公司 超声波成像方法、装置、设备及可读存储介质

Similar Documents

Publication Publication Date Title
Van Sloun et al. Deep learning in ultrasound imaging
CN104042247B (zh) 使用自适应时间阶段的超声arfi位移成像
US11589840B2 (en) Methods for super-resolution ultrasound imaging of microvessels
KR101906916B1 (ko) 지식-기반 초음파 이미지 개선
US8684934B2 (en) Adaptively performing clutter filtering in an ultrasound system
JP6030825B2 (ja) 電磁気ノイズ除去のための超音波法
EP2016905B1 (en) Ultrasound diagnostic apparatus
EP0884606A2 (en) Ultrasonic diagnostic image flash suppression technique
EP1967867A2 (en) Ultrasound system and method of forming ultrasound images
CN100457045C (zh) 超声波诊断装置和图像处理装置
KR101188593B1 (ko) 복수의 3차원 초음파 영상을 제공하는 초음파 시스템 및 방법
US20140050048A1 (en) Harmonic Ultrasound Imaging Using Synthetic Aperture Sequential Beamforming
CN115670520B (zh) 血管内超声成像方法、装置、计算机设备、存储介质
US11826201B2 (en) Ultrasound lesion assessment and associated devices, systems, and methods
CN114903519A (zh) 超声波成像方法、超声成像设备及图像显示装置
KR20170042936A (ko) 초음파 장치 및 초음파 신호 처리 방법
CN104768470A (zh) 超声波诊断装置
Jensen et al. Fast super resolution ultrasound imaging using the erythrocytes
CN112512436B (zh) 用于超声中的肺搏动检测的设备、***和方法
CN114129190A (zh) 全眼微血流的成像方法、装置及电子设备和存储介质
CN110072462A (zh) 用于对含血液的样品成像的方法和相关装置
Wang et al. A general deep learning model for ultrasound localization microscopy
Liu et al. Super-resolution ultrasound imaging by sparse Bayesian learning method
EP2341364B1 (en) Adaptive clutter filtering method and ultrasound system for the same
JP4982393B2 (ja) 画像フィルタリング装置、画像フィルタリングプログラム、画像フィルタリング方法及び超音波診断装置

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