本发明涉及磁共振成像领域，尤其涉及一种动态测量氧摄取率的方法。
背景技术：
氧摄取率(oxygenextractionfraction，oef)描述了血液流经毛细血管床时组织从血液中摄取氧的量，即耗氧量与供氧量的比值。器官活动所需的主要能量来源于有氧代谢，因此oef可以直接反映器官的可持续性和活跃性，是器官健康和功能的重要标志。在临床应用方面，oef被用于研究脑梗，肿瘤，老年痴呆，脑缺血，贫血等与耗氧有关的疾病中。在认知神经科学方面，oef在功能刺激中的研究将有助于理解相关脑区之间的活动过程以及具体的心理阶段，因此对oef动态测量的研究不仅有利于临床应用的发展，同时还为认知神经科学的发展提供有力的技术支持。
基于磁共振成像(magneticresonanceimaging，mri)技术测量oef的方法主要是利用了血液中具有抗磁性的含氧血红蛋白被组织摄取氧后变成了具有顺磁性的脱氧血红蛋白(deoxygenatedhemoglobin，dhb)从而导致mr信号损失的原理。1994年yablonskiy和haacke提出了一种方向随机的、无限长的(相对于血管半径)圆柱形血管网络模型来研究由dhb引起的组织mr信号的变化，从而建立了oef与组织mr信号变化的关系。
测量组织mr信号变化的常用方法是采用基于梯度回波采样自旋回波(gradientechosamplingofspinecho，gesse)的成像技术。该技术因采用一次激发只填写k空间一行的采集方式，即填满k空间所有行则需要若干次激发，这就使得该成像方法的采集时间以k空间采集行数为倍数增加，同时多次激发使图像对运动非常敏感。因此，gesse成像技术并不能满足oef动态测量在时间分辨率上的要求。
技术实现要素：
为了解决采集时间较长的问题，本发明提出了一种动态测量氧摄取率的方法，基于回波平面成像(echoplanarimaging，epi)的采集模式，即一次激发利用读梯度的正反向快速切换所组成的回波链一次性填满k空间的采集方式，并利用基于自旋回波(spinecho，se)的多个时间抵消回波来测量组织mr信号变化的成像技术。由此原理，本项技术命名为多个时间抵消的自旋回波-回波平面成像(multiple-offset-spin-echoecho-planarimaging，mose-epi)。同时，mose-epi还结合了并行成像技术进一步提高了其时间分辨率，使其在时间分辨率上比gesse序列有了极大的提高。
为了实现上述目的，本发明采用以下技术方案：
一种动态测量氧摄取率的方法，包括：利用回波平面成像与多个时间抵消的自旋回波成像相结合的方式采集待测组织中的mr信号并计算mr信号衰减的可逆部分的弛豫率常数r2′以及静脉血体积分数λ，然后根据r2′和λ计算得到待测组织的氧摄取率oef，实现氧摄取率的动态测量。
进一步地，本发明还包括：在计算r2′和λ之前，对采集得到的mr信号进行预处理，所述预处理包括运动校正，图像配准以及平滑。
进一步地，本发明通过以下方法采集待测组织中的mr信号：
用一个90°选择性水激发脉冲激发待测组织特定层内水的信号，将该层的磁化矢量激发到横平面后，在自旋回波时间(echotimeofspinecho，tese)的一半处施加一个180°回聚脉冲，使得失相的磁化矢量又开始慢慢重聚，之后利用包括自旋回波时间(tese处)在内的多个不同时间抵消的自旋回波(以下简称为抵消回波，包括至少一对相对于自旋回波对称的抵消回波以及在自旋回波后的至少四个抵消回波)采集得到多个不同的mr信号。
进一步地，利用回波平面成像与多个时间抵消的自旋回波成像相结合的方式采集到的待测组织中的mr信号随时间变化的关系可写为：
其中s0为平衡态信号；λ为静脉血体积分数；r2′和r2分别为mr信号衰减的可逆部分和不可逆部分的弛豫率常数；δω为dhb引起的频率偏移(可根据前人测量的有关参数值计算得到)。
进一步地，计算mr信号衰减的可逆部分的弛豫率常数r2′包括：
根据公式(2)利用相对于自旋回波对称的第一个和第三个抵消回波计算得到mr信号衰减的不可逆部分的弛豫率常数r2，
其中，tese为自旋回波时间，δte1为回波间隔，δte1＝te2-te1＝te3-te2，te2＝tese，te1为第一个抵消回波时间，te3为第三个抵消回波时间；
接着对自旋回波后的至少四个抵消回波利用线性最小二乘法曲线拟合出r2′。
进一步地，通过公式(3)计算得到λ：
其中sextrapolated(tese)为将t＝tese代入公式(1)中的第三行式子计算得到在tese时的信号，s(tese)为实际采集得到的自旋回波信号。
进一步地，oef通过将r2′和λ代入公式(4)计算得到：
其中γ为旋磁比；δχ0为全含氧血红蛋白与全脱氧血红蛋白的磁敏感系数差，其值为0.18ppm每单位hct；hct为血红细胞比容分数(其中组织hct的典型值为0.357)；b0为主磁场强度。
本发明的有益效果如下：
1、本发明适用于所有基于yablonskiy和haacke理论的模型假设，例如：单组织模型和多组织模型；
2、利用本发明提供的一种基于多个时间抵消的自旋回波-回波平面成像的新技术——mose-epi，解决了目前因测量时间较长而无法对oef进行动态测量的问题。相较于传统测量方法gesse序列，本发明将oef的测量时间从分钟量级缩短到毫秒量级(以采集单层数据而言)，这不仅节约了时间成本，还为认知神经科学的发展提供了有力的技术支持，同时也为研究生理变化提供了新的视角。
3、本发明同时还可以实现对中间变量r2′，r2和λ的动态测量。
附图说明
图1是本发明的mose-epi磁共振序列示意图。
图2是双手曲直运动刺激下得到的oef激活图。
图3是以主要运动区域中的激活区为感兴趣区得到的oef相对变化时序图。
具体实施方式
为使本发明的特点更为浅显易懂，下文采用具有代表性的较佳实施例，并配合所附图进行详细说明。应指出的是，本发明并不局限于所描述的具体结构、功能、器件和方法，也可以具有其他实施方式，或者是其他实施方式的组合。本发明中所描述的元素数目也可以设置为多个。此外，为避免其他例与本发明发生混淆，对于本领域中众所周知的一些技术特征和细节未进行描述。
本发明提供了一种基于多个时间抵消的自旋回波-回波平面成像的新序列，命名为mose-epi(multiple-offset-spin-echoecho-planarimaging)。该序列因其具有较高的时间分辨率使其能够应用于功能磁共振成像(functionalmagneticresonanceimaging，fmri)的研究中。本例采用区块设计(blockdesign)的磁共振实验方法，以双手曲直运动为实验任务，以mose-epi序列为测量手段，详细说明该序列的具体实施方法。值得强调的是，该方法并不局限于运动任务，其他类型实验任务可以按照相似的步骤处理。
图1是本发明的mose-epi序列的示意图。其中，rf为射频脉冲，gx,y为频率编码梯度和相位编码梯度，gz为层面选择梯度，这里以简单的六个回波为例，即一对对称于自旋回波的抵消回波以及四个位于自旋回波后的抵消回波，te2＝tese，δte1和δte2为回波间隔，其中δte1＝te2-te1＝te3-te2，δte2＝te4-te3＝te5-te4＝te6-te5。
在本例中采用的扫描参数：视野大小为26厘米；矩阵大小为64×64；重复时间(repetitiontime，tr)为2000毫秒；层厚为6毫米，以间隔采集的方式采集了10层以覆盖运动区；tese为100毫秒；δte1为35毫秒；δte2为12.4毫秒。对于上述扫描参数的设置可以任意选取，这里只是针对该例子的较佳实施方案。
具体实验任务为12秒的空扫(为了使信号达到稳态)，之后为三个休息块和两个手动任务块交替出现，利用同步刺激仪在实验任务开始时同步开始采集mose-epi信号，根据以上序列参数设置每2秒可得到六个不同回波时间的组织信号(为方便理解，称为一组信号)，本实验共采集了162秒的数据(即81组信号)。
根据yablonskiy和haacke理论模型，上述序列采集到的信号随时间变化的关系为：
其中s0为平衡态信号；λ为静脉血体积分数；r2′和r2分别为mr信号衰减的可逆部分和不可逆部分的弛豫率常数；δω为dhb引起的频率偏移。
将最开始采集到的12秒空扫数据(即六组信号)去掉，剩下的75组数据利用spm12软件对其经过一系列的包括运动校正，图像配准以及平滑等图像预处理，然后针对每一组数据中的六个不同回波信号进行以下计算：
根据公式(2)利用第一个和第三个抵消回波计算得到r2，
具体计算过程为：将t＝tese-δte1＝te1代入公式(1)中的第一行式子，t＝tese+δte1＝te3代入公式(1)中的第三行式子，结合起来即可抵消掉未知参数，从而计算得到r2；
接着利用后四个抵消回波利用线性最小二乘法曲线拟合出r2′，并结合由第二个抵消回波得到的自旋回波信号s(tese)根据公式(3)计算得到λ，
其中sextrapolated(tese)为将t＝tese代入公式(1)中的第三式计算得到在tese时的信号，s(tese)为实际采集得到的自旋回波信号；
最后将r2′和λ代入公式(4)求得每个体素的oef，即一组数据得到一幅oef图，
其中γ为旋磁比；δχ0为全含氧血红蛋白与全脱氧血红蛋白的磁敏感系数差，其值为0.18ppm每单位hct；hct为血红细胞比容分数(其中组织hct的典型值为0.357)；b0为主磁场强度。
图2是双手曲直运动刺激下得到的oef激活图。该图对每个体素的oef时间序列(75个时间点)通过一般线性模型进行激活分析，并采用体素水平p<0.001，团簇水平p<0.05的fwe校正得到的激活图，从图中可以看出激活区都在运动区域，包括初级运动皮层，感觉运动皮层和辅助运动皮层。这与fmri相关研究的结果保持一致，且该激活区比常用的bold(blood-oxygen-level-dependence)技术测量更集中在灰质区域，由此可见本方法是可行的。
图3是以主要运动区域中的激活区为感兴趣区得到的oef相对变化时序图。从图中可以看出oef在任务态下是降低的，这也符合生理机制，且oef的平均相对变化大小也与其他图像技术，如：正电子发射计算机断层显像技术和扩散光学成像得到的结果相同，相比于这些图像技术，磁共振的优势已经是众所周知，在此就不详细讨论了。同时在静息态下，得到的数值大小(0.349)也与其他诸多oef定量测量得到的结果保持一致，由此进一步说明了，本发明的优势及可行性。
综上所述，利用本发明提供的一种基于多个时间抵消的自旋回波-回波平面成像的新技术——mose-epi，解决了目前因测量时间较长而无法对oef进行动态测量的问题。相较于传统测量方法gesse序列，本发明将oef的测量时间从分钟量级缩短到毫秒量级(以采集单层数据而言)，这不仅节约了时间成本，还为认知神经科学的发展提供了有力的技术支持，同时也为研究生理变化提供了新的视角。