non-local 论文及代码复现总结

一、NLmeansfilter 函数

function [output]=NLmeansfilter(input,t,f,h)
 %  input: 处理的图像
 %  t: 搜索窗的半径
 %  f: 相似窗的半径
 %  h: 滤波参数
 [m,n]=size(input);%图像的范围
 output=zeros(m,n);
 input2=padarray(input,[f f],'symmetric');%便于后面控制边界的处理
 kernel=make_kernel(f);
 kernel=kernel/sum(sum(kernel));%高斯滤波器
 for i=1:m
     for j=1:n
         i1=i+f;
         j1=j+f;
         W1=input2(i1-f:i1+f,j1-f:j1+f);%使用2f+1区域
         wmax=0;
         average=0;
         sweight=0;
         rmin=max(i1-t,f+1);
         rmax=min(i1+t,m+f);
         smin=max(j1-t,f+1);%不断变化搜索窗的大小
         smax=min(j1+t,n+f);%之前始终没懂这里是干什么的，用于控制搜索窗口的移动
         %遍历搜索21*21的区域
         for r=rmin:1:rmax
             for s=smin:1:smax
                 if(r==i1&&s==j1)%遍历搜索，若遇到的相同的像素跳过
                     %以r,s的像素为中心计算高斯加权欧式距离