射频系统中心频率确定方法、装置、计算机设备和介质

摘要

本申请涉及一种射频系统中心频率确定方法、装置、计算机设备、存储介质和程序产品，该方法包括采集第一磁共振信号和第二磁共振信号；确定第一磁共振信号对应的第一频谱，以及第二磁共振信号对应的第二频谱；根据第一频谱和第二频谱确定相对弛豫频谱；通过分析相对弛豫频谱确定目标成像核素在目标组织中的特征频率，根据特征频率确定射频系统的中心频率。通过本申请提供的射频系统中心频率确定方法可以得到更加准确的中心频率。

说明书

技术领域

本申请涉及磁共振成像技术领域，特别是涉及一种射频系统中心频率确定方法、装置、计算机设备和存储介质。

背景技术

在临床磁共振设备的成像过程中，通过磁共振信号对应的频谱信息确定磁共振设备中射频系统的中心频率，以使射频系统的中心频率与成像核素在目标组织中的拉莫频率对准，从而达到最佳的成像效果。

传统技术中，通过对采集到的磁共振信号的谱线特征进行分析，即采用单一的谱线分析方法来确定成像核素在目标组织中的中心频率。然而，面对临床扫描的复杂多变场景，上述单一的谱线分析方法难以适应，进而存在由此分析确定的中心频率不准确的问题。

发明内容

基于此，有必要针对上述技术问题，提供一种射频系统中心频率确定方法、装置、计算机设备、存储介质和程序产品。

第一方面，本申请一个实施例提供一种射频系统中心频率确定方法，包括：

采集第一磁共振信号和第二磁共振信号，第一磁共振信号是由第一脉冲序列得到，第二磁共振信号是由第二脉冲序列得到，第一脉冲序列与第二脉冲序列相邻设置；

确定第一磁共振信号对应的第一频谱，以及第二磁共振信号对应的第二频谱；

根据第一频谱和第二频谱，确定相对弛豫频谱；

通过分析相对弛豫频谱确定目标成像核素在目标组织的特征频率，根据特征频率确定射频系统的中心频率。

在其中一个实施例中，第一脉冲序列与第二脉冲序列之间施加有射频脉冲；射频脉冲具有预设的翻转角。

在其中一个实施例中，在射频脉冲之后预设时间内施加有散相梯度。

在其中一个实施例中，根据第一频谱和第二频谱，确定相对弛豫频谱，包括：

确定第一频谱和第二频谱之间的比例关系，根据比例关系确定相对弛豫频谱；

或，计算第一频谱和第二频谱之间的差值，根据差值确定相对弛豫频谱。

在其中一个实施例中，通过分析相对弛豫频谱确定目标成像核素在目标组织的中特征频率，包括：

通过分析相对弛豫频谱，确定不同组织与频率的对应关系；

根据不同组织与频率的对应关系，确定目标成像核素在目标组织中的频率；

根据目标成像核素在目标组织中的频率确定特征频率。

在其中一个实施例中，通过分析相对弛豫频谱，确定不同组织与频率的对应关系，包括：

根据相对弛豫频谱，获取不同频率的幅值强度变化；

根据不同频率幅值强度变化和预设的不同组织的强度变化条件，确定不同组织与频率的对应关系。

在其中一个实施例中，根据目标成像核素在目标组织中的频率确定特征频率，包括：

获取目标成像核素在目标组织中的频谱；

计算频谱中幅值强度的最大值，将幅值强度的最大值对应的频率作为特征频率。

在其中一个实施例中，根据目标成像核素在目标组织中的频率确定的特征频率，包括：

对目标成像核素在目标组织中的频率进行拟合处理，得到特征频率。

第二方面，本申请一个实施例提供一种射频系统中心频率确定装置，包括：

采集模块，用于采集第一磁共振信号和第二磁共振信号，第一磁共振信号是由第一脉冲序列得到，第二磁共振信号由第二脉冲序列得到，第一脉冲序列与第二脉冲序列相邻设置；

第一确定模块，用于确定第一磁共振信号对应的第一频谱，以及第二磁共振信号对应的第二频谱；

第二确定模块，用于根据第一频谱和第二频谱，确定相对弛豫频谱；

第三确定模块，用于通过分析相对弛豫频谱中不同组织的频率确定目标成像核素在目标组织的特征频率，根据所述特征频率确定射频系统的中心频率。

第三方面，本申请一个实施例提供一种计算机设备，包括存储器和处理器，存储器存储由计算机程序，处理器执行计算机程序时实现如上述实施例提供的射频系统中心频率确定方法的步骤。

第四方面，本申请一个实施例提供一种计算机可读存储介质，其上存储由计算机程序，计算机程序被处理器执行时实现如上述实施例提供的射频系统中心频率确定方法的步骤。

第五方面，本申请一个实施例还提供一种计算机程序产品，包括计算机程序，计算机程序被处理器执行时实现如上述实施例提供的射频系统中心频率确定方法的步骤。

本申请实施例提供一种射频系统中心频率确定方法、装置、计算机设备、存储介质和程序产品。该射频系统中心频率确定方法通过采集第一磁共振信号和第二磁共振信号，确定第一磁共振信号对应的第一频谱，以及第二磁共振信号对应的第二频谱；根据第一频谱和第二频谱确定相对弛豫频谱；通过分析相对弛豫频谱确定目标成像核素在目标组织中的特征频率，根据特征频率确定射频系统的中心频率。本申请实施例提供的射频系统中心频率确定方法考虑到弛豫对磁共振信号的影响，即，第二磁共振信号对应的第二频谱会受到弛豫的影响，通过分析第二磁共振信号对应的第二频谱和第一磁共振信号对应的第一频谱之间的差异(相对弛豫频谱)，可以更加准确的查找到目标成像核素在目标组织中的频率，从而能够提高确定目标成像核素在目标组织中的特征频率的准确性，进而利用该特征频率可以对磁共振设备中的射频系统进行准确的频率校准，得到更加准确的射频系统中心频率。

附图说明

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

图1为一个实施例提供的射频系统中心频率确定方法的应用环境图；

图2为一个实施例提供的射频系统中心频率确定方法的步骤流程示意图；

图3为一个实施例提供的第一脉冲序列、射频脉冲、散相梯度、第二脉冲序列、第一磁共振信号和第二磁共振信号的示意图；

图4为另一个实施例提供的第一脉冲序列、射频脉冲、第二脉冲序列、第一磁共振信号和第二磁共振信号的示意图；

图5为另一个实施例提供的射频系统中心频率确定方法的步骤流程示意图；

图6为另一个实施例提供的射频系统中心频率确定方法的步骤流程示意图；

图7为另一个实施例提供的射频系统中心频率确定方法的步骤流程示意图；

图8为一个实施例提供的射频系统中心频率确定装置的结构示意图；

图9为一个实施例提供的计算机设备的结构示意图。

具体实施方式

为使本申请的上述目的、特征和优点能够更加明显易懂，下面结合附图对本申请的具体实施方式做详细的说明。在下面的描述中阐述了很多具体细节以便于充分理解本申请。但是本申请能够以很多不同于在此描述的其它方式来实施，本领域技术人员可以在不违背本申请内涵的情况下做类似改进，因此本申请不受下面公开的具体实施例的限制。

本文中为部件所编序号本身，例如“第一”、“第二”等，仅用于区分所描述的对象，不具有任何顺序或技术含义。

本申请实施例提供的射频系统中心频率确定方法，可以应用于如图1所示的应用环境中，该应用环境包括终端100和磁共振设备200。其中，终端100可以通过网络与磁共振设备200进行通信。终端100可以但不限于是各种个人计算机、笔记本电脑机和平板电脑。本实施例对磁共振设备200的具体种类不作限制。

请参见图2，本申请一个实施例提供一种射频系统中心频率确定方法，以该方法应用于图1中的终端为例进行说明，包括以下步骤：

步骤200、采集第一磁共振信号和第二磁共振信号，第一磁共振信号是由第一脉冲序列得到，第二磁共振信号由第二脉冲序列得到，第一脉冲序列与所述第二脉冲序列相邻设置。

把射频脉冲、梯度场和信号采集时刻等相关各参数的设置及其在时序上的排列称为磁共振成像的脉冲序列。第一脉冲序列和第二脉冲序列可以相同，也可以不同。第一脉冲序列和第二脉冲序列相邻。第一脉冲序列或第二脉冲序列可以是自由感应衰减类序列、自旋回波类序列、反转恢复类序列、梯度回波类序列和平面回波序列中的任意一种。本实施例对第一脉冲序列和第二脉冲序列的种类不作限制，只要能够实现其功能即可。可选地，第一脉冲序列和第二脉冲序列的参数具有相同的设置(中心频率、序列参数和线圈等)，第一脉冲序列和第二脉冲序列相同。

终端采集第一磁共振信号和第二磁共振信号。在磁共振设备中通过第一脉冲序列可以得到待检测对象对应的第一磁共振信号，通过第二脉冲序列可以得到待检测对象对应的第二磁共振信号。第一磁共振信号和第二磁共振信号可以存储在磁共振设备对应的存储设备中，终端在需要时直接在磁共振设备对应的存储设备中获取。第一磁共振信号和第二磁共振信号也可以直接存储在终端的存储器中，在需要时直接在存储器中获取即可。本实施例对采集第一磁共振信号和第二磁共振信号的具体方法不作限制，只要能够实现其功能即可。

步骤210、确定第一磁共振信号对应的第一频谱，以及第二磁共振信号对应的第二频谱。

终端在采集到第一磁共振信号后，对其进行变换，可以得到第一磁共振信号对应的频谱，即第一频谱，第一频谱是指第一磁共振信号中幅值强度随频。率变换的信息。具体地，终端可以对第一磁共振信号进行傅里叶变换得到第一频谱，也可以对第一磁共振信号进行小波变换得到第一频谱。同样的，终端在采集到第二磁共振信号后，对其进行变换，可以得到第二磁共振信号对应的频谱，即第二频谱，第二频谱是指第二磁共振信号中幅值强度随频率变换的信息。具体地，终端可以对第二磁共振信号进行傅里叶变换得到第二频谱，也可以对第二磁共振信号进行小波变换得到第二频谱。本实施例对确定第一频谱和第二频谱的具体方法不作限制，只要能够实现其功能即可。

步骤220、根据第一频谱和第二频谱，确定相对弛豫频谱。

弛豫过程是指在射频脉冲的作用下，原子核发生磁共振达到稳定的高能态后，在射频脉冲消失开始至恢复发生磁共振前的磁矩状态为止的过程。

在磁共振扫描过程中，在第一磁共振信号对应的射频脉冲的作用下，在原子核还没恢复到磁共振前的磁矩状态采集第二磁共振信号，这样得到的磁共振信号会受到弛豫(T1弛豫)的影响，则第二频谱中的幅值强度也会受到弛豫的影响。根据第一频谱和第二频谱可以确定第二频谱中的幅值强度受弛豫影响的信息，即相对弛豫频谱。相对弛豫频谱可以表示在弛豫影响下频率和幅值强度之间的对应关系。本实施例对具体确定相对弛豫频谱的方法不作限制，只要能够实现其功能即可。

步骤230、通过分析相对弛豫频谱确定目标成像核素在目标组织的特征频率，根据特征频率确定射频系统的中心频率。

目标成像核素可以包括氢原子

终端在得到目标成像核素在目标组织中的特征频率后，可以根据该特征频率对射频系统中的初始中心频率进行校准，得到最终的中心频率。本实施例对具体根据特征频率确定射频系统的中心频率的方法不作限制，只要能够实现其功能即可。

本申请实施例提供的射频系统中心频率确定方法通过采集第一磁共振信号和第二磁共振信号；确定第一磁共振信号对应的第一频谱，以及第二磁共振信号对应的第二频谱；根据第一频谱和第二频谱确定相对弛豫频谱；通过分析相对弛豫频谱中确定目标成像核素在目标组织的特征频率，根据特征频率确定射频系统的中心频率。本申请实施例提供的射频系统中心频率确定方法考虑到弛豫对磁共振信号的影响，即，第二磁共振信号对应的频谱会受到弛豫的影响，通过分析第二磁共振信号对应的第二频谱和第一磁共振信号对应的第一频谱之间的差异，可以更加准确的查找到目标成像核素在目标组织中的频率，从而能够提高确定目标成像核素在目标组织中的特征频率的准确性，进而利用该特征频率可以对磁共振设备中的射频系统进行准确的频率校准，得到射频系统的中心频率。另外，与传统技术中采用单一的谱线分析相比，本申请通过分析相对弛豫频谱来确定目标成像核素在目标组织中的特征频率，适用于不同场强、目标成像核素在目标组织中的频率特征不明显、对频率校准结果可靠性要求高的场景中，即，本申请提供的射频系统中心频率确定方法具有较高的适用性。

在一个实施例中，第一脉冲序列与第二脉冲序列之间施加有射频脉冲；射频脉冲具有预设的翻转角。

在根据第一脉冲序列得到第一磁共振信号后，在第二脉冲序列出现前，会向待检测对象施加一个射频脉冲，该射频脉冲具有预设的翻转角。该翻转角可以根据射频特殊吸收率(specific absorption ration，SAR)、采集时间、目标组织的成分、场强等因素确定。在该射频脉冲的激励下，宏观强化强度矢量会偏离静磁场的方向预设的翻转角。具体地，射频脉冲可以是方形脉冲、sinc波形的脉冲、绝热脉冲等，本实施例对射频脉冲的具体形式不作限制，只要能够实现其功能即可。

在本实施例中，通过在第一脉冲序列和第二脉冲序列之间施加具有翻转角的射频脉冲可以提高弛豫对第二频谱中不同组织的幅值强度的影响(第二频谱中不同组织的幅值强度与第一频谱中不同组织的幅值强度相比，具有不同程度的衰减)，从而可以提高根据相对弛豫频谱确定目标成像核素在目标组织的特征频率的准确性，进而提高确定射频系统中心频率的准确性。

在一个可选地实施例中，具有翻转角的射频脉冲可以包括在第一脉冲序列中，也可以包括在第二脉冲序列中。也就是说，若第一脉冲序列或第二脉冲序列中包括射频脉冲的情况下，再施加具有翻转角的射频脉冲，同样可以提高确定射频系统中心频率的准确性。

在一个实施例中，在射频脉冲之后预设时间内施加有散相梯度。

在对待检测对象施加射频脉冲预设时间后，再向待检测对象施加散相梯度，预设时间在施加射频脉冲和第二脉冲序列出现的时间之间。通过在射频脉冲后施加散相梯度可以减少非预期的磁共振信号的影响(例如：可以减少第一脉冲序列产生的自由感应衰减信号强度)，能够提高第二磁共振信号的准确性，从而可以提高确定目标成像核素在目标组织的特征频率，进而提高确定射频系统的中心频率的准确性。本实施例对散相梯度的具体信息不作限制，只要能够实现其功能即可。

在一个可选的实施例中，第一脉冲序列包括方形射频脉冲和第一磁共振信号采集，第二脉冲序列包括方形射频脉冲和第二磁共振信号采集。如图3所示，使用方形射频脉冲产生的自由感应衰减信号为第一磁共振信号和第二磁共振信号。图3中，RF(RadioFrequency，RF)轴上左边的虚线框中是矩形表示第一脉冲序列中的方形射频脉冲，右边的虚线框中的矩形表示第二脉冲序列中的方形射频脉冲，ADC轴上的左边的虚线框中的矩形表示第一磁共振信号采集，右边的虚线框中的矩形表示第二磁共振信号采集；两个虚线框中间的矩形是指施加在第一脉冲序列和第二脉冲序列之间的具有翻转角射频脉冲，两个虚线框中间的梯形是指施加在射频脉冲后的散相梯度。在射频脉冲后施加散相梯度的预设时间小于施加射频脉冲的等效中心与第二脉冲序列的等效中心之间的时间t。图3中的序列图可用于信噪比低时的扫描。

在另一个可选的实施例中，第一脉冲序列包括射频脉冲、梯度和第一磁共振信号采集，第二脉冲序列包括射频脉冲、梯度和第二磁共振信号采集。其中，梯度包括第一方向轴上的梯度G

请参见图5，在一个实施例中，涉及根据第一频谱和第二频谱，确定相对弛豫频谱的一种可能的实现方式，具体步骤包括：

步骤500、确定第一频谱和第二频谱之间的比例关系。

第一频谱和第二频谱之间的比例关系可以是第一频谱除以第二频谱，也可以是第二频谱除以第一频谱。终端在得到第一频谱和第二频谱后，确定第一频谱和第二频谱之间的比例关系。本实施例对确定第一频谱和第二频谱之间的比例关系的具体方法不作限制，只要能够实现其功能即可。

步骤510、根据比例关系确定相对弛豫频谱。

终端在得到第一频谱和第二频谱之间的比例关系后，可以根据该比例关系确定相对弛豫频谱。具体地，第一频谱和第二频谱之间的比值即为相对弛豫频谱。也就是说，使用第一频谱和第二频谱之间的比例关系表征第二频谱中的幅值强度受弛豫影响的信息。

在一个实施例中，涉及根据第一频谱和第二频谱，确定相对弛豫频谱的一种可能的实现方式，具体步骤包括：

计算第一频谱和第二频谱之间的差值，根据差值确定相对弛豫频谱。

终端在得到第一频谱和第二频谱后，计算第一频谱和第二频谱之间的差值，该差值可以是第一频谱减去第二频谱后的差值，也可以是第二频谱减去第一频谱后的差值。本实施例对差值的具体计算方法不作限制，只要能够实现其功能即可。终端在得到差值后，将该差值作为相对弛豫频谱。也就是说，使用第一频谱和第二频谱之间的差值表征第二频谱中的幅值强度受弛豫影响的信息。

在本实施例中，提供了两种确定相对弛豫频谱的方法，且每种方法简单、容易实现，使用者可以根据实际需求自行选择，使得本实施例提供的射频系统中心频率确定方法具有较高的实用性。

请参见图6，在一个实施例中，涉及通过分析相对弛豫频谱确定目标成像核素在目标组织的特征频率，包括：

步骤600、通过分析相对弛豫频谱，确定不同组织与频率的对应关系。

相对弛豫频谱中包括目标成像核素在不同组织中的频谱，终端通过对相对弛豫频谱进行分析，可以确定相对弛豫频谱中不同组织与频率的对应关系。也就是说，终端通过对相对弛豫频谱进行分析，可以得到各个组织对应的频率。

步骤610、根据不同组织与频率的对应关系，确定目标成像核素在目标组织中的频率。

步骤620、根据目标成像核素在目标组织中的频率确定特征频率。

终端在得到不同组织与频率的对应关系后，在不同组织与频率的对应关系中，查找目标组织对应的频率。具体地，终端可以先在各种组织中查找到目标组织，再根据不同组织与频率的对应关系，获取目标组织的频率，即目标成像核素在目标组织中的频率。

终端在得到目标成像核素在目标组织中的频率，根据该频率确定目标成像核素在目标组织中的特征频率。本实施例对确定特征频率的方法不作限制，只要能够实现其功能即可。

在一个实施例中，如图7所示，涉及通过分析相对弛豫频谱，确定不同组织与频率的对应关系的一种可能的实现方式，步骤包括：

步骤700、根据相对弛豫频谱，获取不同频率的幅值强度变化。

相对弛豫频谱表示的是频率和幅值强度之间的对应关系，不同频率的幅值强度受弛豫影响会发生变化。终端通过对相对弛豫频谱进行分析，可以获取相对弛豫频谱中不同频率的幅值强度变化。本实施例对获取不同频率的幅值强度变化的具体方法不作限制，只要能够实现其功能即可。

步骤710、根据不同频率的幅值强度变化和预设的不同组织的强度变化条件，确定不同组织与频率的对应关系。

不同组织的频率对应的幅值强度的变化不同，预设的不同组织的强度变化条件对应于不同组织。终端通过将不同频率的幅值强度的变化分别与预设的不同组织的强度变化条件进行对比，将幅值强度变化满足预设的不同组织的强度变化条件的幅值强度变化对应频率确定为预设组织对应的频率，预设组织为预设的不同组织的强度变化条件对应的组织，从而得到不同组织与频率的对应关系。预设的不同组织的强度变化条件可以是指不同组织的预设强度阈值，若幅值强度变化大于或等于预设强度阈值，表示幅值强度变化满足预设强度变化条件，若幅值强度变化小于预设强度阈值，表示幅值强度变化不满足预设强度变化条件。具体地，不同的组织包括脑白质和脑灰质，假设脑白质对应的强度变化条件为第一强度变化条件，脑灰质对应的强度变化条件为第二强度变化条件，则终端将获取到的幅值强度变化与第一强度变化条件进行对比，将幅值强度变化满足第一强度变化条件的幅值强度变化对应的频率确定为脑白质对应的频率；终端将获取到的幅值强度变化与第二强度变化条件进行对比，将幅值强度变化满足第二强度变化条件的幅值强度变化对应的频率确定为脑灰质对应的频率，从而可以得到脑白质和与其对应的频率之间的对应关系，脑灰质和与其对应的频率之间的对应关系。

本实施例提供的确定目标成像核素在目标组织中的频率的方法快捷易懂，且容易实现。

在一个具体的实施例中，在第二脉冲序列如图3所示，即第二脉冲序列中仅包括一个射频脉冲时，当施加的射频脉冲的等效中心与第二脉冲序列中的射频脉冲的等效中心之间的时间t满足目标组织的纵向磁化矢量刚好恢复到零，即目标组织的磁共振信号在第二频谱中的占比较小。将施加的射频脉冲记为目标射频脉冲。根据反转恢复过程中信号强度随时间变化的公式：

在一个可选的实施例中，涉及通过分析相对弛豫频谱确定目标成像核素在目标组织的特征频率的一种可能的实现方式，包括：

终端可以基于扩展相位图(extended phase graph,EPG)算法，利用目标组织的弛豫时间常数、第一脉冲序列中的射频脉冲的翻转角、第二脉冲序列中的射频脉冲的翻转角、在第一脉冲序列和第二脉冲序列之间施加的射频脉冲的翻转角、以及该射频脉冲的等效中心和第二脉冲序列中的射频脉冲的等效中心之间的时间t计算得到相对弛豫频谱理论值。根据相对弛豫频谱与相对弛豫频谱理论值之间的差异(权重)，确定目标成像核素在目标组织中的特征频率。

具体地，假设目标组织为肌肉，计算得到的相对弛豫频谱理论值为0.26，相对弛豫频谱中各频率对应的幅值强度为0.1，0.05，0.11，0.19，0.23，0.25，0.27，0.24，0.20，0.14，0.11。相对弛豫频谱与相对弛豫频谱理论值之间的差异可以用相对弛豫频谱与相对弛豫频谱理论值之间的比值表示，即，相对弛豫频谱除以相对弛豫频谱理论值后得到的比值，当该比值大于1或小于2时，取2与该比值的差值；当该比值大于2时取0。根据相对弛豫频谱和相对弛豫频谱理论值，可以得到相对弛豫频谱中幅值强度对应的权重0.38，0.19，0.42，0.73，0.88，0.96，0.96，0.92，0.77，0.54，0.42。应用权重和频谱分析方法可以得到目标成像核素在目标组织中的特征频率。具体地，可以权重处于预设范围内的幅值强度对应的频率确定为目标成像核素在目标组织中的频率，然后根据目标成像核素在目标组织中的频率确定特征频率。

在另一个实施例中，相对弛豫频谱和相对弛豫频谱理论值之间的差异(权重)的计算方法还包括：取差值的1范数，即1-|相对弛豫频谱的幅值强度-相对弛豫频谱理论值|，将小于1的值取为零；取差值的2范数，即1-(相对弛豫频谱的幅值强度-相对弛豫频谱理论值)

在一个实施例中，涉及根据目标成像核素在目标组织中的频率确定特征频率的一种可能的实现方式，步骤包括：

获取目标成像核素在目标组织中的频谱；计算频谱中的幅值强度的最大值，将幅值强度的最大值对应的频率作为中心频率。

终端在从相对弛豫频谱中确定目标成像核素在目标组织中的频率后，获取该频率对应的幅值强度；根据该频率和该频率对应的幅值强度，得到目标成像核素在目标组织中的频谱。终端在得到目标成像核素在目标组织中的频谱后，将该频谱中的幅值强度进行对比，找出幅值强度中的最大值，将该幅值强度最大值对应的频率作为目标成像核素在目标组织中的特征频率。

在一个实施例中，涉及根据目标成像核素在目标组织中的频率确定特征频率的一种可能的实现方法，步骤包括：

对目标成像核素在目标组织中的频率进行拟合处理，得到特征频率。

终端在得到目标成像核素在目标组织中的频率后，将所有频率根据拟合算法进行拟合处理，可以得到目标成像核素在目标组织中的特征频率。具体地，拟合算法包括最小二乘法和化矩阵求伪逆等算法。本实施例对目标成像核素在目标组织中的频率进行拟合处理的具体过程不作限制，只要能够实现其功能即可。

应该理解的是，虽然图中的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，图中的至少一部分步骤可以包括多个子步骤或者多个阶段，这些子步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些子步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。

基于同样的发明构思，本申请实施例还提供了一种用于实现上述所涉及的射频系统中心频率确定方法的射频系统中心频率确定装置。该装置所提供的解决问题的实现方案与上述方法中所记载的实现方案相似，故下面所提供的一个或多个射频系统中心频率确定装置实施例中的具体限定可以参见上文中对于射频系统中心频率确定方法的限定，在此不再赘述。

请参见图8，本申请一个实施例提供一种射频系统中心频率确定装置10,该装置包括采集模块11、变换模块12、第一确定模块13和第二确定模块14。其中，

采集模块11用于采集第一磁共振信号和第二磁共振信号，第一磁共振信号是由第一脉冲序列得到，第二磁共振信号由第二脉冲序列得到，第一脉冲序列与第二脉冲序列相邻设置；

第一确定模块12用于确定第一磁共振信号对应的第一频谱，以及第二磁共振信号对应的第二频谱；

第二确定模块13用于根据第一频谱和第二频谱，确定相对弛豫频谱；

第三确定模块14用于通过分析相对弛豫频谱确定目标成像核素在目标组织的特征频率，根据特征频率确定射频系统的中心频率。

在一个实施例中，第一脉冲序列与第二脉冲序列之间施加有射频脉冲；射频脉冲具有预设的翻转角。

在一个实施例中，在射频脉冲之后预设时间内施加有散相梯度。

在一个实施例中，第二确定模块13具体用于确定第一频谱和第二频谱之间的比例关系；根据比例关系确定相对弛豫频谱。

在一个实施例中，第二确定模块13具体还用于计算第一频谱和第二频谱之间的差值，根据差值确定相对弛豫频谱。

在一个实施例中，第三确定模块14包括第一确定单元、第二确定单元和第三确定单元，第一确定单元用于通过分析相对弛豫频谱，确定不同组织与频率的对应关系；第二确定单元用于根据不同组织与频率的对应关系，确定目标成像核素在目标组织中的频率；第三确定单元用于根据目标成像核素在目标组织中的频率确定特征频率。

在一个实施例中，第二确定单元具体用于根据相对弛豫频谱，获取不同频率的幅值强度变化；根据不同频率的幅值强度变化和预设的不同组织的强度变化条件，确定不同组织与频率的对应关系。

在一个实施例中，第三确定单元用于获取目标成像核素在目标组织中的频谱；计算频谱中幅值强度的最大值，将幅值强度的最大值对应的频率作为特征频率。

在一个实施例中，第三确定单元还用于对目标成像核素在目标组织中的频率进行拟合处理，得到特征频率。

上述射频系统中心频率确定装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。

在一个实施例中，提供了一种计算机设备，该计算机设备可以是终端，其内部结构图可以如图9所示。该计算机设备包括通过系统总线连接的处理器、存储器、通信接口、显示屏和输入装置。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统和计算机程序。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的通信接口用于与外部的终端进行有线或无线方式的通信，无线方式可通过WIFI、移动蜂窝网络、NFC(近场通信)或其他技术实现。该计算机程序被处理器执行时以实现一种射频系统中心频率确定方法。该计算机设备的显示屏可以是液晶显示屏或者电子墨水显示屏，该计算机设备的输入装置可以是显示屏上覆盖的触摸层，也可以是计算机设备外壳上设置的按键、轨迹球或触控板，还可以是外接的键盘、触控板或鼠标等。

本领域技术人员可以理解，图9中示出的结构，仅仅是与本申请方案相关的部分结构的框图，并不构成对本申请方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。

在一个实施例中，提供了一种计算机设备，包括存储器和处理器，存储器中存储有计算机程序，处理器执行计算机程序时实现以下步骤：

采集第一磁共振信号和第二磁共振信号，第一磁共振信号是由第一脉冲序列得到，第二磁共振信号由第二脉冲序列得到，第一脉冲序列与第二脉冲序列相邻设置；

确定第一磁共振信号对应的第一频谱，以及第二磁共振信号对应的第二频谱；

根据第一频谱和第二频谱，确定相对弛豫频谱；

通过分析相对弛豫频谱确定目标成像核素在目标组织的特征频率，根据特征频率确定射频系统的中心频率。

在一个实施例中，第一脉冲序列与第二脉冲序列之间施加有射频脉冲；射频脉冲具有预设的翻转角。

在一个实施例中，在射频脉冲之后预设时间内施加有散相梯度。

在一个实施例中，处理器执行计算机程序时还实现以下步骤：确定第一频谱和第二频谱之间的比例关系；根据比例关系确定相对弛豫频谱。

在一个实施例中，处理器执行计算机程序时还实现以下步骤：计算第一频谱和第二频谱之间的差值，根据差值确定相对弛豫频谱。

在一个实施例中，处理器执行计算机程序时还实现以下步骤：通过分析相对弛豫频谱确定不同组织与频率的对应关系；根据不同组织与频率的对应关系，确定目标成像核素在目标组织中的频率；根据目标成像核素在目标组织中的频率确定特征频率。

在一个实施例中，处理器执行计算机程序时还实现以下步骤：根据相对弛豫频谱，获取不同频率的幅值强度变化；根据不同频率的幅值强度变化和预设的不同组织的强度变化条件，确定不同组织与频率的对应关系。

在一个实施例中，处理器执行计算机程序时还实现以下步骤：获取目标成像核素在目标组织中的频谱；计算频谱中幅值强度的最大值，将幅值强度的最大值对应的频率作为特征频率。

在一个实施例中，处理器执行计算机程序时还实现以下步骤：对目标成像核素在目标组织中的频率进行拟合处理，得到特征频率。

在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：

采集第一磁共振信号和第二磁共振信号，第一磁共振信号是由第一脉冲序列得到，第二磁共振信号由第二脉冲序列得到，第一脉冲序列与第二脉冲序列相邻设置；

确定第一磁共振信号对应的第一频谱，以及第二磁共振信号对应的第二频谱；

根据第一频谱和第二频谱，确定相对弛豫频谱；

通过分析相对弛豫频谱确定目标成像核素在目标组织的特征频率，根据特征频率确定射频系统的中心频率。

在一个实施例中，第一脉冲序列与第二脉冲序列之间施加有射频脉冲；射频脉冲具有预设的翻转角。

在一个实施例中，在射频脉冲之后预设时间内施加有散相梯度。

在一个实施例中，计算机程序被处理器执行时还实现以下步骤：确定第一频谱和第二频谱之间的比例关系；根据比例关系确定相对弛豫频谱。

在一个实施例中，计算机程序被处理器执行时还实现以下步骤：计算第一频谱和第二频谱之间的差值，根据差值确定相对弛豫频谱。

在一个实施例中，计算机程序被处理器执行时还实现以下步骤：通过分析相对弛豫频谱确定不同组织与频率的对应关系；根据不同组织与频率的对应关系，确定目标成像核素在目标组织中的频率；根据目标成像核素在目标组织中的频率确定特征频率。

在一个实施例中，计算机程序被处理器执行时还实现以下步骤：根据相对弛豫频谱，获取不同频率的幅值强度变化；根据不同频率的幅值强度变化和预设的不同组织的强度变化条件，确定不同组织与频率的对应关系。

在一个实施例中，计算机程序被处理器执行时还实现以下步骤：获取目标成像核素在目标组织中的频谱；计算频谱中幅值强度的最大值，将幅值强度的最大值对应的频率作为特征频率。

在一个实施例中，计算机程序被处理器执行时还实现以下步骤：对目标成像核素在目标组织中的频率进行拟合处理，得到特征频率。

在一个实施例中，提供了一种计算机程序产品，包括计算机程序，该计算机程序被处理器执行时实现以下步骤：

采集第一磁共振信号和第二磁共振信号，第一磁共振信号是由第一脉冲序列得到，第二磁共振信号由第二脉冲序列得到，第一脉冲序列与第二脉冲序列相邻设置；

确定第一磁共振信号对应的第一频谱，以及第二磁共振信号对应的第二频谱；

根据第一频谱和第二频谱，确定相对弛豫频谱；

通过分析相对弛豫频谱确定目标成像核素在目标组织的特征频率，根据特征频率确定射频系统的中心频率。

在一个实施例中，第一脉冲序列与第二脉冲序列之间施加有射频脉冲；射频脉冲具有预设的翻转角。

在一个实施例中，在射频脉冲之后预设时间内施加有散相梯度。

在一个实施例中，计算机程序被处理器执行时还实现以下步骤：确定第一频谱和第二频谱之间的比例关系；根据比例关系确定相对弛豫频谱。

在一个实施例中，计算机程序被处理器执行时还实现以下步骤：计算第一频谱和第二频谱之间的差值，根据差值确定相对弛豫频谱。

在一个实施例中，计算机程序被处理器执行时还实现以下步骤：通过分析相对弛豫频谱确定不同组织与频率的对应关系；根据不同组织与频率的对应关系，确定目标成像核素在目标组织中的频率；根据目标成像核素在目标组织中的频率确定特征频率。

在一个实施例中，计算机程序被处理器执行时还实现以下步骤：根据相对弛豫频谱，获取不同频率的幅值强度变化；根据不同频率的幅值强度变化和预设的不同组织的强度变化条件，确定不同组织与频率的对应关系。

在一个实施例中，计算机程序被处理器执行时还实现以下步骤：获取目标成像核素在目标组织中的频谱；计算频谱中幅值强度的最大值，将幅值强度的最大值对应的频率作为中心频率。

在一个实施例中，计算机程序被处理器执行时还实现以下步骤：对目标成像核素在目标组织中的频率进行拟合处理，得到特征频率。

本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限，RAM以多种形式可得，诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。

以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。

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