我正在使用奇异值分解技术进行图像压缩。 我在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) );
请帮我如何保存压缩文件。
答案 0 :(得分:1)
正如Lorem Ipsum所说,你没有保存压缩图像,你正在保存重建图像本身(Q是重建图像)。因此,您的文件将与8位.bmp文件一样大......
即使您只保存相关的奇异值和相应的向量,也不应小于.jpg格式文件。 Jpeg已经被大量压缩,技术比SVD更好......