output_of_DCT_in_image_fusion

代码说明：

  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 )