fft_reconstruct.m

用2D_FFT方法重建核磁共振的图像

%%Reconstruct MRI image By 2D FFT Method

%%%%%%%%%%%%%Tian hongjun  09.14.2005 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [img]=FFT_RENCONSTRUCT(B,NewM,NewN)
[M,N]=size(B);


c=fft2(B);
%c= ifft2(B);

%%c = abs(c);

  
 c = [c(1:M,N/2+1:N),c(1:M,1:N/2)];%这里中间必须为逗号！
 %%c = [c(M/2+1:M,1:N);c(1:M/2,1:N)];%%%%%%%%如果中间的分号改成豆好则出错!!!!!!!!!!!!!

for m = 1 : NewM
    for n = 1 : NewN
        img(m,n) = c(m,n+(N-NewN)/2);
    end
end