CN113729738A - 一种多通道肌电特征图像的构建方法 - Google Patents

Abstract

本发明是一种多通道肌电特征图像数据集的构建方法。首先，结合阈值比较法和巴特沃斯滤波器，完成原始肌电信号的预处理。随后，从常用的时域和频域特征中进行2次筛选，选择4种互不冗余的肌电特征进行特征提取。其次，基于一维信号与二维图像的映射关系生成肌电图像。最后，按照图像的拼接方式，完成多通道肌电特征图像的构建。肌电图像数据集通过深度学习网络的训练，可以有效提高手势识别率。多通道肌电特征图像具有更加丰富的特征信息，多种特征可以完成信息之间的互补，最后的识别率会比单通道肌电特征图像高出8％‑9％。

Description

一种多通道肌电特征图像的构建方法

技术领域

本发明涉及信号处理与人机交互领域，特别涉及一种基于肌电信号的双向递归肌电特征选择方法。

背景技术

随着科技的发展，人机交互的研究与应用受到越来越多的学者关注。其中，手势识别作为人机交互的一个主要分支，是人机交互领域的重要研究课题之一。人机交互中的手势对使用者来说具有便捷、交互性强、表达内容丰富等优点，一直是人机交互领域的首要选择。肌电信号是在人体肌肉活动过程中产生的一种生物电信号，是不同的运动单位动作电势在肌肉纤维和细胞组织中传导叠加的综合结果。肌电信号主要是通过提取肌电特征进行手势识别，但基于传统方法的肌电手势分类识别率一直处于停滞的状态。

鉴于此需要结合深度学习的方法进行肌电手势识别，将传统方法与深度学习结合，使用传统方法提取不同表征的肌电特征并生成肌电图像，通过深度学习理解特征图中包含的更抽象的信息，最后完成信息的融合，以期取得更佳的分类性能。

发明内容

本发明提供了一种多通道肌电特征图像的构建方法，可以从多种肌电特征中，通过肌电特征提取和一维信号到二维图像的映射数学关系，完成单通道肌电特征图像的转化，然后基于图像通道的拼接方法，完成多通道肌电特征图像的构建。

为解决上述技术问题，根据发明内容，采用如下技术方案为一种多通道肌电特征图像的构建方法，具体包括以下步骤：

步骤1：获取人体前臂表面肌电信号sEMG，完成活动段提取和滤波降噪预处理工作；

步骤2：选择合适的时域特征和频域特征；

步骤3：基于一维信号与二维图像的映射关系生成单通道肌电特征图像Image_单通道；

步骤4：对比初步筛选特征的图像，利用相似度分析，筛选出4种特征图像；

步骤5：结合图像通道拼接方法，完成多通道肌电特征图像Image_多通道，其中筛选出的4种特征图像为其4个通道；

步骤6：重复步骤3–步骤5，完成信号向图片的转换，构建多通道肌电特征图像数据集。

进一步地，步骤1中利用阈值比较法提取活动段S(n)：

其中，C为sEMG的采集通道，N为采样点的个数，sEMG_c(n)为c通道的第n个采样点的值，

为c通道肌肉放松时sEMG的平均值，TH为设定的阈值。

通过1HZ的一阶巴特沃斯滤波器进行低通滤波，其中，

为振幅，

为频率，ε为比例系数，n为滤波器的阶数，

为截止频率，

为通频带边缘频率。

进一步地，步骤3中所述单通道肌电特征图像Image_单通道为：

Image_单通道∈S^W×H×1

W表示滑动窗口的宽，也是图像的宽，H表示原始肌电信号的通道数，也表示图像或滑动窗口高度。使用滑动窗口完成一维信号到[0，1]区间内的灰度图像映射，映射关系式为：

A_[0,1]＝F(A_(i,j))

F(A_(i,j))＝(A_(i,j)-A_min)/(A_max-A_min)

A为sEMG通过滑动窗口后所形成的sEMG矩阵，A_[0,1]为转换后的肌电图像，F()为傅里叶函数，A_(i,j)为sEMG矩阵中的第i行第j列的值，其中，0<i<H，0<j<W，A_max为sEMG矩阵中最大的值，A_min为sEMG矩阵中最小的值。

进一步地，步骤5中所述多通道肌电特征图像Image_多通道为：

Image_多通道∈S^W×H×4

多通道肌电特征图像的每个通道由4种单通道肌电特征图像组成，每张单通道肌电特征图像尺寸为W×H，C为图像的通道数，也是特征数量，组合在一起构成的多通道肌电特征图像尺寸为W×H×C，C＝4。

与现有技术相比，本发明的有益效果是：

本发明提出的一种多通道肌电特征图像的构建方法，基于信号与图像的映射关系和图像通道拼接原理，将一维肌电信号转换成二维肌电图像，有助于构建肌电特征图像数据集。肌电图像数据集通过深度学习网络的训练，可以有效提高手势识别率。多通道肌电特征图像具有更加丰富的特征信息，多种特征可以完成信息之间的互补，最后的识别率会比单通道肌电特征图像高出8％-9％。

附图说明

图1：为本发明方法的流程图；

图2：为采集肌电信号的手势图；

图3：为多通道肌电特征图像的构建原理图；

图4：为肌电信号采集设备。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

本发明的具体实施方式采用如下技术方案为一种多通道肌电特征图像的构建方法，具体包括以下步骤：

步骤1：通过如图4所示的肌电信号采集设备获取如图2所示动作的前臂表面肌电信号，这9个手势包括手掌闭合(SH)和手掌张开(SK)，手腕运动包括手腕屈曲(NQ)和外屈曲(WQ)，手指运动包括拇指对食指(MS)、中指(MZ)、无名指(MW)和小指(MX)的作用力，除此之外，测试手势动作还包括一个休息动作(RE)。完成肌电信号采集工作后，进行信号活动段提取和滤波降噪等预处理工作。活动段指动作执行过程中的肌肉收缩信号，活动段提取的目的是为了去除不包含肌肉收缩的非活动段；

步骤1中利用阈值比较法提取活动段S(n)：

其中，C为sEMG的采集通道，N为采样点的个数，sEMG_c(n)为c通道的第n个采样点的值，

为c通道肌肉放松时sEMG的平均值，TH为设定的阈值。

通过1HZ的一阶巴特沃斯滤波器进行低通滤波，其中，

为振幅，

为频率，ε为比例系数，n为滤波器的阶数，

为截止频率，

为通频带边缘频率。

步骤2：选择4种肌电特征，分别为：均方根(RMS)、波形长度(WL)、样本熵(SampEn)以及中值幅谱(MAS)；

步骤3：基于一维信号与二维图像的映射关系生成单通道肌电特征图像；

步骤3中所述单通道肌电特征图像Image_单通道为：

Image_单通道∈S^W×H×1

W表示滑动窗口的宽，也是图像的宽，H表示原始肌电信号的通道数，也表示图像或滑动窗口高度。使用滑动窗口完成一维信号到[0，1]区间内的灰度图像映射，映射关系式为：

A_[0,1]＝F(A_(i,j))

F(A_(i,j))＝(A_(i,j)-A_min)/(A_max-A_min)

A为sEMG通过滑动窗口后所形成的sEMG矩阵，A_[0,1]为转换后的肌电图像，F()为傅里叶函数，A_(i,j)为sEMG矩阵中的第i行第j列的值，其中，0<i<H，0<j<W，A_max为sEMG矩阵中最大的值，A_min为sEMG矩阵中最小的值。

步骤4：对比初步筛选特征的图像，利用相似度分析，筛选出4种特征图像；

步骤5：结合图像通道拼接方法，完成多通道肌电特征图像，其中筛选出的4种特征图像为其4个通道。

步骤5中所述多通道肌电特征图像Image_多通道为：

Image_多通道∈S^W×H×4

多通道肌电特征图像的每个通道由4种单通道肌电特征图像组成，每张单通道肌电特征图像尺寸为W×H，C为图像的通道数，也是特征数量，组合在一起构成的多通道肌电特征图像尺寸为W×H×C，C＝4，其生成原理如图3所示。

本说明中未作详细描述的内容属于本领域专业技术人员公知的现有技术，尽管参照前述实施例对本发明进行了详细的说明，对于本领域的技术人员来说，其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (4)

1.一种多通道肌电特征图像的构建方法，其特征在于，所述方法包括如下步骤：

步骤1：获取人体前臂表面肌电信号sEMG，完成活动段提取和滤波降噪预处理工作；

步骤2：选择合适的时域特征和频域特征；

步骤3：基于一维信号与二维图像的映射关系生成单通道肌电特征图像Image_单通道；

步骤4：对比初步筛选特征的图像，利用相似度分析，筛选出4种特征图像；

步骤5：结合图像通道拼接方法，完成多通道肌电特征图像Image_多通道，其中筛选出的4种特征图像为其4个通道；

步骤6：重复步骤3–步骤5，完成信号向图片的转换，构建多通道肌电特征图像数据集。

2.如权利要求1所述的一种多通道肌电特征图像的构建方法，其特征在于，步骤1中利用阈值比较法提取活动段S(n)：

其中，C为sEMG的采集通道，N为采样点的个数，sEMG_c(n)为c通道的第n个采样点的值，

为c通道肌肉放松时sEMG的平均值，TH为设定的阈值；

通过1HZ的一阶巴特沃斯滤波器进行低通滤波，

其中，

为振幅，

为频率，ε为比例系数，n为滤波器的阶数，

为截止频率，

为通频带边缘频率。

3.如权利要求1所述的一种多通道肌电特征图像的构建方法，其特征在于，步骤3中所述单通道肌电特征图像Image_单通道为：

Image_单通道∈S^W×H×1

W表示滑动窗口的宽，也是图像的宽，H表示原始肌电信号的通道数，也表示图像或滑动窗口高度，使用滑动窗口完成一维信号到[0，1]区间内的灰度图像映射，映射关系式为：

A_[0,1]＝F(A_(i,j))

F(A_(i,j))＝(A_(i,j)-A_min)/(A_max-A_min)

A为sEMG通过滑动窗口后所形成的sEMG矩阵，A_[0,1]为转换后的肌电图像，F()为傅里叶函数，A_(i,j)为sEMG矩阵中的第i行第j列的值，其中，0<i<H，0<j<W，A_max为sEMG矩阵中最大的值，A_min为sEMG矩阵中最小的值。

4.如权利要求1所述的一种多通道肌电特征图像的构建方法，其特征在于，步骤5中所述多通道肌电特征图像Image_多通道为：

Image_多通道∈S^W×H×4

多通道肌电特征图像的每个通道由4种单通道肌电特征图像组成，每张单通道肌电特征图像尺寸为W×H，C为图像的通道数，也是特征数量，组合在一起构成的多通道肌电特征图像尺寸为W×H×C，C＝4。

Images