分裂Bregman算法的源代码1

偏微分图像去噪中的经典算法，2008年分裂Bregman算法的源代码！

%% This file demonstrates the Split Bregman method for Total Variation denoising
%
%   SB_ATV.m  Split Bregman Anisotropic Total Variation Denoising
%   SB_ITV.m  Split Bregman Isotropic Total Variation Denoising
%
% Benjamin Trémoulhéac
% University College London
% b.tremoulheac@cs.ucl.ac.uk
% April 2012

clc; clear all;
close all;

N = 512; n = N^2;
f = double（imread（‘Lena512‘‘png‘））;
g = f（:） + 0.09*max（f（:））*randn（n1）;

mu = 20;

g_denoise_atv = SB_ATV（gmu）;
g_denoise_itv = SB_ITV（gmu）;

fprintf（‘ATV Rel.Err = %g\n‘norm（g_denoise_atv（:） - f（:）） / norm（f（:）））;
fprintf（‘ITV Rel.Err = %g\n‘norm（g_denoise_itv（:） - f（:）） / norm（f（:）））;

figure; colormap gray;
subplot（221）; imagesc（f）; axis image; title（‘Original‘）;
subplot（222）; imagesc（reshape（gNN））; axis image; title（‘Noisy‘）;
subplot（223）; imagesc（reshape（g_denoise_atvNN））; axis image; 
title（‘Anisotropic TV denoising‘）;
subplot（224）; imagesc（reshape（g_denoise_itvNN））; axis image; 
title（‘Isotropic TV denoising‘）;

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

     文件       1341  2012-04-19 09:00  SplitBregmanTVdenoising\license.txt

     文件     153785  2011-02-03 12:24  SplitBregmanTVdenoising\SplitBregmanTVdenoising\Lena512.png

     文件        996  2012-04-19 06:45  SplitBregmanTVdenoising\SplitBregmanTVdenoising\main_denoising.m

     文件       1189  2012-04-19 06:34  SplitBregmanTVdenoising\SplitBregmanTVdenoising\SB_ATV.m

     文件       1196  2012-04-19 06:45  SplitBregmanTVdenoising\SplitBregmanTVdenoising\SB_ITV.m

     目录          0  2014-04-30 16:55  SplitBregmanTVdenoising\SplitBregmanTVdenoising

     目录          0  2014-04-30 16:55  SplitBregmanTVdenoising

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

               158507                    7