u can calculate DCT

clc

inp = imread( Im1.jpg )



inp1 = imread( Im2.jpg )



A = double(inp(:,:,1))





B = double(inp1(:,:,1))





A1=double(blkproc(A,[8 8], dct2 ))

B1=double(blkproc(B,[8 8], dct2 ))



[r,c] = size(A1)

r1 = []

for i=1:r

for j=1:c



if A1(i,j) > B1(i,j)

r1(i,j)=A1(i,j)

else

r1(i,j)=B1(i,j)

end



end

end



r1= double(r1(:,:,1))

C=blkproc(r1,[8 8], idct2 )

round(C)

C1=idct2(r1)



figure

imshow(uint8(C))

figure

imshow(uint8(C1))

imwrite(uint8(C), dct8x8min.jpg ) -u can calculate DCT

clc

inp = imread( Im1.jpg )



inp1 = imread( Im2.jpg )



A = double(inp(:,:,1))





B = double(inp1(:,:,1))





A1=double(blkproc(A,[8 8], dct2 ))

B1=double(blkproc(B,[8 8], dct2 ))



[r,c] = size(A1)

r1 = []

for i=1:r

　　for j=1:c

　　　

　　　 if A1(i,j) > B1(i,j)

　　　　　 r1(i,j)=A1(i,j)

　　　 else

　　　　　　r1(i,j)=B1(i,j)

　　　 end

　　　

　　end

end



r1= double(r1(:,:,1))

C=blkproc(r1,[8 8], idct2 )

round(C)

C1=idct2(r1)



figure

imshow(uint8(C))

figure

imshow(uint8(C1))

imwrite(uint8(C), dct8x8min.jpg )