Image denosingppbNakagami

Image denosingppb Nakagami Matalb实现

clear all

disp（‘Data loading and generation...‘）;
ima = double（imread（‘lena.png‘））;
L = 3;
s = zeros（size（ima））;
for k = 1:L
    s = s + abs（randn（size（ima）） + i * randn（size（ima）））.^2 / 2;
end
ima_nse = ima .* sqrt（s / L）;

disp（‘Non-it PPB filter...‘）;
hW = 10;
hD = 3;
alpha = 0.88;
T = inf;
nbit = 1;
ima_fil_1 = ppb_nakagami（ima_nse L ...
                         hW hD ...
                         alpha T ...
                         nbit）;

disp（‘It PPB filter...‘）;
hW = [1 3 5 10];
hD = [0 1 2 3];
alpha = 0.92;
T = 0.2;
nbit = [1 2 3 4];
ima_fil_2 = ppb_nakagami（ima_nse L ...
                         hW hD ...
                         alpha T ...
                         nbit）;

figure
subplot（221） plotimage_sar（ima）;
title（‘Noise-free image‘）;
subplot（222） plotimage_sar（ima_nse ima）;
title（sprintf（‘Noisy image （PSNR %.3f）‘ snr（ima_nse ima）））;
subplot（223） plotimage_sar（ima_fil_1 ima）;
title（sprintf（‘Non-it PPB filter （PSNR %.3f）‘ snr（ima_fil_1 ima）））;
subplot（224） plotimage_sar（ima_fil_2 ima）;
title（sprintf（‘It PPB filter （PSNR %.3f）‘ snr（ima_fil_2 ima）））;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1114  2010-03-23 13:03  ppbNakagami\example_nakagami.m
     文件        1114  2010-03-23 13:03  ppbNakagami\exemple_nakagami.m
     文件      151199  2010-03-23 13:03  ppbNakagami\lena.png
     文件         187  2010-03-23 13:03  ppbNakagami\plotimage.m
     文件         326  2010-03-23 13:03  ppbNakagami\plotimage_sar.m
     文件        7029  2013-05-17 08:21  ppbNakagami\ppb_nakagami.m
     文件      637537  2010-03-23 13:03  ppbNakagami\ppbNakagami.mexa64
     文件      618644  2010-03-23 13:03  ppbNakagami\ppbNakagami.mexglx
     文件       77824  2010-03-23 13:03  ppbNakagami\ppbNakagami.mexw32
     文件      263168  2012-05-26 13:12  ppbNakagami\ppbNakagami.mexw64
     文件         137  2010-03-23 13:03  ppbNakagami\psnr.m
     文件        2624  2012-05-26 13:29  ppbNakagami\README
     文件         129  2010-03-23 13:03  ppbNakagami\snr.m