【发布时间】:2013-11-03 14:59:27
【问题描述】:
我已经按照教授的解释编写了图像压缩代码。它工作正常,但对于某些图像,它显示以下错误:
**Error using +
Matrix dimensions must agree.
Error in idwt2 (line 90)
x = upsconv2(a,{Lo_R,Lo_R},sx,dwtEXTM,shift)+ ... % Approximation.
Error in nw (line 56)
DeCompImg=idwt2(cA1,c,cV1,cD1,'haar');**
该代码适用于少数图像,但对于少数图像显示上述错误。我检查了工作区。仅当原始图像尺寸为奇数(255x255)时才会出现此错误。如果图像尺寸是偶数(256x256),则此代码可以正常工作。如何解决这个问题without using imresize 命令?请帮忙。
我的压缩代码:
%Reading an Image
I=imread('mandrill.jpg');
G=rgb2gray(I);
%Applying Haar Wavelet
[cA1,cH1,cV1,cD1] = dwt2(G,'haar');
%Applying 2x2 window averaing
CompImgH=blockproc(cH1, [2 2], @(x) mean(x.data(:)));
CompImgV=blockproc(cV1, [2 2], @(x) mean(x.data(:)));
CompImgD=blockproc(cD1, [2 2], @(x) mean(x.data(:)));
figure(1),
subplot(2,2,1);
imshow(G);
title('Original Image');
subplot(2,2,2);
imshow(CompImgH);
title('Horizontal Component');
subplot(2,2,3);
imshow(CompImgV);
title('Vertical Component');
subplot(2,2,4);
imshow(CompImgD);
title('Diagonal Component');
%DECOMPRESSION
%Inverse process for 2x2 window averaging
b=CompImgH;
[m,n,colormap]=size(b);
k=1; %Counter for Row and
l=1; %Column replication
%f=input('enter replication factor: ');
f=2; % replication factor
for i=1:m %Loop for reading row and
for t=1:f %Row replication
for j=1:n %Loop for reading column and
for t=1:f %Column replication
c(k,l)=b(i,j);
l=l+1;
end
end
l=1;
k=k+1;
end
end
DeCompImg=idwt2(cA1,c,cV1,cD1,'haar');
DecompressedImage=uint8(DeCompImg);
Orig_Image = im2double(G);%---Convert image to double class
Reconstructed_Image = im2double(DecompressedImage);%---Convert image to double class
[M N] = size(Orig_Image);%---Size of Original Image
err = Orig_Image - Reconstructed_Image;%---Difference between two images
MSE = (sum(sum(err .* err)))/(M * N);
【问题讨论】:
-
dwtmode('status')返回什么?另外,在循环之后,idwt2之前的size(c)是什么? -
您的错误发生在
idwt2。您可能正在舍入/取整某些东西(也许除以 2),因此当您有奇数时会遇到问题。或者你不能在你划分的地方四舍五入/地板,你应该。查找任何除法,并确保在出现奇数时做正确的事情。 -
我尝试使用以下命令来消除此错误。它在 M=N 时有效,但在 M 不等于 N 时无效。%[p q]= size(G); % if mod(p,2) && mod(q,2) % G=padarray(G,[1 1],'pre'); % 其他
标签: matlab image-processing image-compression