彩色血流成像的壁滤波方法

摘要

本申请涉及一种彩色血流成像的壁滤波方法，属于图像处理技术领域，该方法包括：使用同相分量I和正交分量Q信号进行自相关运算得到组织运动信息；基于组织运动信息确定出当前位置符合相位旋转条件的情况下对IQ信号进行相位旋转处理；使用第i等级的壁滤波器对IQ信号进行壁滤波处理，得到滤波后的IQ信号；对滤波后的IQ信号进行自相关运算，得到血流运动信息；基于血流运动信息确定出继续进行壁滤波处理的情况下，令i＝i+1再次执行壁滤波处理的步骤；在不继续进行壁滤波处理的情况下对滤波后的自相关运算结果进行反向相位旋转处理；可以提升超声彩色多普勒血流成像中的壁滤波效率和血流图像性能。

说明书

【技术领域】

本申请涉及一种彩色血流成像的壁滤波方法，属于图像处理技术领域。

【背景技术】

彩色血流图像的获取原理包括：超声设备上探头发射的波束由被检部位反射产生频移后又被探头接收，经过接收放大。波束合成以后对回波信号进行正交解调，解调后的信号经壁滤波、自相关估计和彩色编码得到彩色血流图像。

一般地，来自人体的回波信号包括：来自血流的血流回波信号、以及来自低速运动的血管壁及组织产生的壁回波信号。其中，血流回波信号是彩色血流成像需要的信号，而壁回波信号是成像的干扰信号，一般使用壁滤波器来过滤。

如果壁滤波太弱，则无法有效滤除杂波并呈现血流；如果壁滤波太强，则低速低能量血流可能被滤除，导致血流灵敏度下降。

现有的壁滤波方法通常需要用户调节壁滤波等级，通过图像确认最合适的等级，然后保存为预设。

然而，用户调节壁滤波等级的方式比较主观，并且在不同的应用场合下需要设置不同的壁滤波等级，无法随着不同的应用和组织特征自适应变化，调节效率较低。

【发明内容】

本申请提供了一种彩色血流成像的壁滤波方法，可以解决现有的壁滤波方式无法自适应地壁滤波的问题。本申请提供如下技术方案：

提供一种彩色血流成像的壁滤波方法，所述方法包括：

获取当前回波信号的同相分量I信号和正交分量Q信号；

使用所述同相分量I信号和所述正交分量Q信号进行自相关运算，得到组织运动信息；

基于所述组织运动信息确定当前位置是否符合相位旋转条件；

在确定出符合相位旋转条件的情况下，对所述同相分量I信号和所述正交分量Q信号进行相位旋转处理，得到旋转后的信号；

使用第i等级的壁滤波器对IQ信号进行壁滤波处理，得到滤波后的IQ信号；所述IQ信号为所述同相分量I信号和正交分量Q信号、或者为所述旋转后的信号；所述壁滤波器包括m个等级，所述壁滤波器的等级按照截止频率或阶数由低到高的顺序设置，所述i依次取1至m的整数；所述m为大于1的整数；

对所述滤波后的IQ信号进行自相关运算，得到血流运动信息；

基于所述血流运动信息确定是否继续进行壁滤波处理；

在确定出继续进行壁滤波处理的情况下，令i＝i+1，并再次执行所述使用第i等级的壁滤波器对IQ信号进行壁滤波处理，得到滤波后的IQ信号的步骤；

在确定出不继续进行壁滤波处理的情况下，对所述滤波后的IQ信号的自相关结果进行反向相位旋转处理，得到最终的血流运动信息。

可选地，所述组织运动信息包括组织运动的能量、频移和方差；

所述基于所述组织运动信息确定当前位置是否符合相位旋转条件，包括：

在所述组织运动的能量、频移和/或方差满足第一阈值要求的情况下，确定所述当前位置符合相位旋转条件；在所述组织运动的能量、频移和/或方差不满足第一阈值要求的情况下，确定所述当前位置不符合相位旋转条件。

可选地，所述第一阈值要求包括：所述能量大于第一能量阈值、和/或所述频率小于第一频率阈值、和/或所述方差小于第一方差阈值。

可选地，在确定出符合相位旋转条件的情况下，所述对所述同相分量I信号和所述正交分量Q信号进行相位旋转处理，得到旋转后的信号，通过下式表示：

φ＝wT＝2πf

I'(n)+jQ'(n)＝(I(n)+jQ(n))*e

＝(I(n)+jQ(n))*(cos(nφ)-jsin(nφ))

＝(I(n)*cos(nφ)+Q(n)*sin(nφ))+j(Q(n)*cos(nφ)-I(n)sin(nφ))

其中，I表示所述同相分量I信号，Q表示所述正交分量Q信号，I’表示旋转后的同相分量I信号，Q’表示旋转后的正交分量Q信号，所述旋转后的信号包括所述旋转后的同相分量I信号和所述旋转后的正交分量Q信号，f

可选地，所述血流运动信息包括对所述滤波后的IQ信号进行自相关运算得到的能量、频移和方差；

所述基于所述血流运动信息确定是否继续进行壁滤波处理，包括：

在自相关运算的自相关结果中的的能量、频移和/或方差满足第二阈值要求的情况下，确定继续进行壁滤波处理；在所述自相关运算的自相关结果中的能量、频移和/或方差不满足第二阈值要求的情况下，确定不继续进行壁滤波处理。

可选地，所述能量大于第二能量阈值、和/或所述频移小于第二频移阈值、和/或所述方差小于第二方差阈值。

可选地，所述对所述滤波后的IQ信号进行反向相位旋转处理，得到分量输出信号，通过下式表示：

φ＝wT＝2πf

D'+jN'＝(D+jN)*e

＝(D+jN)*(cos(φ)+jsin(φ))

＝(D*cos(φ)-N*sin(φ))+j(N*cos(φ)+Dsin(φ))

其中，D和N表示所述组织运动信息中的自相关结果；D’和N’表示反向相位旋转处理后的自相关结果；对于进行相位旋转处理的IQ信号，f

可选地，所述自相关运算通过下式表示：

其中，I表示同相分量；Q表示正交分量；PS是对同一扫查线重复扫查的次数；N和D表示自相关结果；T是脉冲重复频率的倒数，

本申请的有益效果至少包括：能够根据局部组织信号特征自适应的选择合适的壁滤波等级，可以提升超声彩色多普勒血流成像中的壁滤波效率和血流图像性能，更有效地抑制血流成像中的低频杂波信号和闪彩噪声，而无需过多主观的壁滤波调节即可达到较优的血流图像。

另外，通过基于组织运动信息确定当前位置是否符合相位旋转条件，对符合相位旋转条件的部分进行相位旋转处理后再进行壁滤波处理，可以让组织的运动近似归零，从而提高壁滤波的效率，有效滤除低频组织运动，减少闪彩伪像。

另外，通过在确定出停止壁滤波后，对当前壁滤波后的自相关结果进行反向相位旋转操作，可以弥补此前因I/Q相位旋转操作去掉的频移或相移，提高分量输出信号的准确性。

上述说明仅是本申请技术方案的概述，为了能够更清楚了解本申请的技术手段，并可依照说明书的内容予以实施，以下以本申请的较佳实施例并配合附图详细说明如后。

【附图说明】

图1是本申请一个实施例提供的彩色血流成像的壁滤波方法的流程图；

图2是本申请一个实施例提供的彩色血流成像过程的流程图；

图3是本申请一个实施例提供的彩色血流成像的壁滤波装置的框图。

【具体实施方式】

下面结合附图和实施例，对本申请的具体实施方式做进一步详细描述。以下实施例用于说明本申请，但不用来限制本申请的范围。

首先，对本申请涉及的若干名词进行介绍。

有限长单位冲激响应(Finite Impulse Response，FIR)滤波器：可以在保证任意幅频特性的同时具有严格的线性相频特性，同时其单位抽样响应是有限长的。

用于血流成像的回波信号经过FIR滤波器的处理后，数据的有效点数有一定的损失，损失的点数为滤波器阶数。当滤波器阶数较低时，过渡带太宽，滤波效果不理想；滤波器阶数较高时(不能超过数据长度)，数据损失点数太多，会增加血流参数估计的误差。

无限脉冲响应(Infinite Impulse Response，IIR)数字滤波器：采用递归型结构，即结构上带有反馈环路。IIR滤波器运算结构通常由延时、乘以系数和相加等基本运算组成。

IIR滤波器的特点是过渡带较窄，在较低的阶数下即可获得较好的滤波效果。但是，IIR滤波器暂态响应较长，因此，在回波信号较短的数据长度下，暂态响应将使血流参数的估计产生较大偏差，直接影响滤波器的性能。决定滤波器暂态响应的除了IIR滤波器的种类，更重要的是滤波器的初始化方法。目前主要有三种初始化方法：

zero-initialized IIR filters：即在寄存器中载入初始化状态矢量为零。

step-initialized IIR filters：即通过假设第一个复采样值在n为无穷时已经存在来设置初始化状态矢量。

projection initialization：即利用状态空间法使暂态响应最小化。首先把滤波器的输出分解为暂态分量和稳态分量，然后使由初始状态产生的暂态响应的投影的幅值最小化。

理想的滤波器设计完成后，将其滤波器系数输入到设计好的初始化程序中，即可分别得到不同初始化条件下的IIR滤波器。IIR滤波器比较容易设计，而且截止频率可以连续动态变化。三种形式的IIR滤波器中，投影初始化方法(projection initialization)的暂态响应最小，不损失数据点数，效果最好。

回归(Regression)滤波器：基本原理是把输入信号作为时域的一个多项式函数，多普勒信号中低频的杂波成分可以用一个给定阶次的多项式来近似拟合(通常采用最小二乘拟合)，然后从多普勒信号中减掉这部分信号，进而提取出血流信号，该滤波器的实质是一个高通滤波器。

其基本原理如下：

假设以x(n)作为输入信号，y(n)为输出信号，c(n)为拟合的低频杂波信号，a

y(n)＝x(n)-c(n)。

一种典型的回归滤波器包括以勒让德多项式作为基向量构建的回归滤波器。勒让德多项式可以通过对多项式{1,n

回归滤波器具有平滑、单调的频率响应，阻带衰减大、过渡带窄，而且不损失数据点数的特点，因此非常适合作为多普勒血流信号的壁滤波器。从回归滤波器的基本原理和频率响应函数可以看出回归滤波器的截止频率取决于数据长度N和多项式最大阶数K。实际应用中，为了兼顾实时性和计算精度，数据长度通常取8～24之间的整数，而在数据长度一定的情况下，滤波器截止频率则完全取决于多项式最大阶数，而多项式阶数只能取到1、2、3、4之类的整数(阶数过大截止频率太高会过滤掉血流信号)，因此无法设计任意截止频率的回归滤波器，这是回归滤波器的一大缺陷。

上述滤波器仅是示意性地，在实际实现时，用于对多普勒血流信号进行滤波的壁滤波器，还有基于特征值分解或奇异值分解的壁滤波方法，这类方法是对信号进行特征值分解或奇异值分解，然后滤除其中的静态成分和低频成分，本申请在此不再一一列举。

传统的超声彩色多普勒血流成像的壁滤波方法，包括：获取输入信号；根据输入信号数据长度获得回归滤波器各阶次的滤波器系数矩阵；将相邻阶次的滤波器系数矩阵进行线性组合构建新的滤波器系数矩阵；根据新的滤波器系数矩阵对输入信号进行滤波得到输出信号。

但是，上述方法仅适用于回归滤波器，依然无法实现自适应地对多普勒血流信号进行壁滤波。

另外，传统的超声彩色多普勒血流成像的壁滤波方法，还包括：获取彩色回波信号，根据所述彩色回波信号得到每个采样位置的彩色特征数据；根据所述彩色特征数据进行DSC数字扫描变换，得到图像显示区域的每个像素点的彩色图像数据；确定满足预设阈值条件的像素点的数据为血流数据，所述预设阈值条件根据所述彩色图像数据确定；根据所述血流数据得到对应像素点的血流速度数据，根据所述血流速度数据对图像进行彩色编码成像。

然而，上述方法也无法自适应地实现壁滤波。

基于传统的壁滤波方式无法实现自适应滤波的问题，本申请提供一种彩色血流成像的壁滤波方法，主要目的是根据组织信号特征实现自适应的壁滤波，从而提高血流成像的灵敏度和表现力，并减轻医务人员的调节负担。

下面对本申请提供的彩色血流成像的壁滤波方法进行详细介绍，本申请以各个实施例提供的彩色血流成像的壁滤波方法用于电子设备中为例进行说明，该电子设备为超声设备、或者与超声设备通信相连的终端或服务器，终端可以为计算机、或者平板电脑等，本实施例不对电子设备的类型作限定。

图1是本申请一个实施例提供的彩色血流成像的壁滤波方法的流程图，该方法至少包括以下几个步骤：

步骤101，获取当前回波信号的同相分量I信号和正交分量Q信号。

对当前回波信号进行正交解调后得到该当前回波信号的同相分量I信号和正交分量Q信号。

其中，当前回波信号用于进行彩色血流成像，当前回波信号也称为多普勒血流信号或者超声回波信号等，本实施例不对当前回波信号的名称作限定。

步骤102，使用同相分量I信号和正交分量Q信号进行自相关运算，得到组织运动信息。

自相关运算通过下式表示：

其中，I表示同相分量(包括本步骤中的同相分量I信号)；Q表示正交分量(包括本步骤中的正交分量Q信号)；PS是对同一扫查线重复扫查的次数；N和D表示自相关结果；T是脉冲重复频率的倒数，

示意性地，为了保证一定的帧频和计算精度，PS的取值范围通常在8～24之间。

根据上述内容可知，将同相分量I信号和正交分量Q信号输入上述自相关运算公式后，得到的组织运动信息包括：组织运动的自相关结果N和D、能量R(0)、频移和方差。

步骤103，基于组织运动信息确定当前位置是否符合相位旋转条件。

示意性地，基于组织运动信息确定当前位置是否符合相位旋转条件，包括：在组织运动的能量、频移和/或方差满足第一阈值要求的情况下，确定当前位置符合相位旋转条件。在组织运动的能量、频移和/或方差不满足第一阈值要求的情况下，确定当前位置不符合相位旋转条件。

由于一般能量较大，频移和方差较小的回波信号会被初步判定为组织，需要做相位旋转(rotation)处理，以便将组织运动的频移归零。而能量较小，频移和方差较大的回波信号会被初步判定为血流，不符合相位旋转条件，不需要做rotation处理，以防止正常血流信号被丢失。基于此，第一阈值要求包括能量大于第一能量阈值、和/或所述频率小于第一频率阈值、和/或所述方差小于第一方差阈值。相应地，判断逻辑通过下式表示：

Rw>prev_r和/或fw

Rwprev_f和/或δw>prev_δ，则确定出不符合相位旋转条件，而是组织执行步骤105。

其中，Rw表示能量、fw表示频移、δw表示方差；prev_r表示第一能量阈值、prev_f表示第一频移阈值、prev_δ表示第一方差阈值。

步骤104，在确定出符合相位旋转条件的情况下，对同相分量I信号和正交分量Q信号进行相位旋转处理，得到旋转后的信号。

本实施例中，对不是血流的回波信号进行相位旋转处理的目的为：让组织的运动近似归零，从而提高壁滤波的效率，有效滤除低频组织运动，减少闪彩伪像。

在确定出符合相位旋转条件的情况下，对同相分量I信号和正交分量Q信号进行相位旋转处理，得到旋转后的信号，通过下式表示：

φ＝wT＝2πf

I'(n)+jQ'(n)＝(I(n)+jQ(n))*e

＝(I(n)+jQ(n))*(cos(nφ)-jsin(nφ))

＝(I(n)*cos(nφ)+Q(n)*sin(nφ))+j(Q(n)*cos(nφ)-I(n)sin(nφ))

其中，I表示同相分量I信号，Q表示正交分量Q信号，I’表示旋转后的同相分量I信号，Q’表示旋转后的正交分量Q信号，旋转后的信号包括旋转后的同相分量I信号和旋转后的正交分量Q信号，f

经过上述相位旋转处理可知，回波信号中原有的平均频移f

步骤105，使用第i等级的壁滤波器对IQ信号进行壁滤波处理，得到滤波后的IQ信号；IQ信号为同相分量I信号和正交分量Q信号、或者为旋转后的信号；壁滤波器包括m个等级，壁滤波器的等级按照截止频率或阶数由低到高的顺序设置，i依次取1至m的整数。

其中，m为大于1的整数。

在当前位置不符合相位旋转条件的情况下IQ信号为同相分量I信号和正交分量Q信号；在当前位置为组织的情况下，IQ信号为或者为旋转后的信号。

其中，m个等级中的壁滤波器中每个等级的壁滤波器由至少一种不同的壁滤波器组成。至少一种不同的壁滤波器包括以下几种：FIR滤波器，无IIR数字滤波器，回归Regression滤波器，基于特征值分解的滤波器，基于奇异值分解的滤波器。

本实施例中，随着壁滤波器等级由低到高，壁滤波器截止频率或壁滤波器阶数依次由低到高，壁滤波器属于高通滤波器，即随着壁滤波器等级由低到高，所滤除的低频成分越来越多。并且，壁滤波处理从最低等级开始调用，且不超过最高等级。

步骤106，对滤波后的IQ信号进行自相关运算，得到血流运动信息。

壁滤波后的IQ信号的自相关运算的自相关结果主要反应血流运动信息。本步骤中的自相关运算参考步骤102中的自相关运算过程，此时，公式中的I表示滤波后的同相分量I信号；Q表示滤波后的正交分量Q信号。

根据步骤102中的公式可知，血流运动信息包括对滤波后的IQ信号进行自相关运算得到的的能量、频移和方差。

步骤107，基于血流运动信息确定是否继续进行壁滤波处理。

示意性地，基于血流运动信息确定是否继续进行壁滤波处理，包括：在自相关运算的自相关结果中的能量、频移和/或方差满足第二阈值要求的情况下，确定继续进行壁滤波处理。

一般能量较大，频移和方差较小的信号点需要做更高等级的壁滤波处理，以便将组织运动的低频成分更好的滤除；而能量较小，频移和方差较大的信号点不需要做更高等级的壁滤波处理，以防止正常血流信号被滤除和增加计算量。基于此，第二阈值要求包括：能量大于第二能量阈值、和/或所述频移小于第二频移阈值、和/或所述方差小于第二方差阈值。相应地，判断逻辑通过下式表示：

R>post_r和/或f

Rpost_f和/或δ>post_δ，则不继续进行壁滤波处理，执行步骤109。

其中，第二阈值要求与壁滤波器的阻带衰减、截止频率等设计参数相关。

步骤108，在确定出继续进行壁滤波处理的情况下，令i＝i+1，并再次执行使用第i等级的壁滤波器对IQ信号进行壁滤波处理，得到滤波后的IQ信号的步骤，即再次执行步骤105。

步骤109，在确定出不继续进行壁滤波处理的情况下，对滤波后的IQ信号的自相关结果进行反向相位旋转处理，得到最终的血流运动信息。

经过步骤107的判断后，一部分信号点中断壁滤波处理；一部分信号点继续更高等级壁滤波处理和自相关运算，直到满足退出壁滤波处理的条件或壁滤波等级达到最高。退出壁滤波处理后，需要对当前壁滤波后的自相关结果进行反向相位旋转(fw)操作，目的是弥补此前因I/Q相位旋转(-fw)操作去掉的频移或相移，对于没有实施相位旋转处理的信号点，fw置为0。

具体地，对滤波后的IQ信号的自相关结果进行反向相位旋转处理，得到最终的血流运动信息，通过下式表示：

φ＝wT＝2πf

D'+jN'＝(D+jN)*e

＝(D+jN)*(cos(φ)+jsin(φ))

＝(D*cos(φ)-N*sin(φ))+j(N*cos(φ)+Dsin(φ))

其中，D和N表示组织运动信息中的自相关结果；D’和N’表示反向相位旋转处理后的自相关结果；对于进行相位旋转处理的IQ信号，f

可选地，参考图2，在进行反向相位旋转处理后，还可以包括步骤21-28：

步骤21，对血流运动信息进行血流能量对数压缩；

血流能量对数压缩的目的是将较大的动态范围压缩到0～255的显示范围；

步骤22，对压缩后的信号进行N/D/R二维空间平滑；

N/D/R二维空间平滑的目的是提高血流的平滑度和表现力；

步骤23，对平滑后的信号进行血流速度、能量、方差计算和仲裁；

血流速度、能量、方差计算和仲裁的目的是将N/D/R转换成血流速度、能量、方差并判断出哪些信号点是血流；

步骤24，对仲裁后的信号进行帧平均；

帧平均的目的是对相邻帧之间进行时间平滑，以提升血流的时间连续性和表现力；

步骤25，对帧平均后的信号进行血流边界处理；

血流边界处理的目的是消除空洞和边界锯齿；

步骤26，对边界处理后的信号进行扫描变换，以将扫描线空间转换到屏幕像素空间；

步骤27，将变换后的信号进行B/C图像混合及仲裁；

B/C图像混合及仲裁的目的是对扫描变换后的B模式和CF模式图像进行alpha混合并判断哪些像素点是血流；

步骤28，将混合后的图像显示。

具体地，混合后的RGB图像送往显示器进行最终显示。

需要补充说明的是，本申请的主要技术方案是彩色多普勒血流成像中的一种自适应的壁滤波方法，其它彩色多普勒血流成像过程不是本申请关注的重点，在此不做深入展开。

值得一提的是，在血流速度、能量、方差计算和仲裁的过程中，仲裁判断中的阈值与本申请的自适应壁滤波方法中涉及的壁滤波器阻带衰减、截止频率等设计参数直接相关，因此也是本申请自适应壁滤波方法的相关部分。

综上所述，本实施例提供的彩色血流成像的壁滤波方法，能根据局部组织信号特征自适应的选择合适的壁滤波等级，可以提升超声彩色多普勒血流成像中的壁滤波效率和血流图像性能，更有效地抑制血流成像中的低频杂波信号和闪彩噪声，而无需过多主观的壁滤波调节即可达到较优的血流图像。

另外，通过基于组织运动信息确定当前位置是否符合相位旋转条件，对符合相位旋转条件的部分进行相位旋转处理后再进行壁滤波处理，可以让组织的运动近似归零，从而提高壁滤波的效率，有效滤除低频组织运动，减少闪彩伪像。

另外，通过在确定出停止壁滤波后，对当前壁滤波后的自相关结果进行反向相位旋转操作，可以弥补此前因I/Q相位旋转操作去掉的频移或相移，提高分量输出信号的准确性。

图3是本申请一个实施例提供的彩色血流成像的壁滤波装置的框图。该装置至少包括以下几个模块：信号获取模块310、第一自相关模块320、血流判断模块330、相位旋转模块340、壁滤波模块350、第二自相关模块360、滤波判断模块370和反向旋转模块380。

信号获取模块310，用于获取当前回波信号的同相分量I信号和正交分量Q信号；

第一自相关模块320，用于使用所述同相分量I信号和所述正交分量Q信号进行自相关运算，得到组织运动信息；

血流判断模块330，用于基于所述组织运动信息确定当前位置是否符合相位旋转条件；

相位旋转模块340，用于在确定出符合相位旋转条件的情况下，对所述同相分量I信号和所述正交分量Q信号进行相位旋转处理，得到旋转后的信号；

壁滤波模块350，用于使用第i等级的壁滤波器对IQ信号进行壁滤波处理，得到滤波后的IQ信号；所述IQ信号为所述同相分量I信号和正交分量Q信号、或者为所述旋转后的信号；所述壁滤波器包括m个等级，所述壁滤波器的等级按照截止频率或阶数由低到高的顺序设置，所述i依次取1至m的整数；

第二自相关模块360，用于对所述滤波后的IQ信号进行自相关运算，得到血流运动信息；

滤波判断模块370，用于基于所述血流运动信息确定是否继续进行壁滤波处理；

所述壁滤波模块350，还用于在确定出继续进行壁滤波处理的情况下，令i＝i+1，并再次执行所述使用第i等级的壁滤波器对IQ信号进行壁滤波处理，得到滤波后的IQ信号的步骤；

反向旋转模块380，用于在确定出不继续进行壁滤波处理的情况下，对所述滤波后的IQ信号的自相关结果进行反向相位旋转处理，得到最终的血流运动信息。

相关细节参考上述方法实施例。

需要说明的是：上述实施例中提供的彩色血流成像的壁滤波装置在进行彩色血流成像的壁滤波时，仅以上述各功能模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能模块完成，即将彩色血流成像的壁滤波装置的内部结构划分成不同的功能模块，以完成以上描述的全部或者部分功能。另外，上述实施例提供的彩色血流成像的壁滤波装置与彩色血流成像的壁滤波方法实施例属于同一构思，其具体实现过程详见方法实施例，这里不再赘述。

可选地，本申请还提供有一种计算机可读存储介质，所述计算机可读存储介质中存储有程序，所述程序由处理器加载并执行以实现上述方法实施例的彩色血流成像的壁滤波方法。

可选地，本申请还提供有一种计算机产品，该计算机产品包括计算机可读存储介质，所述计算机可读存储介质中存储有程序，所述程序由处理器加载并执行以实现上述方法实施例的彩色血流成像的壁滤波方法。

以上所述实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。

以上所述实施例仅表达了本申请的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本申请构思的前提下，还可以做出若干变形和改进，这些都属于本申请的保护范围。因此，本申请专利的保护范围应以所附权利要求为准。