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

目录

  • 一、相位差map重建一般步骤和反正切函数主值范围
  • 二、反正切运算
  • 三、可预期相位误差和伴随场的校正
  • 四、图形变形校正


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

  MRI是一个相敏成像模态MR原始数据傅里叶变换后的复数图像中每个像素值有模和相位。标准模重建是遗弃相位,只用数值建立图像。然而,编码在相位里有有用的信息
  比如相位map可产生B0均匀度信息,可用于匀场;还可产生流体流动信息,可用于相敏法血管造影;相位信息还可用于体内MR温度分布测量以及MR弹性造影。

  采集两个独立的数据集,然后形成一个相位差map,就 可以解决相位map污染问题。而所谓的相位差map,就是用一对相位像的像素逐个相减得到的,目的是突出所希望的相位,同时对消有害的相位
  从相位像中提取物理参数值不是一个简单事情。因为 很多因素对相位像贡献了不需要的相位,例如梯度涡流化学位移磁化率变化以及伴随场等等。这些不需要的贡献使得相位map受到污染,难以解释,所希望的信息往往被这些有害相位淹没。

   相位差重建 的核心是对每个体素进行反正切运算来产生相位图像(map),反正切函数的输出被定义在一个有限的范围内,称为主值范围。在MRI领域,相位范围一般通过一个四象限反正切函数而扩展到(-π,π]。主值范围以外的相位值被主值内的一个值代表,混叠(alias)。相位map中混叠伴随着不连续性(在±π跃变),叫做相位卷绕
  为提高计算效率,并最小化相位卷绕数,一个最佳相位差重建应该每像素只做一次反正切运算。还有,因为不连续性,希望在计算反正切之前先执行比如相位阵列多线圈数据组合运算和伴随场相位校正运算。
在这里插入图片描述

相位差重建步骤流程

二、反正切运算

  假定已经采集了两个K空间数据,想产生一个相位差map。数据集首先符号交替并填零,然后分别进行傅里叶变换,产生两个独立的复数像
  恢复范围(-π,π]内的相位差值到其真实值的过程称为 相位解卷绕
  因为涡流式依赖于MR硬件的设计和校准的,通常难以预期产生的相位误差的准确的空间依赖。
在这里插入图片描述

相位卷绕

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

  相位差像形成后,还经常残留着有害的相位误差。如由梯度涡流造成的相位污染,给相位对比度血管造影造成麻烦。因此此类相位误差经常靠实验拟合来校正
  例如,在应该有零相位差的某个区域,比如在一个相位对比度血管造影图中的静止组织,一个恒定相位和线性相位可用相位差像的多项式拟合来确定。因为拟合的相位大概是全部由系统不完善譬如涡流造成的,用后处理可用将其消除。

  其他相位误差,比如由伴随场产生的那些相位误差,可以精确预测,因为它们是基础物理效应(可以由麦克斯韦方程描写),在计算反正切之前校正这些相位误差是有利的
  在反正切运算之前分别应用图像变形校正到实像和虚像的另一个优点是:被变形校正算法引起的图像强度的任何变换都不会影响到相位差map


四、图形变形校正

  由梯度非线性造成的图像变形的校正运算也可用于相位差像。对于多线圈组合,在反正切运算之前应用图像变形校正更是可取的。这是因为图像变形校正是用内插方法比如立方样条的正投影映射。如果在反正切运算之后应用图像变形校正,将会遭遇相位卷绕,并且在图像中其跳变的不连续还会引起多余的上冲。
在这里插入图片描述
  内插相位像的一个缺陷:如果相位像被内插(例如图像变形校正或矩形视野缩小),在混叠边界处会发生上冲(箭头),这种内插过的相位像就不适合解卷绕处理更可取的方法是在反正切之前按内插实像和虚像,而不要内插相位像

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.hqwc.cn/news/92203.html

如若内容造成侵权/违法违规/事实不符,请联系编程知识网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Spring与Mybatis整合aop整合pageHelper分页插件

前言 Spring与MyBatis整合的意义在于提供了一种结合优势的方式,以便更好地开发和管理持久层(数据库访问)代码。 这里也是总结了几点主要意义 简化配置:Spring与MyBatis整合后,可以通过Spring的配置文件来管理和配置M…

C语言之数组题

目录 1.使用函数实现数组操作 2.冒泡排序 3.三子棋 4.【一维数组】交换数组 5.扫雷 6.概念辨析tips 我又来了,今天是数组题,本人还在补军训真的热!🆗 1.使用函数实现数组操作 2.冒泡排序 3.三子棋 4.【一维数组】交换数组 …

自定义spring-boot-start的jar包被引用时,提示找不到bean

类似这个报错: 重点要看一下我们自定义的start包下的config配置 BeanConditionalOnProperty(prefix "file", value "iSenable", havingValue "true")public FileServiceTemplate fileServiceTemplate(){return new FileServiceTe…

app加固怎么做?

app加固是一种保护应用程序免受攻击的技术,可以防止未经授权的访问、数据泄露和其他安全问题。下面是一些app加固的最佳实践: 1.加密和数据保护:对应用程序的数据进行加密和保护,以防止未经授权的访问和数据泄露。您可以使用AES、…

【ES6】—【新特性】—Symbol详情

一、一种新的原始数据类型 定义:独一无二的字符串 二、 声明方式 1. 无描述声明 let s1 Symbol() let s2 Symbol() console.log(s1, s2) // Symbol() Symbol() console.log(s1 s2) // falsePS: Symbol 声明的值是独一无二的 2. 有描述的声明 let s1 Symb…

玩转 PI 系列-看起来像服务器的 ARM 开发板矩阵-Firefly Cluster Server

前言 基于我个人的工作内容和兴趣,想要在家里搞一套服务器集群,用于容器/K8s 等方案的测试验证。 考虑过使用二手服务器,比如 Dell R730, 还搞了一套配置清单,如下: Dell R7303.5 尺寸规格硬盘CPU: 2686v4*2 内存&a…

大语言模型之六- LLM之企业私有化部署

数据安全是每个公司不得不慎重对待的,为了提高生产力,降本增效又不得不接受新技术带来的工具,私有化部署对于公司还是非常有吸引力的。大语言模型这一工具结合公司的数据可以大大提高公司生产率。 私有化LLM需要处理的问题 企业内私有化LLM…

深入解析Kotlin类与对象:构造、伴生、单例全面剖析

前言 本篇文章将带您了解Kotlin编程中的重要概念:类及构造函数、访问修饰符、伴生对象和单例模式。就像搭积木一样,我们会逐步揭开这些概念的面纱,让您轻松理解它们的作用和用法。无论您是编程新手还是有经验的开发者,本文都将为…

GPU编程(基于Python和CUDA)(一)——零基础安装pycuda

系列文章目录 GPU编程(基于Python和CUDA)(一)——零基础安装pycuda GPU编程(基于Python和CUDA)(二)——显示GPU信息 零基础安装pycuda 系列文章目录前言安装Visual Studio安装C/C环…

【附源码】Python-3.9.5安装教程

软件下载 软件:Python版本:3.9.5语言:英文大小:26.9M安装环境:Win11/Win10/Win8/Win7硬件要求:CPU2.5GHz 内存2G(或更高)下载通道①百度网盘丨64位下载链接:https://pan.baidu.com/…

ConsoleApplication17_2项目免杀(Fiber+VEH Hook Load)

加载方式FiberVEH Hook Load Fiber是纤程免杀,VEH是异常报错,hook使用detours来hook VirtualAlloc和sleep,通过异常报错调用实现主动hook 纤程Fiber的概念:纤程是比线程的更小的一个运行单位。可以把一个线程拆分成多个纤程&#…

promethues监控postgres,emqx,redis

一、监控postgres 1、安装监控 docker pull wrouesnel/postgres_exporter2、执行 docker run -d -p 9187:9187 --name postgres_exporter --nethost -d -e DATA_SOURCE_NAME"postgresql://postgres:123456192.168.12.116:5432/rcc-manage?sslmodedisable" wroues…