CN109363722B - 彩色血流成像中运动伪像的抑制方法及设备 - Google Patents
彩色血流成像中运动伪像的抑制方法及设备 Download PDFInfo
- Publication number
- CN109363722B CN109363722B CN201811444026.XA CN201811444026A CN109363722B CN 109363722 B CN109363722 B CN 109363722B CN 201811444026 A CN201811444026 A CN 201811444026A CN 109363722 B CN109363722 B CN 109363722B
- Authority
- CN
- China
- Prior art keywords
- blood flow
- flow velocity
- color
- frames
- imaging
- 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.)
- Active
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 100
- 238000000034 method Methods 0.000 title claims abstract description 38
- 230000017531 blood circulation Effects 0.000 claims abstract description 294
- 238000004891 communication Methods 0.000 claims description 9
- 238000009825 accumulation Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 abstract description 15
- 230000002401 inhibitory effect Effects 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 6
- 230000001629 suppression Effects 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 210000001519 tissue Anatomy 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 239000000523 sample Substances 0.000 description 3
- 210000004204 blood vessel Anatomy 0.000 description 2
- 238000002592 echocardiography Methods 0.000 description 2
- 210000003743 erythrocyte Anatomy 0.000 description 2
- 210000003205 muscle Anatomy 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000029058 respiratory gaseous exchange Effects 0.000 description 2
- 238000002604 ultrasonography Methods 0.000 description 2
- 238000005311 autocorrelation function Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000002203 pretreatment Methods 0.000 description 1
- 230000003252 repetitive effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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/5269—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
- A61B8/5276—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts due to motion
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/06—Measuring blood flow
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)
- Hematology (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Image Analysis (AREA)
Abstract
本发明实施例提供了一种彩色血流成像中运动伪像的抑制方法及设备。其中，所述方法包括：获取彩色血流成像中的血流速度及相应的灰度图像，将所述血流速度中小于流速阈值的血流速度剔除，得到处理后血流速度，在所述彩色血流成像的相同位置，若在彩色血流成像时刻之前的若干帧图像的血流速度均为零，或仅有一帧图像的血流速度不为零，和/或所述相应的灰度图像的灰度值大于等于灰度直方图均值，则确定所述处理后血流速度为运动伪像；去除所述运动伪像，得到真实的血流速度。本发明实施例提供的彩色血流成像中运动伪像的抑制方法及设备，可以达到抑制彩色血流成像中运动伪像的技术效果。
Description
技术领域
本发明实施例涉及超声成像技术领域，尤其涉及一种彩色血流成像中运动伪像的抑制方法及设备。
背景技术
1842年奥地利物理学家多普勒首先提出一种物理现象，即振动源与接收器之间存在运动时，所接收到的振动频率因运动而发生改变。现在超声检查应用的超声探头就是振动源，被检查人的心脏及血流就是接收器，两者之间存在着振动频率的改变，被称之为“多普勒”效应。在这个基础上科学家以脉冲多普勒技术为基础，用运动目标检测(MTI)、自相关函数计算、数字扫描转换及彩色编码等技术，得到血流的彩色显像，因此将其称之为彩色多普勒血流成像。作为超声诊断技术中非常重要的一支，超声多普勒血流成像能够在较大范围内实时显示人体血流分布情况，并兼有组织结构成像和多普勒流速检测等功能，成为目前高端超声诊断设备中不可缺少的部分。彩色血流成像也常被俗称为“彩超”或“彩色多普勒”。
在理想情况下，我们认为超声回波是完全由红细胞的背向散射信号组成的。但人体中还有许多其他的反射源，如血管壁、肌肉和组织，同时人体的呼吸、心跳以及探头的移动等都可以产生多普勒频移。因此得到的每一个多普勒矢量信号是由三种信号夹杂在一起形成的。
第一种是在人体血液中的大量红细胞的背向散射形成，分布在发射频率附近的窄带信号。也就是计算血流参数的超声多普勒血流信号；第二种是由血管壁、肌肉、软组织等运动速度非常慢的人体器官或组织对超声波的反射形成的信号，这种信号具有频率低、幅度大的特点，被称为杂波(clutter)；第三种就是系统内外都有可能形成的噪声(比如白噪声)，还有探头的移动、人体的呼吸、心跳等产生的伪像，我们称之为运动伪像，主要表现为大片状或块状的不规则彩色信号。
上述第二种和第三种信号都不是我们需要的，因此需要去掉。关于第二种信号的去除，杂波滤波，之前有很多人做了这方面的研究，也取得了一定的效果，但是很难完全去掉第三种信号。另外，调节较高的速度标尺也会对运动伪像有一定的抑制作用，但是同时也会去掉一些低速血流信号。针对运动伪像的研究不是很多，大部分都是从灰度图、血流、能量、方差等多个方面入手，处理起来比较复杂。因此，找到一种能够有效去除彩色血流成像中的各种运动伪像，获取精确的血流图像，就成为业界亟待解决的主要问题。
发明内容
针对现有技术存在的上述问题，本发明实施例提供了一种彩色血流成像中运动伪像的抑制方法及设备。
第一方面，本发明的实施例提供了一种彩色血流成像中运动伪像的抑制方法，包括：获取彩色血流成像中的血流速度及相应的灰度图像，将所述血流速度中小于流速阈值的血流速度剔除，得到处理后血流速度，在所述彩色血流成像的相同位置，若在彩色血流成像时刻之前的若干帧图像的血流速度均为零，或仅有一帧图像的血流速度不为零，和/或所述相应的灰度图像的灰度值大于等于灰度直方图均值，则确定所述处理后血流速度为运动伪像；去除所述运动伪像，得到真实的血流速度。
进一步地，所述灰度直方图均值，包括：
A＝[0,1,...,255]；f＝[f1,f2,f3,...,f256]
其中，hist_avg为灰度直方图均值，A为灰度数据，A(i)为A的第i个元素，f为灰度出现的次数，fi为f的第i个元素。
进一步地，所述在所述彩色血流成像的相同位置，若在彩色血流成像时刻之前的若干帧图像的血流速度均为零，或仅有一帧图像的血流速度不为零，和/或所述相应的灰度图像的灰度值大于等于灰度直方图均值，则确定所述处理后血流速度为运动伪像，相应地，确定真实血流速度包括：
将所述处理后血流速度及若干帧对比血流速度二值化，若二值化后的处理后血流速度为1，且二值化后的若干帧对比血流速度之和大于2，则确定所述处理后血流速度为真实血流速度；其中，所述若干帧对比血流速度为在所述彩色血流成像的相同位置，在彩色血流成像时刻之前的若干帧图像的血流速度，所述二值化为1、0二值化。
进一步地，所述在所述彩色血流成像的相同位置，若在彩色血流成像时刻之前的若干帧图像的血流速度均为零，或仅有一帧图像的血流速度不为零，和/或所述相应的灰度图像的灰度值大于等于灰度直方图均值，则确定所述处理后血流速度为运动伪像，相应地，确定真实血流速度包括：
将所述处理后血流速度及若干帧对比血流速度二值化，若二值化后的处理后血流速度为1，且二值化后的若干帧对比血流速度之和等于2，且所述相应的灰度图像的灰度值小于灰度直方图均值，则确定所述处理后血流速度为真实血流速度；其中，所述若干帧对比血流速度为在所述彩色血流成像的相同位置，在彩色血流成像时刻之前的若干帧图像的血流速度，所述二值化为1、0二值化。
进一步地，获取所述若干帧对比血流速度之和的方法，包括：将二值化后的若干帧对比血流速度进行算术累加，得到所述若干帧对比血流速度之和。
进一步地，获取所述若干帧对比血流速度之和的方法，包括：若所述若干帧对比血流速度中的一个血流速度是真实血流速度，则将计数加1，直到遍历完全部所述若干帧对比血流速度，得到的所述计数的最终值为所述若干帧对比血流速度之和。
第二方面，本发明的实施例提供了一种彩色血流成像中运动伪像的抑制装置，包括：
血流运动伪像确定模块，用于获取彩色血流成像中的血流速度及相应的灰度图像，将所述血流速度中小于流速阈值的血流速度剔除，得到处理后血流速度，在所述彩色血流成像的相同位置，若在彩色血流成像时刻之前的若干帧图像的血流速度均为零，和/或所述相应的灰度图像的灰度值大于等于灰度直方图均值，则确定所述处理后血流速度为运动伪像；
血流运动伪像去除模块，用于去除所述运动伪像，得到真实的血流速度。
第三方面，本发明的实施例提供了一种电子设备，包括：
至少一个处理器；以及
与处理器通信连接的至少一个存储器，其中：
存储器存储有可被处理器执行的程序指令，处理器调用程序指令能够执行第一方面的各种可能的实现方式中任一种可能的实现方式所提供的彩色血流成像中运动伪像的抑制方法。
第四方面，本发明的实施例提供了一种非暂态计算机可读存储介质，非暂态计算机可读存储介质存储计算机指令，计算机指令使计算机执行第一方面的各种可能的实现方式中任一种可能的实现方式所提供的彩色血流成像中运动伪像的抑制方法。
本发明实施例提供的彩色血流成像中运动伪像的抑制方法及设备，通过采用流速阈值对血流速度进行筛选，再将经过筛选的血流速度与之前若干帧血流速度图像进行对比，可以精确定位真实血流速度及血流速度的运动伪像，进而将彩色血流成像中的血流速度的伪像去除，达到抑制彩色血流成像中运动伪像的技术效果。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图做一简单的介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的彩色血流成像中运动伪像的抑制方法流程图；
图2为现有技术提供的彩色血流成像系统的结构示意图；
图3为本发明实施例提供的彩色血流成像系统的结构示意图；
图4为本发明实施例提供的另一彩色血流成像系统的结构示意图；
图5为本发明实施例提供的彩色血流成像中运动伪像的抑制效果示意图；
图6为本发明实施例提供的彩色血流成像中运动伪像的抑制装置结构示意图；
图7为本发明实施例提供的电子设备的实体结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。另外，本发明提供的各个实施例或单个实施例中的技术特征可以相互任意结合，以形成可行的技术方案，但是必须是以本领域普通技术人员能够实现为基础，当技术方案的结合出现相互矛盾或无法实现时，应当认为这种技术方案的结合不存在，也不在本发明要求的保护范围之内。
本发明实施例提供了一种彩色血流成像中运动伪像的抑制方法，参见图1，该方法包括：
101、获取彩色血流成像中的血流速度及相应的灰度图像，将所述血流速度中小于流速阈值的血流速度剔除，得到处理后血流速度，在所述彩色血流成像的相同位置，若在彩色血流成像时刻之前的若干帧图像的血流速度均为零，或仅有一帧图像的血流速度不为零，和/或所述相应的灰度图像的灰度值大于等于灰度直方图均值，则确定所述处理后血流速度为运动伪像；
102、去除所述运动伪像，得到真实的血流速度。
在本发明的各个实施例中，与血流速度相关的数据主要是指频差，血流速度与频差的关系可以表示如下：
其中，f0为发射频率，取决于仪器发射参数，单位是MHz；fd为频差，仪器根据接收信号提取，与f0相关，单位Hz；C为声速，1540m/s；θ为声速与目标运动方向夹角。在此基础上，在另一实施例中流速阈值可以由频差决定，具体可以取9Hz(通常取5Hz-12Hz)。
在上述实施例的基础上，本发明实施例中提供的彩色血流成像中运动伪像的抑制方法，所述灰度直方图均值，包括：
A＝[0,1,...,255]；f＝[f1,f2,f3,...,f256]
其中，hist_avg为灰度直方图均值，A为灰度数据，A(i)为A的第i个元素，f为灰度出现的次数，fi为f的第i个元素。在另一实施例中，上述灰度直方图均值可以取73.4951。
在上述实施例的基础上，本发明实施例中提供的彩色血流成像中运动伪像的抑制方法，所述在所述彩色血流成像的相同位置，若在彩色血流成像时刻之前的若干帧图像的血流速度均为零，或仅有一帧图像的血流速度不为零，和/或所述相应的灰度图像的灰度值大于等于灰度直方图均值，则确定所述处理后血流速度为运动伪像，相应地，确定真实血流速度包括：将所述处理后血流速度及若干帧对比血流速度二值化，若二值化后的处理后血流速度为1，且二值化后的若干帧对比血流速度之和大于2，则确定所述处理后血流速度为真实血流速度；其中，所述若干帧对比血流速度为在所述彩色血流成像的相同位置，在彩色血流成像时刻之前的若干帧图像的血流速度，所述二值化为1、0二值化。
在上述实施例的基础上，本发明实施例中提供的彩色血流成像中运动伪像的抑制方法，所述在所述彩色血流成像的相同位置，若在彩色血流成像时刻之前的若干帧图像的血流速度均为零，或仅有一帧图像的血流速度不为零，和/或所述相应的灰度图像的灰度值大于等于灰度直方图均值，则确定所述处理后血流速度为运动伪像，相应地，确定真实血流速度包括：将所述处理后血流速度及若干帧对比血流速度二值化，若二值化后的处理后血流速度为1，且二值化后的若干帧对比血流速度之和等于2，且所述相应的灰度图像的灰度值小于灰度直方图均值，则确定所述处理后血流速度为真实血流速度；其中，所述若干帧对比血流速度为在所述彩色血流成像的相同位置，在彩色血流成像时刻之前的若干帧图像的血流速度，所述二值化为1、0二值化。
在上述实施例的基础上，本发明实施例中提供的彩色血流成像中运动伪像的抑制方法，获取所述若干帧对比血流速度之和的方法，包括：将二值化后的若干帧对比血流速度进行算术累加，得到所述若干帧对比血流速度之和。具体地，可以采用程序化语言描述为：
for frame＝1:5
end
此处，若干帧对比血流速度具体是5帧，即彩色血流成像针对的血流速度的时刻前5帧的二值化(0、1二值化)结果，然后进行算术累加，得到最终结果flag。其中，C_pre表示当前血流速度之前的血流速度，baseline表示基线值。flag大于2对应于前述的情况是：将所述处理后血流速度及若干帧对比血流速度二值化，若二值化后的处理后血流速度为1，且二值化后的若干帧对比血流速度之和大于2，则确定所述处理后血流速度为真实血流速度；flag等于2对应于前述的情况是：将所述处理后血流速度及若干帧对比血流速度二值化，若二值化后的处理后血流速度为1，且二值化后的若干帧对比血流速度之和等于2，且所述相应的灰度图像的灰度值小于灰度直方图均值，则确定所述处理后血流速度为真实血流速度。
在上述实施例的基础上，本发明实施例中提供的彩色血流成像中运动伪像的抑制方法，获取所述若干帧对比血流速度之和的方法，包括：若所述若干帧对比血流速度中的一个血流速度是真实血流速度，则将计数加1，直到遍历完全部所述若干帧对比血流速度，得到的所述计数的最终值为所述若干帧对比血流速度之和。具体地，可以采用程序化语言描述为：
此处，是对总共5帧数据中的每帧数据判断每个数据点是否是真实的血流速度，若是则加1(血流速度已经0、1二值化)，直到5帧数据遍历完，得到最终的结果，即flag。其中，C_pre表示当前血流速度之前的血流速度，baseline表示基线值。flag大于2对应于前述的情况是：将所述处理后血流速度及若干帧对比血流速度二值化，若二值化后的处理后血流速度为1，且二值化后的若干帧对比血流速度之和大于2，则确定所述处理后血流速度为真实血流速度；flag等于2对应于前述的情况是：将所述处理后血流速度及若干帧对比血流速度二值化，若二值化后的处理后血流速度为1，且二值化后的若干帧对比血流速度之和等于2，且所述相应的灰度图像的灰度值小于灰度直方图均值，则确定所述处理后血流速度为真实血流速度。
为了将本发明各个实施例的技术效果更加清晰的展示出来，下面将本发明的实施例中的技术方案在现实中的实施情况与现有技术中的实施情况(例如，有无运动伪像抑制，在哪些地方抑制等)做出对比。需要说明的是，这些对比仅仅是为了将本发明各个实施例的精神本质进行阐述，并不是对本发明做出了限制，凡是符合本发明精神实质的技术方案，均在本专利的保护范围之内。参见图2，图2中展示了现有技术中彩色血流成像系统的情况，步骤201的超声脉冲发射和接收，向待测对象发射超声脉冲并从待测对象接收超声回波信号。步骤202的正交解调，包括低通滤波和降采样，获得不同深度的距离单元的正交多普勒信号。正交解调后的信号按照深度进行存储，多普勒处理需要在同一扫描线上多次重复发射脉冲的回波上进行处理。步骤203的壁滤波，也叫杂波滤波器，是一个高通滤波器，滤除静止或者缓慢运动的组织回波。步骤204的自相关，是对壁滤波后的多普勒信号进行自相关估计，获得对应的血流速度、方差和能量。步骤205的DSC，是对自相关得到的血流速度数据或者是方差和能量进行数字扫描转换，得到叠加之前用于显示的血流数据C。步骤206的B数据，也就是把经过波束合成后的信号进行脉冲回波包络检测、对数压缩等处理，得到组织灰度数据。步骤207的DSC，是灰度图像数字扫描转换，用于把步骤206得到的数据转换成用于显示的灰度图像数据B。步骤208是B和C的叠加，叠加一般采用优先权编码方式来进行。优先权编码是一组二值数据，为1则显示血流速度数据C，否则显示灰度图像数据B。步骤209的显示器，显示B和C叠加的结果。
与上述现有技术中超声血流多普勒系统相比，本发明中的各个实施例的彩色血流成像系统的实施情况可以参见图3及图4。图3中301至309各个步骤的功能与图2中的现有技术相同，在此不再赘述。与现有技术不同的是，在位于彩色DSC(305)之前增加了运动伪像抑制(3045)。图4中401至409各个步骤的功能与图2中的现有技术相同，在此不再赘述。与现有技术不同的是，在位于彩色DSC(405)之后增加了运动伪像抑制(4058)。
采用本发明的各个实施例的技术方案取得的技术效果可以参见图5，图5中包括：处理前效果501、处理后效果502及运动伪像去除后对比效果503。由图5中可见，采用了本发明的各个实施例的技术方案对处理前效果501进行处理后，得到了处理后效果502，对比501及502，可以明显看到将血液流速的运动伪像进行了去除，具体可参见运动伪像去除后对比效果503。
本发明实施例提供的彩色血流成像中运动伪像的抑制方法，通过采用流速阈值对血流速度进行筛选，再将经过筛选的血流速度与之前若干帧血流速度图像进行对比，可以精确定位真实血流速度及血流速度的运动伪像，进而将彩色血流成像中的血流速度的伪像去除，达到抑制彩色血流成像中运动伪像的技术效果。
本发明各个实施例的实现基础是通过具有处理器功能的设备进行程序化的处理实现的。因此在工程实际中，可以将本发明各个实施例的技术方案及其功能封装成各种模块。基于这种现实情况，在上述各实施例的基础上，本发明的实施例提供了一种彩色血流成像中运动伪像的抑制装置，该装置用于执行上述方法实施例中的彩色血流成像中运动伪像的抑制方法。参见图6，该装置包括：
血流运动伪像确定模块601，用于获取彩色血流成像中的血流速度及相应的灰度图像，将所述血流速度中小于流速阈值的血流速度剔除，得到处理后血流速度，在所述彩色血流成像的相同位置，若在彩色血流成像时刻之前的若干帧图像的血流速度均为零，和/或所述相应的灰度图像的灰度值大于等于灰度直方图均值，则确定所述处理后血流速度为运动伪像；
血流运动伪像去除模块602，用于去除所述运动伪像，得到真实的血流速度。
本发明实施例提供的彩色血流成像中运动伪像的抑制装置，采用血流运动伪像确定模块及血流运动伪像去除模块，通过采用流速阈值对血流速度进行筛选，再将经过筛选的血流速度与之前若干帧血流速度图像进行对比，可以精确定位真实血流速度及血流速度的运动伪像，进而将彩色血流成像中的血流速度的伪像去除，达到抑制彩色血流成像中运动伪像的技术效果。
本发明实施例的方法是依托电子设备实现的，因此对相关的电子设备有必要做一下介绍。基于此目的，本发明的实施例提供了一种电子设备，如图7所示，该电子设备包括：至少一个处理器(processor)701、通信接口(Communications Interface)704、至少一个存储器(memory)702和通信总线703，其中，至少一个处理器701，通信接口704，至少一个存储器702通过通信总线703完成相互间的通信。至少一个处理器701可以调用至少一个存储器702中的逻辑指令，以执行如下方法：获取彩色血流成像中的血流速度及相应的灰度图像，将所述血流速度中小于流速阈值的血流速度剔除，得到处理后血流速度，在所述彩色血流成像的相同位置，若在彩色血流成像时刻之前的若干帧图像的血流速度均为零，或仅有一帧图像的血流速度不为零，和/或所述相应的灰度图像的灰度值大于等于灰度直方图均值，则确定所述处理后血流速度为运动伪像；去除所述运动伪像，得到真实的血流速度。
此外，上述的至少一个存储器702中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。例如包括：获取彩色血流成像中的血流速度及相应的灰度图像，将所述血流速度中小于流速阈值的血流速度剔除，得到处理后血流速度，在所述彩色血流成像的相同位置，若在彩色血流成像时刻之前的若干帧图像的血流速度均为零，或仅有一帧图像的血流速度不为零，和/或所述相应的灰度图像的灰度值大于等于灰度直方图均值，则确定所述处理后血流速度为运动伪像；去除所述运动伪像，得到真实的血流速度。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下，即可以理解并实施。
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件。基于这样的理解，上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (8)
1.一种彩色血流成像中运动伪像的抑制方法，其特征在于，包括：
获取彩色血流成像中的血流速度及相应的灰度图像，将所述血流速度中小于流速阈值的血流速度剔除，得到处理后血流速度，在所述彩色血流成像的相同位置，若在彩色血流成像时刻之前的若干帧图像的血流速度均为零，或仅有一帧图像的血流速度不为零，和/或所述相应的灰度图像的灰度值大于等于灰度直方图均值，则确定所述处理后血流速度为运动伪像；
去除所述运动伪像，得到真实的血流速度；
所述灰度直方图均值，包括：
A＝[0,1,...,255]；f＝[f1,f2,f3,...,f256]
其中，hist_avg为灰度直方图均值，A为灰度数据，A(i)为A的第i个元素，f为灰度出现的次数，fi为f的第i个元素。
2.根据权利要求1所述的彩色血流成像中运动伪像的抑制方法，其特征在于，所述在所述彩色血流成像的相同位置，若在彩色血流成像时刻之前的若干帧图像的血流速度均为零，或仅有一帧图像的血流速度不为零，和/或所述相应的灰度图像的灰度值大于等于灰度直方图均值，则确定所述处理后血流速度为运动伪像，相应地，确定真实血流速度包括：
将所述处理后血流速度及若干帧对比血流速度二值化，若二值化后的处理后血流速度为1，且二值化后的若干帧对比血流速度之和大于2，则确定所述处理后血流速度为真实血流速度；
其中，所述若干帧对比血流速度为在所述彩色血流成像的相同位置，在彩色血流成像时刻之前的若干帧图像的血流速度，所述二值化为1、0二值化。
3.根据权利要求1所述的彩色血流成像中运动伪像的抑制方法，其特征在于，所述在所述彩色血流成像的相同位置，若在彩色血流成像时刻之前的若干帧图像的血流速度均为零，或仅有一帧图像的血流速度不为零，和/或所述相应的灰度图像的灰度值大于等于灰度直方图均值，则确定所述处理后血流速度为运动伪像，相应地，确定真实血流速度包括：
将所述处理后血流速度及若干帧对比血流速度二值化，若二值化后的处理后血流速度为1，且二值化后的若干帧对比血流速度之和等于2，且所述相应的灰度图像的灰度值小于灰度直方图均值，则确定所述处理后血流速度为真实血流速度；
其中，所述若干帧对比血流速度为在所述彩色血流成像的相同位置，在彩色血流成像时刻之前的若干帧图像的血流速度，所述二值化为1、0二值化。
4.根据权利要求2或3所述的彩色血流成像中运动伪像的抑制方法，其特征在于，获取所述若干帧对比血流速度之和的方法，包括：
将二值化后的若干帧对比血流速度进行算术累加，得到所述若干帧对比血流速度之和。
5.根据权利要求2或3所述的彩色血流成像中运动伪像的抑制方法，其特征在于，获取所述若干帧对比血流速度之和的方法，包括：
若所述若干帧对比血流速度中的一个血流速度是真实血流速度，则将计数加1，直到遍历完全部所述若干帧对比血流速度，得到的所述计数的最终值为所述若干帧对比血流速度之和。
6.一种彩色血流成像中运动伪像的抑制装置，其特征在于，包括：
血流运动伪像确定模块，用于获取彩色血流成像中的血流速度及相应的灰度图像，将所述血流速度中小于流速阈值的血流速度剔除，得到处理后血流速度，在所述彩色血流成像的相同位置，若在彩色血流成像时刻之前的若干帧图像的血流速度均为零，和/或所述相应的灰度图像的灰度值大于等于灰度直方图均值，则确定所述处理后血流速度为运动伪像；
血流运动伪像去除模块，用于去除所述运动伪像，得到真实的血流速度；
所述灰度直方图均值，包括：
A＝[0,1,...,255]；f＝[f1,f2,f3,...,f256]
其中，hist_avg为灰度直方图均值，A为灰度数据，A(i)为A的第i个元素，f为灰度出现的次数，fi为f的第i个元素。
7.一种电子设备，其特征在于，包括：
至少一个处理器、至少一个存储器、通信接口和总线；其中，
所述处理器、存储器、通信接口通过所述总线完成相互间的通信；
所述存储器存储有可被所述处理器执行的程序指令，所述处理器调用所述程序指令，以执行如权利要求1至5任一项所述的方法。
8.一种非暂态计算机可读存储介质，其特征在于，所述非暂态计算机可读存储介质存储计算机指令，所述计算机指令使所述计算机执行如权利要求1至5中任一项所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811444026.XA CN109363722B (zh) | 2018-11-29 | 2018-11-29 | 彩色血流成像中运动伪像的抑制方法及设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811444026.XA CN109363722B (zh) | 2018-11-29 | 2018-11-29 | 彩色血流成像中运动伪像的抑制方法及设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109363722A CN109363722A (zh) | 2019-02-22 |
| CN109363722B true CN109363722B (zh) | 2021-08-06 |
Family
ID=65375075
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811444026.XA Active CN109363722B (zh) | 2018-11-29 | 2018-11-29 | 彩色血流成像中运动伪像的抑制方法及设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109363722B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110353727B (zh) * | 2019-06-03 | 2021-11-23 | 浙江聚康生物工程有限公司 | 无线掌上彩色成像设备及成像方法 |
| CN115919362B (zh) * | 2023-03-15 | 2023-05-30 | 深圳英美达医疗技术有限公司 | 超声成像系统的伪像去除方法、装置、设备及存储介质 |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1939223A (zh) * | 2005-09-29 | 2007-04-04 | 深圳迈瑞生物医疗电子股份有限公司 | 消除彩色血流图像中速度异常点的方法 |
| CN101371794A (zh) * | 2008-08-05 | 2009-02-25 | 深圳市蓝韵实业有限公司 | 一种彩色多普勒超声诊断仪编码装置和方法 |
| CN101416885A (zh) * | 2007-10-22 | 2009-04-29 | 深圳迈瑞生物医疗电子股份有限公司 | 用于超声成像系统的帧平均电路 |
| CN102113900A (zh) * | 2010-01-05 | 2011-07-06 | 深圳迈瑞生物医疗电子股份有限公司 | 彩色血流动态帧相关方法和装置 |
| CN102266239A (zh) * | 2010-06-04 | 2011-12-07 | 株式会社东芝 | 超声波诊断装置 |
| CN102429684A (zh) * | 2010-09-28 | 2012-05-02 | 深圳迈瑞生物医疗电子股份有限公司 | 一种多普勒彩色血流成像方法和装置 |
| CN103099642A (zh) * | 2013-02-04 | 2013-05-15 | 左旺孟 | 一种超声血流信号质量实时分析方法 |
| CN103340621A (zh) * | 2013-06-04 | 2013-10-09 | 中国科学院苏州生物医学工程技术研究所 | 运动伪影消除装置及其消除运动伪影的方法 |
| CN103761716A (zh) * | 2014-01-25 | 2014-04-30 | 深圳市恩普电子技术有限公司 | 超声系统中的自适应伪差去除方法及系统 |
| CN105073020A (zh) * | 2013-03-12 | 2015-11-18 | 株式会社东芝 | 超声波诊断装置及超声波图像处理方法 |
| CN105122299A (zh) * | 2013-02-13 | 2015-12-02 | Ge医疗系统环球技术有限公司 | 超声图像显示装置和用于显示超声图像的方法 |
| CN105640592A (zh) * | 2016-03-18 | 2016-06-08 | 深圳安盛生物医疗技术有限公司 | 一种彩色血流帧相关计算方法 |
| US9420997B2 (en) * | 2012-06-14 | 2016-08-23 | Siemens Medical Solutions Usa, Inc. | Motion artifact suppression in ultrasound diagnostic imaging |
-
2018
- 2018-11-29 CN CN201811444026.XA patent/CN109363722B/zh active Active
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1939223A (zh) * | 2005-09-29 | 2007-04-04 | 深圳迈瑞生物医疗电子股份有限公司 | 消除彩色血流图像中速度异常点的方法 |
| CN101416885A (zh) * | 2007-10-22 | 2009-04-29 | 深圳迈瑞生物医疗电子股份有限公司 | 用于超声成像系统的帧平均电路 |
| CN101371794A (zh) * | 2008-08-05 | 2009-02-25 | 深圳市蓝韵实业有限公司 | 一种彩色多普勒超声诊断仪编码装置和方法 |
| CN102113900A (zh) * | 2010-01-05 | 2011-07-06 | 深圳迈瑞生物医疗电子股份有限公司 | 彩色血流动态帧相关方法和装置 |
| CN102266239A (zh) * | 2010-06-04 | 2011-12-07 | 株式会社东芝 | 超声波诊断装置 |
| CN102429684A (zh) * | 2010-09-28 | 2012-05-02 | 深圳迈瑞生物医疗电子股份有限公司 | 一种多普勒彩色血流成像方法和装置 |
| US9420997B2 (en) * | 2012-06-14 | 2016-08-23 | Siemens Medical Solutions Usa, Inc. | Motion artifact suppression in ultrasound diagnostic imaging |
| CN103099642A (zh) * | 2013-02-04 | 2013-05-15 | 左旺孟 | 一种超声血流信号质量实时分析方法 |
| CN105122299A (zh) * | 2013-02-13 | 2015-12-02 | Ge医疗系统环球技术有限公司 | 超声图像显示装置和用于显示超声图像的方法 |
| CN105073020A (zh) * | 2013-03-12 | 2015-11-18 | 株式会社东芝 | 超声波诊断装置及超声波图像处理方法 |
| CN103340621A (zh) * | 2013-06-04 | 2013-10-09 | 中国科学院苏州生物医学工程技术研究所 | 运动伪影消除装置及其消除运动伪影的方法 |
| CN103761716A (zh) * | 2014-01-25 | 2014-04-30 | 深圳市恩普电子技术有限公司 | 超声系统中的自适应伪差去除方法及系统 |
| CN105640592A (zh) * | 2016-03-18 | 2016-06-08 | 深圳安盛生物医疗技术有限公司 | 一种彩色血流帧相关计算方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109363722A (zh) | 2019-02-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11589840B2 (en) | Methods for super-resolution ultrasound imaging of microvessels | |
| Aysal et al. | Rayleigh-maximum-likelihood filtering for speckle reduction of ultrasound images | |
| CN114375178B (zh) | 对微血管进行高空间和时间分辨率超声成像的方法 | |
| Narayanan et al. | A view on despeckling in ultrasound imaging | |
| Yu et al. | Ultrasound speckle reduction by a SUSAN-controlled anisotropic diffusion method | |
| Czerwinski et al. | Line and boundary detection in speckle images | |
| JP6963774B2 (ja) | 流体流路の高解像度超音波画像の撮像方法 | |
| EP3048980B1 (en) | Method and apparatus for detecting a gas pocket using ultrasound | |
| WO2012051216A1 (en) | Direct echo particle image velocimetry flow vector mapping on ultrasound dicom images | |
| CN103202714B (zh) | 超声波诊断装置、医用图像处理装置以及医用图像处理方法 | |
| US10456115B2 (en) | Ultrasound system and clutter filtering method thereof | |
| CN114072068A (zh) | 利用分离的超声数据子集进行超分辨率微血管成像 | |
| CN111265246B (zh) | 超声彩色成像处理方法及装置 | |
| EP4059442B1 (en) | Image processing method and apparatus based on ultrasound contrast image | |
| CN109363722B (zh) | 彩色血流成像中运动伪像的抑制方法及设备 | |
| JP2020517363A (ja) | 改善されたクラッタ抑制を含むパワードップラー・イメージング・システム及び方法 | |
| WO2018000359A1 (zh) | 一种增强超声造影图像的方法、系统及超声造影成像设备 | |
| US8500646B2 (en) | Color Doppler mode image processing in an ultrasound system | |
| Tay et al. | A wavelet thresholding method to reduce ultrasound artifacts | |
| Lee et al. | Probability Mass Function-Based Adaptive Median Filtering for Acoustic Radiation Force Impulse Imaging: A Feasibility Study | |
| US20060030777A1 (en) | T-statistic method for suppressing artifacts in blood vessel ultrasonic imaging | |
| Zhang et al. | Evaluation of segmentation algorithms for vessel wall detection in echo particle image velocimetry | |
| Shiney et al. | An improved speckle noise reduction scheme using switching and flagging of noisy data for pre-processing of ultrasonograms in detection of down syndrome during first and second trimesters | |
| Randhawa et al. | Investigation of performance of Savitzky-Golay filter for speckle reduction in ultrasound images | |
| Tay et al. | A transform method to remove ultrasound artifacts |
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 |