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

最新推荐文章于 2024-10-26 22:57:55 发布

CLARKYY

最新推荐文章于 2024-10-26 22:57:55 发布

阅读量501

点赞数 1

文章标签： matlab

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

版权

该博客详细介绍了使用MATLAB进行2D傅立叶变换的过程。通过加载'magni2D'和'phase2D'数据，计算复数乘法生成k空间数据，然后应用快速傅立叶变换(fft2c)。结果存储在变量ksp中，并保存为'ksp.mat'。此外，还展示了如何从k空间数据重构图像并比较原始图像与重构图像的幅度和相位差异。

摘要由CSDN通过智能技术生成

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