9/13/2021 2-DCT Image compression - MATLAB Answers - MATLAB Central

MATLAB Answers

2-DCT Image compression

 Follow 108 views (last 30 days) Show older comments

Amarnath R
on 3 Nov 2018  Vote 0
Commented: nor saziana ariani sazali
on 25 Jun 2021

Problem:

I tried implementing Discrete Cosine Transformation compression using matlab. Input image would a
jpg image (Lena) having a size 512 X 512.

There are two stages namely compression and decompression.

Compression and Quantization:

The input image is converted to YCbCr component. Then Y component is taken up for compression.
Further DCT will quantized.

Quantization and Decompression:

The quantized image is undergoes dequantization for decompression.

Issues:

Rectangular boxes are spotted in the decompressed version of the image. Is anything wrong with the
code? For your inference, below are the sample input and output images and followed by the matlab
code.

Input image:

https://in.mathworks.com/matlabcentral/answers/427708-2-dct-image-compression 1/9
9/13/2021 2-DCT Image compression - MATLAB Answers - MATLAB Central

MATLAB Answers

Y Component in YCbCr:

https://in.mathworks.com/matlabcentral/answers/427708-2-dct-image-compression 2/9
9/13/2021 2-DCT Image compression - MATLAB Answers - MATLAB Central

MATLAB Answers

Output image:

https://in.mathworks.com/matlabcentral/answers/427708-2-dct-image-compression 3/9
9/13/2021 2-DCT Image compression - MATLAB Answers - MATLAB Central

MATLAB Answers

Code:

 
clc;
clear all;
close all;
I = imread('lena512.jpg');
figure, imshow(I);
% Y = I;
YCbCr = rgb2ycbcr(I);
figure, imshow(YCbCr);
Y = YCbCr(:,:, 1);
figure, imshow(Y);
[h, w] = size(Y);
r = h/8;
c = w/8;

s = 1;
q50 = [16 11 10 16 24 40 51 61;
12 12 14 19 26 58 60 55;
14 13 16 24 40 57 69 56;
14 17 22 29 51 87 80 62;
18 22 37 56 68 109 103 77;
24 35 55 64 81 104 113 92;
49 64 78 87 103 121 120 101;

https://in.mathworks.com/matlabcentral/answers/427708-2-dct-image-compression 4/9
9/13/2021 2-DCT Image compression - MATLAB Answers - MATLAB Central

72 92 95 98 112 100 103 99];

MATLAB Answers
% COMPRESSION
for i=1:r
e = 1;
for j=1:c
block = Y(s:s+7,e:e+7);
cent = double(block) - 128;
for m=1:8
for n=1:8
if m == 1
u = 1/sqrt(8);
else
u = sqrt(2/8);
end
if n == 1
v = 1/sqrt(8);
else
v = sqrt(2/8);
end
comp = 0;
for x=1:8
for y=1:8
comp = comp + cent(x, y)*(cos((((2*(x-1))+1)*(m-1)*pi)/16))
end
end

F(m, n) = v*u*comp;
end
end

for x=1:8
for y=1:8
cq(x, y) = round(F(x, y)/q50(x, y));
end
end
Q(s:s+7,e:e+7) = cq;
e = e + 8;
end
s = s + 8;
end

% % % % % % % % % % % % % % %
% % DECOMPRESSION
% % % % % % %
s = 1;
for i=1:r
e = 1;
for j=1:c
cq = Q(s:s+7,e:e+7);
https://in.mathworks.com/matlabcentral/answers/427708-2-dct-image-compression 5/9
9/13/2021 2-DCT Image compression - MATLAB Answers - MATLAB Central

for x=1:8
for y=1:8
MATLAB Answers
DQ(x, y) = q50(x, y)*cq(x, y);
end
end
for m=1:8
for n=1:8
if m == 1
u = 1/sqrt(8);
else
u = sqrt(2/8);
end
if n == 1
v = 1/sqrt(8);
else
v = sqrt(2/8);
end
comp = 0;
for x=1:8
for y=1:8
comp = comp + u*v*DQ(x, y)*(cos((((2*(x-1))+1)*(m-1)*pi)/16
end
end

bf(m, n) = round(comp)+128;
end

end

Org(s:s+7,e:e+7) = bf;
e = e + 8;
end
s = s + 8;
end

imwrite(Y, 'F:\workouts\phd\jpeg\input.jpg');
imwrite(uint8(Org), 'F:\workouts\phd\jpeg\output.jpg');
return;

   0 Comments
Sign in to comment.

Sign in to answer this question.


Answers (1)

https://in.mathworks.com/matlabcentral/answers/427708-2-dct-image-compression 6/9
9/13/2021 2-DCT Image compression - MATLAB Answers - MATLAB Central

Pham Minh Hanh

on 26 Mar 2021  Vote 0  Link
MATLAB Answers
Edited: Pham Minh Hanh
on 26 Mar 2021

Your decompression code is not right. IDCT is defined by:

so that u and v have to be inside the 2nd loop, not outside. Here is my code. The right image
is the input image, the left one is output. The decompressed image is similar to the input one,
so DCT works well as expected.

 
for x = 1:8
for y = 1:8
comp = 0;
for m = 1:8
for n = 1:8
if m == 1
u = 1/sqrt(2);
else
u = 1;
end
if n == 1
v = 1/sqrt(2);
https://in.mathworks.com/matlabcentral/answers/427708-2-dct-image-compression 7/9
9/13/2021 2-DCT Image compression - MATLAB Answers - MATLAB Central

else
v = 1; MATLAB Answers
end

comp = comp + u*v*DQ(m, n)*(cos((((2*(x-1))+1)*(m-1)*pi

end
end
bf(x, y) = round((1/4) *comp + 128);
end

   5 Comments Show 4 older comments

nor saziana ariani sazali

on 25 Jun 2021 

ohh i see..thankyou for explanation Walter

Sign in to comment.

Sign in to answer this question.

See Also
MATLAB Answers
image compression should bitmap image?
1 Answer

Why is my image not recombining properly after a DCT?

2 Answers

good evening , iam getting error as (im2jpeg(x, quality) Error using im2jpeg (line 21) The input must be a UINT8 image) at
line ...
1 Answer

Entire Website
Obfuscated MATLAB Code
Blogs

RGB2YCC(seq)
File Exchange

dct(y)
File Exchange

Categories
Image Processing and Computer Vision
> Image Processing Toolbox
> Image Filtering and Enhancement

Tags
dct discrete cosine tr... jpeg compression decompression

Products
MATLAB

https://in.mathworks.com/matlabcentral/answers/427708-2-dct-image-compression 8/9
9/13/2021 2-DCT Image compression - MATLAB Answers - MATLAB Central

Release
 MATLAB Answers
R2014a

Poll

What’s your beverage of choice while programming in MATLAB?

Water

Tea

Coffee

Juice

Something harder

Nothing, I'm good

11962 votes
26 comments
2863 views

Vote

mathworks.com
© 1994-2021 The MathWorks, Inc. MATLAB and Simulink are registered trademarks of The MathWorks, Inc. See mathworks.com/trademarks for a list of additional
trademarks. Other product or brand names may be trademarks or registered trademarks of their respective holders.

https://in.mathworks.com/matlabcentral/answers/427708-2-dct-image-compression 9/9