磁共振系统频率漂移量的确定方法、装置和计算机设备与流程

1.本技术涉及磁共振技术领域，特别是涉及一种磁共振系统频率漂移量的确定方法、装置和计算机设备。


背景技术：

2.在磁共振成像过程中，系统频率的准确性直接影响图像质量。在进行长时间的序列扫描时，磁共振系统会因为电流不断切换引起的振动而发热，进而导致系统频率发生漂移。
3.相关技术中，为了避免系统频率漂移降低图像质量，通常会对系统频率进行实时校准。但是，在扫描乳腺或者颈部的时候，常常需要用户手动对系统频率进行调整。这时，如果仍然按照之前的校准方式进行校准，会将系统频率的手动调整量当成频率漂移，从而导致系统频率校准错误。


技术实现要素：

4.基于此，有必要针对上述技术问题，提供一种能够避免将系统频率的手动调整量当成频率漂移，从而避免系统频率校准错误的磁共振系统频率漂移量的确定方法、装置和计算机设备。
5.一种磁共振系统频率漂移量的确定方法，该方法包括：
6.采用预设扫描序列进行第一扫描以获取第一组回波信号；预设扫描序列为梯度双回波序列；
7.在第一扫描的预设时长后，采用预设扫描序列进行第二扫描以获取第二组回波信号；
8.分别根据第一组回波信号和第二组回波计算，得到第一导航频率和第二导航频率；
9.根据第一导航频率和第二导航频率确定磁共振系统频率漂移量。
10.在其中一个实施例中，上述每组回波信号均包括第一回波信号和第二回波信号，分别根据第一组回波信号和第二组回波计算，得到第一导航频率和第二导航频率，包括：
11.针对每一组回波信号，分别对第一回波信号和第二回波信号进行频域转换，得到转换后的第三回波信号和第四回波信号；
12.根据第三回波信号和第四回波信号计算导航频率。
13.在其中一个实施例中，上述根据第三回波信号和第四回波信号计算导航频率，包括：
14.将第三回波信号与第四回波信号进行整合得到系统回波信号；系统回波信号采用振幅和相位角表示；
15.对系统回波信号进行相位拟合，得到第一相位角；
16.根据频率与相位角之间的对应关系和第一相位角计算导航频率。
17.在其中一个实施例中，上述根据第三回波信号和第四回波信号计算导航频率，包括：
18.根据第三回波信号、第四回波信号和预先设置的相位角函数进行计算，得到第二相位角；
19.根据频率与相位角之间的对应关系和第二相位角计算导航频率。
20.在其中一个实施例中，上述频率与相位角之间的对应关系包括频率与相位角成正比，与第一回波信号和第二回波信号之间的时间差成反比。
21.在其中一个实施例中，第一组回波信号包括多对回波信号，每对回波信号包括两个回波信号，上述分别根据第一组回波信号和第二组回波计算，得到第一导航频率和第二导航频率，包括：
22.针对第一组回波信号，根据多对回波信号计算导航频率，得到多个候选导航频率；
23.将多个候选导航频率的平均值确定为第一导航频率。
24.在其中一个实施例中，上述根据第一导航频率和第二导航频率确定磁共振系统频率漂移量，包括：
25.计算第一导航频率和第二导航频率之间的频率差值；
26.对频率差值进行解缠绕处理，得到磁共振系统频率漂移量。
27.在其中一个实施例中，在上述根据第一导航频率和第二导航频率确定磁共振系统频率漂移量之后，该方法还包括：
28.根据磁共振系统频率漂移量修正磁共振系统频率。
29.在其中一个实施例中，在上述根据第一导航频率和第二导航频率确定磁共振系统频率漂移量之后，该方法还包括：
30.根据磁共振系统频率漂移量修正发射线圈的频率和/或接收线圈的频率。
31.一种磁共振系统频率漂移量的确定装置，该装置包括：
32.第一序列扫描模块，用于采用预设扫描序列进行第一扫描以获取第一组回波信号；预设扫描序列为梯度双回波序列；
33.第二序列扫描模块，用于在第一扫描的预设时长后，采用预设扫描序列进行第二扫描以获取第二组回波信号；
34.导航频率计算模块，用于分别根据第一组回波信号和第二组回波计算，得到第一导航频率和第二导航频率；
35.漂移量确定模块，用于根据第一导航频率和第二导航频率确定磁共振系统频率漂移量。
36.在其中一个实施例中，每组回波信号均包括第一回波信号和第二回波信号，上述导航频率计算模块包括：
37.频域转换子模块，用于针对每一组回波信号，分别对第一回波信号和第二回波信号进行频域转换，得到转换后的第三回波信号和第四回波信号；
38.导航频率计算子模块，用于根据第三回波信号和第四回波信号计算导航频率。
39.在其中一个实施例中，上述导航频率计算子模块，用于将第三回波信号与第四回波信号进行整合得到系统回波信号；系统回波信号采用振幅和相位角表示；对系统回波信号进行相位拟合，得到第一相位角；根据频率与相位角之间的对应关系和第一相位角计算
导航频率。
40.在其中一个实施例中，上述导航频率计算子模块，用于根据第三回波信号、第四回波信号和预先设置的相位角函数进行计算，得到第二相位角；根据频率与相位角之间的对应关系和第二相位角计算导航频率。
41.在其中一个实施例中，上述频率与相位角之间的对应关系包括频率与相位角成正比，与第一回波信号和第二回波信号之间的时间差成反比。
42.在其中一个实施例中，上述第一组回波信号包括多对回波信号，每对回波信号包括两个回波信号，上述第一导航频率计算模块，具体用于根据多对回波信号计算导航频率，得到多个候选导航频率；将多个候选导航频率的平均值确定为第一导航频率。
43.在其中一个实施例中，上述漂移量确定模块，具体用于计算第一导航频率和第二导航频率之间的频率差值；对频率差值进行解缠绕处理，得到磁共振系统频率漂移量。
44.在其中一个实施例中，该装置还包括：
45.导航频率修正模块，用于根据磁共振系统频率漂移量修正磁共振系统频率。
46.在其中一个实施例中，该装置还包括：
47.线圈频率修改模块，用于根据磁共振系统频率漂移量修正发射线圈的频率和/或接收线圈的频率。
48.一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现以下步骤：
49.采用预设扫描序列进行第一扫描以获取第一组回波信号；预设扫描序列为梯度双回波序列；
50.在第一扫描的预设时长后，采用预设扫描序列进行第二扫描以获取第二组回波信号；
51.分别根据第一组回波信号和第二组回波计算，得到第一导航频率和第二导航频率；
52.根据第一导航频率和第二导航频率确定磁共振系统频率漂移量。
53.一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现以下步骤：
54.采用预设扫描序列进行第一扫描以获取第一组回波信号；预设扫描序列为梯度双回波序列；
55.在第一扫描的预设时长后，采用预设扫描序列进行第二扫描以获取第二组回波信号；
56.分别根据第一组回波信号和第二组回波计算，得到第一导航频率和第二导航频率；
57.根据第一导航频率和第二导航频率确定磁共振系统频率漂移量。
58.上述磁共振系统频率漂移量的确定方法、装置和计算机设备，采用预设扫描序列进行第一扫描以获取第一组回波信号；在第一扫描的预设时长后，采用预设扫描序列进行第二扫描以获取第二组回波信号；分别根据第一组回波信号和第二组回波计算，得到第一导航频率和第二导航频率；根据第一导航频率和第二导航频率确定磁共振系统频率漂移量。通过本技术实施例，在用户手动调整系统频率之后，采用预设扫描序列进行第一扫描并
获取第一组回波信号，经过预设时长再采用预设扫描序列进行第二扫描并获取第二组回波信号，分别根据第一组回波信号和第二组回波信号计算第一导航频率和第二导航频率，由于第一导航频率和第二导航频率中均包含系统频率的手动调整量，因此可以通过计算将系统频率的手动调整量抵消掉，从而避免将系统频率的手动调整量当成频率漂移而导致系统频率校准错误的问题。
附图说明
59.图1为一个实施例中磁共振系统频率漂移量的确定方法的应用环境图；
60.图2为一个实施例中磁共振系统频率漂移量的确定方法的流程示意图；
61.图3为一个实施例中的波形图；
62.图4为一个实施例中根据回波信号计算导航频率步骤的流程示意图；
63.图5为另一个实施例中磁共振系统频率漂移量的确定方法的流程示意图；
64.图6为采用现有技术方法获得的多期相的乳腺磁共振图像；
65.图7为一个实施例中对磁共振系统频率修正后获得的多期相的乳腺磁共振图像；
66.图8为一个实施例中磁共振系统频率漂移量的确定装置的结构框图；
67.图9为一个实施例中计算机设备的内部结构图。
具体实施方式
68.为了使本技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本技术进行进一步详细说明。应当理解，此处描述的具体实施例仅仅用以解释本技术，并不用于限定本技术。
69.本技术提供的磁共振系统频率漂移量的确定方法，可以应用于如图1所示的应用环境中。该应用环境为磁共振系统，该磁共振系统至少包括梯度系统101、射频系统102和控制终端103。其中，梯度系统101包括梯度线圈等设备，用于产生梯度脉冲；射频系统102包括射频线圈等设备，射频线圈包括发射线圈和接收线圈，射频系统102用于产生射频脉冲信号以及检测对象核自旋产生的磁共振信号；控制终端103与梯度系统101和射频系统102通过网络进行通信，用于控制梯度系统101和射频系统102。上述控制终端103可以但不限于是各种个人计算机、笔记本电脑和平板电脑。
70.在一个实施例中，如图2所示，提供了一种磁共振系统频率漂移量的确定方法，以该方法应用于图1中的控制终端为例进行说明，包括以下步骤：
71.步骤201，采用预设扫描序列进行第一扫描以获取第一组回波信号；预设扫描序列为梯度双回波序列。
72.由于在扫描乳腺或者颈部的时候，常常需要用户手动对系统频率进行调整，因此采用实时校准的方式，会将系统频率的手动调整量当成频率漂移，从而导致系统频率校准错误。为了避免这种情况，本技术实施例在检测到针对系统频率的调整操作，即用户手动调整系统频率之后，不再进行实时频率校准，而是采用动态频率校准。具体地，先采用预设扫描序列进行第一扫描，其中预设扫描序列可以是梯度双回波序列，见图3所示的波形图，rf为射频脉冲，g
ss
为选层方向梯度，g
ro
为频率编码方向梯度，adc表示模数转换器接收的信号。图3所示的该梯度双回波序列只有频率编码梯度，没有相位编码梯度。并且，为了降低预
设扫描序列对磁场的影响，在满足信噪比情况下，激发翻转角尽量小，各个逻辑轴的梯度净面积为零。
73.在采用预设扫描序列进行第一扫描的同时，采集第一组回波信号，见图3所示的波形图，沿g
ss
方向与射频脉冲相对应的施加选层梯度和该选层梯度后紧邻施加的散相梯度；沿g
ro
方向施加连续频率编码梯度，以避免梯度的频繁切换引起梯度涡流；adc采样得到第一回波信号e0和第二回波信号e1，该两个回波所对应的g
ro
方向的梯度矩相同；回波时间te(echo time)是射频脉冲中点到第一回波信号e0的时间间隔，δte为第一回波信号e0与第二回波信号e1的回波时间差值。在此实施例中，预设扫描序列不包含相位编码方向的梯度，也即本技术实施例中采样得到第一回波信号e0和第二回波信号e1是未经过相位编码处理的。在一实施例中，te的值满足保证检测对象的水信号和脂肪信号处于同相位。在其中一个实施例中，在用户手动调整系统频率之前，先进行一次系统频率的校准，以便在确定磁共振系统频率漂移量之后，可以更加准确地修正系统频率。其中，系统频率为拉莫尔频率，即磁性原子核的进动频率或主磁场频率。
74.步骤202，在第一扫描的预设时长后，采用预设扫描序列进行第二扫描以获取第二组回波信号。
75.本技术实施例中，在获取第一组回波信号之后，系统频率经过一段时间可能发生了漂移，因此采用预设扫描序列进行第二扫描，同时获取第二组回波信号。
76.例如，采用预设扫描序列进行一次扫描，5分钟后采用预设扫描序列再进行一次扫描，并获取第二组回波信号。本技术实施例对预设时长不作详细限定，可以根据实际情况进行设置。
77.步骤203，分别根据第一组回波信号和第二组回波计算，得到第一导航频率和第二导航频率。
78.本技术实施例中，根据第一组回波信号计算得到第一导航频率，根据第二组回波信号计算得到第二导航频率。具体地，对每组回波信号进行频域转换，然后根据转换后的回波信号计算相位角，再根据相位角和频率之间的对应关系计算导航频率。本技术实施例对具体计算方式不作详细限定，可以根据实际情况进行设置。
79.步骤204，根据第一导航频率和第二导航频率确定磁共振系统频率漂移量。
80.本技术实施例中，由于第一导航频率中包含系统频率的手动调整量，第二导航频率中也包含系统频率的手动调整量，因此在得到第一导航频率和第二导航频率之后，可以计算第一导航频率和第二导航频率之间的频率差值，从而将系统频率的手动调整量抵消掉，进而得到磁共振系统频率漂移量。
81.上述磁共振系统频率漂移量的确定方法中，采用预设扫描序列进行第一扫描以获取第一组回波信号；在第一扫描的预设时长后，采用预设扫描序列进行第二扫描以获取第二组回波信号；分别根据第一组回波信号和第二组回波计算，得到第一导航频率和第二导航频率；根据第一导航频率和第二导航频率确定磁共振系统频率漂移量。本技术实施例中，第一导航频率和第二导航频率中均包含系统频率的手动调整量，因此可以通过计算将系统频率的手动调整量抵消掉，从而避免将系统频率的手动调整量当成频率漂移而导致系统频率校准错误的问题。
82.在一个实施例中，如图4所示，涉及根据回波信号计算导航频率的一种可选过程。
在上述实施例的基础上，每组回波信号均包括第一回波信号和第二回波信号，因此计算第一导航频率和第二导航频率，均可以采用如下步骤：
83.步骤301，针对每一组回波信号，分别对第一回波信号和第二回波信号进行频域转换，得到转换后的第三回波信号和第四回波信号。
84.本技术实施例中，采用梯度双回波序列进行扫描，可以采集到两个回波信号，即第一回波信号和第二回波信号。之后，对第一回波信号进行频域转换得到第三回波信号，对第二回波信号进行频域转换得到第四回波信号。具体地，频域转换可以采用傅里叶变换，如公式(1)(2)：
85.e0＝fft(e0)
-------------------------------------------
(1)
86.e1＝fft(e1)
--------------------------------------------
(2)
87.其中，e0为第一回波信号，e1为第二回波信号，e0为第三回波信号，e1为第四回波信号。
88.步骤302，根据第三回波信号和第四回波信号计算导航频率。
89.本技术实施例中，得到频域下的第三回波信号和第四回波信号之后，可以采用以下两种方式计算导航频率。
90.方式一：将第三回波信号与第四回波信号进行整合得到系统回波信号；系统回波信号采用振幅和相位角表示；对系统回波信号进行相位拟合，得到第一相位角；根据频率与相位角之间的对应关系和第一相位角计算导航频率。
91.具体地，先将两个回波信号整合为一个系统回波信号，其中系统回波信号采用振幅和相位角表示。整合过程可以如公式(3)：
[0092][0093]
然后将e(i)记作公式(4)：
[0094][0095]
其中，e0(i)*e1(i)为第三回波信号e0和第四回波信号e1共轭相乘；a
i
为系统回波信号e(i)的幅值，为系统回波信号e(i)的相位角，n表示所采集得到的梯度双回波的数量，n≥1。
[0096]
在得到系统回波信号之后，可以采用最小二乘法进行相位拟合，如公式(5)：
[0097][0098]
其中，θ1为拟合得到的第一相位角。
[0099]
之后，根据频率与相位角之间的对应关系，以及拟合得到的第一相位角计算导航频率，具体如公式(6)：
[0100]
[0101]
上述频率与相位角之间的对应关系包括：频率与相位角成正比，与第一回波信号和第二回波信号之间的时间差成反比。其中，f为计算出的导航频率，δte为第一回波信号和第二回波信号之间的时间差。
[0102]
方式二：根据第三回波信号、第四回波信号和预先设置的相位角函数进行计算，得到第二相位角；根据频率与相位角之间的对应关系和第二相位角计算导航频率。
[0103]
具体地，可以直接根据相位角函数进行计算，得到第二相位角，如公式(7)：
[0104]
θ2＝angle(∑
i
(e0(i)*e1(i)))
-----------------------
(7)
[0105]
其中，θ1为第二相位角。
[0106]
之后，根据频率与相位角之间的对应关系，以及计算得到的第二相位角计算导航频率，具体如公式(8)：
[0107][0108]
上述频率与相位角之间的对应关系包括：频率与相位角成正比，与第一回波信号和第二回波信号之间的时间差成反比。其中，f为计算出的导航频率，δte为第一回波信号和第二回波信号之间的时间差。
[0109]
上述根据回波信号计算导航频率的过程中，针对每组回波信号，分别对第一回波信号和第二回波信号进行频域转换，得到转换后的第三回波信号和第四回波信号；根据第三回波信号和第四回波信号计算导航频率。通过本技术实施例，先将时域表示的回波信号转换为频域表示的回波信号，由频域表示的回波信号可以计算出相位角，再根据频率与相位角之间的对应关系可以计算出导航频率，从而确定第一导航频率和第二导航频率，进而得到磁共振系统频率漂移量。
[0110]
在一个实施例中，如图5所示，提供了一种磁共振系统频率漂移量的确定方法，以该方法应用于图1中的控制终端为例进行说明，包括以下步骤：
[0111]
步骤401，采用预设扫描序列进行第一扫描以获取第一组回波信号；预设扫描序列为梯度双回波序列。
[0112]
步骤402，在第一扫描的预设时长后，采用预设扫描序列进行第二扫描以获取第二组回波信号。
[0113]
步骤403，分别根据第一组回波信号和第二组回波计算，得到第一导航频率和第二导航频率。
[0114]
本技术实施例中，为了使第一导航频率更加准确，针对第一组回波信号，可以根据多对回波信号计算导航频率，得到多个候选导航频率；将多个候选导航频率的平均值确定为第一导航频率。具体地，采集多对回波信号，每对回波信号均包括两个回波信号；针对每对回波信号进行导航频率计算，得到对应的候选导航频率，这样就可以得到多个候选导航频率；最后，计算多个候选导航频率的平均值，将平均值作为第一导航频率。
[0115]
在其中一个实施例中，针对每对回波信号计算导航频率，得到对应的候选导航频率可以包括：分别对第一回波信号和第二回波信号进行频域转换，得到转换后的第三回波信号和第四回波信号；根据第三回波信号和第四回波信号计算导航频率。
[0116]
在其中一个实施例中，根据第三回波信号和第四回波信号计算导航频率，包括：将第三回波信号与第四回波信号进行整合得到系统回波信号；系统回波信号采用振幅和相位
角表示；对系统回波信号进行相位拟合，得到第一相位角；根据频率与相位角之间的对应关系和第一相位角计算导航频率。
[0117]
在其中一个实施例中，根据第三回波信号和第四回波信号计算导航频率，包括：根据第三回波信号、第四回波信号和预先设置的相位角函数进行计算，得到第二相位角；根据频率与相位角之间的对应关系和第二相位角计算导航频率。
[0118]
在其中一个实施例中，频率与相位角之间的对应关系包括频率与相位角成正比，与第一回波信号和第二回波信号之间的时间差成反比。
[0119]
第二组回波信号包括一对回波信号，因此计算第二导航频率可以参照上述实施例，在此不再赘述。
[0120]
步骤404，计算第一导航频率和第二导航频率之间的频率差值；对频率差值进行解缠绕处理，得到磁共振系统频率漂移量。
[0121]
本技术实施例中，在得到第一导航频率和第二导航频率之后，计算第一导航频率和第二导航频率之间的频率差值，抵消系统频率的手动调整量。但是，由于相位缠绕问题，频率差值不能直接反应磁共振系统频率漂移量，因此需要对频率差值进行解缠绕处理。具体如公式(9)：
[0122][0123]
其中，δf为磁共振系统频率漂移量，f
k
为第二导航频率，f
r
为第一导航频率。并且，设置频率阈值为n为正整数并取适当值，使得
[0124]
得到磁共振系统频率漂移量之后，可以采用以下任意一个步骤。
[0125]
步骤405，根据磁共振系统频率漂移量修正磁共振系统频率。
[0126]
本技术实施例中，在得到磁共振系统频率漂移量之后，可以根据磁共振系统频率漂移量修正磁共振系统频率。
[0127]
在实际应用中，如图6所示为采用现有技术方法获得的多期相的乳腺磁共振图像，乳腺6期相扫描，如果不采用上述步骤计算磁共振系统频率漂移量进行系统频率修正，图中第一排从左至右依次对应第一期相、第二期相和第三期相的乳腺磁共振图像，在第三期相的乳腺磁共振图像存在由于脂肪压制不均匀而产生的白色伪影区域(图中箭头所示)，第二排从左至右依次对应第四期相、第五期相和第六期相的乳腺磁共振图像，三个期相的存在由于脂肪压制不均匀而产生的明显白色伪影区域，即随着期相的增加，压脂效果逐渐变差。而如图7所示为采用如图2所述的方法对磁共振系统频率修正后获得的多期相的乳腺磁共振图像，同样是乳腺6期相扫描，采用上述步骤计算出的磁共振系统频率漂移量进行系统频率修正，不同期相的压脂效果基本一致。可见，采用本技术实施例的方式计算磁共振系统频率漂移量进行系统频率修正，可以提高系统频率的准确性，进而提高成像效果。
[0128]
步骤406，根据磁共振系统频率漂移量修正发射线圈的频率和/或接收线圈的频率。
[0129]
本技术实施例中，在得到磁共振系统频率漂移量之后，可以修正发射线圈的频率，或者修正接收线圈的频率；或者同时修正发射线圈的频率和接收线圈的频率，从而达到修
正系统频率的效果。
[0130]
上述磁共振系统频率漂移量的确定方法中，采用预设扫描序列进行第一扫描以获取第一组回波信号；在第一扫描的预设时长后，采用预设扫描序列进行第二扫描以获取第二组回波信号；分别根据第一组回波信号和第二组回波计算，得到第一导航频率和第二导航频率；计算第一导航频率和第二导航频率之间的频率差值；对频率差值进行解缠绕处理，得到磁共振系统频率漂移量；可以根据磁共振系统频率漂移量修正导航频率；也可以根据磁共振系统频率漂移量修正发射线圈的频率和/或接收线圈的频率。通过本技术实施例，根据第一导航频率和第二导航频率确定磁共振系统频率漂移量，并根据磁共振系统频率漂移量对系统频率进行修正，可以提高系统频率的准确性，进而提高成像效果。
[0131]
应该理解的是，虽然图2-图5的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，图2-图5中的至少一部分步骤可以包括多个步骤或者多个阶段，这些步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤中的步骤或者阶段的至少一部分轮流或者交替地执行。
[0132]
在一个实施例中，如图8所示，提供了一种磁共振系统频率漂移量的确定装置，包括：
[0133]
第一序列扫描模块501，用于采用预设扫描序列进行第一扫描以获取第一组回波信号；预设扫描序列为梯度双回波序列；
[0134]
第二序列扫描模块502，用于在第一扫描的预设时长后，采用预设扫描序列进行第二扫描以获取第二组回波信号；
[0135]
导航频率计算模块503，用于分别根据第一组回波信号和第二组回波计算，得到第一导航频率和第二导航频率；
[0136]
漂移量确定模块504，用于根据第一导航频率和第二导航频率确定磁共振系统频率漂移量。
[0137]
在其中一个实施例中，每组回波信号均包括第一回波信号和第二回波信号，上述导航频率计算模块503包括：
[0138]
频域转换子模块，用于针对每一组回波信号，分别对第一回波信号和第二回波信号进行频域转换，得到转换后的第三回波信号和第四回波信号；
[0139]
导航频率计算子模块，用于根据第三回波信号和第四回波信号计算导航频率。
[0140]
在其中一个实施例中，上述导航频率计算子模块，用于将第三回波信号与第四回波信号进行整合得到系统回波信号；系统回波信号采用振幅和相位角表示；对系统回波信号进行相位拟合，得到第一相位角；根据频率与相位角之间的对应关系和第一相位角计算导航频率。
[0141]
在其中一个实施例中，上述导航频率计算子模块，用于根据第三回波信号、第四回波信号和预先设置的相位角函数进行计算，得到第二相位角；根据频率与相位角之间的对应关系和第二相位角计算导航频率。
[0142]
在其中一个实施例中，上述频率与相位角之间的对应关系包括频率与相位角成正
比，与第一回波信号和第二回波信号之间的时间差成反比。
[0143]
在其中一个实施例中，上述第一组回波信号包括多对回波信号，每对回波信号包括两个回波信号，上述导航频率计算模块503，具体用于根据多对回波信号计算导航频率，得到多个候选导航频率；将多个候选导航频率的平均值确定为第一导航频率。
[0144]
在其中一个实施例中，上述漂移量确定模块504，具体用于计算第一导航频率和第二导航频率之间的频率差值；对频率差值进行解缠绕处理，得到磁共振系统频率漂移量。
[0145]
在其中一个实施例中，该装置还包括：
[0146]
系统频率修正模块，用于根据磁共振系统频率漂移量修正磁共振系统频率。
[0147]
在其中一个实施例中，该装置还包括：
[0148]
线圈频率修改模块，用于根据磁共振系统频率漂移量修正发射线圈的频率和/或接收线圈的频率。
[0149]
关于磁共振系统频率漂移量的确定装置的具体限定可以参见上文中对于磁共振系统频率漂移量的确定方法的限定，在此不再赘述。上述磁共振系统频率漂移量的确定装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。
[0150]
在一个实施例中，提供了一种计算机设备，该计算机设备可以是终端，其内部结构图可以如图9所示。该计算机设备包括通过系统总线连接的处理器、存储器、通信接口、显示屏和输入装置。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统和计算机程序。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的通信接口用于与外部的终端进行有线或无线方式的通信，无线方式可通过wifi、运营商网络、nfc(近场通信)或其他技术实现。该计算机程序被处理器执行时以实现一种磁共振系统频率漂移量的确定方法。该计算机设备的显示屏可以是液晶显示屏或者电子墨水显示屏，该计算机设备的输入装置可以是显示屏上覆盖的触摸层，也可以是计算机设备外壳上设置的按键、轨迹球或触控板，还可以是外接的键盘、触控板或鼠标等。
[0151]
本领域技术人员可以理解，图9中示出的结构，仅仅是与本技术方案相关的部分结构的框图，并不构成对本技术方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
[0152]
在一个实施例中，提供了一种计算机设备，包括存储器和处理器，存储器中存储有计算机程序，该处理器执行计算机程序时实现以下步骤：
[0153]
采用预设扫描序列进行第一扫描以获取第一组回波信号；预设扫描序列为梯度双回波序列；
[0154]
在第一扫描的预设时长后，采用预设扫描序列进行第二扫描以获取第二组回波信号；
[0155]
分别根据第一组回波信号和第二组回波计算，得到第一导航频率和第二导航频率；
[0156]
根据第一导航频率和第二导航频率确定磁共振系统频率漂移量。
[0157]
在一个实施例中，上述每组回波信号均包括第一回波信号和第二回波信号，该处
理器执行计算机程序时实现以下步骤：
[0158]
针对每一组回波信号，分别对第一回波信号和第二回波信号进行频域转换，得到转换后的第三回波信号和第四回波信号；
[0159]
根据第三回波信号和第四回波信号计算导航频率。
[0160]
在一个实施例中，该处理器执行计算机程序时实现以下步骤：
[0161]
将第三回波信号与第四回波信号进行整合得到系统回波信号；系统回波信号采用振幅和相位角表示；
[0162]
对系统回波信号进行相位拟合，得到第一相位角；
[0163]
根据频率与相位角之间的对应关系和第一相位角计算导航频率。
[0164]
在一个实施例中，该处理器执行计算机程序时实现以下步骤：
[0165]
根据第三回波信号、第四回波信号和预先设置的相位角函数进行计算，得到第二相位角；
[0166]
根据频率与相位角之间的对应关系和第二相位角计算导航频率。
[0167]
在一个实施例中，上述频率与相位角之间的对应关系包括频率与相位角成正比，与第一回波信号和第二回波信号之间的时间差成反比。
[0168]
在一个实施例中，第一组回波信号包括多对回波信号，每对回波信号包括两个回波信号，该处理器执行计算机程序时实现以下步骤：
[0169]
针对第一组回波信号，根据多对回波信号计算导航频率，得到多个候选导航频率；
[0170]
将多个候选导航频率的平均值确定为第一导航频率。
[0171]
在一个实施例中，该处理器执行计算机程序时实现以下步骤：
[0172]
计算第一导航频率和第二导航频率之间的频率差值；
[0173]
对频率差值进行解缠绕处理，得到磁共振系统频率漂移量。
[0174]
在一个实施例中，该处理器执行计算机程序时实现以下步骤：
[0175]
根据磁共振系统频率漂移量修正磁共振系统频率。
[0176]
在一个实施例中，该处理器执行计算机程序时实现以下步骤：
[0177]
根据磁共振系统频率漂移量修正发射线圈的频率和/或接收线圈的频率。
[0178]
在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：
[0179]
采用预设扫描序列进行第一扫描以获取第一组回波信号；预设扫描序列为梯度双回波序列；
[0180]
在第一扫描的预设时长后，采用预设扫描序列进行第二扫描以获取第二组回波信号；
[0181]
分别根据第一组回波信号和第二组回波计算，得到第一导航频率和第二导航频率；
[0182]
根据第一导航频率和第二导航频率确定磁共振系统频率漂移量。
[0183]
在一个实施例中，上述每组回波信号均包括第一回波信号和第二回波信号，计算机程序被处理器执行时实现以下步骤：
[0184]
针对每一组回波信号，分别对第一回波信号和第二回波信号进行频域转换，得到转换后的第三回波信号和第四回波信号；
[0185]
根据第三回波信号和第四回波信号计算导航频率。
[0186]
在一个实施例中，计算机程序被处理器执行时实现以下步骤：
[0187]
将第三回波信号与第四回波信号进行整合得到系统回波信号；系统回波信号采用振幅和相位角表示；
[0188]
对系统回波信号进行相位拟合，得到第一相位角；
[0189]
根据频率与相位角之间的对应关系和第一相位角计算导航频率。
[0190]
在一个实施例中，计算机程序被处理器执行时实现以下步骤：
[0191]
根据第三回波信号、第四回波信号和预先设置的相位角函数进行计算，得到第二相位角；
[0192]
根据频率与相位角之间的对应关系和第二相位角计算导航频率。
[0193]
在一个实施例中，上述频率与相位角之间的对应关系包括频率与相位角成正比，与第一回波信号和第二回波信号之间的时间差成反比。
[0194]
在一个实施例中，第一组回波信号包括多对回波信号，每对回波信号包括两个回波信号，计算机程序被处理器执行时实现以下步骤：
[0195]
针对第一组回波信号，根据多对回波信号计算导航频率，得到多个候选导航频率；
[0196]
将多个候选导航频率的平均值确定为第一导航频率。
[0197]
在一个实施例中，计算机程序被处理器执行时实现以下步骤：
[0198]
计算第一导航频率和第二导航频率之间的频率差值；
[0199]
对频率差值进行解缠绕处理，得到磁共振系统频率漂移量。
[0200]
在一个实施例中，计算机程序被处理器执行时实现以下步骤：
[0201]
根据磁共振系统频率漂移量修正磁共振系统频率。
[0202]
在一个实施例中，计算机程序被处理器执行时实现以下步骤：
[0203]
根据磁共振系统频率漂移量修正发射线圈的频率和/或接收线圈的频率。
[0204]
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本技术所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和易失性存储器中的至少一种。非易失性存储器可包括只读存储器(read-only memory，rom)、磁带、软盘、闪存或光存储器等。易失性存储器可包括随机存取存储器(random access memory，ram)或外部高速缓冲存储器。作为说明而非局限，ram可以是多种形式，比如静态随机存取存储器(static random access memory，sram)或动态随机存取存储器(dynamic random access memory，dram)等。
[0205]
以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
[0206]
以上所述实施例仅表达了本技术的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本技术构思的前提下，还可以做出若干变形和改进，这些都属于本技术的保护范围。因此，本技术专利的保护范围应以所附权利要求为准。