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

目录

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


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

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

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

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

相位差重建步骤流程

二、反正切运算

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

相位卷绕

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

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

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


四、图形变形校正

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

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

如若内容造成侵权/违法违规/事实不符,请联系我们进行投诉反馈qq邮箱809451989@qq.com,一经查实,立即删除!

相关文章

Mysql高级语句

高级语句 1.按关键字排序 SELECT column1, column2, ... FROM table_name ORDER BY column1, column2, ... ASC|DESC ASC 是按照升序进行排序的,是默认的排序方式,即 ASC 可以省略。 SELECT 语句中如果没有指定具体的排序方式,则默认按 ASC…

【数据结构】——查找、散列表的相关习题

目录 一、选择填空判断题题型一(顺序、二分查找的概念)题型二(分块查找的概念)题型三(关键字比较次数) 二、应用题题型一(二分查找判定树) 一、选择填空判断题 题型一(顺…

Leetcode328 奇偶链表

思路:分别处理奇偶,保存奇偶的第一个和最后一个节点,注意最后链接的时候需要把偶数的next去掉再拼接不然就成环了 class Solution:def oddEvenList(self, head: ListNode) -> ListNode:if not head or not head.next or not head.next.ne…

YOLOv5算法改进(10)— 替换主干网络之GhostNet

前言:Hello大家好,我是小哥谈。GhostNet是一种针对计算机视觉任务的深度神经网络架构,它于2020年由中国科学院大学的研究人员提出。GhostNet的设计目标是在保持高精度的同时,减少模型的计算和存储成本。GhostNet通过引入Ghost模块…

阿里云架构

负载均衡slb 分类以及应用场景 负载均衡slb clb 传统的负载均衡(原slb) 支持4层和7层(仅支持对uri(location),域名进行转发) 一般使用slb(clb) alb 应用负载均衡 只支持7层,整合了nginx负载均衡的各种功能,可以根据用户请求头,响应头 如果需要详细处理用户请求(浏…

【Dart】学习使用(二):基本类型

前言 基本类型是语言的基础。 Dart 语言支持以下基础类型:Numbers(int、double), 整形Strings(String), 字符串Booleans(bool) , 布尔型Records((value1,value2)) 记录Lists(List ) 数组Sets(Set) 集合Maps(Map) 映射Runes(Runes,通常由 characters AP…

springBoot--终

目录 Web定制SpringMVC定制某些mvc配置定制mvc核心组件完全自定义mvc 静态资源静态资源映射静态资源缓存示例自定义静态资源规则配置方式代码方式 欢迎页Favicon(网站图标)路径匹配策略内容协商制定返回json类型数据制定返回xml类型数据制定返回自定义类…

2023年高教社杯 国赛数学建模思路 - 案例:ID3-决策树分类算法

文章目录 0 赛题思路1 算法介绍2 FP树表示法3 构建FP树4 实现代码 建模资料 0 赛题思路 (赛题出来以后第一时间在CSDN分享) https://blog.csdn.net/dc_sinor?typeblog 1 算法介绍 FP-Tree算法全称是FrequentPattern Tree算法,就是频繁模…

本地部署体验LISA模型(LISA≈图像分割基础模型SAM+多模态大语言模型LLaVA)

GitHub地址:https://github.com/dvlab-research/LISA 该项目论文paper reading:https://blog.csdn.net/Transfattyacids/article/details/132254770 在GitHub上下载源文件,进入下载的文件夹,打开该地址下的命令控制台,…

Vue安装过程的困惑解答——nodejs和vue关系、vue的项目结构

文章目录 一、为什么在使用vue前要下载nodejs?二、为什么安装nodejs后就能使用NPM包管理工具?三、为什么是V8引擎并且使用C实现?四、为什么会安装淘宝镜像?五、什么是webpack模板,为什么需要他?六、vue项目…

深入探索C语言自定义类型:打造你的编程世界

一、什么是自定义类型 C语言提供了丰富的内置类型,常见的有int, char, float, double, 以及各种指针。 除此之外,我们还能自己创建一些类型,这些类型称为自定义类型,如数组,结构体,枚举类型和联合体类型。 …

Unreal5(虚幻5)学习记录 快捷键

虚幻5学习记录 快捷键 世界场景中漫游(镜头移动): 按住鼠标右键 键盘的W(前) S(后) A(左) D(右) E(上) Q(下)键 透视 透视 ALTG 上部分 ALTJ 底视图ALTSHIFTJ 左视图 ALTK 右视图 ALTSHIFTK 前视图 ALTH 后视图 ALTSHIFTH 内容浏览器 Ctrl Space 内容浏览器…

stm32CubeMX HAL W5500芯片介绍 第一章

W5500芯片介绍 文章目录 W5500芯片介绍简单简绍以太网以太网分五层:第一层物理层:第二层:数据链路层:第三层:网络层:第四层:传输层:第五层:应用层:以太网应用…

Scala的特质trait与java的interface接口的区别,以及Scala特质的自身类型和依赖注入

1. Scala的特质trait与java接口的区别 Scala中的特质(trait)和Java中的接口(interface)在概念和使用上有一些区别: 默认实现:在Java中,接口只能定义方法的签名,而没有默认实现。而在…

Matlab(结构化程式和自定义函数)

目录 1.脚本编辑器 2.脚本流 2.1 控制流 2.2 关系(逻辑)操作符 3.脚本与函数 1.脚本编辑器 Matlab的命名规则: 常用功能: 智能缩进: 在写代码的时候,有的时候代码看起来并不是那么美观(可读性…

单片机通用学习-​什么是寄存器?​

什么是寄存器? 寄存器是一种特殊的存储器,主要用于存储和检查微机的状态。CPU寄存器用于存储和检查CPU的状态,具体包括计算中途数据、程序因中断或子程序分支时的返回地址、计算结果为零时的负值、计算结果为零时的信息、进位值等。 由于CP…

怎么提取视频中的音乐保存到本地?其实方法很简单

当你想要使用视频中的音乐时,你可以考虑将它从视频中提取出来。这可以用于制作音频样本集,制作铃声或其他音频素材,或者向其他人展示视频的音乐部分而无需显示视频本身。如果你是一位音乐制作人员,你可能会需要一些特定类型的音效…

HTML之VSCode简单配置与创建

目录 插件下载 然后输入源码&#xff1a; 使用 效果 插件下载 下载这个插件后可以直接运行&#xff1a; 然后创建一个文件&#xff1a; 然后输入源码&#xff1a; <!DOCTYPE html> <html lang"en"> <head><meta charset"UTF-8"…

Vue3列表竖向滚动(包含使用swiper的翻页效果)

一、使用element-plus表格进行滚动&#xff1a; 可以满足的需求&#xff1a;表格一行一行竖向滚动&#xff0c;类似走马灯。 不能满足的需求&#xff1a;表格分页竖向滚动&#xff0c;有翻页的效果。 代码&#xff1a; <template><el-table:data"tableData"…

AMEYA360:兆易创新获得ISO 26262 ASIL D流程认证, 汽车功能安全管理体系再上新台阶

中国北京(2023年8月29日) —— 业界半导体器件供应商兆易创新GigaDevice(股票代码 603986)今日宣布&#xff0c;获得由国际公认的测试、检验和认证机构通标标准技术服务有限公司(以下简称SGS)授予的ISO 26262:2018汽车功能安全最高等级ASIL D流程认证证书&#xff0c;这标志着兆…