基于外插磁化转移信号的CEST数据拟合方法、装置及介质

文档序号:26050680发布日期:2021-07-27 15:25阅读:来源:国知局

技术特征:

1.一种基于外插磁化转移信号的cest数据拟合方法,其特征在于,针对cest图像中的每一个体素,依次执行s1~s3,得到每一个体素的消除背景干扰信号的cest效应值;其中s1~s3为:

s1、对获得的原始z谱进行主磁场频率b0偏移校正,得到校正后z谱;

s2、取校正后z谱中仅含有mt效应的部分数据点作为拟合样本数据,对预先设置有待定参数初值和上下限的两池模型进行拟合,得到待定参数的拟合值;

所述两池模型表示为

式中:exp表示以自然常数e为底的指数函数;δt表示时刻t的增量;分别为t+δt时刻、t时刻的磁化矢量其中矩阵m=[mxa,mya,mza,mzb]t,mxa、mya、mza、mzb分别为水集团的x方向磁化强度分量、水集团的y方向磁化强度分量、水集团的z方向磁化强度分量、mt集团的z方向磁化强度分量;系数矩阵其中矩阵向量c=[0,0,r1am0a,r1bm0b]t,δωa表示外界施加的射频脉冲的频率和水的频率差,ω1表示所施加的射频脉冲强度,kab表示从水集团交换到mt集团的速率,kba表示从mt集团交换到水集团的速率,rrfb为mt集团对射频脉冲的吸收速率,δωb表示外界施加的射频脉冲的频率和mt的频率差,r1a=1/t1a和r2a=1/t2a分别为水集团的纵向弛豫速率和横向弛豫速率,r1b=1/t1b和r2b=1/t2b分别为mt集团的纵向弛豫速率和横向弛豫速率,m0a和m0b分别为水集团和mt集团的稳态磁化强度;

所述两池模型中的待定参数为t1a,t2a,t2b,m0b,kba;

s3、基于代入所述拟合值的两池模型得到包含mt效应的拟合曲线,将所述拟合曲线与所述校正后z谱上感兴趣频率处的差值,作为消去背景干扰信号的cest效应值。

2.如权利要求1所述的基于外插磁化转移信号的cest数据拟合方法,其特征在于,所述s2中,两池模型中待定参数的拟合具有两轮,第一轮拟合时对5个待定参数在其预先设定的上下限范围内均利用所述拟合样本数据进行参数优化,获得第一轮拟合值;第二轮拟合时,以第一轮拟合值为中心,选择相对于第一轮拟合上下限范围更窄的上下限,重新使用和第一轮拟合中相同的拟合样本数据对所述两池模型进行第二轮精度更高的拟合,得到待定参数的最终拟合值;

优选的,在进行第二轮拟合时,在所述两池模型的5个待定参数中仅选择部分参数调整其上下限,剩余参数的上下限与第一轮拟合保持一致;

进一步优选的,在进行第二轮拟合时,需改变上下限的参数组合为t1a,t2b,m0b或t1a,t2b或t2b,m0b。

3.如权利要求1所述的基于外插磁化转移信号的cest数据拟合方法,其特征在于,所述s2中,仅含有mt效应的数据点为频率位于80~6ppm范围内的数据点,优选为频率位于80~20ppm范围内的数据点。

4.如权利要求2所述的基于外插磁化转移信号的cest数据拟合方法,其特征在于,所述s2中,进行第一轮拟合时,待定参数t1a,t2a,t2b,m0b,kba预设的上限是[1.5,0.25,12e-6,0.16,40],预设的下限是[0.6,0.065,8e-6,0.001,20]。

5.如权利要求2所述的基于外插磁化转移信号的cest数据拟合方法,其特征在于,所述s2中,在进行第二轮拟合时,对于需改变上下限的参数,其上限值改为该参数的第一轮拟合值的120%~150%,下限值改为该参数的第一轮拟合值的80%~50%。

6.如权利要求2所述的基于外插磁化转移信号的cest数据拟合方法,其特征在于,所述s2中,mt集团对射频脉冲的吸收速率其中g(·)是一个线形函数,包括洛伦兹型、高斯型或超洛伦兹型。

7.如权利要求6所述的基于外插磁化转移信号的cest数据拟合方法,其特征在于,所述s3中,所述感兴趣频率为3.5ppm或-3.5ppm,3.5ppm的感兴趣频率对应的cest效应值为apt#值,-3.5ppm的感兴趣频率对应的cest效应值为noe#

8.一种计算机可读存储介质,其特征在于,所述存储介质上存储有计算机程序,当所述计算机程序被处理器执行时,实现如权利要求1~7任一项所述的基于外插磁化转移信号的cest数据拟合方法。

9.一种基于外插磁化转移信号数值拟合的cest数据拟合装置,其特征在于,包括处理器和存储介质,所述存储介质上存储有计算机程序,当所述计算机程序被处理器执行时,实现如权利要求1~7任一项所述的基于外插磁化转移信号的cest数据拟合方法。

10.一种去除cest成像干扰效应的磁共振成像装置,其特征在于,包括磁共振扫描器以及控制单元,所述磁共振扫描器用于通过磁共振cest成像获取cest图像;所述控制单元能获取所述cest图像且控制单元中存储有计算机程序,当所述计算机程序被执行时,用于实现如权利要求1~7任一项所述的基于外插磁化转移信号的cest数据拟合方法,输出每一个体素消除背景干扰信号的cest效应值。


技术总结
本发明公开了一种基于外插磁化转移信号的CEST数据拟合方法、装置及介质,属于磁共振成像领域。本发明提供的基于外插磁化转移信号的CEST数据拟合方法,将MT效应的线型融入到Bloch‑McConnell方程中,通过对该方程代表的模型的拟合,可以得出MT效应,将拟合得出的曲线与实验采集到的曲线作差,就可以去除干扰效应。本发明的基于Bloch‑McConnell的拟合方法,可以考虑实验中实际使用的射频脉冲波形,并且在系统非稳态下也可以使用。

技术研发人员:张祎;雍兴旺;吴丹
受保护的技术使用者:浙江大学
技术研发日:2021.04.29
技术公布日:2021.07.27
当前第2页1 2 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1