CN113030814B - 一种磁共振k空间数据校正的方法 - Google Patents

一种磁共振k空间数据校正的方法 Download PDF

Info

Publication number
CN113030814B
CN113030814B CN202110253687.XA CN202110253687A CN113030814B CN 113030814 B CN113030814 B CN 113030814B CN 202110253687 A CN202110253687 A CN 202110253687A CN 113030814 B CN113030814 B CN 113030814B
Authority
CN
China
Prior art keywords
data
phase
ceil
gain
difference factor
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
Application number
CN202110253687.XA
Other languages
English (en)
Other versions
CN113030814A (zh
Inventor
吴林
刘杭
余洁
张涛
尧德中
解玺洁
胡霞飞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN202110253687.XA priority Critical patent/CN113030814B/zh
Publication of CN113030814A publication Critical patent/CN113030814A/zh
Application granted granted Critical
Publication of CN113030814B publication Critical patent/CN113030814B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/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/58—Calibration of imaging systems, e.g. using test probes, Phantoms; Calibration objects or fiducial markers such as active or passive RF coils surrounding an MR active material

Landscapes

  • Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本发明公开一种磁共振K空间数据校正的方法，应用于磁共振K空间数据处理领域，针对现有技术缺少K空间数据拼接的具体实现方法，本发明对拼接后的K空间数据进行增益与相位补偿，具体的：将三个不同增益参数下的K空间数据通过增益差校正与补偿，归一化到同一个接收增益水平；再对增益补偿后的K空间数据进行相位差校正与补偿。

Description

一种磁共振K空间数据校正的方法
技术领域
本发明属于计算机技术领域，特别涉及一种磁共振K空间数据处理技术。
背景技术
在MRIK空间中所有相位编码行的回波信号中，每个相位编码行的回波信号幅值存在明显差异，较大的相位编码行比较小的相位编码行信号幅值小，因此在增益设置过程中，在保证整个K空间所有相位编码行回波幅值最大的点不溢出的前提下设置固定的接收增益，会致使整体回波信号达不到最大的信号动态范围。Mark A等人和C.H.OH等人发现在不改动磁共振接收机电路的情况下，将K空间不同区域的相位编码行设置为不同的接收增益参数，既保证了相位编码较小行所在区域的回波信号不溢出，同时又使得较大相位编码行所在区域的回波信号达到更大的动态范围。
1998年，Mark A等人设置了6个接收增益参数，每个接收增益参数按照6dB递增顺序设置，在每个接收增益参数下扫描一次序列，每次得到一组K空间数据，然后将得到的6个K空间拼接成一个K空间数据，通过傅里叶变换得到一个3D图像，通过对拼接K空间数据重建得到的图像与最低增益下得到的K空间数据重建的图像对比，图像质量得到提升。该文献仅提及到需要一个校正系统来校正不同增益参数的差异，但是未详细阐述接收增益差校正的实现方法。2010年，C.H.OH等人提出通过优化K空间各相位编码行的接收增益参数来实现对MRI图像信噪比的提升，该方法首先要从所有相位编码行的回波信号中搜索出的最大值M(n)，并确定每个相位编码行的最佳接收增益
Figure BDA0002966938180000011
其中Sm是整个K空间信号的最大值，n为相位编码行的行号。通过仿真发现，接收机在固定接收增益下，ADC的最佳采样位宽最少需要19位，而在可变接收增益下，ADC的最佳采样位宽只需要大于15位即可满足，对ADC采样位宽需求明显降低。由于该文进行实验的MRI平台不支持在序列运行期间动态更改接收增益参数，因此该文实际实验范式设计参考了Mark A等人的思路，在每个接收增益参数设置下得到一组K空间数据，然后将不同接收增益参数下得到的K空间数据进行拼接，在拼接时提到需要将不同接收增益参数下的K空间数据归一化到同一接收接收增益参数下，但是也未提及具体的实现方法。
鉴于目前大部分磁共振平台在序列运行期间，接收增益参数不能动态更新，要使K空间整体达到最大的信号动态范围，如果沿用上述技术提出的接收增益参数设置方法和序列多次扫描的方式，即在每个接收增益参数下扫描一次序列，得到一组K空间原始数据，然后将不同接收增益参数下得到的K空间数据拼接成一幅K空间数据，在拼接时需要将不同接收增益参数下的K空间数据归一化到同一接收增益参数下。但是上述技术未详细阐述K空间数据拼接方法，也没有详细阐述磁共振接收增益差和相位差的校正问题。
发明内容
为解决上述技术问题，本发明提出一种磁共振K空间数据校正的方法，基于增益差校正和补偿方法将K空间不同接收增益参数下得到的数据归一化到最大接收增益参数下的幅值。
本发明采用的技术方案为：一种磁共振K空间数据校正的方法，包括：
将K空间的相位编码行行号表示为m，K空间一共A个相位编码行；
分别在增益参数为：τ1、τ2、τ3下得到各自对应的K空间数据：
Figure BDA0002966938180000021
τ1＞τ2＞τ3；
将
Figure BDA0002966938180000022
分成3个区域，
Figure BDA0002966938180000023
的第一区域数据记为
Figure BDA0002966938180000024
Figure BDA0002966938180000025
的第二区域数据记为
Figure BDA0002966938180000026
Figure BDA0002966938180000027
的第三区域数据记为
Figure BDA0002966938180000028
将
Figure BDA0002966938180000029
分成3个区域，
Figure BDA00029669381800000210
的第一区域数据记为
Figure BDA00029669381800000211
Figure BDA00029669381800000212
的第二区域数据记为
Figure BDA00029669381800000213
Figure BDA00029669381800000214
的第三区域数据记为
Figure BDA00029669381800000215
将
Figure BDA00029669381800000216
分成3个区域，
Figure BDA00029669381800000217
的第一区域数据记为
Figure BDA00029669381800000218
Figure BDA00029669381800000219
的第二区域数据记为
Figure BDA00029669381800000220
Figure BDA00029669381800000221
的第三区域数据记为
Figure BDA00029669381800000222
其中：
Figure BDA00029669381800000223
Figure BDA00029669381800000224
应的相位编码行行号为：m＝1到m＝A/6，m＝5A/6+1到m＝A；
Figure BDA00029669381800000225
对应的相位编码行行号为：m＝A/6+1到m＝A/3，m＝2A/3+1到m＝5A/6；
Figure BDA00029669381800000226
对应的相位编码行行号为：m＝A/3+1到m＝2A/3；
拼接后的数据空间为K，记K的第1区域数据为K1，K的第2区域数据为K2，K的第3区域数据为K3；将
Figure BDA00029669381800000227
的
Figure BDA00029669381800000228
作为K的K1，
Figure BDA00029669381800000229
的
Figure BDA00029669381800000230
作为K的K2，
Figure BDA00029669381800000231
的
Figure BDA00029669381800000232
作为K的K3；
通过增益差校正与补偿，将三组K空间数据归一化到同一个接收增益水平；具体的：让K2每行每一个数据幅值与第一增益差因子相乘得到K2*，让K3每行每一个数据幅值与第二增益差因子相乘得到K3，则K1、K2*、K3*构成增益补偿后的K；
对增益补偿后的K进行相位差校正和补偿。
对增益补偿后的K进行相位差校正和补偿，具体的：在增益补偿后，K1数据点相位保持不变；K2*每行数据点的相位加第一相位差因子得到K2*@；K3*每行数据点的相位加第二相位差因子得到K3*@。
第一增益差因子计算过程为：在
Figure BDA0002966938180000031
中，每次顺序取出相同4行的数据，取相位编码行行号为ceil(m/4)*4的中间N个数据的模定义为
Figure BDA0002966938180000032
令：
Figure BDA0002966938180000033
Figure BDA0002966938180000034
Figure BDA0002966938180000035
即为第一增益差因子。
第二增益差因子计算过程为：在
Figure BDA0002966938180000036
与
Figure BDA0002966938180000037
中，每次顺序取出相同4行的数据，取相位编码行行号为ceil(m/4)*4的中间N个数据的模定义为
Figure BDA0002966938180000038
令：
Figure BDA0002966938180000039
Figure BDA00029669381800000310
Figure BDA00029669381800000311
即为
Figure BDA00029669381800000312
和
Figure BDA00029669381800000313
的对应相位编码行的增益差因子ΔG2(ceil(m/4)*4)；ΔG1(ceil(m/4)*4)*ΔG2(ceil(m/4)*4)则第二增益差因子。
第一增益差因子计算过程为：在
Figure BDA00029669381800000314
中，每次顺序取出相同4行的数据，取相位编码行行号为ceil(m/4)*4的中间N个数据的模定义为
Figure BDA00029669381800000315
令：
Figure BDA00029669381800000316
Figure BDA00029669381800000317
Figure BDA00029669381800000318
为
Figure BDA00029669381800000319
和
Figure BDA00029669381800000320
对应相位编码行的增益差因子ΔG1(ceil(m/4)*4)，在
Figure BDA00029669381800000321
Figure BDA00029669381800000322
中，将这些区域所有行的增益差因子ΔG1(ceil(m/4)*4)求和再取平均值得到avg[ΔG1(ceil(m/4)*4)]，avg[ΔG1(ceil(m/4)*4)]为第一增益差因子。
第二增益差因子的计算过程为：
在
Figure BDA00029669381800000323
与
Figure BDA00029669381800000324
中，每次顺序取出相同4行的数据，取相位编码行行号为ceil(m/4)*4的中间N个数据的模定义为
Figure BDA00029669381800000325
令：
Figure BDA0002966938180000041
Figure BDA0002966938180000042
Figure BDA0002966938180000043
即为
Figure BDA0002966938180000044
与
Figure BDA0002966938180000045
对应相位编码行的增益差因子ΔG2(ceil(m/4)*4)，将这些区域所有行的ΔG2(ceil(m/4)*4)求和再取平均值得到avg[ΔG2(ceil(m/4)*4)]；avg[ΔG2(ceil(m/4)*4)]即为第二增益差因子。
第一相位差因子的计算过程为：在
Figure BDA0002966938180000046
中，每次顺序取出相同4行的数据，取相位编码行行号为ceil(m/4)*4的中间N个数据的相位定义为
Figure BDA0002966938180000047
Figure BDA0002966938180000048
令：
Figure BDA0002966938180000049
对γ1i进行出现频次统计，出现频次最多的值记为第一相位差因子ΔΦ1(ceil(m/4)*4)。
第二相位差因子的计算过程为：在
Figure BDA00029669381800000410
中，每次顺序取出相同4行的数据，取相位编码行行号为ceil(m/4)*4的中间N个数据的相位定义为
Figure BDA00029669381800000411
令：
Figure BDA00029669381800000412
对γ2i进行出现频次统计，出现频次最多的值则为
Figure BDA00029669381800000413
和
Figure BDA00029669381800000414
对应相位编码行的相位差因子ΔΦ2(ceil(m/4)*4)；ΔΦ1(ceil(m/4)*4)+ΔΦ2(ceil(m/4)*4)记为第二相位差因子。
本发明的有益效果：本发明分别在三个不同接收增益下得到的三组数据，通过增益差和相位差的校正和补偿，将三组K空间数据归一化到同一个接收增益水平；然后再同一接收增益水平下进行拼接，从而得到拼接后的K空间数据；
虽然本发明专利基于目前大部分磁共振平台在序列运行期间接收增益参数不能动态更新的背景下提出的，但是本发明专利提出的方法同样也适用于序列运行时接收增益参数能实时更新的磁共振平台，序列扫描期间，基于正常的扫描矩阵参数，在K空间不同区域实时动态地更新接收增益参数，扫描得到一幅完整的K空间数据，然后基于本发明提出的增益差校正和补偿方法将K空间不同接收增益参数下得到的数据归一化到最大接收增益参数下的幅值；采用本发明的方法可以使K空间整体达到最大的信号动态范围。
附图说明
图1位本发明实施例提供的系统执行流程图；
图2为本发明实施例提供的方案一的K空间数据拼接及增益校正示意图；
图3为本发明实施例提供的方案二的K空间数据拼接及增益校正示意图。
具体实施方式
为便于本领域技术人员理解本发明的技术内容，下面结合附图对本发明内容进一步阐释。
在序列扫描数据采集中，每个接收增益控制参数调整一次，则执行一次序列扫描，获得到一组K空间原始数据，获得接收增益参数分别为τdB、τ-6dB、τ-12dB的三组K空间数据Kτ、Kτ-6、Kτ-12。通过分析这三组数据的增益差和相位差，构建一个拼接的图像的数据空间，将拼接后的数据空间通过傅里叶变换可以获得一幅图像，系统执行流程如图1所示。
本发明提出了一种将不同接收增益参数下得到的K空间数据拼接成一幅K空间数据的方法，以及将不同接收增益参数下的K空间数据归一化到同一接收增益参数下的方法。虽然本发明基于目前大部分磁共振平台在序列运行期间接收增益参数不能动态更新的背景下提出的，但是本发明提出的方法同样也适用于序列运行时接收增益参数能实时更新的磁共振平台，序列扫描期间，基于正常的扫描矩阵参数，在K空间不同区域实时动态地更新接收增益参数，扫描得到一幅完整的K空间数据，然后基于本发明提出的增益差校正和补偿方法将K空间不同接收增益参数下得到的数据归一化到最大接收增益参数下的幅值。
本实施例以增益参数τdB、τ-6dB、τ-12dB为例进行说明：
分别在接收增益参数为τdB、τ-6dB、τ-12dB下得到的三组K空间数据Kτ、Kτ-6、Kτ-12。
将K空间的相位编码行行号表示为m，K空间一共A个相位编码行。为了在同一幅图像中同时描述3个接收增益的表达值，首先将Kτ分成3个区域分别用Kτ 1，Kτ 2，Kτ 3来表示，其中：Kτ 1为：m＝1到m＝A/6，m＝5A/6+1到m＝A；Kτ 2为：m＝A/6+1到m＝A/3，m＝2A/3+1到m＝5A/6；Kτ 3为：m＝A/3+1到m＝2A/3。利用同样的方式，对Kτ-6、Kτ-12也分成3个区域。
K空间所有区域在同一接收增益下，回波信号的幅值满足如下准则：
K() 1<K() 2<K() 3
其中，K() 1、K() 2、K() 3分别表示K()的第1、2、3区域回波信号幅值；这里的K()用于代表某个接收增益，比如本实施例中的Kτ或Kτ-6或Kτ-12。
拼接后数据空间，K() 1的回波信号在最大增益下得到最大的动态范围，同时K() 3的回波信号不溢出，并且K() 3使用最小增益下的数据填充。
为了将三个K空间的数据归一化到同一个接收增益水平，需要对数据进行增益差和相位差的校正。针对增益差和相位差的校正和补偿问题，本发明提出两种可选的实现方案，分别如下文方案一和方案二所述。
方案一：
如图2所示为方案一的K空间数据拼接及增益校正示意图，具体包括以下内容：
数据拼接方法：
拼接后的K空间数据的第1区域数据由Kτ第1区域数据填充；拼接后的K空间数据的第2区域数据由Kτ-6第2区域数据填充；拼接后的K空间数据的第3区域数据由Kτ-12第3区域数据填充。
增益校正：
在Kτ的Kτ 2和Kτ 3与Kτ-6的Kτ-6 2和Kτ-6 3中，每次顺序取出相同4行的数据，取相位编码行行号为ceil(m/4)*4的中间N个数据的模定义为δτi，δ(τ-6)i，i＝1,2,…,N。令：
Figure BDA0002966938180000061
Figure BDA0002966938180000062
本领域的人员应注意α1i在这里仅作为便于计算的中间变量，无具体的物理含义。
βτ,τ-6即为Kτ和Kτ-6对应相位编码行的增益差因子ΔG1(ceil(m/4)*4)；在Kτ-6的Kτ-6 3与Kτ-12的Kτ-12 3中，每次顺序取出相同4行的数据，取相位编码行行号为ceil(m/4)*4的中间N个数据的模定义为δ(τ-6)i，δ(τ-12)i，i＝1,2,…,N。令：
Figure BDA0002966938180000063
Figure BDA0002966938180000064
本领域的人员应注意α2i在这里仅作为便于计算的中间变量，无具体的物理含义。
βτ-6,τ-12即为Kτ-6 3和Kτ-12 3的对应相位编码行的增益差因子ΔG2(ceil(m/4)*4)。ΔG1(ceil(m/4)*4)*ΔG2(ceil(m/4)*4)则为在Kτ 3与Kτ-12 3中对应相位编码行的增益差因子。
增益补偿：
假定拼接后的数据空间为K，将Kτ的Kτ 1作为K的K1。Kτ-6的Kτ-6 2作为K的K2，让K2每行每一个数据幅值与ΔG1(ceil(m/4)*4)相乘得到K2*。Kτ-12的Kτ-12 3作为K的K3，让K3每行每一个数据幅值与ΔG1(ceil(m/4)*4)*ΔG2(ceil(m/4)*4)相乘得到K3*。则K1，K2*，K3*构成的K空间为补偿后的K空间。
相位校正：
在Kτ的Kτ 2和Kτ 3与Kτ-6的Kτ-6 2和Kτ-6 3中，每次顺序取出相同4行的数据，取相位编码行行号为ceil(m/4)*4的中间N个数据的相位定义为ετi，ε(τ-6)i，i＝1,2,…,N令：
γ1i＝ετi-ε(τ-6)i,i＝1,2,…,N
对γ1i进行出现频次统计，出现频次最多的值则是Kτ和Kτ-6对应相位编码行的相位差因子ΔΦ1(ceil(m/4)*4)；在Kτ-6的Kτ-6 3与Kτ-12的Kτ-12 3中，每次顺序取出相同4行的数据，取相位编码行行号为ceil(m/4)*4的中间N个数据的相位定义为ε(τ-6)i，ε(τ-12)i，i＝1,2,…,N令：
γ2i＝ε(τ-6)i-ε(τ-12)i,i＝1,2,…,N
对γ2i进行出现频次统计，出现频次最多的值则为Kτ-6 3和Kτ-12 3对应相位编码行的相位差因子ΔΦ2(ceil(m/4)*4)；ΔΦ1(ceil(m/4)*4)+ΔΦ2(ceil(m/4)*4)则为在Kτ 3与Kτ-12 3中相位差因子。
相位补偿：
在增益补偿后，K1数据点相位保持不变；K2*每行数据点的相位加
Figure BDA0002966938180000071
得到K2*@；K3*每行数据点的相位加ΔΦ1(ceil(m/4)*4)+ΔΦ2(ceil(m/4)*4)得到K3*@。
如图3所示为方案一的K空间数据拼接及增益校正示意图，具体包括以下内容：
方案二：
数据拼接方法：
拼接后的K空间数据的第1区域数据由Kτ第1区域数据填充；拼接后的K空间数据的第2区域数据由Kτ-6第2区域数据填充；拼接后的K空间数据的第3区域数据由Kτ-12第3区域数据填充。
增益校正：
在Kτ的Kτ 2和Kτ 3与Kτ-6的Kτ-6 2和Kτ-6 3中，每次顺序取出相同4行的数据，取相位编码行行号为ceil(m/4)*4的中间N个数据的模定义为δτi，δ(τ-6)i，i＝1,2,…,N。令：
Figure BDA0002966938180000081
Figure BDA0002966938180000082
βτ,τ-6即为Kτ和Kτ-6对应相位编码行的增益差因子ΔG1(ceil(m/4)*4)，在Kτ的Kτ 2和Kτ 3与Kτ-6的Kτ-6 2和Kτ-6 3中，将这些区域所有行的增益差因子ΔG1(ceil(m/4)*4)求和再取平均值得到avg[ΔG1(ceil(m/4)*4)]。avg[ΔG1(ceil(m/4)*4)]即为在Kτ的Kτ 2和Kτ 3与Kτ-6的Kτ-6 2和Kτ-6 3中的增益差因子。
在Kτ-6的Kτ-6 3与Kτ-12的Kτ-12 3中，每次顺序取出相同4行的数据，取相位编码行行号为ceil(m/4)*4的中间N个数据的模定义为δ(τ-6)i，δ(τ-12)i，i＝1,2,…,N。令：
Figure BDA0002966938180000083
Figure BDA0002966938180000084
βτ-6,τ-12即为Kτ-6 3与Kτ-12 3对应相位编码行的增益差因子ΔG2(ceil(m/4)*4)，将这些区域所有行的ΔG2(ceil(m/4)*4)求和再取平均值得到avg[ΔG2(ceil(m/4)*4)]。avg[ΔG2(ceil(m/4)*4)]即为在Kτ-6 3与Kτ-12 3中的增益差因子。
avg[ΔG1(ceil(m/4)*4)]*avg[ΔG2(ceil(m/4)*4)]则为Kτ 3与Kτ-12 3中增益差因子.
增益补偿：
假定拼接后的数据空间为K，将Kτ的Kτ 1作为K的K1。Kτ-6的Kτ-6 2作为K的K2，让K2每行每一个数据幅值与avg[ΔG1(ceil(m/4)*4)]相乘得到K2*。Kτ-12的Kτ-12 3作为K的K3，让K3每行每一个数据幅值与avg[ΔG1(ceil(m/4)*4)]*avg[ΔG2(ceil(m/4)*4)]相乘得到K3*。则K1，K2*，K3*构成的K空间为补偿后的K空间。
由于方案二的相位校正和相位补偿跟方案一相同，不再重复阐述。
图2、3中的area1表示第一区域，area2表示第二区域，area3表示第三区域。
本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的权利要求范围之内。

Claims (8)

1.一种磁共振K空间数据校正的方法，其特征在于，将三个不同增益参数下的K空间数据通过增益差校正与补偿，归一化到同一个接收增益水平；再对增益补偿后的K空间数据进行相位差校正与补偿；具体包括以下过程：
将K空间的相位编码行行号表示为m，K空间一共A个相位编码行；
分别在增益参数为：τ1、τ2、τ3下得到各自对应的K空间数据：
Figure FDA0003351458690000011
将
Figure FDA0003351458690000012
分成3个区域，
Figure FDA0003351458690000013
的第一区域数据记为
Figure FDA0003351458690000014
Figure FDA0003351458690000015
的第二区域数据记为
Figure FDA0003351458690000016
Figure FDA0003351458690000017
的第三区域数据记为
Figure FDA0003351458690000018
将
Figure FDA0003351458690000019
分成3个区域，
Figure FDA00033514586900000110
的第一区域数据记为
Figure FDA00033514586900000111
Figure FDA00033514586900000112
的第二区域数据记为
Figure FDA00033514586900000113
Figure FDA00033514586900000114
的第三区域数据记为
Figure FDA00033514586900000115
将
Figure FDA00033514586900000116
分成3个区域，
Figure FDA00033514586900000117
的第一区域数据记为
Figure FDA00033514586900000118
Figure FDA00033514586900000119
的第二区域数据记为
Figure FDA00033514586900000120
Figure FDA00033514586900000121
的第三区域数据记为
Figure FDA00033514586900000122
其中：
Figure FDA00033514586900000123
Figure FDA00033514586900000124
应的相位编码行行号为：m＝1到m＝A/6，m＝5A/6+1到m＝A；
Figure FDA00033514586900000125
对应的相位编码行行号为：m＝A/6+1到m＝A/3，m＝2A/3+1到m＝5A/6；
Figure FDA00033514586900000126
对应的相位编码行行号为：m＝A/3+1到m＝2A/3；
拼接后的数据空间为K，记K的第1区域数据为K1，K的第2区域数据为K2，K的第3区域数据为K3；将
Figure FDA00033514586900000127
的
Figure FDA00033514586900000128
作为K的K1，
Figure FDA00033514586900000129
的
Figure FDA00033514586900000130
作为K的K2，
Figure FDA00033514586900000131
的
Figure FDA00033514586900000132
作为K的K3；
通过增益差校正与补偿，将三组K空间数据归一化到同一个接收增益水平；具体的：让K2每行每一个数据幅值与第一增益差因子相乘得到K2*，让K3每行每一个数据幅值与第二增益差因子相乘得到K3，则K1、K2*、K3*构成增益补偿后的K；
对增益补偿后的K进行相位差校正和补偿。
2.根据权利要求1所述的一种磁共振K空间数据校正的方法，其特征在于，对增益补偿后的K进行相位差校正和补偿，具体的：在增益补偿后，K1数据点相位保持不变；K2*每行数据点的相位加第一相位差因子得到K2*@；K3*每行数据点的相位加第二相位差因子得到K3*@。
3.根据权利要求2所述的一种磁共振K空间数据校正的方法，其特征在于，第一增益差因子计算过程为：在
Figure FDA00033514586900000133
中，每次顺序取出相同4行的数据，取相位编码行行号为ceil(m/4)*4的中间N个数据的模定义为
Figure FDA00033514586900000134
令：
Figure FDA0003351458690000021
Figure FDA0003351458690000022
Figure FDA0003351458690000023
即为第一增益差因子。
4.根据权利要求3所述的一种磁共振K空间数据校正的方法，其特征在于，第二增益差因子计算过程为：在
Figure FDA0003351458690000024
与
Figure FDA0003351458690000025
中，每次顺序取出相同4行的数据，取相位编码行行号为ceil(m/4)*4的中间N个数据的模定义为
Figure FDA0003351458690000026
令：
Figure FDA0003351458690000027
Figure FDA0003351458690000028
Figure FDA0003351458690000029
即为
Figure FDA00033514586900000210
和
Figure FDA00033514586900000211
的对应相位编码行的增益差因子ΔG2(ceil(m/4)*4)；ΔG1(ceil(m/4)*4)*ΔG2(ceil(m/4)*4)则第二增益差因子。
5.根据权利要求2所述的一种磁共振K空间数据校正的方法，其特征在于，第一增益差因子计算过程为：在
Figure FDA00033514586900000212
中，每次顺序取出相同4行的数据，取相位编码行行号为ceil(m/4)*4的中间N个数据的模定义为
Figure FDA00033514586900000213
令：
Figure FDA00033514586900000214
Figure FDA00033514586900000215
Figure FDA00033514586900000216
为
Figure FDA00033514586900000217
和
Figure FDA00033514586900000218
对应相位编码行的增益差因子ΔG1(ceil(m/4)*4)，在
Figure FDA00033514586900000219
Figure FDA00033514586900000220
中，将这些区域所有行的增益差因子ΔG1(ceil(m/4)*4)求和再取平均值得到avg[ΔG1(ceil(m/4)*4)]，avg[ΔG1(ceil(m/4)*4)]为第一增益差因子。
6.根据权利要求5所述的一种磁共振K空间数据校正的方法，其特征在于，第二增益差因子的计算过程为：
在
Figure FDA00033514586900000221
与
Figure FDA00033514586900000222
中，每次顺序取出相同4行的数据，取相位编码行行号为ceil(m/4)*4的中间N个数据的模定义为
Figure FDA00033514586900000223
令：
Figure FDA00033514586900000224
Figure FDA0003351458690000031
Figure FDA0003351458690000032
即为
Figure FDA0003351458690000033
与
Figure FDA0003351458690000034
对应相位编码行的增益差因子ΔG2(ceil(m/4)*4)，将这些区域所有行的ΔG2(ceil(m/4)*4)求和再取平均值得到avg[ΔG2(ceil(m/4)*4)]；avg[ΔG2(ceil(m/4)*4)]即为第二增益差因子。
7.根据权利要求4或6所述的一种磁共振K空间数据校正的方法，其特征在于，第一相位差因子的计算过程为：在
Figure FDA0003351458690000035
中，每次顺序取出相同4行的数据，取相位编码行行号为ceil(m/4)*4的中间N个数据的相位定义为
Figure FDA0003351458690000036
令：
Figure FDA0003351458690000037
对γ1i进行出现频次统计，出现频次最多的值记为第一相位差因子ΔΦ1(ceil(m/4)*4)。
8.根据权利要求7所述的一种磁共振K空间数据校正的方法，其特征在于，第二相位差因子的计算过程为：在
Figure FDA0003351458690000038
中，每次顺序取出相同4行的数据，取相位编码行行号为ceil(m/4)*4的中间N个数据的相位定义为
Figure FDA0003351458690000039
令：
Figure FDA00033514586900000310
对γ2i进行出现频次统计，出现频次最多的值则为
Figure FDA00033514586900000311
和
Figure FDA00033514586900000312
对应相位编码行的相位差因子ΔΦ2(ceil(m/4)*4)；ΔΦ1(ceil(m/4)*4)+ΔΦ2(ceil(m/4)*4)记为第二相位差因子。
CN202110253687.XA 2021-03-09 2021-03-09 一种磁共振k空间数据校正的方法 Active CN113030814B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110253687.XA CN113030814B (zh) 2021-03-09 2021-03-09 一种磁共振k空间数据校正的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110253687.XA CN113030814B (zh) 2021-03-09 2021-03-09 一种磁共振k空间数据校正的方法

Publications (2)

Publication Number Publication Date
CN113030814A CN113030814A (zh) 2021-06-25
CN113030814B true CN113030814B (zh) 2021-12-21

Family

ID=76467130

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110253687.XA Active CN113030814B (zh) 2021-03-09 2021-03-09 一种磁共振k空间数据校正的方法

Country Status (1)

Country Link
CN (1) CN113030814B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5451876A (en) * 1993-10-18 1995-09-19 General Electric Company MRI system with dynamic receiver gain
CN1575747A (zh) * 2003-07-02 2005-02-09 Ge医疗系统环球技术有限公司 增益调整方法和磁共振成像设备
CN106251313A (zh) * 2016-08-15 2016-12-21 上海联影医疗科技有限公司 医学成像方法及系统
JP2017003336A (ja) * 2015-06-06 2017-01-05 国立大学法人 岡山大学 磁場計測装置及びこの磁場計測装置を用いた非破壊検査装置
CN108324276A (zh) * 2018-01-11 2018-07-27 上海东软医疗科技有限公司 磁共振成像方法和装置
CN110720047A (zh) * 2017-06-09 2020-01-21 皇家飞利浦有限公司 双回波Dixon型水/脂肪分离MR成像

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103747728B (zh) * 2012-08-04 2016-06-29 株式会社东芝 磁共振成像装置以及磁共振成像方法
CN111521963B (zh) * 2020-05-06 2021-02-19 电子科技大学 一种梯度波形速率可自动适配输出的方法及磁共振系统

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5451876A (en) * 1993-10-18 1995-09-19 General Electric Company MRI system with dynamic receiver gain
CN1575747A (zh) * 2003-07-02 2005-02-09 Ge医疗系统环球技术有限公司 增益调整方法和磁共振成像设备
JP2017003336A (ja) * 2015-06-06 2017-01-05 国立大学法人 岡山大学 磁場計測装置及びこの磁場計測装置を用いた非破壊検査装置
CN106251313A (zh) * 2016-08-15 2016-12-21 上海联影医疗科技有限公司 医学成像方法及系统
CN110720047A (zh) * 2017-06-09 2020-01-21 皇家飞利浦有限公司 双回波Dixon型水/脂肪分离MR成像
CN108324276A (zh) * 2018-01-11 2018-07-27 上海东软医疗科技有限公司 磁共振成像方法和装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
利用动态接收增益提高磁共振成像信噪比的方法;孙宏宇等;《波谱学杂志》;20141231;第31卷(第4期);515-521 *

Also Published As

Publication number Publication date
CN113030814A (zh) 2021-06-25

Similar Documents

Publication Publication Date Title
US10489943B2 (en) System and method for sparse image reconstruction
CN110490832B (zh) 一种基于正则化深度图像先验方法的磁共振图像重建方法
CN101682686B (zh) 图像处理设备、图像处理方法
CN101357063B (zh) 一种磁共振快速自旋回波成像方法
JPH0646409B2 (ja) 像を補正する装置と方法
CN110895320B (zh) 基于深度学习的无变形单次激发平面回波成像方法及装置
US10261158B2 (en) Method and apparatus for eliminating motion artifact in magnetic resonance imaging
CN108287324B (zh) 磁共振多对比度图像的重建方法和装置
CN111443318A (zh) 磁共振图像处理方法、装置、存储介质和磁共振成像系统
CN111352055B (zh) 一种前瞻性相位校正平面回波成像方法
CN113030814B (zh) 一种磁共振k空间数据校正的方法
CN109544652B (zh) 基于深度生成对抗神经网络的核磁共振多加权成像方法
CN114841872A (zh) 一种基于多智能体深度强化学习的数字半色调处理方法
CN116934755A (zh) 基于直方图均衡化的肺结核ct图像增强系统
US10852380B2 (en) Magnetic resonance image reconstruction
JPH09117427A (ja) ＭＲデータラインを得るためにＭＲパルスシーケンスを使用してＭＲデータを得て、その様なデータラインでｋ空間マトリックスを充填する方法
JPH0738746A (ja) 画像処理装置および画像形成装置
WO2008063339A1 (en) Method and system for reducing auto-moire artifacts in periodic halftone screening
CN112070676B (zh) 一种双通道多感知卷积神经网络的图片超分辨率重建方法
US20020008879A1 (en) Image processing method
CN101998030A (zh) 信号处理设备、再现设备、信号处理方法和程序
CN115685032B (zh) 一种幅条k空间下的快速自旋回波的校正方法及系统
CN112684393B (zh) 一种增强磁共振谱仪动态范围的方法
CN108333544B (zh) 平面回波成像方法和装置
JP2000354163A (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