双边滤波的基本特性:
普通的时空域的低通滤波器,在像素空间完成滤波以后,导致图像的边缘部分也变得不那么明显,
整张图像都变得同样的模糊,图像边缘细节丢失。双边滤波器(Bilateral Filter)可以很好的保
留边缘的同时消除噪声。双边滤波器能做到这些原因在于它不像普通的高斯/卷积低通滤波,只考
虑了位置对中心像素的影响,它还考虑了卷积核中像素与中心像素之间相似程度的影响,根据位置
影响与像素值之间的相似程度生成两个不同的权重表(WeightTable),在计算中心像素的时候加以
考虑这两个权重,从而实现双边低通滤波。据说AdobePhotoshop的高斯磨皮功能就是应用了
双边低通滤波算法实现。
数学公式:
其中为输出影像,
为输入影像,函数的值代表给定像素上的色彩资讯。 我们将输入影像
和输出影像
写成向量的原因是,操作的影像不限定是只有单一通道的灰阶影像,他们也可以是多通道的彩色影像。 式子中的
代表影像中的任意一点,而
代表在
邻近的点。 而
为归一化函数:
双边滤波器的目的是,对于影像中的每一个点,都要对邻近的所有点进行加权平均,以达到平滑化的目的。 而对此,双边滤波器在邻近的点
采取了不同的权重, 策略包含了两个核心函数
和
, 这两个式子分别代表了
和
在空间几何关系上的差异和光度/色彩差异上的差异。 以下将针对这两个核心函数进行说明。
matlab代码:
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Implements bilateral filter for color images.
function B = bfltColor(A,w,sigma_d,sigma_r)
% Convert input sRGB image to CIELab color space.
if exist('applycform','file')
A = applycform(A,makecform('srgb2lab'));
else
A = colorspace('Lab<-RGB',A);
end
% Pre-compute Gaussian domain weights.
[X,Y] = meshgrid(-w:w,-w:w);
G = exp(-(X.^2+Y.^2)/(2*sigma_d^2));
% Rescale range variance (using maximum luminance).
sigma_r = 100*sigma_r;
% Create waitbar.
h = waitbar(0,'Applying bilateral filter...');
set(h,'Name','Bilateral Filter Progress');
% Apply bilateral filter.
dim = size(A);
B = zeros(dim);
for i = 1:dim(1)
for j = 1:dim(2)
% Extract local region.
iMin = max(i-w,1);
iMax = min(i+w,dim(1));
jMin = max(j-w,1);
jMax = min(j+w,dim(2));
I = A(iMin:iMax,jMin:jMax,:);
% Compute Gaussian range weights.
dL = I(:,:,1)-A(i,j,1);
da = I(:,:,2)-A(i,j,2);
db = I(:,:,3)-A(i,j,3);
H = exp(-(dL.^2+da.^2+db.^2)/(2*sigma_r^2));
% Calculate bilateral filter response.
F = H.*G((iMin:iMax)-i+w+1,(jMin:jMax)-j+w+1);
norm_F = sum(F(:));
B(i,j,1) = sum(sum(F.*I(:,:,1)))/norm_F;
B(i,j,2) = sum(sum(F.*I(:,:,2)))/norm_F;
B(i,j,3) = sum(sum(F.*I(:,:,3)))/norm_F;
end
waitbar(i/dim(1));
end
% Convert filtered image back to sRGB color space.
if exist('applycform','file')
B = applycform(B,makecform('lab2srgb'));
else
B = colorspace('RGB<-Lab',B);
end
% Close waitbar.
close(h);
使用代码如下:I=imread('lena.png');
I=double(I)/255;
I1= bfltColor(I,1,5,0.1);
subplot(1,2,1);
imshow(I);
subplot(1,2,2);
imshow(I1)
处理结果如下: