www.gusucode.com > 《matlab图像处理与界面编程宝典》秦襄培 编著，每章的MATLAB源代码程序 > 第20章/代码20-3.txt

    
clear;                            % 清理工作空间     
load wbarb;                       % 装载原始图像               
figure;                           % 新建窗口           
image(X);                         % 显示图像
colormap(map);                    % 设置色彩索引图
title('原图');                    % 设置图像标题
axis square;                      % 设置显示比例
%生成含噪图像
init=2055615866;                  % 初始值
randn('seed',init)                % 随机值    
x=X+38*randn(size(X));            % 添加随机噪声
figure;                           % 新建窗口
image(x);                         % 显示图像
colormap(map);                    % 设置色彩索引图
title('含噪声图像');              % 设置图像标题
axis square;                      % 设置显示比例    
% 下面对图像消噪，首先利用sym4函数对图像进行2层分解
[c,s]=wavedec2(x,2,'sym4');       % 分解  
a1=wrcoef2('a',c,s,'sym4',1);     % 图像第一层的重构逼近信号　
figure;                           % 新建窗口
image(a1);                        % 显示图像
colormap(map);                    % 设置色彩索引图
title('第一层重构图像');          % 设置图像标题
axis square;                      % 设置显示比例
a2=wrcoef2('a',c,s,'sym4',2);     % 图像第二层的重构逼近信号
figure;                           % 新建窗口
image(a2);                        % 显示图像
colormap(map);                    % 设置色彩索引图
title('第二层重构图像');          % 设置图像标题
axis square;                      % 设置显示比例