CN113625209B - 磁共振***频率漂移量的确定方法、装置和计算机设备 - Google Patents
磁共振***频率漂移量的确定方法、装置和计算机设备 Download PDFInfo
- Publication number
- CN113625209B CN113625209B CN202010387615.XA CN202010387615A CN113625209B CN 113625209 B CN113625209 B CN 113625209B CN 202010387615 A CN202010387615 A CN 202010387615A CN 113625209 B CN113625209 B CN 113625209B
- Authority
- CN
- China
- Prior art keywords
- frequency
- echo
- echo signals
- navigation
- echo signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000004590 computer program Methods 0.000 claims description 28
- 238000006243 chemical reaction Methods 0.000 claims description 24
- 238000004364 calculation method Methods 0.000 claims description 17
- 230000006870 function Effects 0.000 claims description 10
- 238000002592 echocardiography Methods 0.000 claims description 5
- 230000005284 excitation Effects 0.000 claims description 5
- 230000000694 effects Effects 0.000 abstract description 5
- 238000003384 imaging method Methods 0.000 abstract description 3
- 210000000481 breast Anatomy 0.000 description 12
- 238000010586 diagram Methods 0.000 description 7
- 238000004891 communication Methods 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000001629 suppression Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000000593 degrading effect Effects 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
- 238000004804 winding Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
- G01R33/38—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field
- G01R33/385—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field using gradient magnetic field coils
Landscapes
- Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本申请涉及一种磁共振***频率漂移量的确定方法、装置和计算机设备。所述方法包括：采用预设扫描序列进行第一扫描以获取第一组回波信号；所述预设扫描序列为梯度双回波序列；在所述第一扫描的预设时长后，采用所述预设扫描序列进行第二扫描以获取第二组回波信号；分别根据所述第一组回波信号和所述第二组回波计算，得到第一导航频率和第二导航频率；根据所述第一导航频率和所述第二导航频率确定磁共振***频率漂移量。采用本方法能够得到准确的磁共振***频率漂移量，从而提高***频率的准确性，提高成像效果。
Description
技术领域
本申请涉及磁共振技术领域，特别是涉及一种磁共振***频率漂移量的确定方法、装置和计算机设备。
背景技术
在磁共振成像过程中，***频率的准确性直接影响图像质量。在进行长时间的序列扫描时，磁共振***会因为电流不断切换引起的振动而发热，进而导致***频率发生漂移。
相关技术中，为了避免***频率漂移降低图像质量，通常会对***频率进行实时校准。但是，在扫描乳腺或者颈部的时候，常常需要用户手动对***频率进行调整。这时，如果仍然按照之前的校准方式进行校准，会将***频率的手动调整量当成频率漂移，从而导致***频率校准错误。
发明内容
基于此，有必要针对上述技术问题，提供一种能够避免将***频率的手动调整量当成频率漂移，从而避免***频率校准错误的磁共振***频率漂移量的确定方法、装置和计算机设备。
一种磁共振***频率漂移量的确定方法，该方法包括：
采用预设扫描序列进行第一扫描以获取第一组回波信号；预设扫描序列为梯度双回波序列；
在第一扫描的预设时长后，采用预设扫描序列进行第二扫描以获取第二组回波信号；
分别根据第一组回波信号和第二组回波计算，得到第一导航频率和第二导航频率；
根据第一导航频率和第二导航频率确定磁共振***频率漂移量。
在其中一个实施例中，上述每组回波信号均包括第一回波信号和第二回波信号，分别根据第一组回波信号和第二组回波计算，得到第一导航频率和第二导航频率，包括：
针对每一组回波信号，分别对第一回波信号和第二回波信号进行频域转换，得到转换后的第三回波信号和第四回波信号；
根据第三回波信号和第四回波信号计算导航频率。
在其中一个实施例中，上述根据第三回波信号和第四回波信号计算导航频率，包括：
将第三回波信号与第四回波信号进行整合得到***回波信号；***回波信号采用振幅和相位角表示；
对***回波信号进行相位拟合，得到第一相位角；
根据频率与相位角之间的对应关系和第一相位角计算导航频率。
在其中一个实施例中，上述根据第三回波信号和第四回波信号计算导航频率，包括：
根据第三回波信号、第四回波信号和预先设置的相位角函数进行计算，得到第二相位角；
根据频率与相位角之间的对应关系和第二相位角计算导航频率。
在其中一个实施例中，上述频率与相位角之间的对应关系包括频率与相位角成正比，与第一回波信号和第二回波信号之间的时间差成反比。
在其中一个实施例中，第一组回波信号包括多对回波信号，每对回波信号包括两个回波信号，上述分别根据第一组回波信号和第二组回波计算，得到第一导航频率和第二导航频率，包括：
针对第一组回波信号，根据多对回波信号计算导航频率，得到多个候选导航频率；
将多个候选导航频率的平均值确定为第一导航频率。
在其中一个实施例中，上述根据第一导航频率和第二导航频率确定磁共振***频率漂移量，包括：
计算第一导航频率和第二导航频率之间的频率差值；
对频率差值进行解缠绕处理，得到磁共振***频率漂移量。
在其中一个实施例中，在上述根据第一导航频率和第二导航频率确定磁共振***频率漂移量之后，该方法还包括：
根据磁共振***频率漂移量修正磁共振***频率。
在其中一个实施例中，在上述根据第一导航频率和第二导航频率确定磁共振***频率漂移量之后，该方法还包括：
根据磁共振***频率漂移量修正发射线圈的频率和/或接收线圈的频率。
一种磁共振***频率漂移量的确定装置，该装置包括：
第一序列扫描模块，用于采用预设扫描序列进行第一扫描以获取第一组回波信号；预设扫描序列为梯度双回波序列；
第二序列扫描模块，用于在第一扫描的预设时长后，采用预设扫描序列进行第二扫描以获取第二组回波信号；
导航频率计算模块，用于分别根据第一组回波信号和第二组回波计算，得到第一导航频率和第二导航频率；
漂移量确定模块，用于根据第一导航频率和第二导航频率确定磁共振***频率漂移量。
在其中一个实施例中，每组回波信号均包括第一回波信号和第二回波信号，上述导航频率计算模块包括：
频域转换子模块，用于针对每一组回波信号，分别对第一回波信号和第二回波信号进行频域转换，得到转换后的第三回波信号和第四回波信号；
导航频率计算子模块，用于根据第三回波信号和第四回波信号计算导航频率。
在其中一个实施例中，上述导航频率计算子模块，用于将第三回波信号与第四回波信号进行整合得到***回波信号；***回波信号采用振幅和相位角表示；对***回波信号进行相位拟合，得到第一相位角；根据频率与相位角之间的对应关系和第一相位角计算导航频率。
在其中一个实施例中，上述导航频率计算子模块，用于根据第三回波信号、第四回波信号和预先设置的相位角函数进行计算，得到第二相位角；根据频率与相位角之间的对应关系和第二相位角计算导航频率。
在其中一个实施例中，上述频率与相位角之间的对应关系包括频率与相位角成正比，与第一回波信号和第二回波信号之间的时间差成反比。
在其中一个实施例中，上述第一组回波信号包括多对回波信号，每对回波信号包括两个回波信号，上述第一导航频率计算模块，具体用于根据多对回波信号计算导航频率，得到多个候选导航频率；将多个候选导航频率的平均值确定为第一导航频率。
在其中一个实施例中，上述漂移量确定模块，具体用于计算第一导航频率和第二导航频率之间的频率差值；对频率差值进行解缠绕处理，得到磁共振***频率漂移量。
在其中一个实施例中，该装置还包括：
导航频率修正模块，用于根据磁共振***频率漂移量修正磁共振***频率。
在其中一个实施例中，该装置还包括：
线圈频率修改模块，用于根据磁共振***频率漂移量修正发射线圈的频率和/或接收线圈的频率。
一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现以下步骤：
采用预设扫描序列进行第一扫描以获取第一组回波信号；预设扫描序列为梯度双回波序列；
在第一扫描的预设时长后，采用预设扫描序列进行第二扫描以获取第二组回波信号；
分别根据第一组回波信号和第二组回波计算，得到第一导航频率和第二导航频率；
根据第一导航频率和第二导航频率确定磁共振***频率漂移量。
一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现以下步骤：
采用预设扫描序列进行第一扫描以获取第一组回波信号；预设扫描序列为梯度双回波序列；
在第一扫描的预设时长后，采用预设扫描序列进行第二扫描以获取第二组回波信号；
分别根据第一组回波信号和第二组回波计算，得到第一导航频率和第二导航频率；
根据第一导航频率和第二导航频率确定磁共振***频率漂移量。
上述磁共振***频率漂移量的确定方法、装置和计算机设备，采用预设扫描序列进行第一扫描以获取第一组回波信号；在第一扫描的预设时长后，采用预设扫描序列进行第二扫描以获取第二组回波信号；分别根据第一组回波信号和第二组回波计算，得到第一导航频率和第二导航频率；根据第一导航频率和第二导航频率确定磁共振***频率漂移量。通过本申请实施例，在用户手动调整***频率之后，采用预设扫描序列进行第一扫描并获取第一组回波信号，经过预设时长再采用预设扫描序列进行第二扫描并获取第二组回波信号，分别根据第一组回波信号和第二组回波信号计算第一导航频率和第二导航频率，由于第一导航频率和第二导航频率中均包含***频率的手动调整量，因此可以通过计算将***频率的手动调整量抵消掉，从而避免将***频率的手动调整量当成频率漂移而导致***频率校准错误的问题。
附图说明
图1为一个实施例中磁共振***频率漂移量的确定方法的应用环境图；
图2为一个实施例中磁共振***频率漂移量的确定方法的流程示意图；
图3为一个实施例中的波形图；
图4为一个实施例中根据回波信号计算导航频率步骤的流程示意图；
图5为另一个实施例中磁共振***频率漂移量的确定方法的流程示意图；
图6为采用现有技术方法获得的多期相的乳腺磁共振图像；
图7为一个实施例中对磁共振***频率修正后获得的多期相的乳腺磁共振图像；
图8为一个实施例中磁共振***频率漂移量的确定装置的结构框图；
图9为一个实施例中计算机设备的内部结构图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本申请进行进一步详细说明。应当理解，此处描述的具体实施例仅仅用以解释本申请，并不用于限定本申请。
本申请提供的磁共振***频率漂移量的确定方法，可以应用于如图1所示的应用环境中。该应用环境为磁共振***，该磁共振***至少包括梯度***101、射频***102和控制终端103。其中，梯度***101包括梯度线圈等设备，用于产生梯度脉冲；射频***102包括射频线圈等设备，射频线圈包括发射线圈和接收线圈，射频***102用于产生射频脉冲信号以及检测对象核自旋产生的磁共振信号；控制终端103与梯度***101和射频***102通过网络进行通信，用于控制梯度***101和射频***102。上述控制终端103可以但不限于是各种个人计算机、笔记本电脑和平板电脑。
在一个实施例中，如图2所示，提供了一种磁共振***频率漂移量的确定方法，以该方法应用于图1中的控制终端为例进行说明，包括以下步骤：
步骤201，采用预设扫描序列进行第一扫描以获取第一组回波信号；预设扫描序列为梯度双回波序列。
由于在扫描乳腺或者颈部的时候，常常需要用户手动对***频率进行调整，因此采用实时校准的方式，会将***频率的手动调整量当成频率漂移，从而导致***频率校准错误。为了避免这种情况，本申请实施例在检测到针对***频率的调整操作，即用户手动调整***频率之后，不再进行实时频率校准，而是采用动态频率校准。具体地，先采用预设扫描序列进行第一扫描，其中预设扫描序列可以是梯度双回波序列，见图3所示的波形图，RF为射频脉冲，GSS为选层方向梯度，GRO为频率编码方向梯度，ADC表示模数转换器接收的信号。图3所示的该梯度双回波序列只有频率编码梯度，没有相位编码梯度。并且，为了降低预设扫描序列对磁场的影响，在满足信噪比情况下，激发翻转角尽量小，各个逻辑轴的梯度净面积为零。
在采用预设扫描序列进行第一扫描的同时，采集第一组回波信号，见图3所示的波形图，沿GSS方向与射频脉冲相对应的施加选层梯度和该选层梯度后紧邻施加的散相梯度；沿GRO方向施加连续频率编码梯度，以避免梯度的频繁切换引起梯度涡流；ADC采样得到第一回波信号e0和第二回波信号e1，该两个回波所对应的GRO方向的梯度矩相同；回波时间TE(echo time)是射频脉冲中点到第一回波信号e0的时间间隔，ΔTE为第一回波信号e0与第二回波信号e1的回波时间差值。在此实施例中，预设扫描序列不包含相位编码方向的梯度，也即本申请实施例中采样得到第一回波信号e0和第二回波信号e1是未经过相位编码处理的。在一实施例中，TE的值满足保证检测对象的水信号和脂肪信号处于同相位。在其中一个实施例中，在用户手动调整***频率之前，先进行一次***频率的校准，以便在确定磁共振***频率漂移量之后，可以更加准确地修正***频率。其中，***频率为拉莫尔频率，即磁性原子核的进动频率或主磁场频率。
步骤202，在第一扫描的预设时长后，采用预设扫描序列进行第二扫描以获取第二组回波信号。
本申请实施例中，在获取第一组回波信号之后，***频率经过一段时间可能发生了漂移，因此采用预设扫描序列进行第二扫描，同时获取第二组回波信号。
例如，采用预设扫描序列进行一次扫描，5分钟后采用预设扫描序列再进行一次扫描，并获取第二组回波信号。本申请实施例对预设时长不作详细限定，可以根据实际情况进行设置。
步骤203，分别根据第一组回波信号和第二组回波计算，得到第一导航频率和第二导航频率。
本申请实施例中，根据第一组回波信号计算得到第一导航频率，根据第二组回波信号计算得到第二导航频率。具体地，对每组回波信号进行频域转换，然后根据转换后的回波信号计算相位角，再根据相位角和频率之间的对应关系计算导航频率。本申请实施例对具体计算方式不作详细限定，可以根据实际情况进行设置。
步骤204，根据第一导航频率和第二导航频率确定磁共振***频率漂移量。
本申请实施例中，由于第一导航频率中包含***频率的手动调整量，第二导航频率中也包含***频率的手动调整量，因此在得到第一导航频率和第二导航频率之后，可以计算第一导航频率和第二导航频率之间的频率差值，从而将***频率的手动调整量抵消掉，进而得到磁共振***频率漂移量。
上述磁共振***频率漂移量的确定方法中，采用预设扫描序列进行第一扫描以获取第一组回波信号；在第一扫描的预设时长后，采用预设扫描序列进行第二扫描以获取第二组回波信号；分别根据第一组回波信号和第二组回波计算，得到第一导航频率和第二导航频率；根据第一导航频率和第二导航频率确定磁共振***频率漂移量。本申请实施例中，第一导航频率和第二导航频率中均包含***频率的手动调整量，因此可以通过计算将***频率的手动调整量抵消掉，从而避免将***频率的手动调整量当成频率漂移而导致***频率校准错误的问题。
在一个实施例中，如图4所示，涉及根据回波信号计算导航频率的一种可选过程。在上述实施例的基础上，每组回波信号均包括第一回波信号和第二回波信号，因此计算第一导航频率和第二导航频率，均可以采用如下步骤：
步骤301，针对每一组回波信号，分别对第一回波信号和第二回波信号进行频域转换，得到转换后的第三回波信号和第四回波信号。
本申请实施例中，采用梯度双回波序列进行扫描，可以采集到两个回波信号，即第一回波信号和第二回波信号。之后，对第一回波信号进行频域转换得到第三回波信号，对第二回波信号进行频域转换得到第四回波信号。具体地，频域转换可以采用傅里叶变换，如公式(1)(2)：
E0＝FFT(e0)-------------------------------------------(1)
E1＝FFT(e1)--------------------------------------------(2)
其中，e0为第一回波信号，e1为第二回波信号，E0为第三回波信号，E1为第四回波信号。
步骤302，根据第三回波信号和第四回波信号计算导航频率。
本申请实施例中，得到频域下的第三回波信号和第四回波信号之后，可以采用以下两种方式计算导航频率。
方式一：将第三回波信号与第四回波信号进行整合得到***回波信号；***回波信号采用振幅和相位角表示；对***回波信号进行相位拟合，得到第一相位角；根据频率与相位角之间的对应关系和第一相位角计算导航频率。
具体地，先将两个回波信号整合为一个***回波信号，其中***回波信号采用振幅和相位角表示。整合过程可以如公式(3)：
然后将E(i)记作公式(4)：
其中，E0(i)*E1(i)为第三回波信号E0和第四回波信号E1共轭相乘；Ai为***回波信号E(i)的幅值，为***回波信号E(i)的相位角，n表示所采集得到的梯度双回波的数量，n≥1。
在得到***回波信号之后，可以采用最小二乘法进行相位拟合，如公式(5)：
其中，θ1为拟合得到的第一相位角。
之后，根据频率与相位角之间的对应关系，以及拟合得到的第一相位角计算导航频率，具体如公式(6)：
上述频率与相位角之间的对应关系包括：频率与相位角成正比，与第一回波信号和第二回波信号之间的时间差成反比。其中，f为计算出的导航频率，ΔTE为第一回波信号和第二回波信号之间的时间差。
方式二：根据第三回波信号、第四回波信号和预先设置的相位角函数进行计算，得到第二相位角；根据频率与相位角之间的对应关系和第二相位角计算导航频率。
具体地，可以直接根据相位角函数进行计算，得到第二相位角，如公式(7)：
θ2＝angle(∑i(E0(i)*E1(i)))-----------------------(7)
其中，θ1为第二相位角。
之后，根据频率与相位角之间的对应关系，以及计算得到的第二相位角计算导航频率，具体如公式(8)：
上述频率与相位角之间的对应关系包括：频率与相位角成正比，与第一回波信号和第二回波信号之间的时间差成反比。其中，f为计算出的导航频率，ΔTE为第一回波信号和第二回波信号之间的时间差。
上述根据回波信号计算导航频率的过程中，针对每组回波信号，分别对第一回波信号和第二回波信号进行频域转换，得到转换后的第三回波信号和第四回波信号；根据第三回波信号和第四回波信号计算导航频率。通过本申请实施例，先将时域表示的回波信号转换为频域表示的回波信号，由频域表示的回波信号可以计算出相位角，再根据频率与相位角之间的对应关系可以计算出导航频率，从而确定第一导航频率和第二导航频率，进而得到磁共振***频率漂移量。
在一个实施例中，如图5所示，提供了一种磁共振***频率漂移量的确定方法，以该方法应用于图1中的控制终端为例进行说明，包括以下步骤：
步骤401，采用预设扫描序列进行第一扫描以获取第一组回波信号；预设扫描序列为梯度双回波序列。
步骤402，在第一扫描的预设时长后，采用预设扫描序列进行第二扫描以获取第二组回波信号。
步骤403，分别根据第一组回波信号和第二组回波计算，得到第一导航频率和第二导航频率。
本申请实施例中，为了使第一导航频率更加准确，针对第一组回波信号，可以根据多对回波信号计算导航频率，得到多个候选导航频率；将多个候选导航频率的平均值确定为第一导航频率。具体地，采集多对回波信号，每对回波信号均包括两个回波信号；针对每对回波信号进行导航频率计算，得到对应的候选导航频率，这样就可以得到多个候选导航频率；最后，计算多个候选导航频率的平均值，将平均值作为第一导航频率。
在其中一个实施例中，针对每对回波信号计算导航频率，得到对应的候选导航频率可以包括：分别对第一回波信号和第二回波信号进行频域转换，得到转换后的第三回波信号和第四回波信号；根据第三回波信号和第四回波信号计算导航频率。
在其中一个实施例中，根据第三回波信号和第四回波信号计算导航频率，包括：将第三回波信号与第四回波信号进行整合得到***回波信号；***回波信号采用振幅和相位角表示；对***回波信号进行相位拟合，得到第一相位角；根据频率与相位角之间的对应关系和第一相位角计算导航频率。
在其中一个实施例中，根据第三回波信号和第四回波信号计算导航频率，包括：根据第三回波信号、第四回波信号和预先设置的相位角函数进行计算，得到第二相位角；根据频率与相位角之间的对应关系和第二相位角计算导航频率。
在其中一个实施例中，频率与相位角之间的对应关系包括频率与相位角成正比，与第一回波信号和第二回波信号之间的时间差成反比。
第二组回波信号包括一对回波信号，因此计算第二导航频率可以参照上述实施例，在此不再赘述。
步骤404，计算第一导航频率和第二导航频率之间的频率差值；对频率差值进行解缠绕处理，得到磁共振***频率漂移量。
本申请实施例中，在得到第一导航频率和第二导航频率之后，计算第一导航频率和第二导航频率之间的频率差值，抵消***频率的手动调整量。但是，由于相位缠绕问题，频率差值不能直接反应磁共振***频率漂移量，因此需要对频率差值进行解缠绕处理。具体如公式(9)：
其中，Δf为磁共振***频率漂移量，fk为第二导航频率，fr为第一导航频率。并且，设置频率阈值为n为正整数并取适当值，使得/>
得到磁共振***频率漂移量之后，可以采用以下任意一个步骤。
步骤405，根据磁共振***频率漂移量修正磁共振***频率。
本申请实施例中，在得到磁共振***频率漂移量之后，可以根据磁共振***频率漂移量修正磁共振***频率。
在实际应用中，如图6所示为采用现有技术方法获得的多期相的乳腺磁共振图像，乳腺6期相扫描，如果不采用上述步骤计算磁共振***频率漂移量进行***频率修正，图中第一排从左至右依次对应第一期相、第二期相和第三期相的乳腺磁共振图像，在第三期相的乳腺磁共振图像存在由于脂肪压制不均匀而产生的白色伪影区域(图中箭头所示)，第二排从左至右依次对应第四期相、第五期相和第六期相的乳腺磁共振图像，三个期相的存在由于脂肪压制不均匀而产生的明显白色伪影区域，即随着期相的增加，压脂效果逐渐变差。而如图7所示为采用如图2所述的方法对磁共振***频率修正后获得的多期相的乳腺磁共振图像，同样是乳腺6期相扫描，采用上述步骤计算出的磁共振***频率漂移量进行***频率修正，不同期相的压脂效果基本一致。可见，采用本申请实施例的方式计算磁共振***频率漂移量进行***频率修正，可以提高***频率的准确性，进而提高成像效果。
步骤406，根据磁共振***频率漂移量修正发射线圈的频率和/或接收线圈的频率。
本申请实施例中，在得到磁共振***频率漂移量之后，可以修正发射线圈的频率，或者修正接收线圈的频率；或者同时修正发射线圈的频率和接收线圈的频率，从而达到修正***频率的效果。
上述磁共振***频率漂移量的确定方法中，采用预设扫描序列进行第一扫描以获取第一组回波信号；在第一扫描的预设时长后，采用预设扫描序列进行第二扫描以获取第二组回波信号；分别根据第一组回波信号和第二组回波计算，得到第一导航频率和第二导航频率；计算第一导航频率和第二导航频率之间的频率差值；对频率差值进行解缠绕处理，得到磁共振***频率漂移量；可以根据磁共振***频率漂移量修正导航频率；也可以根据磁共振***频率漂移量修正发射线圈的频率和/或接收线圈的频率。通过本申请实施例，根据第一导航频率和第二导航频率确定磁共振***频率漂移量，并根据磁共振***频率漂移量对***频率进行修正，可以提高***频率的准确性，进而提高成像效果。
应该理解的是，虽然图2-图5的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，图2-图5中的至少一部分步骤可以包括多个步骤或者多个阶段，这些步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤中的步骤或者阶段的至少一部分轮流或者交替地执行。
在一个实施例中，如图8所示，提供了一种磁共振***频率漂移量的确定装置，包括：
第一序列扫描模块501，用于采用预设扫描序列进行第一扫描以获取第一组回波信号；预设扫描序列为梯度双回波序列；
第二序列扫描模块502，用于在第一扫描的预设时长后，采用预设扫描序列进行第二扫描以获取第二组回波信号；
导航频率计算模块503，用于分别根据第一组回波信号和第二组回波计算，得到第一导航频率和第二导航频率；
漂移量确定模块504，用于根据第一导航频率和第二导航频率确定磁共振***频率漂移量。
在其中一个实施例中，每组回波信号均包括第一回波信号和第二回波信号，上述导航频率计算模块503包括：
频域转换子模块，用于针对每一组回波信号，分别对第一回波信号和第二回波信号进行频域转换，得到转换后的第三回波信号和第四回波信号；
导航频率计算子模块，用于根据第三回波信号和第四回波信号计算导航频率。
在其中一个实施例中，上述导航频率计算子模块，用于将第三回波信号与第四回波信号进行整合得到***回波信号；***回波信号采用振幅和相位角表示；对***回波信号进行相位拟合，得到第一相位角；根据频率与相位角之间的对应关系和第一相位角计算导航频率。
在其中一个实施例中，上述导航频率计算子模块，用于根据第三回波信号、第四回波信号和预先设置的相位角函数进行计算，得到第二相位角；根据频率与相位角之间的对应关系和第二相位角计算导航频率。
在其中一个实施例中，上述频率与相位角之间的对应关系包括频率与相位角成正比，与第一回波信号和第二回波信号之间的时间差成反比。
在其中一个实施例中，上述第一组回波信号包括多对回波信号，每对回波信号包括两个回波信号，上述导航频率计算模块503，具体用于根据多对回波信号计算导航频率，得到多个候选导航频率；将多个候选导航频率的平均值确定为第一导航频率。
在其中一个实施例中，上述漂移量确定模块504，具体用于计算第一导航频率和第二导航频率之间的频率差值；对频率差值进行解缠绕处理，得到磁共振***频率漂移量。
在其中一个实施例中，该装置还包括：
***频率修正模块，用于根据磁共振***频率漂移量修正磁共振***频率。
在其中一个实施例中，该装置还包括：
线圈频率修改模块，用于根据磁共振***频率漂移量修正发射线圈的频率和/或接收线圈的频率。
关于磁共振***频率漂移量的确定装置的具体限定可以参见上文中对于磁共振***频率漂移量的确定方法的限定，在此不再赘述。上述磁共振***频率漂移量的确定装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。
在一个实施例中，提供了一种计算机设备，该计算机设备可以是终端，其内部结构图可以如图9所示。该计算机设备包括通过***总线连接的处理器、存储器、通信接口、显示屏和输入装置。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作***和计算机程序。该内存储器为非易失性存储介质中的操作***和计算机程序的运行提供环境。该计算机设备的通信接口用于与外部的终端进行有线或无线方式的通信，无线方式可通过WIFI、运营商网络、NFC(近场通信)或其他技术实现。该计算机程序被处理器执行时以实现一种磁共振***频率漂移量的确定方法。该计算机设备的显示屏可以是液晶显示屏或者电子墨水显示屏，该计算机设备的输入装置可以是显示屏上覆盖的触摸层，也可以是计算机设备外壳上设置的按键、轨迹球或触控板，还可以是外接的键盘、触控板或鼠标等。
本领域技术人员可以理解，图9中示出的结构，仅仅是与本申请方案相关的部分结构的框图，并不构成对本申请方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
在一个实施例中，提供了一种计算机设备，包括存储器和处理器，存储器中存储有计算机程序，该处理器执行计算机程序时实现以下步骤：
采用预设扫描序列进行第一扫描以获取第一组回波信号；预设扫描序列为梯度双回波序列；
在第一扫描的预设时长后，采用预设扫描序列进行第二扫描以获取第二组回波信号；
分别根据第一组回波信号和第二组回波计算，得到第一导航频率和第二导航频率；
根据第一导航频率和第二导航频率确定磁共振***频率漂移量。
在一个实施例中，上述每组回波信号均包括第一回波信号和第二回波信号，该处理器执行计算机程序时实现以下步骤：
针对每一组回波信号，分别对第一回波信号和第二回波信号进行频域转换，得到转换后的第三回波信号和第四回波信号；
根据第三回波信号和第四回波信号计算导航频率。
在一个实施例中，该处理器执行计算机程序时实现以下步骤：
将第三回波信号与第四回波信号进行整合得到***回波信号；***回波信号采用振幅和相位角表示；
对***回波信号进行相位拟合，得到第一相位角；
根据频率与相位角之间的对应关系和第一相位角计算导航频率。
在一个实施例中，该处理器执行计算机程序时实现以下步骤：
根据第三回波信号、第四回波信号和预先设置的相位角函数进行计算，得到第二相位角；
根据频率与相位角之间的对应关系和第二相位角计算导航频率。
在一个实施例中，上述频率与相位角之间的对应关系包括频率与相位角成正比，与第一回波信号和第二回波信号之间的时间差成反比。
在一个实施例中，第一组回波信号包括多对回波信号，每对回波信号包括两个回波信号，该处理器执行计算机程序时实现以下步骤：
针对第一组回波信号，根据多对回波信号计算导航频率，得到多个候选导航频率；
将多个候选导航频率的平均值确定为第一导航频率。
在一个实施例中，该处理器执行计算机程序时实现以下步骤：
计算第一导航频率和第二导航频率之间的频率差值；
对频率差值进行解缠绕处理，得到磁共振***频率漂移量。
在一个实施例中，该处理器执行计算机程序时实现以下步骤：
根据磁共振***频率漂移量修正磁共振***频率。
在一个实施例中，该处理器执行计算机程序时实现以下步骤：
根据磁共振***频率漂移量修正发射线圈的频率和/或接收线圈的频率。
在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：
采用预设扫描序列进行第一扫描以获取第一组回波信号；预设扫描序列为梯度双回波序列；
在第一扫描的预设时长后，采用预设扫描序列进行第二扫描以获取第二组回波信号；
分别根据第一组回波信号和第二组回波计算，得到第一导航频率和第二导航频率；
根据第一导航频率和第二导航频率确定磁共振***频率漂移量。
在一个实施例中，上述每组回波信号均包括第一回波信号和第二回波信号，计算机程序被处理器执行时实现以下步骤：
针对每一组回波信号，分别对第一回波信号和第二回波信号进行频域转换，得到转换后的第三回波信号和第四回波信号；
根据第三回波信号和第四回波信号计算导航频率。
在一个实施例中，计算机程序被处理器执行时实现以下步骤：
将第三回波信号与第四回波信号进行整合得到***回波信号；***回波信号采用振幅和相位角表示；
对***回波信号进行相位拟合，得到第一相位角；
根据频率与相位角之间的对应关系和第一相位角计算导航频率。
在一个实施例中，计算机程序被处理器执行时实现以下步骤：
根据第三回波信号、第四回波信号和预先设置的相位角函数进行计算，得到第二相位角；
根据频率与相位角之间的对应关系和第二相位角计算导航频率。
在一个实施例中，上述频率与相位角之间的对应关系包括频率与相位角成正比，与第一回波信号和第二回波信号之间的时间差成反比。
在一个实施例中，第一组回波信号包括多对回波信号，每对回波信号包括两个回波信号，计算机程序被处理器执行时实现以下步骤：
针对第一组回波信号，根据多对回波信号计算导航频率，得到多个候选导航频率；
将多个候选导航频率的平均值确定为第一导航频率。
在一个实施例中，计算机程序被处理器执行时实现以下步骤：
计算第一导航频率和第二导航频率之间的频率差值；
对频率差值进行解缠绕处理，得到磁共振***频率漂移量。
在一个实施例中，计算机程序被处理器执行时实现以下步骤：
根据磁共振***频率漂移量修正磁共振***频率。
在一个实施例中，计算机程序被处理器执行时实现以下步骤：
根据磁共振***频率漂移量修正发射线圈的频率和/或接收线圈的频率。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和易失性存储器中的至少一种。非易失性存储器可包括只读存储器(Read-Only Memory，ROM)、磁带、软盘、闪存或光存储器等。易失性存储器可包括随机存取存储器(Random Access Memory，RAM)或外部高速缓冲存储器。作为说明而非局限，RAM可以是多种形式，比如静态随机存取存储器(Static Random Access Memory，SRAM)或动态随机存取存储器(Dynamic Random Access Memory，DRAM)等。
以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本申请构思的前提下，还可以做出若干变形和改进，这些都属于本申请的保护范围。因此，本申请专利的保护范围应以所附权利要求为准。
Claims (7)
1.一种磁共振***频率漂移量的确定方法，其特征在于，所述方法包括：
采用预设扫描序列进行第一扫描以获取第一组回波信号；所述预设扫描序列为梯度双回波序列，所述第一组回波信号包括多对回波信号，每对回波信号包括同次射频脉冲激发后先后采集的两个回波信号；
在所述第一扫描的预设时长后，采用所述预设扫描序列进行第二扫描以获取第二组回波信号，所述第二组回波信号包括多对回波信号，每对回波信号包括同次射频脉冲激发后先后采集的两个回波信号；
分别根据所述第一组回波信号和所述第二组回波信号计算，得到第一导航频率和第二导航频率；
计算所述第一导航频率和所述第二导航频率之间的频率差值；
对所述频率差值进行解缠绕处理，得到所述磁共振***频率漂移量；
其中，每对回波信号包括第一回波信号和第二回波信号，所述分别根据所述第一组回波信号和所述第二组回波计算，得到第一导航频率和第二导航频率，包括：
针对每一组回波信号，分别对所述第一回波信号和所述第二回波信号进行频域转换，得到转换后的第三回波信号和第四回波信号；
根据所述第三回波信号、所述第四回波信号和预先设置的相位角函数确定相位角，并根据所述相位角、所述第一回波信号和所述第二回波信号的时间差计算导航频率。
2.根据权利要求1所述的方法，其特征在于，频率与相位角成正比，与所述第一回波信号和所述第二回波信号之间的时间差成反比。
3.根据权利要求1所述的方法，其特征在于，所述分别根据所述第一组回波信号和所述第二组回波计算，得到第一导航频率和第二导航频率，包括：
针对所述第一组回波信号，根据所述多对回波信号计算导航频率，得到多个候选导航频率；
将多个所述候选导航频率的平均值确定为所述第一导航频率。
4.根据权利要求1所述的方法，其特征在于，在所述根据所述第一导航频率和所述第二导航频率确定磁共振***频率漂移量之后，所述方法还包括：
根据所述磁共振***频率漂移量修正磁共振***频率；或者，
根据所述磁共振***频率漂移量修正发射线圈的频率和/或接收线圈的频率。
5.一种磁共振***频率漂移量的确定装置，其特征在于，所述装置包括：
第一序列扫描模块，用于采用预设扫描序列进行第一扫描以获取第一组回波信号；所述预设扫描序列为梯度双回波序列，所述第一组回波信号包括多对回波信号，每对回波信号包括同次射频脉冲激发后先后采集的两个回波信号；
第二序列扫描模块，用于在所述第一扫描的预设时长后，采用所述预设扫描序列进行第二扫描以获取第二组回波信号，所述第二组回波信号包括多对回波信号，每对回波信号包括同次射频脉冲激发后先后采集的两个回波信号；
导航频率计算模块，用于分别根据所述第一组回波信号和所述第二组回波计算，得到第一导航频率和第二导航频率；
漂移量确定模块，用于计算所述第一导航频率和所述第二导航频率之间的频率差值；对所述频率差值进行解缠绕处理，得到所述磁共振***频率漂移量；
其中，每对回波信号包括第一回波信号和第二回波信号，所述导航频率计算模块，具体用于针对每一组回波信号，分别对所述第一回波信号和所述第二回波信号进行频域转换，得到转换后的第三回波信号和第四回波信号；
根据所述第三回波信号、所述第四回波信号和预先设置的相位角函数确定相位角，并根据所述相位角、所述第一回波信号和所述第二回波信号的时间差计算导航频率。
6.一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至4中任一项所述的方法的步骤。
7.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1至4中任一项所述的方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010387615.XA CN113625209B (zh) | 2020-05-09 | 2020-05-09 | 磁共振***频率漂移量的确定方法、装置和计算机设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010387615.XA CN113625209B (zh) | 2020-05-09 | 2020-05-09 | 磁共振***频率漂移量的确定方法、装置和计算机设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113625209A CN113625209A (zh) | 2021-11-09 |
| CN113625209B true CN113625209B (zh) | 2024-02-27 |
Family
ID=78377490
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010387615.XA Active CN113625209B (zh) | 2020-05-09 | 2020-05-09 | 磁共振***频率漂移量的确定方法、装置和计算机设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113625209B (zh) |
Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4728906A (en) * | 1985-10-21 | 1988-03-01 | Wiltron Measurements Limited | Tuning and calibration circuits for CW and sweep frequency signal generators |
| CN1216242A (zh) * | 1997-10-17 | 1999-05-12 | 日立医药株式会社 | 磁共振成像方法及实施该方法的装置 |
| CN1693885A (zh) * | 2004-04-30 | 2005-11-09 | 西门子公司 | 在磁共振光谱学试验中动态检测磁共振频率的方法 |
| CN101034149A (zh) * | 2007-04-19 | 2007-09-12 | 新奥博为技术有限公司 | 一种基于预扫描测量校正磁共振信号相位偏移的方法 |
| CN102078196A (zh) * | 2009-11-27 | 2011-06-01 | 西门子公司 | 借助磁共振采集相应于对象的呼吸运动的信号数据的方法 |
| CN102232849A (zh) * | 2010-04-30 | 2011-11-09 | 西门子(深圳)磁共振有限公司 | 一种动态频率漂移校正方法 |
| CN102654574A (zh) * | 2012-05-21 | 2012-09-05 | 西安电子科技大学 | 基于fpga的多普勒中心估计方法 |
| CN102866372A (zh) * | 2011-07-08 | 2013-01-09 | 西门子公司 | 对磁共振设备进行频率校准的方法及相应的磁共振设备 |
| CN102944875A (zh) * | 2012-11-21 | 2013-02-27 | 西安电子科技大学 | Isar图像距离单元选择横向定标方法 |
| CN105247382A (zh) * | 2013-03-29 | 2016-01-13 | 皇家飞利浦有限公司 | 具有重影伪影降低的磁共振成像***和其操作的方法 |
| CN105997074A (zh) * | 2016-04-28 | 2016-10-12 | 华东师范大学 | 一种磁共振定量磁化率成像多回波相位拟合方法 |
| CN108375746A (zh) * | 2017-01-18 | 2018-08-07 | 上海联影医疗科技有限公司 | 一种相位反卷折方法与设备 |
| CN108627834A (zh) * | 2018-06-07 | 2018-10-09 | 北京城建勘测设计研究院有限责任公司 | 一种基于地基InSAR的地铁路基结构监测方法及装置 |
| CN108742623A (zh) * | 2018-02-05 | 2018-11-06 | 苏州朗润医疗***有限公司 | 基于三维梯度双回波的动态匀场方法 |
| CN110068780A (zh) * | 2019-05-07 | 2019-07-30 | 上海东软医疗科技有限公司 | 一种磁共振成像的自适应方法、装置和磁共振成像*** |
| CN110907872A (zh) * | 2018-09-18 | 2020-03-24 | 西门子医疗有限公司 | 产生磁共振图像的方法、程序、数据载体和磁共振设备 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| IL196487A (en) * | 2009-01-13 | 2016-03-31 | Aspect Imaging Ltd | Means for buying sharp resolution mri |
| DE102011006230B4 (de) * | 2011-03-28 | 2013-01-24 | Siemens Aktiengesellschaft | Pixelweise Korrektur einer Phaseninformation in MR-Bildern unter Verwendung eines Navigatorsignals |
| JP5925529B2 (ja) * | 2011-03-31 | 2016-05-25 | 株式会社東芝 | 磁気共鳴イメージング装置 |
| WO2015017764A1 (en) * | 2013-08-02 | 2015-02-05 | The General Hosptial Corporation | System and method for real-time frequency correction for magnetic resonance imaging |
-
2020
- 2020-05-09 CN CN202010387615.XA patent/CN113625209B/zh active Active
Patent Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4728906A (en) * | 1985-10-21 | 1988-03-01 | Wiltron Measurements Limited | Tuning and calibration circuits for CW and sweep frequency signal generators |
| CN1216242A (zh) * | 1997-10-17 | 1999-05-12 | 日立医药株式会社 | 磁共振成像方法及实施该方法的装置 |
| CN1693885A (zh) * | 2004-04-30 | 2005-11-09 | 西门子公司 | 在磁共振光谱学试验中动态检测磁共振频率的方法 |
| CN101034149A (zh) * | 2007-04-19 | 2007-09-12 | 新奥博为技术有限公司 | 一种基于预扫描测量校正磁共振信号相位偏移的方法 |
| CN102078196A (zh) * | 2009-11-27 | 2011-06-01 | 西门子公司 | 借助磁共振采集相应于对象的呼吸运动的信号数据的方法 |
| CN102232849A (zh) * | 2010-04-30 | 2011-11-09 | 西门子(深圳)磁共振有限公司 | 一种动态频率漂移校正方法 |
| CN102866372A (zh) * | 2011-07-08 | 2013-01-09 | 西门子公司 | 对磁共振设备进行频率校准的方法及相应的磁共振设备 |
| CN102654574A (zh) * | 2012-05-21 | 2012-09-05 | 西安电子科技大学 | 基于fpga的多普勒中心估计方法 |
| CN102944875A (zh) * | 2012-11-21 | 2013-02-27 | 西安电子科技大学 | Isar图像距离单元选择横向定标方法 |
| CN105247382A (zh) * | 2013-03-29 | 2016-01-13 | 皇家飞利浦有限公司 | 具有重影伪影降低的磁共振成像***和其操作的方法 |
| CN105997074A (zh) * | 2016-04-28 | 2016-10-12 | 华东师范大学 | 一种磁共振定量磁化率成像多回波相位拟合方法 |
| CN108375746A (zh) * | 2017-01-18 | 2018-08-07 | 上海联影医疗科技有限公司 | 一种相位反卷折方法与设备 |
| CN108742623A (zh) * | 2018-02-05 | 2018-11-06 | 苏州朗润医疗***有限公司 | 基于三维梯度双回波的动态匀场方法 |
| CN108627834A (zh) * | 2018-06-07 | 2018-10-09 | 北京城建勘测设计研究院有限责任公司 | 一种基于地基InSAR的地铁路基结构监测方法及装置 |
| CN110907872A (zh) * | 2018-09-18 | 2020-03-24 | 西门子医疗有限公司 | 产生磁共振图像的方法、程序、数据载体和磁共振设备 |
| CN110068780A (zh) * | 2019-05-07 | 2019-07-30 | 上海东软医疗科技有限公司 | 一种磁共振成像的自适应方法、装置和磁共振成像*** |
Non-Patent Citations (2)
| Title |
|---|
| 基于动态频率漂移校正的二项式水激发磁共振成像;修峰;《中国优秀硕士学位论文全文数据库 (医药卫生科技辑)》(第5期);第4章 * |
| 定量磁化率成像多回波相位拟合算法研究;赵欣欣;薄斌仕;刘田;王乙;李建奇;;波谱学杂志(第04期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113625209A (zh) | 2021-11-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9835700B2 (en) | Magnetic resonance imaging apparatus | |
| US9046590B2 (en) | Magnetic resonance imaging apparatus phase correction using one or more prescans with identical readout and slice gradients | |
| US7323873B2 (en) | Phase error measuring method, MR imaging method, and MRI system | |
| US8854037B2 (en) | Method for compensating for eddy current fields in magnetic resonance images | |
| JP2012205897A (ja) | 磁気共鳴画像内の位相情報の補正方法および磁気共鳴設備 | |
| JPWO2008152937A1 (ja) | 磁気共鳴イメージング装置及び傾斜磁場に起因する誤差補正方法 | |
| US7446527B2 (en) | Methods and apparatus for reducing artifacts in MRI images | |
| JP2012005823A (ja) | 磁気共鳴イメージング装置及び方法 | |
| US9971007B2 (en) | Method and apparatus for accelerated magnetic resonance imaging | |
| US9846215B2 (en) | MRI embodiments for controlling an arrangement order of multiple echoes in a k-space | |
| US10705172B2 (en) | Magnetic resonance apparatus and method for dynamic adjustment thereof with multiple adjustment parameters | |
| CN111190131B (zh) | 用于磁共振图像重建的方法和磁共振*** | |
| US11624793B2 (en) | Method for acquiring measurement data using a magnetic resonance system with correction of k-space trajectories | |
| US10473743B2 (en) | Method and magnetic resonance apparatus for determining a scan sequence based on a representation of a pulse response in k-space of the gradient system | |
| US9417306B2 (en) | Magnetic resonance trajectory correcting with GRAPPA operator gridding | |
| CN113625209B (zh) | 磁共振***频率漂移量的确定方法、装置和计算机设备 | |
| US11402454B2 (en) | Correction of distorted diffusion-weighted magnetic resonance image data | |
| KR101778139B1 (ko) | 자기 공명 영상 장치 및 자기 공명 영상 장치의 자장 보정 방법 | |
| CN115586473B (zh) | 磁场匀场方法、装置、磁共振成像设备及可读存储介质 | |
| JP3938717B2 (ja) | 磁気共鳴撮像装置 | |
| CN116990736B (zh) | 一种核磁共振成像方法、装置、设备及存储介质 | |
| JP2002306454A (ja) | 磁気共鳴映像装置 | |
| US20220206097A1 (en) | Magnetic resonance image processing method and device, and computer readable storage medium | |
| US10823796B2 (en) | Magnetic resonance imaging apparatus and pulse setting method | |
| JP4901627B2 (ja) | 磁気共鳴撮影装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |