splitting Bregman

splitting Bergman 对应matlab程序

%% 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‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件      153785  2011-02-03 17:24  SplitBregmanTVdenoising\Lena512.png
     文件        1189  2012-04-19 10:34  SplitBregmanTVdenoising\SB_ATV.m
     文件         996  2012-04-19 10:45  SplitBregmanTVdenoising\main_denoising.m
     文件        1196  2012-04-19 10:45  SplitBregmanTVdenoising\SB_ITV.m
     文件        1510  2014-02-12 13:58  license.txt