现在来将之前2个像素成像的例子扩展到6个像素,看看2维傅里叶变换是如何解析出每个像素各自贡献的信号。下图所示是一幅脑扫描像,在中央标识了2x3的像素阵列(A-F),以及频率编码和相位编码方向的定义。
在成像执行过程中,我们仍然进行2次相位编码步。相编步1时,相编梯度关闭,则6个体素的正弦波信号均保持了同相位。而相编步2时,我们假设第2排的像素(B,D,F)相对于第1排,都积累到了相位增益量180度。并在每一相编步后,都采集一次合成的磁共振信号,是6个像素一次总的贡献。那么也就是说,这里对于一个断层面,2次相编步就会采集两个磁共振信号。
相编步1采到的信号是对应于两排像素间没有相位差时的,对该信号做傅里叶变换会产生一个频谱结果,包含有3个频率成分(w1,w2,w3),各自的幅值为(A+B,C+D,E+F),正是频率编码方向上每个列的像素信号值之和。到这儿,只是知道每频率分量列上像素的值之和,具体各个像素的值有待进一步解析。
当相编步2后又采到同一层面6个像素的一次和信号,如前所述,此中的第2排每像素都增益了180度相位,再做傅里叶变换频谱也随之发生改变。同样地,解析出(w1,w2,w3)3个频率成分,但这次3个的幅值分别变成了(A-B)、(C-D)和(E-F)。还是同样的处境,只靠单步的测量信号采集,不能确定出每个像素的各自贡献。步1是两排像素间的和值,步2是差值。
由傅里叶变换已经解析出两次叠加采集信号的频谱,有3个频率成分/列,现在综合起来、各个像素的相位信息就能够计算得到。像之前文章中2像素的例子一样,用简单的代数运算基于2步的磁共振叠加信号,像素值可以一一确定。
(S0+S1)=(A+B)+(C+D)+(E+F)+(A-B)+(C-D)+(E-F)=2(A+C+E)
(S0-S1)=(A+B)+(C+D)+(E+F)-(A-B)-(C-D)-(E-F)=2(B+D+F)
(注,真实的正交解调应该是对任意的同频正弦波S0=A0sin(wt+phi0),解析得到A0和phi0共同决定像素值S0)
一幅扫描的磁共振图像可能有几百个像素点,因此需要执行上百步的相位编码梯度,但是背后的基本重建原理都是一样的。
http://mriquestions.com