CN104749540A - 一种指示组织间区分度的方法与装置 - Google Patents

Abstract

本发明提供了一种指示组织间区分度的方法。对多种组织中的每种组织，根据随机扫描序列，生成对应所述每种组织的一组磁化矢量；根据所述多种组织中每种组织所对应的磁化矢量，计算得到所述多种组织两两之间的区分度指示值。这样，医生可以根据本发明提供的方法，选择适合的随机扫描序列，生成包括多种组织对应的多条亮度曲线的磁共振图像，使得在该磁共振图像中多种组织分别对应的亮度曲线的趋势相差较大。本发明还提供了一种指示组织间区分度的装置。

Description

一种指示组织间区分度的方法与装置

技术领域

本发明涉及磁共振领域，特别是一种指示组织间区分度的方法与装置。

背景技术

磁共振技术是一种强大、多用途的测量技术。磁共振成像(MRI)技术通过利用磁场和无线电波脉冲来生成身体组织和结构的图像，可以定量的对被该磁共振技术所介入的一组特征进行观察。磁共振指纹法(magnetic resonance fingerprinting，MRF)是一种MRI新方法，可以早期常规筛查某些特异的癌症、多发性硬化症、心脏病及其他疾病。每个身体组织和疾病都具有一种独特的指纹，可用于快速诊断问题。通过利用随机扫描序列，MRF技术可以同时扫描多种具有不同的物理特性的组织，生成包括多种组织的亮度曲线的磁共振图像。这样，医生可以利用该包括了多种组织亮度曲线的磁共振图像进行疾病的诊断。

理论上，当磁共振图像中的多条亮度曲线的趋势相差较大时，医生能够在一幅图像中，准确地分辨出多种组织的特性，因而，能够在该图像中同时观察、对比、分析该多种组织。反之，当磁共振图像中的多条亮度曲线的趋势相近时，则无法从该图像中分离多个组织的特征，这样，该磁共振图像便无法辅助医生进行疾病的诊断。然而，针对不同的组织，采用不同的随机扫描序列，通过MRF方法所生成的一幅磁共振图像中组织的多条亮度曲线的趋势相差也不相同。目前尚无法预先得知针对某些组织，采用哪个随机扫描序列能够生成包括这些组织所对应的亮度曲线的趋势相差较大的磁共振图像。只能针对多种组织使用多个不同的随机序列，生成多幅磁共振图像，从生成的多幅磁共振图像中选择该多种组织对应的亮度曲线趋势相差较大的图像。这样，将浪费大量的处理资源。

发明内容

有鉴于此，本发明提出了一种指示组织区分度的方法，用以指示利用随机扫描序列生成在磁共振图像中，多种组织分别对应的多条亮度曲线的区分度。

一种指示组织间区分度的方法，包括：

对多种组织中的每种组织，根据随机扫描序列，生成对应所述每种组织的一组磁化矢量；

根据所述多种组织中每种组织所对应的磁化矢量，计算得到所述多种组织两两之间的区分度指示值。

优选地，所述对多种组织中的每种组织，根据随机扫描序列，生成对应所述组织的一组磁化矢量，包括：

根据所述随机扫描序列，生成每种组织的射频脉冲旋转矩阵R_RF；

根据所述随机扫描序列，生成每种组织的弛豫矩阵R_relax；

判断所述多种组织中是否具有偏共振特性df(Hz)的组织；

若是，则根据所述随机扫描序列，生成每种具有偏共振特性的组织的偏共振旋转矩阵R_off，将所述射频脉冲旋转矩阵R_RF、所述弛豫矩阵R_relax、所述偏共振旋转矩阵R_off作用于初始磁化矢量，通过迭代计算，生成每种组织的磁化矢量；

否则，将所述射频脉冲旋转矩阵R_RF、所述弛豫矩阵R_relax作用于初始磁化矢量，通过迭代计算，生成每种组织的所述磁化矢量。

优选地，所述随机扫描序列包括如下随机扫描序列中的至少一种：翻转角随机扫描序列和重复时间随机扫描序列。

优选地，根据所述随机扫描序列，生成所述组织的射频脉冲旋转矩阵R_RF，包括：

判断所述随机扫描序列是否包括所述翻转角随机扫描序列；

若所述随机扫描序列包括所述翻转角随机扫描序列，则根据所述翻转角随机扫描序列，生成每种组织的所述射频脉冲旋转矩阵R_RF；

若所述随机扫描序列不包括翻转角随机扫描序列，则根据预设的翻转角非随机扫描序列，生成每种组织的所述射频脉冲旋转矩阵R_RF。

优选地，根据所述随机扫描序列，生成所述组织的弛豫矩阵R_relax，包括：

判断所述随机扫描序列是否包括重复时间随机扫描序列；

若所述随机扫描序列包括重复时间随机扫描序列，则根据所述重复时间随机扫描序列，生成每种组织的一个弛豫矩阵R_relax；

若所述随机扫描序列不包括重复时间随机扫描序列，则根据预设的重复时间非随机扫描序列，生成每种组织的一个弛豫矩阵R_relax。

优选地，所述根据所述随机扫描序列，生成每种具有偏共振特性的组织的偏共振旋转矩阵R_off，包括：

若所述随机扫描序列包括所述重复时间随机扫描序列，则根据每种具有偏共振特性的组织的偏共振特性df(Hz)与所述重复时间随机扫描序列，得到每种具有偏共振特性的组织的偏共振旋转矩阵R_off；

若所述随机扫描序列不包括所述重复时间随机扫描序列，则根据每种具有偏共振特性的组织的偏共振特性df(Hz)与所述预设的重复时间非随机扫描序列，得到每种具有偏共振特性的组织的偏共振旋转矩阵R_off。

优选地，所述根据所述多种组织中每种组织所对应的磁化矢量，计算得到所述多种组织的区分度指示值，包括：

根据每种组织所对应的磁化矢量，获得所述组织的一组织磁共振指纹演化矢量，每种组织磁共振指纹演化矢量包括所述每种组织对应的磁化矢量中每个元素分别投影在XOY平面上矢量的模值；

将所述多种组织中的任意两种组织的所述组织磁共振指纹演化矢量进行互相关计算，得到任意两种组织的区分度指示值。

优选地，本方法进一步包括：

根据所述多种组织的区分度指示值，生成多种组织的区分度指示矩阵，所述多种组织的区分度指示矩阵中i行j列的元素为第i种组织与第j种组织的区分度指示值。

优选地，本方法还包括，根据所述区分度指示值，确定利用所述随机扫描序列生成的磁共振图像中所述多种组织分别对应的亮度曲线的两两之间的区分度。

优选地，所述根据所述区分度指示值，确定利用所述随机扫描序列生成的磁共振图像中所述多种组织分别对应的亮度曲线的两两之间的区分度，包括：

判断所述多种组织中的任意两种组织的区分度指示值是否小于设定的阈值；

若所述任意两种组织的区分度指示值小于设定的阈值，则确定在一幅图像中，所述任意两种组织分别对应的两条亮度曲线的区分度高。

一种指示组织间区分度的装置，包括：

磁化矢量生成模块，用于对多种组织中的每种组织，根据随机扫描序列，生成对应每种组织的一组磁化矢量；

组织区分度指示获得模块，用于根据所述多种组织中每种组织所对应的磁化矢量，计算得到所述多种组织两两之间的区分度指示值；

优选地，磁化矢量生成模块包括第一旋转矩阵生成子模块、第二旋转矩阵生成子模块，第三旋转矩阵生成子模块、判断子模块、磁矢量生成子模块：

所述第一旋转矩阵生成子模块，用于根据所述随机扫描序列，生成每种组织的射频脉冲旋转矩阵R_RF；

所述第二旋转矩阵生成子模块，用于根据所述随机扫描序列，生成每种组织的弛豫矩阵R_relax；

所述判断子模块，用于判断所述多种组织中是否具有偏共振特性df(Hz)的组织；

所述第三旋转矩阵生成子模块，用于根据所述随机扫描序列，生成每种具有偏共振特性的组织的偏共振旋转矩阵R_off；

所述磁化矢量生成子模块，用于若所述判断子模块判断出所述多种组织中是否具有偏共振特性df(Hz)的组织，则将所述射频脉冲旋转矩阵R_RF、所述弛豫矩阵R_relax、所述偏共振旋转矩阵R_off作用于初始磁化矢量，通过迭代计算，生成每种组织的磁化矢量；否则，将所述射频脉冲旋转矩阵R_RF、所述弛豫矩阵R_relax作用于初始磁化矢量，通过迭代计算，生成每种组织的所述磁化矢量。

优选地，所述随机扫描序列包括如下随机扫描序列中的至少一种：翻转角随机扫描序列、重复时间随机扫描序列。

优选地，所述判断子模块进一步用于判断所述随机扫描序列是否包括翻转角随机扫描序列；

所述第一旋转矩阵生成子模块，用于若所述判断模块判断出所述随机扫描序列包括翻转角随机扫描序列，则根据所述翻转角随机扫描序列，生成每种组织的所述射频脉冲旋转矩阵R_RF；若判断模块判断出所述随机扫描序列不包括翻转角随机扫描序列，则根据预设的翻转角非随机扫描序列，生成每种组织的所述射频脉冲旋转矩阵R_RF。

优选地，所述判断子模块进一步用于判断所述随机扫描序列是否包括重复时间随机扫描序列；

所述第二旋转矩阵生成子模块，用于若所述判断模块判断出所述随机扫描序列包括所述重复时间随机扫描序列，则根据所述重复时间随机扫描序列，生成每种组织的一个弛豫矩阵R_relax；若所述判断模块判断出所述随机扫描序列不包括重复时间随机扫描序列，则根据预设的重复时间非随机扫描序列，生成每种组织的一个弛豫矩阵R_relax。

优选地，所述第三旋转矩阵生成子模块用于若所述判断子模块判断出所述随机扫描序列包括所述重复时间随机扫描序列，则根据每种具有偏共振特性的组织的偏共振特性df(Hz)与所述重复时间随机扫描序列，得到每种具有偏共振特性的组织的偏共振旋转矩阵R_off；若所述判断子模块判断出所述随机扫描序列不包括所述重复时间随机扫描序列，则根据每种具有偏共振特性的组织的偏共振特性df(Hz)与所述预设的重复时间非随机扫描序列，得到每种具有偏共振特性的组织的偏共振旋转矩阵R_off。

优选地，所述组织区分度指示获得模块包括：组织磁共振指纹演化矢量计算子模块和相关计算子模块；其中，

所述组织磁共振指纹演化矢量计算子模块用于根据每种组织所对应的磁化矢量，获得每种组织的一组织磁共振指纹演化矢量，所述组织磁共振指纹演化矢量包括所述组织对应的磁化矢量中每个元素分别投影在XOY平面上矢量的模值；

所述相关计算子模块用于将所述多种组织中的任意两种组织的所述磁共振指纹演化矢量进行互相关计算，得到任意两种组织的区分度指示值。

优选地，所述组织区分度指示获得模块进一步用于根据所述多种组织的区分度指示值，生成多种组织的区分度指示矩阵，所述多种组织的区分度指示矩阵中i行j列的元素为第i种组织与第j种组织的区分度指示值。

优选地，所述装置进一步包括组织区分度确定模块，用于根据所述区分度指示值，确定利用所述随机扫描序列生成的磁共振图像中所述多种组织分别对应的亮度曲线的两两之间的区分度。

优选地，所述组织区分度确定模块用于判断所述多种组织两两之间的区分度指示值是否小于设定的阈值；若所述任意两种组织的区分度指示值小于设定的阈值，则确定在一幅图像中，所述任意两种组织分别对应的两条亮度曲线的区分度高。

从上述方案中可以看出，由于对多种组织中的每种组织，根据随机扫描序列，生成对应所述每种组织的一组磁化矢量；根据所述多种组织中每种组织所对应的磁化矢量，计算得到所述多种组织两两之间的区分度指示值。这样，医生可以根据本发明提供的方法，选择适合的随机扫描序列，生成包括多种组织对应的多条亮度曲线的磁共振图像，使得在该磁共振图像中多种组织分别对应的亮度曲线的趋势相差较大。这样，医生就能够在一幅图像中，准确地分辨出多种组织的特性，能够在该图像中同时观察、对比、分析多种组织。

附图说明

下面将通过参照附图详细描述本发明的优选实施例，使本领域的普通技术人员更清楚本发明的上述及其它特征和优点，附图中：

图1为根据本发明实施例提供的一种指示组织间区分度的方法的一流程图；

图2为根据本发明实施例提供的一种指示组织间区分度的方法的又一流程图；

图3为根据本发明实施例提供的计算偏共振旋转矩阵作用后的磁化矢量示意图；

图4为根据本发明实施例提供的另一种计算偏共振旋转矩阵作用后的磁化矢量示意图；

图5为根据本发明实施例提供的多种组织的磁共振图像；

图6为根据本发明实施例提供的一种指示组织间区分度的装置的结构示意图；

图7为根据本发明实施例提供的磁化矢量生成模块的结构示意图；

图8为根据本发明实施例提供的组织区分度指示获得的结构示意图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚，以下举实施例对本发明进一步详细说明。

图1示出了本发明实施例提供的指示组织区分度的方法，该方法可以由计算设备执行。该计算设备可以是服务器、PC机等。该方法过程详述如下：

步骤101：对多种组织中的每种组织，根据随机扫描序列，生成对应所述每种组织的一组磁化矢量；

根据所述随机扫描序列，生成每种组织的射频脉冲旋转矩阵R_RF。根据所述随机扫描序列，生成每种组织的弛豫矩阵R_relax。判断所述多种组织中是否具有偏共振特性df(Hz)的组织。若是，则根据所述随机扫描序列，生成每种具有偏共振特性的组织的偏共振旋转矩阵R_off，将所述射频脉冲旋转矩阵R_RF、所述弛豫矩阵R_relax、所述偏共振旋转矩阵R_off作用于初始磁化矢量，通过迭代计算，生成每种组织的磁化矢量。否则，将所述射频脉冲旋转矩阵R_RF、所述弛豫矩阵R_relax作用于初始磁化矢量，通过迭代计算，生成每种组织的所述磁化矢量。

在一个实施方式中，所述随机扫描序列可包括翻转角随机扫描序列和重复时间随机扫描序列中的至少一种。

其中，判断所述随机扫描序列是否包括所述翻转角随机扫描序列；若所述随机扫描序列包括所述翻转角随机扫描序列，则根据所述翻转角随机扫描序列，生成每种组织的所述射频脉冲旋转矩阵R_RF；若所述随机扫描序列不包括翻转角随机扫描序列，则根据预设的翻转角非随机扫描序列，生成每种组织的所述射频脉冲旋转矩阵R_RF。

其中，判断所述随机扫描序列是否包括重复时间随机扫描序列；若所述随机扫描序列包括重复时间随机扫描序列，则根据所述重复时间随机扫描序列，生成每种组织的一个弛豫矩阵R_relax；若所述随机扫描序列不包括重复时间随机扫描序列，则根据预设的重复时间非随机扫描序列，生成每种组织的一个弛豫矩阵R_relax。

其中，若所述随机扫描序列包括所述重复时间随机扫描序列，则根据每种具有偏共振特性的组织的偏共振特性df(Hz)与所述重复时间随机扫描序列，得到每种具有偏共振特性的组织的偏共振旋转矩阵R_off；若所述随机扫描序列不包括所述重复时间随机扫描序列，则根据每种具有偏共振特性的组织的偏共振特性df(Hz)与所述预设的重复时间非随机扫描序列，得到每种具有偏共振特性的组织的偏共振旋转矩阵R_off。

步骤102：根据所述多种组织中每种组织所对应的磁化矢量，计算得到所述多种组织两两之间的区分度指示值。

其中，针对所述多种组织中的每一种组织，可根据所述组织所对应的磁化矢量，获得所述组织的一个组织MRF演化矢量，所述组织MRF演化矢量包括所述组织对应的磁化矢量中每个元素分别投影在XOY平面上矢量的模值；之后可将所述多种组织中的任意两种组织的组织MRF演化矢量进行互相关计算，得到任意两种组织的区分度指示值。

进一步地，可根据所述多种组织的区分度指示值，生成多种组织的区分度指示矩阵，所述多种组织的区分度指示矩阵中i行j列的元素为第i种组织与第j种组织的区分度指示值。

进一步地，根据所述区分度指示值，确定利用所述随机扫描序列生成的磁共振图像中所述多种组织分别对应的亮度曲线的两两之间的区分度。

具体地，可判断所述多种组织区分度指示值是否小于一设定的阈值；若多种组织区分度指示值小于该设定的阈值，则可确定在一幅图像中，所述多种组织分别对应的多条亮度曲线的区分度高；否则，可确定在一幅图像中，所述多种组织分别对应的多条亮度曲线的区分度低。

下面结合一个具体示例，对本发明实施方式中的指示组织间区分度的方法进行详细描述。图2为根据本发明实施方式的又一指示组织间区分度的方法的流程示意图。该方法可以由计算设备执行。该计算设备可以是服务器、PC机等。

如图2所示，该方法可具体包括如下步骤：

步骤201：判断随机扫描序列是否包括翻转角随机扫描序列；

该随机扫描序列包括如下随机扫描序列中的至少一种：翻转角随机扫描序列、重复时间随机扫描序列。

若该随机扫描序列包括翻转角随机扫描序列，则执行步骤202。否则，执行步骤206。

步骤202：针对多种组织中的每一种组织，根据翻转角(Flip Angle，FA)随机扫描序列FA(i)(i＝1，2，…，N)，生成该组织的射频脉冲旋转矩阵R_RF；

具体的，根据该翻转角随机序列，获得射频脉冲的翻转角度α、相位φ。将该相位φ沿着z轴进行旋转，得到一个旋转矩阵Rot_Z(φ)。将翻转角度α沿着x轴进行旋转，得到一个旋转矩阵Rot_X(α)。并且，将-φ沿着z轴进行旋转，得到一个旋转矩阵Rot_Z(-φ)。根据Rot_Z(φ)、Rot_X(α)、Rot_Z(-φ)，得到该组织对应的射频脉冲旋转矩阵R_RF。将该射频脉冲旋转矩阵R_RF作用于初始磁化矢量，通过迭代方法，生成该组织的R_RF作用后的一组磁化矢量。

在一个具体的实施中，如图3所示，可通过下式(1)，根据翻转角随机扫描序列FA(i)，获得该射频脉冲的翻转角度α：

α=abs(FA(i))，i＝1，2，…，N；   (1)

可通过如下式(2)，根据翻转角随机扫描序列FA(i)，获得相位φ：

φ=angLe(FA(i))，i＝1，2，…，N；   (2)

可通过如下式(3)，将相位φ沿着z轴进行旋转，得到旋转矩阵Rot_Z(φ)：

${\mathrm{Rot}}_Z\left(\mathrm{\φ}\right)=\left(\begin{array}{ccc}\cos \left(\mathrm{\φ}\right)& \sin \left(\mathrm{\φ}\right)& 0\\ -\sin \left(\mathrm{\φ}\right)& \cos \left(\mathrm{\φ}\right)& 0\\ 0& 0& 1\end{array}\right)---\left(3\right)$

可通过如下式(4)，将射频脉冲的翻转角度α沿着x轴进行旋转，得到旋转矩阵Rot_X(α)；

${\mathrm{Rot}}_X\left(\mathrm{\α}\right)=\left(\begin{array}{ccc}1& 0& 0\\ 0& \cos \left(\mathrm{\α}\right)& \sin \left(\mathrm{\α}\right)\\ 0& -\sin \left(\mathrm{\α}\right)& \cos \left(\mathrm{\α}\right)\end{array}\right)---\left(4\right)$

可通过如下式(5)，将-φ沿着z轴进行旋转，得到旋转矩阵Rot_Z(-φ)：

${\mathrm{Rot}}_Z(-\mathrm{\φ})=\left(\begin{array}{ccc}\cos (-\mathrm{\φ})& \sin (-\mathrm{\φ})& 0\\ -\sin (-\mathrm{\φ})& \cos (-\mathrm{\φ})& 0\\ 0& 0& 1\end{array}\right)---\left(5\right)$

可通过如下式(6)，根据Rot_Z(φ)、Rot_X(α)、Rot_Z(-φ)，得到射频脉冲旋转矩阵R_RF：

R_RF＝Rot_Z(φ)Rot_X(α)Rot_Z(-φ)   (6)

在一种具体实时方式中，如图4所示，场不均匀相位β也可以用于射频脉冲旋转矩阵R_RF的计算。具体的，根据频率偏移△ω、磁旋比γ、射频场B1，分别计算有效射频脉冲的翻转角度α′与场不均匀相位β。根据该翻转角随机序列，计算获得射频脉冲的相位φ。将射场不均匀相位β沿着y轴进行旋转，得到旋转矩阵Rot_Y(β)。将射场不均匀相位-β沿着y轴进行旋转，得到旋转矩阵Rot_Y(-β)。有效射频脉冲的翻转角度α′沿着x轴进行旋转，得到旋转矩阵Rot_X(α′)。将相位φ沿着z轴进行旋转，得到旋转矩阵Rot_Z(φ)。可通过上述式(5)，将-φ沿着z轴进行旋转，得到旋转矩阵Rot_Z(-φ)。根据Rot_Z(φ)、Rot_X(α′)、Rot_Z(-φ)、Rot_Y(β)、Rot_Y(-β)，得到射频脉冲旋转矩阵R_RF。

在一个具体实施中，可通过如下式(7)，根据频率偏移△ω、磁旋比γ、射频场B1，计算有效射频脉冲的翻转角度α′：

${\mathrm{\α}}^{\′}\underset{t_0}{\overset{t_0+\mathrm{\τ}}{\∫}}\sqrt{{\left(\mathrm{\Δ\ω}\right)}^2+{\left({\mathrm{\γB}}_1\right)}^2}\mathrm{dt}---\left(7\right)$

其中，特定空间位置、时间t上的频率偏移△ω可以通过下式(8)通过磁场强度B，设定的背景场强度B0，磁旋比γ，z方向的位置矢量，计算得到的：

$\mathrm{\Δ\ω}(\stackrel{\→}{r},t)=\mathrm{\γ}(B_0-B(\stackrel{\→}{r},t))\stackrel{\→}{z}---\left(8\right)$

可通过上式(2)，根据翻转角随机扫描序列FA(i)，获得相位φ。

可通过如下(9)，计算场不均匀相位β：

$\mathrm{\β}=\arctan \left(\frac{\mathrm{\Δ\ω}}{{\mathrm{\γB}}_1}\right)---\left(9\right)$

可通过如下式(10)，将射场不均匀相位β沿着y轴进行旋转，得到旋转矩阵Rot_Y(β)；

${\mathrm{Rot}}_Y\left(\mathrm{\β}\right)=\left(\begin{array}{ccc}\cos \left(\mathrm{\β}\right)& 0& \sin \left(\mathrm{\β}\right)\\ 0& 1& 0\\ -\sin \left(\mathrm{\β}\right)& 0& \cos \left(\mathrm{\β}\right)\end{array}\right)---\left(10\right)$

可通过如下式(11)，将射场不均匀相位-β沿着y轴进行旋转，得到旋转矩阵Rot_Y(-β)；

${\mathrm{Rot}}_Y\left(\mathrm{\β}\right)=\left(\begin{array}{ccc}\cos (-\mathrm{\β})& 0& \sin (-\mathrm{\β})\\ 0& 1& 0\\ -\sin (-\mathrm{\β})& 0& \cos (-\mathrm{\β})\end{array}\right)---\left(11\right)$

可通过如下式(12)，将有效射频脉冲的翻转角度α′沿着x轴进行旋转，得到旋转矩阵Rot_X(α′)；

${\mathrm{Rot}}_X\left({\mathrm{\α}}^{\′}\right)=\left(\begin{array}{ccc}1& 0& 0\\ 0& \cos \left({\mathrm{\α}}^{\′}\right)& \sin \left({\mathrm{\α}}^{\′}\right)\\ 0& -\sin \left({\mathrm{\α}}^{\′}\right)& \cos \left({\mathrm{\α}}^{\′}\right)\end{array}\right)---\left(12\right)$

可通过上述式(3)，将相位φ沿着z轴进行旋转，得到旋转矩阵Rot_Z(φ)。可通过上述式(5)，将-φ沿着z轴进行旋转，得到旋转矩阵Rot_Z(-φ)。

可通过如下式(13)，根据Rot_Z(φ)、Rot_X(α′)、Rot_Z(-φ)、Rot_Y(β)、Rot_Y(-β)，得到射频脉冲旋转矩阵R_RF：

R_RF=Rot_Z(φ)Rot_Y(β)Rot_X(α′)Rot_Y(-β)Rot_Z(-φ)   (13)

步骤203：判断随机扫描序列是否包括重复时间随机扫描序列；

若该随机扫描序列包括重复时间随机扫描序列，则执行步骤204。否则，执行步骤207。

步骤204：针对多种组织中的每一种组织，根据重复时间(Time ofRepeat，TR)随机扫描序列TR(i)(i＝1，2，…，N)，生成该组织的弛豫矩阵R_relax；

具体的，可根据重复时间随机序列、组织的纵向弛豫时间常数T1，组织的横向弛豫时间常数T2，获得该组织的弛豫矩阵R_relax。将该该组织的弛豫矩阵R_relax作用于磁化矢量m^-，可得到该组织的弛豫矩阵R_relax作用后的磁化矢量m⁺。其中，T1、T2反映了组织的特性。

在一个具体的实现中，可通过如下式(14)(15)，根据重复时间TR，获得弛豫矩阵R_relax：

△t=TR(i)i＝1，2，…，N   (14)

$R_{\mathrm{relax}}=\left(\begin{array}{ccc}{e}^{-\frac{\mathrm{\Δt}}{T_2}}& 0& 0\\ 0& e^{-\frac{\mathrm{\Δt}}{T_2}}& 0\\ 0& 0& 1-e^{-\frac{\mathrm{\Δt}}{T_1}}\end{array}\right)---\left(15\right)$

其中，T1是纵向弛豫时间常数，T2是横向弛豫时间常数。

步骤205：判断所述多种组织中是否包括具有偏共振特性df(Hz)的组织。若所述多种组织中包括具有偏共振特性df(Hz)的组织，则针对多种组织中具有偏共振特性df(Hz)的每一种组织，根据该组织的偏共振特性df(Hz)与重复时间随机扫描序列，声称该组织的偏共振旋转矩阵R_off；

具体的，根据重复时间随机序列TR(i)与该组织的偏共振特性df(Hz)，可生成偏共振相位角(单位rad)。根据该偏共振相位角可生成偏共振旋转矩阵R_off。将该偏共振旋转矩阵R_off作用于磁化矢量m^-，可得到R_off作用后的磁化矢量m⁺。

在一个具体实现中，可通过如下式(16)，根据重复时间随机序列TR(i)与该组织的偏共振特性df(Hz)，生成偏共振相位角

通过如下式(17)，根据该偏共振相位角生成偏共振旋转矩阵R_off：

执行步骤205后，跳转到执行步骤209。

步骤206：针对多种组织中的每一种组织，根据预设的翻转角非随机扫描序列，生成该组织的射频脉冲旋转矩阵R_RF；

在一种实施方式中，该翻转角非随机扫描序列可以是翻转角正弦序列，或翻转角余弦序列等。

根据翻转角非随机扫描序列，获得该组织的射频脉冲旋转矩阵R_RF作用后的一组磁化矢量的方法与步骤202中方法相同。其中，频脉冲的翻转角度α、相位φ是根据翻转角非随机扫描序列所获得的。

执行步骤206后，跳转执行步骤204。

步骤207：针对多种组织中的每一种组织，根据预设的重复时间非随机扫描序列，生成该组织弛豫矩阵R_relax；

在一种实施方式中，该重复时间非随机扫描序列可以是重复时间正弦序列，或重复时间余弦序列等。

根据预设的重复时间非随机扫描序列，获得该组织弛豫矩阵R_relax作用后的一组磁化矢量的方法与步骤204中方法相同。其中，该组织的弛豫矩阵R_relax可根据重复时间非随机序列、组织的纵向弛豫时间常数T1，组织的横向弛豫时间常数T2获得。

步骤208：判断所述多种组织中是否包括具有偏共振特性df(Hz)的组织。若所述多种组织中包括具有偏共振特性df(Hz)的组织，则针对多种组织中具有偏共振特性df(Hz)的每一种组织，根据该组织的偏共振特性df(Hz)与重复时间非随机扫描序列，生成该组织偏共振旋转矩阵R_off；

根据该组织的偏共振特性df(Hz)与重复时间非随机扫描序列，可得到该组织偏共振旋转矩阵R_off作用的一组磁化矢量的方法与步骤205中的方法相同。其中，偏共振相位角根据重复时间随机序列TR(i)与该组织的偏共振特性df(Hz)生成。

步骤209：若判断出所述多种组织中具有偏共振特性df(Hz)的组织，则将射频脉冲旋转矩阵R_RF、所述弛豫矩阵R_relax、所述偏共振旋转矩阵R_off作用于初始磁化矢量，通过迭代计算，生成每种组织的磁化矢量；否则，将所述射频脉冲旋转矩阵R_RF、所述弛豫矩阵R_relax作用于初始磁化矢量，通过迭代计算，生成每种组织的所述磁化矢量；

可通过迭代方法，通过如下式(18)，计算得到射频脉冲旋转矩阵R_RF作用后的磁化矢量m⁺：

m⁺=R_RF·m^-   (18)

其中，m^-是本次射频脉冲旋转矩阵R_RF作用之前的磁化矢量，在一个具体的实现中，初始磁化矢量可设置为m₀=[0，0,-1]^T。当场不均匀相位β也可以用于射频脉冲旋转矩阵R_RF的计算时，初始磁化向量可设置为有效射频场如图4所示，通过上述式(8)，计算有效射频场可通过如下式(19)，计算有效射频场

${\stackrel{\→}{B}}_{\mathrm{eff}}={\stackrel{\→}{B}}_1+\left(\frac{\mathrm{\Δ\ω}}{\mathrm{\γ}}\right)\·\stackrel{\→}{z}---\left(19\right)$

其中，B1是射频场，是场偏置。

可通过如下式(20)，将该组织的弛豫矩阵R_relax作用于磁化矢量m^-，得到组织的弛豫矩阵R_relax作用后的磁化矢量m⁺：

m⁺＝R_relax·m^-   (20)

其中，m^-为R_relax作用之前的磁化矢量，在一个具体的实现中，初始磁化矢量可以设置为m₀=[0，0,-1]^T。

通过如下式(21)，将该组织的偏共振旋转矩阵R_off作用于磁化矢量m^-，可得到该组织的偏共振旋转矩阵R_off作用后的磁化矢量m⁺：

m⁺=R_off·m^-   (21)。

其中，m^-为R_off作用之前的磁化矢量，在一个具体的实现中，若在计算初始磁化矢量可以设置为m₀=[0，0,-1]^T。

步骤210：根据所述多种组织中每种组织所对应的磁化矢量，计算得到所述多种组织的区分度指示值；

其中，多组组织的区分度指示值能够指示在一幅图像中所述多种组织分别对应的多条亮度曲线的区分度。

在一种具体实施中，可以通过如下方法，获得多种组织的区分度指示值：

针对所述多种组织中的每一种组织，根据该组织所对应的磁化矢量，获得该组织的一个组织MRF演化矢量。该组织MRF演化矢量包括所述组织对应的磁化矢量中每个元素分别投影在XOY平面上矢量的模值。将多种组织中的任意两种组织的组织MRF演化矢量进行互相关计算，得到任意两种组织的区分度指示值。

例如，通过如下方法，可得到组织Ti的组织MRF演化矢量：根据如上步骤201—203，以初始磁化矢量为初始值，分别进行N次迭代运算，得到组织Ti的N个磁化矢量；将该N个磁化矢量在XOY平面上投影的向量作为该组织MRF演化矢量的N个向量。例如，可将初始磁化矢量m₀作为初始值，根据如上步骤201-203迭代方法，可依次获得组织Ti的磁化矢量(m₁～m_N)。可分别计算该磁化矢量(m₁～m_N)在XOY平面上投影向量的模值S(1)～S(N)，得S_Ti＝[|S(1)|，|S(2)|，…，|S(N)|]。

在一种具体实施中，可通过如下式(22)，将组织Ti、Tj的组织MRF演化矢量进行互相关，得到组织Ti、Tj区分度指示值R_Ti，TJ：

$R_{\mathrm{Ti},\mathrm{Tj}}=\frac{E\left\{\right[S_{\mathrm{Ti}}-{\mathrm{\μ}}_{\mathrm{Ti}}]\·[S_{\mathrm{Tj}}-{\mathrm{\μ}}_{\mathrm{Tj}}\left]\right\}}{\sqrt{E\left\{{\left|S_{\mathrm{Ti}}\right|}^2\right\}\·E\left\{{\left|S_{\mathrm{Tj}}\right|}^2\right\}}}---\left(22\right)$

其中，S_Ti为组织Ti的组织向量；

S_Tj为组织Tj的组织向量；

μ_Ti为组织Ti的亮度均值；μ_Tj为组织Tj的亮度均值

在另一种具体实施中，还可以通过关联度评估，模式识别，机器学习等方法，获得多种组织的区分度指示值。

进一步地，还可以根据所述多种组织的区分度指示值，生成多种组织的区分度指示矩阵。该多种组织的区分度指示矩阵中i行j列的元素为第i种组织与第j种组织的区分度指示值。具体的，该矩阵可以表示为：

其中，该矩阵具有M列、M行，分别代表M种组织。该矩阵中位于i行j列的元素C_ij为两种组织Ti、Tj区分度指示值。

步骤211：根据所述多种组织的区分度指示值，确定利用所述随机扫描序列生成的磁共振图像中所述多种组织分别对应的亮度曲线的区分度。

具体的，可判断获得的多种组织区分度指示值是否小于设定的阈值。若多种组织区分度指示值是否小于设定的阈值，则在一幅图像中，所述多种组织分别对应的多条亮度曲线的区分度较高。

在一种具体实施中，可从该矩阵中分别得到该多种组织中任意两种组织区分度指示值。若获得的任意两种组织区分度指示值均小于设定的阈值，则使用该翻转角随机扫描序列与重复时间随机扫描序列能够得到对该多组织区分度较高的图像，即则在一幅图像中，所述多种组织分别对应的多条亮度曲线的趋势相差较大。

在另一种具体实施中，分别得到该多种组织中任意两种组织区分度指示值。该指示值可以通过两种组织的组织MRF演化矢量进行互相关的取反运算获得。若获得的任意两种组织区分度指示值均大于设定的阈值，则使用该翻转角随机扫描序列与重复时间随机扫描序列能够得到对该多组织区分度较高的图像，即则在一幅图像中，所述多种组织分别对应的多条亮度曲线的趋势相差较大。

示范性地，当翻转角随机扫描序列重复时间随机扫描序列TR_t=U(4)+10.0时，可以通过如上方法，确定利用上述随机扫描序列生成的磁共振图像中多种组织(CSF、脂肪Fat、白质White Matter、Gray Matter、白质-30Hz WhiteMatter-30Hz)分别对应的亮度曲线的区分度。其中，White Matter-30Hz指具有-30Hz偏振特性的白质。其每种组织特性参数如下表所示：

表1：组织特性参数

通过图2所示的方法，得到任意两种组织的区分度指示值如下表所示：

表2：组织区分度指示值

如表2所示，任意两种组织交叉位置所示为该两种组织区分度指示值。例如，WhiteMatter与Fat区分度指示值为0.9226，CSF与Fat的区分度指示值为0.4579。

根据表2所示的组织区分度指示值，若将阈值设置为0.7，则由于White Matter与Fat区分度指示值、CSF与Fat区分度指示值、WM(-30Hz)与CSF区分度指示值都小于设定的阈值，则在一幅图像中，White Matter、Fat、CSF、WM(-30Hz)分别对应的多条亮度曲线的区分度较高。

图5为使用本实施例中的翻转角随机扫描序列和重复时间随机扫描序列生成的磁共振图像。该图像中包括CSF、Fat、White Matters、Gray Matters、WM(-30Hz)分别对应的亮度曲线。由图5可以看出，组织区分度指示值小于设定阈值的组织所对应的亮度曲线趋势相差较大，例如CSF与Fat。组织区分度指示值于设定阈值的组织所对应的亮度曲线趋势相差较小，例如White Matters与Fat。

从上述方案中可以看出，由于本发明对多种组织中的每种组织，根据随机扫描序列，生成对应所述每种组织的一组磁化矢量；根据所述多种组织中每种组织所对应的磁化矢量，计算得到所述多种组织两两之间的区分度指示值。这样，医生可以根据本发明提供的方法，选择适合的随机扫描序列，生成包括多种组织对应的多条亮度曲线的磁共振图像，使得在该磁共振图像中多种组织分别对应的亮度曲线的趋势相差较大。这样，医生就能够在一幅图像中，准确地分辨出多种组织的特性，能够在该图像中同时观察、对比、分析多种组织。

本发明实施方式还提出了一种指示组织区分度的装置。图3为根据本发明实施例提供的一种指示组织区分度装置的结构示意图。如图3所示，该装置包括磁化矢量生成模块301、组织区分度指示获得模块302。

磁化矢量生成模块301，用于对多种组织中的每种组织，根据随机扫描序列，生成对应每种组织的一组磁化矢量；

组织区分度指示获得模块302，用于根据所述多种组织中每种组织所对应的磁化矢量，计算得到所述多种组织两两之间的区分度指示值。

随机扫描序列包括如下随机扫描序列中的至少一种：翻转角随机扫描序列、重复时间随机扫描序列。

具体的，如图7所示，磁化矢量生成模块包括第一旋转矩阵生成子模块3011、第二旋转矩阵生成子模块3012，第三旋转矩阵生成子模块3013、判断子模块3014、磁矢量生成子模块3015。

第一旋转矩阵生成子模块3011，用于根据所述随机扫描序列，生成每种组织的射频脉冲旋转矩阵R_RF。第二旋转矩阵生成子模块3012，用于根据所述随机扫描序列，生成每种组织的弛豫矩阵R_relax。判断子模块3014，用于判断所述多种组织中是否具有偏共振特性df(Hz)的组织。第三旋转矩阵生成子模块3013用于根据所述随机扫描序列，生成每种具有偏共振特性的组织的偏共振旋转矩阵R_off。磁化矢量生成子模块3015，用于若所述判断子模块3014判断出所述多种组织中是否具有偏共振特性df(Hz)的组织，则将所述射频脉冲旋转矩阵R_RF、所述弛豫矩阵R_relax、所述偏共振旋转矩阵R_off作用于初始磁化矢量，通过迭代计算，生成每种组织的磁化矢量；否则，将所述射频脉冲旋转矩阵R_RF、所述弛豫矩阵R_relax作用于初始磁化矢量，通过迭代计算，生成每种组织的所述磁化矢量。

判断子模块3014进一步用于判断所述随机扫描序列是否包括翻转角随机扫描序列。第一旋转矩阵生成子模块3011，用于若所述判断模块判断出所述随机扫描序列包括翻转角随机扫描序列，则根据所述翻转角随机扫描序列，生成每种组织的所述射频脉冲旋转矩阵R_RF；若判断模块3014判断出所述随机扫描序列不包括翻转角随机扫描序列，则根据预设的翻转角非随机扫描序列，生成每种组织的所述射频脉冲旋转矩阵R_RF。

判断子模块3014进一步用于判断所述随机扫描序列是否包括重复时间随机扫描序列。第二旋转矩阵生成子模块3012，用于若所述判断模块3014判断出所述随机扫描序列包括所述重复时间随机扫描序列，则根据所述重复时间随机扫描序列，生成每种组织的一个弛豫矩阵R_relax；若所述判断模块3014判断出所述随机扫描序列不包括重复时间随机扫描序列，则根据预设的重复时间非随机扫描序列，生成每种组织的一个弛豫矩阵R_Relax。

第三旋转矩阵生成子模块3013用于若所述判断子模块3014判断出所述随机扫描序列包括所述重复时间随机扫描序列，则根据每种具有偏共振特性的组织的偏共振特性df(Hz)与所述重复时间随机扫描序列，得到每种具有偏共振特性的组织的偏共振旋转矩阵R_off；若所述判断子模块3014判断出所述随机扫描序列不包括所述重复时间随机扫描序列，则根据每种具有偏共振特性的组织的偏共振特性df(Hz)与所述预设的重复时间非随机扫描序列，得到每种具有偏共振特性的组织的偏共振旋转矩阵R_off。

这里，磁化矢量生成模块301所使用的方法与步骤201-209中的方法一致，这里不再赘述。

具体的，如图8所示，所述组织区分度指示获得模块302包括：组织MRF演化矢量计算子模块3021和相关计算子模块3022。其中，组织MRF演化矢量计算子模块3021用于针对所述多种组织中的每一种组织，根据所述组织所对应的磁化矢量，获得所述组织的一个组织MRF演化矢量，组织MRF演化矢量包括所述组织对应的磁化矢量中每个元素分别投影在XOY平面上矢量的模值。相关计算子模块3022用于将所述多种组织中的任意两种组织的组织MRF演化矢量进行互相关计算，得到任意两种组织的区分度指示值。

所述组织区分度指示获得模块302进一步用于根据所述多种组织的区分度指示值，生成多种组织的区分度指示矩阵，所述多种组织的区分度指示矩阵中i行j列的元素为第i种组织与第j种组织的区分度指示值。

这里，组织区分度指示获得模块302所使用的方法与步骤210中的方法一致，这里不再赘述。

该装置进一步还包括组织区分度确定模块，用于根据所述区分度指示值，确定利用所述随机扫描序列生成的磁共振图像中所述多种组织分别对应的亮度曲线的两两之间的区分度。

具体的，所述组织区分度确定模块用于判断所述多种组织区分度指示值是否小于设定的阈值；若多种组织区分度指示值小于设定的阈值，则确定在一幅图像中，所述多种组织分别对应的多条亮度曲线的区分度高。

这里，组织区分度确定模块所使用的方法与步骤211中的方法一致，这里不再赘述。

从上述装置中可以看出，由于对多种组织中的每种组织，根据随机扫描序列，生成对应所述每种组织的一组磁化矢量；根据所述多种组织中每种组织所对应的磁化矢量，计算得到所述多种组织两两之间的区分度指示值。这样，医生可以根据本发明提供的装置，选择适合的随机扫描序列，生成包括多种组织对应的多条亮度曲线的磁共振图像，使得在该磁共振图像中多种组织分别对应的亮度曲线的趋势相差较大。这样，医生就能够在一幅图像中，准确地分辨出多种组织的特性，能够在该图像中同时观察、对比、分析多种组织。

本发明还提供了一种机器可读的存储介质，存储用于使一机器执行如本文所述的一种指示组织间区分度的方法的指令。具体地，可以提供配有存储介质的系统或者装置，在该存储介质上存储着实现上述实施例中任一实施例的功能的软件程序代码，且使该系统或者装置的计算机(或CPU或MPU)读出并执行存储在存储介质中的程序代码。

在这种情况下，从存储介质读取的程序代码本身可实现上述实施例中任何一项实施例的功能，因此程序代码和存储程序代码的存储介质构成了本发明的一部分。

用于提供程序代码的存储介质实施例包括软盘、硬盘、磁光盘、光盘(如CD—ROM、CD-R、CD—RW、DVD—ROM、DVD—RAM、DVD—RW、DVD+RW)、磁带、非易失性存储卡和ROM。可选择地，可以由通信网络从服务器计算机上下载程序代码。

此外，应该清楚的是，不仅可以通过执行计算机所读出的程序代码，而且可以通过基于程序代码的指令使计算机上操作的操作系统等来完成部分或者全部的实际操作，从而实现上述实施例中任意一项实施例的功能。

此外，可以理解的是，将由存储介质读出的程序代码写到插入计算机内的扩展板中所设置的存储器中或者写到与计算机相连接的扩展单元中设置的存储器中，随后基于程序代码的指令使安装在扩展板或者扩展单元上的CPU等来执行部分和全部实际操作，从而实现上述实施例中任一实施例的功能。

以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (20)

1.一种指示组织间区分度的方法，包括：

对多种组织中的每种组织，根据随机扫描序列，生成对应所述每种组织的一组磁化矢量；

根据所述多种组织中每种组织所对应的磁化矢量，计算得到所述多种组织两两之间的区分度指示值。

2.根据权利要求1所述的方法，其特征在于，所述对多种组织中的每种组织，根据随机扫描序列，生成对应所述组织的一组磁化矢量，包括：

根据所述随机扫描序列，生成每种组织的射频脉冲旋转矩阵R_RF；

根据所述随机扫描序列，生成每种组织的弛豫矩阵R_relax；

判断所述多种组织中是否具有偏共振特性df(Hz)的组织；

若是，则根据所述随机扫描序列，生成每种具有偏共振特性的组织的偏共振旋转矩阵R_off，将所述射频脉冲旋转矩阵R_RF、所述弛豫矩阵R_relax、所述偏共振旋转矩阵R_off作用于初始磁化矢量，通过迭代计算，生成每种组织的磁化矢量；

否则，将所述射频脉冲旋转矩阵R_RF、所述弛豫矩阵R_relax作用于初始磁化矢量，通过迭代计算，生成每种组织的所述磁化矢量。

3.根据权利要求2所述的方法，其特征在于，所述随机扫描序列包括如下随机扫描序列中的至少一种：翻转角随机扫描序列和重复时间随机扫描序列。

4.根据权利要求3所述的方法，其特征在于，根据所述随机扫描序列，生成所述组织的射频脉冲旋转矩阵R_RF，包括：

判断所述随机扫描序列是否包括所述翻转角随机扫描序列；

若所述随机扫描序列包括所述翻转角随机扫描序列，则根据所述翻转角随机扫描序列，生成每种组织的所述射频脉冲旋转矩阵R_RF；

若所述随机扫描序列不包括翻转角随机扫描序列，则根据预设的翻转角非随机扫描序列，生成每种组织的所述射频脉冲旋转矩阵R_RF。

5.根据权利要求3或4所述的方法，其特征在于，根据所述随机扫描序列，生成所述组织的弛豫矩阵R_relax，包括：

判断所述随机扫描序列是否包括重复时间随机扫描序列；

若所述随机扫描序列包括重复时间随机扫描序列，则根据所述重复时间随机扫描序列，生成每种组织的一个弛豫矩阵R_relax；

若所述随机扫描序列不包括重复时间随机扫描序列，则根据预设的重复时间非随机扫描序列，生成每种组织的一个弛豫矩阵R_relax。

6.根据权利要求5所述的方法，其特征在于，所述根据所述随机扫描序列，生成每种具有偏共振特性的组织的偏共振旋转矩阵R_off，包括：

若所述随机扫描序列包括所述重复时间随机扫描序列，则根据每种具有偏共振特性的组织的偏共振特性df(Hz)与所述重复时间随机扫描序列，得到每种具有偏共振特性的组织的偏共振旋转矩阵R_off；

若所述随机扫描序列不包括所述重复时间随机扫描序列，则根据每种具有偏共振特性的组织的偏共振特性df(Hz)与所述预设的重复时间非随机扫描序列，得到每种具有偏共振特性的组织的偏共振旋转矩阵R_off。

7.根据权利要求1所述的方法，其特征在于，所述根据所述多种组织中每种组织所对应的磁化矢量，计算得到所述多种组织的区分度指示值，包括：

根据每种组织所对应的磁化矢量，获得所述组织的一组织磁共振指纹演化矢量，每种组织磁共振指纹演化矢量包括所述每种组织对应的磁化矢量中每个元素分别投影在XOY平面上矢量的模值；

将所述多种组织中的任意两种组织的所述组织磁共振指纹演化矢量进行互相关计算，得到任意两种组织的区分度指示值。

8.根据权利要求1所述的方法，其特征在于，进一步包括：

根据所述多种组织的区分度指示值，生成多种组织的区分度指示矩阵，所述多种组织的区分度指示矩阵中i行j列的元素为第i种组织与第j种组织的区分度指示值。

9.根据权利要求7所述的方法，其特征在于，还包括，根据所述区分度指示值，确定利用所述随机扫描序列生成的磁共振图像中所述多种组织分别对应的亮度曲线的两两之间的区分度。

10.根据权利要求9所述的方法，其特征在于，所述根据所述区分度指示值，确定利用所述随机扫描序列生成的磁共振图像中所述多种组织分别对应的亮度曲线的两两之间的区分度，包括：

判断所述多种组织中的任意两种组织的区分度指示值是否小于设定的阈值；

若所述任意两种组织的区分度指示值小于设定的阈值，则确定在一幅图像中，所述任意两种组织分别对应的两条亮度曲线的区分度高。

11.一种指示组织间区分度的装置，包括：

磁化矢量生成模块，用于对多种组织中的每种组织，根据随机扫描序列，生成对应每种组织的一组磁化矢量；

组织区分度指示获得模块，用于根据所述多种组织中每种组织所对应的磁化矢量，计算得到所述多种组织两两之间的区分度指示值。

12.根据权利要求11所述的装置，其特征在于，磁化矢量生成模块包括第一旋转矩阵生成子模块、第二旋转矩阵生成子模块，第三旋转矩阵生成子模块、判断子模块、磁矢量生成子模块：

所述第一旋转矩阵生成子模块，用于根据所述随机扫描序列，生成每种组织的射频脉冲旋转矩阵R_RF；

所述第二旋转矩阵生成子模块，用于根据所述随机扫描序列，生成每种组织的弛豫矩阵R_relax；

所述判断子模块，用于判断所述多种组织中是否具有偏共振特性df(Hz)的组织；

所述第三旋转矩阵生成子模块，用于根据所述随机扫描序列，生成每种具有偏共振特性的组织的偏共振旋转矩阵R_off；

所述磁化矢量生成子模块，用于若所述判断子模块判断出所述多种组织中是否具有偏共振特性df(Hz)的组织，则将所述射频脉冲旋转矩阵R_RF、所述弛豫矩阵R_relax、所述偏共振旋转矩阵R_off作用于初始磁化矢量，通过迭代计算，生成每种组织的磁化矢量；否则，将所述射频脉冲旋转矩阵R_RF、所述弛豫矩阵R_relax作用于初始磁化矢量，通过迭代计算，生成每种组织的所述磁化矢量。

13.根据权利要求12所述的装置，其特征在于，所述随机扫描序列包括如下随机扫描序列中的至少一种：翻转角随机扫描序列、重复时间随机扫描序列。

14.根据权利要求13所述的装置，其特征在于，所述判断子模块进一步用于判断所述随机扫描序列是否包括翻转角随机扫描序列；

所述第一旋转矩阵生成子模块，用于若所述判断模块判断出所述随机扫描序列包括翻转角随机扫描序列，则根据所述翻转角随机扫描序列，生成每种组织的所述射频脉冲旋转矩阵R_RF；若判断模块判断出所述随机扫描序列不包括翻转角随机扫描序列，则根据预设的翻转角非随机扫描序列，生成每种组织的所述射频脉冲旋转矩阵R_RF。

15.根据权利要求13或14所述的装置，其特征在于，所述判断子模块进一步用于判断所述随机扫描序列是否包括重复时间随机扫描序列；

所述第二旋转矩阵生成子模块，用于若所述判断模块判断出所述随机扫描序列包括所述重复时间随机扫描序列，则根据所述重复时间随机扫描序列，生成每种组织的一个弛豫矩阵R_relax；若所述判断模块判断出所述随机扫描序列不包括重复时间随机扫描序列，则根据预设的重复时间非随机扫描序列，生成每种组织的一个弛豫矩阵R_relax。

16.根据权利要求15所述的装置，其特征在于，所述第三旋转矩阵生成子模块用于若所述判断子模块判断出所述随机扫描序列包括所述重复时间随机扫描序列，则根据每种具有偏共振特性的组织的偏共振特性df(Hz)与所述重复时间随机扫描序列，得到每种具有偏共振特性的组织的偏共振旋转矩阵R_off；若所述判断子模块判断出所述随机扫描序列不包括所述重复时间随机扫描序列，则根据每种具有偏共振特性的组织的偏共振特性df(Hz)与所述预设的重复时间非随机扫描序列，得到每种具有偏共振特性的组织的偏共振旋转矩阵R_off。

17.根据权利要求11所述的装置，其特征在于，所述组织区分度指示获得模块包括：组织磁共振指纹演化矢量计算子模块和相关计算子模块；其中，

所述组织磁共振指纹演化矢量计算子模块用于根据每种组织所对应的磁化矢量，获得每种组织的一组织磁共振指纹演化矢量，所述组织磁共振指纹演化矢量包括所述组织对应的磁化矢量中每个元素分别投影在XOY平面上矢量的模值；

所述相关计算子模块用于将所述多种组织中的任意两种组织的所述磁共振指纹演化矢量进行互相关计算，得到任意两种组织的区分度指示值。

18.根据权利要求11所述的装置，其特征在于，所述组织区分度指示获得模块进一步用于根据所述多种组织的区分度指示值，生成多种组织的区分度指示矩阵，所述多种组织的区分度指示矩阵中i行j列的元素为第i种组织与第j种组织的区分度指示值。

19.根据权利要求17所述的装置，其特征在于，所述装置进一步包括组织区分度确定模块，用于根据所述区分度指示值，确定利用所述随机扫描序列生成的磁共振图像中所述多种组织分别对应的亮度曲线的两两之间的区分度。

20.根据权利要求19所述的装置，其特征在于，所述组织区分度确定模块用于判断所述多种组织两两之间的区分度指示值是否小于设定的阈值；若所述任意两种组织的区分度指示值小于设定的阈值，则确定在一幅图像中，所述任意两种组织分别对应的两条亮度曲线的区分度高。