[发明专利]基于多线性奇异值张量分解核磁共振FID信号噪声抑制方法

说明书

本发明提供了一种基于多线性奇异值张量分解核磁共振FID信号噪声抑制方法，包括：获取多通道信号：通过延时采样，每间隔时间Δt，使能一个通道开始采样，得到多通道信号；将每个通道信号变换为Hankel矩阵，构成三阶张量对三阶张量进行Tucker分解并处理，得到三阶张量将三阶张量恢复为多通道信号，对该多通道信号构成的二阶张量X进行CP张量分解处理，最后得到多通信信号融合后的高信噪比信号x_new。本发明有益效果是：有效克服了当前算法在强噪声干扰下的局限性，同时也保证了算法的实时性与普适性，可适用于质子磁力仪、核磁共振找水仪等仪器中，提高其频域测量的精度、水文参数的准确度，有效提高相关地学仪器的精度与准确度。

技术领域

本发明涉及信号处理领域，尤其涉及基于多线性奇异值张量分解核磁共振FID信号噪声抑制方法及方法。

背景技术

自由感应衰减(Free induction decay,FID)信号是核磁共振技术中产生的信号，在利用核磁共振技术的地体物理探测领域的相关仪器中(如质子旋进类磁力仪、核磁共振找水仪等)都是通过对核磁共振FID信号的参数进行提取、频率测量、水文参数反演等手段来获得所需的地磁场值、水文参数等，但是获得的FID信号都极其微弱，FID信号很容易淹没在环境噪声和人为噪声中。所以如何让抑制各种类别的噪声成为核磁共振技术在相关地学仪器中的主要限制。因此，为了提高地学仪器的测量精度和准确度，需要对获取到的FID信号进行噪声抑制，以得到更准确的频率、水文参数等信息。

当前阶段的FID信号噪声抑制手段主要是通过算法来进行噪声抑制，常用的算法有SVD、小波变换、PCA、EMD等算法，但是这些算法在实际应用中仍存在以下局限性：1)数据量要求较高，在采样点数很少的情况下，其噪声抑制效果会大打折扣；2)在时频域进行分析处理，其运算量较大，实时性无法保障；3)普适性较弱，这些算法针对特定类别的噪声进行噪声抑制，但在实际情况中会有多种噪声同时存在。因此当获得的FID信号噪声较大且噪声类别较多时，这些算法已近无法进行有效的噪声抑制，使得仪器无法正常工作。

发明内容

有鉴于此，针对以上技术缺陷，本发明提出了一种基于多线性奇异值张量分解核磁共振FID信号噪声抑制方法，包括以下步骤：

S101：获取多通道信号：通过延时采样，每间隔时间Δt，使能一个通道开始采样，得到多通道信号；

S102：基于多线性奇异值张量分解的噪声抑制：

将每个通道信号变换为Hankel矩阵，构成三阶张量

对三阶张量进行Tucker分解并处理，得到三阶张量

将三阶张量恢复为多通道信号，对该多通道信号构成的二阶张量X进行CP张量分解处理，最后得到多通信信号融合后的高信噪比信号x_new。

进一步地，步骤S102中获得的三阶张量进行Tucker分解处理的具体过程为：

对张量进行Tucker分解，得到一个核心张量与多个伴随矩阵U⁽ⁿ⁾；

对核心张量中数值相对较小的奇异值进行置零处理，得到处理后的核心张量将处理后的核心张量与上一步得到的多个伴随矩阵U⁽ⁿ⁾进行张量模乘运算，得到三阶张量即完成一部分的噪声抑制。

进一步地，对三阶张量恢复为多通道信号X，并进行CP张量分解的具体过程为：

将三阶张量通过Hankel逆变换为多通道信号矩阵X(即二阶张量)；利用CP张量分解将二阶张量X分解为两组一维向量外积的形式，得到a，b向量，对a，b向量进行运算，得到融合后的去噪FID信号，完成全部噪声抑制。

进一步地，对a，b向量进行运算的过程如下：