CN110824398B - 化学交换特征定量方法及设备 - Google Patents

Abstract

本申请适用于生物医学工程领域，提供了一种化学交换特征定量方法及设备，所述方法对于多种CEST对比剂的CEST信号的信号强度实测值，基于各CEST对比剂的理论化学交换特征定量范围，分别计算在与信号强度实测值相同的饱和脉冲参数、共振频率下对应的信号强度理论值，将匹配度最高的理论值对应的理论化学交换特征定量作为CEST对比剂的化学交换特征定量结果。本申请的化学交换特征定量过程无需选择参照信号，因此不易受到其他化学交换效应的影响，由此使得匹配结果对应的理论化学交换特征定量能更准确地代表实际的化学交换特定量，有效地提高了化学交换特定量的准确率。

Description

化学交换特征定量方法及设备

技术领域

本申请涉及生物医学工程领域，具体涉及一种化学交换特定量方法及设备。

背景技术

化学交换饱和转移(Chemical Exchange Saturation Transfer,CEST)成像是一种磁共振分子影像手段，基于内源性或外源性CEST对比剂中可交换质子与水分子的化学交换作用，可无创获取生物体组织能量代谢、酸碱环境、对比剂含量等微观层面信息，在疾病的识别、诊断、评估中具有重要的研究价值和应用潜力。CEST对比剂，包括机体自身代谢物质及外源性对比剂；其中，代表性的内源性CEST对比剂包括蛋白质、肌酸、葡萄糖等，可反映细胞构成、酸碱环境、能量代谢及摄取等重要信息；外源性CEST对比剂主要来源于人工合成的化合物。

在CEST对比剂成像过程中，由于自由水的直接饱和效应和磁化转移(MT)效应相对较强，干扰CEST对比剂的信号，导致CEST对比剂的化学交换特征(如化学交换速率、浓度等)同步定量结果准确率相对较低。针对该问题现有技术中在CEST对比剂成像时常用的方法包括非对称性分析法、三点法、高阶多项式拟合法等。然而，这些方法在计算CEST效应时需要寻找参照信号，而参照信号易包含其他CEST效应，故在计算反映CEST效应强弱的CEST信号时易受到其他CEST效应的影响，从而导致多种CEST对比剂化学交换特征定量准确率较低。

发明内容

本申请的实施例提供一种化学交换特征定量方法及设备，可以解决现有技术中CEST对比剂的化学交换特征定量准确率低的问题。

为了解决上述技术问题，本申请实施例采用的技术方案是：

第一方面本申请实施例提供了一种化学交换特征定量方法，包括：采集目标物体中待探测的多种化学交换饱和转移CEST对比剂各自的CEST信号的信号强度实测值；其中，CEST信号在对CEST对比剂的共振频率上按预设条件施加饱和脉冲时产生；基于每种CEST对比剂对应的理论化学交换特征定量范围，计算每种CEST对比剂各自的CEST信号对应的多个信号强度理论值；其中，每个CEST信号的信号强度理论值在计算时对应的饱和脉冲参数、共振频率与其信号强度实测值对应的饱和脉冲参数、共振频率相同；将每种CEST对比剂各自对应的CEST信号的信号强度实测值与其对应的多个信号强度理论值进行匹配，将匹配度最高的信号强度理论值对应的理论化学交换特征定量作为CEST对比剂的化学交换特征定量结果。

可选地，采集目标物体中待探测的多种化学交换饱和转移CEST对比剂各自的CEST信号的信号强度实测值之前，还包括：

根据各CEST对比剂的化学交换特征，确定各CEST对比剂的有效成像范围；有效成像范围中的参数为饱和脉冲参数；根据所有CEST对比剂的有效成像范围，确定饱和脉冲的饱和脉冲参数范围；其中，饱和脉冲参数范围中的饱和脉冲参数属于各CEST对比剂的有效成像范围。

示例性的，预设条件为：

依次施加至少一个饱和脉冲，每个饱和脉冲参数在设定的饱和脉冲参数范围内随机化。

在第一方面的一种可能实现的方式中，采集目标物体中待探测的多种化学交换饱和转移CEST对比剂各自的CEST信号的信号强度实测值，包括：根据多种CEST对比剂的共振频率，确定待施加的饱和脉冲的频率范围；其中，饱和脉冲的频率范围中的频率涵盖各CEST对比剂的共振频率；从饱和脉冲的频率范围中选择多个频率，多个频率包含各CEST对比剂的共振频率；在多个频率中的每个频率上按预设条件施加饱和脉冲时，采集每个频率对应的CEST信号的信号强度实测值；从每个频率对应的CEST信号的信号强度实测值中，提取每个CEST对比剂的共振频率对应的CEST信号的信号强度实测值作为CEST对比剂的CEST信号的信号强度实测值。

进一步地，基于每种CEST对比剂对应的理论化学交换特征定量范围，计算每种CEST对比剂各自的CEST信号对应的多个信号强度理论值之前，还包括：

量化目标物体的磁化转移效应，得到磁化转移效应的化学交换特征定量结果；基于每种CEST对比剂对应的理论化学交换特征定量范围，计算每种CEST对比剂各自的CEST信号对应的多个信号强度理论值，包括：获取每种CEST对比剂的CEST信号对应的第一生成条件，第一生成条件包括CEST对比剂的共振频率，以及施加饱和脉冲时对应的饱和脉冲参数；确定每种CEST对比剂的理论化学交换特征定量范围；基于理论化学交换特征定量范围中不同的理论化学交换特征定量，利用磁化转移效应的化学交换特征定量结果计算在第一生成条件下生成的CEST信号所对应的多个信号强度理论值。

在第一方面的一种可选的实现方式中，量化目标物体的磁化转移效应，得到磁化转移效应的化学交换特征定量结果，包括：

获取在远离自由水和所有CEST对比剂的共振频率上，按照预设条件施加饱和脉冲时，采集的磁化转移信号的信号强度实测值；基于磁化转移信号的理论化学交换特征定量范围，计算磁化转移信号对应的多个信号强度理论值；其中，在计算磁化转移信号的信号强度理论值时所对应的饱和脉冲参数以及共振频率，与生成磁化转移信号的信号强度实测值时所对应的饱和脉冲参数、共振频率相同；将磁化转移信号的信号强度实测值与其对应的多个信号强度理论值进行匹配，将匹配度最高的信号强度理论值对应的理论化学交换特征定量作为目标物体的磁化转移效应的化学交换特征定量结果。

进一步地，基于磁化转移信号的理论化学交换特征定量范围，计算磁化转移信号对应的多个信号强度理论值，包括：

获取磁化转移信号对应的第二生成条件，第二生成条件包括远离自由水的信息和所有CEST对比剂的共振频率，以及施加饱和脉冲时对应的饱和脉冲参数；确定磁化转移信号的理论化学交换特征定量范围；基于理论化学交换特征定量范围中不同的理论化学交换特征定量，计算在第二生成条件下生成的磁化转移信号所对应的多个信号强度理论值。

可选地，饱和脉冲参数包括：脉冲强度、持续时间、翻转角、占空比中的任一种或一种以上的组合。

第二方面，本申请实施例提供了一种化学交换特征定量装置，包括：

第一采集模块，用于采集目标物体中待探测的多种化学交换饱和转移CEST对比剂各自的CEST信号的信号强度实测值；其中，CEST信号在对CEST对比剂的共振频率上按预设条件施加饱和脉冲时产生；第一计算模块，用于基于每种CEST对比剂对应的理论化学交换特征定量范围，计算每种CEST对比剂各自的CEST信号对应的多个信号强度理论值；其中，每个CEST信号的信号强度理论值在计算时对应的饱和脉冲参数、共振频率与其信号强度实测值对应的饱和脉冲参数、共振频率相同；第一匹配模块，用于将每种CEST对比剂各自对应的CEST信号的信号强度实测值与其对应的多个信号强度理论值进行匹配，将匹配度最高的信号强度理论值对应的理论化学交换特征定量作为CEST对比剂的化学交换特征定量结果。

可选地，化学交换特征定量装置还包括：

第一设定模块，用于根据各CEST对比剂的化学交换特征，确定各CEST对比剂的有效成像范围；有效成像范围中的参数为饱和脉冲参数；第二设定模块，用于根据所有CEST对比剂的有效成像范围，确定饱和脉冲的饱和脉冲参数范围；其中，饱和脉冲参数范围中的饱和脉冲参数属于各CEST对比剂的有效成像范围。

进一步地，预设条件为：依次施加至少一个饱和脉冲，每个饱和脉冲参数在设定的饱和脉冲参数范围内随机化。

可选地，第一采集模块包括：

频率范围设定单元，用于根据多种CEST对比剂的共振频率，确定待施加的饱和脉冲的频率范围；其中，饱和脉冲的频率范围中的频率涵盖各CEST对比剂的共振频率；频率选择单元，用于从饱和脉冲的频率范围中选择多个频率，多个频率包含各CEST对比剂的共振频率；第一采集子单元，用于在多个频率中的每个频率上按预设条件施加饱和脉冲时，采集每个频率对应的CEST信号的信号强度实测值；提取单元，用于从每个频率对应的CEST信号的信号强度实测值中，提取每个CEST对比剂的共振频率对应的CEST信号的信号强度实测值作为CEST对比剂的CEST信号的信号强度实测值。

进一步地，化学交换特征定量装置还包括：

量化模块，用于量化目标物体的磁化转移效应，得到磁化转移效应的化学交换特征定量结果；

第一计算模块包括：

第一获取单元，用于获取每种CEST对比剂的CEST信号对应的第一生成条件，第一生成条件包括CEST对比剂的共振频率，以及施加饱和脉冲时对应的饱和脉冲参数；第一确定单元，用于确定每种CEST对比剂的理论化学交换特征定量范围；第一计算子单元，用于基于理论化学交换特征定量范围中不同的理论化学交换特征定量，利用磁化转移效应的化学交换特征定量结果计算在第一生成条件下生成的CEST信号所对应的多个信号强度理论值。

可选地，量化模块包括：

第二采集模块，用于获取在远离自由水和所有CEST对比剂的共振频率上，按照预设条件施加饱和脉冲时，采集的磁化转移信号的信号强度实测值；第二计算模块，用于基于磁化转移信号的理论化学交换特征定量范围，计算磁化转移信号对应的多个信号强度理论值；其中，在计算磁化转移信号的信号强度理论值时所对应的饱和脉冲参数以及共振频率，与生成磁化转移信号的信号强度实测值时所对应的饱和脉冲参数、共振频率相同；第二匹配模块，用于将磁化转移信号的信号强度实测值与其对应的多个信号强度理论值进行匹配，将匹配度最高的信号强度理论值对应的理论化学交换特征定量作为目标物体的磁化转移效应的化学交换特征定量结果。

进一步地，第二计算模块包括：

第二获取单元，用于获取磁化转移信号对应的第二生成条件，第二生成条件包括远离自由水的信息和所有CEST对比剂的共振频率，以及施加饱和脉冲时对应的饱和脉冲参数；第二确定单元，用于确定磁化转移信号的理论化学交换特征定量范围；第二计算子单元，用于基于理论化学交换特征定量范围中不同的理论化学交换特征定量，计算在第二生成条件下生成的磁化转移信号所对应的多个信号强度理论值。

第三方面，本申请提供了一种化学交换特征定量设备，包括存储器、处理器以及存储在存储器中并可在处理器上运行的计算机程序，处理器执行计算机程序时实现第一方面所述的化学交换特征定量方法。

第四方面，本申请提供了一种计算机可读存储介质，计算机可读存储介质存储有计算机程序，计算机程序被处理器执行时实现第一方面所述的化学交换特征定量方法。

第五方面，本申请提供了一种计算机程序产品，当计算机程序产品在化学交换特征定量设备上运行时，使得化学交换特征定量设备执行第一方面所述的化学交换特征定量方法的步骤。

本申请实施例与现有技术相比存在的有益效果是：

本申请的实施例中针对于每一个CEST对比剂的CEST信号的信号强度实测值，在信号强度理论值进行计算时与信号强度实测值在采集时的采集环境，即饱和脉冲参数、共振频率相同，因此基于理论化学交换特征定量范围，可以计算出多个与信号强度实测值对应的信号强度理论值；信号强度理论值是在与信号强度实测值采集时相同的采集环境下计算得到的多个值，采用一对多的匹配方式使得信号强度实测值能有更准确的匹配结果。本申请的化学交换特征定量过程无需选择参照信号，不易受到其他化学交换效应的影响，由此使得匹配结果对应的理论化学交换特征定量能更准确地代表实际的化学交换特定量，故本申请有效地提高了化学交换特定量的准确率。

可以理解的是，上述第二方面至第五方面的有益效果可以参见上述第一方面中的相关描述，在此不再赘述。

附图说明

为了更清楚地说明本申请实施例中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

图1是本申请一个实施例提供的化学交换特征定量方法流程示意图；

图2是施加饱和脉冲后产生CEST信号的示意图；

图3是本申请的一个实施例提供的S101的流程示意图；

图4是本申请的一个实施例提供的S102的流程示意图；

图5是本申请的一个实施例提供的量化目标物体的磁化转移效应的流程示意图；

图6是本申请的一个实施例提供的计算磁化转移信号对应的多个信号强度理论值的流程示意图；

图7是本申请的一个实施例中计算得到的CEST信号的多个信号强度实测值随饱和脉冲序号的变化曲线；

图8是本申请的一个实施例采集的CEST信号的信号强度实测值随饱和脉冲序号的变化曲线；

图9是图8中的变化曲线在图7中多个变化曲线的匹配结果；

图10是本申请实施例提供的一种化学交换特征定量装置的结构框图；

图11是本申请的一个实施例中第一采集模块的结构框图；

图12是本申请的一个实施例中第一计算模块的结构框图；

图13是本申请的一个实施例中量化模块的结构框图；

图14是本申请的一个实施例提供的化学交换特征定量设备的示意图。

具体实施方式

以下描述中，为了说明而不是为了限定，提出了诸如特定系统结构、技术之类的具体细节，以便透彻理解本申请实施例。然而，本领域的技术人员应当清楚，在没有这些具体细节的其它实施例中也可以实现本申请。在其它情况中，省略对众所周知的系统、装置、电路以及方法的详细说明，以免不必要的细节妨碍本申请的描述。

应当理解，当在本申请说明书和所附权利要求书中使用时，术语“包括”指示所描述特征、整体、步骤、操作、元素和/或组件的存在，但并不排除一个或多个其它特征、整体、步骤、操作、元素、组件和/或其集合的存在或添加。

还应当理解，在本申请说明书和所附权利要求书中使用的术语“和/或”是指相关联列出的项中的一个或多个的任何组合以及所有可能组合，并且包括这些组合。

如在本申请说明书和所附权利要求书中所使用的那样，术语“如果”可以依据上下文被解释为“当...时”或“一旦”或“响应于确定”或“响应于检测到”。类似地，短语“如果确定”或“如果检测到[所描述条件或事件]”可以依据上下文被解释为意指“一旦确定”或“响应于确定”或“一旦检测到[所描述条件或事件]”或“响应于检测到[所描述条件或事件]”。

在本申请说明书中描述的参考“一个实施例”或“一些实施例”等意味着在本申请的一个或多个实施例中包括结合该实施例描述的特定特征、结构或特点。由此，在本说明书中的不同之处出现的语句“在一个实施例中”、“在一些实施例中”、“在其他一些实施例中”、“在另外一些实施例中”等不是必然都参考相同的实施例，而是意味着“一个或多个但不是所有的实施例”，除非是以其他方式另外特别强调。术语“包括”、“包含”、“具有”及它们的变形都意味着“包括但不限于”，除非是以其他方式另外特别强调。

为了说明本申请技术方案，下面通过具体实施例来进行说明。

应理解，下述实施例中各步骤的序号的大小并不意味着执行顺序的先后，各过程的执行顺序应以其功能和内在逻辑确定，而不应对本申请实施例的实施过程构成任何限定。

化学交换饱和转移(Chemical Exchange Saturation Transfer,CEST)定量是一种利用饱和脉冲对待成像的目标物体中CEST对比剂，例如蛋白质、葡萄糖、肌酸等进行预饱和，CEST对比剂会以一定的交换速率和周围的水分子发生化学交换，从而影响自由水的信号强度；通过采用成像序列采集反映CEST效应强弱的CEST信号，即可进行CEST对比剂化学交换特征定量。

本申请的实施例中对目标物体中存在的多种CEST对比剂施加饱和脉冲进行CEST对比剂化学交换特征定量，目标物体是指内部含有CEST对比剂的物体，例如生物体组织。示例性的，需要对生物体的大脑内的CEST对比剂进行化学交换特征定量时，生物体的大脑即为目标物体。

如图1所示，本申请的一个实施例提供了一种化学交换特征定量方法，本实施例中的化学交换特征定量方法的执行主体为化学交换特征定量设备，化学交换特征定量设备可以是服务器，也可以是终端设备，化学交换特征定量方法包括：

S101，采集目标物体中待探测的多种化学交换饱和转移CEST对比剂各自的CEST信号的信号强度实测值；其中，CEST信号在对CEST对比剂的共振频率上按预设条件施加饱和脉冲时产生。

本申请的实施例中，利用饱和脉冲对CEST对比剂进行预饱和以采集CEST信号的信号强度实测值。在每个CEST对比剂对应的共振频率上按预设条件施加饱和脉冲，即可采集到该CEST对比剂的CEST信号的信号强度实测值。其中，CEST对比剂可以是内源性CEST对比剂或外源性CEST对比剂。施加饱和脉冲的预设条件与采集到的CEST信号的强弱存在对应关系，当依次施加多个不同预设条件的饱和脉冲时，可对应采集到多个CEST信号的信号强度实测值。分别针对每一个CEST对比剂在其对应的共振频率上按照预设条件施加饱和脉冲，当施加多个饱和脉冲时，可得到针对于每一个CEST对比剂的CEST信号的多个信号强度实测值。

请一并参阅图2，图2展示出了施加饱和脉冲后产生CEST信号的示意图；在施加饱和脉冲后，采集CEST信号，即可得到CEST信号的信号强度实测值。

在一种可选的实现方式中，施加饱和脉冲进行CEST对比剂的预饱和之前，首先设定好饱和脉冲施加的预设条件。预设条件是饱和脉冲施加所遵循的条件，可以包括饱和脉冲施加的次数和/或饱和脉冲参数等。其中，在CEST定量中有多个重要的饱和脉冲参数，包括脉冲强度、持续时间、翻转角、占空比等；可以根据实际需要选择其中一种或一种以上的组合。例如，本实施例中选择脉冲强度ω₁、持续时间T_s这两个饱和脉冲参数。

可选地，在一个具体的示例中，预设条件为依次施加至少一个饱和脉冲，每个饱和脉冲参数在设定的饱和脉冲参数范围内随机化。

即，在每一个CEST对比剂的共振频率上，依次施加至少一个饱和脉冲，并采集每次施加饱和脉冲之后探测到的CEST信号的信号强度实测值。本实施例中具体的饱和脉冲参数选取利用随机化的方法，即先设定好饱和脉冲参数范围，施加饱和脉冲之前先在饱和脉冲参数范围中随机选取饱和脉冲参数的值，按照选取的饱和脉冲参数的值施加饱和脉冲。在其他的示例中，饱和脉冲参数的值也可以根据实际需求设定为固定值。

例如在一个CEST对比剂的共振频率选择施加n个饱和脉冲，则CEST信号的信号强度实测值的采集过程为：

照预设条件中的饱和脉冲参数范围首先随机产生第一个(当饱和脉冲参数为一种时)或第一组(当饱和脉冲参数为一种以上时)饱和脉冲参数，记录按照该饱和脉冲参数施加第一个饱和脉冲后探测到的CEST信号的信号强度实测值；然后随机产生第二个或第二组饱和脉冲参数并按照该饱和脉冲参数施加饱和脉冲，记录探测到的CEST信号的信号强度实测值；以此类推，在每个饱和脉冲施加之前，先在选定的饱和脉冲参数范围内随机选择饱和脉冲参数，然后再按照该饱和脉冲参数施加饱和脉冲，直至n个饱和脉冲全部施加完，即可得到在该CEST对比剂的共振频率上及设定的饱和脉冲参数范围内的n个CEST信号的信号强度实测值。

S102，基于每种CEST对比剂对应的理论化学交换特征定量范围，计算每种CEST对比剂各自的CEST信号对应的多个信号强度理论值；其中，每个CEST信号的信号强度理论值在计算时对应的饱和脉冲参数、共振频率与其信号强度实测值对应的饱和脉冲参数、共振频率相同。

CEST对比剂的理论化学交换特征定量范围，是指该CEST对比剂可能的化学交换特征定量范围，该范围可通过先验知识确定，例如通过以往的试验数据，或通过查询资料的方式来确定。其中，化学交换特征可以是化学交换速率和/或浓度等；在化学交换特征选择化学交换速率和浓度时，理论化学交换特征定量范围指化学交换速率范围和浓度范围。

基于CEST对比剂对应的理论化学交换特征定量范围，在与CEST信号的信号强度实测值采集时相同的饱和脉冲参数、共振频率条件下，可计算多个与CEST信号对应的信号强度理论值。例如，可利用理论化学交换特征定量范围中不同的理论化学交换特征定量进行计算，针对于一个理论化学交换特征定量，通过理论计算的方式可得到一个信号强度理论值；那么将范围中的多个理论化学交换特征定量分别计算信号强度理论值，即可得到对应于CEST信号的多个信号强度理论值。由于针对一个CEST信号可探测到一个信号强度实测值，因此信号强度理论值与信号强度实测值存在多对一的对应关系。

当针对于每个CEST对比剂施加多个饱和脉冲而采集多个CEST信号时，在计算每个CEST信号的对应的多个信号强度理论值过程中，首先从CEST对比剂的理论化学交换特征定量范围中选取多个理论化学交换特征定量，按照取出的化学交换特征定量，依次计算在每个化学交换特征定量下，每个CEST信号的信号强度理论值。

例如，对于一个CEST对比剂，其共振频率为F1；预设条件中，设定施加n个饱和脉冲，设定的饱和脉冲参数范围是：脉冲强度ω₀～ω_a、持续时间T_s0～T_sb。该CEST对比剂的理论化学交换特征定量范围为：化学交换速率范围K，浓度范围c。在理论化学交换特征定量范围中选择多个理论化学交换特征定量，例如在化学交换速率范围K中等间隔选取p个化学交换速率K₁、K₂…K_p，在浓度范围c中等间隔选取q个浓度c₁、c₂…c_q，则化学交换速率、浓度共可以构成p×q个组合(K₁,c₁)、(K₁,c₂)、(K₁,c₃)…(K_p,c_q)，也即p×q个理论化学交换特征定量。

在采集该CEST对比剂的CEST信号的信号强度实测值时：

在共振频率F1上施加第1至第n个饱和脉冲，随机选择的饱和脉冲参数(脉冲强度,持续时间)分别为：(ω₁，T_s1)、(ω₂，T_s2)…(ω_n，T_sn)，其中，ω₁,ω₂,…ω_n∈(ω₀～ω_a)，T_s1,T_s2,…T_sn∈(T_s0～T_sb)。则采集到的第1个至第n个CEST信号的信号强度实测值分别为CEST_acq1、CEST_acq2…CEST_acqn。

在计算该CEST对比剂的CEST信号的多个信号强度理论值时：

对于第1个理论化学交换特征定量，本例中为第1个化学交换速率和第1个浓度(K₁,c₁)，计算第1～n个CEST信号的信号强度理论值时，分别将CEST对比剂的共振频率F1，以及采集1～n个CEST信号的信号强度实测值时的饱和脉冲参数(ω₁，T_s1)、(ω₂，T_s2)…(ω_n，T_sn)作为计算条件，计算得到第1～n个信号强度理论值CEST¹ _ideal1、CEST¹ _ideal2…CEST¹ _idealn。即，第1个信号强度理论值CEST¹ _ideal1在计算时所用的参量为：[(K₁,c₁)、F1、(ω₁，T_s1)]；第2个信号强度理论值CEST¹ _ideal2在计算时所用的参量为：[(K₁,c₁)、F1、(ω₂，T_s2)]…第n个信号强度理论值CEST¹ _idealn在计算时所用的参量为：[(K₁,c₁)、F1、(ω_n，T_sn)]。

通过以上方法，计算得到了针对于第1个理论化学交换特征定量的n个CEST信号的n个信号强度理论值。

对于第2个理论化学交换特征定量(K₁,c₂)，计算第1～n个CEST信号的信号强度理论值的参量分别为[(K₁,c₂)、F1、(ω₁，T₁)]、[(K₁,c₂)、F1、(ω₂，T₂)]…[(K₁,c₂)、F1、(ω_n，T_n)]，对应计算出1～n个信号强度理论值CEST² _ideal1、CEST² _ideal2…CEST² _idealn。

按照相同的方法，可依次计算每一个理论化学交换特征定量对应第1至第n个CEST信号的信号强度理论值。由于每一个理论化学交换特征定量均可对应计算出一个信号强度理论值，则从理论化学交换特征定量范围中选取多个理论化学交换特征定量后，针对于一个CEST对比剂通过理论计算可得到与其信号强度实测值对应的多个信号强度理论值。

S103，将每种CEST对比剂各自对应的CEST信号的信号强度实测值与其对应的多个信号强度理论值进行匹配，将匹配度最高的信号强度理论值对应的理论化学交换特征定量作为CEST对比剂的化学交换特征定量结果。

可选地，利用点乘最大化原则等算法，将每一个CEST对比剂的CEST信号的信号强度实测值与对应的多个信号强度理论值进行匹配；例如，施加n个饱和脉冲后记录CEST信号的n个信号强度实测值为CEST_acq1、CEST_acq2、CEST_acq3…、CEST_acqn，在对应于这n个信号强度实测值的多个信号强度理论值中，寻找与实测值CEST_acq1、CEST_acq2、CEST_acq3…、CEST_acqn匹配度最高的n个理论值，这n个理论值在计算过程中对应的理论化学交换特征定量即为目标物体的磁化转移效应的化学交换特征定量。

例如，在S102计算时所采用的理论化学交换特征定量中，第m个理论化学交换特征定量(K_i,c_j)∈[(K₁,c₁)、(K₁,c₂)、(K₁,c₃)…(K_q,c_q)]所对应的n个理论值CEST^m _ideal1、CEST^m _ideal2、CEST^m _ideal3…CEST^m _idealn与n个信号强度实测值CEST_ideal1、CEST_ideal2、CEST_ideal3、…、CEST_idealn的匹配程度最高，则理论化学交换特征定量(K_i,c_j)即为CEST对比剂的化学交换特征定量结果。

对每个CEST对比剂的CEST信号的信号强度实测值均与对应的多个信号强度理论值进行匹配，匹配结果对应的理论化学交换特征定量即为每个CEST对比剂的化学交换特征定量结果。

本申请的实施例针对于每一个CEST对比剂的CEST信号的信号强度实测值，在信号强度理论值进行计算时与信号强度实测值在采集时的采集环境，即饱和脉冲参数、共振频率相同，因此基于理论化学交换特征定量范围，可以计算出多个与信号强度实测值对应的信号强度理论值；由于信号强度理论值是在与信号强度实测值采集时相同的采集环境下计算得到的多个值，采用一对多的匹配方式使得信号强度实测值能有更准确的匹配结果。本申请的化学交换特征定量过程无需选择参照信号，因此不易受到其他化学交换效应的影响，由此使得匹配结果对应的理论化学交换特征定量能更准确地代表实际的化学交换特定量，故有效地提高了化学交换特定量的准确率。

在一个示例中，S101之前还可以包括S100-1和S100-2：

S100-1，根据各CEST对比剂的化学交换特征，确定各CEST对比剂的有效成像范围；所述有效成像范围中的参数为饱和脉冲参数。

示例性的，CEST对比剂的化学交换特征可以是例如化学交换速率和/或浓度等。有效成像范围即为针对于该CEST对比剂的饱和脉冲参数施加的一个有效范围，在该范围内，CEST对比剂具有较好的成像效果。各CEST对比剂的有效成像范围可通过实验来确定，即针对于不同的饱和脉冲参数，分别进行CEST对比剂的成像实验，然后根据成像效果来选定一个具有较好成像效果的饱和脉冲参数范围，作为该CEST对比剂的有效成像范围。还可以通过先验知识进行确定，例如通过查询已有资料中该CEST对比剂成像时惯用的饱和脉冲参数范围作为有效成像范围。

S100-2，根据所有所述CEST对比剂的有效成像范围，确定所述饱和脉冲的饱和脉冲参数范围；其中，所述饱和脉冲参数范围中的饱和脉冲参数属于各CEST对比剂的有效成像范围。

在S100-2中，选定的饱和脉冲参数范围中的饱和脉冲参数，应该同时属于所有CEST对比剂的有效成像范围，这样在施加饱和脉冲时，随机选择的饱和脉冲可在不同的CEST对比剂的共振频率上满足该CEST对比剂的成像要求。

需要说明的是，当饱和脉冲参数选择一种以上的组合时，例如选择脉冲强度ω₁、持续时间T_s时，所对应的饱和脉冲参数范围也是两个，即对应于脉冲强度的饱和脉冲参数范围以及对应于持续时间的饱和脉冲参数范围。

以饱和脉冲参数选择两种，即选择脉冲强度、持续时间这两种饱和脉冲参数为例，如针化学交换速率较慢(10-100Hz)的CEST对比剂，其有效成像范围可以是：饱和脉冲强度为0.5-4.0μT，持续时间为30-100ms，代表性的CEST对比剂有-NH等；而针对交换速率较快(500-1500Hz)的CEST对比剂，其有效成像范围可以是：饱和脉冲强度>2.0μT，持续时间<50ms，代表性的CEST对比剂有-OH等。该示例中，饱和脉冲参数范围可以选择：饱和脉冲强度0.5～5μT，持续时间10～100ms。

现有技术中，CEST成像中饱和脉冲参数的选取往往由经验而来，且在数据采集过程中其参数保持恒定，每次仅能对单一CEST对比剂实施有效定量，一方面较难获取最佳CEST效应，另一方面，又难以同时对具备不同化学交换特性的多种CEST对比剂实施有效定量。本申请的实施例中，根据CEST对比剂的化学交换特征，首先确定各CEST对比剂的有效成像范围，在此基础上确定饱和脉冲参数范围。在进行CEST定量时，以设定的步长改变饱和脉冲频率，可依次在每个CEST对比剂的共振频率上进行饱和脉冲的施加，由此可对不同CEST对比剂进行同步定量；且依据CEST对比剂化学交换特征确定的饱和脉冲参数范围内采集的CEST数据能更好地获取CEST效应。

请参阅图3，在一种可能的实现方式中，S101可以包括S1011～S1014，该示例中采用设定频率范围并从中选取多个频率，在每个频率上均施加饱和脉冲，从探测到的所有的CEST信号的信号强度实测值中提取每个CEST对比剂对应的CEST信号的信号强度实测值，具体如下：

S1011，根据多种CEST对比剂的共振频率，确定待施加的饱和脉冲的频率范围；其中，饱和脉冲的频率范围中的频率涵盖各CEST对比剂的共振频率。

可选地，将各CEST对比剂的共振频率对应的频率范围与主磁场的频率偏移程度范围求并集，将结果作为饱和脉冲的频率范围，由此确定的频率范围中即涵盖各CEST对比剂的共振频率。作为一种通用的表述习惯，常采用频率偏移对频率进行说明。频率偏移和频率之间存在换算关系，具体为：

频率偏移＝(f-f₀)×10⁶/f₀ppm 式1

式1中，f表示CEST对比剂的共振频率或施加饱和脉冲的频率，f₀为参考频率，例如采用水的共振频率作为参考频率。

例如，内源性CEST对比剂相对于水分子的共振频率偏移通常在-4～4ppm，若主磁场的频率偏移在-1～1ppm范围内，则饱和脉冲的频率偏移范围可以选择-5～5ppm。

本实施例中，以同时探测生物体组织内的三种CEST对比剂：氨基(-NH，3.5ppm)、胺基(-NH₂，2ppm)、谷氨酸盐(3ppm)为例，这三种CEST对比剂的共振频率偏移对应的频率偏移范围为2～3.5ppm，在主磁场频率偏移为-1～1ppm时，饱和脉冲的频率偏移范围可选择1～4.5ppm。确定频率偏移范围之后，通过式1将频率偏移范围转换成频率，即得到频率范围。

S1012，从饱和脉冲的频率范围中选择多个频率，多个频率包含各CEST对比剂的共振频率。

在确定饱和脉冲的频率范围后，以合适的步长对频率范围进行划分。其中步长的选取应该满足使得划分后的多个频率包含各CEST对比剂的共振频率。以频率偏移为例进行说明：

在饱和脉冲的频率偏移范围1～4.5ppm内，选择适当的频率偏移步长，例如待探测的CEST对比剂为氨基(-NH，3.5ppm)、胺基(-NH₂，2ppm)、谷氨酸盐(3ppm)三种，共振频率偏移分别为3.5ppm、2ppm、3ppm，则可以选择频率偏移步长为0.5ppm；即频率偏移步长满足被待探测的多种CEST对比剂的共振频率偏移整除时，即可满足划分后的多个频率偏移满足涵盖各个CEST对比剂共振频率的要求。将偏移步长、频率偏移按照式1转换成频率，即得到步长以及划分后的频率。

S1013，在多个频率中的每个频率上按预设条件施加饱和脉冲时，采集每个频率对应的CEST信号的信号强度实测值。

在每个频率上按照预设条件施加饱和脉冲后，均可以探测到对应于该频率的CEST信号的信号强度实测值。以将频率范围划分成多个频率[f1,f2,f3,…]，每个频率上施加n组饱和脉冲进行说明。

首先在第一个频率f1上施加饱和脉冲，施加之前先在设定的饱和脉冲参数范围内随机化选择饱和脉冲参数，按照选择的饱和脉冲参数施加饱和脉冲，即可探测到在该频率上的第一个CEST信号的信号强度实测值。同样是在该频率上，再次随机化选择饱和脉冲参数，按照随机化后的饱和脉冲参数施加饱和脉冲，探测到该频率上的第二个CEST信号的信号强度实测值；按照同样的方法，直至n个饱和脉冲均施加完毕之后，得到所有针对于第一个频率f1的n个CEST信号的信号强度实测值。在其他的频率f2,f3,…都依照此方法，即可得到每个频率对应的CEST信号的信号强度实测值。

S1014，从每个频率对应的CEST信号的信号强度实测值中，提取每个CEST对比剂的共振频率对应的CEST信号的信号强度实测值作为CEST对比剂的CEST信号的信号强度实测值。

由于设定的频率范围中的频率涵盖每个CEST对比剂的共振频率，所以针对一个CEST对比剂，通过提取出与其共振频率相同的饱和脉冲施加频率时所探测到的CEST信号的信号强度实测值，即为该CEST对比剂的CEST信号的信号强度实测值。

因为每个频率上施加的饱和脉冲是n组，故提取出的CEST对比剂的CEST信号的实测值CEST_acq也有n个，这n个实测值反映出了该CEST对比剂的CEST信号的信号强度随n个饱和脉冲参数的变化情况。

由于在CEST信号的信号强度实测值采集时可能存在主磁场偏移，实际采集的频率点会有偏差，故本实施例采用上述首先设定频率范围，然后将在频率范围内选择多个频率，在每个频率上分别采集信号强度实测值，在从所有信号强度实测值中提取每个CEST对比剂的CEST信号的信号强度实测值的方法，该方法可在每个CEST对比剂的共振频率附近采集多个频率点的数据，然后采用水饱和移位(water saturation shift referencing,WASSR)校正法等方法进行主磁场校正以获取频率校正后各CEST对比剂的CEST信号的信号强度实测值；校正后的实测值具有更高的准确性，后续的处理过程采用的是校正后的CEST信号的信号强度实测值。

在进行CEST对比剂化学交换特征定量过程中，自由水的磁化转移(MagnetizationTransfer，MT)效应相对较强，干扰CEST信号。

在一个可选的示例中，为避免MT效应对CEST信号强度的影响，S102之前，还包括：

S1-102，量化所述目标物体的磁化转移效应，得到磁化转移效应的化学交换特征定量结果。

其中，S1-102在S101之后执行。

基于上述磁化转移效应的化学交换特征定量结果，如图4所示，S102可以具体包括S1021～S1023：

S1021，获取每种CEST对比剂的CEST信号对应的第一生成条件，第一生成条件包括CEST对比剂的共振频率，以及施加饱和脉冲时对应的饱和脉冲参数。

第一生成条件，即S101中采集CEST信号的信号强度实测值时对应的条件。在S101中采集CEST信号的信号强度实测值时，需要在CEST对比剂的共振频率上按照预设条件施加饱和脉冲。其中预设条件中包含饱和脉冲参数范围，采集信号强度实测值时饱和脉冲对应的饱和脉冲参数在设定的饱和脉冲参数范围内随机化；例如上述示例中，对于第1个饱和脉冲施加后采集到的CEST信号，第一生成条件为CEST对比剂的共振频率F₁，以及饱和脉冲参数：脉冲强度ω₁、持续时间T_s。

S1022，确定每种CEST对比剂的理论化学交换特征定量范围。

可选地，采用以往的实验数据来确定每种CEST对比剂的理论化学交换特征定量范围。例如，根据以往的实验数据，CEST对比剂氨基(-NH，3.5ppm)的理论化学交换特征定量范围可以是：化学交换速率范围50-500Hz，浓度范围100-1000mM。

S1023，基于理论化学交换特征定量范围中不同的理论化学交换特征定量，利用磁化转移效应的化学交换特征定量结果计算在第一生成条件下生成的CEST信号所对应的多个信号强度理论值。

可选地，在一个实施例中，计算CEST信号对应的多个信号强度理论值时利用三化学交换池Bloch方程。其中，三化学交换池Bloch方程是指自由水、磁化转移MT Bloch方程；例如对于第1个饱和脉冲施加后采集的CEST信号，通过该方程，利用磁化转移效应的化学交换特征定量结果、第一生成条件中的共振频率F1，以及饱和脉冲参数：脉冲强度ω₁、持续时间T_s，可计算得到对于CEST对比剂的不同理论化学交换特征定量时的多个信号强度理论值。

三化学交换池Bloch方程考虑了自由水和磁化转移效应，因此在CEST信号的信号强度理论值时可有效地排除自有水的直接饱和效应和磁化转移效应，使得理论值的计算结果更加精准。

可选地，根据计算得到的CEST信号的信号强度的理论值，可进一步建立每一个CEST对比剂的CEST信号字典集以保存每一个CEST信号的多个信号强度实测值，以便于后续的匹配过程。

请参阅图5，作为一个具体的示例，如图5所示，S1-102还可以包括：S1-1021～S1-1023：

S1-1021，获取在远离自由水和所有CEST对比剂的共振频率上，按照预设条件施加饱和脉冲时，采集的磁化转移信号的信号强度实测值。

其中，远离自由水和所有CEST对比剂的共振频率，例如选择共振频率偏移为10ppm处。自由水的共振频率偏移为0ppm，本实施例中列举的三种CEST对比剂：氨基(-NH，3.5ppm)、胺基(-NH₂，2ppm)、谷氨酸盐(3ppm)的共振频率偏移范围为2～3.5ppm，因此10ppm的共振频率偏移相对于自由水和CEST对比剂的共振频率已足够远，故在10ppm的共振频率偏移处所对应的共振频率上施加饱和脉冲可以满足获取磁化转移效应的化学交换特征定量的需求。10ppm处的共振频率已经是一个比较远的共振频率，满足一般情况下多数CEST对比剂进行磁化转移效应的化学交换特征定量的需求。需要说明的是，此处的列举只是示意性的，不是限制性的；也可以选择8ppm、9ppm、11ppm等其他共振频率进行饱和脉冲的施加。“远离”，可以理解成是在所有CEST对比剂的共振频率中，选择最大的一个共振频率值，在这个值对应的共振频率偏移的数值上加f，得到的共振频率偏移所对应的共振频率即可理解为远离自由水和CEST对比剂的共振频率；其中f可以是一个自然数，例如其取值>2ppm。

在确定好上述远离自由水和所有CEST对比剂的共振频率后，按照预设条件对目标物体施加饱和脉冲。这里的预设条件与S101中预设条件相同，即依次施加至少一个饱和脉冲，每个饱和脉冲参数在设定的饱和脉冲参数范围内随机化。其中，饱和脉冲参数范围也可以与S101中相同。

在确定的共振频率上施加饱和脉冲，每个饱和脉冲的饱和脉冲参数在设定的饱和脉冲参数范围内随机化，采集施加饱和脉冲后产生的磁化转移信号的信号强度实测值MT_acq。饱和脉冲的数量与采集的磁化转移信号的数量存在对应关系，当依次施加多个饱和脉冲时，可对应采集到多个磁化转移信号的信号强度实测值。本实施例中，采集磁化转移信号时与采集CEST信号施加饱和脉冲数量一致，即最终采集到n个磁化转移信号的信号强度实测值。

S1-1022，基于磁化转移信号的理论化学交换特征定量范围，计算磁化转移信号对应的多个信号强度理论值；其中，在计算磁化转移信号的信号强度理论值时所对应的饱和脉冲参数以及共振频率，与生成磁化转移信号的信号强度实测值时所对应的饱和脉冲参数、共振频率相同。

该步骤中，通过理论计算的方法确定磁化转移信号在不同理论化学交换特征定量时的多个信号强度理论值。其中，化学交换特征可以是化学交换速率和/或浓度等，在化学交换特征选择化学交换速率和浓度时，理论化学交换特征定量范围指化学交换速率范围和浓度范围。磁化转移信号的理论化学交换特征定量范围可以通过先验知识确定，例如通过以往的试验数据，或通过查询资料的方式来确定。

基于磁化转移信号的理论化学交换特征定量范围，在与磁化转移信号的信号强度实测值采集时相同的饱和脉冲参数、共振频率条件下，可计算多个与磁化转移信号对应的信号强度理论值。例如，可利用理论化学交换特征定量范围中不同的理论化学交换特征定量进行计算，针对于一个磁化转移信号的理论化学交换特征定量，通过理论计算的方式可得到一个信号强度理论值；那么将范围中的多个理论化学交换特征定量分别计算信号强度理论值，即可得到对应于磁化转移信号的多个信号强度理论值。由于针对一磁化转移信号可探测到一个信号强度实测值，因此信号强度理论值与信号强度实测值存在多对一的对应关系。

当针对于目标物体施加多个饱和脉冲而采集了多个磁化转移信号时，在计算每个磁化转移信号的对应的多个信号强度理论值过程中，首先从磁化转移信号的理论化学交换特征定量范围中选取多个理论化学交换特征定量，按照取出的化学交换特征定量，依次计算在每个化学交换特征定量下，每个磁化转移信号的信号强度理论值。

例如，在采集磁化转移信号的信号强度实测值时确定的远离自由水和所有CEST对比剂的共振频率为F2，预设条件中，设定施加n个饱和脉冲，设定的饱和脉冲参数范围是：脉冲强度ω₀～ω_a、持续时间T_s0～T_sb。磁化转移信号的理论化学交换特征定量范围为：化学交换速率范围K’，浓度范围c’。在理论化学交换特征定量范围中选择多个理论化学交换特征定量，例如在化学交换速率范围K’中等间隔选取a个化学交换速率K’₁、K’₂…K’_a，在浓度范围c’中等间隔选取b个浓度c’₁、c’₂…c’_b，则化学交换速率、浓度共可以构成a×b个组合(K’₁,c’₁)、(K’₁,c’₂)、(K’₁,c’₃)…(K’_a,c’_b)，也即a×b个理论化学交换特征定量。

在采集磁化转移信号的信号强度实测值时：

在共振频率F2上施加第1至第n个饱和脉冲，随机选择的饱和脉冲参数(脉冲强度,持续时间)分别为：(ω’₁，T’_s1)、(ω’₂，T’_s2)…(ω’_n，T’_sn)，其中，ω’₁,ω’₂,…ω’_n∈(ω₀～ω_a)，T’_s1,T’_s2,…T’_sn∈(T_s0～T_sb)。则采集到的第1个至第n个磁化转移信号的信号强度实测值分别为MT_acq1、MT_acq2…MT_acqn。

在计算磁化转移信号对应的多个信号强度理论值时：

对于第1个理论化学交换特征定量，本例中为第1个化学交换速率和第1个浓度(K’₁,c’₁)，计算第1～n个磁化转移信号的信号强度理论值时，分别将共振频率F2，以及采集1～n个磁化转移信号的信号强度实测值时的饱和脉冲参数(ω’_1，T’_s1)、(ω’_2，T’_s2)…(ω’_n，T’_sn)作为计算条件，计算得到第1～n个信号强度理论值MT¹ _ideal1、MT¹ _ideal2…MT¹ _idealn。即，第1个信号强度理论值MT¹ _ideal1在计算时所用的参量为：[(K’₁,c’₁)、F2、(ω’₁，T’_s1)]；第2个信号强度理论值MT¹ _ideal2在计算时所用的参量为：[(K’₁,c’₁)、F2、(ω’₂，T’_s2)]…第n个信号强度理论值MT¹ _idealn在计算时所用的参量为：[(K’₁,c’₁)、F2、(ω’_n，T’_sn)]。

通过以上方法，计算得到了针对于第1个理论化学交换特征定量的n个磁化转移信号的n个信号强度理论值。

对于第2个理论化学交换特征定量(K’₁,c’₂)，计算第1～n个CEST信号的信号强度理论值的参量分别为[(K’₁,c’₂)、F2、(ω’₁，T’_s1)]、[(K’₁,c’₂)、F2、(ω’₂，T’_s2)]…[(K’₁,c’₂)、F2、(ω’_n，T’_sn)]，对应计算出1～n个信号强度理论值MT² _ideal1、MT² _ideal2…MT² _idealn。

按照相同的方法，可依次计算每一个理论化学交换特征定量对应第1至第n个磁化转移信号的信号强度理论值。由于每一个理论化学交换特征定量均可对应计算出一个信号强度理论值，则从理论化学交换特征定量范围中选取多个理论化学交换特征定量后，针对于一个磁化转移信号通过理论计算可得到与其信号强度实测值对应的多个信号强度理论值。

S1-1023，将磁化转移信号的信号强度实测值与其对应的多个信号强度理论值进行匹配，将匹配度最高的信号强度理论值对应的理论化学交换特征定量作为目标物体的磁化转移效应的化学交换特征定量结果。

可选地，利用点乘最大化原则等算法，将磁化转移信号的信号强度实测值与对应的多个信号强度理论值进行匹配；例如，施加n个饱和脉冲后记录磁化转移信号的n个信号强度实测值为MT_acq1、MT_acq2、MT_acq3…、MT_acqn，在磁化转移信号字典集中存储的对应于这n个信号强度实测值的多个信号强度理论值中，寻找与实测值MT_acq1、MT_acq2、MT_acq3…、MT_acqn匹配度最高的n个理论值，这n个理论值在计算过程中对应的理论化学交换特征定量即为目标物体的磁化转移效应的化学交换特征定量。

例如，在S1-1022计算时采用的理论化学交换特征定量中，第k个理论化学交换特征定量(K’_i,c’_j)∈[(K’₁,c’₁)、(K’₁,c’₂)、(K’₁,c’₃)…(K’_a,c’_b)]所对应的n个理论值MT^k _ideal1、MT^k _ideal2、MT^k _ideal3…MT^k _idealn与n个信号强度实测值MT_acq1、MT_acq2、MT_acq3…、MT_acqn匹配度最高，则理论化学交换特征定量(K’_i,c’_j)即为目标物体的磁化转移效应的化学交换特征定量结果。

本申请的实施例中，采用计算磁化转移信号的信号强度理论值的方法，使得对于每一个磁化转移信号的信号强度实测值，都有多个对应的信号强度理论值；由于这些理论值是基于理论化学交换特征定量范围、并且在与实测值采集时相同的共振频率和饱和脉冲参数条件下计算得出的，故实测值可精确地找到匹配度最高的理论值，从而使得磁化转移效应定量结果更加精确。可以在理论化学交换特征定量范围中选择更多的定量值计算更多的信号强度理论值，以进一步提升磁化转移效应定量结果的精确程度。

在一个可选的示例中，如图6所示，S1-1022还可以包括S1-10221～S1-10223：

S1-10221，获取磁化转移信号对应的第二生成条件，第二生成条件包括远离自由水的信息和所有CEST对比剂的共振频率，以及施加饱和脉冲时对应的饱和脉冲参数。

第二生成条件，即S1-10221中采集磁化转移信号的信号强度实测值时对应的条件。在S1-10221中采集磁化转移信号的信号强度实测值时，需要确定在远离自由水和所有CEST对比剂的共振频率，并按照预设条件施加饱和脉冲。其中预设条件中包含饱和脉冲参数范围，采集信号强度实测值时饱和脉冲对应的饱和脉冲参数在设定的饱和脉冲参数范围内随机化；例如上述示例中，对于第1个饱和脉冲施加后采集到的磁化转移信号，第二生成条件为远离自由水和所有CEST对比剂的共振频率为F2，以及饱和脉冲参数：脉冲强度ω’₁、持续时间T’_s。

S1-10222，确定磁化转移信号的理论化学交换特征定量范围；

该示例中，磁化转移信号的理论化学交换特征定量范围通过以往的试验数据得到，示例性的，磁化转移信号的理论化学交换特征定量范围中，化学交换速率K范围设置为5-100Hz，浓度c范围2-30M。

S1-10223，基于理论化学交换特征定量范围中不同的理论化学交换特征定量，计算在第二生成条件下生成的磁化转移信号所对应的多个信号强度理论值。

可选地，在一个实施例中，计算磁化转移信号对应的多个信号强度理论值时利用双化学交换池Bloch方程。其中，双化学交换池Bloch方程是指自由水、磁化转移MT Bloch方程；例如对于第1个饱和脉冲施加后采集的磁化转移信号，通过该方程，利用第二生成条件中的共振频率为F2，以及饱和脉冲参数：脉冲强度ω’₁、持续时间T’_s，可计算得到对于磁化转移信号的不同理论化学交换特征定量时的多个信号强度理论值。

可选地，根据计算得到磁化转移信号的信号强度的理论值，可进一步建立磁化转移信号字典集以保存每一个磁化转移信号的多个信号强度实测值，以便于后续的匹配过程。

如图7所示，图7中展示出了一个CEST对比剂的CEST信号的多个信号强度理论值。为了更加直观地看到信号强度理论值的变化情况，图7中以饱和脉冲组的序号作为横坐标轴，将不同饱和脉冲对应的CEST信号强度的理论值依次连接，即得到了不同理论化学交换特征定量所对应的信号强度理论值随饱和脉冲序号的变化曲线。图7的示例中，共有8条变化曲线，即有8个理论化学交换特定量对应的信号强度理论值变化曲线。如图8所示，为CEST信号的实测值随饱和脉冲序号的变化曲线。如图9所示，将图8中的曲线在图7中进行匹配，匹配度最高的一条变化曲线中各个信号强度理论值对应的理论化学交换特征定量即CEST对比剂的化学交换特征定量。从图9的匹配结果可以看出，本申请实施例的方法有很高的匹配程度。

本申请的实施例中对磁化转移效应进行了定量，在进行CEST信号的信号强度理论值计算过程中，消除了磁化转移效应对计算环境的影响，从而使得计算得到的信号强度理论值是在更为贴近目标物体真实环境下得出的，相比于传统方法依赖参考信号的准确性的局限，本申请中由于能够对磁化转移效应进行准确量化，因此能有效地去除磁化转移效应，进一步地提升了定量结果的准确率。

将每个CEST对比剂的化学交换特征定量结果以伪彩图等方式制图，即可得到各CEST对比剂的同步定量成像结果。

图10示出了本申请提供的一种化学交换特征定量装置，本实施例中的化学交换特征定量装置包含的各模块用于执行上述方法实施例中图1、图3～图6各自对应的实施例，具体实现过程可参阅上述实施例的相关描述，此次不赘述。化学交换特征定量装置包括：

第一采集模块1，用于采集目标物体中待探测的多种化学交换饱和转移CEST对比剂各自的CEST信号的信号强度实测值；其中，CEST信号在对CEST对比剂的共振频率上按预设条件施加饱和脉冲时产生；

第一计算模块2，用于基于每种CEST对比剂对应的理论化学交换特征定量范围，计算每种CEST对比剂各自的CEST信号对应的多个信号强度理论值；其中，每个CEST信号的信号强度理论值在计算时对应的饱和脉冲参数、共振频率与其信号强度实测值对应的饱和脉冲参数、共振频率相同；

第一匹配模块3，用于将每种CEST对比剂各自对应的CEST信号的信号强度实测值与其对应的多个信号强度理论值进行匹配，将匹配度最高的信号强度理论值对应的理论化学交换特征定量作为CEST对比剂的化学交换特征定量结果。

可选地，化学交换特征定量装置还包括：

第一设定模块，用于根据各CEST对比剂的化学交换特征，确定各CEST对比剂的有效成像范围；有效成像范围中的参数为饱和脉冲参数；

第二设定模块，用于根据所有CEST对比剂的有效成像范围，确定饱和脉冲的饱和脉冲参数范围；其中，饱和脉冲参数范围中的饱和脉冲参数属于各CEST对比剂的有效成像范围。

进一步地，预设条件为：依次施加至少一个饱和脉冲，每个饱和脉冲参数在设定的饱和脉冲参数范围内随机化。

可选地，如图11所示，第一采集模块1包括：

频率范围设定单元11，用于根据多种CEST对比剂的共振频率，确定待施加的饱和脉冲的频率范围；其中，饱和脉冲的频率范围中的频率涵盖各CEST对比剂的共振频率；

频率选择单元12，用于从饱和脉冲的频率范围中选择多个频率，多个频率包含各CEST对比剂的共振频率；

第一采集子单元13，用于在多个频率中的每个频率上按预设条件施加饱和脉冲时，采集每个频率对应的CEST信号的信号强度实测值；

提取单元14，用于从每个频率对应的CEST信号的信号强度实测值中，提取每个CEST对比剂的共振频率对应的CEST信号的信号强度实测值作为CEST对比剂的CEST信号的信号强度实测值。

进一步地，化学交换特征定量装置还包括：

量化模块4，用于量化目标物体的磁化转移效应，得到磁化转移效应的化学交换特征定量结果；

如图12所示，第一计算模块2包括：

第一获取单元21，用于获取每种CEST对比剂的CEST信号对应的第一生成条件，第一生成条件包括CEST对比剂的共振频率，以及施加饱和脉冲时对应的饱和脉冲参数；

第一确定单元22，用于确定每种CEST对比剂的理论化学交换特征定量范围；

第一计算子单元23，用于基于理论化学交换特征定量范围中不同的理论化学交换特征定量，利用磁化转移效应的化学交换特征定量结果计算在第一生成条件下生成的CEST信号所对应的多个信号强度理论值。

可选地，如图13所示，量化模块4包括：

第二采集模块41，用于获取在远离自由水和所有CEST对比剂的共振频率上，按照预设条件施加饱和脉冲时，采集的磁化转移信号的信号强度实测值；

第二计算模块42，用于基于磁化转移信号的理论化学交换特征定量范围，计算磁化转移信号对应的多个信号强度理论值；其中，在计算磁化转移信号的信号强度理论值时所对应的饱和脉冲参数以及共振频率，与生成磁化转移信号的信号强度实测值时所对应的饱和脉冲参数、共振频率相同；

第二匹配模块43，用于将磁化转移信号的信号强度实测值与其对应的多个信号强度理论值进行匹配，将匹配度最高的信号强度理论值对应的理论化学交换特征定量作为目标物体的磁化转移效应的化学交换特征定量结果。

进一步地，第二计算模块包括：

第二获取单元，用于获取磁化转移信号对应的第二生成条件，第二生成条件包括远离自由水的信息和所有CEST对比剂的共振频率，以及施加饱和脉冲时对应的饱和脉冲参数；

第二确定单元，用于确定磁化转移信号的理论化学交换特征定量范围；

第二计算子单元，用于基于理论化学交换特征定量范围中不同的理论化学交换特征定量，计算在第二生成条件下生成的磁化转移信号所对应的多个信号强度理论值。

需要说明的是，上述模块/单元之间的信息交互、执行过程等内容，由于与本申请方法实施例基于同一构思，其具体功能及带来的技术效果，具体可参见方法实施例部分，此处不再赘述。

请参阅图14，本申请实施例进一步提供一种化学交换特征定量设备5，如图14所示，包括存储器52、处理器51以及存储在存储器52中并可在处理器上运行的计算机程序53，处理器51执行计算机程序53时实现上述化学交换特征定量方法，例如，图1所示的S101至S103。

计算机程序53也可以被分割成一个或多个模块/单元，一个或者多个模块/单元被存储在存储器52中，并由处理器51执行，以完成本申请。一个或多个模块/单元可以是能够完成特定功能的一系列计算机程序指令段，该指令段用于描述计算机程序53在化学交换特征定量设备5中的执行过程，例如，计算机程序53可以被分割为第一采集模块、第一计算模块、第一匹配模块，各模块的功能参见前述系统中的描述，不再赘述。

本申请的实施提供了一种计算机可读存储介质，计算机可读存储介质存储有计算机程序，计算机程序被处理器执行时实现上述化学交换特征定量方法，例如，图1所示的S101至S103。

集成的模块/单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本申请实现上述实施例方法中的全部或部分流程，也可以通过计算机程序来指令相关的硬件来完成，计算机程序可存储于一计算机可读存储介质中，该计算机程序在被处理器执行时，可实现上述各个方法实施例的步骤。其中，计算机程序包括计算机程序代码，计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。计算机可读介质可以包括：能够携带计算机程序代码的任何实体或装置、记录介质、U盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、电载波信号、电信信号以及软件分发介质等。需要说明的是，计算机可读介质包含的内容可以根据司法管辖区内立法和专利实践的要求进行适当的增减，例如在某些司法管辖区，根据立法和专利实践，计算机可读介质不包括电载波信号和电信信号。

本申请的实施例还提供了一种计算机程序产品，当计算机程序产品在化学交换特征定量设备上运行时，使得化学交换特征定量设备执行上述化学交换特征定量方法，例如，图1所示的S101至S103。

以上实施例仅用以说明本申请的技术方案，而非对其限制；尽管参照前述实施例对本申请进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本申请各实施例技术方案的精神和范围，均应包含在本申请的保护范围之内。

Claims (10)

1.一种化学交换特征定量方法，其特征在于，包括：

采集目标物体中待探测的多种化学交换饱和转移CEST对比剂各自的CEST信号的信号强度实测值；其中，所述CEST信号在对所述CEST对比剂的共振频率上按预设条件施加饱和脉冲时产生；

基于每种所述CEST对比剂对应的理论化学交换特征定量范围，计算每种所述CEST对比剂各自的CEST信号对应的多个信号强度理论值；其中，每个CEST信号的信号强度理论值在计算时对应的饱和脉冲参数、共振频率与其信号强度实测值对应的饱和脉冲参数、共振频率相同；

将每种所述CEST对比剂各自对应的CEST信号的信号强度实测值与其对应的多个信号强度理论值进行匹配，将匹配度最高的信号强度理论值对应的理论化学交换特征定量作为所述CEST对比剂的化学交换特征定量结果。

2.根据权利要求1所述的化学交换特征定量方法，其特征在于，所述采集目标物体中待探测的多种化学交换饱和转移CEST对比剂各自的CEST信号的信号强度实测值之前，还包括：

根据各CEST对比剂的化学交换特征，确定各CEST对比剂的有效成像范围；所述有效成像范围中的参数为饱和脉冲参数；

根据所有所述CEST对比剂的有效成像范围，确定所述饱和脉冲的饱和脉冲参数范围；其中，所述饱和脉冲参数范围中的饱和脉冲参数属于各CEST对比剂的有效成像范围。

3.根据权利要求1或2所述的化学交换特征定量方法，其特征在于，所述预设条件为：

依次施加至少一个饱和脉冲，每个所述饱和脉冲参数在设定的饱和脉冲参数范围内随机化。

4.根据权利要求1所述的化学交换特征定量方法，其特征在于，所述采集目标物体中待探测的多种化学交换饱和转移CEST对比剂各自的CEST信号的信号强度实测值，包括：

根据多种所述CEST对比剂的共振频率，确定待施加的饱和脉冲的频率范围；其中，所述饱和脉冲的频率范围中的频率涵盖各CEST对比剂的共振频率；

从所述饱和脉冲的频率范围中选择多个频率，所述多个频率包含各CEST对比剂的共振频率；

在所述多个频率中的每个频率上按所述预设条件施加饱和脉冲时，采集每个频率对应的CEST信号的信号强度实测值；

从所述每个频率对应的CEST信号的信号强度实测值中，提取每个CEST对比剂的共振频率对应的CEST信号的信号强度实测值作为所述CEST对比剂的CEST信号的信号强度实测值。

5.根据权利要求1所述的化学交换特征定量方法，其特征在于，所述基于每种所述CEST对比剂对应的理论化学交换特征定量范围，计算每种所述CEST对比剂各自的CEST信号对应的多个信号强度理论值之前，还包括：

量化所述目标物体的磁化转移效应，得到磁化转移效应的化学交换特征定量结果；

所述基于每种所述CEST对比剂对应的理论化学交换特征定量范围，计算每种所述CEST对比剂各自的CEST信号对应的多个信号强度理论值，包括：

获取每种所述CEST对比剂的CEST信号对应的第一生成条件，所述第一生成条件包括所述CEST对比剂的共振频率，以及施加饱和脉冲时对应的饱和脉冲参数；

确定每种所述CEST对比剂的理论化学交换特征定量范围；

基于所述理论化学交换特征定量范围中不同的理论化学交换特征定量，利用所述磁化转移效应的化学交换特征定量结果计算在所述第一生成条件下生成的CEST信号所对应的多个信号强度理论值。

6.根据权利要求5所述的化学交换特征定量方法，其特征在于，所述量化目标物体的磁化转移效应，得到磁化转移效应的化学交换特征定量结果，包括：

获取在远离自由水和所有CEST对比剂的共振频率上，按照所述预设条件施加饱和脉冲时，采集的磁化转移信号的信号强度实测值；

基于所述磁化转移信号的理论化学交换特征定量范围，计算所述磁化转移信号对应的多个信号强度理论值；其中，在计算磁化转移信号的信号强度理论值时所对应的饱和脉冲参数以及共振频率，与生成磁化转移信号的信号强度实测值时所对应的饱和脉冲参数、共振频率相同；

将所述磁化转移信号的信号强度实测值与其对应的多个信号强度理论值进行匹配，将匹配度最高的信号强度理论值对应的理论化学交换特征定量作为所述目标物体的磁化转移效应的化学交换特征定量结果。

7.根据权利要求6所述的化学交换特征定量方法，其特征在于，所述基于所述磁化转移信号的理论化学交换特征定量范围，计算所述磁化转移信号对应的多个信号强度理论值，包括：

获取所述磁化转移信号对应的第二生成条件，所述第二生成条件包括所述远离自由水的信息和所有CEST对比剂的共振频率，以及施加饱和脉冲时对应的饱和脉冲参数；

确定所述磁化转移信号的理论化学交换特征定量范围；

基于所述理论化学交换特征定量范围中不同的理论化学交换特征定量，计算在所述第二生成条件下生成的磁化转移信号所对应的多个信号强度理论值。

8.根据权利要求1、2、4至7中任一权利要求所述的化学交换特征定量方法，其特征在于，所述饱和脉冲参数包括：脉冲强度、持续时间、翻转角、占空比中的一种或多种。

9.一种化学交换特征定量设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，其特征在于，处理器执行计算机程序时实现如权利要求1-8中任一项所述的化学交换特征定量方法的步骤。

10.一种计算机可读存储介质，计算机可读存储介质存储有计算机程序，其特征在于，计算机程序被处理器执行时实现如权利要求1-8中任一项所述的化学交换特征定量方法的步骤。

Images