快速非局部均值滤波算法

clear all;
close all;
ima =  double（imread（‘barbara.png‘））;
hD = 3;
hW = 10;
D=2*hD;
W=2*hW+1;
D1=D+1;

sigma=20;
% h =6*sigma;
% h=h^2;
h = quantile_gauss（D1 0.88） * （4* sigma^2） .* D1.^2;
randn（‘seed‘0）;
ima_nse=ima+sigma*randn（size（ima））;
imshow（ima_nse[]）;
[widthheight] = size（ima_nse）;
ima_nse2 = padarray（ima_nse[hD hD]‘symmetric‘）; 
ima_nse3 = padarray（ima_nse[hW+hD hW+hD]‘symmetric‘）; 
output =zeros（widthheight）;
Z =zeros（widthheight）;

k=1;
w = cell（W*W）;
for m = -hW:hW
    for n = -hW:hW
        if（m==0 && n==0） continue; end;    %保证搜索窗比较窗不重叠

        Sd = （ima_nse2-ima_nse3（hW+1+m:hW+m+width+DhW+1+n:hW+n+height+D））.^2;
        Sd = cumsum