书签 分享 收藏 举报 版权申诉 / 18
立即下载 加入VIP,免费下载
当前位置： 首页 > 物理 > 测量；测试 > 磁共振成像的方法和系统及其用途.pdf

磁共振成像的方法和系统及其用途.pdf

  • 上传人：1***
  • 文档编号：1696079
  • 上传时间：2018-07-05
  • 格式：PDF
  • 页数：18
  • 大小：523.90KB
    • 摘要
    申请专利号：

    CN201080021664.6

    		 申请日：

    2010.04.22
    公开号：

    CN102428383A

    		 公开日：

    2012.04.25
    当前法律状态：

    授权

    		 有效性：

    有权
    法律详情：

    授权|||实质审查的生效IPC(主分类):G01R 33/563申请日:20100422|||公开
    IPC分类号：

    G01R33/563; A61B5/055; G01N24/08; G01R33/46

    		 主分类号：

    G01R33/563
    申请人：

    CR发展公司
    发明人：

    丹尼尔·托普加德
    地址：

    瑞典伦德
    优先权：

    2009.05.22 SE 0950363-2; 2009.05.22 US 61/180,752
    专利代理机构：

    北京安信方达知识产权代理有限公司 11262

    		 代理人：

    周靖;郑霞
    PDF完整版下载： PDF下载
    内容摘要

    本发明涉及用于磁共振成像或者核磁共振光谱的方法，包括：向受到磁场影响的物体发射射频和梯度脉冲序列，其中所述物体包括具有核自旋非0的原子的分子；编码、检测以及获取来自所述物体的对应于所述发射的射频和梯度脉冲序列的磁共振信号，其中所述射频和梯度脉冲序列包括第一加权块、具有持续时间tm的混合块和第二加权块；通过增加所述第一加权块和所述第二加权块中的至少一个的强度，将编码、检测以及获取来自所述物体的磁共振信号限制到信号强度I的初始衰减，其中分析所述初始信号衰减速率随tm的变化来得到表观交换速率AXR。

    权利要求书

    1： 一种用于磁共振成像或者核磁共振光谱的方法， 包括 ： 向受到磁场影响的物体发射射频和梯度脉冲序列， 其中所述物体包括具有核自旋非 0 的原子的分子 ； 编码、 检测以及获取来自所述物体的对应于所述发射的射频和梯度脉冲序列的磁共振 信号， 其中所述射频和梯度脉冲序列包括第一加权块、 具有持续时间 tm 的混合块和第二加 权块 ； 其特征在于， 通过增加所述第一加权块和所述第二加权块中的至少一个的强度， 将编码、 检测以及 获取来自所述物体的磁共振信号限制到所述信号的强度 I 的初始衰减， 其中分析所述初始 信号衰减速率随 tm 的变化来得到表观交换速率 AXR。
    2： 如权利要求 1 所述的方法， 其中 所述第一加权块是具有弛豫加权 τ1 的第一弛豫加权块， 并且所述第二加权块是具有 弛豫加权 τ2 的第二弛豫块， 并且其中通过增加 τ1 或者 τ2 将编码、 检测以及获取来自所 述物体的磁共振信号限制到所述信号的强度 I 的初始衰减来得到表观弛豫速率 R ； 或者其中 所述第一加权块是具有扩散加权 b1 的第一扩散加权块， 并且所述第二加权块是具有扩 散加权 b2 的第二扩散块， 并且其中通过增加 b1 或 b2 将编码、 检测以及获取来自所述物体的 磁共振信号限制到所述信号的强度 I 的初始衰减来得到表观扩散系数 ADC ； 并且其中， 其后分析所述表观弛豫速率 R 或者所述表观扩散系数 ADC 随 tm 的变化来得 到所述表观交换速率 AXR。
    3： 如权利要求 1 或 2 所述的方法， 其中 AXR 用于产生图像对比度。
    4： 如权利要求 1-3 任一项所述的方法， 其中所述第一加权块是具有扩散加权 b1 的第一 扩散加权块， 并且所述第二加权块是具有扩散加权 b2 的第二扩散块， 并且其中根据下式限 制所述磁共振信号的编码、 检测以及获取 ： 其中， 变量 b 指的是 b1 或 b2， 并且 是表观扩散系数。
    5： 如权利要求 4 所述的方法， 其中根据下式 随 tm， n 的变化依赖于交换速率 k ： 其中 是在平衡处的 ， 并且 是在 tm， n ＝ 0 处的 。
    6： 如权利要求 1-5 任一项所述的方法， 其中进行所述编码、 检测以及获取以用于多位 置交换。
    7： 如权利要求 6 所述的方法， 其中根据下式计算用于多位置交换的所述表观交换速率 (AXR) ： In(b) ＝ I0， nexp{-ADC·[1-σexp(-AXR·tm， n)]·b} 其中 ADC ＝ ， 并且 2 其中 σ 限制在范围 0 ＜ σ ＜ 1 内。
    8： 如上述任一项权利要求所述的方法， 其中所述物体是沸石、 脂质体、 泡囊或者生物细 胞。
    9： 如上述任一项权利要求所述的方法， 其中所述第一加权块是具有扩散加权 b1 的第一 扩散加权块， 并且所述第二加权块是具有扩散加权 b2 的第二扩散块， 并且其中下述中的一 个是有效的 ： -b1 处于开或关模式， b2 固定且 tm 固定 ； -b2 处于开或关模式， b1 固定且 tm 固定 ； -b1 处于开或关模式， b2 固定且 tm 变化 ； -b2 处于开或关模式， b1 固定且 tm 变化 ； -b1 固定， b2 变化且 tm 变化， 其由 b1 等于零、 tm 固定且 b2 变化的系列补充 ； 或者 -b2 固定， b1 变化且 tm 变化， 其由 b2 等于零、 tm 固定且 b1 变化的系列补充。
    10： 如上述任一项权利要求所述的方法， 其中使用非对称的脉冲对或梯度调制来提高 滤波器效率。
    11： 如上述任一项权利要求所述的方法， 其中进行数值优化来提高滤波器效率。
    12： 一种用于磁共振成像或者核磁共振光谱的系统， 包括 ： 射频和梯度脉冲单元， 其用于向受到磁场影响的物体发射射频和梯度脉冲序列， 其中 所述物体包括具有核自旋非 0 的原子的分子 ； 检测器单元， 其用于检测来自所述物体的对应于所述发射的射频和梯度脉冲序列的磁 共振信号， 其中所述射频和梯度脉冲序列包括第一加权块、 具有持续时间 tm 的混合块和第 二加权块 ； 其特征在于， 所述检测器单元被设置为用于检测， 通过增加所述第一加权块和所述第二加权块中的 至少一个的强度而将所述检测限制到所述信号的强度 I 的初始衰减， 以便可以分析所述初 始信号衰减速率随 tm 的变化来得到表观交换速率 AXR。
    13： 如权利要求 12 所述的系统， 其中 所述第一加权块是具有弛豫加权 τ1 的第一弛豫加权块， 并且所述第二加权块是具有 弛豫加权 τ2 的第二弛豫块， 并且其中通过增加 τ1 或者 τ2 将编码、 检测以及获取来自所 述物体的磁共振信号限制到所述信号的强度 I 的初始衰减来得到表观弛豫速率 R ； 或者其中 所述第一加权块是具有扩散加权 b1 的第一扩散加权块， 并且所述第二加权块是具有扩 散加权 b2 的第二扩散块， 并且其中通过增加 b1 或 b2 将编码、 检测以及获取来自所述物体的 磁共振信号限制到所述信号的强度 I 的初始衰减来得到表观扩散系数 ADC ； 以便其后可以分析所述表观弛豫速率 R 或者所述表观扩散系数 ADC 随 tm 的变化来得 到所述表观交换速率 AXR。
    14： 如权利要求 12 或 13 所述的系统， 还包括图像生成单元， 该图像生成单元被配置成 基于所述磁共振信号生成图像。
    15： 一种医疗工作站， 包括用于执行根据所述权利要求 1-11 中的任一项的方法的装 置。 3
    16： 如权利要求 1-11 中的任一项的方法、 如权利要求 12-14 中的任一项的系统、 或者如 权利要求 15 的医疗工作站的作为诊断疾病或疾患的诊断工具的用途。
    17： 如权利要求 1-11 中的任一项的方法、 如权利要求 12-14 中的任一项的系统、 或者如 权利要求 15 的医疗工作站的用于研究体内活细胞的新陈代谢的用途。
    18： 如权利要求 1-11 中的任一项的方法、 如权利要求 12-14 中的任一项所述的系统、 或 者如权利要求 15 的医疗工作站的用于研究医疗药物通过细胞膜的跨膜扩散的用途。

    说明书


    磁共振成像的方法和系统及其用途

        发明领域 本发明大体上涉及核磁共振和磁共振成像领域。
         发明背景
         扩散核磁共振 (NMR) 已经使用了 40 多年来确定自扩散系数， 其可以从集料粒度、 分子通过其移动的介质的渗透性、 以及发生在扩散物质和较大分子或多孔基体之间的结合 事件 (binding events) 的角度来解释。最常见的扩散 NMR 技术依赖于扩散编码块， 其包 括用于标记 NMR 无线频率信号在脉冲之间的时间发生的位移的成对的磁场梯度脉冲。扩 散 NMR 技术和分析方法不仅适用于体外而且还适用于病理状态的检测的医疗磁共振成像 (MRI) 的背景， 该病理状态例如缺血性卒中、 脱髓鞘作用紊乱、 以及肿瘤。在一些情况下， 尤 其是卒中情况， 基于扩散的图像对比度比基于核弛豫率 R1 和 R2 的对比度的更加传统的方 式提供更多信息。
         在具有不同的弛豫 / 扩散特性的部分之间的水交换速率是用于组织内的病理
         状态的潜在有用的标记。作为扩散的功能执行的扩散 NMR 实验可以利用 进行分析来估计交换速率模型J.， H.Pfeifer， 和 W.Heink.1988.Principles andapplications of self-diffusion measurements by nuclear magnetic resonance.Adv. Magn.Reson.12 ： 1-89)。该分析受到 NMR 信号对交换速率的弱依赖性的限制。 出 版 物 PT.Callaghan 和 I.Furó ， Diffusion-diffusion correlation and exchange as a signature for local order and dynamics.J.Chem. Phys.120(2004)4032-4038 以及 K.E.Washburn， 和 PT.Callaghan， Tracking pore to pore exchange using relaxation exchange spectroscopy.Phys.Rev.Lett.97(2006)175502. 公开了用于检查交换过程的二维核磁共振实验。该方法——被知晓为扩散交换光谱和弛 豫交换光谱——利用被混合时间分离的两个独立增加的弛豫 / 扩散加权块、 以及使用二维 拉普拉斯逆变换的信号分析。但是， 严重的缺点是为获得二维拉普拉斯逆变换分析所需要 的大量数据而对仪器时间的过度要求， 因此使该方法对于耐心有限的人类主体的研究不实 用。
         总的来说， 直到现在， 目前用于估计交换速率的可用的扩散 NMR 方法非常消耗时 间 (Callaghan) 或者依赖仅在评估参数和实验数据中的信息之间的弱依赖性的曲线拟合
         。其他已知的能够用于获得交换时间的方法是有创性方法， 例如通过光散射、 显 微镜检查、 吸收光谱和 X 射线的方式观察标记分子的扩散。这不仅仅由于毒理学风险难以 在体内使用而且人们永远无法保证组织和体液不会受到引入的标记的影响。 但是， PCT 申请 WO2008/147326 公开了一种方法， 其解决了上面揭露的问题。根据 PCT 申请 WO2008/147326 的方法包括向受到磁场影响的物体发射无线频率脉冲序列 ( 其中 该物体包括具有不为 0 的核自旋的原子的分子 )， 向该物体发射梯度脉冲序列， 检测该物体 的对应于发射的无线频率脉冲序列的磁共振信号， 并且处理该磁共振信号。
         进一步的， 根据 PCT 申请 WO2008/147326 的方法的特征为该梯度脉冲序列包括配 置成降低该物体的扩散系数的至少一部分分布的滤波块 ( 第一扩散加权块 )， 和在发射该
         滤波块后的预定时间发生的扩散编码块 ( 第二扩散加权块 )， 并且该过程包括将磁共振信 号的一部分与预定的磁共振信号的一部分比较， 产生比较信号， 其中该预定的磁共振信号 的一部分或者由用户定义或者由之前应用的梯度脉冲序列产生。
         根据 PCT 申请 WO2008/147326 的方法具有几个优点。 根据 PCT 申请 WO2008/147326 的发明的一般解决方案是它使用梯度脉冲的序列作为扩散实验的滤波器。 因而相同的分子 可以基于它们的扩散是如何被限制的而被单独地分析并且区分。 在各个部分之间的交换速 率是重要的参数， 其可以根据一些具体实施方式利用该发明得到。进一步的， 根据 PCT 申请 WO2008/147326 的发明提供用于材料 ( 例如组织 ) 的 MRI 研究的新的对比模式， 其中交换速 率作为位置的函数而变化。
         另外， 根据 PCT 申请 WO2008/147326 的发明具有的其他优点是所需要的总的实 验持续时间的缩短， 就这一点而言其使该发明能够用于体内， 例如， 诸如用于磁共振成像 (MRI) 的对比度的装置， 在一些情况下， 为了获得可依赖的交换速率结果， 不存在从其他实 验获取背景信息 ( 例如， 所研究的分子的形状或者扩散系数 ) 的需要， 以及提供其对比度依 赖于交换速率的不同的图像的可能性。
         但 是， 在 PCT 申 请 WO2008/147326 公 开 的 方 法 和 协 议 中 以 及 卡 拉 汉 协 议 (Callaghan’ s protocol) 中也存在问题。 一个问题是这些协议不能够在所有的 MRI 仪器上使用的事实。目前使用的一些标 准的 MRI 扫描器不能够获取足够的数据来允许例如根据 WO2008/147326 的方法的应用， 例 如用于整体双部分拟合 (global two component fit) 或 ILT 分析。
         临床 MRI 的另一个问题总的来说是高噪声水平。根据 WO2008/147326 由于大量的 数据获取， 事实上噪声可能是大问题。
         本发明的一个目的是提供用于 MRI 的方法， 该方法能够应用于非常大的范围的 MRI 扫描器， 例如目前使用的标准的临床全身 MRI 扫描器， 本发明的另一个目的是提供用于 MRI 的方法， 其受噪声影响不大。
         发明概述
         以上公开的目的通过一种用于磁共振成像或者核磁共振光谱的方法实现， 该方法 包括向受到磁场影响的物体发射射频和梯度脉冲序列， 其中所述物体包括具有核自旋非 0 的原子的分子 ；
         编码、 检测以及获取来自所述物体的对应于所述发射的射频和梯度脉冲序列的磁 共振信号， 其中所述射频和梯度脉冲序列包括第一加权块、 具有持续时间 tm 的混合块和第 二加权块 ； 其中
         通过增加所述第一加权块和所述第二加权块中的至少一个的强度， 将编码、 检测 以及获取来自所述物体的磁共振信号限制到信号的强度 I 的初始衰减， 其中分析所述初始 信号衰减速率随 tm 的变化来得到表观交换速率 AXR。
         与在 WO2008/147326 中说明的之前的协议相比， 需要较小的信号衰减来获得关于 交换的信息。对于扩散加权的情况， 较小的信号衰减意味着相应地可以使用扩散加权参数 b 和梯度幅值 G 的较小的值。该协议因此可以在任何临床 MRI 扫描器而不仅在具有非常高 的最大梯度强度 ( ＞ 40mT/m) 的最现代的 MRI 扫描器上实施。较小的信号衰减还意味着该 信号受噪声的影响较小。
         附图的简要说明
         图 1 示出用于扩散加权的核磁共振 (NMR) 或者磁共振成像 (MRI) 的众所周知的通 常的脉冲序列。
         图 2 示出了用于编码 NMR 或 MRI 的脉冲序列的示意图， 该 NMR 或 MRI 用于在具有 快扩散和慢扩散的部分之间的分子交换。
         图 3 示出了利用根据本发明的方法在不同溶液和细胞内执行的实验数据。
         本发明的一般基础
         本发明的一般基础在 PCT 申请 WO2008/147326 中公开， 该申请的全部通过引用被 并入。
         以下， 将简要说明该一般基础。
         扩散加权的 NMR 和 MRI
         图 1 示出了用于扩散加权的核磁共振 (NMR) 或者磁共振成像 (MRI) 的众所周知的 通常的脉冲序列。在检测块之前插入扩散加权 (DW) 块， 因而导致在 NMR 的每一个频道或者 在 MRI 的每一个像素的强度依赖分子运动的衰减。一般来说， DW 块由一对梯度脉冲 ( 最常 见的是矩形、 梯形、 或者正弦形状 ) 组成， 通过回聚 180°射频 (RF) 脉冲分离， 该 180°射频 (RF) 脉冲反转由之前应用的梯度脉冲导致的相位移动。第一和第二梯度脉冲分别表示散 焦脉冲和回聚脉冲。在图 1 的展开图中示出最简单有效的梯度形状。用于 NMR 的检测块通 常包括记录在自由旋进过程中的时域信号， 其在傅里叶变换 (FT) 以后产生 NMR 光谱。对于 MRI 应用， 检测块包括单个或者系列梯度或者 RF 回波 (echo)， 其基于 FT 产生 1D、 2D、 或者 3D 图像。通常的检测方案包括但不限于平面回波成像、 快速自旋回波、 螺旋成像、 螺旋桨成 像等。 DW 块的信号衰减通常使用扩散感应变量 (diffusion sensitizing variable)b 和 自扩散系数 D 通过式 (1) 报告， 其中 -bD
         I ＝ I0e (1)
         其中， I 是检测的信号强度并且 I0 是在零梯度力度处的信号强度。b 的值通过式 (2) 计算 ：
         其对于矩形梯度脉冲的估值为
         b ＝ (γGδ)2(Δ-δ/3)， (3)
         其中 γ 是磁旋比， G 是梯度脉冲的幅值， δ 是梯度脉冲的持续时间， 并且 Δ 是梯 度脉冲的前沿之间的间隔。
         D 通过下式与均方位移 相关
         ＝ 2Dt (4)
         其中 t 是测量扩散的时间间隔。
         在短梯度脉冲的限制中， 其中与在脉冲和系统的结构长度尺寸之间的时间发生的 位移比较， 在脉冲期间的分子位移是不显著的， 可以根据由 DW 块导致的信号衰减使用下式 2 估计 ，
         其中， q 是由散焦梯度脉冲导致的磁化螺旋的波矢量。q 的值通过下式由散焦梯度 脉冲的面积提供
         其对于矩形梯度脉冲等同于对于具有高斯扩散的系统， 如果估计的 指的是在有效扩散时间 td 内发生的位 移， 式 5 仍然是有效的， 对于矩形梯度脉冲 td 由下式提供 td ＝ Δ-δ/3。 (8) 2
         甚至对于具有非高斯扩散的系统， 可以根据初始值、 低 q、 在上面定义的短梯 度脉冲条件下的信号衰减估计。 对于梯度脉冲具有有限长度的情况， 人们可以根据初始值、 2 低 G、 E ＝ I/I0 的衰减定义表观均方位移 和对应的表观扩散系数 D(δ， Δ) ：
         和在球形细胞内的扩散
         对于限制在具有半径 r 的球形腔内的具有较大的扩散系数 D0 的流体， 可以根据 WO2008/147326 示出如下
         (11)其中， am 是下式的根通过对在式 (11) 中的指数进行一系列展开， 得到下面的极限操作 ：
         8 ＝ 2D0Δ(15)102428383 A CN 102428390
         说明书5/12 页受限情况和非受限情况在短 td 时一致。对于受限情况， 在长 td 和短 δ 时达到上 限。当 DW 块作为滤波器使用来消除来自非受限部分的信号时， 应当选择 δ 和 Δ 来最大化 2 1/2 在自由部分和受限部分之间的 的差， 同时保持 δ 和 Δ 比用于在部分之间 交换的特征时间要短得多。上面的式子使 DW 滤波器的合理设计成为可能。
         在细胞外和细胞内的部分之间的分子交换
         在细胞内和细胞外的部分内的分子数量 ni 和 ne 之间的比是均衡常数 K ：
         在细胞内和细胞外的流体内的自扩散以扩散系数 Di， 0 和 De 发生。为了方便起见， 在 De 内包括由球形细胞的出现产生的细胞外流体的阻力效果。根据式 (10) 和 (11)， 在细 胞内部分的有效扩散系数 Di 依赖于 δ 和 Δ 的值。在细胞膜之间的分子交换以渗透系数 P 发生。通过下式提供向外和向内的交换速率 ki 和 ke
         ki ＝ 3P/r (19)
         和
         交换速率与在细胞内和细胞外阶段的平均滞留时间 τi 和 τe 相关， τi 和 τe 经 由下式得出。
         τi ＝ 1/ki τe ＝ 1/ke。(21)
         使用宏观方法， 在 DW 块的信号衰减可以通过求解以下微分方程计算
         对于在细胞内和细胞外阶段的磁化强度 Mi 和 Me， 使用初始条件其中， M0 是在均分状况的总的磁化强度。式 (22) 假设在两个阶段的高斯扩散具有扩散 系数 Di 和 De。在两个阶段之间的交换以速率常数 ki 和 ke 发生。在式 (22) 内， 没有涉及系 统的微观几何结构。在短的梯度脉冲近似内， 式 (22) 的细胞内和细胞外的磁化强度 Mi， 1 和 Me， 1 的解在时间 t1 处是
         其中总的 NMR 信号与 Mi 和 Me 的和成比例， 并且比例 I/I0 等于 (Mi， 1+Me， 1)/(Mi， 0+Me， 0)。
         执行式 (24) 的以下修改来使它对于在球形细胞和连续的介质之间的交换及具有 有限长度的梯度脉冲也有效 ： 利用式 (10) 和 (11) 计算 Di 并且利用 td 替换指数中的 Δ。 只 要 δ ＜＜ τi， τe， 这个修改被预期是准确的。推导的表达式可以例如准确地描述在大范 围的 δ 和 Δ 取值内在酵母细胞的包围的沉积物内水的扩散加权。使式 (24) 适合于实验 数据产生了对用于交换的时标的估计， 但是由于实验观察到的信号强度对交换时间的依赖 性非常差， 因此这个方法对临床应用不实用。
         DW 滤波器 ： 消除细胞外的信号
         使用上面公开的模型和对 Di， De、 r、 K 和 P 的取值的合理假设， 对 {δ， Δ， G} 参 0、 数集做出明智的选择是可能的， 该参数集被设计来消除由细胞外部分产生的信号， 同时保 持来自细胞内部分的信号。更有效地， 可以使用迭代数值程序找到对于细胞内部分的给定 衰减最小化来自细胞外部分的信号的 {δ， Δ， G} 集合。
         当条件 δ ＜＜ τi， τe 没有被满足或者如果梯度调制比矩形调制更加复杂， 可以 根据 WO2008/147326 使用数值方法来求解式 (22)。
         扩散 - 扩散交换
         在图 2 中展示了用于编码 NMR 或 MRI 的脉冲序列的示意图， 该 NMR 或 MRI 用于在 具有快扩散和慢扩散的部分之间的分子交换。两个扩散加权块 DW1 和 DW2 位于信号检测之 前被具有持续时间 tm 的混合块分离。在图 2 中的每一个 DW 块遵循相同的限制， 并且可以 以与图 1 所示的 DW 块类似的方式分析。时间点 t0、 t1、 t2 和 t3 在图中标明。对于这些时间 点 tn 的每一个， 人们可以分别估计细胞内的磁化强度、 细胞外的磁化强度和总的磁化强度 的幅值 Mi，n、 Me，n 以及 Mn ＝ Mi，n+Me，n。在 t0 和 t1 之间的时间内磁化强度的改变可以通过式 (24) 提供。
         忽略核弛豫过程， 混合块的效果是重新分布在细胞内和细胞外部分之间的磁化强 度， 同时保持总的磁化强度的整体幅值， 即 M2 ＝ M1。在混合时间以后对总的磁化强度的相 对作用可以示出为
         在 DW2 块， 即 t2 和 t3 之间的时间的磁化强度 Mi 和 Me 内的演进再次由式 (24) 提供。 因此， 可以应用式 (24)、 然后应用式 (26)、 以及最后再次应用式 (24) 计算 Mi， 3 和 Me， 3。对于 新协议的合理设计， 用于 DW1 混合 DW2 序列以后的信号的理论表达式可能是关键的。
         推广
         为了简便， 上面的分析集中于在两个部分之间的交换 ： 一个部分自由扩散和一个
         部分以渗透膜受限在球形细胞内扩散。本领域技术人员应当理解， 上面的分析以及下面呈 现的协议可以推广到弛豫而不是扩散、 以及其他几何结构、 多个部分、 和梯度调制方案。
         卡拉汉协议
         卡拉汉之前介绍过与在图 2 中示出的相同的一般的脉冲序列 ( 卡拉汉， Furó. J.Chem.Phys.2004， 120， 4032)。以下面的方式执行根据卡拉汉的方法 ：
         1) 保持 δ1 ＝ δ2， Δ1 ＝ Δ2， 并且 tm 为常数， 独立地改变 G1 和 G2( 通常通过 16×16 ＝ 256 个分开的步骤 (separate steps))， 并且执行 2D 拉普拉斯逆变换。在这样得 到的扩散 - 扩散交换 2D 相关频谱内的 “交叉峰” 的出现表明在时间标度 tm 的交换。
         2) 对于一系列 tm( 通常为 4 或 8， 这样产生 16×16×8 ＝ 2048 个分开的步骤 ) 重 复在 1) 中描述的协议。通过对作为 tm 的函数的交叉峰的量的变化的分析得到 (ki+ke) 的 数值估计。
         根据 WO2008/147326 的新协议与卡拉汉的比较
         根据 WO2008/147326 的协议与卡拉汉介绍的协议在以下方面不同 ： 描述每一个 DW 块的参数是变化的， 以及分析数据的方式是不同的， 因而导致对于相同的信息内容， 更 短数量级的实验时间。在实验时间上的减小对在临床设置上的实际实现是关键的。在 WO2008/147326 中提出下面的协议 ： 1) 找到一组 {δ1， Δ1， G1} 来尽可能多地减少细胞外的部分而不影响细胞内的部 分 ( 例如， 通过基于知识的猜测或者数值方案 )。找到一组 {δ2， Δ2， G2} 来充分减少细胞 外的部分同时保持尽可能多的细胞内的部分。 基于预期的交换速率和核弛豫时间选择合理 的混合时间 tm。使用参数 {δ1， Δ 1， G 1， tm， δ2， Δ2， G2} 记录图像 1， 并利用 {δ1， Δ1， G1 ＝ 0， tm， δ2， Δ2， G2} 记录图像 2。只有在由实验变量定义的时标上存在分子交换时， 通过从图 像 2 中减去图像 1 得到的差图像产生信号强度。该协议在超过 1/100 的更短实验时间内提 供与上面的卡拉汉 1) 相同的信息。
         2) 对于一系列 tm 重复在 1) 中描述的协议。通过分析作为 tm 的函数的信号强度 的变化得到 (ki+ke) 的数值估计。再一次， 该协议在超过 1/100 的更短实验时间内提供与上 面的卡拉汉 2) 相同的信息。
         3) 对于一系列 G2 和 tm( 通常 16×7 ＝ 116 个分开的步骤 ) 利用 {δ1， Δ1， G1， tm， δ2， Δ2， G2} 重复协议。利用在 tm 为最小值时的一系列 G2 进行补充并且使用 G1 ＝ 0( 通常 16 个间距 )。根据 WO2008/147326， 这个后续系列改进了估计参数的准确性。根据下面部分 的分析产生 km、 ke 和 K 的估值。
         数据估计
         尽管协议 1) 和 2) 要求的数据估计是琐碎的——分别采用在两个图像之间的差异 和适配于一系列差图像的指数——协议 3) 的估计不是那么简单。代替新协议 1) 和 2) 以 及卡拉汉协议的更多唯象参数， 更高级的分析不只是通过估计与细胞微观结构和细胞动力 学直接相关的参数来进行合理的证明。
         例如， 对于在酵母细胞沉积物内的水， 当 DW1 块被关闭 (G1 ＝ 0) 时， 在 NMR 信号与 b( 在公式 (2) 内定义 ) 的关系的图中， 较慢 ( 细胞内 ) 的部分和较快 ( 细胞外 ) 的部分之 间存在明显差别。 该数据系列等效于利用无限长的 tm 所得到的结果。 因为信号将由于核弛 豫过程而降低到噪声水平以下， 这样的测量在实际中是不可能执行的。当打开 DW1 块时 (G1
         ＝ 0.30Tm-1)， 较快的部分消失。增加 tm 导致较快的部分再现而消耗较慢的部分。该观察明 确表明了在细胞内和细胞外的部分之间的分子交换。
         出于数据匹配的目的， 根据 WO2008/147326 便于将式 (24) 到 (26) 改写为
         其中并且
         Xe， (29) n ＝ 1-Xi， n。
         在上面的式 (27)-(29) 以及下面的所有式子中， 用于不同 tm 的数据系列的观测量 和变量 In(b)、 Xi，n、 Xe，n、 Pn(D)、 tm，n 以索引 n 标识。如上所证明的， n ＝ 0 时 G1 ＝ 0 的 系列被作为 tm ＝∞的系列处理。以 n ＝ 1 标识的数据系列指的是 tm ＝ 0。在式 (27) 中， b 涉及 DW2 块。用式 (27) 中的双指数替换公式 (24) 是近似， 只要在每一个 DW 块内 τi 和 τe 比 δ 和 Δ 长得多， 该近似预期没有问题。使用 Di、 De、 k、 Xi， Xi， 0、 1 和一组 I0， n 作为可调节参 数， 式 (27) 和式 (28) 和 (29) 适合于实验数据的整个集合。系统参数 K 和 ki 与匹配参数 k 和 Xi， 0 通过下式相关
         并且如上面描述的， 根据 WO2008/147326 的整体拟合产生所估计的参数的最准确的结果。 对于比酵母细胞沉积物更复杂的系统， 在式 (27) 中的双指数函数可以用其他多 指数表达式代替。可选择地， 人们可以通过下式将信号 In(b) 与扩散系数分布 Pn(D) 关联
         使用拉普拉斯逆变换 (ILT) 算法可以根据实验获得的 In(b) 估计 Pn(D)。可以使 用与式 (28) 类似的式子分析在得到的 Pn(D) 中的作为 tm 函数的多个部分的幅值的变化以 用于交换。不幸的是， ILT 算法由于其不稳定性而声名狼藉， 导致在 Pn(D) 中波峰的位置有 时大幅波动。 根据 WO2008/147326， 通过自定义的 ILT 算法解决该问题， 其中， 波峰位置 ( 但 不是幅值 ) 被迫使对具有 tm 的所有系列为常数。对所有的 Pn(D) 强加固定波峰位置的限制 改进了估计参数的准确性。 对于每一个 Pn(D)， 通过积分对应于细胞内和细胞外部分的波峰 的面积估计 Xi，n。在后续步骤中， 通过将式 (28) 适合于数据来估计 k、 Xi， 0 和 Xi， 1。估计的 参数与之前说明的整体拟合程序的结果有利地一致。 ILT 分析方法比整体拟合更加通用， 但 是为了改进数值稳定性， 需要更多的数据点， 因此需要更长的实验时间。
         根据本发明的具体实施方式的说明
         如同可根据上面的发明概述注意到的， 根据本发明的方法针对通过限制到信号强 度 I 的初始衰减来编码、 检测并且获取来自所述物体的磁共振信号。根据本发明， 这可以通 过在混合块之前和之后使用不同类型的加权块 ( 即， 第一和第二加权块 ) 实现。因而， 根据
         本发明的一个具体实施方式， 提供了一种根据本发明的方法， 其中第一加权块是具有弛豫 加权 τ1 的第一弛豫加权块， 并且第二加权块是具有弛豫加权 τ2 的第二弛豫块， 并且其中 通过增加 τ1 或者 τ2 将编码、 检测和获取来自所述物体的磁共振信号限制到信号强度 I 的 初始衰减来得到表观弛豫率 R ； 或者其中
         第一加权块是具有扩散加权 b1 的第一扩散加权块， 并且第二加权块是具有扩散加 权 b2 的第二扩散块， 并且其中通过增加 b1 或者 b2 将编码、 检测和获取来自所述物体的磁共 振信号限制到信号强度 I 的初始衰减来得到表观扩散系数 ADC ；
         并且其中， 其后分析表观弛豫速率 R 或者表观扩散系数 ADC 随 tm 的变化来得到表 观交换速率 AXR。
         在根据本发明的弛豫加权块的情况下， 这些可以通过弛豫加权变量 τ1 和 τ2 而不 是在扩散加权块的情况下的 b1 或 b2 来限定。弛豫率 R 可以指众所周知的核弛豫率 R1( 径 向或自旋 - 晶格 )、 R2( 横向或自旋 - 自旋 )、 或者 R1， ρ( 在旋转帧的自旋 - 晶格 ) 中的任一 种。
         根据本发明， 所得到的表观交换速率 AXR 可以用于产生图像对比度。重要的是要 理解不仅主参数 AXR， 而是还有下面解释的 σ( 西格玛 )、 ADC 和 I0， 其单独或者结合 n 集合， 可以用于产生灰度或彩色图像。作为一个例子， AXR 可以提供红色水平， ADC 提供绿色水平 并且西格玛提供蓝色水平。
         由于在临床 MRI 扫描器上有限的梯度长度， 难以获取足够高的 b 值的数据来允许 如上所述的两部分整体拟合或者 ILT 分析。临床 MRI 的另一个问题通常是高噪声水平。当 使用根据本发明的扩散加权块时， 在这种情况下可以根据下式将限制对 In(b) 的初始斜度 的分析。
         其中 给出有效衰减速率， Pn(D) 的平均值为 ：
         根据下式双位置交换 以指数方式接近 的平衡值在式 (35) 中的表达式由式 (28) 推导。
         因而， 根据本发明的一个具体实施方式， 当使用扩散加权块时， 第一加权块是具有 扩散加权 b1 的第一扩散加权块， 并且第二加权块是具有扩散加权 b2 的第二扩散块， 并且其 中根据下式取编码、 检测以及获取磁共振信号的极限 ：
         其中， 变量 b 指的是 b1 或 b2， 并且 是表观扩散系数。
         进一步的， 根据本发明的再一具体实施方式， 根据下式通过扩散系数分布 Pn(D) 的 平均值给出自扩散系数 ：
         根据本发明的再一实施方式， 随 tm， n 的变化根据下式依赖于交换速率 k ：其中 是在平衡处的 并且 是在 tm， n ＝ 0 处的 。
         甚至对于多位置交换， 式 (36) 和式 (38) 是对在 tm 的 NMR 信号的演进的好的近似。 与在传统的 DW MRI 中使用的表观扩散系数 (ADC) 的概念类似， 在多位置情况下 k 的值应当 被认为是表观交换速率 (AXR)， 其具有潜力成为在 MRI 内的对比度的有价值的模型。 为了分 析实验数据， 式 (33) 和式 (35) 被重写为
         In(b) ＝ I0， (39) nexp{-ADC·[1-σexp(-AXR·tm， n)]·b}
         其中 ADC ＝ 并且
         σ 被限制在范围 0 ＜ σ ＜ 1 内。对于宽 P(D) 和有效的 DW 滤波器得到 σ 的较大值。 根据上述公开内容， 根据本发明的一个具体实施方式， 对多位置交换进行编码、 检 测和获取。根据另一个具体实施方式， 对于多位置交换的表观交换速率 (AXR) 根据下式计 算：
         In(b) ＝ I0， (41) nexp{-ADC·[1-σexp(-AXR·tm， n)]·b}
         其中 ADC ＝ 并且
         其中 σ 被限制在范围 0 ＜ σ ＜ 1 内。
         本发明可以用于不同的物体， 诸如例如沸石、 脂质体、 泡囊或者生物细胞。
         进一步的， 当使用扩散加权块时， 限定这些块的参数可以以不同的配置组织。 根据 本发明的一个具体实施方式， 第一加权块是具有扩散加权 b1 的第一扩散加权块， 并且第二 加权块是具有扩散加权 b2 的第二扩散块， 并且其中下述中的一个是有效的 ：
         -b1 处于开或关模式， b2 固定且 tm 固定 ；
         -b2 处于开或关模式， b1 固定且 tm 固定 ；
         -b1 处于开或关模式， b2 固定且 tm 变化 ；
         -b2 处于开或关模式， b1 固定且 tm 变化 ；
         -b1 固定， b2 变化且 tm 变化， 其由 b1 等于零， tm 固定且 b2 变化的系列补充 ； 或者
         -b2 固定， b1 变化且 tm 变化， 其由 b2 等于零， tm 固定且 b1 变化的系列补充。
         术语 “处于开或关模式” 在本文意思是处于开或关模式的参数可以在这样的情况 下只有两种值， 或者零， 即关， 或别的值， 即开。在 b1 处于开或关模式、 b2 和 tm 固定的情况 下， 通过根据 (2×1×1) 进行一系列尝试 ( 即两个尝试 )， 根据本发明实现差图像和由此获
         得关于是否存在交换的信息是可能的。
         作为另一个例子， 在 b1 处于开或关模式、 b2 固定并且 tm 变化的情况下， 根据本发明 实现一系列差图像和由此获得交换速率 (ki+ke) 是可能的。
         作为再一个例子， 在 b1 固定、 b2 变化并且 tm 变化的情况下， 和然后 b1 ＝ 0、 tm 固定 并且 b2 变化的补充系列， 根据本发明从提供交换速率 (ki+ke) 和部分 Xi， 0 的全双峰匹配实现 交换速率 ki 是可能的。
         根据本发明的另一个具体实施方式， 提供了一种方法， 其中使用非对称的脉冲对 或梯度调制来提高滤波器效率。但是在扩散加权块的结尾， 有效梯度调制的时间积分应当 估计为零。
         根据本发明的一个具体实施方式， 进行数值优化来提高滤波器效率。
         根据本发明的再一个具体实施方式， 提供一种用于磁共振成像或者核磁共振光谱 的系统， 包括 ：
         射频和梯度脉冲单元， 其用于向受到磁场影响的物体发射射频和梯度脉冲序列， 其中所述物体包括具有核自旋非 0 的原子的分子 ；
         检测器单元， 其用于检测来自所述物体的对应于所述发射的射频和梯度脉冲序列 的磁共振信号 ； 其中所述射频和梯度脉冲序列包括第一加权块、 具有持续时间 tm 的混合块 和第二加权块 ；
         其中
         所述检测器单元被提供用于检测， 通过增加所述第一加权块和所述第二加权块中 的至少一个的强度将所述检测限制到所述信号强度 I 的初始衰减， 以便可以分析初始信号 衰减速率随 tm 的变化来得到表观交换速率 AXR。
         根据另一个具体的实施方式， 根据本发明提供了一种系统， 其中第一加权块是具 有弛豫加权 τ1 的第一弛豫加权块， 并且所述第二加权块是具有弛豫加权 τ2 的第二弛豫 块， 并且其中通过增加 τ1 或者 τ2 将编码、 检测以及获取来自所述物体的磁共振信号限制 到所述信号强度 I 的初始衰减来得到表观弛豫速率 R ；
         或者其中
         所述第一加权块是具有扩散加权 b1 的第一扩散加权块， 并且所述第二加权块是具 有扩散加权 b2 的第二扩散块， 并且其中通过增加 b1 或 b2 将编码、 检测以及获取来自所述物 体的磁共振信号限制到所述信号强度 I 的初始衰减来得到表观扩散系数 ADC ；
         以便以后可以分析所述表观弛豫速率 R 或者所述表观扩散系数 ADC 随 tm 的变化 来得到所述表观交换速率 AXR。
         进一步的， 根据另一具体实施方式， 根据以上所述的系统还包括图像生成单元， 其 被配置成基于所述磁共振信号生成图像。
         根据另一个实施方式， 提供了一种医疗工作站， 其包括用于执行根据本发明的方 法的装置。
         进一步的， 根据具体实施方式， 提供了根据本发明的方法、 系统或者医疗工作站的 应用， 所述医疗工作站作为诊断疾病或者疾患的诊断工具， 用于研究在体内的活细胞的新 陈代谢或者用于研究医疗药物通过细胞膜的跨膜扩散。
         虽然以上已经结合具体实施方式说明了本发明， 但是不旨在将其限制到本文阐述的具体形式。 相反地， 本发明只通过所附的权利要求来限制， 并且在这些所附的权利要求的 范围内， 上面的具体实施方式以外的其他实施方式同样是可能的。
         附图的详细说明
         图 1 示出编码用于分子扩散的 NMR 光谱或者 MR 图像的脉冲序列的示意图。信号 强度通过位于信号检测前的扩散加权块 DW 衰减。DW 块包括一对具有持续时间 σ 和幅值 G 的梯度脉冲， 其具有相反的有效极性。Δ 是在梯度脉冲的起始点之间的时间。在 DW 块的起 始和末尾的时间点被分别标记为 t0 和 t1。扩散加权 b 通过式 (3) 提供。
         图 2 示出编码用于在具有慢扩散和快扩散的部分之间的分子交换的 NMR 光谱或者 MR 图像的脉冲序列的示意图。两个扩散加权块 DW1 和 DW2 通过具有持续时间 tm 的混合块分 离。每一个 DW 块与在图 1 中示出的展开图类似。使用式 (3) 可以计算用于每一个块的扩 散加权 b。利用式 (24)、 (26)、 和 (24) 可以分别计算在时间点 t0-t1、 t1-t2 和 t2-t3 之间的 细胞内和细胞外的磁化强度 Mi 和 Me 的演进。
         图 3 示出利用根据本发明的方法在不同的溶液和细胞内执行的实验数据。在所有 的图形 a-f 中示出了归一化的 MR 信号强度与 DW2 扩散加权块的扩散加权 b 之间的关系。实 验数据被示出为具有下面的含义的符号 ： 圆形 ( 数据系列 n ＝ 0 ： DW1 扩散加权 b1 ＝ 0， tm ＝ 9 -2 29.0ms， 其对应于 tm ＝∞ )， 正方形 ( 数据系列 n ＝ 1 ： b1 ＝ 2.76·10 sm ， tm ＝ 29.0ms)， 三角形 ( 数据系列 n ＝ 2 ： b1 ＝ 2.76· 109sm-2， tm ＝ 128.0ms)， 以及菱形 ( 数据系列 n ＝ 3 ： 9 -2 b1 ＝ 2.76· 10 sm ， tm ＝ 328.0ms)。线表示产生 ADC、 AXR 和 σ 的估计的式 (36) 的全模型 拟合的结果。调查的样本是 (a) 聚 ( 乙二醇 ) 水溶液， (b) 酵母细胞， (c)MCF-10A 健康的乳 房细胞， (d) 和 (e)MCF-7 具有雌激素受体的癌变的乳腺细胞， 以及 (f)SK-BR-3 没有雌激素 受体的癌变的乳腺细胞。
         与在 WO2008/147326 中说明的之前的协议相比， 需要较小的信号衰减来获得关于 交换的信息， 其意味着可以使用对应的较低的 b 和 G 值。该协议因此可以在任何临床 MRI 扫描器上实施。较小的信号衰减还意味着该信号受到噪声的影响较小。

    举报
    举报     版权申诉
    版权申诉     PDF格式文档无特别注明外均可编辑修改；预览文档经过压缩，下载后原文更清晰！ 立即下载
    关 键  词：
    磁共振 成像 方法 系统 及其 用途
      专利查询网所有文档均是用户自行上传分享，仅供网友学习交流，未经上传用户书面授权，请勿作他用。
    0条评论

    还可以输入200字符

    暂无评论，赶快抢占沙发吧。

    关于本文
    本文标题：磁共振成像的方法和系统及其用途.pdf
    链接地址：https://www.zhuanlichaxun.net/p-1696079.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    copyright@ 2017-2018 zhuanlichaxun.net网站版权所有
    经营许可证编号:粤ICP备2021068784号-1