

ksp1= importdata('MRI_kspace_data.txt');

ksp2= ksp1(:,1) + ksp1(:,2)*j;

ksp = reshape(ksp2,256,256);

figure(2)
image(1e6*abs(ksp))
colormap(gray)


figure(3)
imagesc(angle(ksp))
colormap(gray)

figure(4)
image(1e6*real(ksp))
colormap(gray)

figure(5)
image(1e6*imag(ksp))
colormap(gray)



im= fftshift( ifft2(ksp) );
figure(6)
imagesc(imag(im))
colormap(gray)


