"Gradient Domain Guided Image Filtering"论文中边缘权重函数matlab实现代码

"Kou F, Chen W, Wen C, et al. Gradient Domain Guided Image Filtering[J]. Image Processing, IEEE Transactions on, 2015, 24(11): 4528-4539."论文中权重函数matlab实现的代码。

clear all;

g = double(imread('tulips.png'))/255.0;
[m, n, d]=size(g);
if(size(g,3)==3) 
    g=rgb2gray(g); 
end

%% weight 1
% epsilon=(0.001*(max(g(:))-min(g(:))))^2;
% r=1;
% 
% N = boxfilter(ones(m, n), r); % the size of each local patch; N=(2r+1)^2 except for boundary pixels.
% mean_I = boxfilter(g, r) ./ N;
% mean_II = boxfilter(g.*g, r) ./ N;
% var_I = mean_II - mean_I .* mean_I;
% 
% weight1=(var_I+epsilon)/(mean(var_I(:))+epsilon);       

%% weight 2
epsilon=(0.001*(max(g(:))-min(g(:))))^2;
r1=1;
r2=16;

N1 = boxfi