MRI多任务技术及应用

news2024/11/26 14:58:03

目录

  • 一、定量心血管磁共振成像(CMR)的改进方法
  • 二、磁共振多任务
  • 三、磁共振多任务的成像框架
  • 四、磁共振多任务的图像模型和采样和重建策略
  • 五、利用MR多任务进行快速三维稳态CEST(ss-CEST)成像
    • 5.1 利用MR多任务进行快速三维稳态CEST(ss-CEST)成像介绍
    • 5.2 利用MR多任务进行快速三维稳态CEST(ss-CEST)成像
    • 5.3 利用MR多任务进行快速三维稳态CEST(ss-CEST)成像的讨论


一、定量心血管磁共振成像(CMR)的改进方法

  可以将生理运动和其他动态过程概念化为多个时间维度,通过 低秩张量(LRT)成像解决运动伪影,实现多达四个时间维度的运动解决定量成像。这种连续采集的方法,称之为CMR多任务。使用捕捉运动而非避免运动的方法不需要使用心电图触发或病人憋气,以有效地进行CMR定量成像。

  CMR在各种重叠的动态过程中成像,包括生理(如心脏跳动、呼吸运动)和物理(如T1、T2弛豫)动态,使得成像过程复杂化。传统处理心血管成像中的动态重叠的策略使用复杂的心电图(ECG)控制、病人屏气和/或短暂的采集“冻结”尽可能多的动态,但这意味着放弃了其余动态的有用的时间,而且在采集中的每个动态都需要采用不同动态所对应的冻结机制组合

  对于那些有心律失常或憋气困难的不健康受试者来说,CMR检查的时间是冗长的。


二、磁共振多任务

  磁共振多任务(multitasking)是一个连续采集框架,可以 同时解决定量心脑血管成像中涉及的许多重叠的动态。磁共振多任务将图像动态的不同来源概念化为不同的时间维度,利用多任务来捕捉(而非避免)运动、松弛和其他动态,进行时间分辨的T1映射,纠正了动态对比度增强(DCE)成像中的信号饱和问题

  LRT图像模型允许磁共振多任务,利用不同时间点的图像之间的相关性来减少采样要求(也可通过压缩、磁共振指纹减少采样),打破了“扫描时间随维度的增加而呈指数级增长”的维度诅咒。LRT的自由度及其所需的扫描时间与维度的数量呈近似线性关系,适合于多维成像。
在这里插入图片描述

图为非心电图、自由呼吸的原生心肌T1图的多个时间维度的说明和分析

  a.不同图像在三个时间维度空间的位置。不同的T,加权位于反转时间轴(水平)不同的心动期位于心动时间轴(深度)不同的呼吸期位于呼吸时间轴(垂直)
  b.描述图像张量每个维度的三个最重要的基函数,由1分钟的数据重建
  c.来自12.3分钟的高阶svD的奇异值曲线的原始子空间训练数据(即有足够的数据来覆盖所有的运动状态和对比度组合),表明对于原始数据张量的所有展开,奇异值都会迅速衰减

  磁共振多任务处理解决了定量CMR中许多长期的限制
  (1)有效地处理了运动,消除了对心电图控制和/或屏气的依赖,为心律失常患者进行定量CMR提供了潜在的手段。
  (2)纠正了从T1加权成像中量化造影剂浓度的不准确性,允许从单一造影剂中进行量化
  (3)简化了工作流程,使用单一的免设置的扫描来产生共同登记(?)的、运动分辨的参数图而不是一连串的错位扫描,每个扫描都涉及一个复杂的设置过程,以确定心电图触发延迟时间、呼吸门控窗口中心和宽带,和/或选择适当图像对比度的时间参数。
  该方法允许灵活的采样有效的因子张量重建。与磁共振指纹不同的是,除了核磁共振松弛外,它还能 对运动和DCE进行成像。此外,由于CMR多任务与LRT成像可以扫描广泛的“自然”成像对比(即与磁共振指纹的随机图像对比相反),多任务产生的图像也可能反过来验证它们的准确性。
  LRT成像不同于其他多为成像方法,如XD-GRASP通过隐含地假设沿每个时间维度横向片状恒定演化来利用图像的“局部”相似性,而LRT成像则全面利用了在整个多维时间空间中的图像的横向和对角线上都有的相关性


三、磁共振多任务的成像框架

  CMR多任务将一组心脑血管图像表示为一个多为张量(或阵列),其中一个维度索引体素位置(即结合空间维度),其他维度索引N个不同的时间维度每个维度对应于要成像的不同“任务”。通过对这个张量进行低秩的建模,我们描述并利用图像的相关性(沿着每个时间维度和整个多维时间维度的对角线)。
  通过将稀疏采样的图像数据与经常对K空间子集进行采样的辅助子空间训练数据交错获取数据,为图像重建开发了一种记忆和时间效率高的因子方法,其中核心张量和N个时间因子矩阵是从子空间训练数据中估计出来的,空间因子矩阵是通过将核心张量和时间因子矩阵与测量数据的其余部分进行拟合而恢复的。
在这里插入图片描述

CMR多任务技术用于非心电图、自由呼吸的原生心肌T1图

  a.CMR多任务技术产生了沿反转时间维度的精细分辨的对比度变化
  b.多任务技术产生的心脏分辨的T1图


四、磁共振多任务的图像模型和采样和重建策略

  CMR多任务心血管图像表示为空间位置的多维函数I(x,t1,t2,…tN)和N个时间维度t1,t2,…tN每个时间维度对应要解决的不同“任务”,例如时间维度对应于心脏运动、呼吸运动、磁化准备后的时间(沿途的信号演变取决于组织特性参数,如T1和T2)和成像开始后的时间(对描述通过的时间很有用)。
  尽管LRT模型释放了从维度诅咒中获得的采样要求,但N+1路张量(或称多维数组)A的大小仍然是指数级的增长。因此,通过单值分解(SVD)阈值化解决方案涉及的相同大小的辅助变量,以未压缩的形式存储图像I用离散形式表示为的N+1路张量(或称多维数组)A对内存的要求是非常高的。此外,图像重建算法的每一次迭代都可能涉及到对A的每一列(其数量可能达到数十万)的操作,以及大矩阵的多个SVD,所有这些都需要大量的计算时间。
  训练数据比K空间的其他位置需要更频繁的采样。这个数据子集(“子空间训练数据”)包含有限的空间信息,但有大量适合确定Φ的时间信息。收集的数据要足够频繁,以解决最精细的生理学问题。
  如果有序的均匀采样与时间基础函数一致,例如定期采样与呼吸周期、心动周期或磁化准备期同步,这时有序的均匀采样的采样性能会得到心脏和呼吸频率变化的加持。磁化准备时间表通常正好是周期性的,但准备期和采样期的选择可以避免相互同步。为了避免这些问题,可以采用黄金角径向取样或随机笛卡尔取样的采样方案,即使在周期性运动的情况下也能提供不连贯性。
在这里插入图片描述

CMR多任务处理非心电图,第一道心肌灌注T绘图

  a.对比剂动态被捕捉到的收缩期和舒张期心脏阶段
  b.这是由于该方法能够解决心脏运动
  c.经过的时间维度(用于描绘对比剂动态)和饱和时间维度相结合,产生二维信号强度表面,而不是传统的一维信号强度曲线
  d.这些信号强度表面被用来绘制R(t),它考虑了信号饱和度,并在线性转换后直接得到Gd浓度


五、利用MR多任务进行快速三维稳态CEST(ss-CEST)成像

5.1 利用MR多任务进行快速三维稳态CEST(ss-CEST)成像介绍

  化学交换饱和转移是―种非对比性的MRI技术,通过在不同的频率偏移下的预饱和,间接地检测出pool中的可交换元化学交换饱和转移核磁共振成像提供了一种新的对比机制来成像重要的生理学信息,如在交换质子池和水池之间的PH值和me-tabolite浓度。在不同的饱和频率下收集图像,产生所谓的Z谱,它反映了在给定饱和功率下采样频率偏移的稳态信号
  Z谱的广泛和对称覆盖允许多池分析,同时揭示不同的CEST效应(如酰胺质子转移(APT)、中继核Overhauser增强(rNOE)效应),以及其他特定的应用效应(如糖化CEST、CrCEST和糖化NOE)。为了实现可靠的多池分析,通常要对宽Z谱进行密集采样。考虑到通常采样的几十个频率偏移,每个频率偏移的采集时间(包括长的饱和模块)最好限制在几秒钟,以保持临床实践中可接受的扫描时间。这种时间限制通常只允许对饱和模块后的每个频率偏移进行单次K空间采集。这种单次拍摄方案使得快速、高质量的三维CEST成像变得困难
在这里插入图片描述
(A)序列设计:每个模块(TR=70 ms)包含一个单叶高斯饱和脉冲(txat = 30 ms,翻转角=500°),然后是一个扰流梯度和八个FLASH读出线(一个训练线加七个成像线;翻转角=5°)
该模块在一个特定的频率偏移(Nw=80)处重复,然后切换到另一个
(B)k-空间采样模式:在每个模块中﹐首先获得“训练线”(中心k空间线),然后是七条"成像线"(在xy和ga方向具有高斯分布的伪随机采样线)


5.2 利用MR多任务进行快速三维稳态CEST(ss-CEST)成像

  通过优化K空间采样效率,如使用螺旋中心重排序的K空间采集进行快照梯度回波(GRE)读出,或使用CAIPIRINHA下采样的三维EPI读出,缩短了扫描时间。它可以提供1.7×1.7毫米的2面内分辨率,FOV为220×180×54毫米,使用三维GRE读出,每次偏移7秒获得,使用三维EPI读出,每次偏移4.3秒获得1.8毫米的同向分辨率,FOV为256×224×156毫米
  一种潜在的更快的方法,即稳态CEST(ss-CEST)方法以重复模块的交错模式进行预饱和和k空间采样以确保饱和交换稳态在每个频率偏移的大部分时间内保持不变,交错模式在序列设计和可能的加速方面提供了更多的灵活性。
  但是最初的ss-CEST方法需要12分钟以上的时间来获取整个Z谱,这对于实际使用来说仍然太长。我们可将径向读出与多线奇异值分解相结合,进一步减少总的扫描时间,使其少于5分钟。与以前的ss-CEST方法相比,每个频率偏移的采集时间从10多秒减少到7.6秒。
  磁共振多任务是一种低秩张量的成像策略。通过低秩张量建模,利用在不同频率偏移下获得的图像之间的相关性,以及在接近稳态期间的数据之间的相关性,以减少扫描时间并提高图像质量。这样就可以在5.5分钟内获得覆盖全脑的1.7×1.7×3.0毫米的Z谱的空间分辨率。
在这里插入图片描述

二维单次FLASH法(A)和拟议的多任务ss-CEST法(B)所产生的map的比较

  注意二维单次FLASH CEST方法中使用的切片厚度为10毫米,大于拟议的三维多任务ss-CEST方法中的3毫米。因此,尽管(A)和(B)之间的切片中心相匹配,但空间覆盖率并不完全相同。


5.3 利用MR多任务进行快速三维稳态CEST(ss-CEST)成像的讨论

  与单次采集方法(或假稳态方法)相比,ss-CEST方法确保化学交换过程的稳态在大部分采集时间内得到保持。与单次拍摄方法相比﹐它允许连续采集,以便更有效地穿越K空间。但在原始ss-cEST实现中,仅在平行成像中拥有额外加速源
  而MR-Multitasking中的低秩张量模相对于两个单独的序列参数维度:频率偏移和频率增加后的时间,强制执行综合空间维度的低秩性(因此没有对空间结构进行假设)。此外,多任务系统在图像重建过程中采用这种模式来加速采集。
  除了加速采集,多作务ss-CEST方法还有两个额外的优点:
  (1)在每个频率偏移处的稳态方法被建模,稳态之前的信号被排除允许使用未被破坏的稳态值进行量化
  (2)Z谱在低秩约束下被自动去噪。考虑到足够的时空关系,加速的唯一明确限制因素为在每个频率偏移处达到稳定状态所的时间
  还可以进一步优化多作务ss-CEST方法。首先,我们可以利用先进的K-空间采样轨迹,例如非笛卡尔轨迹(如螺旋线),有可能比笛卡尔采集提高采样效率和不连贯性进一步减少扫描时间和提高图像质量
  其次,可以 进一步优化特定频率偏移的采样模式,以减少采样频率的总数,这样就可以减少总的扫描时间,同时保持多pool拟合数据的可靠性和稳健性。

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

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

嵌入式学习笔记(14)位置有关编码

位置无关编码(PIC,position independent code):汇编源文件被编码成二进制可执行程序时编码方式与位置(内存地址)无关。在我们写程序时,必须给链接器指定地址。将来的程序被执行时必须放在当时链接时给定的地…

iPhone照片导入电脑的绝佳教程,赶快get起来!

为什么都说苹果手机拍照比安卓手机好看?因为苹果手机拍照更加真实,其镜头和硬件能让照片更好地接近原始色彩。手机拍照好看就会让人忍不住多拍几张,久而久之手机内存就满了。有什么方法能将照片保存到电脑上吗?iPhone照片导入电脑…

《vue3实战》运用push()方法实现电影评价系统的添加功能

目录 前言 电影评价系统的添加功能是什么? 电影评价系统的添加功能有什么作用? 一、push()方法是什么?它有什么作用? 含义: 作用: 二、功能实现 这段是添加开始时点击按钮使…

私募证券基金动态-23年8月报

成交量:8月日均8,252.00亿元 8月A股两市日均成交8,252.00亿元,环比下降12.23%、同比下降18.11%。8月整体23个交易日,仅有3个交易日单日成交金额过万亿,且成交量起伏较大,单日成交金额最低仅有6805.32亿元(…

C# 实现电子签名

本项目基于Emgu.CV(C#下OpenCv的封装)开发的,编译器最新版Vs2022,编译环境x86 直接看效果图 1.主页面 2.我们先看手写的方式: 点击确认就到主界面,如下 : 点击自动适配-,再点击生成…

windows安装向量数据库milvus

本文介绍windows下安装milvus的方法。 一.Docker安装 1.1docker下载 首先到Docker官网上下载docker:Docker中文网 官网 1.2.安装前前期准备 先使用管理员权限打开windows powershell 然后在powershell里面输入下面那命令,启用“适用于 Linux 的 Windows 子系统”…

JMeter压力测试入门教程

Apache JMeter是Apache组织开发的基于Java的压力测试工具。用于对软件做压力测试,它最初被设计用于Web应用测试但后来扩展到其他测试领域。 它可以用于测试静态和动态资源例如静态文件、Java小服务程序、CGI脚本、Java 对象、数据库, FTP服务器, 等等。J…

JavaScript中的深拷贝和浅拷贝

聚沙成塔每天进步一点点 ⭐ 专栏简介⭐ 浅拷贝(Shallow Copy):⭐深拷贝(Deep Copy):⭐ 写在最后 ⭐ 专栏简介 前端入门之旅:探索Web开发的奇妙世界 记得点击上方或者右侧链接订阅本专栏哦 几何带…

数据结构day7栈-顺序栈的实现

用指针比用数组好&#xff0c;这样用户可以自己指定空间的大小&#xff0c;有参与感。 全部代码: main.c #include <stdio.h> #include <string.h> #include "sqstack.h"int main(int argc, char *argv[]) {sqstack *s;int i;s stack_create(100);if(…

数据结构 - 单链表

文章目录 目录 文章目录 一、什么是链表? 线性表: 顺序表: 二、链表的分类和实现 分类: 实现: 1.创建节点类 2.创建单链表 1.addTail(尾增) 2.删除节点值为key的第一个节点 3.插入节点(在指定位置) 4.获取链表长度 总结 前言 大家好,这篇博客给大家讲一下什么是…

《代码随想录》刷题笔记——数组篇【java实现】

*二分查找 题目链接 https://leetcode.cn/problems/binary-search/ 左闭右闭区间实现 时间复杂度&#xff1a;O(log n)空间复杂度&#xff1a;O(1) /*** 左闭右闭写法** param nums* param target* return*/ public static int search1(int[] nums, int target) {if (nums…

关于工信部发布的app备案以及小程序备案流程

一、相关政策 通知&#xff1a;https://beian.miit.gov.cn/#/Integrated/lawStatute 腾讯备案&#xff1a;网站备案 首次备案-网站备案-文档中心-腾讯云 阿里备案&#xff1a;网站备案_ICP备案_备案迁移_备案-阿里云 二、遇到的问题 APP备案 安卓获取平台公钥方法&#xf…

Vue + Element UI 前端篇(九):接口格式定义

接口请求格式定义 前台显示需要后台数据&#xff0c;我们这里先把前后端交互接口定义好&#xff0c;没有后台的时候&#xff0c;也方便用mock模拟。 接口定义遵循几个规范&#xff1a; 1. 接口按功能模块划分。 系统登录&#xff1a;登录相关接口 用户管理&#xff1a;用户…

Redis-1.4-过期策略

1 设置带过期时间的 key # 时间复杂度&#xff1a;O&#xff08;1&#xff09;&#xff0c;最常用方式 expire key seconds# 字符串独有方式 setex(String key, int seconds, String value)除了string独有设置过期时间的方法&#xff0c;其他类型都需依靠expire方法设置时间&a…

说说MySQL回表查询与覆盖索引

分析&回答 什么是回表查询&#xff1f; 通俗的讲就是&#xff0c;如果索引的列在 select 所需获得的列中&#xff08;因为在 mysql 中索引是根据索引列的值进行排序的&#xff0c;所以索引节点中存在该列中的部分值&#xff09;或者根据一次索引查询就能获得记录就不需要…

【Linux】Qt Remote之Remote开发环境搭建填坑小记

总体思路 基于WSL2&#xff08;Ubuntu 22.04 LTS&#xff09;原子Alpha开发板进行Qt开发实验&#xff0c;基于Win11通过vscode remote到WSL2&#xff0c;再基于WSL2通过Qt 交叉编译&#xff0c;并通过sshrsync远程到开发板&#xff0c;构建起开发工具链。 Step1 基于Win11通过…

Mock接口测试

什么是mock? 测试桩&#xff0c;模拟被测对象的返回&#xff0c;用于测试 通常意义的mock指的就是mock server, 模拟服务端返回的接口数据&#xff0c;用于前端开发&#xff0c;第三方接口联调 为什么要mock? 1. 解决依赖问题&#xff1a;当我们测试一个接口或者功能模块…

Python Opencv实践 - 轮廓特征(最小外接圆,椭圆拟合)

import cv2 as cv import numpy as np import matplotlib.pyplot as pltimg cv.imread("../SampleImages/stars.PNG") plt.imshow(img[:,:,::-1])#轮廓检测 img_gray cv.cvtColor(img, cv.COLOR_BGR2GRAY) ret,thresh cv.threshold(img_gray, 127, 255, 0) contou…

Unity汉化一个插件 制作插件汉化工具

我是编程一个菜鸟&#xff0c;英语又不好&#xff0c;有的插件非常牛&#xff01;我想学一学&#xff0c;页面全是英文&#xff0c;完全不知所措&#xff0c;我该怎么办啊...尝试在Unity中汉化一个插件 效果&#xff1a; 思路&#xff1a; 如何在Unity中把一个自己喜欢的插件…

Flutter 最优秀动画库「完全商业化」,Rive 2 你全面了解过吗?

说到 rive &#xff0c;非 Flutter 开发者可能会感觉比较陌生&#xff0c;而做过 Flutter 开发的可能对 rive 会有所耳闻&#xff0c;因为 rive 在一开始叫 flare &#xff0c;是 2dimensions 公司的开源动画产品&#xff0c;在发布之初由于和 Flutter 团队有深入合作&#xff…