nlmeans算法

nlmeans算法 matlab code

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（‘average‘）;
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‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        2507  2008-08-01 02:28  NLmeansfilter.m
     文件         567  2008-08-01 02:51  demo.m
     文件        1348  2009-05-25 11:58  license.txt