由幅度和相位数据，计算K空间数据

最新推荐文章于 2024-05-28 10:41:10 发布

CLARKYY

最新推荐文章于 2024-05-28 10:41:10 发布

阅读量448

点赞数 1

文章标签： matlab

本文链接：https://blog.csdn.net/CLARKYY/article/details/123927890

版权

load magni2D.mat
load phase2D.mat

[x,y,z] = size(magni2D);
ksp = zeros(x,y,z);


for i=1:z
    m = magni2D(:,:,i);
    p = phase2D(:,:,i);
    x = m.*cos(p)+1j*m.*sin(p);%generate k-space data
    y = fft2c(x);
    ksp(:,:,i)=y;
end

% save ksp.mat ksp

% test
slice = 17;
mrec = abs(ifft2c(ksp(:,:,slice))); mraw = magni2D(:,:,slice);
figure;imshow(mrec,[]);
figure;imshow(mraw,[]);
figure;imshow(mrec-mraw,[]);
prec = angle(ifft2c(ksp(:,:,slice))); praw = phase2D(:,:,slice);
figure;imshow(prec,[]);
figure;imshow(praw,[]);
figure;imshow(prec-praw,[]);

fft2c.m and ifft2c.m 来自软件包：
http://people.eecs.berkeley.edu/~mlustig/Software.html