一种ADC-T2二维图谱的测量方法、装置、计算机设备及非均匀场磁共振系统与流程

imaging，se-epi)，即自旋回波序列(se)进行扩散梯度编码，平面回波序列(epi)进行信号读出。在非均匀磁场磁共振系统(也简称非均匀场磁共振系统)中，类似的扩散加权成像技术被引入，用于测量物质的表观扩散系数。几种典型的测量表观扩散系数的脉冲序列如图1所示。图1中的(a)为自旋回波-cpmg序列(spin echo-cpmg sequence，se-cpmg)，即基于自旋回波进行扩散梯度编码，然后用超快速的cpmg序列进行信号读出。图1中的(b)为双自旋回波-cpmg序列(dual spin echo-cpmg sequence，dse-cpmg)，即基于双回波序列进行扩散梯度编码，同样用超快速的cpmg序列进行信号读出，该方法能降低低速液体流动带来的影响。图1中的(c)为受激回波-cpmg序列(stimulated echo-cpmg sequence，ste-cpmg)，即基于受激回波序列进行扩散梯度编码，该方法能降低t1恢复的影响，当被检测物体的t1/t2比较小时，用该序列测量表观扩散系数(apparent diffusion coefficient，adc)可以提升测量准确度。
[0006]
现有技术中，adc测量脉冲序列都是由扩散梯度编码模块和信号读出模块组成。由于在非均匀场磁共振系统中，梯度磁场非常大，通常比常规磁共振成像系统(magnetic resonance imaging，mri)的梯度场高2～3个数量级，且无法控制该梯度磁场在信号读出阶段的改变(mri系统中的dwi技术，在信号读出阶段可以控制梯度场降低，而在非均匀磁体中，由于静态梯度磁场g往往非常大，且为恒定不可控梯度场)，而该梯度磁场在信号读出阶段仍然会起到扩散编码的作用，使得会影响到表观扩散系数adc测量的准确性。因此需要采用超快速信号读出模块，才能降低读出信号过程中扩散效应的影响。即现有技术都是采用超短回波间隔的cpmg序列进行信号读出，例如使用的回波间隔为40us。这样对核磁共振系统的谱仪设备、射频功放和射频线圈等硬件要求都非常高。


技术实现要素：

[0007]
为了解决在具有极不均匀磁场的核磁共振系统中需采用超快速信号读出模块才能准确测量adc和t2的二维图谱分布，导致对谱仪设备、射频功放和射频线圈等硬件要求都非常高的问题，本发明目的在于提供一种横向磁化矢量衰减时间常数的测量方法、装置、计算机设备、非均匀场磁共振系统及计算机可读存储介质。
[0008]
第一方面，本发明提供了一种adc-t2二维图谱的测量方法，包括：
[0009]
获取基于cpmg核磁共振序列采集得到的多组回波信号，其中，所述cpmg核磁共振序列用于进行扩散梯度编码和各组回波信号的读出，并针对各组回波信号采用有不同的回波间隔；
[0010]
对所述多组回波信号进行如下积分公式的拟合：
[0011][0012]
式中，s
′
(i,j)表示在所述多组回波信号中与变量i和j对应的回波信号，i＝1,2,
…
,m，j＝1,2,
…
,n，τ
i
表示在第i组回波信号采集过程中所采用的回波间隔，m表示所述多组回波信号的总组数，n
j
表示在单组回波信号中第j个回波信号，n表示单组回波信号的回波链长度，d
max
表示表观扩散系数adc的最大值边界，d
min
表示表观扩散系数adc的最小值边界，t
2-max
表示横向磁化矢量衰减时间常数t2的最大值边界，t
2-min
表示横向磁化矢量衰减时间常数t2的最小值边界，s(d,t2)表示待求解的adc-t2二维图谱，d表示表观扩散系数
adc，t2表示横向磁化矢量衰减时间常数t2，e表示自然对数的底数；
[0013]
将所述待求解的adc-t2二维图谱的维度设定为p
×
q，然后将所述积分公式转换为如下的矩阵公式：
[0014]
s
′
＝k
·
s
[0015]
式中，s
′
＝vect(s
′
)，s
′
表示s
′
(i,j)，vect(s
′
)表示s
′
(i,j)的向量形式，k表示由m
·
n
×
p
·
q个元素k
a,b
排列成的矩阵，a＝1,2,
…
,m
·
n，b＝1,2,
…
,p
·
q，τ
a
＝k
τ
(a)＝vect(repmat(a,1,n))，n
a
＝k
n
(a)＝vect(repmat(a,m,1))，d
b
＝k
d
(b)＝vect(repmat(b,1,q))，repmat()表示用于复制和平铺矩阵的函数，vect(repmat())表示repmat()的向量形式，s＝vect(s)，s表示s(d,t2)，vect(s)表示s(d,t2)的向量形式；
[0016]
根据所述矩阵公式，将求解s(d,t2)的问题转化为如下公式的求解最优化问题：
[0017][0018]
式中，f(s)表示关于变量s的函数，|| ||2表示向量的二范数，λ表示约束项系数；
[0019]
根据所述求解最优化问题求解得到最小的且大于0的拟合向量s；
[0020]
根据求解得到的拟合向量s，还原得到所述adc-t2二维图谱的求解结果。
[0021]
基于上述发明内容，可以在具有极不均匀的磁场或者是无法实现极短的回波时间的核磁共振系统中，通过获取的多个采用不同回波间隔的cpmg序列采集的回波信号，从多组回波信号中拟合出adc系数和t2值，从而能测量出adc-t2图谱，使得不再需要复杂的扩散加权序列，具有算法简单且对系统要求低的优点，可以降低对对谱仪设备、射频功放和射频线圈等硬件系统的成本。同时所述测量方法还具有算法稳定的特点，不易受流动液体影响，对于t1/t2较小的物质同样适用。即测量序列简单快速，无需扩散梯度编码模块，相对于现有技术更容易实现；对于极不均匀的磁场，或者是无法实现极短的回波时间的核磁共振系统，该方法仍然能测量出adc-t2图谱；基于前述两点优势，在设计用于测量adc-t2的非均匀场核磁共振系统时可以降低系统硬件要求，从而降低硬件成本；所述adc-t2测量方法，不易受流动液体影响，对于t1/t2较小的物质同样适用；可以较为准确得区分被检物体中的水和脂肪成分。
[0022]
在一个可能的设计中，所述扩散敏感系数b值与(t
ee
)3成正比，其中，t
ee
表示在单组回波信号采集过程中首个90
°
激发脉冲的中心点时刻至首个回波信号的中心点时刻的时长。
[0023]
在一个可能的设计中，根据所述求解最优化问题求解得到最小的且大于0的拟合向量s，包括：
[0024]
对矩阵k进行如下的且基于奇异值分解的降维处理：
[0025][0026]
式中，u
m
·
n
×
m
·
n
表示m
·
n
×
m
·
n阶酉矩阵，表示p
·
q
×
p
·
q阶酉矩阵v
p
·
q
×
p
·
q
的共轭转置矩阵，∑
m
·
n
×
p
·
q
对角线上的元素即为奇异值，k
′
表示降维后的奇异值个数，u
m
·
n
×
k
′
表示u
m
·
n
×
m
·
n
的前k
′
列矩阵；
[0027]
按照如下公式将向量s
′
从m
×
n维度降维到k
′
维度：
[0028][0029]
式中，s
s
′
vd
表示将向量s
′
从m
×
n维度降维到k
′
维度时的向量，为矩阵u
m
·
n
×
k
′
的转置矩阵；
[0030]
将所述求解最优化问题转换为如下形式后，再求解得到所述拟合向量s：
[0031][0032]
式中，表示将矩阵k从m
×
n维度降维到k
′
维度时的矩阵。
[0033]
在一个可能的设计中，在对所述多组回波信号进行积分公式拟合前，还包括对所述多组回波信号进行如下的预处理：
[0034]
在单次读出数据的采样点数维度上，对所述多组回波信号进行快速傅立叶变换，得到频域数据，然后保留低于预设频率阈值的低频部分并进行均值化处理，得到在所述采样点数维度上无差异的多组回波信号。
[0035]
在一个可能的设计中，在对所述多组回波信号进行积分公式拟合前，还包括对所述多组回波信号进行如下的预处理：
[0036]
在平均次数维度上，对所述多组回波信号进行均值化处理，得到在所述平均次数维度上无差异的多组回波信号。
[0037]
在一个可能的设计中，在对所述多组回波信号进行积分公式拟合前，还包括对所述多组回波信号进行如下的预处理：
[0038]
在回波链长度维度上，对所述多组回波信号进行基于奇异值分解的滤波处理，得到在所述回波链长度维度上无噪声影响的多组回波信号。
[0039]
第二方面，本发明提供了一种adc-t2二维图谱的测量装置，包括依次通信连接的回波信号获取模块、积分公式拟合模块、矩阵公式转换模块、求解问题转化模块、拟合向量求解模块和求解结果还原模块；
[0040]
所述回波信号获取模块，用于获取基于cpmg核磁共振序列采集得到的多组回波信号，其中，所述cpmg核磁共振序列用于进行扩散梯度编码和各组回波信号的读出，并针对各组回波信号采用有不同的回波间隔；
[0041]
所述积分公式拟合模块，用于对所述多组回波信号进行如下积分公式的拟合：
[0042][0043]
式中，s
′
(i,j)表示在所述多组回波信号中与变量i和j对应的回波信号，i＝1,2,
…
,m，j＝1,2,
…
,n，τ
i
表示在第i组回波信号采集过程中所采用的回波间隔，m表示所述多组回波信号的总组数，n
j
表示在单组回波信号中第j个回波信号，n表示单组回波信号的回波链长度，d
max
表示表观扩散系数adc的最大值边界，d
min
表示表观扩散系数adc的最小值边界，t
2-max
表示横向磁化矢量衰减时间常数t2的最大值边界，t
2-min
表示横向磁化矢量衰减时间常数t2的最小值边界，s(d,t2)表示待求解的adc-t2二维图谱，d表示表观扩散系数adc，t2表示横向磁化矢量衰减时间常数t2，e表示自然对数的底数；
[0044]
所述矩阵公式转换模块，用于将所述待求解的adc-t2二维图谱的维度设定为p
×
q，然后将所述积分公式转换为如下的矩阵公式：
[0045]
s
′
＝k
·
s
[0046]
式中，s
′
＝vect(s
′
)，s
′
表示s
′
(i,j)，vect(s
′
)表示s
′
(i,j)的向量形式，k表示由m
·
n
×
p
·
q个元素k
a,b
排列成的矩阵，a＝1,2,
…
,m
·
n，b＝1,2,
…
,p
·
q，τ
a
＝k
τ
(a)＝vect(repmat(a,1,n))，n
a
＝k
n
(a)＝vect(repmat(a,m,1))，d
b
＝k
d
(b)＝vect(repmat(b,1,q))，repmat()表示用于复制和平铺矩阵的函数，vect(repmat())表示repmat()的向量形式，s＝vect(s)，s表示s(d,t2)，vect(s)表示s(d,t2)的向量形式；
[0047]
所述求解问题转化模块，用于根据所述矩阵公式，将求解s(d,t2)的问题转化为如下公式的求解最优化问题：
[0048][0049]
式中，f(s)表示关于变量s的函数，|| ||2表示向量的二范数，λ表示约束项系数；
[0050]
所述拟合向量求解模块，用于根据所述求解最优化问题求解得到最小的且大于0的拟合向量s；
[0051]
所述求解结果还原模块，用于根据求解得到的拟合向量s，还原得到所述adc-t2二维图谱的求解结果。
[0052]
第三方面，本发明提供了一种计算机设备，包括通信相连的存储器和处理器，其中，所述存储器用于存储计算机程序，所述处理器用于读取所述计算机程序，执行如第一方面或第一方面中任意一种可能设计所述的测量方法。
[0053]
第四方面，本发明提供了一种非均匀场磁共振系统，包括有控制台、核磁共振谱仪、磁体和射频子系统；
[0054]
所述控制台，通信连接所述核磁共振谱仪，用于向所述核磁共振谱仪发送指令，以便控制测量序列的参数选择和感兴趣区域定位，并接收由所述核磁共振谱仪采集到的磁共振信号，完成数据处理，其中，所述数据处理包括执行如第一方面或第一方面中任意一种可能设计所述的测量方法；
[0055]
所述核磁共振谱仪，通信连接所述射频子系统，用于执行来自所述控制台的指令，并通过所述射频子系统发射测量序列的射频激励信号和接收所述磁共振信号；
[0056]
所述磁体，用于布置在被检物体和所述射频子系统中信号收发部件的正上方；
[0057]
所述射频子系统，用于在所述核磁共振谱仪的控制下，发射测量序列的射频激励信号和接收所述磁共振信号。
[0058]
在一个可能的设计中，所述射频子系统包括有射频功率放大器、前置放大器、收发转换开关和射频线圈；
[0059]
所述射频功率放大器的信号输入端电连接所述核磁共振谱仪的信号输出端，所述射频功率放大器的信号输出端电连接所述收发转换开关的第一切换端；
[0060]
所述前置放大器的信号输入端电连接所述收发转换开关的第二切换端，所述前置放大器的信号输出端电连接所述核磁共振谱仪的信号输入端；
[0061]
所述收发转换开关的受控端通信连接所述核磁共振谱仪的输出端，所述收发转换
开关的切换公共端电连接所述射频线圈；
[0062]
所述射频线圈，用于作为所述射频子系统的信号收发部件，向被检物体发射所述射频激励信号和接收来自被检物体的磁共振信号。
[0063]
第五方面，本发明提供了一种计算机可读存储介质，所述计算机可读存储介质上存储有指令，当所述指令在计算机上运行时，执行如上第一方面或第一方面中任意一种可能设计的所述测量方法。
[0064]
第六方面，本发明提供了一种包含指令的计算机程序产品，当所述指令在计算机上运行时，使所述计算机执行如上第一方面或第一方面中任意一种可能设计的所述测量方法。
附图说明
[0065]
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0066]
图1是现有技术在非均匀场磁共振系统中所采用的adc测量脉冲序列的示例图，其中，(a)采用se-cpmg序列；(b)采用dse-cpmg序列；(c)采用ste-cpmg序列。
[0067]
图2是本发明提供的测量方法的流程示意图。
[0068]
图3是本发明提供的在非均匀场磁共振系统中所采用的adc-t2测量脉冲序列的示例图。
[0069]
图4是本发明提供的当被检物质为0.5mmol/l mncl2溶液和花生油时针对不同回波间隔τ＝te所获取的cpmg测量数据示例图。
[0070]
图5是本发明提供的当被检物质为0.5mmol/l mncl2溶液和花生油时测量得到的adc-t2二维图谱示例图。
[0071]
图6是本发明提供的测量装置的结构示意图。
[0072]
图7是本发明提供的计算机设备的结构示意图。
[0073]
图8是本发明提供的非均匀场磁共振系统的结构示意图。
具体实施方式
[0074]
下面结合附图及具体实施例来对本发明作进一步阐述。在此需要说明的是，对于这些实施例方式的说明虽然是用于帮助理解本发明，但并不构成对本发明的限定。本文公开的特定结构和功能细节仅用于描述本发明的示例实施例。然而，可用很多备选的形式来体现本发明，并且不应当理解为本发明限制在本文阐述的实施例中。
[0075]
应当理解，尽管本文可能使用术语第一、第二等等来描述各种单元，但是这些单元不应当受到这些术语的限制。这些术语仅用于区分一个单元和另一个单元。例如可以将第一单元称作第二单元,并且类似地可以将第二单元称作第一单元，同时不脱离本发明的示例实施例的范围。
[0076]
应当理解，对于本文中可能出现的术语“和/或”，其仅仅是一种描述关联对象的关联关系，表示可以存在三种关系，例如，a和/或b，可以表示：单独存在a，单独存在b，同时存
在a和b三种情况；对于本文中可能出现的术语“/和”，其是描述另一种关联对象关系，表示可以存在两种关系，例如，a/和b，可以表示：单独存在a，单独存在a和b两种情况；另外，对于本文中可能出现的字符“/”，一般表示前后关联对象是一种“或”关系。
[0077]
应当理解，在本文中若将单元称作与另一个单元“连接”、“相连”或“耦合”时，它可以与另一个单元直相连接或耦合，或中间单元可以存在。相対地，在本文中若将单元称作与另一个单元“直接相连”或“直接耦合”时，表示不存在中间单元。另外，应当以类似方式来解释用于描述单元之间的关系的其他单词(例如，“在
……
之间”对“直接在
……
之间”,“相邻”对“直接相邻”等等)。
[0078]
应当理解，本文使用的术语仅用于描述特定实施例，并不意在限制本发明的示例实施例。若本文所使用的，单数形式“一”、“一个”以及“该”意在包括复数形式，除非上下文明确指示相反意思。还应当理解，若术语“包括”、“包括了”、“包含”和/或“包含了”在本文中被使用时,指定所声明的特征、整数、步骤、操作、单元和/或组件的存在性,并且不排除一个或多个其他特征、数量、步骤、操作、单元、组件和/或他们的组合存在性或增加。
[0079]
应当理解，还应当注意到在一些备选可能设计中，所出现的功能/动作可能与附图出现的顺序不同。例如,取决于所涉及的功能/动作,实际上可以实质上并发地执行,或者有时可以以相反的顺序来执行连续示出的两个图。
[0080]
应当理解，在下面的描述中提供了特定的细节,以便于对示例实施例的完全理解。然而,本领域普通技术人员应当理解可以在没有这些特定细节的情况下实现示例实施例。例如可以在框图中示出系统，以避免用不必要的细节来使得示例不清楚。在其他实例中，可以不以非必要的细节来示出众所周知的过程、结构和技术，以避免使得示例实施例不清楚。
[0081]
如图2～5所示，本实施例第一方面提供的所述adc-t2二维图谱的测量方法,可以但不限于由在具有极不均匀的磁场或者是无法实现极短回波时间的核磁共振系统中的控制台(其与核磁共振谱仪连接，以便控制测量序列的参数选择和感兴趣区域定位，并接收由所述核磁共振谱仪采集到的磁共振信号，完成数据处理)执行，以便得到准确的adc-t2二维图谱。所述adc-t2二维图谱的测量方法，包括但不限于有如下步骤s101～s106。
[0082]
s101.获取基于cpmg核磁共振序列采集得到的多组回波信号，其中，所述cpmg核磁共振序列用于进行扩散梯度编码和各组回波信号的读出，并针对各组回波信号采用有不同的回波间隔。
[0083]
在所述步骤s101中，如图3所示，即采用一个典型的θ-2θ-2θ-2θ
……
射频脉冲序列：首个激发脉冲翻转角为θ，其后跟随若干个回聚脉冲，翻转角为2θ；首个激发脉冲和首个回聚脉冲之间的相位差为90度，首个激发脉冲与首个回聚脉冲之间的时间间隔为τ/2,首个回聚脉冲到首个采样窗之间的时间间隔为τ/2；回聚脉冲之间的时间间隔均为τ，称为回波间隔。在核磁共振系统中，恒定梯度场为磁体的天然梯度场，不需要控制。此外，虽然一次激发即可采集n个回波信号(其中，n为正整数)，但也需要重复多次，以便通过平均信号提升信噪比，以及为了后续能够拟合得到adc-t2二维图谱，需要通过改变回波间隔τ来采集多组回波信号。
[0084]
s102.对所述多组回波信号进行如下积分公式的拟合：
[0085]
[0086]
式中，s
′
(i,j)表示在所述多组回波信号中与变量i和j对应的回波信号，i＝1,2,
…
,m，j＝1,2,
…
,n，τ
i
表示在第i组回波信号采集过程中所采用的回波间隔，m表示所述多组回波信号的总组数，n
j
表示在单组回波信号中第j个回波信号，n表示单组回波信号的回波链长度，d
max
表示表观扩散系数adc的最大值边界，d
min
表示表观扩散系数adc的最小值边界，t
2-max
表示横向磁化矢量衰减时间常数t2的最大值边界，t
2-min
表示横向磁化矢量衰减时间常数t2的最小值边界，s(d,t2)表示待求解的adc-t2二维图谱，d表示表观扩散系数adc，t2表示横向磁化矢量衰减时间常数t2，e表示自然对数的底数。
[0087]
在所述步骤s102中，采集获取的所述多组回波信号可用一个4维数组s(m,n,a,q)来表示，其中，第一维数组s(m)基于不同的回波间隔，长度为正整数m，即对应m个不同长度的回波间隔τ；第二维数组s(n)基于不同的回波链长度，长度为正整数n；第三维数组s(a)基于在进行重复采集时的不同平均次数，长度为正整数a；第四维数组s(q)基于单次读出数据的不同采样点数，长度为正整数q。
[0088]
在所述步骤s102中，考虑所述多组回波信号是一个4维数组s(m,n,a,q)，维度过多且存在噪声影响，会增加不必要的积分计算量以及带来误差，由此有必要先进行降维和/或去噪处理，即在对所述多组回波信号进行积分公式拟合前，还包括但不限于对所述多组回波信号进行如下方式(a)至(c)中任意一项或任意组合的预处理：(a)在单次读出数据的采样点数维度上，对所述多组回波信号进行快速傅立叶变换，得到频域数据，然后保留低于预设频率阈值的低频部分并进行均值化处理，得到在所述采样点数维度上无差异的多组回波信号；(b)在平均次数维度上，对所述多组回波信号进行均值化处理，得到在所述平均次数维度上无差异的多组回波信号；(c)在回波链长度维度上，对所述多组回波信号进行基于奇异值分解(singular value decompositionm,简称svd，是在机器学习领域应用较为广泛的算法之一，也是学习机器学习算法绕不开的基石之一。svd算法主要用在降维算法中的特征分解、推荐系统、自然语言处理计算机视觉等领域，通俗一点讲就是将一个线性变换分解为两个线性变换，一个线性变换代表旋转，一个线性变换代表拉伸)的滤波处理，得到在所述回波链长度维度上无噪声影响的多组回波信号。由此通过前述方式(a)和(b)的处理，可将四维数组s(m,n,a,q)降维成二维数组s(m,n)，利于减少后续计算。
[0089]
s103.将所述待求解的adc-t2二维图谱的维度设定为p
×
q，然后将所述积分公式转换为如下的矩阵公式：
[0090]
s
′
＝k
·
s
[0091]
式中，s
′
＝vect(s
′
)，s
′
表示s
′
(i,j)，vect(s
′
)表示s
′
(i,j)的向量形式，k表示由m
·
n
×
p
·
q个元素k
a,b
排列成的矩阵，a＝1,2,
…
,m
·
n，b＝1,2,
…
,p
·
q，τ
a
＝k
τ
(a)＝vect(repmat(a,1,n))，n
a
＝k
n
(a)＝vect(repmat(a,m,1))，d
b
＝k
d
(b)＝vect(repmat(b,1,q))，repmat()表示用于复制和平铺矩阵的函数，vect(repmat())表示repmat()的向量形式，s＝vect(s)，s表示s(d,t2)，vect(s)表示s(d,t2)的向量形式。
[0092]
s104.根据所述矩阵公式，将求解s(d,t2)的问题转化为如下公式的求解最优化问题：
[0093][0094]
式中，f(s)表示关于变量s的函数，|| ||2表示向量的二范数，λ表示约束项系数。
[0095]
s105.根据所述求解最优化问题求解得到最小的且大于0的拟合向量s。
[0096]
在所述步骤s105中，根据所述求解最优化问题求解所述拟合向量s的具体过程为现有常规求解方式。优化的，为了降低在求解过程中的计算量，可以但不限于按照如下步骤s501～s503进行降维求解。
[0097]
s501.对矩阵k进行如下的且基于奇异值分解的降维处理：
[0098][0099]
式中，u
m
·
n
×
m
·
n
表示m
·
n
×
m
·
n阶酉矩阵，表示p
·
q
×
p
·
q阶酉矩阵v
p
·
q
×
p
·
q
的共轭转置矩阵，∑
m
·
n
×
p
·
q
对角线上的元素即为奇异值，k
′
表示降维后的奇异值个数，u
m
·
n
×
k
′
表示u
m
·
n
×
m
·
n
的前k
′
列矩阵；
[0100]
s502.按照如下公式将向量s
′
从m
×
n维度降维到k
′
维度：
[0101][0102]
式中，s
s
′
vd
表示将向量s
′
从m
×
n维度降维到k
′
维度时的向量，为矩阵u
m
·
n
×
k
′
的转置矩阵。
[0103]
s503.将所述求解最优化问题转换为如下形式后，再求解得到所述拟合向量s：
[0104][0105]
式中，表示将矩阵k从m
×
n维度降维到k
′
维度时的矩阵。
[0106]
s106.根据求解得到的拟合向量s，还原得到所述adc-t2二维图谱的求解结果。
[0107]
在所述步骤s106中，由于拟合向量s＝vect(s)＝vect(s(d,t2))，因此可通过常规方式还原得到所述adc-t2二维图谱的求解结果。即考虑非均匀场中存在的强梯度场，在cpmg序列中既可以起频率编码的作用，同时也始终存在扩散编码的作用。也就是说如果cpmg序列的回波间隔越大，则信号越低，其中，扩散效应对cpmg信号的影响，可以用如下公式描述：
[0108][0109]
式中，γ表示原子核的磁旋比，d
′
表示物质的adc系数，g表示梯度磁场大小，τ表示cpmg序列的回波间隔。由该公式可以看出，cpmg信号和物质的adc系数、梯度场和回波间隔有关，由此可以通过采用不同回波间隔的cpmg信号估计出adc系数和t2值，得到所述adc-t2二维图谱。
[0110]
如图4～5所示，下面以在非均匀场核磁共振系统上对0.5mmol/l mncl2溶液和花生油进行检测为例，来具体说明本实施例的实验结果。实验参数主要如下：感兴趣区域(region of interest，roi)的恒定磁场为0.07t，梯度场为110gauss/cm，cpmg序列扩散梯度编码的回波间隔τ分别为0.5ms、1.2ms、1.5ms、1.7ms、1.85ms、2.0ms、2.5ms、6.0ms、7.5ms、8.5ms、9.0ms和10.0ms共12个(即m＝12)，cpmg序列的回波链长度为256(即n＝256)，平均次数为64次(即a＝64)，回波采样点数为64(即q＝64)。
[0111]
根据图4可见具有不同回波间隔的cpmg序列采集信号衰减程度不一样，这是由于分子扩散的影响；根据图5可见0.5mmol/l mncl2溶液在adc-t2二维图谱上为一均匀峰，而花生油在adc-t2二维图谱为双峰，两种物质的adc系数和t2值均与理论值接近。此外，还可以看到，mncl2溶液和花生油在图谱上的峰区分明显：mncl2溶液的adc较大，而花生油具有更小的adc。
[0112]
由此通过上述步骤s101～s106所详细描述的测量方案，可以在具有极不均匀的磁场或者是无法实现极短的回波时间的核磁共振系统中，通过获取的多个采用不同回波间隔的cpmg序列采集的回波信号，从多组回波信号中拟合出adc系数和t2值，从而能测量出adc-t2图谱，使得不再需要复杂的扩散加权序列，具有算法简单且对系统要求低的优点，可以降低对对谱仪设备、射频功放和射频线圈等硬件系统的成本。同时所述测量方法还具有算法稳定的特点，不易受流动液体影响，对于t1/t2较小的物质同样适用。即测量序列简单快速，无需扩散梯度编码模块，相对于现有技术更容易实现；对于极不均匀的磁场，或者是无法实现极短的回波时间的核磁共振系统，该方法仍然能测量出adc-t2图谱；基于前述两点优势，在设计用于测量adc-t2的非均匀场核磁共振系统时可以降低系统硬件要求，从而降低硬件成本；所述adc-t2测量方法，不易受流动液体影响，对于t1/t2较小的物质同样适用；可以较为准确得区分被检物体中的水和脂肪成分。
[0113]
如图6所示，本实施例第二方面提供了一种实现第一方面或第一方面中任意一种可能设计的所述测量方法的虚拟装置，包括依次通信连接的回波信号获取模块、积分公式拟合模块、矩阵公式转换模块、求解问题转化模块、拟合向量求解模块和求解结果还原模块；
[0114]
所述回波信号获取模块，用于获取基于cpmg核磁共振序列采集得到的多组回波信号，其中，所述cpmg核磁共振序列用于进行扩散梯度编码和各组回波信号的读出，并针对各组回波信号采用有不同的回波间隔；
[0115]
所述积分公式拟合模块，用于对所述多组回波信号进行如下积分公式的拟合：
[0116][0117]
式中，s
′
(i,j)表示在所述多组回波信号中与变量i和j对应的回波信号，i＝1,2,
…
,m，j＝1,2,
…
,n，τ
i
表示在第i组回波信号采集过程中所采用的回波间隔，m表示所述多组回波信号的总组数，n
j
表示在单组回波信号中第j个回波信号，n表示单组回波信号的回波链长度，d
max
表示表观扩散系数adc的最大值边界，d
min
表示表观扩散系数adc的最小值边界，t
2-max
表示横向磁化矢量衰减时间常数t2的最大值边界，t
2-min
表示横向磁化矢量衰减时间常数t2的最小值边界，s(d,t2)表示待求解的adc-t2二维图谱，d表示表观扩散系数adc，t2表示横向磁化矢量衰减时间常数t2，e表示自然对数的底数；
[0118]
所述矩阵公式转换模块，用于将所述待求解的adc-t2二维图谱的维度设定为p
×
q，然后将所述积分公式转换为如下的矩阵公式：
[0119]
s
′
＝k
·
s
[0120]
式中，s
′
＝vect(s
′
)，s
′
表示s
′
(i,j)，vect(s
′
)表示s
′
(i,j)的向量形式，k表示由m
·
n
×
p
·
q个元素k
a,b
排列成的矩阵，a＝1,2,
…
,m
·
n，b＝1,
2,
…
,p
·
q，τ
a
＝k
τ
(a)＝vect(repmat(a,1,n))，n
a
＝k
n
(a)＝vect(repmat(a,m,1))，d
b
＝k
d
(b)＝vect(repmat(b,1,q))，repmat()表示用于复制和平铺矩阵的函数，vect(repmat())表示repmat()的向量形式，s＝vect(s)，s表示s(d,t2)，vect(s)表示s(d,t2)的向量形式；
[0121]
所述求解问题转化模块，用于根据所述矩阵公式，将求解s(d,t2)的问题转化为如下公式的求解最优化问题：
[0122][0123]
式中，f(s)表示关于变量s的函数，|| ||2表示向量的二范数，λ表示约束项系数；
[0124]
所述拟合向量求解模块，用于根据所述求解最优化问题求解得到最小的且大于0的拟合向量s；
[0125]
所述求解结果还原模块，用于根据求解得到的拟合向量s，还原得到所述adc-t2二维图谱的求解结果。
[0126]
本实施例第二方面提供的前述装置的工作过程、工作细节和技术效果，可以参见第一方面或第一方面中任意一种可能设计所述的测量方法，于此不再赘述。
[0127]
如图7所示，本实施例第三方面提供了一种执行第一方面或第一方面中任意一种可能设计的所述测量方法的计算机设备，包括通信相连的存储器和处理器，其中，所述存储器用于存储计算机程序，所述处理器用于读取所述计算机程序，执行如第一方面或第一方面中任意一种可能设计所述的测量方法。具体举例的，所述存储器可以但不限于包括随机存取存储器(random-access memory，ram)、只读存储器(read-only memory，rom)、闪存(flash memory)、先进先出存储器(first input first output，fifo)和/或先进后出存储器(first input last output，filo)等等；所述处理器可以不限于采用型号采用stm32f105系列的微处理器。此外，所述计算机设备还可以但不限于包括有电源模块、显示屏和其它必要的部件。
[0128]
本实施例第三方面提供的前述计算机设备的工作过程、工作细节和技术效果，可以参见第一方面或第一方面中任意一种可能设计所述的测量方法，于此不再赘述。
[0129]
如图8所示，本实施例第四方面提供了一种执行第一方面或第一方面中任意一种可能设计的所述测量方法的非均匀场磁共振系统，包括有控制台、核磁共振谱仪、磁体和射频子系统；
[0130]
所述控制台，通信连接所述核磁共振谱仪，用于向所述核磁共振谱仪发送指令，以便控制测量序列的参数选择和感兴趣区域定位，并接收由所述核磁共振谱仪采集到的磁共振信号，完成数据处理，其中，所述数据处理包括执行如第一方面或第一方面中任意一种可能设计所述的测量方法；
[0131]
所述核磁共振谱仪，通信连接所述射频子系统，用于执行来自所述控制台的指令，并通过所述射频子系统发射测量序列的射频激励信号和接收所述磁共振信号；
[0132]
所述磁体，用于布置在被检物体和所述射频子系统中信号收发部件的正上方；
[0133]
所述射频子系统，用于在所述核磁共振谱仪的控制下，发射测量序列的射频激励信号和接收所述磁共振信号。
[0134]
如图8所示，在所述非均匀场磁共振系统的具体结构中，所述磁体一般设计为永磁
体，例如单边永磁体，使得在感兴趣区域内仍然具有高度不均匀的磁场。
[0135]
在一种可能设计中，所述射频子系统包括有射频功率放大器、前置放大器、收发转换开关和射频线圈；
[0136]
所述射频功率放大器的信号输入端电连接所述核磁共振谱仪的信号输出端，所述射频功率放大器的信号输出端电连接所述收发转换开关的第一切换端；
[0137]
所述前置放大器的信号输入端电连接所述收发转换开关的第二切换端，所述前置放大器的信号输出端电连接所述核磁共振谱仪的信号输入端；
[0138]
所述收发转换开关的受控端通信连接所述核磁共振谱仪的输出端，所述收发转换开关的切换公共端电连接所述射频线圈；
[0139]
所述射频线圈，用于作为所述射频子系统的信号收发部件，向被检物体发射所述射频激励信号和接收来自被检物体的磁共振信号。
[0140]
如图8所示，在所述射频子系统的具体结构中，所述射频功率放大器用于放大待发射的所述射频激励信号；所述前置放大器用于放大接收的所述磁共振信号。所述收发转换开关用于通过切换控制，使所述射频线圈即可发射所述射频激励信号，也能接收所述磁共振信号。
[0141]
本实施例第四方面提供的前述非均匀场磁共振系统的工作过程、工作细节和技术效果，可以参见第一方面或第一方面中任意一种可能设计所述的测量方法，于此不再赘述。
[0142]
本实施例第五方面提供了一种存储包含第一方面或第一方面中任意一种可能设计的所述测量方法的指令的计算机可读存储介质，即所述计算机可读存储介质上存储有指令，当所述指令在计算机上运行时，执行如第一方面或第一方面中任意一种可能设计所述的测量方法。其中，所述计算机可读存储介质是指存储数据的载体，可以但不限于包括软盘、光盘、硬盘、闪存、优盘和/或记忆棒(memory stick)等，所述计算机可以是通用计算机、专用计算机、计算机网络、或者其他可编程装置。
[0143]
本实施例第五方面提供的前述计算机可读存储介质的工作过程、工作细节和技术效果，可以参见第一方面或第一方面中任意一种可能设计所述的测量方法，于此不再赘述。
[0144]
本实施例第六方面提供了一种包含指令的计算机程序产品，当所述指令在计算机上运行时，使所述计算机执行如第一方面或第一方面中任意一种可能设计所述的测量方法。其中，所述计算机可以是通用计算机、专用计算机、计算机网络、或者其他可编程装置。
[0145]
以上所描述的实施例仅仅是示意性的，若涉及到作为分离部件说明的单元，其可以是或者也可以不是物理上分开的；若涉及到作为单元显示的部件，其可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下，即可以理解并实施。
[0146]
以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换。而这些修改或者替换，并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围。
[0147]
最后应说明的是，本发明不局限于上述可选的实施方式，任何人在本发明的启示下都可得出其他各种形式的产品。上述具体实施方式不应理解成对本发明的保护范围的限
制，本发明的保护范围应当以权利要求书中界定的为准，并且说明书可以用于解释权利要求书。