CN113219391B - 一种磁共振成像***中加速涡流矫正的方法和设备 - Google Patents
一种磁共振成像***中加速涡流矫正的方法和设备 Download PDFInfo
- Publication number
- CN113219391B CN113219391B CN202110587406.4A CN202110587406A CN113219391B CN 113219391 B CN113219391 B CN 113219391B CN 202110587406 A CN202110587406 A CN 202110587406A CN 113219391 B CN113219391 B CN 113219391B
- Authority
- CN
- China
- Prior art keywords
- eddy current
- time
- gradient
- axis
- delay time
- 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
- 238000012937 correction Methods 0.000 title claims abstract description 78
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000002595 magnetic resonance imaging Methods 0.000 title claims abstract description 22
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 20
- 238000005457 optimization Methods 0.000 claims abstract description 17
- 238000012360 testing method Methods 0.000 claims description 38
- 238000001514 detection method Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 8
- 238000005070 sampling Methods 0.000 claims description 8
- 230000006870 function Effects 0.000 claims description 7
- 230000005284 excitation Effects 0.000 claims description 6
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 claims description 3
- 229910052739 hydrogen Inorganic materials 0.000 claims description 3
- 239000001257 hydrogen Substances 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 5
- 238000013461 design Methods 0.000 description 19
- 238000010586 diagram Methods 0.000 description 4
- 230000001133 acceleration Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 239000004020 conductor Substances 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000005415 magnetization Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
Images
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/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
- G01R33/56518—Correction of image distortions, e.g. due to magnetic field inhomogeneities due to eddy currents, e.g. caused by switching of the gradient magnetic field
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公开了一种磁共振成像***加速涡流矫正的方法和设备，包括：设置一组初始延迟时间；其中，一组初始延迟时间不具有时间轴上的连续性和完备性；利用涡流矫正序列采集初始延迟时间的涡流信号；计算涡流信号的时间补偿参数；基于时间补偿参数判断残余涡流是否小于涡流阈值；若小于，完成涡流矫正；否则利用延迟时间优化算法，在初始延迟时间和时间补偿参数的基础上，每次迭代计算一组新的延迟时间；每次基于一组新的延迟时间对当前轴进行涡流矫正，直至当次的残余涡流小于所述涡流阈值。本发明能够在保证测量精度的同时大幅降低涡流矫正的时间。
Description
技术领域
本发明属于磁共振成像技术领域，具体涉及一种磁共振成像***中加速涡流矫正的方法和设备。
背景技术
在磁共振成像(magnetic resonance imaging，MRI)过程中，梯度线圈中的电流随着时间的切换，会在周围导体结构中产生涡流，涡流会产生一个空间和时间上不断变化的磁场，使得成像区域内的磁场梯度发生畸变，影响MRI的图像质量。因此需要对磁共振成像***进行涡流矫正，现有技术中常见的涡流矫正方法是通过对脉冲序列发出的梯度波形进行预矫正，使实际产生的梯度场接近于理想波形的涡流矫正方法。具体如下：
如图1所示，其中RF(Radio Frequency)为射频脉冲时间轴，α为RF轴上射频脉冲的翻转角，Gdir为矫正涡流方向的梯度时间轴，G_test为Gdir轴上测量涡流的测试梯度，TE(Time of Echo)是回波时间，数据采集窗位于数据采集时间轴RX上。在测试梯度产生涡流后，通过改变产生涡流到序列接收涡流信号的时间，来检测不同延迟时间的涡流信号，进而确定当前涡流的幅度常数和时间常数数据，在经过多次延迟时间的参数扫描后，得到涡流场模型方程的时间常数和幅度常数，从而对涡流进行一次矫正。再次进行信号检测后，判断梯度实际波形与理想波形的差距达到要求，否则重复进行第二次矫正，直到涡流场补偿到要求状态。并采用同样的方式对梯度的其他轴进行矫正，以完成涡流的整体矫正。
然而，上述方案由于测试梯度一般调节到磁共振***的最大幅值，同时为了不引入额外的涡流导致之后处理的涡流模型复杂，进而使得涡流补偿参数估计不准确，测试梯度的上升时间和保持时间通常都需要在5秒左右，整个TR(Time of Repetition,重复时间)通常在10秒以上。则单个轴一次涡流的完整时间检测扫描时间将达到数小时，且每个轴上的涡流矫正通常需要重复多次，导致整个涡流矫正过程十分耗时。
发明内容
本发明的目的在于提供一种磁共振成像***中加速涡流矫正的方法和设备，用于解决现有技术中磁共振成像***的涡流矫正过程十分耗时的问题。
为了实现上述目的，本发明采用以下技术方案：
第一方面，本发明提供一种磁共振成像***中加速涡流矫正的方法，包括：
设置一组初始延迟时间；其中，一组所述初始延迟时间不具有时间轴上的连续性和完备性；
利用涡流矫正序列采集所述初始延迟时间的涡流信号；
计算所述涡流信号的时间补偿参数；
基于所述时间补偿参数判断残余涡流是否小于涡流阈值；
若小于，完成涡流矫正；否则利用延迟时间优化算法，在所述初始延迟时间和所述时间补偿参数的基础上，每次迭代计算一组新的延迟时间；
每次基于一组新的延迟时间对当前轴进行涡流矫正，直至当次的所述残余涡流小于所述涡流阈值。
在一种可能的设计中，相邻两个延迟时间之间的差值至少大于三倍的被扫描模体的纵向弛豫时间。
在一种可能的设计中，利用延迟时间优化算法，在所述初始延迟时间和所述时间补偿参数的基础上，每次迭代计算一组新的延迟时间，包括：
利用计算时间复杂度的log多分段原理和完备递归采样值原理，在所述初始延迟时间和所述时间补偿参数的基础上，迭代计算一组新的延迟时间，计算公式如下：
Γ＝{U(log(*,T))} (1)
其中，Tk表示第k次计算得到的一组延迟时间，i表示涡流拟合的项数，
表示关于时间补偿参数的函数关系，Γ表示延迟时间的优化算法，log(*,T)表示关于延迟时间任意数乘积的log优化算子，U表示在完整延迟时间集上的算子集合。
在一种可能的设计中，利用涡流矫正序列采集所述初始延迟时间的涡流信号，包括：
在一个重复时间TR内，利用自相关涡流矫正序列施加一次测试梯度Gtest和检测信号梯度在第一方向梯度轴Gdir上，并进行多次激发，以在数据采集时间轴RX采集到一组初始延时时间的涡流信号。
在一种可能的设计中，利用涡流矫正序列采集所述初始延迟时间的涡流信号，还包括：
在一个重复时间TR内，利用互相关涡流矫正序列分别施加测试梯度Gtest和检测信号梯度在第二方向梯度轴Gdir1和第三方向梯度轴Gdir2上，并进行多次激发，以对所述第二方向梯度轴Gdir1和所述第三方向梯度轴Gdir2轴的涡流进行互相关矫正。
在一种可能的设计中，当对所述测试梯度Gtest和检测信号梯度进行两次激发时，所述自相关涡流矫正序列和所述互相关涡流矫正序列的射频脉冲时间轴RF上的翻转角α0和翻转角α1的角度分别为90°和180°。
在一种可能的设计中，在对所述第一方向梯度轴Gdir、所述第二方向梯度轴Gdir1或所述第三方向梯度轴Gdir2进行涡流检测时，至少需要采集两个重复时间TR内的涡流信号，且两个重复时间TR内施加的测试梯度Gtest的极性相反。
在一种可能的设计中，利用涡流矫正序列采集所述初始延迟时间的涡流信号之后，还包括：
从涡流信号中分离计算所述测试梯度Gtest产生的梯度涡流信号，公式如下:
其中，Be中是涡流场强度，Φ0是所述自相关涡流矫正序列采集时的相位值；γ是氢核旋磁比，t表示时间，xs表示扫描点在第一方向梯度轴Gdir上的位置，Φ1(t)、Φ2(t)分别表示两次扫描得到的总相位值，Φ(t)表示涡流场产生的相位值。
在一种可能的设计中，计算所述涡流信号的时间补偿参数，包括：
计算在第一方向梯度轴Gdir上的涡流信号的梯度强度，计算公式如下：
其中，Gx(t)，Gy(t)，Gz(t)分别为第一方向梯度轴Gdir、第二方向梯度轴Gdir1和第三方向梯度轴Gdir2的理想波形，Gx'(t)为第一方向梯度轴Gdir上受涡流影响的实际波形，
为第i个时间常数对应的第一方向梯度轴Gdir自相关涡流补偿卷积核函数，
和
为第i个交叉项涡流补偿卷积核函数，即第二方向梯度轴Gdir1和第三方向梯度轴Gdir2对第一方向梯度轴Gdir的影响；其中，所述卷积核函数为：
其中，
和
分别是X轴梯度上涡流矫正的第i项幅度补偿参数和时间补偿参数。
第二方面，本发明提供一种计算机设备，包括：依次通信相连的存储器、处理器和收发器，其中，所述存储器用于存储计算机程序，所述收发器用于收发消息，所述处理器用于读取所述计算机程序，执行如第一方面所述的磁共振成像***中加速涡流矫正的方法。
第三方面，本发明提供了一种可读存储介质，所述可读存储介质上存储有指令，当所述指令在计算机上运行时，执行如上第一方面或第一方面中任意一种可能设计的所述磁共振成像***中加速涡流矫正的方法。
第四方面，本发明提供了一种包含指令的计算机程序产品，当所述指令在计算机上运行时，使所述计算机执行如上第一方面或第一方面中任意一种可能设计的所述磁共振成像***中加速涡流矫正的方法。
有益效果：
1.本发明通过设置一组初始延迟时间；然后利用涡流矫正序列采集所述初始延迟时间的涡流信号；通过计算所述涡流信号的时间补偿参数；然后基于所述时间补偿参数判断残余涡流是否小于涡流阈值；若小于，完成当前轴的涡流矫正；否则利用延迟时间优化算法，在所述初始延迟时间和所述时间补偿参数的基础上，迭代计算一组新的延迟时间；基于一组新的延迟时间对当前轴进行涡流矫正，直至所述残余涡流小于所述涡流阈值。其中，由于初始延迟时间是不连续和不完备的，且延迟时间优化算法，在所述初始延迟时间和所述时间补偿参数的基础上，迭代计算一组新的延迟时间，则每一组得到的延迟时间参数均是不连续和不完备的，则无需在时间轴上检测连续的涡流信号，能够在保证涡流的测量精度的同时，大幅降低涡流的矫正时间。
2.本发明通过应用计算时间复杂度的log多分段原理和完备递归采样值原理，在初始延迟时间和拟合得到的时间补偿参数的基础上，计算一组新的延迟时间；其中，基于log多分段原理和完备递归采样值原理获得的新的延迟时间本身也是不连续的，则通过迭代计算新的延迟时间并每次更新应用到新的一次涡流矫正中，能够进一步避免每次检测都需要在连续完备的时间上检测，使得涡流矫正时间大幅降低。
附图说明
图1为现有技术中的涡流矫正序列图；
图2为本发明提供的涡流检测序列的示意图；
图3为本发明提供的磁共振成像***中加速涡流矫正的方法的流程图；
图4(a)为本发明提供的加速涡流矫正方法的涡流信号与拟合信号的示意图；
图4(b)为现有涡流矫正方法的涡流信号与拟合信号的示意图。
具体实施方式
为使本说明书实施例的目的、技术方案和优点更加清楚，下面将结合本说明书实施例中的附图，对本说明书实施例中的技术方案进行清楚、完整地描述。基于本说明书中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本说明保护的范围。
实施例
为了解决现有的磁共振成像***的涡流矫正过程十分耗时的问题，如图2-4所示，第一方面，本发明提供一种磁共振成像***中加速涡流矫正的方法，包括但不限于由步骤S101～S106实现。具体如下：
步骤S101.设置一组初始延迟时间；其中，一组所述初始延迟时间不具有时间轴上的连续性和完备性；
在步骤S101中，需要说明的是，一组所述初始延迟时间中包含若干个由小到大的延迟时间，优选的，包含最大延迟时间(对应于长时涡流)和最小延迟时间(对应短时涡流)，优选的，相邻两个延迟时间之间的差值至少大于三倍的被扫描模体的纵向弛豫时间T1；其中，所述纵向弛豫时间T1是指纵向磁化强度恢复的时间常数。
步骤S102.利用涡流矫正序列采集所述初始延迟时间的涡流信号；
如图2所示，其中，需要说明的是，所述涡流矫正序列包括自相关涡流矫正序列和互相关涡流矫正序列，所述自相关涡流矫正序列和所述互相关涡流矫正序列均包括射频脉冲时间轴RF和数据采集时间轴RX，所述自相关涡流矫正序列还包括第一方向梯度轴Gdir，所述互相关涡流矫正序列还包括第二方向梯度轴Gdir1和第三方向梯度轴Gdir2。
其中，在一种可能的设计中，利用涡流矫正序列采集所述初始延迟时间的涡流信号包括：
在一个重复时间TR内，利用自相关涡流矫正序列施加一次测试梯度Gtest和检测信号梯度在第一方向梯度轴Gdir上，并进行多次激发，以在数据采集时间轴RX上采集到一组初始延时时间的涡流信号。
在一种可能的设计中，利用涡流矫正序列采集所述初始延迟时间的涡流信号，还包括：
在一个重复时间TR内，利用互相关涡流矫正序列分别施加测试梯度Gtest和检测信号梯度在第二方向梯度轴Gdir1和第三方向梯度轴Gdir2上，并进行多次激发，以对所述第二方向梯度轴Gdir1和所述第三方向梯度轴Gdir2轴的涡流进行互相关矫正。
其中，在一种可能的设计中，当对所述测试梯度Gtest和检测信号梯度进行两次激发时，所述自相关涡流矫正序列和所述互相关涡流矫正序列的射频脉冲时间轴RF上的翻转角α0和翻转角α1的角度分别为90°和180°。
在一种可能的设计中，在对所述第一方向梯度轴Gdir、所述第二方向梯度轴Gdir1或所述第三方向梯度轴Gdir2进行涡流检测时，至少需要采集两个重复时间TR内的涡流信号，且两个重复时间TR内施加的测试梯度Gtest的极性相反。
其中，在一种可能的设计中，利用涡流矫正序列采集所述初始延迟时间的涡流信号之后，还包括：
从涡流信号中分离计算所述测试梯度Gtest产生的梯度涡流信号，公式如下:
其中，Be中是涡流场强度，Φ0是所述自相关涡流矫正序列采集时的相位值；γ是氢核旋磁比，t表示时间，xs表示扫描点在第一方向梯度轴Gdir上的位置，Φ1(t)、Φ2(t)分别表示两次扫描得到的总相位值，Φ(t)表示涡流场产生的相位值。
步骤S103.计算所述涡流信号的时间补偿参数；
在一种可能的设计中，计算所述涡流信号的时间补偿参数，包括：
计算在第一方向梯度轴Gdir上的涡流信号的梯度强度，计算公式如下：
其中，Gx(t)，Gy(t)，Gz(t)分别为第一方向梯度轴Gdir、第二方向梯度轴Gdir1和第三方向梯度轴Gdir2的理想波形，Gx'(t)为第一方向梯度轴Gdir上受涡流影响的实际波形，
为第i个时间常数对应的第一方向梯度轴Gdir自相关涡流补偿卷积核函数，
和
为第i个交叉项涡流补偿卷积核函数，即第二方向梯度轴Gdir1和第三方向梯度轴Gdir2对第一方向梯度轴Gdir的影响；其中，所述卷积核函数为：
其中，
和
分别是X轴梯度上涡流矫正的第i项幅度补偿参数和时间补偿参数。
步骤S104.基于所述时间补偿参数判断残余涡流是否小于涡流阈值；
在步骤104中，需要说明的是，基于所述时间补偿参数可以对所述残余涡流进行描述，当所述残余涡流对磁场均匀度的影响小于涡流阈值对磁场均匀度的影响时，则说明满足涡流矫正的精度要求。
步骤S105.若小于，完成涡流矫正；否则利用延迟时间优化算法，在所述初始延迟时间和所述时间补偿参数的基础上，每次迭代计算一组新的延迟时间；
在一种可能的设计中，利用延迟时间优化算法，在所述初始延迟时间和所述时间补偿参数的基础上，迭代计算一组新的延迟时间，包括：
利用计算时间复杂度的log多分段原理和完备递归采样值原理，在所述初始延迟时间和所述时间补偿参数的基础上，迭代计算一组新的延迟时间，计算公式如下：
Γ＝{U(log(*,T))} (1)
其中，Tk表示第k次计算得到的一组延迟时间，i表示涡流拟合的项数，
表示关于时间补偿参数的函数关系，Γ表示延迟时间的优化算法，log(*,T)表示关于延迟时间任意数乘积的log优化算子，U表示在完整延迟时间集上的算子集合。
其中，应用计算时间复杂度log多分段原理和完备递归采样值原理的理由如下：
由于算法的时间复杂度是一个函数，它定性描述该算法的运行时间。常见的算法时间复杂度由小到大依次为：O(1)<O(log2n)<O(n)<O(nlog2n)<O(n2)<O(n3)<…<O(2n)<O(n！)；其中，O(1)表示基本语句的执行次数是一个常数的时间，O(log2n)、O(n)、O(nlog2n)、O(n2)和O(n3)称为多项式时间，而O(2n)和O(n！)称为指数时间。计算机科学家普遍认为越靠前的是效率越高的算法，因此简单来说算法中对数时间(log)级别的时间复杂度使用了分治思想，这个底数直接由分治的复杂度决定。例如如果采用二分法，那么就会以2为底数，三分法就会以3为底数，其他亦然，该算法的log级别的渐进意义，时间复杂度的增长与处理数据多少的增长的关系都是一样的。
因此，由于在涡流矫正中需要计算连续的延迟时间上的涡流信号，将计算机的时间复杂度的log多分段原理和完备递归采样值原理用来计算涡流检测加速序列的延迟时间组，能够得到不连续的多段延迟时间组，每次更新应用在涡流的多次矫正中，避免了每次检测都需要在连续完备的时间上检测，使得涡流矫正时间大幅降低。
步骤S106.每次基于一组新的延迟时间进行涡流矫正，直至当次的残余涡流小于所述涡流阈值。
其中，需要说明的是，利用延迟时间优化算法，可以逐次迭代计算出多组新的延迟时间，将每次计算得出的一组新的延迟时间重复步骤S102～步骤S104逐次进行涡流矫正，直到基于当前新的延迟时间得到的残余涡流小于涡流阈值，则完成涡流矫正。
基于上述公开的内容，通过设置一组初始延迟时间；然后利用涡流矫正序列采集所述初始延迟时间的涡流信号；通过计算所述涡流信号的时间补偿参数；然后基于所述时间补偿参数判断残余涡流是否小于涡流阈值；若小于，完成涡流矫正；否则利用延迟时间优化算法，在所述初始延迟时间和所述时间补偿参数的基础上，每次迭代计算一组新的延迟时间；每次基于一组新的延迟时间进行涡流矫正，直至当次的残余涡流小于所述涡流阈值。其中，由于初始延迟时间是不连续和不完备的，且利用延迟时间优化算法，在所述初始延迟时间和所述时间补偿参数的基础上，每次迭代计算一组新的延迟时间，则每一组得到的延迟时间参数均是不连续和不完备的，因此无需在时间轴上检测连续的涡流信号，能够在保证涡流的测量精度的同时，大幅降低涡流的矫正时间。
作为本实施例的一个应用例，以四项的涡流矫正因子为例(即公式(4)中，N＝4)。
测量涡流在时间0.02ms到4000ms的涡流大小，测量序列选择两次加速，首次设定测量的延迟时间参数组分为五对，分别为[(20，1000000)；(80000，1800000)；(160000，2500000)；(300000，3200000)；(600000，4000000)]，单位为us；利用如第一方面所述的检测涡流的加速序列，施加测试梯度Gtest和检测信号梯度在同一个梯度轴上，每一个TR时间内采集到一对涡流延迟时间的信号，在五个TR后，反转测试梯度的极性，再次采集；计算得到测试梯度的涡流信号，经过涡流场的描述公式(4)和(5)计算得到拟合的涡流补偿参数，其中四个时间补偿参数是[29999.199219，499989，626878.375，5068739.5]，这时计算残余涡流对磁场均匀度的影响是否小于0.2个ppm，第一次矫正后残余涡流对磁场均匀度的影响是2.667个ppm，不满足要求，则将这次涡流补偿参数的时间补偿参数以及上述的首次设定的延迟时间参数，应用计算时间复杂度的log多分段原理和完备递归采样值原理，根据上一次迭代设置的延迟时间参数以及当前拟合得到的时间补偿参数，计算一组新的延迟时间参数：[20.0，201062.542990901；11356.2145054558，795817.460312703；37880.4042463233，1604419.81682821；739603.954582226，2239864.86401948；110199.233558885，4000000]；利用基于新的延迟时间参数的加速序列再次采集信号。重复上述步骤，可得知在第四次矫正时，达到涡流补偿后的磁场均匀度0.177ppm，满足精度要求，则完成了当前轴的涡流自相关矫正。
然后，将测试涡流梯度分别施加在其余两个轴上，进行其他轴的涡流互相关矫正。通过在三个梯度轴上重复上述步骤，直到完成整体的涡流矫正。
如图4(a)和图4(b)所示，采用本实施例的加速涡流矫正方法，完成单个轴上的涡流校正时间在140个TR时间，约1500秒左右；如果采用现有的加速涡流矫正方法，在完整延迟时间上采集，经过涡流场的描述公式(4)和(5)计算得到拟合的涡流补偿参数，如果重复矫正次数同为四次时满足精度要求，则完成单个轴上的涡流校正时间最少需要1200个TR时间，约12850秒左右。可见，采用本实施例的加速涡流矫正方法进行涡流矫正能够大幅降低涡流矫正时间，提高矫正效率。
第二方面，本发明提供一种计算机设备，包括：依次通信相连的存储器、处理器和收发器，其中，所述存储器用于存储计算机程序，所述收发器用于收发消息，所述处理器用于读取所述计算机程序，执行如第一方面所述的磁共振成像***中加速涡流矫正的方法。
具体举例的，所述存储器可以但不限于包括随机存取存储器(Random-AccessMemory，RAM)、只读存储器(Read-Only Memory，ROM)、闪存(Flash Memory)、先进先出存储器(First Input First Output，FIFO)和/或先进后出存储器(First Input Last Output，FILO)等等；所述收发器可以但不限于为WiFi(无线保真)无线收发器、蓝牙无线收发器、GPRS(General Packet Radio Service，通用分组无线服务技术)无线收发器和/或ZigBee(紫蜂协议，基于IEEE802.15.4标准的低功耗局域网协议)无线收发器等；所述处理器可以不限于采用型号为STM32F105系列的微处理器。此外，所述网关设备还可以但不限于包括有电源模块、显示屏和其它必要的部件。
第三方面，本发明提供了一种可读存储介质，所述可读存储介质上存储有指令，当所述指令在计算机上运行时，执行如上第一方面或第一方面中任意一种可能设计的所述磁共振成像***中加速涡流矫正的方法。其中，所述可读存储介质是指存储数据的载体，可以但不限于包括软盘、光盘、硬盘、闪存、优盘和/或记忆棒(Memory Stick)等，所述计算机可以是通用计算机、专用计算机、计算机网络、或者其他可编程装置。
第四方面，本发明提供了一种包含指令的计算机程序产品，当所述指令在计算机上运行时，使所述计算机执行如上第一方面或第一方面中任意一种可能设计的所述磁共振成像***中加速涡流矫正的方法。其中，所述计算机可以是通用计算机、专用计算机、计算机网络、或者其他可编程装置
最后应说明的是：以上所述仅为本发明的优选实施例而已，并不用于限制本发明的保护范围。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (8)
1.一种磁共振成像***中加速涡流矫正的方法，其特征在于，包括：
设置一组初始延迟时间；其中，一组所述初始延迟时间不具有时间轴上的连续性和完备性；
利用涡流矫正序列采集所述初始延迟时间的涡流信号；
计算所述涡流信号的时间补偿参数；
基于所述时间补偿参数判断残余涡流是否小于涡流阈值；
若小于，完成涡流矫正；否则利用延迟时间优化算法，在所述初始延迟时间和所述时间补偿参数的基础上，每次迭代计算一组新的延迟时间；
每次基于一组新的延迟时间进行涡流矫正，直至当次的残余涡流小于所述涡流阈值，完成涡流矫正；
相邻两个延迟时间之间的差值至少大于三倍的被扫描模体的纵向弛豫时间；
在所述初始延迟时间和所述时间补偿参数的基础上，每次迭代计算一组新的延迟时间，包括：
利用计算时间复杂度的log多分段原理和完备递归采样值原理，在所述初始延迟时间和所述时间补偿参数的基础上，每次迭代计算一组新的延迟时间，计算公式如下：
Γ＝{U(log(*,T))} (1)
其中，Tk表示第k次计算得到的一组延迟时间，i表示涡流拟合的项数，
表示关于时间补偿参数的函数关系，Γ表示延迟时间的优化算法，log(*,T)表示关于延迟时间任意数乘积的log优化算子，U表示在完整延迟时间集上的算子集合。
2.根据权利要求1所述的方法，其特征在于，利用涡流矫正序列采集所述初始延迟时间的涡流信号，包括：
在一个重复时间TR内，利用自相关涡流矫正序列施加一次测试梯度Gtest和检测信号梯度在第一方向梯度轴Gdir上，并进行多次激发，以在数据采集时间轴RX上采集到一组初始延时时间的涡流信号。
3.根据权利要求2所述的方法，其特征在于，利用涡流矫正序列采集所述初始延迟时间的涡流信号，还包括：
在一个重复时间TR内，利用互相关涡流矫正序列分别施加测试梯度Gtest和检测信号梯度在第二方向梯度轴Gdir1和第三方向梯度轴Gdir2上，并进行多次激发，以对所述第二方向梯度轴Gdir1和所述第三方向梯度轴Gdir2轴的涡流进行互相关矫正。
4.根据权利要求3所述的方法，其特征在于，当对所述测试梯度Gtest和检测信号梯度进行两次激发时，所述自相关涡流矫正序列和所述互相关涡流矫正序列的射频脉冲时间轴RF上的翻转角α0和翻转角α1的角度分别为90°和180°。
5.根据权利要求3所述的方法，其特征在于，在对所述第一方向梯度轴Gdir、所述第二方向梯度轴Gdir1或所述第三方向梯度轴Gdir2进行涡流检测时，至少需要采集两个重复时间TR内的涡流信号，且两个重复时间TR内施加的测试梯度Gtest的极性相反。
6.根据权利要求5所述的方法，其特征在于，利用涡流矫正序列采集所述初始延迟时间的涡流信号之后，还包括：
从涡流信号中分离计算所述测试梯度Gtest产生的梯度涡流信号，公式如下:
其中，Be中是涡流场强度，Φ0是所述自相关涡流矫正序列采集时的相位值；γ是氢核旋磁比，t表示时间，xs表示扫描点在第一方向梯度轴Gdir上的位置，Φ1(t)、Φ2(t)分别表示两次扫描得到的总相位值，Φ(t)表示涡流场产生的相位值。
7.根据权利要求3所述的方法，其特征在于，计算所述涡流信号的时间补偿参数，包括：
计算在第一方向梯度轴Gdir上的涡流信号的梯度强度，计算公式如下：
其中，Gx(t)，Gy(t)，Gz(t)分别为第一方向梯度轴Gdir、第二方向梯度轴Gdir1和第三方向梯度轴Gdir2的理想波形，Gx'(t)为第一方向梯度轴Gdir上受涡流影响的实际波形，
为第i个时间常数对应的第一方向梯度轴Gdir自相关涡流补偿卷积核函数，
和
为第i个交叉项涡流补偿卷积核函数，即第二方向梯度轴Gdir1和第三方向梯度轴Gdir2对第一方向梯度轴Gdir的影响；其中，所述卷积核函数为：
其中，
和
分别是X轴梯度上涡流矫正的第i项幅度补偿参数和时间补偿参数。
8.一种计算机设备，其特征在于，包括：依次通信相连的存储器、处理器和收发器，其中，所述存储器用于存储计算机程序，所述收发器用于收发消息，所述处理器用于读取所述计算机程序，执行如权利要求1-7任意一项所述的磁共振成像***中加速涡流矫正的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110587406.4A CN113219391B (zh) | 2021-05-27 | 2021-05-27 | 一种磁共振成像***中加速涡流矫正的方法和设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110587406.4A CN113219391B (zh) | 2021-05-27 | 2021-05-27 | 一种磁共振成像***中加速涡流矫正的方法和设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113219391A CN113219391A (zh) | 2021-08-06 |
| CN113219391B true CN113219391B (zh) | 2021-12-14 |
Family
ID=77098930
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110587406.4A Active CN113219391B (zh) | 2021-05-27 | 2021-05-27 | 一种磁共振成像***中加速涡流矫正的方法和设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113219391B (zh) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102298129A (zh) * | 2011-05-30 | 2011-12-28 | 苏州安科医疗***有限公司 | 一种用于核磁共振成像***的涡流测量及补偿方法 |
| CN102866373A (zh) * | 2011-12-12 | 2013-01-09 | 中国科学院深圳先进技术研究院 | 磁共振温度成像中温度测量的矫正方法和*** |
| EP3640661A1 (en) * | 2018-10-16 | 2020-04-22 | Siemens Healthcare GmbH | Method and control unit for active compensation of eddy current induced magnetic fields in magnetic resonance imaging |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0323841A (ja) * | 1989-06-20 | 1991-01-31 | Sanyo Electric Co Ltd | Ｍｒｉ装置の勾配磁場電源回路 |
| CN102096054B (zh) * | 2010-12-22 | 2013-03-13 | 华东师范大学 | 一种快速自旋回波脉冲序列的梯度***延时校正方法 |
| CN106483482B (zh) * | 2015-08-25 | 2019-08-23 | 上海联影医疗科技有限公司 | 用于磁共振成像***的梯度涡流校正方法和装置 |
| CN108279393A (zh) * | 2017-01-05 | 2018-07-13 | 上海康达卡勒幅医疗科技有限公司 | 一种磁共振成像全自动预加重校正涡流的方法 |
| CN109799472B (zh) * | 2019-03-08 | 2022-06-28 | 沈阳工业大学 | 一种基于深度学习的磁共振涡流补偿方法 |
-
2021
- 2021-05-27 CN CN202110587406.4A patent/CN113219391B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102298129A (zh) * | 2011-05-30 | 2011-12-28 | 苏州安科医疗***有限公司 | 一种用于核磁共振成像***的涡流测量及补偿方法 |
| CN102866373A (zh) * | 2011-12-12 | 2013-01-09 | 中国科学院深圳先进技术研究院 | 磁共振温度成像中温度测量的矫正方法和*** |
| EP3640661A1 (en) * | 2018-10-16 | 2020-04-22 | Siemens Healthcare GmbH | Method and control unit for active compensation of eddy current induced magnetic fields in magnetic resonance imaging |
Non-Patent Citations (1)
| Title |
|---|
| 《永磁磁共振扩散加权平面回波成像涡流伪影校正》;王艳飞等;《东北大学学报(自然科学版)》;20200715 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113219391A (zh) | 2021-08-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20110304333A1 (en) | Means And Methods For Providing High Resolution MRI | |
| US9046590B2 (en) | Magnetic resonance imaging apparatus phase correction using one or more prescans with identical readout and slice gradients | |
| US9316713B2 (en) | Magnetic resonance imaging apparatus and RF pulse for navigator and imaging sequence applying method | |
| JP2000037367A (ja) | 磁場のシミング | |
| WO2021248752A1 (zh) | 一种用于磁共振成像的梯度涡流补偿方法及*** | |
| CN106872920B (zh) | 磁共振成像***射频校准方法及装置 | |
| CN110215210A (zh) | 磁共振成像方法、***、计算机设备和可读存储介质 | |
| CN110320485B (zh) | 测量梯度延时和一阶场不均匀度的装置、方法及存储介质 | |
| CN113219391B (zh) | 一种磁共振成像***中加速涡流矫正的方法和设备 | |
| US6470203B2 (en) | MR imaging method, phase error measuring method and MRI system | |
| CN115629347B (zh) | 一种磁共振成像***中获得梯度轨迹的方法、装置和介质 | |
| CN104155621B (zh) | 一种准确测量静磁场b0分布的方法 | |
| CN107479015B (zh) | 磁共振校准扫描序列配置和图像获取的方法及*** | |
| US10429471B2 (en) | MRI with variable density sampling | |
| CN110558984B (zh) | 射频特定吸收率的确定方法、检测装置、磁共振***及扫描控制方法 | |
| CN115905817B (zh) | 一种多通道非均匀场信号降噪方法、装置和存储介质 | |
| CN111830450A (zh) | 确定和消除磁共振设备中射频脉冲与选层梯度之间的时延的方法 | |
| JP4456009B2 (ja) | Ｎｍｒ測定方法および装置 | |
| CN110068780A (zh) | 一种磁共振成像的自适应方法、装置和磁共振成像*** | |
| US20020050816A1 (en) | MR imaging method, phase error measuring method, and MRI apparatus | |
| Andris et al. | Optimized measurement of magnetic field maps using nuclear magnetic resonance (NMR) | |
| CN115951281B (zh) | 磁共振梯度延迟和梯度切换率同时检测方法、装置和介质 | |
| JPH05237067A (ja) | 磁気共鳴イメージング装置 | |
| CN114167331A (zh) | 一种非均匀磁场下测量纵向弛豫时间的方法和设备 | |
| CN106597335A (zh) | 电压校准方法、磁共振成像方法及*** |
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 | ||
| CP03 | Change of name, title or address |
Address after: 610000 No. 104, floor 1, building 2, No. 88, Keyuan South Road, high tech Zone, Chengdu, Sichuan Patentee after: Chengdu Mingshi Junzhi Technology Co.,Ltd. Country or region after: China Address before: 610000 No. 104, floor 1, building 2, No. 88, Keyuan South Road, high tech Zone, Chengdu, Sichuan Patentee before: Chengdu mingshijunzhi Medical Technology Co.,Ltd. Country or region before: China |
|
| CP03 | Change of name, title or address |