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