基于迭代卡尔曼滤波的核磁共振数据噪声抑制方法技术

本发明专利技术属于核磁共振数据处理领域，具体来讲为一种基于迭代卡尔曼滤波的核磁共振数据噪声抑制方法，包括：对采集得到的核磁共振回波数据y(t)进行分段处理，时间t＝2nτ，τ为半回波间隔时间，n为回波个数，n＝1,2,3,

【技术实现步骤摘要】
基于迭代卡尔曼滤波的核磁共振数据噪声抑制方法


[0001]本专利技术属于核磁共振数据处理领域，具体来讲为一种基于迭代卡尔曼滤波的核磁共振数据噪声抑制方法。

技术介绍

[0002]核磁共振(nuclear magnetic resonance，NMR)方法基于背景场强度，可分为高场核磁共振(大于0.5T)和低场核磁共振。而低场核磁共振具有成本低、原位测试等优点，因此被广泛应用于检测物质的物理性质，主要包含食品科学、高分子材料和石油勘探等。低场核磁共振通过施加射频脉冲(CPMG序列)，激发介质中的氢质子发生共振,利用接收线圈即可采集得到NMR回波曲线。
[0003]专利ZL202010574811.8专利技术了一种核磁共振信号的去噪方法及装置，通过对核磁共振信号分离出的多个分量信号进行信噪分离，得到核磁共振信号中的噪声信号，进而实现核磁共振信号的去噪，该无需预先设置信噪分离阈值，保证了信噪分离准确性。
[0004]专利202110794109.7专利技术了一种适用于地学核磁共振仪器的核磁共振噪声抑制方法，该放大基于多线性奇异值张量分本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.基于迭代卡尔曼滤波的核磁共振数据噪声抑制方法，其特征在于，该方法包括，步骤1，对采集得到的核磁共振回波数据y(t)进行分段处理，时间t＝2nτ，τ为半回波间隔时间，n为回波个数，n＝1,2,3,
…
N，分段方案为t＝2τ～2Nτ之间等对数间隔截取K份，每段NMR回波数据记为y
k
，k＝1,2,
…
,K；步骤2，对第一段至第K段NMR回波数据，进行扩展卡尔曼滤波处理得到y
k
段回波数据对应的模型值M
k
和对应的时间t
k
；步骤3，将步骤2输出的y
k
段回波数据对应的模型值M
k
和对应的时间t
k
，k＝1,2,
…
,K代入NMR回波正演函数，得到抑制噪声后的2.按照权利要求1所述的基于迭代卡尔曼滤波的核磁共振数据噪声抑制方法，其特征在于，所述步骤2具体包括：步骤A，初始化，令k＝1，初始化M0，P0，M0＝[A0,T
20
]，A0的范围为1至回波曲线最大值M
max
间的任意实数，T
20
的范围为0.001至1之间的任意实数，模型的先验误差估计P0为M0的协方差矩阵；步骤B，构建NMR回波数据的模型空间和测量空间，分别如式(1)和式(2)所示：M
k+1
＝F(M
k
)+θ
k
ꢀꢀꢀ
(1)y
k
＝G(M
k
‑1)+ω
k
ꢀꢀꢀꢀ
(2)其中M
k
为y
k
段回波数据对应的模型值，M
k
＝[A
k
,T

【专利技术属性】
技术研发人员：矫德强，马海涛，陈楠，李岩，邵向鑫，
申请(专利权)人：长春工业大学，
类型：发明
国别省市：