matlab mri的k空间,磁共振K空间解析

关键词 ：磁共振 MRI K空间 matlab

以上是我们要用的工具，一台电脑，matlab软件(为了方便起见且单位有授权，用了7.0版本)，一张小姐姐的图像。

假设我们有一台B0和B1都结对均匀的磁共振，且扫描物中只有水，且在TE=inphase的时间扫描(以上所有的设定是为了假设扫描出的图像的相位角都为0)，扫描物是上面的小姐姐，那么我们就得到了小姐姐的灰白图像。

然后我们往回推，我们把这张灰度图进行FFT转换(FFT和IFFT得到的结果都是一样的)，这样我们就得到了在图像重建之前和扫描之后的K空间数据。

i=imread('aql.bmp');

figure(1);

imshow(i);

title('安琪拉');

t=rgb2gray(i);

j=fft2(t);%进行FFT变换

f=fftshift(j);%进行K空间的修正(可忽略)

f=log(abs(f)+1);%取模值

K空间的像素中填充的都是向量，我们人类可能看到四维图，因此我们取K空间的模值(绝对值)，就得到了一个可视的图像。

以往的教材中经常出现下面一张图，表示K空间可转换为图像。

但模值不能代表K空间的真实状态，取模值的时候损失了向量的夹角(相位)，因此使用能量谱不能变换出想要的图像来，因此以往的教材中的示意图只能用于示意，而非真正的k空间。

K空间中的单元都是向量，高中数学在学复数时，老师告诉我们，向量都可以用

表示，m代表长度，θ代表向量转动的角度，在MRI中，θ是和梯度的施加有关的，θ=γGt,  γ是常数，G是所处位置的梯度强度，t是梯度施加时间。

如果把K空间中的单元进行实部(频率编码)和虚部(相位编码)两个方向进行二次积分，便得到了一个意想不到的答案。

上图是我们得到的答案

这个是百度百科里傅里叶变换的公式，我们观察我们自己得到的公式，

如果把里面的变量表换一下，就是这个公式。

以上。

这个还是几年前做的研究，发出来以便能让搜索引擎搜到，让想探索磁共振成像原理的人少走一些弯路。

笔者将在后面推出初代Dixon的算法详解，敬请期待。