一、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的像素为中心计算高斯加权欧式距离