本文主要实现多通道脑部MRI图片的加速重建,使用GRAPPA方法。
目录
加载满采数据
数据欠采样
GRAPPA重建
完整数据代码下载
加载满采数据
load brain_8ch
DATA = DATA/max(max(max(abs(ifft2c(DATA))))) + eps;
数据欠采样
mask=ones(sy,sx);
for i=1:2:sy
if 94<i && i<106
continue
end
mask(i,:)=zeros(sx,1);
end
mask = repmat(mask,[1,1,8]);
GRAPPA重建
disp('Performing GRAPPA reconstruction ... slow in Matlab! ')
tic; reskGRAPPA = GRAPPA(DATAc,calib,ksize,0.01);toc
resGRAPPA = ifft2c(reskGRAPPA);
完整数据代码下载
MATLAB -传统 GRAPPA MRI 重建 完整代码和数据