一种图像去噪方法技术

本发明专利技术提供一种图像去噪方法，包括步骤：S1)提取待处理图像的边缘信息E和纹理信息T，并存储；S2)对待处理图像进行二维经验模式分解，得到图像的内禀模式函数IMF和最终的大趋势项R；S3)对高频IMF进行基于边缘和纹理特征的去噪；S4)将去噪后的高频IMF与低频IMF以及大趋势项R合成处理，得到最终的增强图像。本发明专利技术通过对原有二维经验模式分解方法的改进，在二维经验模式分解得到高频分量后，参考原始图像的纹理特征和边缘信息，进行改进的非局部均值去噪，能较好地去除高斯噪声并保留图像的细节信息，解决了图像变得模糊的问题，较原来的算法取得更好的效果。

【技术实现步骤摘要】
一种图像去噪方法
本专利技术涉及图像处理
，特别是涉及一种图像去噪方法。
技术介绍
在数字图像处理中，图像去噪一直是一个非常关键的环节，去噪效果的好坏直接影响到后续的处理以及视觉效果。目前的图像去噪算法多种多样，基本可分为基于空域的滤波、基于小波变换的滤波、基于模糊理论滤波、基于各向异性扩散(ADE)滤波及基于图像分解的滤波等等。基于二维经验模式分解的图像增强已经取得一些较好的效果。经验模式分解(EMD)方法是1998年Huang提出的一种新的信号处理方法，在非平稳信号分析方面有较好性能。该方法将复杂的非平稳信号分解为若干具有不同尺度的平稳数据分量和趋势项的叠加，具有自适应性。EMD在一维信号处理上取得了很好效果，国内学者将其推广到二维，提出二维经验模式分解(BEMD)方法。然而，二维经验模式分解BEMD应用在去噪方面，如果单纯将分解后的高频成分去除，则会损失大量细节信息，图像变得模糊，图像变得模糊究其原因是没有很好的保留图像的边缘和纹理等细节信息。因而，需要将二维经验模式分解(BEMD)方法与其他去噪方法相结合才会取得较理想的效果。
技术实现思路
鉴于以上所述现有技术的缺点，本专利技术的目的在于提供一种图像去噪方法，用于解决现有技术中二维经验模式分解方法单纯将分解后的高频成分去除，导致图像大量细节信息的损失、图像变得模糊的问题。为实现上述目的，本专利技术提供一种图像去噪方法，将二维经验模式分解和纹理边缘等形态学去噪方法结合起来，在去除噪声的同时，尽量保留图像细节信息，以得到较好的去噪效果。本专利技术是采用以下方案来实现的：一种图像去噪方法，包括以下步骤：1)提取原始图像f(x,y)的边缘信息E和纹理信息T，并存储；2)对原始图像f(x,y)进行二维经验模式分解，得到图像的内禀模式函数IMF和最终的大趋势项R；3)对高频内禀模式函数IMF进行基于边缘和纹理特征的非局部均值去噪；4)将去噪后的高频内禀模式函数IMF与低频内禀模式函数IMF以及大趋势项R合成处理，得到最终的增强图像。于本专利技术的一实施方式中，在所述步骤1)中，通过边缘检测算子和滤波器组提取原始图像f(x,y)的边缘信息E和纹理信息T。于本专利技术的一实施方式中，在所述步骤2)中，对原始图像f(x,y)进行二维经验模式分解的步骤包括：21)提取原始图像f(x,y)的局部极大值和极小值；22)利用插值方法建立极大值层和极小值层，确保原始图像f(x,y)的所有点在这两层之间；23)计算极大值层和极小值层的均值；24)计算原始图像f(x,y)和所述均值的差；25)重复步骤21)～24)，直到标准差SD比预先定义的阈值小时，得到第一个内禀模式函数IMF1；26)将第一个内禀模式函数IMF1从f(x,y)中提取出来，将余项作为新数据，重复步骤21)～25)n次，直到没有更多的内禀模式函数能够被提取出来则停止；27)最终的原始图像表示成：其中，f(x,y)为原始图像数据，cj(x,y)为第j个内禀模式函数，是分解得到的较小尺度细节信息，rn(x,y)是第n个余项，亦即最终得到的大尺度趋势项。于本专利技术的一实施方式中，在所述步骤25)中，预先定义的标准差SD的阈值范围为0.2～0.3。于本专利技术的一实施方式中，在步骤3)中进行非局部均值去噪时，包括以下步骤：31)计算原始图像的纹理信息T和边缘信息E，并归一化为Tnorm和Enorm，存储待查；32)通过Tnorm和Enorm计算权值wte＝f(Tnorm,Enorm)，用来选择搜索窗口和比较窗口大小，若像素i有更多纹理信息，则搜索窗口和比较窗口都减小；若像素i在边缘区域，则扩大搜索窗口，减小比较窗口；33)根据确定的比较窗口和搜索窗口大小，计算像素i和j的相似程度：其中，表示像素i和像素j之间以它们为中心的子块之间高斯距离平方，u(Ni)表示i周围局部子块像素集合，h为滤波参数，h越大，滤波效果越强，0≤w(i,j)≤1，且满足：∑w(i,j)＝1；34)设含噪声的图像为u(i)，去除噪声后的图像为NL(u)(i)，对每个像素的值通过加权得到：NL(u)(i)＝∑w(i,j)u(j)。如上所述，本专利技术的图像去噪方法，具有以下有益效果：在基于高频分量去噪时，依据提取的图像的边缘和纹理信息，自适应的调节去噪系数，较好地保留了图像的细节；且根据各个IMF分量所包含的噪声强度越来越弱的特点，只针对前几个高频IMF分量去噪而保留低频IMF效果会更好。附图说明图1显示为本专利技术图像去噪方法的流程图。图2显示为本专利技术图像去噪方法中二维经验模式分解的流程图。元件标号说明S1～S4步骤S21～S27步骤具体实施方式以下通过特定的具体实例说明本专利技术的实施方式，本领域技术人员可由本说明书所揭露的内容轻易地了解本专利技术的其他优点与功效。本专利技术还可以通过另外不同的具体实施方式加以实施或应用，本说明书中的各项细节也可以基于不同观点与应用，在没有背离本专利技术的精神下进行各种修饰或改变。需说明的是，在不冲突的情况下，以下实施例及实施例中的特征可以相互组合。需要说明的是，以下实施例中所提供的图示仅以示意方式说明本专利技术的基本构想，虽图示中仅显示与本专利技术中有关的组件而非按照实际实施时的组件数目、形状及尺寸绘制，其实际实施时各组件的型态、数量及比例可为一种随意的改变，且其组件布局型态也可能更为复杂。请参阅图1，为本专利技术图像去噪方法的流程图，所述图像去噪方法包括以下步骤：S1)提取待处理图像的边缘和纹理信息，并存储待查；S2)对待处理图像进行二维经验模式分解，得到图像的内禀模式函数IMF和最终的大趋势项R；S3)对高频IMF进行基于边缘和纹理特征的去噪；S4)将去噪后的高频IMF与低频IMF以及大趋势项R合成处理，得到最终的增强图像。作为示例，在所述步骤S1)中，通过边缘检测算子和滤波器组提取原始图像f(x,y)的边缘和纹理信息。作为示例，在所述步骤S2)中，对待处理图像进行二维经验模式分解，得到图像的内禀模式函数IMF和最终的大趋势项R。需要注意的是，所述内禀模式函数(也称固有模态函数)需要满足两个条件：一、在整个数据集合上，极值点和零交叉点个数相等或者相差1；二、在任意时间点，由极大值定义的包络和极小值定义的包络均值为0。作为示例，请参阅图2，在所述步骤S2)中，对原始图像f(x,y)进行二维经验模式分解的步骤包括：S21)、提取待处理图像f(x,y)的局部极大值点和极小值点，并进行插值处理，得到极大值包络为Emax(x,y)，极小值包络为Emin(x,y)。S22)、计算平均包络，即极大值包络Emax(x,y)和极小值包络Emin(x,y)的代数均值，记为e1(x,y)：S23)、计算待处理图像f(x,y)与平均包络e1(x,y)的差值，记为d1(x,y)：d1(x,y)＝f(x,y)-e1(x,y)。S24)、判断是否满足内禀模式函数条件，若不满足，则重复步骤S21)～S23)过程k次，直到d1k(x,y)为內禀模式函数：d1(x,y)-e11(x,y)＝d11(x,y)......d1(k-1)(x,y)-e1k(x,y)＝d1k(x,y)；定义c1(x,y)＝d1k(x,y)，满足内禀模式函数IMF条件，即为从原始数据中分离出的第一个内禀模式函数I本文档来自技高网...

【技术保护点】
一种图像去噪方法，其特征在于，包括以下步骤：1)提取原始图像f(x,y)的边缘信息E和纹理信息T，并存储；2)对原始图像f(x,y)进行二维经验模式分解，得到图像的内禀模式函数IMF和最终的大趋势项R；3)对高频内禀模式函数IMF进行基于边缘和纹理特征的非局部均值去噪；4)将去噪后的高频内禀模式函数IMF与低频内禀模式函数IMF以及大趋势项R合成处理，得到最终的增强图像。

【技术特征摘要】
1.一种图像去噪方法，其特征在于，包括以下步骤：1)提取原始图像f(x,y)的边缘信息E和纹理信息T，并存储；2)对原始图像f(x,y)进行二维经验模式分解，得到图像的内禀模式函数IMF和最终的大趋势项R；3)对高频内禀模式函数IMF进行基于边缘和纹理特征的非局部均值去噪；4)将去噪后的高频内禀模式函数IMF与低频内禀模式函数IMF以及大趋势项R合成处理，得到最终的增强图像。2.根据权利要求1所述的图像去噪方法，其特征在于，在所述步骤1)中，通过边缘检测算子和滤波器组提取原始图像f(x,y)的边缘信息E和纹理信息T。3.根据权利要求1所述的图像去噪方法，其特征在于，在所述步骤2)中，对原始图像f(x,y)进行二维经验模式分解的步骤包括：21)提取原始图像f(x,y)的局部极大值和极小值；22)利用插值方法建立极大值层和极小值层，确保原始图像f(x,y)的所有点在这两层之间；23)计算极大值层和极小值层的均值；24)计算原始图像f(x,y)和所述均值的差；25)重复步骤21)～24)，直到标准差SD比预先定义的阈值小时，得到第一个内禀模式函数IMF1；26)将第一个内禀模式函数IMF1从f(x,y)中提取出来，将余项作为新数据，重复步骤21)～25)n次，直到没有更多的内禀模式函数能够被提取出来则停止；27)最终的原始图像表示成：其...

【专利技术属性】
技术研发人员：刘立庄，张浩，赵丹，
申请(专利权)人：中国科学院上海高等研究院，
类型：发明
国别省市：上海,31