program_16_23.m

来自「这是基于MATLAB系统上开发的,应用于小波分析方面,很有用的程序模块,书名为"」· M 代码 · 共 29 行

%装载并显示原始图像
load bust;
%加入含噪
init=2055615866;
randn('seed',init);
X1=X+20*randn(size(X));
subplot(221);image(X1);
%colormap(map);
title('原始图像');axis square;
%用小波db4对图像X进行一层小波包分解
T=wpdec2(X1,1,'db4');
%重构图像近似部分
A=wprcoef(T,[1 0]);
subplot(222);image(A);
title('图像的近似部分');
axis square;
%边缘检测
%%原图像的边缘检测
BW1 = edge(A,'sobel');
subplot(223);imshow(BW1);
title('原图像的边缘');
axis square;
%图像近似部分的边缘检测
BW2= edge(X1,'sobel');
subplot(224);imshow(BW2);
title('图像近似部分的边缘');
axis square;