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

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

一、部分傅里叶重建

  在部分傅里叶采集中,数据并不是绕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空间)算法处理


相关文章
|
存储 缓存 安全
C++数组全解析:从基础知识到高级应用,领略数组的魅力与技巧
C++数组全解析:从基础知识到高级应用,领略数组的魅力与技巧
589 1
|
机器学习/深度学习 自然语言处理 PyTorch
PyTorch应用实战六:利用LSTM实现文本情感分类
PyTorch应用实战六:利用LSTM实现文本情感分类
631 0
|
安全 Linux Shell
structure needs cleaning结构需要清理解决方案
structure needs cleaning结构需要清理解决方案
4789 0
|
机器学习/深度学习 传感器 算法
【虚警检测】基于CA-CFAR(单元平均恒虚警)、GO-CFAR(最大选择恒虚警)、SO-CFAR(最小选择恒虚警)算法的恒虚警检测附Matlab代码
【虚警检测】基于CA-CFAR(单元平均恒虚警)、GO-CFAR(最大选择恒虚警)、SO-CFAR(最小选择恒虚警)算法的恒虚警检测附Matlab代码
|
11月前
|
缓存 前端开发 Android开发
安卓开发中的自定义视图:从零到英雄
【10月更文挑战第42天】 在安卓的世界里,自定义视图是一块画布,让开发者能够绘制出独一无二的界面体验。本文将带你走进自定义视图的大门,通过深入浅出的方式,让你从零基础到能够独立设计并实现复杂的自定义组件。我们将探索自定义视图的核心概念、实现步骤,以及如何优化你的视图以提高性能和兼容性。准备好了吗?让我们开始这段创造性的旅程吧!
167 1
|
6月前
|
安全 Java 程序员
巧用Optional之优雅规避NPE问题
本文探讨了Java中常见的NullPointerException问题及其解决方案,重点介绍了Optional类的使用。通过实例代码分析,展示了如何用Optional替代传统的空值检查,使代码更简洁、优雅。文章详细讲解了Optional的创建方法(如of、ofNullable、empty)及常用方法(如get、orElse、map、flatMap、filter),并通过实战案例演示了其在实际开发中的应用,帮助开发者有效避免NPE问题,提升代码质量。
157 2
巧用Optional之优雅规避NPE问题
|
搜索推荐 Linux Android开发
如何根据自己的开发板型号下载和配置交叉编译链
【8月更文挑战第24天】本指南详细介绍了为特定开发板下载及配置交叉编译链的过程。首先,需明确开发板型号与架构,通过查阅文档了解其处理器架构和支持的操作系统。其次,根据开发板架构及目标操作系统确定所需的交叉编译链类型。下载环节推荐三种途径:在线搜索、访问官方站点以及开源社区。安装阶段涉及解压文件并设置环境变量,以确保能在终端直接调用交叉编译工具。最后,通过检查版本信息及编译测试程序验证交叉编译链是否安装正确。整个过程中应注意选择合适的版本、遵循安装指导并妥善处理遇到的问题。
302 3
|
Linux pouch 容器
CentOS7部署阿里巴巴开源的pouch容器管理工具实战
关于如何在CentOS 7.6操作系统上安装和使用阿里巴巴开源的Pouch容器管理工具的实战教程。
315 2
CentOS7部署阿里巴巴开源的pouch容器管理工具实战
Anaconda——添加清华源
Anaconda——添加清华源
593 0
|
机器学习/深度学习 存储 算法
【类脑计算】突触可塑性模型之Hebbian学习规则和STDP
本文介绍了突触可塑性中的Hebbian学习规则和STDP(Spike-Timing Dependent Plasticity),两种基于神经元活动调节突触强度的机制,其中Hebbian规则强调同时活动的神经元间的连接增强,而STDP则考虑了脉冲时间差异对突触强度的调节作用。
919 2