CN119375799A - 一种超低场磁共振成像方法、系统、电子设备及产品 - Google Patents
一种超低场磁共振成像方法、系统、电子设备及产品 Download PDFInfo
- Publication number
- CN119375799A CN119375799A CN202411500525.1A CN202411500525A CN119375799A CN 119375799 A CN119375799 A CN 119375799A CN 202411500525 A CN202411500525 A CN 202411500525A CN 119375799 A CN119375799 A CN 119375799A
- Authority
- CN
- China
- Prior art keywords
- gradient
- spin echo
- sequence
- readout
- magnetic resonance
- 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.)
- Granted
Links
- 238000002595 magnetic resonance imaging Methods 0.000 title claims abstract description 46
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000000264 spin echo pulse sequence Methods 0.000 claims description 114
- 238000012937 correction Methods 0.000 claims description 26
- 238000004590 computer program Methods 0.000 claims description 12
- 230000001351 cycling effect Effects 0.000 claims description 6
- 230000006870 function Effects 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 4
- 238000004891 communication Methods 0.000 claims description 3
- 239000010410 layer Substances 0.000 description 21
- 238000003384 imaging method Methods 0.000 description 11
- 230000000694 effects Effects 0.000 description 7
- 238000013461 design Methods 0.000 description 6
- 238000000819 phase cycle Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 230000006698 induction Effects 0.000 description 4
- 239000000243 solution Substances 0.000 description 4
- 238000012935 Averaging Methods 0.000 description 2
- 230000004087 circulation Effects 0.000 description 2
- 230000001934 delay Effects 0.000 description 2
- 230000008030 elimination Effects 0.000 description 2
- 238000003379 elimination reaction Methods 0.000 description 2
- 238000000685 Carr-Purcell-Meiboom-Gill pulse sequence Methods 0.000 description 1
- 238000005481 NMR spectroscopy Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000002872 contrast media Substances 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 239000008155 medical solution Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000013421 nuclear magnetic resonance imaging Methods 0.000 description 1
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method 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/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/56572—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of a gradient magnetic field, e.g. non-linearity of a gradient magnetic field
-
- 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/4818—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
-
- 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/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
Landscapes
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General 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)
- Nonlinear Science (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明属于磁共振成像技术领域，其目的在于提供一种超低场磁共振成像方法、系统、电子设备及产品。本发明通过迭代的方式对不同辐条的各梯度轴进行读出预散相梯度Grp的面积及读出梯度Gr的梯度延迟的补偿，使采用Radial采集方法读出的两梯度轴梯度的延迟差尽量小，使k空间中心尽量对齐而减少偏移，以减少因梯度性能差异带来的影响，从而提高图像质量；同时，本发明改进了成像序列，通过对快速自旋回波序列进行修正，可在超低场能获的更好的成像效果，使得平均次数最少可设置为1，从而可以采集更密的辐条，在提高分辨率同时不需要两次循环，由此可在确保图像分辨率和质量满足要求的前提下，缩短扫描时长。
Description
技术领域
本发明属于磁共振成像技术领域，具体涉及一种超低场磁共振成像方法、系统、电子设备及产品。
背景技术
磁共振成像技术是利用待测对象(如人体组织)中某种原子核的核磁共振现象，将所得的磁共振信号经过电子计算机处理，以重建出人体某一层面的图像的诊断技术。目前，在核磁共振成像技术的应用中，超低场磁共振成像因其成本低、安全性高、便携性强和对屏蔽要求低等优势，在需要经济高效医疗解决方案、特殊人群成像、安全使用对比剂以及便携医疗设备等场景中得以广泛应用。
由于低磁场强度和环境干扰等限制，在对待测对象进行超低场磁共振成像时，主磁场(B0)的均匀性通常较差，而这会影响最终生成图像的质量。为此，现有技术中，通常使用快速自旋回波类序列来获得质量更好的图像，同时，在进行超低场磁共振成像数据采集时，还结合Radial采集方法(一种非笛卡尔坐标系的采样方法，其通过径向轨迹来获取k空间数据)增加对运动的不敏感性，以降低图像中因待测对象运动而出现的伪影。
但是，在使用现有技术过程中，发明人发现现有技术中至少存在如下问题：
在采用Radial采集方法时，由于读出梯度会随着辐条的旋转角度而有所变化，不同梯度轴和相同梯度轴的不同极性下的梯度性能(如梯度延迟和剩磁)都有所差异，使得在超低场磁共振成像中应用Radial序列进行数据采集成像的效果不佳；另外，常规的自旋回波类序列通常利用CPMG(Carr-Purcell-Meiboom-Gill，一种脉冲序列)条件来减弱振铃和FID(Free Induction Decay，自由感应衰减)造成的伪影，需要至少进行两次平均才能达到理想的成像效果，导致扫描耗时较长。
发明内容
本发明旨在至少在一定程度上解决上述技术问题，本发明提供了一种超低场磁共振成像方法、系统、电子设备及产品。
为了实现上述目的，本发明采用以下技术方案：
第一方面，本发明提供了一种超低场磁共振成像方法，包括：
采用指定快速自旋回波序列对待测对象进行预扫描；其中，对所述待测对象进行预扫描时，不开启层编码梯度和层编码梯度重绕梯度；
获取由所述待测对象产生的各梯度轴的多个回波信号，并根据各梯度轴的多个回波信号，得到与多个所述回波信号一一对应的多个回波辐条图像；
从各梯度轴的多个回波辐条图像中分别获取各梯度轴的第一回波辐条图像和第二回波辐条图像，并根据各梯度轴的第一回波辐条图像和第二回波辐条图像，得到各梯度轴的读出预散相梯度的面积补偿量及读出梯度的梯度延迟补偿量；
根据各梯度轴的读出预散相梯度的面积补偿量及读出梯度的梯度延迟补偿量，对所述指定快速自旋回波序列进行修正，得到修正后快速自旋回波序列；
将所述修正后快速自旋回波序列作为新的指定快速自旋回波序列，以便重新采用所述指定快速自旋回波序列对待测对象进行预扫描，直到各梯度轴的读出预散相梯度的面积补偿量及读出梯度的梯度延迟补偿量均低于预设阈值时，将得到的修正后快速自旋回波序列作为最终快速自旋回波序列；
采用所述最终快速自旋回波序列对待测对象进行扫描，以便得到所述待测对象的磁共振图像；其中，采用所述最终快速自旋回波序列对待测对象进行扫描时，开启层编码梯度和层编码梯度重绕梯度。
在一个可能的设计中，所述指定快速自旋回波序列采用3D-FSE序列，所述指定快速自旋回波序列的ETL方向为层编码方向。
在一个可能的设计中，采用指定快速自旋回波序列对待测对象进行预扫描时，将两次不同的相位循环序列应用于相邻的两个辐条。
在一个可能的设计中，采用Radial采集方法获取所述待测对象产生的由多个回波信号，且读出预散相梯度的幅值与读出梯度的幅值保持一致。
在一个可能的设计中，任一梯度轴的读出预散相梯度的面积补偿量均包括正极面积补偿量和负极面积补偿量，任一梯度轴的读出梯度的梯度延迟补偿量包括正极梯度延迟补偿量和负极梯度延迟补偿量；所述任一梯度轴为X轴时，其正极面积补偿量为：
式中，Δx为+X轴的图像分辨率；N为+X方向的总读出点数；phase(·)为取相位操作函数；C1(n)为+X方向读出点n对应的第一回波辐条图像，n＝1,2,3,…,N-1；C2(n)为+X方向读出点n对应的第二回波辐条图像；
所述任一梯度轴为X轴时，其正极梯度延迟补偿量为：
式中，dw为+X方向的读出点采样间隔。
在一个可能的设计中，根据各梯度轴的读出预散相梯度的面积补偿量及读出梯度的梯度延迟补偿量，对所述指定快速自旋回波序列进行修正，得到修正后快速自旋回波序列，包括：
获取所述指定快速自旋回波序列中的辐条旋转角度；
根据所述指定快速自旋回波序列中的辐条旋转角度、各梯度轴的读出预散相梯度的面积补偿量和读出梯度的梯度延迟补偿量，得到各梯度轴的修正后读出预散相梯度和修正后梯度延迟补偿量；
根据各梯度轴的修正后读出预散相梯度和修正后梯度延迟补偿量对所述指定快速自旋回波序列进行修正，得到修正后快速自旋回波序列。
在一个可能的设计中，根据各梯度轴的修正后读出预散相梯度和修正后梯度延迟补偿量对所述指定快速自旋回波序列进行修正，得到修正后快速自旋回波序列，包括：
根据各梯度轴的修正后读出预散相梯度和修正后梯度延迟补偿量对所述指定快速自旋回波序列进行初始修正，得到初始修正后快速自旋回波序列；
对所述初始修正后快速自旋回波序列进行二次修正，得到修正后快速自旋回波序列。
第二方面，本发明提供了一种超低场磁共振成像系统，用于实现如上述任意一项所述的超低场磁共振成像方法；所述超低场磁共振成像系统包括：
磁共振序列修正模块，用于采用指定快速自旋回波序列对待测对象进行预扫描；其中，对所述待测对象进行预扫描时，不开启层编码梯度和层编码梯度重绕梯度；
所述磁共振序列修正模块，还用于获取由所述待测对象产生的各梯度轴的多个回波信号，并根据各梯度轴的多个回波信号，得到与多个所述回波信号一一对应的多个回波辐条图像；
所述磁共振序列修正模块，还用于从各梯度轴的多个回波辐条图像中分别获取各梯度轴的第一回波辐条图像和第二回波辐条图像，并根据各梯度轴的第一回波辐条图像和第二回波辐条图像，得到各梯度轴的读出预散相梯度的面积补偿量及读出梯度的梯度延迟补偿量；
所述磁共振序列修正模块，还用于根据各梯度轴的读出预散相梯度的面积补偿量及读出梯度的梯度延迟补偿量，对所述指定快速自旋回波序列进行修正，得到修正后快速自旋回波序列；
所述磁共振序列修正模块，还用于将所述修正后快速自旋回波序列作为新的指定快速自旋回波序列，以便重新采用所述指定快速自旋回波序列对待测对象进行预扫描，直到各梯度轴的读出预散相梯度的面积补偿量及读出梯度的梯度延迟补偿量均低于预设阈值时，将得到的修正后快速自旋回波序列作为最终快速自旋回波序列；
磁共振成像模块，与所述磁共振序列修正模块通信连接，用于采用所述最终快速自旋回波序列对待测对象进行扫描，以便得到所述待测对象的磁共振图像；其中，采用所述最终快速自旋回波序列对待测对象进行扫描时，开启层编码梯度和层编码梯度重绕梯度。
第三方面，本发明提供了一种电子设备，包括：
存储器，用于存储计算机程序指令；以及，
处理器，用于执行所述计算机程序指令从而完成如上述任意一项所述的超低场磁共振成像方法的操作。
第四方面，本发明提供了一种计算机程序产品，包括计算机程序或指令，所述计算机程序或所述指令在被计算机执行时实现如上述任意一项所述的超低场磁共振成像方法。
本发明的有益效果为：
本发明公开了一种超低场磁共振成像方法、系统、电子设备及产品，可提高成像质量，同时可缩短磁共振扫描耗时。具体地，本发明在实施过程中，通过预扫描来了解不同辐条在各梯度轴上的梯度性能差异，并基于预扫描结果，通过迭代的方式对不同辐条的各梯度轴进行读出预散相梯度Grp的面积及读出梯度Gr的梯度延迟的补偿，使采用Radial采集方法读出的两梯度轴梯度的延迟差尽量小，使k空间中心尽量对齐而减少偏移，以减少因梯度性能差异带来的影响，从而提高图像质量；同时，本发明改进了成像序列，通过对快速自旋回波序列进行修正，可在超低场能获的更好的成像效果，使得平均次数最少可设置为1，从而可以采集更密的辐条，在提高分辨率同时不需要两次循环，由此可在确保图像分辨率和质量满足要求的前提下，缩短扫描时长。
本发明的其他有益效果将在具体实施方式中进一步进行说明。
附图说明
图1是实施例中一种超低场磁共振成像方法的流程图；
图2是实施例中一种超低场磁共振成像系统的模块框图。
具体实施方式
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将结合附图和实施例或现有技术的描述对本发明作简单地介绍，显而易见地，下面关于附图结构的描述仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在此需要说明的是，对于这些实施例方式的说明用于帮助理解本发明，但并不构成对本发明的限定。
实施例1：
本实施例公开了一种超低场磁共振成像方法，可以但不限于由具有一定计算资源的计算机设备或虚拟机执行，例如由个人计算机、智能手机、个人数字助理或可穿戴设备等电子设备执行，或者由虚拟机执行。
如图1所示，一种超低场磁共振成像方法，可以但不限于包括有如下步骤：
S1.采用指定快速自旋回波序列对待测对象进行预扫描；其中，对所述待测对象进行预扫描时，不开启层编码梯度Gse和层编码梯度重绕梯度Gsr，也即在对待测对象进行预扫描的过程中没有层间编码的梯度，由此可确保扫描只专注于读出方向的优化，而不受层编码的影响。
本实施例中，在层选方向(指通过施加梯度场来选择成像层面的方向，用于确定成像层面的位置和厚度，根据不同的梯度场方向，可以实现不同类型的层面选择)为Z轴，读出方向(也即相位编码梯度的方向，相位编码梯度可实现对图像的空间分辨率和质量的控制)为X轴和Y轴时，读出预散相梯度Grp和读出梯度Gr向X轴和Y轴的分配方式为：
Grp_x＝Grp*cos(θ)，Grp_y＝Grp*sin(θ)；
Gr_x＝Gr*cos(θ)，Gr_y＝Gr*sin(θ)；
式中，θ为指定的辐条旋转角度；Grp_x为读出预散相梯度Grp向X轴的分配量；Grp_y为读出预散相梯度Grp向Y轴的分配量；Gr_x为读出梯度Gr向X轴的分配量；Gr_y为读出梯度Gr向Y轴的分配量。
步骤S1中，所述指定快速自旋回波序列采用3D-FSE(3D-Fast Spin Echo，三维快速自旋回波)序列，所述指定快速自旋回波序列的ETL(Echo Train Length，回波链长度)方向为层编码方向，也即回波链的方向与层编码方向一致，以提高成像效率和质量。需要说明的是，本实施例中，所述指定快速自旋回波序列采用3D-FSE序列，可利于在较短的时间内获取高质量的图像。
本实施例改进了成像序列，通过对快速自旋回波序列相位循环方式进行改进，可在超低场能获的更好的成像效果，使得平均次数最少可设置为1，从而可以采集更密的辐条，在提高分辨率同时不需要两次循环，由此可在确保图像分辨率和质量满足要求的前提下，缩短扫描时长。具体地，步骤S1中，采用指定快速自旋回波序列对待测对象进行预扫描时，将两次不同的相位循环序列应用于相邻的两个辐条。
传统的CPMG序列的相位循环序列可以是[+x,+y,+y,+y,…]和[+x,-y,-y,-y,…]，通过两次的相位循环序列能够消除回波脉冲中的自由感应衰减(FID)信号和振铃信号，从而减少图像伪影，然而，采用此种相位循环方式需要至少2次平均(即两次相位循环)，采集时间较长。
本实施例中，通过优化所述指定快速自旋回波序列中相位循环的设置，将两次不同的相位循环序列(如[+x,+y,+y,+y,…]和[+x,-y,-y,-y,…])分别应用到相邻的两个辐条，使得本实施例可以在仅进行一次平均的情况下消除伪影，也即可采用相位循环技术消除回聚脉冲的自由感应衰减信号以及振铃信号，仍能保证伪影消除的效果，同时可减少采集时间，使采集时间缩短一半。只要保证相邻的两个辐条有不同的相位循环序列，就能够实现伪影的消除，而不必重复多次采集。
S2.获取由所述待测对象产生的各梯度轴的多个回波信号，并根据各梯度轴的多个回波信号，得到与多个所述回波信号一一对应的多个回波辐条图像；应当理解的是，多个回波信号可称之为回波链，每个回波信号对应一个时间点的k空间数据，多个回波信号会填充k空间中的数据。
步骤S2中，采用Radial采集方法获取所述待测对象产生的由多个回波信号，且读出预散相梯度Grp的幅值与读出梯度Gr的幅值保持一致，以确保相位编码的一致性和减少伪影。
需要说明的是，在步骤S1中，采用指定快速自旋回波序列对待测对象进行预扫描时，分别在多个特定的梯度方向进行扫描，多个特定的梯度方向如包括+X方向(沿X轴正方向)、X方向(沿X轴负方向)、+Y方向(沿Y轴正方向)和Y方向(沿Y轴负方向)四个梯度方向，此时每个梯度方向都会产生一组回波数据；以+X方向为例，对+X方向上产生的一组回波数据进行傅里叶逆变换，即可以得到一组对应于该方向的回波辐条图像。对应地，步骤S2中，得到的由所述待测对象产生的任一回波信号中包括多个梯度方向的回波数据，该任一回波信号对应的回波图像即包括多个梯度方向的回波辐条图像。
需要说明的是，对多个所述回波信号分别进行傅里叶变换处理，即可得到与多个所述回波信号一一对应的多个回波辐条图像，其中，每一回波辐条图像代表一个特定时间点的特定梯度方向的磁共振图像，这些数据共同参与了回波图像的重建。
S3.从各梯度轴的多个回波辐条图像中分别获取各梯度轴的第一回波辐条图像和第二回波辐条图像，并根据各梯度轴的第一回波辐条图像和第二回波辐条图像，得到各梯度轴的读出预散相梯度Grp的面积补偿量及读出梯度Gr的梯度延迟补偿量；需要说明的是，所述第一回波辐条图像和所述第二回波辐条图像分别为基于多个回波信号中的第一个回波信号和第二个回波信号生成的回波辐条图像，其中，所述第一回波信号是在最初的时间点采集的，其通常具有最高的信噪比和最低的衰减效应，所述第二回波信号是在接近第一个回波后短时间内采集的，通常也具有较高的图像质量，所述第一回波辐条图像和所述第二回波辐条图像的线性相位差和图像质量通常较好，二者的线性相位差反映了读出方向磁场不均匀、读出预散相梯度编码面积不准确及其他系统误差带来的相位变化。
本实施例中，各梯度轴的读出预散相梯度Grp的面积补偿量可用于补偿由于各梯度轴的读出梯度Gr所导致的信号相位误差，从而提高最终图像的质量。梯度延迟是指信号在不同时间点到达检测器的时间差异，通过计算这些延迟，可以补偿由于各梯度轴造成的时间误差，减少成像伪影。
具体地，步骤S3中，任一梯度轴的读出预散相梯度Grp的面积补偿量均包括正极面积补偿量和负极面积补偿量，任一梯度轴的读出梯度Gr的梯度延迟补偿量包括正极梯度延迟补偿量和负极梯度延迟补偿量；所述任一梯度轴为X轴时，其正极面积补偿量为：
式中，Δx为+X轴的图像分辨率；N为+X方向的总读出点数；phase(·)为取相位操作函数；C1(n)为+X方向读出点n对应的第一回波辐条图像，n＝1,2,3,…,N-1；C2(n)为+X方向读出点n对应的第二回波辐条图像；
所述任一梯度轴为X轴时，其正极梯度延迟补偿量为：
式中，dw为+X方向的读出点采样间隔。
需要说明的是，基于上述步骤，可以计算出+X方向、-X方向、+Y方向和-Y方向时的读出预散相梯度Grp的面积补偿量Δkrpx+,Δkrpx-,Δkrpy+和Δkrpy-以及其读出梯度Gr的梯度延迟补偿量Δτx+,Δτx-,Δτy+和Δτy-。
S4.根据各梯度轴的读出预散相梯度Grp的面积补偿量及读出梯度Gr的梯度延迟补偿量，对所述指定快速自旋回波序列进行修正，得到修正后快速自旋回波序列。
步骤S4中，根据各梯度轴的读出预散相梯度Grp的面积补偿量及读出梯度Gr的梯度延迟补偿量，对所述指定快速自旋回波序列进行修正，得到修正后快速自旋回波序列，包括：
S401.获取所述指定快速自旋回波序列中的辐条旋转角度；
S402.根据所述指定快速自旋回波序列中的辐条旋转角度、各梯度轴的读出预散相梯度Grp的面积补偿量和读出梯度Gr的梯度延迟补偿量，得到各梯度轴的修正后读出预散相梯度和修正后梯度延迟补偿量；
具体地，所述梯度轴包括X轴和Y轴，对应地，X轴的修正后读出预散相梯度面积为：
Y轴的修正后读出预散相梯度为：
X轴的修正后梯度延迟补偿量为：
Y轴的修正后梯度延迟补偿量为：
式中，θ为指定辐条旋转角度；Arp为理想的读出预散相梯度面积；sign(·)为符号函数，自变量为正时取1，为负时取-1。
S403.根据各梯度轴的修正后读出预散相梯度和修正后梯度延迟补偿量对所述指定快速自旋回波序列进行修正，得到修正后快速自旋回波序列。
本实施例的步骤S403中，根据各梯度轴的修正后读出预散相梯度和修正后梯度延迟补偿量对所述指定快速自旋回波序列进行修正，得到修正后快速自旋回波序列，包括：
根据各梯度轴的修正后读出预散相梯度和修正后梯度延迟补偿量对所述指定快速自旋回波序列进行初始修正，得到初始修正后快速自旋回波序列；
对所述初始修正后快速自旋回波序列进行二次修正，得到修正后快速自旋回波序列。
具体地，本实施例中，对所述初始修正后快速自旋回波序列进行二次修正时，先保持无编码状态，采用所述初始修正后快速自旋回波序列对待测对象进行预扫描，以得到所有辐条角度下的数据，而后再次计算得到各梯度轴的读出预散相梯度Grp的面积补偿量及读出梯度Gr的梯度延迟补偿量，并将计算得到的各梯度轴的读出预散相梯度Grp的面积补偿量直接按照余弦、正弦关系分配给X轴和Y轴，不再区分极性，将计算得到的读出梯度Gr的梯度延迟补偿量直接同等地赋给X轴和Y轴，而不考虑极性，由此得到修正后快速自旋回波序列。
S5.将所述修正后快速自旋回波序列作为新的指定快速自旋回波序列，以便重新采用所述指定快速自旋回波序列对待测对象进行预扫描，直到各梯度轴的读出预散相梯度Grp的面积补偿量及读出梯度Gr的梯度延迟补偿量均低于预设阈值时，将得到的修正后快速自旋回波序列作为最终快速自旋回波序列。
S6.采用所述最终快速自旋回波序列对待测对象进行扫描，以便得到所述待测对象的磁共振图像；其中，采用所述最终快速自旋回波序列对待测对象进行扫描时，开启层编码梯度Gse和层编码梯度重绕梯度Gsr。
本实施例可提高成像质量，同时可缩短磁共振扫描耗时。具体地，本实施例在实施过程中，通过预扫描来了解不同辐条在各梯度轴上的梯度性能差异，并基于预扫描结果，通过迭代的方式对不同辐条的各梯度轴进行读出预散相梯度Grp的面积及读出梯度Gr的梯度延迟的补偿，使采用Radial采集方法读出的两梯度轴梯度的延迟差尽量小，使k空间中心尽量对齐而减少偏移，以减少因梯度性能差异带来的影响，从而提高图像质量；同时，本实施例改进了成像序列，通过对快速自旋回波序列相位循环方式进行改进，可在超低场能获的更好的成像效果，使得平均次数最少可设置为1，从而可以采集更密的辐条，在提高分辨率同时不需要两次循环，由此可在确保图像分辨率和质量满足要求的前提下，缩短扫描时长。
实施例2：
本实施例公开了一种超低场磁共振成像系统，用于实现实施例1中超低场磁共振成像方法；如图2所示，所述超低场磁共振成像系统包括：
磁共振序列修正模块，用于采用指定快速自旋回波序列对待测对象进行预扫描；其中，对所述待测对象进行预扫描时，不开启层编码梯度和层编码梯度重绕梯度；
所述磁共振序列修正模块，还用于获取由所述待测对象产生的各梯度轴的多个回波信号，并根据各梯度轴的多个回波信号，得到与多个所述回波信号一一对应的多个回波辐条图像；
所述磁共振序列修正模块，还用于从各梯度轴的多个回波辐条图像中分别获取各梯度轴的第一回波辐条图像和第二回波辐条图像，并根据各梯度轴的第一回波辐条图像和第二回波辐条图像，得到各梯度轴的读出预散相梯度的面积补偿量及读出梯度的梯度延迟补偿量；
所述磁共振序列修正模块，还用于根据各梯度轴的读出预散相梯度的面积补偿量及读出梯度的梯度延迟补偿量，对所述指定快速自旋回波序列进行修正，得到修正后快速自旋回波序列；
所述磁共振序列修正模块，还用于将所述修正后快速自旋回波序列作为新的指定快速自旋回波序列，以便重新采用所述指定快速自旋回波序列对待测对象进行预扫描，直到各梯度轴的读出预散相梯度的面积补偿量及读出梯度的梯度延迟补偿量均低于预设阈值时，将得到的修正后快速自旋回波序列作为最终快速自旋回波序列；
磁共振成像模块，与所述磁共振序列修正模块通信连接，用于采用所述最终快速自旋回波序列对待测对象进行扫描，以便得到所述待测对象的磁共振图像；其中，采用所述最终快速自旋回波序列对待测对象进行扫描时，开启层编码梯度和层编码梯度重绕梯度。
需要说明的是，本实施例2提供的超低场磁共振成像系统的工作过程、工作细节和技术效果，可以参见实施例1，于此不再赘述。
实施例3：
在实施例1或2的基础上，本实施例公开了一种电子设备，该设备可以是智能手机、平板电脑、笔记本电脑或者台式电脑等。电子设备可能被称为用户终端、便携式终端、台式终端等，电子设备包括：
存储器，用于存储计算机程序指令；以及，
处理器，用于执行所述计算机程序指令从而完成如实施例1中任一所述的超低场磁共振成像方法的操作。
实施例4：
在实施例1至3任意一项实施例的基础上，本实施例公开了一种计算机程序产品，包括计算机程序或指令，所述计算机程序或所述指令在被计算机执行时实现如实施例1中任意一项所述的超低场磁共振成像方法。
显然，本领域的技术人员应该明白，上述的本发明的各模块或各步骤可以用通用的计算装置来实现，它们可以集中在单个的计算装置上，或者分布在多个计算装置所组成的网络上，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它们存储在存储装置中由计算装置来执行，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样，本发明不限制于任何特定的硬件和软件结合。
最后应说明的是，以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换。而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (10)
1.一种超低场磁共振成像方法，其特征在于，包括：
采用指定快速自旋回波序列对待测对象进行预扫描；其中，对所述待测对象进行预扫描时，不开启层编码梯度和层编码梯度重绕梯度；
获取由所述待测对象产生的各梯度轴的多个回波信号，并根据各梯度轴的多个回波信号，得到与多个所述回波信号一一对应的多个回波辐条图像；
从各梯度轴的多个回波辐条图像中分别获取各梯度轴的第一回波辐条图像和第二回波辐条图像，并根据各梯度轴的第一回波辐条图像和第二回波辐条图像，得到各梯度轴的读出预散相梯度的面积补偿量及读出梯度的梯度延迟补偿量；
根据各梯度轴的读出预散相梯度的面积补偿量及读出梯度的梯度延迟补偿量，对所述指定快速自旋回波序列进行修正，得到修正后快速自旋回波序列；
将所述修正后快速自旋回波序列作为新的指定快速自旋回波序列，以便重新采用所述指定快速自旋回波序列对待测对象进行预扫描，直到各梯度轴的读出预散相梯度的面积补偿量及读出梯度的梯度延迟补偿量均低于预设阈值时，将得到的修正后快速自旋回波序列作为最终快速自旋回波序列；
采用所述最终快速自旋回波序列对待测对象进行扫描，以便得到所述待测对象的磁共振图像；其中，采用所述最终快速自旋回波序列对待测对象进行扫描时，开启层编码梯度和层编码梯度重绕梯度。
2.根据权利要求1所述的一种超低场磁共振成像方法，其特征在于，所述指定快速自旋回波序列采用3D-FSE序列，所述指定快速自旋回波序列的ETL方向为层编码方向。
3.根据权利要求1所述的一种超低场磁共振成像方法，其特征在于，采用指定快速自旋回波序列对待测对象进行预扫描时，将两次不同的相位循环序列应用于相邻的两个辐条。
4.根据权利要求1所述的一种超低场磁共振成像方法，其特征在于，采用Radial采集方法获取所述待测对象产生的由多个回波信号，且读出预散相梯度的幅值与读出梯度的幅值保持一致。
5.根据权利要求1所述的一种超低场磁共振成像方法，其特征在于，任一梯度轴的读出预散相梯度的面积补偿量均包括正极面积补偿量和负极面积补偿量，任一梯度轴的读出梯度的梯度延迟补偿量包括正极梯度延迟补偿量和负极梯度延迟补偿量；所述任一梯度轴为X轴时，其正极面积补偿量为：
；
式中，为+X轴的图像分辨率；N为+X方向的总读出点数；为取相位操作函数；为+X方向读出点n对应的第一回波辐条图像，n=1,2,3,…,N-1；为+X方向读出点n对应的第二回波辐条图像；
所述任一梯度轴为X轴时，其正极梯度延迟补偿量为：
；
式中，为+X方向的读出点采样间隔。
6.根据权利要求1所述的一种超低场磁共振成像方法，其特征在于，根据各梯度轴的读出预散相梯度的面积补偿量及读出梯度的梯度延迟补偿量，对所述指定快速自旋回波序列进行修正，得到修正后快速自旋回波序列，包括：
获取所述指定快速自旋回波序列中的辐条旋转角度；
根据所述指定快速自旋回波序列中的辐条旋转角度、各梯度轴的读出预散相梯度的面积补偿量和读出梯度的梯度延迟补偿量，得到各梯度轴的修正后读出预散相梯度和修正后梯度延迟补偿量；
根据各梯度轴的修正后读出预散相梯度和修正后梯度延迟补偿量对所述指定快速自旋回波序列进行修正，得到修正后快速自旋回波序列。
7.根据权利要求6所述的一种超低场磁共振成像方法，其特征在于，根据各梯度轴的修正后读出预散相梯度和修正后梯度延迟补偿量对所述指定快速自旋回波序列进行修正，得到修正后快速自旋回波序列，包括：
根据各梯度轴的修正后读出预散相梯度和修正后梯度延迟补偿量对所述指定快速自旋回波序列进行初始修正，得到初始修正后快速自旋回波序列；
对所述初始修正后快速自旋回波序列进行二次修正，得到修正后快速自旋回波序列。
8.一种超低场磁共振成像系统，其特征在于，用于实现如权利要求1至7中任意一项所述的超低场磁共振成像方法；所述超低场磁共振成像系统包括：
磁共振序列修正模块，用于采用指定快速自旋回波序列对待测对象进行预扫描；其中，对所述待测对象进行预扫描时，不开启层编码梯度和层编码梯度重绕梯度；
所述磁共振序列修正模块，还用于获取由所述待测对象产生的各梯度轴的多个回波信号，并根据各梯度轴的多个回波信号，得到与多个所述回波信号一一对应的多个回波辐条图像；
所述磁共振序列修正模块，还用于从各梯度轴的多个回波辐条图像中分别获取各梯度轴的第一回波辐条图像和第二回波辐条图像，并根据各梯度轴的第一回波辐条图像和第二回波辐条图像，得到各梯度轴的读出预散相梯度的面积补偿量及读出梯度的梯度延迟补偿量；
所述磁共振序列修正模块，还用于根据各梯度轴的读出预散相梯度的面积补偿量及读出梯度的梯度延迟补偿量，对所述指定快速自旋回波序列进行修正，得到修正后快速自旋回波序列；
所述磁共振序列修正模块，还用于将所述修正后快速自旋回波序列作为新的指定快速自旋回波序列，以便重新采用所述指定快速自旋回波序列对待测对象进行预扫描，直到各梯度轴的读出预散相梯度的面积补偿量及读出梯度的梯度延迟补偿量均低于预设阈值时，将得到的修正后快速自旋回波序列作为最终快速自旋回波序列；
磁共振成像模块，与所述磁共振序列修正模块通信连接，用于采用所述最终快速自旋回波序列对待测对象进行扫描，以便得到所述待测对象的磁共振图像；其中，采用所述最终快速自旋回波序列对待测对象进行扫描时，开启层编码梯度和层编码梯度重绕梯度。
9.一种电子设备，其特征在于，包括：
存储器，用于存储计算机程序指令；以及，
处理器，用于执行所述计算机程序指令从而完成如权利要求1至7中任意一项所述的超低场磁共振成像方法的操作。
10.一种计算机程序产品，包括计算机程序或指令，其特征在于，所述计算机程序或所述指令在被计算机执行时实现如权利要求1至7中任意一项所述的超低场磁共振成像方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202411500525.1A CN119375799B (zh) | 2024-10-25 | 2024-10-25 | 一种超低场磁共振成像方法、系统、电子设备及产品 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202411500525.1A CN119375799B (zh) | 2024-10-25 | 2024-10-25 | 一种超低场磁共振成像方法、系统、电子设备及产品 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN119375799A true CN119375799A (zh) | 2025-01-28 |
| CN119375799B CN119375799B (zh) | 2025-04-15 |
Family
ID=94335386
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202411500525.1A Active CN119375799B (zh) | 2024-10-25 | 2024-10-25 | 一种超低场磁共振成像方法、系统、电子设备及产品 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN119375799B (zh) |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1058123A1 (en) * | 1999-06-03 | 2000-12-06 | General Electric Company | Fast spin echo phase correction for MRI system |
| US20090072826A1 (en) * | 2007-09-13 | 2009-03-19 | The Board Of Trustees Of The Leland Stanford Junior University | Magnetic resonance imaging with bipolar multi-echo sequences |
| CN107024670A (zh) * | 2016-01-29 | 2017-08-08 | 上海联影医疗科技有限公司 | 磁共振系统的校正方法及装置 |
| US20180188343A1 (en) * | 2016-12-30 | 2018-07-05 | Uih America, Inc. | System and method for magnetic resonance imaging |
| WO2018184056A1 (en) * | 2017-04-05 | 2018-10-11 | The University Of Queensland | Magnetic resonance imaging method and apparatus |
| CN111352055A (zh) * | 2020-03-17 | 2020-06-30 | 无锡鸣石峻致医疗科技有限公司 | 一种前瞻性相位校正平面回波成像技术 |
| CN113093077A (zh) * | 2021-04-08 | 2021-07-09 | 无锡鸣石峻致医疗科技有限公司 | 3d-dess序列的优化方法、3d-dess优化序列和磁共振成像的匀场方法 |
| US20220206095A1 (en) * | 2019-04-12 | 2022-06-30 | Baycrest Centre For Geriatric Care | Dual gradient echo and spin echo magnetic resonance fingerprinting for simultaneous estimation of t1, t2, and t2* with integrated b1 correction |
| CN115639510A (zh) * | 2022-10-09 | 2023-01-24 | 深圳市联影高端医疗装备创新研究院 | 磁共振成像方法、波谱成像方法、装置、设备和存储介质 |
| US20230168324A1 (en) * | 2021-11-30 | 2023-06-01 | GE Precision Healthcare LLC | Magnetic resonance imaging system and method, and computer-readable storage medium |
-
2024
- 2024-10-25 CN CN202411500525.1A patent/CN119375799B/zh active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1058123A1 (en) * | 1999-06-03 | 2000-12-06 | General Electric Company | Fast spin echo phase correction for MRI system |
| US20090072826A1 (en) * | 2007-09-13 | 2009-03-19 | The Board Of Trustees Of The Leland Stanford Junior University | Magnetic resonance imaging with bipolar multi-echo sequences |
| CN107024670A (zh) * | 2016-01-29 | 2017-08-08 | 上海联影医疗科技有限公司 | 磁共振系统的校正方法及装置 |
| US20180188343A1 (en) * | 2016-12-30 | 2018-07-05 | Uih America, Inc. | System and method for magnetic resonance imaging |
| WO2018184056A1 (en) * | 2017-04-05 | 2018-10-11 | The University Of Queensland | Magnetic resonance imaging method and apparatus |
| US20220206095A1 (en) * | 2019-04-12 | 2022-06-30 | Baycrest Centre For Geriatric Care | Dual gradient echo and spin echo magnetic resonance fingerprinting for simultaneous estimation of t1, t2, and t2* with integrated b1 correction |
| CN111352055A (zh) * | 2020-03-17 | 2020-06-30 | 无锡鸣石峻致医疗科技有限公司 | 一种前瞻性相位校正平面回波成像技术 |
| CN113093077A (zh) * | 2021-04-08 | 2021-07-09 | 无锡鸣石峻致医疗科技有限公司 | 3d-dess序列的优化方法、3d-dess优化序列和磁共振成像的匀场方法 |
| US20230168324A1 (en) * | 2021-11-30 | 2023-06-01 | GE Precision Healthcare LLC | Magnetic resonance imaging system and method, and computer-readable storage medium |
| CN115639510A (zh) * | 2022-10-09 | 2023-01-24 | 深圳市联影高端医疗装备创新研究院 | 磁共振成像方法、波谱成像方法、装置、设备和存储介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN119375799B (zh) | 2025-04-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108375746B (zh) | 一种相位反卷折方法与设备 | |
| CN102772210B (zh) | 弥散加权磁共振成像 | |
| Stehning et al. | Fast isotropic volumetric coronary MR angiography using free‐breathing 3D radial balanced FFE acquisition | |
| US7696753B2 (en) | MRI apparatus and method with B0 corrections using reference and object average units | |
| CN112014782B (zh) | 磁共振成像方法、装置及计算机存储介质 | |
| WO2020019412A1 (zh) | 磁共振成像方法、装置、设备及存储介质 | |
| CN103260510B (zh) | 磁共振成像装置以及对比度增强图像取得方法 | |
| US8538115B2 (en) | Coil compression for three dimensional autocalibrating parallel imaging with cartesian sampling | |
| CN103477238A (zh) | 采用来自先验采集的约束的压缩感测mr图像重建 | |
| CN108459289B (zh) | 一种基于数据一致性的多次激发弥散加权磁共振成像方法 | |
| CN108020796B (zh) | 一种磁共振扩散加权成像方法和装置 | |
| US11194003B2 (en) | MRI apparatus and MRI method | |
| CN113391251A (zh) | 磁共振图像重建方法、装置和设备 | |
| CN112272535A (zh) | 磁共振弥散成像的方法和系统 | |
| CN115639510A (zh) | 磁共振成像方法、波谱成像方法、装置、设备和存储介质 | |
| CN110109036B (zh) | 二维时空编码多扫磁共振成像非笛卡尔采样及重建方法 | |
| Saju et al. | Suppressing image blurring of PROPELLER MRI via untrained method | |
| CN110226100A (zh) | 用于磁共振成像的系统和方法 | |
| CN111257811B (zh) | 磁共振成像方法、装置及系统 | |
| CN119375799B (zh) | 一种超低场磁共振成像方法、系统、电子设备及产品 | |
| CN109696646B (zh) | 一种快速自旋回波脉冲序列中射频脉冲的优化方法 | |
| CN118799235B (zh) | 一种磁共振图像畸变矫正方法、装置、设备及介质 | |
| US10823808B2 (en) | Image processing apparatus, magnetic resonance imaging apparatus, and image generating method | |
| JP2020175194A (ja) | 磁気共鳴イメージング装置及び磁気共鳴イメージング方法 | |
| CN117503100A (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 |