NLM非局部均值滤波

NLM的源代码。使用NLM，对高斯噪声，乘性噪声，椒盐噪声进行性能测试

clcclear allclose all;
[filename pathname] = uigetfile（{‘*.jpg‘; ‘*.bmp‘; ‘*.gif‘} ‘选择图片‘）;
if filename==0
    return;
end
I = imread（[pathname filename]）;
sigma=10;

[Ay Ax dim] = size（I）;

%转换为灰度图
if dim>1
    I = rgb2gray（I）;
end

x=（10^2）/（255^2）;
% add  noise


%I1=imnoise（I‘gaussian‘0x）;

%I2=imnoise（I‘salt & pepper‘0.01）;
I3=imnoise（I‘speckle‘0.01）;
%I1=double（I1）;
%I2=double（I2）;
I3=double（I3）;

% denoise
%if dim>2
 %   for i=1:dim     
 %     II3（::i）=NLmeansfilter（I3（::i）52sigma）;
%    end
%end

      
%II1=NLmeansfilter（I1217sigma）;
%II2=NLmeansfilter（I2217sigma）;

II3=NLmeansfilter（I3217sigma）;

 PSNR=0;
 for i=1:Ax
     for j=1:Ay
         PSNR=PSNR+（II3（ij）-I3（ij））^2;
     end
 end
 PSNR=PSNR/Ax/A

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        258  2018-11-22 18:23  make_kernel.m

     文件       1612  2018-11-28 22:27  NLmeansfilter.m

     文件       1225  2018-11-29 17:12  go.m

----------- ---------  ---------- -----  ----

                 3095                    3