clear
clc
clf
colormap（gray）

%?create?example?image
ima=100*ones（100）;
ima（50:100:）=50;
ima（:50:100）=2*ima（:50:100）;
fs=fspecial（‘a(chǎn)verage‘）;
ima=imfilter（imafs‘symmetric‘）;

%?add?some?noise
sigma=10;
rima=ima+sigma*randn（size（ima））;

%?show?it
imagesc（rima）
drawnow

%?denoise?it
fima=NLmeansfilter（ima52sigma）;

%?show?results
clf
subplot（221）imagesc（ima）title（‘original‘）;
subplot（222）imagesc（rima）title（‘noisy‘）;
subplot（223）imagesc（fima）title（‘filtered‘）;
subplot（224）imagesc（rima-fima）title（‘residuals‘）;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件????????2507??2008-08-01?02:28??NLmeansfilter.m
?????文件?????????567??2008-08-01?02:51??demo.m
?????文件????????1348??2009-05-25?11:58??license.txt