1

我正在使用奇异值分解技术进行图像压缩。我在 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) );

请帮助我如何保存压缩文件。

4

1 回答 1

1

正如 Lorem Ipsum 所说，您不是在保存压缩图像，而是在保存重建图像本身（Q 是重建图像）。因此，您的文件将与 8 位 .bmp 文件一样大...

而且即使只保存相关的奇异值和对应的向量，它也不应该小于 .jpg 格式的文件。Jpeg已经被高度压缩，具有比 SVD 更好的技术......

于 2013-03-13T08:56:29.200 回答