【核磁共振成像】相位差重建

简介: 【核磁共振成像】相位差重建

一、相位差map重建一般步骤和反正切函数主值范围

  MRI是一个相敏成像模态MR原始数据傅里叶变换后的复数图像中每个像素值有模和相位。标准模重建是遗弃相位,只用数值建立图像。然而,编码在相位里有有用的信息

  比如相位map可产生B0均匀度信息,可用于匀场;还可产生流体流动信息,可用于相敏法血管造影;相位信息还可用于体内MR温度分布测量以及MR弹性造影。

  采集两个独立的数据集,然后形成一个相位差map,就 可以解决相位map污染问题。而所谓的相位差map,就是用一对相位像的像素逐个相减得到的,目的是突出所希望的相位,同时对消有害的相位

  从相位像中提取物理参数值不是一个简单事情。因为 很多因素对相位像贡献了不需要的相位,例如梯度涡流化学位移磁化率变化以及伴随场等等。这些不需要的贡献使得相位map受到污染,难以解释,所希望的信息往往被这些有害相位淹没。

   相位差重建 的核心是对每个体素进行反正切运算来产生相位图像(map),反正切函数的输出被定义在一个有限的范围内,称为主值范围。在MRI领域,相位范围一般通过一个四象限反正切函数而扩展到(-π,π]。主值范围以外的相位值被主值内的一个值代表,混叠(alias)。相位map中混叠伴随着不连续性(在±π跃变),叫做相位卷绕

  为提高计算效率,并最小化相位卷绕数,一个最佳相位差重建应该每像素只做一次反正切运算。还有,因为不连续性,希望在计算反正切之前先执行比如相位阵列多线圈数据组合运算和伴随场相位校正运算。

相位差重建步骤流程


二、反正切运算

  假定已经采集了两个K空间数据,想产生一个相位差map。数据集首先符号交替并填零,然后分别进行傅里叶变换,产生两个独立的复数像

  恢复范围(-π,π]内的相位差值到其真实值的过程称为 相位解卷绕

  因为涡流式依赖于MR硬件的设计和校准的,通常难以预期产生的相位误差的准确的空间依赖。

相位卷绕


三、可预期相位误差和伴随场的校正

  相位差像形成后,还经常残留着有害的相位误差。如由梯度涡流造成的相位污染,给相位对比度血管造影造成麻烦。因此此类相位误差经常靠实验拟合来校正

  例如,在应该有零相位差的某个区域,比如在一个相位对比度血管造影图中的静止组织,一个恒定相位和线性相位可用相位差像的多项式拟合来确定。因为拟合的相位大概是全部由系统不完善譬如涡流造成的,用后处理可用将其消除。

  其他相位误差,比如由伴随场产生的那些相位误差,可以精确预测,因为它们是基础物理效应(可以由麦克斯韦方程描写),在计算反正切之前校正这些相位误差是有利的

  在反正切运算之前分别应用图像变形校正到实像和虚像的另一个优点是:被变形校正算法引起的图像强度的任何变换都不会影响到相位差map


四、图形变形校正

  由梯度非线性造成的图像变形的校正运算也可用于相位差像。对于多线圈组合,在反正切运算之前应用图像变形校正更是可取的。这是因为图像变形校正是用内插方法比如立方样条的正投影映射。如果在反正切运算之后应用图像变形校正,将会遭遇相位卷绕,并且在图像中其跳变的不连续还会引起多余的上冲。

  内插相位像的一个缺陷:如果相位像被内插(例如图像变形校正或矩形视野缩小),在混叠边界处会发生上冲(箭头),这种内插过的相位像就不适合解卷绕处理更可取的方法是在反正切之前按内插实像和虚像,而不要内插相位像


相关文章
|
编解码 算法
【核磁共振成像】部分傅里叶重建
【核磁共振成像】部分傅里叶重建
|
数据采集 存储
【核磁共振成像】观共享重建
【核磁共振成像】观共享重建
|
6月前
R语言分析蛋白质组学数据:飞行时间质谱(MALDI-TOF)法、峰值检测、多光谱比较
R语言分析蛋白质组学数据:飞行时间质谱(MALDI-TOF)法、峰值检测、多光谱比较
|
6月前
|
机器学习/深度学习 算法 前端开发
高速公路表面图像裂缝检测程序
高速公路表面图像裂缝检测程序
|
6月前
|
编解码 算法 安全
图像识别在轻轨轨道的损坏程度的分析
图像识别在轻轨轨道的损坏程度的分析
54 0
|
6月前
|
机器人
[贴装专题] 基于多目视觉的手眼标定
[贴装专题] 基于多目视觉的手眼标定
66 0
|
机器学习/深度学习 数据采集 编解码
【核磁共振成像】并行采集MRI
【核磁共振成像】并行采集MRI
【核磁共振成像】并行采集MRI
|
传感器 机器学习/深度学习 算法
【姿态解算】基于扩展卡尔曼滤波九轴传感器姿态解算研究附代码
【姿态解算】基于扩展卡尔曼滤波九轴传感器姿态解算研究附代码
|
存储 数据采集 编解码
【核磁共振成像】傅里叶重建
【核磁共振成像】傅里叶重建
|
存储 算法 NoSQL
【核磁共振成像】方格化重建
【核磁共振成像】方格化重建