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(A

clc close all clear all I=imread( desert.jpg ) [r c]=size(I) bs=16 Block Size (8x8) nob=(r/bs)*(c/bs) Total number of 8x8 Blocks Dividing the image into 8x8 Blocks kk=0 for i=1:(r/bs) for j=1:(c/bs) Block(:,:,kk+j)=I