我的工作是使用奇异值分解技术进行图像压缩。我用Matlab为它写了代码。我压缩了一个文件大小为8.15KB的图像(255*255*3)。当我保存压缩图像( jpg格式)时,它的文件大小超过了原始图像文件的大小。
B=imread('lena.jpg');
figure,
imshow(B), title( sprintf('size=%d',numel(B)) ) // displaying the original image
A=im2double(B);
A1=A(:,:,1);
A2=A(:,:,2);
A3=A(:,:,3);
tic;
// applying svd for each layer
[U1,S1,V1]=svd(A1);
[U2,S2,V2]=svd(A2);
[U3,S3,V3]=svd(A3);
// reconstuctin compressed image
p=100;
U1p=U1(:,1:p);
V1p=V1(:,1:p);
S1p=diag(S1(1:p,1:p));
C1=U1p * diag(S1p) * V1p';
C1=255*C1;
C1=uint8(C1);
U2p=U2(:,1:p);
V2p=V2(:,1:p);
S2p=diag(S2(1:p,1:p));
C2=U2p * diag(S2p) * V2p';
C2=255*C2;
C2=uint8(C2);
U3p=U3(:,1:p);
V3p=V3(:,1:p);
S3p=diag(S3(1:p,1:p));
C3=U3p * diag(S3p) * V3p';
C3=255*C3;
C3=uint8(C3);
Q(:,:,1)=C1;
Q(:,:,2)=C2;
Q(:,:,3)=C3;
// finding size and error of the compressed image
sz = (3*(numel(U1p) + numel(V1p) + numel(S1p)));
err = mean( abs(B(:)-Q(:)) );
toc;
t=toc;
// displying the compressed image
figure,
imshow(Q)
title( sprintf('p=%d, size=%d,err=%d', p, sz,err) );请告诉我如何保存压缩文件。
发布于 2013-03-13 16:56:29
正如Lorem Ipsum所说,您不是在保存压缩图像,而是在保存重建图像本身(Q是重建图像)。所以你的文件将和一个8位的.bmp文件一样大。
即使您只保存相关的奇异值和相应的向量,这也不应小于.jpg格式的文件。Jpeg已经被高度压缩,拥有比奇异值分解更好的技术...
https://stackoverflow.com/questions/15377437
复制相似问题