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

这里使用了双线性插值来处理极坐标转换到直角坐标下时，点不对称的问题，不过最后生成的图有一点问题，本来应该是一圈一圈的圆环，但是转化后的只有部分圆环了 rMin = 0.1 rMax =0.9 im2 = imread( ivus.bmp ) figure(4) imshow(im2) imR2 = PolarToIm(im2, rMin, rMax, 300, 300) figure(5) imshow(imR2, [0 255]) 然后直接调用附件的funct

Here is a code for gray to RGB conversion. You must use gray image for im1 & RGB image for im2