一种基于参数选择的图像自适应指导滤波方法
【技术领域】
[0001] 本发明涉及一种图像自适应指导滤波方法，尤其是一种基于参数选择的图像自适 应指导滤波方法，属于图像处理技术领域。
【背景技术】
[0002] 在计算机视觉和计算机图像学中，滤波技术被应用到大部分领域，其主要作用是 抑制或提取图像的某些内容。经典的平均滤波和高斯滤波属于线性滤波，采用的是加权平 均的思想，其核函数具有空间不变性且独立于图像的内容，在抑制噪声的同时，也模糊了图 像的边缘。双边滤波是一种典型的非线性滤波，考虑了像素的空间邻近关系，但传统双边滤 波存在的非线性估计不足问题。因此，指导滤波应运而生。指导滤波是一种局部的线性模 型，在滤波过程中加入一个额外的指导图像并显式地表现在滤波器核函数之中，能快速实 现图像平滑滤波和保留边缘的功能，在图像去噪、图像增强、图像去雾等领域有着广泛的应 用前景。
[0003] 图像的指导滤波是一个线性移可变的滤波过程，包括指导图像G，输入图像I和输 出图像〇。其中指导图像G可以根据应用事先设定，也可以直接取为输入图像I。
[0004] 对于输出图像中的第i个像素 O1，其计算方法为：
[0005]
:⑴：
[0006] 其中i，j为像素位置，Wu为指导滤波器的核函数。其定义为：
[0007]
[0008] 式中，Wk为第k个核函数窗口，一般来讲，各核函数窗口内的像素个数相等，均为 |w|，μ^Ρ 4是指导图像G在第k个核函数窗口内灰度值的均值和方差，ε是平滑因子。 k: (i，j)表示第k个核函数窗口的像素位置，&和G j是该核函数窗口内第i个像素和第j 个像素的灰度值。
[0009] 在上述的指导滤波的原理中，主要涉及两个参数的选择：滤波窗口 w和平滑因子 ε。传统指导滤波未考虑窗口内像素灰度值偏差对滤波的意义，也就是缺少对窗口内最大 值和最小值相对均值的偏差分析。为了改进不足，有研究者引入偏移函数I 1，带偏移函数 ξ i的指导滤波器定义如下：
[0010]
[0011] 其中，偏移函数ξ』勺选择是：
[0012]
[0013] 式中，偏移量A1= G1-Pk。此方法的偏移函数ξ ^艮据各核函数窗口内像素灰度 值的最大值MAX (wk)和最小值MIN(Wk)确定。当核函数窗口内第i个像素的灰度值Gi大于 窗口内灰度均值μ 3寸，偏移函数ξ 1取窗口内最大值MAX (Wk)与该点灰度值的差；同理，当 核函数窗口内第i个像素的灰度值G1小于窗口内灰度均值μ ,时，偏移函数ξ 1取该核函 数窗口内像素灰度值的最小值MIN(Wk)与该点灰度值的差；当该核函数窗口内第i个像素 的灰度值G 1等于该核函数窗口内像素灰度值的均值μ ,时，偏移函数ξ i取零，该像素的灰 度值不做任何处理。
[0014] 在指导滤波的核函数之中加入偏移函数，在边缘保持方面效果有所提高。但偏移 函数1的选择相对片面，在选择偏移函数ξ i之前，缺少对每一个核函数窗口内各像素灰 度值的具体分析，未考虑到窗口内各像素灰度值之间的联系，造成滤波后图像边缘存在伪 影。而且这种方法并没有考虑到指导滤波自身参数平滑因子的选择对滤波结果的影响，整 幅图像的平滑因子ε只选用单一的固定值。
[0015] 当指导图像G的当前窗口方差较大时，说明图像在此窗口内变化较多，包含边缘 信息或突出信息。此时窗口的方差将大于平滑因子ε，那么应选择一个较小的ε，这样经 过核函数滤波之后，滤波图像将更接近指导图像的变化特征。当指导图像G的当前窗口的 方差较小时，说明图像在此窗口内变化较小，相对平滑，此时窗口的方差将小于平滑因子 ε，那么应选择一个较大的ε，使滤波后的图像保持原有指导图像的平滑特征。
[0016] 因此，平滑因子ε的选择不应是单一固定值，应考虑指导图像滤波窗口的方差。 在上述方法中，都没有考虑选择不同的平滑因子来进一步对滤波窗口内的边缘和细节信息 进行处理。
【发明内容】
[0017] 本发明要解决的技术问题是提供一种基于参数选择的图像自适应指导滤波方法。
[0018] 本发明采用下述技术方案：
[0019] -种基于参数选择的图像自适应指导滤波方法，包括以下步骤：
[0020] 步骤1 :选择输入图像I和指导图像G，设定滤波窗口的大小w = (2*r+l)*(2*r+l)，r为窗口半径，|w|代表窗口内像素点的个数；
[0021] 步骤2 :分别以各像素为滤波窗口的中心像素，使用所述核函数滤波输入图像，包 括以下具体步骤：
[0022] 步骤2-1 :判断指导图像G的当前滤波窗口中各像素点的灰度值Gpl < i < |w 和中心像素点的灰度值Gni是否相似，如果不相似，将所述像素点的灰度值G i置为0,否则所 述像素点的灰度值G1保持原值；
[0023] 步骤2-2 :统计当前滤波窗口中和中心像素点的灰度值Gni不相似的像素点的个 数；判断其是否超过预设阈值，如果是，设置偏移函数ξ 〇 ;否则，设置偏移函数为：
[0024] \；
[0025] 其中偏移量Λ FGfyk, MIN(Wk)、MAX(Wk)、别为指导图像G的当前滤波窗 口内的最小像素灰度值、最大像素灰度值和像素灰度均值；G i, I < i < |w|为指导图像G当 前滤波窗口内的第i个像素的灰度值；
[0026] 步骤2-3 :确定平滑因子ε，逐一计算各滤波窗口中像素灰度值的方差并保存在 方差矩阵Var ;
[0027] 步骤2-4 :依次判断各滤波窗口的方差是否大于等于方差矩阵Var的均值，如果 是，定义为大方差滤波窗口，否则，定义为小方差滤波窗口；大方差滤波窗口内的平滑因子 小于小方差滤波窗口内的平滑因子，按如下公式选择：
[0028]
[0029] 其中i和j代表像素坐标，μ _代表方差矩阵Var的均值；Var (i，j)代表方差矩 阵Var在（i，j)点的值；
[0030] 步骤2-5 :将平滑因子设置为核函数的参数，使用核函数计算滤波后像素的灰度 值，输出滤波后的图像。
[0031] 所述步骤2-1中判断所述指导图像G在滤波窗口中各像素点Gpl彡i彡|w|和中 心像素点Gni的灰度值是否相似的方法为判断所述像素点G 灰度值与中心像素点G "的灰 度值的差值的绝对值是否小于所述指导图像G在滤波窗口中像素灰度值标准差σ ,的1/3。
[0032] 所述步骤2-3中，大方差滤波窗口中的平滑因子和小方差滤波窗口中的平滑因子 ε根据相似性判断后的方差矩阵Var自适应选择。
[0033] 采用上述技术方案所产生的有益效果在于：
[0034] 1、本发明综合考虑偏移函数I1和平滑因子ε对滤波结果的影响，在二者选择前 加入灰度相似性判断，实现了偏移函数ξ i的全面选择，并依据滤波窗口内的方差大小自动 确定平滑因子ε，共同实现二者的自适应选择。
[0035] 2、本发明有区别地处理窗口内的局部边缘特征，因此滤波后的图像在边缘和细节 方面具有良