MRI图像重建，添加所有列切片

Question

我有一个包含32个图像集的MRI数据集。 现在，我设法逐列重建。 如何全部添加 或其他任何重建图像的好方法。

源文件：https://gist.github.com/saran08/1a85237851b049c9a7f2b78a55fc37d9

Nx = 64; % Resolution (assumes equal phase and frequency resolutions) 
data = readtable('5data.txt','ReadVariableNames',false); 
data = table2array(data(2:end,2:3)); 
data = strrep(data,',','.'); 
data = str2double(data); 
data = data(:,1) + i*data(:,2); 
data = vec2mat(data,Nx); 
% Data is now in an Nx X Nx array 'data' of complex numbers 


for i=1:32 
s= data(:,i); 
im(:,i)=fftshift(ifft2(ifftshift(squeeze(s(:,:))))); 
im(:,i)=abs(im(:,i)); 
figure((i));imshow(im(:,i),[]),title(['reconstructed Image' num2str(i)]) 
end