傅里叶成像与图像重建.ppt

傅里叶成像与图像重建 纽迈科技 f Amplitude A: 261.63Hz B: 392Hz C: 523.25Hz 傅里叶变换 A B C 傅里叶变换 时域中无法区分的信号可以利用频率域区分 在脉冲NMR中，90°脉冲后，FID信号产生出来，但不检测，此时加宽度为t1的脉冲梯度Gx，在Gx作用下，t1时间段内，FID信号自然发展，t1结束时，关闭Gx，马上打开Gz，Gz将持续t2时间。在t2开始时开始检测FID信号，相等间隔取样，N个取样点。要建N×N的图像需要采到N ×N 的数据矩阵。因此需要采N个不同的FID信号，实验需要重复N次，每次改变Gx脉冲强度，保持其他参数不变。如此采集到N ×N的数据在时域表象中形成等距的规格，如下图。 傅里叶成像方案 采样点在时域表象中的分布 这里二维时域是一种虚构，t1，t2只是两个真实的时间段，在形式上画成正交二维时域的形式。时域表象也可变换到K空间表象： Kx=?gxt1 Kz=?gzt2 测量的Ｎ个FID可用S(t1, t2)表示，对其进行二维傅里叶变化可得到S(ωx, ωz)，由下式得到自旋密度ρ(x, z)的分布图像： X=?x/?gx Z= ?z/?gz 时域S(t1, t2) 混合域S(t1, ω2) 频域S(ω1, ω2) F F-1 F F2 F-2 F-1 时域、混合域和频域相互间的关系，F表示正傅里叶变换，F-1表示反傅里叶变换，F2及F-2表示正、反二维傅里叶变换 S(t1, t2)经过一维傅里叶变换后得到混合域上的信号S(t1, ω2)这相当于一维NMR谱即一维频域信号，但由于t1不同，信号幅度和相位受到了不同的调制，具体说，采集到的Ｎ个FID信号都包含频率成分 ωz=γzgz 所以经过一维傅里叶变换后，可显示出ωz谱，而在t2后FID信号具有不同的初相位： Φx= γxgxt1 可见gxt1影响被检测信号的初相位。初相位Φx是ｘ的函数，沿ｘ轴不同位置核自旋在旋转坐标系中具有不同的初相位。在固定ｘ位置看， Φx又是gx的函数，Ｎ个FID对应不同的t1，可见ｘ轴位置信息已包含在FID相位中，因此t1称为相位编码时间，gx称为相位编码梯度，gz称为频率编码梯度，t2称为频率编码时间。 如下图中第一个FID ，由于t1=0，采集到的FID是完整的，一维傅里叶变换后，给出最大幅度的两个峰，由于初相位为0，两个峰几乎等高，随FID序号的增加，峰高减小，由于相位调制，峰型也有相应的变化。 对S(t1, ω2)再做一次傅里叶变化，就得到S(ω1, ω2) ，再得到ρ(x, z)。 在线性梯度gx，gz作用下的部分FID信号及其一维傅里叶变换，左边号码是原FID中的编号 为了应用离散快速傅里叶变化FFT，取样数Ｎ，相位编码步数Ｍ通常选择为２的整数次幂，为了得到优良的傅里叶变换表象，一个常用的步骤是在被变换的阵列中添零。 Kumar等人于1975年提出原始的傅里叶成像，到1980年Edelstein等人提出一个重要修改，克服了原始傅里叶成像方法的缺点，被普遍接受流传至今。修改主要是将修改t1为固定时间段τ，而让gx为步进变量。第n步相位编码为 Φxn=ωxnt1=γgxnxt1=γτxn*Δgx 修改后，弛豫衰减对每次扫描都保持相同；修改后梯度磁场暂态非线性对每次扫描近似相同，从而使相位误差减小。 ρ(x,z); T1 weighted (x,z); T2 weighted (x,z) … S(t1,t2)→S(ω1,ω2) →S(x,z) 二维傅里叶变换 梯度空间编码 谢 谢！ 上海市普陀区中江路889号-1517 （200333） 苏州高新区苏州科技城苏高新软件园2-101 （215163） E-mail: support@ Tel: 0218008；0512