在利用二维小波变换给图像进行小波去噪时,用下面的算法总是出现错误
%读入原始图像并显示
i=imread('lena.jpg');
subplot(2,2,1);
imshow(i);
title('原始图像');
axis square;
%生成含噪图像并显示
j=imnoise(i,'gaussian',0,0.01);
subplot(2,2,2);
imshow(j);
title('含噪图像');
axis square;
%用sym4小波函数对j进行2层分解
[c,l]=wavedec2(j,2,'sym4');
%实现低通滤波消噪
a1=uint8(wrcoef2('a',c,l,'sym4',2));
%用coif2小波函数对j进行2层分解
[gc,gl]=wavedec2(a1,2,'coif2');
n=[1,2];%设置尺度向量
p=[10.28,24.08];%设置阈值向量
%对三个高频系数进行阈值处理
nc=wthcoef2('h',gc,gl,n,p,'s');
nc=wthcoef2('v',gc,gl,n,p,'s');
nc=wthcoef2('d',gc,gl,n,p,'s');
mc=wthcoef2('h',gc,gl,n,p,'s');
mc&#