CN103519809A - 氧代谢参数估测方法和系统 - Google Patents

Abstract

本发明涉及一种氧代谢参数估测方法和系统。所述方法，包括以下步骤：图像获取步骤，获取磁敏感度加权成像图或定量的磁敏感度图；相位差及磁敏感度变化计算步骤，从所述磁敏感度加权成像图或定量的磁敏感度图确定血管位置，计算血管与周围组织的相位差，并根据所述血管与周围组织的相位差计算组织之间局部的磁敏感度变化；氧饱和估测步骤，根据所述组织之间局部的磁敏感度变化估测出氧饱和。上述氧代谢参数估测方法和系统，通过相位差可求得磁敏感度变化，根据磁敏感度变化可求取氧饱和参数，该估测方法无侵害，无任何副作用，提高了氧饱和参数估测的安全性。

Description

氧代谢参数估测方法和系统

技术领域

本发明涉及生理参数检测领域，特别是涉及一种氧代谢参数估测方法和系统。

背景技术

脑的氧代谢对维持大脑正常功能运转起着关键的作用，它是大脑神经功能的主要能量源。为了了解氧代谢的情况，传统的方式有采用正电子放射断层造影术测量氧摄取分数（Oxygen Extraction Fraction，OEF）、静脉氧饱和（Yv）、大脑的氧代谢率（Cerebral Metabolic Rate Of Oxygen，CMRO₂）等参数，然而，需要对被试者注射外源性对比剂，且可能需要多次注射，因外源性对比剂一般为放射性元素，对人体会产生不良影响，影响了人体的健康。

发明内容

基于此，有必要针对现有的氧代谢参数检查存在安全隐患的问题，提供一种安全性高的氧代谢参数估测方法。

此外，还有必要提供一种安全性高的氧代谢参数估测系统。

一种氧代谢参数估测方法，包括以下步骤：

图像获取步骤，获取磁敏感度加权成像图或定量的磁敏感度图；

相位差及磁敏感度变化计算步骤，从所述磁敏感度加权成像图或定量的磁敏感度图确定血管位置，计算血管与周围组织的相位差，并根据所述血管与周围组织的相位差计算组织之间局部的磁敏感度变化；

氧饱和估测步骤，根据所述组织之间局部的磁敏感度变化估测出氧饱和。

在其中一个实施例中，所述图像获取步骤中获取磁敏感度加权成像图的步骤包括：

采集步骤，采用梯度回波成像，得到原始相位图和原始幅值图；

校正步骤，将原始相位图除以低通滤波后的k空间数据，得到校正后的相位图；

标准化步骤，将校正后的相位图中相位值进行标准化处理，得到相位掩膜图；

合成步骤，将所述相位掩膜图作预设值的次幂后和原始幅值图相乘得到磁敏感度加权成像图。

在其中一个实施例中，在所述合成步骤之后，还包括：

对得到的磁敏感度加权成像图采用最小强度投影进行处理。

在其中一个实施例中，所述相位差及磁敏感度变化计算步骤中从所述磁敏感度加权成像图确定血管位置，计算血管与周围组织的相位差，并根据所述血管与周围组织的相位差计算组织之间局部的磁敏感度变化的公式包括：

其中，为血管与周围组织的相位差，γ为旋磁比，ΔB为组织间的磁场变化，TE为回波时间，Δχ为组织之间局部的磁敏感度变化，B₀为主磁场场强；

所述氧代谢参数估测方法还包括：

将所述血管与周围组织的相位差的计算公式简化处理得到：

其中，θ为血管与主磁场B₀之间的夹角。

在其中一个实施例中，所述相位差及磁敏感度变化计算步骤中从所述定量的磁敏感度图确定血管位置，计算血管与周围组织的相位差，并根据所述血管与周围组织的相位差计算组织之间局部的磁敏感度变化的公式包括：

δ＝F^-1DFΔχ;

$D=\frac{1}{3}-\frac{{k_z}^2}{k^2};$

其中，δ为预设的标准场图，F为傅里叶变化，D是k空间的敏感度核，k_z表示z轴数据，为血管与周围组织的相位差，γ为旋磁比，TE为回波时间，Δχ为组织之间局部的磁敏感度变化，B₀为主磁场场强。

在其中一个实施例中，所述氧饱和估测步骤中根据所述组织之间局部的磁敏感度变化估测出氧饱和的公式包括：

Δχ＝Δχ₀·Hct·(1-Y)

其中，Δχ为组织之间局部的磁敏感度变化，Δχ₀为每单元的血细胞比容中所有去氧血和含氧血的磁敏感度变化，Hct为分数比容，Y为氧饱和。

一种氧代谢参数估测系统，包括：

图像获取模块，用于获取磁敏感度加权成像图或定量的磁敏感度图；

相位差及磁敏感度变化计算模块，用于从所述磁敏感度加权成像图或定量的磁敏感度图确定血管位置，计算血管与周围组织的相位差，并根据所述血管与周围组织的相位差计算组织之间局部的磁敏感度变化；

氧饱和估测模块，用于根据所述组织之间局部的磁敏感度变化估测出氧饱和。

在其中一个实施例中，所述图像获取模块包括：

采集子模块，用于采用梯度回波成像，得到原始相位图和原始幅值图；

校正子模块，用于将原始相位图除以低通滤波后的k空间数据，得到校正后的相位图；

标准化子模块，用于将校正后的相位图中相位值进行标准化处理，得到相位掩膜图；

合成子模块，用于将所述相位掩膜图作预设值的次幂后和原始幅值图相乘得到磁敏感度加权成像图。

在其中一个实施例中，所述图像获取模块还包括：

调整子模块，用于对得到的磁敏感度加权成像图采用最小强度投影进行处理。

在其中一个实施例中，所述相位差及磁敏感度变化计算模块从所述磁敏感度加权成像图确定血管位置，计算血管与周围组织的相位差，并根据所述血管与周围组织的相位差计算组织之间局部的磁敏感度变化的公式包括：

其中，为血管与周围组织的相位差，γ为旋磁比，ΔB为组织间的磁场变化，TE为回波时间，Δχ为组织之间局部的磁敏感度变化，B₀为主磁场场强；

所述氧代谢参数估测方法还包括：

将所述血管与周围组织的相位差的计算公式简化处理得到：

其中，θ为血管与主磁场B₀之间的夹角。

在其中一个实施例中，所述相位差及磁敏感度变化计算模块从所述定量的磁敏感度图确定血管位置，计算血管与周围组织的相位差，并根据所述血管与周围组织的相位差计算组织之间局部的磁敏感度变化的公式包括：

δ＝F^-1DFΔχ;

$D=\frac{1}{3}-\frac{{k_z}^2}{k^2};$

其中，δ为预设的标准场图，F为傅里叶变化，D是k空间的敏感度核，k_z表示z轴数据，

为血管与周围组织的相位差，γ为旋磁比，TE为回波时间，Δχ为组织之间局部的磁敏感度变化，B₀为主磁场场强。

在其中一个实施例中，所述氧饱和估测模块根据所述组织之间局部的磁敏感度变化估测出氧饱和的公式包括：

Δχ＝Δχ₀·Hct·(1-Y)

其中，Δχ为组织之间局部的磁敏感度变化，Δχ₀为每单元的血细胞比容中所有去氧血和含氧血的磁敏感度变化，Hct为分数比容，Y为氧饱和。

上述氧代谢参数估测方法和系统，在磁敏感度加权成像图或定量的磁敏感度图中使得T2*衰减增强，导致幅值图中信号减弱，使得组织与血管相位差明显，通过相位差可求得磁敏感度变化，根据磁敏感度变化可求取氧饱和参数，该估测方法无侵害，无任何副作用，提高了氧饱和参数估测的安全性。

附图说明

图1为一个实施例中氧代谢参数估测方法的流程图；

图2为获取磁敏感度加权成像图的流程图；

图3a的原始幅值图；

图3b的原始相位图；

图4为对原始相位图校正的示意图；

图5为负的相位掩膜图；

图6为合成SWI图像的示意图；

图7为通过MIP处理后的图像；

图8为通过相位变化表示组织和血管之间的磁敏感度变化；

图9a为SWI的相位图；

图9b为图9a中相位曲线图；

图10为一个实施例中氧代谢参数估测系统的结构框图；

图11为图10中图像获取模块的内部结构框图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

如图1所示，为一个实施例中氧代谢参数估测方法的流程图。该氧代谢参数估测方法，包括以下步骤：

步骤102，图像获取步骤，获取磁敏感度加权成像图或定量的磁敏感度图；

步骤104，相位差及磁敏感度变化计算步骤，从该磁敏感度加权成像图或定量的磁敏感度图确定血管位置，计算血管与周围组织的相位差，并根据该血管与周围组织的相位差计算组织之间局部的磁敏感度变化；

步骤106，氧饱和估测步骤，根据该组织之间局部的磁敏感度变化估测出氧饱和。

上述氧代谢参数估测方法中，在磁敏感度加权成像图或定量的磁敏感度图中使得T2*衰减增强，导致幅值图中信号减弱，使得组织与血管相位差明显，通过相位差可求得磁敏感度变化，根据磁敏感度变化可求取氧饱和参数，该估测方法无侵害，无任何副作用，提高了氧饱和参数估测的安全性。

在一个实施例中，如图2所示，该图像获取步骤中获取磁敏感度加权成像图的步骤包括：

步骤202，采集步骤，采用梯度回波成像，得到原始相位图和原始幅值图。

具体的，磁敏感加权成像（Susceptibility Weighted Imaging，SWI）是一个高分辨率的相位对比增强的梯度回波成像方法，其在三个成像方向上采用流速补偿。SWI对静脉氧，出血以及铁过量较为敏感。本实施例中，采用梯度回波序列成像，得到原始相位图和原始幅值图，如图3a的原始幅值图和图3b的原始相位图。该梯度回波时间TE为18至25毫秒。

步骤204，校正步骤，将原始相位图除以低通滤波后的k空间数据，得到校正后的相位图。

具体的，对原始相位图施加一个低通滤波器，然后在复数域中用原始相位图除以低通滤波后的k空间数据，去除由于背景磁场不均匀造成的低频扰动，得到高通滤过的相位图像，即校正后的相位图。如图4所示，原始的相位图402除以低通滤波的k空间数据图404，得到高通滤过的相位图像406。

步骤206，标准化步骤，将校正后的相位图中相位值进行标准化处理，得到相位掩膜图。

具体的，将校正后的相位图中不同组织的相位值进行标准化处理。本实施例中，采用负的相位掩膜，将小于0的相位值

转化到（0，1）区间，大于0的设为1，如图5所示为负的相位掩膜图。因在相位图中，顺磁物质如静脉的相位信号表现为明显的负值，而脑实质（如大部分脑实质及脑脊液等）相位值通常为正值或较小的负值，因此静脉等顺磁物质在相位图中的相位值转化为（0，1）。

步骤208，合成步骤，将该相位掩膜图作预设值的次幂后和原始幅值图相乘得到磁敏感度加权成像图。

如图6所示，相位掩膜图602作n次幂后和原始幅值图604相乘得到SWI图像606。其中，n一般取3～5可得到信噪比较高的图像。其中，相位掩膜图的函数表示为

原始幅值图函数表示为ρ(x)。

步骤210，对得到的磁敏感度加权成像图采用最小强度投影进行处理。

具体的，采用最小强度投影（Minimum Intensity Project，MIP）使得分散在各个层面的静脉信号连续化，显示连续的静脉血管结构，如图7所示。通过SWI图像进行处理提高了磁矩图像的对比，对静脉血、出血和铁沉积高度敏感，甚至可检测到一个体素的静脉血管。

在其他实施例中，步骤210可省略。

因静脉血管中的脱氧血红蛋白与周围组织存在着磁敏感度差异，会产生局部的磁场变化ΔB，局部的磁场变化ΔB决定了组织内部和外部的空间相位变化，该种场的变化引起两种反应：（1）T2*变短；（2）存在可测量的相位变化值，即相位差。SWI幅值图可反映增强T2*和引起相位变化的敏感性差异，故此，在一个实施例中，对于一个给定的血管位置，从该磁敏感度加权成像图确定血管位置，计算血管与周围组织的相位差，并根据该血管与周围组织的相位差计算组织之间局部的磁敏感度变化的公式包括：

其中，

为血管与周围组织的相位差，γ为旋磁比，ΔB为组织间的磁场变化，TE为回波时间，Δχ为组织之间局部的磁敏感度变化，B₀为主磁场场强；γ为2.678×10⁸rad/s/T（弧度/秒/特斯拉）。

假设血管是一个有限的圆柱体，该氧代谢参数估测方法还包括：

将该血管与周围组织的相位差的计算公式简化处理得到：

其中，θ为血管与主磁场B₀之间的夹角。

该氧饱和估测步骤中根据该组织之间局部的磁敏感度变化估测出氧饱和的公式包括：

Δχ＝Δχ₀·Hct·(1-Y)    （3）

其中，Δχ为组织之间局部的磁敏感度变化，Δχ₀为每单元的血细胞比容中所有去氧血和含氧血的磁敏感度变化，Hct为分数比容（在给定的大脑组织中红细胞的百分比），Y为氧饱和。其中，Hct一般为0.4，Δχ₀为0.18ppm/per（百万分比）unit Hct

随着时间变化的相对氧饱和变化ΔY的公式如下：

其中，小标1和0表示不同时刻的MRI（Magnetic Resonance Imaging，磁共振成像）扫描，

和为两个不同时刻的相同血管与周围组织的相位差，

可通过公式（5）计算得到：

其中，表示血管一边的尖峰位置a的相位值，

表示血管另一边的尖峰位置b的相位值，

表示血管中相位值最小的位置c的相位值。如图8所示，为通过相位变化表示组织和血管之间的磁敏感度变化。当出现多个尖峰时，公式（5）可求取均值，即

表示血管一边的所有尖峰位置a的相位均值，

表示血管另一边的所有尖峰位置b的相位均值。

图9a为SWI的相位图；图9b为图9a中相位曲线图。

在一个实施例中，定量的磁敏感度图中，磁敏感度变化采用矩阵表示即可定义一个未知方程，

δ＝F^-1DFΔχ    （6）

其中，δ为预设的标准场图，F为傅里叶变化，D是k空间的敏感度核，Δχ为组织之间局部的磁敏感度变化。

$D=\frac{1}{3}-\frac{{k_z}^2}{k^2}---\left(8\right)$

其中，D是k空间的敏感度核，k_z表示z轴数据，

为血管与周围组织的相位差，γ为旋磁比，TE为回波时间，B₀为主磁场场强。

该相位差及磁敏感度变化计算步骤中从该定量的磁敏感度图确定血管位置，计算血管与周围组织的相位差，并根据该血管与周围组织的相位差计算组织之间局部的磁敏感度变化的公式包括公式（6）（7）和（8）。

通过公式（6）（7）和（8）计算得到磁敏感度变化后，再通过公式（3）计算得到氧饱和Y。

如图10所示，为一个实施例中氧代谢参数估测系统的结构框图。该氧代谢参数估测系统，包括图像获取模块1020、相位差及磁敏感度变化计算模块1040和氧饱和估测模块1060。其中：

图像获取模块1020，用于获取磁敏感度加权成像图或定量的磁敏感度图。

相位差及磁敏感度变化计算模块1040，用于从该磁敏感度加权成像图或定量的磁敏感度图确定血管位置，计算血管与周围组织的相位差，并根据该血管与周围组织的相位差计算组织之间局部的磁敏感度变化。

氧饱和估测模块1060，用于根据该组织之间局部的磁敏感度变化估测出氧饱和。

上述氧代谢参数估测系统，在磁敏感度加权成像图或定量的磁敏感度图中使得T2*衰减增强，导致幅值图中信号减弱，使得组织与血管相位差明显，通过相位差可求得磁敏感度变化，根据磁敏感度变化可求取氧饱和参数，该估测方法无侵害，无任何副作用，提高了氧饱和参数估测的安全性。

如图11所示，该图像获取模块1020包括采集子模块1022、校正子模块1024、标准化子模块1026、合成子模块1028和调整子模块1029。其中：

采集子模块1022，用于采用梯度回波成像，得到原始相位图和原始幅值图；

校正子模块1024，用于将原始相位图除以低通滤波后的K空间数据，得到校正后的相位图；

标准化子模块1026，用于将校正后的相位图中相位值进行标准化处理，得到相位掩膜图；

合成子模块1028，用于将该相位掩膜图作预设值的次幂后和原始幅值图相乘得到磁敏感度加权成像图。

调整子模块1029，用于对得到的磁敏感度加权成像图采用最小强度投影进行处理。

在一个实施例中，该相位差及磁敏感度变化计算模块1040从该磁敏感度加权成像图确定血管位置，计算血管与周围组织的相位差，并根据该血管与周围组织的相位差计算组织之间局部的磁敏感度变化的公式包括：

其中，

为血管与周围组织的相位差，γ为旋磁比，ΔB为组织间的磁场变化，TE为回波时间，Δχ为组织之间局部的磁敏感度变化，B₀为主磁场场强；

该氧代谢参数估测方法还包括：

将该血管与周围组织的相位差的计算公式简化处理得到：

其中，θ为血管与主磁场B₀之间的夹角。

在一个实施例中，该相位差及磁敏感度变化计算模块1040从该定量的磁敏感度图确定血管位置，计算血管与周围组织的相位差，并根据该血管与周围组织的相位差计算组织之间局部的磁敏感度变化的公式包括：

δ＝F^-1DFΔχ;

$D=\frac{1}{3}-\frac{{k_z}^2}{k^2};$

其中，δ为预设的标准场图，F为傅里叶变化，D是k空间的敏感度核，k_z表示z轴数据，

为血管与周围组织的相位差，γ为旋磁比，TE为回波时间，Δχ为组织之间局部的磁敏感度变化，B₀为主磁场场强。

该氧饱和估测模块1060根据该组织之间局部的磁敏感度变化估测出氧饱和的公式包括：

Δχ＝Δχ₀·Hct·(1-Y)

其中，Δχ为组织之间局部的磁敏感度变化，Δχ₀为每单元的血细胞比容中所有去氧血和含氧血的磁敏感度变化，Hct为分数比容，Y为氧饱和。

本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述的存储介质可为磁碟、光盘、只读存储记忆体（Read-Only Memory，ROM）或随机存储记忆体（Random Access Memory，RAM）等。

以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。

Claims (12)

1.一种氧代谢参数估测方法，包括以下步骤：

图像获取步骤，获取磁敏感度加权成像图或定量的磁敏感度图；

相位差及磁敏感度变化计算步骤，从所述磁敏感度加权成像图或定量的磁敏感度图确定血管位置，计算血管与周围组织的相位差，并根据所述血管与周围组织的相位差计算组织之间局部的磁敏感度变化；

氧饱和估测步骤，根据所述组织之间局部的磁敏感度变化估测出氧饱和。

2.根据权利要求1所述的氧代谢参数估测方法，其特征在于，所述图像获取步骤中获取磁敏感度加权成像图的步骤包括：

采集步骤，采用梯度回波成像，得到原始相位图和原始幅值图；

校正步骤，将原始相位图除以低通滤波后的k空间数据，得到校正后的相位图；

标准化步骤，将校正后的相位图中相位值进行标准化处理，得到相位掩膜图；

合成步骤，将所述相位掩膜图作预设值的次幂后和原始幅值图相乘得到磁敏感度加权成像图。

3.根据权利要求2所述的氧代谢参数估测方法，其特征在于，在所述合成步骤之后，还包括：

对得到的磁敏感度加权成像图采用最小强度投影进行处理。

4.根据权利要求1所述的氧代谢参数估测方法，其特征在于，所述相位差及磁敏感度变化计算步骤中从所述磁敏感度加权成像图确定血管位置，计算血管与周围组织的相位差，并根据所述血管与周围组织的相位差计算组织之间局部的磁敏感度变化的公式包括：

其中，

为血管与周围组织的相位差，γ为旋磁比，ΔB为组织间的磁场变化，TE为回波时间，Δχ为组织之间局部的磁敏感度变化，B₀为主磁场场强；

所述氧代谢参数估测方法还包括：

将所述血管与周围组织的相位差的计算公式简化处理得到：

其中，θ为血管与主磁场B₀之间的夹角。

5.根据权利要求1所述的氧代谢参数估测方法，其特征在于，所述相位差及磁敏感度变化计算步骤中从所述定量的磁敏感度图确定血管位置，计算血管与周围组织的相位差，并根据所述血管与周围组织的相位差计算组织之间局部的磁敏感度变化的公式包括：

δ＝F^-1DFΔχ;

$D=\frac{1}{3}-\frac{{k_z}^2}{k^2};$

其中，δ为预设的标准场图，F为傅里叶变化，D是k空间的敏感度核，k_z表示z轴数据，为血管与周围组织的相位差，γ为旋磁比，TE为回波时间，Δχ为组织之间局部的磁敏感度变化，B₀为主磁场场强。

6.根据权利要求1所述的氧代谢参数估测方法，其特征在于，所述氧饱和估测步骤中根据所述组织之间局部的磁敏感度变化估测出氧饱和的公式包括：

Δχ＝Δχ₀·Hct·(1-Y)

其中，Δχ为组织之间局部的磁敏感度变化，Δχ₀为每单元的血细胞比容中所有去氧血和含氧血的磁敏感度变化，Hct为分数比容，Y为氧饱和。

7.一种氧代谢参数估测系统，其特征在于，包括：

图像获取模块，用于获取磁敏感度加权成像图或定量的磁敏感度图；

相位差及磁敏感度变化计算模块，用于从所述磁敏感度加权成像图或定量的磁敏感度图确定血管位置，计算血管与周围组织的相位差，并根据所述血管与周围组织的相位差计算组织之间局部的磁敏感度变化；

氧饱和估测模块，用于根据所述组织之间局部的磁敏感度变化估测出氧饱和。

8.根据权利要求7所述的氧代谢参数估测系统，其特征在于，所述图像获取模块包括：

采集子模块，用于采用梯度回波成像，得到原始相位图和原始幅值图；

校正子模块，用于将原始相位图除以低通滤波后的k空间数据，得到校正后的相位图；

标准化子模块，用于将校正后的相位图中相位值进行标准化处理，得到相位掩膜图；

合成子模块，用于将所述相位掩膜图作预设值的次幂后和原始幅值图相乘得到磁敏感度加权成像图。

9.根据权利要求8所述的氧代谢参数估测系统，其特征在于，所述图像获取模块还包括：

调整子模块，用于对得到的磁敏感度加权成像图采用最小强度投影进行处理。

10.根据权利要求7所述的氧代谢参数估测系统，其特征在于，所述相位差及磁敏感度变化计算模块从所述磁敏感度加权成像图确定血管位置，计算血管与周围组织的相位差，并根据所述血管与周围组织的相位差计算组织之间局部的磁敏感度变化的公式包括：

其中，为血管与周围组织的相位差，γ为旋磁比，ΔB为组织间的磁场变化，TE为回波时间，Δχ为组织之间局部的磁敏感度变化，B₀为主磁场场强；

所述氧代谢参数估测方法还包括：

将所述血管与周围组织的相位差的计算公式简化处理得到：

其中，θ为血管与主磁场B₀之间的夹角。

11.根据权利要求7所述的氧代谢参数估测系统，其特征在于，所述相位差及磁敏感度变化计算模块从所述定量的磁敏感度图确定血管位置，计算血管与周围组织的相位差，并根据所述血管与周围组织的相位差计算组织之间局部的磁敏感度变化的公式包括：

δ＝F^-1DFΔχ;

其中，δ为预设的标准场图，F为傅里叶变化，D是k空间的敏感度核，k_z表示z轴数据，

为血管与周围组织的相位差，γ为旋磁比，TE为回波时间，Δχ为组织之间局部的磁敏感度变化，B₀为主磁场场强。

12.根据权利要求7所述的氧代谢参数估测系统，其特征在于，所述氧饱和估测模块根据所述组织之间局部的磁敏感度变化估测出氧饱和的公式包括：

Δχ＝Δχ₀·Hct·(1-Y)

其中，Δχ为组织之间局部的磁敏感度变化，Δχ₀为每单元的血细胞比容中所有去氧血和含氧血的磁敏感度变化，Hct为分数比容，Y为氧饱和。

Images