【核磁共振成像】部分傅里叶重建

简介: 【核磁共振成像】部分傅里叶重建

一、部分傅里叶重建

  在部分傅里叶采集中,数据并不是绕K空间中心对称收集的,而是K空间的一半是完全填充的,另一半只收集了一小部分数据

  部分傅里叶采集所依据的原理是如果物体是实的(real),其傅里叶变换是厄米的,意味着绕K空间中心实部是对称的,虚部是反对称的

  在频率相位方向部分傅里叶采集的主要优点分别降低了回波时间(TE)和扫描时间。在频率编码方向,部分傅里叶采集也减少了沿该轴的梯度矩,从而降低了流动和运动伪影。由K空间最大范围决定的空间分辨率不受影响,与全K空间采集情况等价,只是SNR有所降低,也可能引进一些图像伪影。

  然而,由诸多因素例如运动、共振频率偏移、硬件群时延、涡流、接受线圈B1场不均匀等引起的有害相移,使被建图像失真。因此,在图中显示的未完全填充的半K空间中的额外数据(有时叫”过扫描数据”)被用来克服这些问题。

  部分傅里叶采集也可用于3D扫描的层面选择方向(即第二个相位编码方向),但不普遍。因为层面的相位编码步数一般不是很大,考虑到还需要一些过扫描数据,使得时间的节省并不特别显著

Kx和Ky分别代表频率编码和相位编码方向

  (a)是沿频率编码方向覆盖部分K空间(部分Kx或部分回波)

  (b)是沿相位编码方向覆盖部分K空间(部分Ky或部分Nex)


二、部分傅里叶重建算法

2.1 填零

  填零就是 未测量的K空间数据用零来代替,然后用传统重建以得到图像。对于全K空间采集原始数据填零被普遍用来内插图像并降低部分体积效应;对于部分傅里叶采集填零可以用来代替未测量的数据,如果希望图像内插还可以补额外的零。填零之后,可用基于标准傅里叶变换的全K空间重建

  填零通常导致靠近锐利边缘处有一些Gibbs跳动伪影,是 由于K空间数据截断引起的填零的优点是能给出物体低空间频率过扫描范围内比较真实的图像。对此低空间频率范围,相位信息被保留。因此,大结构的相位是精确的,允许填零由于相敏重建。合理的相位精度通常需要采集比较高分数的K空间,例如0.75或更高。

Kx和Ky分别代表频率编码和相位编码方向

  (a)是沿频率编码方向覆盖部分K空间(部分Kx或部分回波)

  (b)是沿相位编码方向覆盖部分K空间(部分Ky或部分Nex)


2.2 零差处理

  零差处理从自身数据产生的低空间频率相位map去校正由不完全K空间数据重建产生的相位误差零差处理利用的是K空间数据的厄米共轭对称性。零差不适合需要图像相位的那些应用,如匀场、相位对比度(用相位差重建的)和相敏热成像;然而,如果复数差运算在K空间执行,则复数差相位对比度可以用零差成像

零差高通滤波器

  电感阻止高频信号通过而允许低频信号通过,电容的特性却相反。信号能够通过电感的滤波器、或者通过电容连接到地的滤波器对于低频信号的衰减要比高频信号小,称为低通滤波器

  最简单的高通滤波器是“一阶高通滤波器”,它的的特性一般用一阶线性微分方程表示,它的左边与一阶低通滤波器完全相同,仅右边是激励源的导数而不是激励源本身当较低的频率通过该系统时,没有或几乎没有什么输出,而当较高的频率通过该系统时,将会受到较小的衰减。实际上,对于极高的频率而言,电容器相当于“短路”一样,这些频率,基本上都可以在电阻两端获得输出。换言之,这个系统适宜于通过高频率而对低频率有较大的阻碍作用,是一个最简单的“高通滤波器”。

  零差方法的缺点是 方程中用的低频相位map不能准确地描写急速变换的相位。为解决此问题,可用迭代部分傅里叶方法。该方法 用零差重建估计模像,同时从低频相位map估计相位。结合估计的模像和相位像给出一个复数像,对此复数像经傅里叶变换可得到估计的K空间数据。

平滑混合K空间所有点信号值的估计,迭代零差重建由方程给出的结合函数W(k)的曲线

  如果部分傅里叶采集用在两个正交方向,假定部分采集分数大于0.5,则一个方向可用零差重建处理,而另一个方向只能用填零处理。即如果部分傅里叶采集只用在一个方向,其他方向必须首先用正常(即全K空间)算法处理


相关文章
|
机器学习/深度学习 传感器 算法
【信号去噪】基于Sage-Husa自适应卡尔曼滤波器实现海浪磁场噪声抑制及海浪磁场噪声的产生附matlab代码
【信号去噪】基于Sage-Husa自适应卡尔曼滤波器实现海浪磁场噪声抑制及海浪磁场噪声的产生附matlab代码
【信号去噪】基于Sage-Husa自适应卡尔曼滤波器实现海浪磁场噪声抑制及海浪磁场噪声的产生附matlab代码
【核磁共振成像】相位差重建
【核磁共振成像】相位差重建
|
机器学习/深度学习 传感器 算法
使用显著性检测的可见光和红外图像的两尺度图像融合(Matlab代码实现)
使用显著性检测的可见光和红外图像的两尺度图像融合(Matlab代码实现)
|
存储 数据采集 编解码
【核磁共振成像】傅里叶重建
【核磁共振成像】傅里叶重建
|
存储 算法 NoSQL
【核磁共振成像】方格化重建
【核磁共振成像】方格化重建
|
机器学习/深度学习 传感器 算法
【红外与可见光图像融合】离散平稳小波变换域中基于离散余弦变换和局部空间频率的红外与视觉图像融合方法(Matlab代码实现)
【红外与可见光图像融合】离散平稳小波变换域中基于离散余弦变换和局部空间频率的红外与视觉图像融合方法(Matlab代码实现)
|
机器学习/深度学习 传感器 算法
【红外图像增强】基于引力和侧向抑制网络的红外图像增强模型(Matlab代码实现)
【红外图像增强】基于引力和侧向抑制网络的红外图像增强模型(Matlab代码实现)
|
编解码
超声非线性成像-谐波成像基本知识
超声成像因其相对较低的价格、无电磁辐射、无创伤、成像方式多样而在医学成像中得到了广泛的应用。早期的超声成像技术都是基于线性声学原理。然而在实际的医学超声成像应用中,声波在生物组织的非线性传播和造影剂的非线性振动能够产生非线性声信号。在接下来的报告中讲解了如何使用非线性信号中的谐波信号进行成像,以及相关的谐波信号提取方法。最后使用一些代表性的应用实例进行验证,得出谐波成像在分辨率和对比度上确实优于传统基波成像的结论。
353 0
|
机器学习/深度学习 传感器 编解码
【图像重建】基于小波变换图像分解重建(PSNR对比)附matlab代码
【图像重建】基于小波变换图像分解重建(PSNR对比)附matlab代码
|
编解码 算法 Java
基于Gabor-小波滤波深度图表面法线的特征提取算法【通过正常Gabor-小波的直方图进行2D或3D特征提取】研究(Matlab代码实现)
基于Gabor-小波滤波深度图表面法线的特征提取算法【通过正常Gabor-小波的直方图进行2D或3D特征提取】研究(Matlab代码实现)
115 0