CN102306377A

CN102306377A - 一种超声图像去噪的方法和装置

Info

Publication number: CN102306377A
Application number: CN201110281240A
Authority: CN
Inventors: 凌锋; 齐保春
Original assignee: Edan Instruments Inc
Current assignee: Edan Instruments Inc
Priority date: 2011-09-21
Filing date: 2011-09-21
Publication date: 2012-01-04
Anticipated expiration: 2031-09-21
Also published as: CN102306377B

Abstract

本发明涉及到图像处理技术，特别涉及超声成像中的图像数据处理技术，尤其是一种超声图像去噪的方法和装置，本发明包括：读取超声图像数据，以各个像素点为中心选取一个邻域，计算邻域内像素点各个方向的方差均值比，根据所述的方差均值比计算出判别因子，根据所述的判别因子分别区分出像素点邻域为边缘区、非边缘区和半边缘区，对所述不同的边缘区分别进行滤波处理，输出处理后的超声图像数据。本发明所提供的技术方案可以同时实现对图像的边缘增强与斑点抑制，并且具有算法简单、自适应性，易于硬件实现，可实时处理，实用性强等优点。

Description

一种超声图像去噪的方法和装置

技术领域

本发明涉及到图像处理技术，特别涉及超声成像中的图像数据处理技术，尤其是一种超声图像去噪的方法和装置。

背景技术

医学超声成像中，由于成像器官或组织结构的不均匀性，一些微小的结构不能被分辨，加上声波信号的干涉现象，在超声图像上形成了特有的斑点。降低了图像质量，目前常用的去噪算法主要有分辨率增强、叠加平均以及后处理方法。

分辨率增强算法主要是增强超声系统的分辨率，达到降噪的目的。这类方法有提高探头频率、编码激励(Coded-excitation)、阵列探头(Matrix-array transducer)以及谐波(Harmonics)，但探头频率过高在实际应用中并不可行。叠加平均法对多帧去相关图像进行叠加平均，由于去相关后的各帧图像具有不同的speckle噪声分布，同时具有相同的特征信息，故对其进行叠加平均可减少speckle噪声并对特征信息进行增强。但是，叠加平均法只能有限获得speckle抑制效果(1/n，n为帧数)。而且，叠加平均法中所取的多帧图像会降低超声系统的帧率，限制了此项技术的实际应用。后处理方法一般分为四类：中值滤波、维纳滤波、基于扩散方程的滤波以及基于小波变换的滤波方法。中值滤波方法根据图像的局部统计特征来自动选取滤波窗口内的各点权值，或者自动选取滤波窗口的大小和形状(AWMF，自适应加权中值滤波)，维纳滤波方法中，计算一定邻域大小的方差和均值。它是一种自适应的滤波方法，当邻域内方差大的时候，平滑就弱一些，反之，平滑就强以一些。尽管这些方法在保留图像细节方面取得了一定的效果，但这类基于窗口平滑的滤波方法对于窗口形状和大小均十分敏感。基于小波变换的滤波方法将超声图像变换到小波域，利用小波阈值处理将某些尺寸的元素丢弃，再进行逆变换以去除Speckle。但在该方法中，很难选择合适的尺度，如果丢弃小尺度信号，则不能有效去除Speckle；如果丢弃大尺度信号，则有可能丢弃有用信号。基于扩散方程的滤波方法是近年来出现的一种新的超声图像滤波技术，即求解初始值为输入图像的非线性热扩散方程。在扩散方程中，通过引入图像特征，设计合适的扩散系数来控制扩散方程的扩散行为，使得在平滑图像的同时能够保留甚至增强图像的特征信息。

美国专利US6208763和US6592523中公开了一种基于梯度判别因子的图像去噪边缘增强算法，该算法依据图像像素点的梯度幅值大小及方向信息，把图像分为边缘区域及非边缘区域，对于边缘区域采用的是各向异性边缘增强处理，对于非边缘区域采用的是各向同性平滑处理。

中国专利CN200510035913.8中提出一种超声图像边缘锐化与斑点抑制方法，该算法依据图像像素点的梯度信息、方差信息、灰度值信息，把图像分为边缘区域及非边缘区域，对边缘点先进行方向性滤波，再做方向性增强，得到增强的边缘，对于非边缘点求出该点为中心像素点的一个邻域内的像素数据均值，以及各向同性增强，再求两者的加权平均值作为所述非边缘点的灰度均值。

上述专利方法的不足之处在于：根据一定邻域像素点的数学统计特性，简单的将像素点分成边缘区域与非边缘区域，这种分法会造成图像信息量的丢失，且计算量大，难以实际应用。

发明内容

为克服上述缺陷，本发明的目的即在于一种超声图像去噪的方法和装置。

本发明的目的是通过以下技术方案来实现的：

本发明一种超声图像去噪的方法和装置，其中所述的一种超声图像去噪的装置，主要包括：

图像数据读取模块、邻域选取模块、方差均值比计算模块、判断因子计算模块、区域判断模块、滤波处理模块和输出模块，

所述的图像数据读取模块，用于读取超声图像的数据；

所述的邻域选取模块，与所述的图像数据读取模块相连接，用于对所述读取的超声图像数据的每一个像素点，选取以其为中心的一定大小的区域为所述像素点对应的一个邻域。

所述的方差均值比计算模块，与所述的邻域选取模块相连接，用于计算所述邻域内像素点各个方向的方差均值比；

所述的判断因子计算模块，与所述的方差均值比计算模块相连接，用于根据各个方向上的方差均值比计算出判别因子；

所述的区域判断模块，与所述的判断因子计算模块相连接，用于根据判别因子分别区分出像素点邻域为边缘区、非边缘区和半边缘区；

所述的滤波处理模块，与所述的区域判断模块相连接，用于对不同的边缘区分别进行滤波处理；

所述的输出模块，与所述的滤波处理模块相连接，用于输出处理后的超声图像数据。

作为本发明进一步的技术方案，所述的一种超声图像去噪的装置还包括处理完成判断模块，与所述的邻域选取模块、滤波处理模块和输出模块相连接，用于根据所述滤波处理模块的处理结果判断是否处理完所述读取超声图像的数据中的每一个像素点，如果完成则将结果发送到所述的输出模块，如果未完成，则返回到所述的邻域选取模块继续进行像素点的处理。

作为本发明更进一步的技术方案，所述的滤波处理模块，包括了，

边缘区滤波单元，用于对所述的邻域为边缘区进行增强型方向性滤波；

半边缘区滤波单元，用于邻域为半边缘区进行方向性滤波、方向性增强并加权平均；

非边缘区滤波单元，用于邻域为非边缘区进行各向同性平滑滤波。

本发明一种超声图像去噪的方法，包括如下步骤，

步骤1，读取超声图像数据；

步骤2，以各个像素点为中心选取一个邻域；

步骤3，计算邻域内像素点各个方向的方差均值比；

步骤4，根据所述的方差均值比计算出判别因子；

步骤5，根据所述的判别因子分别区分出像素点邻域为边缘区、非边缘区和半边缘区；

步骤6，对所述不同的边缘区分别进行滤波处理；

步骤7，输出处理后的超声图像数据。

作为本发明进一步的技术方案，所述的步骤3，进一步包括，

步骤31，计算邻域内像素点四个方向的方差均值比；

步骤32，计算出各个方差均值比的最大值和最小值。

作为本发明更进一步的技术方案，所述的步骤4，进一步包括，根据所述方差均值比的最大值和最小值之差与最大值的比值，计算出所述的判别因子。

作为本发明更进一步的技术方案，所述的步骤5，进一步包括，当判别因子趋近于1时，所述邻域为图像的边缘区；当判别因子趋近于0时，所述邻域为图像的非边缘区；当判别因子在0与1之间时，所述区域为图像的半边缘区。

作为本发明更进一步的技术方案，所述的步骤6，进一步包括，

步骤61，邻域为边缘区进行增强型方向性滤波；

步骤62，邻域为半边缘区进行方向性滤波、方向性增强并加权平均；

步骤63，邻域为非边缘区进行各向同性平滑滤波。

作为本发明更进一步的技术方案，所述的步骤7之前还包括，判断是否处理完所述读取超声图像的数据中的每一个像素点，如果完成对每一个像素点的处理则进入所述的步骤7，如果没有则继续从所述步骤2开始对未处理的像素点进行相应的计算处理。

作为本发明更进一步的技术方案，所述的步骤7，进一步包括，将经过处理后的最终图像数据进行输出，包括显示和或存储和或打印。

本发明所提供的技术方案采用了对边缘区域、半边缘区域、非边缘区域的分别判别和处理，且判别因子

可根据邻域计算得到，因此，可以同时实现对图像的边缘增强与斑点抑制，并且具有算法简单、自适应性，易于硬件实现，可实时处理，实用性强等优点。

附图说明

为了易于说明，本发明由下述的较佳实施例及附图作以详细描述。

图1为超声成像系统的结构框图；

图2为本发明一种超声图像去噪的方法的一种实施例流程图；

图3为本发明一种超声图像去噪的方法的另一种实施例流程图；

图4为本发明一种超声图像去噪的装置的一种实施例结构图；

图5为本发明一种超声图像去噪的装置的另一种实施例结构图；

图6为像素点四个计算方向的示意图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

如图1所示，为典型的超声成像系统的结构框图，图中所示由主控制器、探头、波束合成、图像处理、数字扫描变换(DSC, digital scan conversion)、显示器组成。在主控制器的控制下，探头发射一定频率的超声波，经过一定时间的延迟后，接收由组织反射回来的声波信号，经过聚焦延时、加权求和，形成一条条扫描线，再经过滤波、检波等一系列的数字信号处理获得扫描回波的包络信号，形成图像数据，再把图像数据送往显示器进行显示。这其中就要涉及到斑点噪声的去除，一般的噪声去除、图像增强处理方法是在DSC之后，也就是形成图像数据后，再对图像数据进行处理，但为了得到更真实的数据，也可以位于DSC之前，即获得扫描线数据(也是超声图像数据)后，就进行处理。

本发明的一种具体实施方式如图2所示：

101.读取超声图像数据；

在数字扫描变换之前或者之后，获得超声图像的数据；

102.以各个像素点为中心选取一个邻域；

以任何一个点为中心的任何开区间称为该点的邻域，对于每一个像素点，将以其为中心的一定大小(如3×3，或者5×5)的区域选取为其对应的一个邻域；

103.计算邻域内像素点各个方向的方差均值比；

如图6所示，本实施例优选计算邻域内所述像素点四个方向上的方差与均值的比值，即为方差均值比；

104.根据方差均值比计算出判别因子；

根据上述计算得到的四个方向上的四个方差均值比，本发明优选根据其最大值和最小值之差与最大值的比值，计算出判别因子；另外，还可以有其他的方式，比如直接根据方差均值比的大小来作为判别因子判断区域，方差均值比值大时为边缘区域，相反小的时候为非边缘区域，中间值时为半边缘区域。

105.根据判别因子分别区分出像素点邻域为边缘区、非边缘区和半边缘区；

当判别因子趋近于1时，所述邻域为图像的边缘区；当判别因子趋近于0时，所述邻域为图像的非边缘区；当判别因子在0与1之间时，所述区域为图像的半边缘区；

106.对不同的边缘区分别进行滤波处理；

本发明实施例中，优选对于边缘区域像素点的处理为增强型方向性滤波，在增强边缘的同时抑制斑点噪声；优选对于半边缘区域像素点的处理为先进行方向性滤波，再采用方向性增强，并对两计算结果加权平均计算得到最终结果；优选对于非边缘区域像素点的处理是各向同性平滑滤波。

107.输出处理后的超声图像数据；

将经过上述处理后的最终图像数据进行输出，包括显示和或存储和或打印；

为了更好的理解本发明一种超声图像去噪的方法，本发明的另一种实施方式如图3所示：

101.读取超声图像数据；

在数字扫描变换之前或者之后，获得超声图像的数据；

102.以各个像素点为中心选取一个邻域；

1031.计算邻域内像素点四个方向的方差均值比；

计算得出方差均值比

分别为：，

，

Figure 201110281240X100002DEST_PATH_IMAGE005

，；

1032.计算出各个方差均值比的最大值和最小值；

得到最大值

Figure 201110281240X100002DEST_PATH_IMAGE007

和最小值

；

以上的1031和1032两步完成了对邻域内像素点各个方向的方差均值比的整个计算过程。

104.根据方差均值比的最大值和最小值计算出判别因子λ；

在本步骤中，优选的计算方式为： ；

105.根据判别因子λ分别区分出像素点邻域为边缘区、非边缘区和半边缘区；

1061.邻域为边缘区进行增强型方向性滤波；

以3×3大小的邻域来说明其计算过程：

对每一个中心像素点，在邻域内计算四个方向上(方向如图6所示)的均值，注意此均值的计算不包括中心像素点，即计算四个方向上两个像素点的均值；

再分别计算中心像素点的灰度值(灰度值是计算机图像中的术语，用来描述生成的图像所能包含的颜色数)与四个方向上均值的差值，比较得到最小差值的方向；

将最小差值方向上的均值赋给中心像素点。

1062.邻域为半边缘区进行方向性滤波、方向性增强并加权平均；

以3×3大小的邻域来说明其计算过程：

对每一个中心像素点，在邻域内计算四个方向上(方向如图6所示)的方差，注意，在计算方差时所要进行的均值计算包括中心像素点，即涉及到每个方向上三个像素点的计算。

再比较得到最小方差的方向，并将最小方差方向上的均值赋给中心像素点，此时中心像素点的灰度值为。

为进一步增强边缘，系统对于半边缘区域像素点还包括方向性增强处理，具体过程如图7所示，以3×3大小的邻域来说明其计算过程：

在邻域内计算四个方向上(方向如图6所示)的一维laplacian值，可以设置一维laplacian模板为{-1,2,-1}，比较得到四个方向上最大的laplacian幅值

Figure 201110281240X100002DEST_PATH_IMAGE011

，并将

赋给中心像素点。

中心像素点处理后的最终灰度值

：

Figure 201110281240X100002DEST_PATH_IMAGE013

也就是说，当

趋近于1时，对所述邻域进行边缘增强；当

趋近于0时，对所述邻域进行方向性平滑滤波；当

在0与1之间时，对所述区域同时进行边缘增强和方向性平滑滤波，

为前面计算出来的区域判别因子，

为、

的加权平均，也就是说对于半边缘区域像素点的处理为：先计算方向性滤波值

，再方向性增强值

，并对这两个值求加权平均，得到最终的结果(

)，加权因子为

(如上式所示)。

上述中的一维laplacian计算方法为，以水平方向上为例，设置一维laplacian模板为{-1,2,-1}，中心像素点为,三个像素点为

Figure 201110281240X100002DEST_PATH_IMAGE015

为例，计算为：

1063.邻域为非边缘区进行各向同性平滑滤波；

所述的各向同性平滑滤波处理，以3×3大小的邻域来说明其计算过程：计算邻域9个像素点的均值，并将该均值赋给中心像素点，另外各向同性滤波可以包括很多的算法，高斯平滑，均值滤波等。

如上所述的步骤1061、1062、1063为对不同的边缘区分别进行滤波处理的详细步骤，当然，进行边缘滤波处理的方法很多，还包括业界悉知的小波处理、高斯-拉普拉斯处理等等。

108.判断是否处理完每一个像素点；

如果完成对每一个像素点的处理则进入下一步骤，如果没有则继续从所述步骤102开始对未处理的像素点进行相应的计算处理。

107.输出处理后的超声图像数据；

为了更好的理解本发明一种超声图像去噪的装置，本发明的装置的一种实施方式如图4所示：

图像数据读取模块301；

用于读取超声图像的数据。

邻域选取模块302；

与所述的图像数据读取模块301相连接，用于对于每一个像素点，将以其为中心的一定大小(如3×3，或者5×5)的区域选取为其对应的一个邻域。

方差均值比计算模块303；

与所述的邻域选取模块302相连接，用于计算所述邻域内像素点各个方向的方差均值比。

判断因子计算模块304；

与所述的方差均值比计算模块303相连接，用于根据各个方向上的方差均值比的最大值和最小值之差与最大值的比值，计算出判别因子。

区域判断模块305；

与所述的判断因子计算模块304相连接，用于根据判别因子分别区分出像素点邻域为边缘区、非边缘区和半边缘区。

滤波处理模块306；

与所述的区域判断模块305相连接，用于对不同的边缘区分别进行滤波处理。

输出模块307；

与所述的滤波处理模块306相连接，用于输出处理后的超声图像数据。

为了更好的理解本发明一种超声图像去噪的装置，本发明的装置的另一种实施方式如图5所示，还包括了：

处理完成判断模块308；

与所述的邻域选取模块302、滤波处理模块306和输出模块307相连接，用于根据所述滤波处理模块306的处理结果判断是否处理完每一个像素点，如果完成则将结果发送到所述的输出模块307，如果未完成，则返回到所述的邻域选取模块302继续进行像素点的处理。

所述的滤波处理模块306，进一步包括：

边缘区滤波单元309，用于对所述的邻域为边缘区进行增强型方向性滤波；

半边缘区滤波单元310，用于邻域为半边缘区进行方向性滤波、方向性增强并加权平均；

非边缘区滤波单元311，用于邻域为非边缘区进行各向同性平滑滤波；

以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims

1.一种超声图像去噪的装置，其特征在于，包括，图像数据读取模块、邻域选取模块、方差均值比计算模块、判断因子计算模块、区域判断模块、滤波处理模块和输出模块，

所述的图像数据读取模块，用于读取超声图像的数据；

所述的邻域选取模块，与所述的图像数据读取模块相连接，用于对所述读取的超声图像数据的每一个像素点，选取以其为中心的一定大小的区域为所述像素点对应的一个邻域；

所述的方差均值比计算模块，与所述的邻域选取模块相连接，用于计算所述邻域内像素点各个方向的方差均值比；

所述的判断因子计算模块，与所述的方差均值比计算模块相连接，用于根据各个方向上的方差均值比计算出判别因子；

所述的区域判断模块，与所述的判断因子计算模块相连接，用于根据判别因子分别区分出像素点邻域为边缘区、非边缘区和半边缘区；

所述的滤波处理模块，与所述的区域判断模块相连接，用于对不同的边缘区分别进行滤波处理；

所述的输出模块，与所述的滤波处理模块相连接，用于输出处理后的超声图像数据。

2.根据权利要求1所述的一种超声图像去噪的装置，其特征在于，还包括处理完成判断模块，与所述的邻域选取模块、滤波处理模块和输出模块相连接，用于根据所述滤波处理模块的处理结果判断是否处理完所述读取超声图像的数据中的每一个像素点，如果完成则将结果发送到所述的输出模块，如果未完成，则返回到所述的邻域选取模块继续进行像素点的处理。

3.根据权利要求1或2所述的一种超声图像去噪的装置，其特征在于，所述的滤波处理模块，包括了，

边缘区滤波单元，用于对所述的邻域为边缘区进行增强型方向性滤波；

半边缘区滤波单元，用于邻域为半边缘区进行方向性滤波、方向性增强并加权平均；

非边缘区滤波单元，用于邻域为非边缘区进行各向同性平滑滤波。

4.一种超声图像去噪的方法，其特征在于，包括如下步骤，

步骤1，读取超声图像数据；

步骤2，以各个像素点为中心选取一个邻域；

步骤3，计算邻域内像素点各个方向的方差均值比；

步骤4，根据所述的方差均值比计算出判别因子；

步骤5，根据所述的判别因子分别区分出像素点邻域为边缘区、非边缘区和半边缘区；

步骤6，对所述不同的边缘区分别进行滤波处理；

步骤7，输出处理后的超声图像数据。

5.根据权利要求4所述的一种超声图像去噪的方法，其特征在于，所述的步骤3，进一步包括，

步骤31，计算邻域内像素点四个方向的方差均值比；

步骤32，计算出各个方差均值比的最大值和最小值。

6.根据权利要求4所述的一种超声图像去噪的方法，其特征在于，所述的步骤4，进一步包括，根据所述方差均值比的最大值和最小值之差与最大值的比值，计算出所述的判别因子。

7.根据权利要求4所述的一种超声图像去噪的方法，其特征在于，所述的步骤5，进一步包括，当判别因子趋近于1时，所述邻域为图像的边缘区；当判别因子趋近于0时，所述邻域为图像的非边缘区；当判别因子在0与1之间时，所述区域为图像的半边缘区。

8.根据权利要求4所述的一种超声图像去噪的方法，其特征在于，所述的步骤6，进一步包括，

步骤61，邻域为边缘区进行增强型方向性滤波；

步骤62，邻域为半边缘区进行方向性滤波、方向性增强并加权平均；

步骤63，邻域为非边缘区进行各向同性平滑滤波。

9.根据权利要求4所述的一种超声图像去噪的方法，其特征在于，所述的步骤7之前还包括，判断是否处理完所述读取超声图像的数据中的每一个像素点，如果完成对每一个像素点的处理则进入所述的步骤7，如果没有则继续从所述步骤2开始对未处理的像素点进行相应的计算处理。

10.根据权利要求4所述的一种超声图像去噪的方法，其特征在于，所述的步骤7，进一步包括，将经过处理后的最终图像数据进行输出，包括显示和或存储和或打印。