www.gusucode.com > Diabetic retinopathy > get_convert_matrix.m
function a = get_convert_matrix(M,N) a=zeros(M*N,M*N); B2=cosine_dct_basis(M*N); B=compute_cosine_modified_basis(M,N); a=B2*inv(B); for u=1:M*N for j=1:M*N if(((u-1)==0) && (mod((j-1),N)~=0)) a(u,j)=(1/sqrt(2))*a(u,j); end if(((u-1)~=0) && (mod((j-1),N)==0)) a(u,j)=sqrt(2)*a(u,j); end a(u,j)=(1/sqrt(M))*a(u,j); end end return