一种基于通道加权鲁棒性主成分分析的三维脉波影像去噪方
法及***
技术领域
本发明属于电子医疗仪器领域,特别涉及一种基于通道加权鲁棒性主 成分分析的三维脉波影像去噪方法及***。
背景技术
“望,闻,问,切”是中医诊断的四大技术。其中“切诊”,即脉诊,是 其中最具中医特色而又最难掌握的。中医诊脉,要求诊者仅凭指端触觉, 感知出脉压搏动的形态细节及其变化趋势,进而准确辨识出几十种脉象, 为中医的辩证施治提供正确的症状信息。
中医诊脉仪正是效仿中医诊脉方法开发的电子医用仪器。使用时,按 照中医取脉的位置设置传感探头,将各脉位处脉搏跳动引起的压力变化转 换成电信号,以曲线的形式显示或记录下来。然后参照西医心电图的方法 对脉搏曲线做波形分析,以寻找和确定脉搏波形与中医脉象之间的对应关 系。但是脉象信号预处理的效果将直接影响后续临床分析的结果。电源干 扰、在传感器上细微的手腕运动、呼吸引起的基线漂移,以及可能由于混沌现象引起的周期偏差,都会成为脉象信号的噪声。如果这些噪声没有被 去除,将干扰脉象信号的后续分析。
相对于采集扰动脉某个点的压力脉波形信号,采集某个区域内多个点 所产生的三维脉象,测量到的噪声干扰较中央的点明显得多。这是由于干 扰由原来的线转变成了动态的面。特别地,在动态面上分布较为边缘的 点。当使用机器手指进行切脉时,传感器除了能接收到桡动脉的脉搏信号 外,还会受到其他皮下组织的干扰。去噪的目标就将脉搏产生的动态压保 留下来,将干扰去除;对边缘点而言,提取淹没在噪声干扰中的动态压尤 为困难。
为解决上述问题,公开号为CN109363646的专利文献公开了一种三维 脉波影像的基面预处理方法及***,所述方法包括如下步骤:(1)数据 获取:三维脉波影像通过安装在机械手指上的传感器获得。(2)静态压 去除:将每个通道采集到的信号都减去所属通道的最小值,即可完成这个 处理。(3)经验模态分解:经验模态分解法的核心思想是假设信号分解 出来的每个模态都处于不同的中心频率附近,然后将这个问题转化为约束 问题来求解出每个模态。脉搏信号一般处于0.8到40Hz内,因此只需保留 合适的模态即可得到处理后的结果。但这些流程是对单传感器(等同于矩 阵式传感器中的单通道)采集的信号进行处理的,整个流程并不包含在 “面”水平上的去噪处理。尽管在数字信号处理的角度看来,这些方法的去 噪效果良好;然而,对于三维脉波影像,动态面的去噪效果却有待提高, 特别是,在存在个别传感器收集的信号被噪声淹没的动态面上,把脉触感 的立体表现严重受影响。
鲁棒主成分分析(Robustprincipalcomponentanalysis,RPCA)被广泛应 用于图像及视频去噪,现有研究成果表明,RPCA可以利用于三维脉波影 像的去噪
[LuJ,XieX,PengB,etal.InterferenceReductionbyUsingRPCAandVariationalMode Decompositionin3DPulseImages[C].InProceedingsofthe4thInternationalConferenceonMultimediaSystemsandSignalProcessing,2019],RPCA能将一系列图 像,分解为无噪声图像矩阵和噪声矩阵之和。其中,无噪声图像矩阵具有 低秩性,而噪声矩阵具有稀疏性。利用矩阵低秩性可以较好复原真实数据 矩阵,通过低秩矩阵来恢复原动态面,从而达到去除噪声的目的。然而, 在现实应用场景中,受测者的呼吸和运动产生的噪声,常常同时以相同的 规律出现在多个通道,并混杂于原始信号,有时甚至淹没原始信号。这 样,在低秩分解的过程中,RPCA将无法对这些通道的脉搏波信号及带有 一些规律的噪声进行正确的区分,也就难以完成信号提取和去除干扰的工 作。
综上,当前医学界亟待需要一种安全可靠的三维脉波影像去噪方法。
发明内容
发明的目的在于针对现有技术的不足之处,提供一种基于通道加权鲁 棒性主成分分析的三维脉波影像去噪方法。
本发明基于通道加权鲁棒性主成分分析的三维脉波影像去噪方法包 括:
获取三维脉波影像数据
其中,m为三维脉波影像的帧 数,n为采集三维脉波影像的矩阵式传感器的通道总个数,Xraw(t,j)则表 示第j个通道在第t帧测量到的压力值;
判断每个压力值的信号质量,并根据判断结果为每个压力值分配权 重,得到被加权的数据矩阵X=XrawWc;
对X=XrawWc分解得到一个低秩矩阵L和一个稀疏噪声矩阵S,即 X=L+S;
根据低秩矩阵
并结合双调和样条插值规则建立三维脉波 影像的曲面信息。
可选的,
步骤所述判断每个压力值的信号质量,并根据判断结果为每个压力值 分配权重,得到被加权的数据矩阵X=XrawWc包括:
使用四项高斯函数模拟在每个通道的每个心跳周期上,拟合压力振幅 值与时间的关系,并使用指标均方根误差描述拟合的效果,通过所述拟合 效果指代信号的质量。
可选的,
使用四项高斯函数模拟在每个通道的每个心跳周期上,拟合压力振幅 值与时间的关系,并使用指标均方根误差描述拟合的效果,通过所述拟合 效果指代信号的质量包括:
使用的四项高斯函数模型及RMSE的计算公式为:
其中C
r(t)为采集信号在单通道单周期上的幅值-时 间关系,而C(t)为对此进行拟合的结果函数;
其中RMSE(p,j)表示第j号通道在第p个周期上拟合结果的RMSE值,而 wc(j)则是计算出来赋予第j号通道的权值。
本发明还提供了一种基于通道加权鲁棒性主成分分析的三维脉波影像 去噪***,包括:
存储器、处理器和程序;
所述存储器用于存储所述程序;
所述处理器用于执行所述程序;
所述程序被配置为:
获取三维脉波影像数据
其中,m为三维脉波影像的帧 数,n为采集三维脉波影像的矩阵式传感器的通道总个数,Xraw(t,j)则表 示第j个通道在第t帧测量到的压力值;
判断每个压力值的信号质量,并根据判断结果为每个压力值分配权 重,得到被加权的数据矩阵X=XrawWc;
对X=XrawWc分解得到一个低秩矩阵L和一个稀疏噪声矩阵S,即 X=L+S;
根据低秩矩阵
并结合双调和样条插值规则建立三维脉波影 像的曲面信息。
可选的,
所述程序被配置为:
获取三维脉波影像数据
其中,m为三维脉波影像的帧 数,n为采集三维脉波影像的矩阵式传感器的通道总个数,Xraw(t,j)则表 示第j个通道在第t帧测量到的压力值;
使用四项高斯函数模拟在每个通道的每个心跳周期上,拟合压力振幅 值与时间的关系,并使用指标均方根误差描述拟合的效果,通过所述拟合 效果指代信号的质量;
对X=XrawWc分解得到一个低秩矩阵L和一个稀疏噪声矩阵S,即 X=L+S;
根据低秩矩阵
并结合双调和样条插值规则建立三维脉波 影像的曲面信息。
可选的,
所述程序被配置为:
获取三维脉波影像数据其中,m为三维脉波影像的帧 数,n为采集三维脉波影像的矩阵式传感器的通道总个数,Xraw(t,j)则表 示第j个通道在第t帧测量到的压力值;
使用的四项高斯函数模型及RMSE的计算公式为:
其中C
r(t)为采集信号在单通道单周期上的幅值-时 间关系,而C(t)为对此进行拟合的结果函数;
其中RMSE(p,j)表示第j号通道在第p个周期上拟合结果的RMSE值,而 wc(j)则是计算出来赋予第j号通道的权值对X=XrawWc分解得到一个低秩矩 阵L和一个稀疏噪声矩阵S,即X=L+S;
根据低秩矩阵
并结合双调和样条插值规则建立三维脉波 影像的曲面信息。
本发明相比于原有的三维脉波影像的预处理流程,添加在三维动态图 像层次的信号去噪的方法。相比于鲁棒性主成分分析(RPCA),本方法考 虑到各通道信号质量的差别,加入了通道加权,改善了去噪的效果。
具体实施方式
下面将结合本发明实施例,对本发明实施例中的技术方案进行清楚、 完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是 全部的实施例。在不脱离本发明的原理情况下,对本发明实施例做出变 化、修改、替换和变形都属于本发明保护的范围。
本发明基于通道加权鲁棒性主成分分析的三维脉波影像去噪方法实施 例具体包括:
获取三维脉波影像数据
其中,m为三维脉波影像的帧 数,n为采集三维脉波影像的矩阵式传感器的通道总个数,Xraw(t,j)则表 示第j个通道在第t帧测量到的压力值;
判断每个压力值的信号质量,并根据判断结果为每个压力值分配权 重,得到被加权的数据矩阵X=XrawWc;
对X=XrawWc分解得到一个低秩矩阵L和一个稀疏噪声矩阵S,即 X=L+S;
根据低秩矩阵
并结合双调和样条插值规则建立三维脉波影 像的曲面信息。
本发明的技术方案具体是这样实现的:
步骤1:输入的三维脉波影像数据为
其中,m为三维脉 波影像的帧数,n为采集三维脉波影像的矩阵式传感器的通道总个数 (n=12),Xraw(t,j)则表示第j个通道在第t帧测量到的压力值。
步骤2:通过右乘一个对角矩阵对输入数据Xraw进行列(通道)加权, 得到被通道加权的数据矩阵:
X=XrawWc (1)
很显然,这里有一个关键问题是权重矩阵Wc如何选择。这里我们借助 四项高斯函数模型确定权重值的大小:使用四项高斯函数模拟在每个通道 的每个心跳周期上,拟合压力振幅值与时间的关系;并使用指标均方根误 差(Root Mean Squared error,RMSE)描述拟合的效果。拟合效果实际可 以一定程度上指代信号的质量。
使用的四项高斯函数模型及RMSE的计算公式为:
其中Cr(t)为采集信号在单通道单周期上的幅值-时间关系,而C(t)为对此进行拟合的结果函数。
计算权重矩阵C每个通道的权重值公式为:
Wc=diag(wc(1),wc(2),...,wc(n)) (5)
其中RMSE(p,j)表示第j号通道在第p个周期上拟合结果的RMSE值,而wc(j)则是 计算出来赋予第j号通道的权值。
步骤3:利用RPCA方法,将上一步的输入数据
分解为一个低秩 矩阵L和一个稀疏噪声矩阵S,即X=L+S。通过解决凸优化问题,矩阵X就 能够被分解为低秩矩阵L及稀疏噪声矩阵S,该凸优化问题的拉格朗日乘项 表达如下,
式中:||S||1——l1范数,即矩阵元素绝对值之和;λ——平衡因子,在 合适的平衡因子λ的作用下,我们希望能够精确地恢复出变量L。
对式(6)构建增广拉格朗日函数如下:
式中:Y——拉格朗日算子,μ>0为常数。
式(7)可以使用交替方向乘子法(ADMM)进行优化求解,输出低秩 矩阵L及稀疏噪声矩阵S。
步骤4:对步骤3输出的二维低秩矩阵
根据步骤2的加权步 骤进行逆向消除:
步骤5:将步骤4输出的二维低秩矩阵
使用双调和样条插值法 建立三维脉波影像的曲面信息,最后的图像就是去噪的结果。
本发明还提供了一种基于通道加权鲁棒性主成分分析的三维脉波影像 去噪***,包括:
存储器、处理器和程序;
所述存储器用于存储所述程序;
所述处理器用于执行所述程序;
所述程序被配置为:
获取三维脉波影像数据
其中,m为三维脉波影像的帧 数,n为采集三维脉波影像的矩阵式传感器的通道总个数,Xraw(t,j)则表 示第j个通道在第t帧测量到的压力值;
判断每个压力值的信号质量,并根据判断结果为每个压力值分配权 重,得到被加权的数据矩阵X=XrawWc;
对X=XrawWc分解得到一个低秩矩阵L和一个稀疏噪声矩阵S,即 X=L+S;
根据低秩矩阵并结合双调和样条插值规则建立三维脉波影 像的曲面信息。
使用峰值信噪比(Peak Signal to Noise Ratio,PSNR)来对各个方法进 行衡量(值越大,效果越好):
其中C为原始三维脉波影像,T为经各种方法去噪后的结果,m,n分别为 三维脉波影像的帧数,及传感器个数。
表1三种不同的方法对观测信号进行处理后的结果峰值信噪比
从表1中可以看出,通道加权鲁棒性主成分分析CW-RPCA去噪效果 最好。
以上所述,仅为本发明较优实施例之一,因此,在不脱离本发明的原 理情况下,对本发明实施例做出变化、修改、替换和变形都属于本发明保 护的范围。