CN112057080A - 一种基于分阶段特征提取的冻结步态检测方法和系统 - Google Patents

一种基于分阶段特征提取的冻结步态检测方法和系统 Download PDF

Info

Publication number
CN112057080A
CN112057080A CN202010798417.2A CN202010798417A CN112057080A CN 112057080 A CN112057080 A CN 112057080A CN 202010798417 A CN202010798417 A CN 202010798417A CN 112057080 A CN112057080 A CN 112057080A
Authority
CN
China
Prior art keywords
gait
sample
feature
acceleration
jerk
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202010798417.2A
Other languages
English (en)
Other versions
CN112057080B (zh
Inventor
赵金
任康
刘西华
周洋
凌云
陈仲略
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huazhong University of Science and Technology
Original Assignee
Huazhong University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Huazhong University of Science and Technology filed Critical Huazhong University of Science and Technology
Priority to CN202010798417.2A priority Critical patent/CN112057080B/zh
Publication of CN112057080A publication Critical patent/CN112057080A/zh
Application granted granted Critical
Publication of CN112057080B publication Critical patent/CN112057080B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
    • A61B5/112—Gait analysis
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235—Details of waveform analysis
    • A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • A61B5/7267—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
    • G06F18/2135—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on approximation criteria, e.g. principal component analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N20/00—Machine learning
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/02—Preprocessing
    • G06F2218/04—Denoising
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/08—Feature extraction
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/12—Classification; Matching

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Pathology (AREA)
  • Software Systems (AREA)
  • Veterinary Medicine (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Public Health (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Physiology (AREA)
  • General Health & Medical Sciences (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Evolutionary Biology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Signal Processing (AREA)
  • Psychiatry (AREA)
  • Fuzzy Systems (AREA)
  • Computing Systems (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

本发明公开了一种基于分阶段特征提取的冻结步态检测方法和系统，属于机器学习领域。包括：基于用户行走过程中采集到的原始加速度数据，构建有标签的样本集，样本为加窗处理后的加速度数据序列，标签表示样本属于正常步态还是冻结步态；对每个样本进行分阶段特征提取；使用PCA对步态特征集进行特征变换，得到低维新步态特征集，对低维新步态特征集进行特征选择，得到最优步态特征子集；使用最优步态特征子集训练基于机器学习的冻结步态检测模型；提取待测样本的分阶段特征，输入至训练好的冻结步态检测模型，得到冻结步态检测结果。本发明提取加速度信号中运动分量及冻结带，并对运动信号进行合成与分解，以充分发挥原始数据的潜在特征。

Description

一种基于分阶段特征提取的冻结步态检测方法和系统
技术领域
本发明属于机器学习领域，更具体地，涉及一种基于分阶段特征提取的冻结步态检测方法和系统。
背景技术
帕金森病(Parkinson’sdisease，PD)是一种常见的神经系统变性疾病，目前我国PD患者已高达290万，并以约每年10万新发患者的速度持续增长。随着人口老龄化比率不断升高，预计到2030年，我国的PD患者将达到500万。冻结步态(Freezing Of Gait，FoG)是中晚期PD很常见的致残性症状，严重影响病人的活动能力。FoG事件的早期检测作为干预措施实施的基础，能帮助发病患者减少跌倒风险，恢复步行及正常活动；作为后续病情评估的依据可提供相关症状信息，对冻结步态的研究与治疗有着重要的指导作用。
现有技术中关于帕金森患者冻结步态的检测方法主要大部分是直接对预处理后的信号进行特征提取，然后将特征集用于分类器的训练，存在特征不明显、鲁棒性低、泛化能力不强等问题。
发明内容
针对现有技术的缺陷和改进需求，本发明提供了一种基于分阶段特征提取的冻结步态检测方法和系统，其目的在于通过提取加速度信号中对运动状态有决定性作用的运动分量及表征FoG的冻结带，并对运动信号进行合成与分解，以充分发挥原始数据的潜在特征，提高检测效率和精度，增强算法鲁棒性和泛化能力。
为实现上述目的，按照本发明的第一方面，提供了一种基于分阶段特征提取的冻结步态检测方法，该方法包括以下步骤：
S1.基于用户行走过程中采集到的原始加速度数据，构建有标签的样本集，其中，样本为加窗处理后的加速度数据序列，标签表示样本属于正常步态还是冻结步态；
S2.对每个样本进行分阶段特征提取：将样本加速度分解为运动分量、重力分量和冻结带；基于样本加速度计算第一急动度；基于运动分量计算第二急动度；基于冻结带加速度计算第三急动度；将样本加速度、第一急动度、第二急动度和第三急动度的三轴数据分别进行合成，得到合成加速度、第一合成急动度、第二合成急动度和第三合成急动度；分别提取合成加速度、第一合成急动度、第二合成急动度和第三合成急动度的时域特征和频域特征，得到步态特征集；
S3.使用主成分分析法对步态特征集进行特征变换，得到代表原始步态特征集的低维新步态特征集，采用基于相关系数或者最大信息系数法对低维新步态特征集进行特征选择，得到最优步态特征子集；
S5.提取待测样本的分阶段特征，输入至训练好的冻结步态检测模型，得到冻结步态检测结果。
优选地，步骤S1中，构建样本集之前，对传感器获取的原始数据进行数据清洗、滤波降噪和动静态分离；
所述滤波降噪，首先使用巴特沃斯低通滤波消除数据清洗处理后信号中的尖峰噪声，再使用滑动平均滤波消除随机噪声；
所述动静态分离，首先计算滤波降噪处理后信号的急动度，再根据急动度的大小，剔除单次试验前后的非运动态加速度数据，得到运动态加速度数据。
优选地，所述根据急动度的大小，剔除单次试验前后的非运动态加速度数据，得到运动态加速度数据，具体包括以下步骤：
I.将滤波降噪处理后加速度数据划分为M个长度为N的连续序列，M与N满足以下关系：
Figure BDA0002625816920000031
其中，Lacc为同一传感器所采集加速度经过滤波降噪处理后的序列长度；
II.计算每个序列的特征值T，计算公式如下：
Figure BDA0002625816920000032
其中，Jxi、Jyi、Jzi分别表示采样点i处x，y，z轴的急动度；
III.对于每个序列的特征值T，设置阈值Tsh，将T＞Tsh的第一个序列记为运动步态开始序列，T＞Tsh的最后一个序列记为运动步态结束序列。
优选地，对运动态加速度数据进行加窗处理，将每个窗口的数据序列
Figure BDA0002625816920000033
作为一个样本，采样点序号sw＝1+(w-1)×S，窗口序号w＝{1，…，m}，S为步长，m为划分窗口总数量，若窗口内正常步态所占比例大于0.5，则样本标签为正常步态，否则，样本标签为冻结步态。
优选地，步骤S1中，构建样本集时，采用SMOTE过采样算法增加冻结步态样本的数量，保证正常步态样本数与冻结步态样本数之比为2∶1。
优选地，运动加速度的运动分量计算过程如下：
(1)采用4阶巴特沃斯高通滤波器分别对降噪后每个维度的加速度滤波，得到高频分量{ahx，ahy，ahz}，截止频率ωn＝0.5Hz；
(2)采用4阶巴特沃斯低通滤波器分别对降噪后每个维度的加速度滤波，得到低频分量{alx，aly，alz}，截止频率ωn＝0.5Hz；
(3)分别计算高频分量的欧几里得范数HFEN和低频分量的欧几里得范数LFEN
Figure BDA0002625816920000041
Figure BDA0002625816920000042
(4)按照公式MAcc＝HFEN+(LFEN-1)计算运动分量。
优选地，所述时域特征包括：均值、标准差、窗口内最大值、窗口内最小值、窗口内最大值与最小值的差值、窗口内超过均值点的个数、第三四分位数与第一四分位数的差距、窗口内数据平方的均值、XY两轴之间的相关性、XZ两轴之间的相关性、YZ两轴之间的相关性和信号幅值面积；所述频域特征包括：直流分量、功率谱密度的幅值统计特征、功率谱密度的形状统计特征和冻结指数FI。
优选地，基于相关系数的特征选择步骤为：
(1)计算低维新步态特征集各特征间的皮尔逊相关系数，获得相关性矩阵；
Figure BDA0002625816920000043
其中，r为特征
Figure BDA0002625816920000044
与
Figure BDA0002625816920000045
间皮尔逊相关系数，特征
Figure BDA0002625816920000046
与
Figure BDA0002625816920000047
对应样本向量分别为
Figure BDA0002625816920000048
(2)设置相关性选择阈值k，k∈[0，1]；
(3)计算各特征的相关性总和，即对相关性矩阵按列求和；
(4)筛选出相关性矩阵中取值大于k的特征对，删除其中相关性总和较大的特征量；
(5)剩余特征组成新的特征集。
优选地，基于最大信息系数法的特征选择步骤为：
(1)计算所有特征Xl与步态类别Y之间的MIC值，选择MIC最大值对应的特征为特征子集的第一个特征；
(2)选择J最大值对应的特征加入特征子集；
Figure BDA0002625816920000051
其中，Xl为候选特征，XF为已选择的特征集；
(3)重复步骤(2)直至所选特征个数达到预设值。
为实现上述目的，按照本发明的第二方面，提供了一种基于分阶段特征提取的冻结步态检测系统，包括：计算机可读存储介质和处理器；
所述计算机可读存储介质用于存储可执行指令；
所述处理器用于读取所述计算机可读存储介质中存储的可执行指令，执行第一方面所述的基于分阶段特征提取的冻结步态检测方法。
总体而言，通过本发明所构思的以上技术方案，能够取得以下有益效果：
(1)针对现有基于机器学习的FoG检测研究通常直接从加速度信号中提取时频域特征用于分类器训练，这种方法没有挖掘加速度信号特征及FoG特性，本发明提出种分阶段特征提取，提取加速度信号中对运动状态有决定性作用的运动分量及表征FoG的冻结带，并对运动信号进行合成与分解，以充分发挥原始数据的潜在特征，并进一步设计了二十种时频域特征用于充分表征FoG。
(2)针对过大的特征维度或计算步骤十分复杂的特征可能会导致无法接受的延迟问题，本发明将一种特征变换和两种特征选择方法用于从较大维度特征中选择最佳特征子集，并对相关系数法和最大信息系数法两种特征选择方法进行改进，将改进后的特征选择算法及主成分分析法用于从较大维度特征中选择最佳特征子集，从而降低计算压力。
(3)本发明将FoG检测问题设定为基于加速度信号的有监督的二分类学习问题，需要构建有标签的原始加速度数据集，通过对传感器采集的信号进行了简单的处理，包括时间戳对齐、缺失值补齐等数据清洗过、滤波降噪过程和动静态分离过程，在此基础上进行通过划分时间窗和正负样本均衡构建原始样本集，从而进一步提高步态检测精度和速度。
附图说明
图1为本发明提供的一种基于分阶段特征提取的冻结步态检测方法流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
如图1所示，本发明提供一种基于分阶段特征提取的冻结步态检测方法，该方法包括以下步骤：
步骤S1.基于用户行走过程中采集到的原始加速度数据，构建有标签的样本集，其中，样本为加窗处理后的加速度数据序列，标签表示样本属于正常步态还是冻结步态。
优选地，步骤S1中，将7个三轴加速度传感器分别固定于用户腰部、左右大腿膝盖上方、左右小腿脚踝上方和左右脚鞋子侧外面，用于采集用户行走过程中的原始加速度数据。
优选地，步骤S1中，构建样本集之前，对传感器获取的原始数据进行数据清洗、滤波降噪和动静态分离，所述数据清洗包括：
(1)时间戳对齐：将最后与数据接收装置建立连接的传感器采集数据的起始时间作为整体起始时间，删除其他传感器在整体起始时间之前的数据。
(2)缺失值补齐：采用邻近点均值法进行加速度数据缺失值的填补。
(3)异常值处理：采用基于距离的离群点检测方法检测离群数据，对检测到的离群数据采用热卡填充进行重构。
所述滤波降噪，首先使用巴特沃斯低通滤波消除数据清洗处理后信号中的尖峰噪声，再使用滑动平均滤波消除随机噪声。
巴特沃斯滤波适用于处理尖峰噪声。由于99％人类运动步态频率低于20Hz，本实施例采用三阶巴特沃斯低通滤波器滤除20Hz以上高频部分。
Figure BDA0002625816920000071
其中，n为滤波器阶数，ωc为截止频率，ωp为通频带边缘频率，
Figure BDA0002625816920000072
为通频带边缘的数值。
滑动平均滤波平滑度高，具有中和随机噪声的作用。其原理为将小段连续序列视为接近平稳，对该序列作局部平均，以减小随机误差造成的起伏。取固定队列长度为N，第t个序列内动态采集数据实际值为{xt-(N-1)，xt-(N-2)，…，xt}根据先进先出原则，每次将新数据放入队尾，并扔掉原来队首的数据，对队列中的N个数据计算算数平均值，得到下一期的预测值
Figure BDA0002625816920000073
但此种算法会丢失t时刻后的数据，因而在此基础上进行了滑动平均滤波算法的改进，取预测时刻前后各M个样本(即序列长度为2M+1)，此时第t个序列内动态采集数据实际值为{xt-M，…，xt-1，xt，xt+1，…，xt+M}根，t时刻预测值
Figure BDA0002625816920000074
所述动静态分离，首先计算滤波降噪处理后信号的急动度，再根据急动度的大小，剔除单次试验前后的非运动态加速度数据，得到运动态加速度数据。
患者运动时加速度值发生较大突变，而静止时加速度值的变化率较小，可以通过计算加速度的变化率来分离动静态。加速度的变化率也称为急动度，根据急动度的大小，将单次试验前后的静态数据与动态数据分离开，便于后续对动态数据进行特征提取。
急动度反应的是加速度的变化率，根据采集到的加速度数据，可由
Figure BDA0002625816920000081
计算其急动度，其中，t2，t1为两个连续采样点对应时刻，Acct2，Acct1分别为t2，t1时刻的加速度值，Δt为采样间隔时间。
优选地，所述根据急动度的大小，剔除单次试验前后的非运动态加速度数据，得到运动态加速度数据，具体包括以下步骤：
I.将滤波降噪处理后加速度数据划分为M个长度为N的连续序列，M与N满足以下关系：
Figure BDA0002625816920000082
其中，Lacc为同一传感器所采集加速度经过滤波降噪处理后的序列长度；
II.计算每个序列的特征值T，计算公式如下：
Figure BDA0002625816920000083
其中，Jxi、Jyi、Jzi分别表示采样点i处x，y，z轴的急动度；
III.对于每个序列的特征值T，设置阈值Tsh，将T＞Tsh的第一个序列记为运动步态开始序列，T＞Tsh的最后一个序列记为运动步态结束序列。
优选地，对运动态加速度数据进行加窗处理，将每个窗口的数据序列
Figure BDA0002625816920000084
作为一个样本，采样点序号sw＝1+(w-1)×S，窗口序号w＝{1，…，m}，S为步长，m为划分窗口总数量，若窗口内正常步态所占比例大于0.5，则样本标签为正常步态，否则，样本标签为冻结步态。
窗口长度N的选择没有限制，由于输入信号中采样点的数量将决定输出频率窗口的分辨率，N越大，对计算涉及离散傅里叶变换的频域特征越有利。但长窗口会降低时间分辨率，且在窗口内无法检测短事件，而FoG持续事件往往较短。另外，窗口内数据量越大，处理时间越长，可能会在数据采集和分类之间引起不必要的延迟，从而影响检测效果。
优选地，步骤S1中，构建样本集时，采用SMOTE过采样算法增加冻结步态样本的数量，保证正常步态样本数与冻结步态样本数之比为2∶1。
划分滑窗后，存在正常步态样本与FoG样本比例不均衡的现象，这将导致分类算法过多的关注正常步态，而使FoG样本的分类性能下降。为了保证数据集的无偏性，本发明采用SMOTE过采样算法增加FoG样本的数量来平衡数据集。使用SMOTE算法将正常步态样本数与FoG样本数之比填充至2∶1，正负样本不均衡问题得到了很好的解决。
步骤S2.对每个样本进行分阶段特征提取：将样本加速度分解为运动分量、重力分量和冻结带；基于样本加速度计算第一急动度；基于运动分量计算第二急动度；基于冻结带加速度计算第三急动度；将样本加速度、第一急动度、第二急动度和第三急动度的三轴数据分别进行合成，得到合成加速度、第一合成急动度、第二合成急动度和第三合成急动度；分别提取合成加速度、第一合成急动度、第二合成急动度和第三合成急动度的时域特征和频域特征，得到步态特征集。
本发明提出一种分阶段的特征提取算法，提取加速度信号中对运动状态有决定性作用的运动分量及表征FoG的冻结带，并对运动信号进行合成与分解，以充分发挥原始数据的潜在特征。
加速度的运动分量是反应受试者运动状态的关键，本发明提出一项用于度量运动加速度分量的指标。优选地，运动加速度的运动分量计算过程如下：
(1)采用4阶巴特沃斯高通滤波器分别对降噪后每个维度的加速度滤波，得到高频分量{ahx，ahy，ahz}，截止频率ωn＝0.5Hz；
由于与大多数人体日常运动相关的加速度频率大于0.5Hz，将滤波器的截止频率设置为0.5Hz。
(2)采用4阶巴特沃斯低通滤波器分别对降噪后每个维度的加速度滤波，得到低频分量{alx，aly，alz}，截止频率ωn＝0.5Hz；
(3)分别计算高频分量的欧几里得范数HFEN和低频分量的欧几里得范数LFEN
Figure BDA0002625816920000101
Figure BDA0002625816920000102
(4)按照公式MAcc＝HFEN+(LFEN-1)计算运动分量。
如果没有旋转运动，则低频分量的欧几里得范数(LFEN)等于重力的大小(1g)，而旋转运动过程中，由于分离不完善，LFEN可能不为1g，当运动感加速度(离心加速度)产生低于截止频率的低频分量(荡秋千时的向心力)时，LFEN＞1g，当部分重力分量转移至高频分量中(垂直平面上的旋转运动)时，LFEN＜1g，因而将低频分量LFEN与1g的差值添加到HFEN中。
冻结步态会在加速度中产生3-8Hz的谐波，故在提取运动加速度分量后，采用巴特沃斯带通滤波器分离出3-8Hz冻结带加速度。
考虑到加速度计采集xyz三轴数据，分别表示三个方向的加速度，而在FoG检测过程中，没有需要特定单轴来辨别方向的动作，故将三轴加速度(急动度)合成为一个加速度(急动度)值，称为合成加速度(急动度)，以合成加速度(急动度)为初始样本进行特征提取和分类，能在保证精度的同时，减少计算的复杂性。根据
Figure BDA0002625816920000103
计算合成加速度与合成急动度ai，axi，ayi，azi分别表示第i个采样点处x，y，z轴的加速度或急动度，其中，i∈{1，2，…，n}，n为采样点个数。
优选地，所述时域特征包括：均值、标准差、窗口内最大值、窗口内最小值、窗口内最大值与最小值的差值、窗口内超过均值点的个数、第三四分位数与第一四分位数的差距、窗口内数据平方的均值、XY两轴之间的相关性、XZ两轴之间的相关性、YZ两轴之间的相关性和信号幅值面积；所述频域特征包括：直流分量、功率谱密度的幅值统计特征、功率谱密度的形状统计特征和冻结指数FI。
本实施例选取了如表1所示的10种时域特征和如表2所示的10种频域特征，这些特征旨在表现FoG与正常步态的不同特点。对于时域特征，均值、方差、标准差、最大值、最小值、过零点个数、四分位差、范围等用于表征运动量及统计力矩，信号幅值面积反应行走速度，信号熵可以体现信息的复杂度且可以较好地衡量非线性信号，此外，相关性用于表征各个运动轴之间的相关性，需要用到各个维度的加速度数据。频域特征用来发现信号中的周期性信息，FoG和正常步态都是典型的周期性运动，可以很好的被频域特征表征。频域信号通过对时域加速度信号进行快速傅里叶变换后得到。频域特征主要包括直流分量和功率谱密度，对功率谱密度进一步提取了幅度统计特征和形状统计特征。此外，还计算了冻结指数FI，其定义为加速度中0.5-3Hz的频率分量，冻结频带定义为3-8Hz的频率分量，将时刻t的冻结指数(freezingindex,FI)定义为以时间t为中心的6s数据窗口中冻结频带功率谱下面积的平方除以运动频带功率谱下面积的平方。
表1
Figure BDA0002625816920000111
Figure BDA0002625816920000121
表2
Figure BDA0002625816920000122
Figure BDA0002625816920000131
步骤S3.使用主成分分析法对步态特征集进行特征变换，得到代表原始步态特征集的低维新步态特征集，采用基于相关系数或者最大信息系数法对低维新步态特征集进行特征选择，得到最优步态特征子集。
经过分阶段提取步态特征后，共获取了高维运动特征，而典型的机器学习算法无法有效的处理如此大量的输入特征，因此需要降低特征维度。本发明选用一种特征变换和两种特征选择方法进行FoG特征降维，其中，特征变换方法为主成分分析法，考虑到FoG检测对算法实时性的要求，且经过分阶段提取后的特征维度较高，采用过滤法进行特征选择。FoG特征均为时间序列下的连续特征，故本发明选择适用于连续特征分类的基于相关系数的特征选择和最大信息系数法，并分别对两种方法进行FoG适应性改进。
优选地，基于相关系数的特征选择步骤为：
(1)计算低维新步态特征集各特征间的皮尔逊相关系数，获得相关性矩阵；
Figure BDA0002625816920000141
其中，r为特征
Figure BDA0002625816920000142
与
Figure BDA0002625816920000143
间皮尔逊相关系数，特征
Figure BDA0002625816920000144
与
Figure BDA0002625816920000145
对应样本向量分别为
Figure BDA0002625816920000146
(2)设置相关性选择阈值k，k∈[0，1]；
(3)计算各特征的相关性总和，即对相关性矩阵按列求和；
(4)筛选出相关性矩阵中取值大于k的特征对，删除其中相关性总和较大的特征量；
(5)剩余特征组成新的特征集。
本实施例中相关性选择阈值k＝0.8。
优选地，基于最大信息系数法的特征选择步骤为：
(1)计算所有特征Xl与步态类别Y之间的MIC值，选择MIC最大值对应的特征为特征子集的第一个特征；
将特征之间的归一化互信息最大的网格分辨率作为MIC的度量值，不同特征之间MIC值越大，说明这两个特征相关性越高。
(2)选择J最大值对应的特征加入特征子集；
Figure BDA0002625816920000147
其中，Xl为候选特征，XF为已选择的特征集；
(3)重复步骤(2)直至所选特征个数达到预设值。
本实施例从高维特征中选出100个最有效特征构建最优步态特征子集。
步骤S4.使用最优步态特征子集训练基于机器学习的冻结步态检测模型，训练过程中，冻结步态检测模型接收样本进行学习后对样本标签做出预测，将预测结果与实际标签进行比较后修正冻结步态检测模型。
为实现FoG的高效精确检测，本发明选用LigthGBM、支持向量机或者随机森林三种分类器进行正常步态与FoG的划分，模型训练过程中采用十折交叉验证法对超参数进行选择和评估。通过对模型测试结果的分析，评估不同窗口大小、特征选择方法和分类算法对模型效果的影响。
步骤S5.提取待测样本的分阶段特征，输入至训练好的冻结步态检测模型，得到冻结步态检测结果。
本发明采用十折交叉验证法对各分类器的超参数进行了测试。表3列出了本实施例各分类器用于十折交叉验证评估的超参数。
表3
Figure BDA0002625816920000151
Figure BDA0002625816920000161
经实验验证，窗口大小为256、使用相关系数法进行特征选择、LightGBM为分类器时，获得最佳分类效果，其精确率为94.66％，召回率为94.66％，准确率为94.66％、FI指数为94.66％、AUC为98.48％，精确率与召回率有较好的平衡。
本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种基于分阶段特征提取的冻结步态检测方法，其特征在于，该方法包括以下步骤：
S1.基于用户行走过程中采集到的原始加速度数据，构建有标签的样本集，其中，样本为加窗处理后的加速度数据序列，标签表示样本属于正常步态还是冻结步态；
S2.对每个样本进行分阶段特征提取：将样本加速度分解为运动分量、重力分量和冻结带；基于样本加速度计算第一急动度；基于运动分量计算第二急动度；基于冻结带加速度计算第三急动度；将样本加速度、第一急动度、第二急动度和第三急动度的三轴数据分别进行合成，得到合成加速度、第一合成急动度、第二合成急动度和第三合成急动度；分别提取合成加速度、第一合成急动度、第二合成急动度和第三合成急动度的时域特征和频域特征，得到步态特征集；
S3.使用主成分分析法对步态特征集进行特征变换，得到代表原始步态特征集的低维新步态特征集，采用基于相关系数或者最大信息系数法对低维新步态特征集进行特征选择，得到最优步态特征子集；
S4.使用最优步态特征子集训练基于机器学习的冻结步态检测模型，训练过程中，冻结步态检测模型接收样本进行学习后对样本标签做出预测，将预测结果与实际标签进行比较后修正冻结步态检测模型；
S5.提取待测样本的分阶段特征，输入至训练好的冻结步态检测模型，得到冻结步态检测结果。
2.如权利要求1所述的方法，其特征在于，步骤S1中，构建样本集之前，对传感器获取的原始数据进行数据清洗、滤波降噪和动静态分离；
所述滤波降噪，首先使用巴特沃斯低通滤波消除数据清洗处理后信号中的尖峰噪声，再使用滑动平均滤波消除随机噪声；
所述动静态分离，首先计算滤波降噪处理后信号的急动度，再根据急动度的大小，剔除单次试验前后的非运动态加速度数据，得到运动态加速度数据。
3.如权利要求2所述的方法，其特征在于，所述根据急动度的大小，剔除单次试验前后的非运动态加速度数据，得到运动态加速度数据，具体包括以下步骤：
I.将滤波降噪处理后加速度数据划分为M个长度为N的连续序列，M与N满足以下关系：
Figure FDA0002625816910000021
其中，Lacc为同一传感器所采集加速度经过滤波降噪处理后的序列长度；
II.计算每个序列的特征值T，计算公式如下：
Figure FDA0002625816910000022
其中，Jxi、Jyi、Jzi分别表示采样点i处x，y，z轴的急动度；
III.对于每个序列的特征值T，设置阈值Tsh，将T＞Tsh的第一个序列记为运动步态开始序列，T＞Tsh的最后一个序列记为运动步态结束序列。
4.如权利要求2或3所述的方法，其特征在于，对运动态加速度数据进行加窗处理，将每个窗口的数据序列
Figure FDA0002625816910000023
作为一个样本，采样点序号sw＝1+(w-1)×S，窗口序号w＝{1，…，m}，S为步长，m为划分窗口总数量，若窗口内正常步态所占比例大于0.5，则样本标签为正常步态，否则，样本标签为冻结步态。
5.如权利要求1至4任一项所述的方法，其特征在于，步骤S1中，构建样本集时，采用SMOTE过采样算法增加冻结步态样本的数量，保证正常步态样本数与冻结步态样本数之比为2∶1。
6.如权利要求1至5任一项所述的方法，其特征在于，运动加速度的运动分量计算过程如下：
(1)采用4阶巴特沃斯高通滤波器分别对降噪后每个维度的加速度滤波，得到高频分量{ahx，ahy，ahz}，截止频率为0.5Hz；
(2)采用4阶巴特沃斯低通滤波器分别对降噪后每个维度的加速度滤波，得到低频分量{alx，aly，alz}，截止频率为0.5Hz；
(3)分别计算高频分量的欧几里得范数HFEN和低频分量的欧几里得范数LFEN
Figure FDA0002625816910000031
Figure FDA0002625816910000032
(4)按照公式MAcc＝HFEN+(LFEN-1)计算运动分量。
7.如权利要求1至6任一项所述的方法，其特征在于，所述时域特征包括：均值、标准差、窗口内最大值、窗口内最小值、窗口内最大值与最小值的差值、窗口内超过均值点的个数、第三四分位数与第一四分位数的差距、窗口内数据平方的均值、XY两轴之间的相关性、XZ两轴之间的相关性、YZ两轴之间的相关性和信号幅值面积；所述频域特征包括：直流分量、功率谱密度的幅值统计特征、功率谱密度的形状统计特征和冻结指数FI。
8.如权利要求1至7任一项所述的方法，其特征在于，基于相关系数的特征选择步骤为：
(1)计算低维新步态特征集各特征间的皮尔逊相关系数，获得相关性矩阵；
Figure FDA0002625816910000033
其中，r为特征
Figure FDA0002625816910000034
与
Figure FDA0002625816910000035
间皮尔逊相关系数，特征
Figure FDA0002625816910000036
与
Figure FDA0002625816910000037
对应样本向量分别为
Figure FDA0002625816910000041
(2)设置相关性选择阈值k，k∈[0，1]；
(3)计算各特征的相关性总和，即对相关性矩阵按列求和；
(4)筛选出相关性矩阵中取值大于k的特征对，删除其中相关性总和较大的特征量；
(5)剩余特征组成新的特征集。
9.如权利要求1至7任一项所述的方法，其特征在于，基于最大信息系数法的特征选择步骤为：
(1)计算所有特征Xl与步态类别Y之间的MIC值，选择MIC最大值对应的特征为特征子集的第一个特征；
(2)选择J最大值对应的特征加入特征子集；
Figure FDA0002625816910000042
其中，Xl为候选特征，XF为已选择的特征集；
(3)重复步骤(2)直至所选特征个数达到预设值。
10.一种基于分阶段特征提取的冻结步态检测系统，其特征在于，包括：计算机可读存储介质和处理器；
所述计算机可读存储介质用于存储可执行指令；
所述处理器用于读取所述计算机可读存储介质中存储的可执行指令，执行权利要求1至9任一项所述的基于分阶段特征提取的冻结步态检测方法。
CN202010798417.2A 2020-08-10 2020-08-10 一种基于分阶段特征提取的冻结步态检测系统 Active CN112057080B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010798417.2A CN112057080B (zh) 2020-08-10 2020-08-10 一种基于分阶段特征提取的冻结步态检测系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010798417.2A CN112057080B (zh) 2020-08-10 2020-08-10 一种基于分阶段特征提取的冻结步态检测系统

Publications (2)

Publication Number Publication Date
CN112057080A true CN112057080A (zh) 2020-12-11
CN112057080B CN112057080B (zh) 2021-07-02

Family

ID=73662622

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010798417.2A Active CN112057080B (zh) 2020-08-10 2020-08-10 一种基于分阶段特征提取的冻结步态检测系统

Country Status (1)

Country Link
CN (1) CN112057080B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112729331A (zh) * 2020-12-21 2021-04-30 卡哈科技(深圳)有限公司 基于三轴加速度计的双带通滤波计步方法及电子设备
CN113408585A (zh) * 2021-05-21 2021-09-17 上海师范大学 一种基于人工智能的智能印章移动检测方法
CN113790722A (zh) * 2021-08-20 2021-12-14 北京自动化控制设备研究所 一种基于惯性数据时频域特征提取的行人步长建模方法
CN114287921A (zh) * 2021-12-23 2022-04-08 常州信息职业技术学院 一种步态双侧相似性分析方法、装置及系统
EP4062827A1 (en) * 2021-03-27 2022-09-28 Tata Consultancy Services Limited Patient invariant model for freezing of gait detection based on empirical wavelet decomposition
CN115171886A (zh) * 2022-07-25 2022-10-11 北京戴来科技有限公司 一种基于随机森林算法的冻结步态检测方法、装置、和存储介质

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009083032A1 (en) * 2007-12-28 2009-07-09 Fundacion Fatronik Method to detect physical blocking episodes on an individual activity
US20150099945A1 (en) * 2013-10-07 2015-04-09 Wahoo Fitness Llc Activity monitoring computing device and system
EP2387357B1 (en) * 2008-12-23 2016-04-20 Oxford Brookes University Gait monitor
CN106456058A (zh) * 2014-02-04 2017-02-22 特拉维夫医疗中心医疗研究、基础设施与健康服务基金会 用在使用身体固定感测器提供帕金森氏病的诊断或预后的方法及系统
WO2017060132A1 (en) * 2015-10-08 2017-04-13 Koninklijke Philips N.V. Treatment apparatus and method for treating a gait irregularity of a person
CN107928679A (zh) * 2017-12-07 2018-04-20 中国科学院宁波材料技术与工程研究所 帕金森患者冻结步态的诊断装置
CN108095729A (zh) * 2018-01-24 2018-06-01 深圳市臻络科技有限公司 一种冻结步态识别方法和装置
CN108731664A (zh) * 2018-05-18 2018-11-02 深圳清创新科技有限公司 机器人状态估计方法、装置、计算机设备和存储介质
CN108836344A (zh) * 2018-04-26 2018-11-20 深圳市臻络科技有限公司 步长步频估算方法和装置及步态检测仪
CN110638458A (zh) * 2019-08-26 2020-01-03 广东省人民医院(广东省医学科学院) 一种基于步态数据的康复训练效果评估方法、装置
CN110916984A (zh) * 2019-12-03 2020-03-27 上海交通大学医学院附属第九人民医院 一种预防冻结步态的穿戴设备及其实现方法
EP3683114A1 (en) * 2019-01-15 2020-07-22 Toyota Jidosha Kabushiki Kaisha Driving assistance apparatus

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009083032A1 (en) * 2007-12-28 2009-07-09 Fundacion Fatronik Method to detect physical blocking episodes on an individual activity
EP2387357B1 (en) * 2008-12-23 2016-04-20 Oxford Brookes University Gait monitor
US20150099945A1 (en) * 2013-10-07 2015-04-09 Wahoo Fitness Llc Activity monitoring computing device and system
CN106456058A (zh) * 2014-02-04 2017-02-22 特拉维夫医疗中心医疗研究、基础设施与健康服务基金会 用在使用身体固定感测器提供帕金森氏病的诊断或预后的方法及系统
WO2017060132A1 (en) * 2015-10-08 2017-04-13 Koninklijke Philips N.V. Treatment apparatus and method for treating a gait irregularity of a person
CN107928679A (zh) * 2017-12-07 2018-04-20 中国科学院宁波材料技术与工程研究所 帕金森患者冻结步态的诊断装置
CN108095729A (zh) * 2018-01-24 2018-06-01 深圳市臻络科技有限公司 一种冻结步态识别方法和装置
CN108836344A (zh) * 2018-04-26 2018-11-20 深圳市臻络科技有限公司 步长步频估算方法和装置及步态检测仪
CN108731664A (zh) * 2018-05-18 2018-11-02 深圳清创新科技有限公司 机器人状态估计方法、装置、计算机设备和存储介质
EP3683114A1 (en) * 2019-01-15 2020-07-22 Toyota Jidosha Kabushiki Kaisha Driving assistance apparatus
CN110638458A (zh) * 2019-08-26 2020-01-03 广东省人民医院(广东省医学科学院) 一种基于步态数据的康复训练效果评估方法、装置
CN110916984A (zh) * 2019-12-03 2020-03-27 上海交通大学医学院附属第九人民医院 一种预防冻结步态的穿戴设备及其实现方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
ANDEAS KUHNER,ET AL.: "Motion Biomarkers Showing Maximum Contrast Between Healthy Subjects and Parkinson"s Disease Patients Treated With Deep Brain Stimulatin of the Subthalamic Nucleus.A Pilot Study", 《FRONTIERS IN NEUROSCIENCE》 *
施翼: "基于运动感知的帕金森患者冻结步态", 《中国优秀硕士学位论文全文数据库医药卫生科技辑》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112729331A (zh) * 2020-12-21 2021-04-30 卡哈科技(深圳)有限公司 基于三轴加速度计的双带通滤波计步方法及电子设备
CN112729331B (zh) * 2020-12-21 2024-03-08 卡哈科技(深圳)有限公司 基于三轴加速度计的双带通滤波计步方法及电子设备
EP4062827A1 (en) * 2021-03-27 2022-09-28 Tata Consultancy Services Limited Patient invariant model for freezing of gait detection based on empirical wavelet decomposition
CN113408585A (zh) * 2021-05-21 2021-09-17 上海师范大学 一种基于人工智能的智能印章移动检测方法
CN113790722A (zh) * 2021-08-20 2021-12-14 北京自动化控制设备研究所 一种基于惯性数据时频域特征提取的行人步长建模方法
CN113790722B (zh) * 2021-08-20 2023-09-12 北京自动化控制设备研究所 一种基于惯性数据时频域特征提取的行人步长建模方法
CN114287921A (zh) * 2021-12-23 2022-04-08 常州信息职业技术学院 一种步态双侧相似性分析方法、装置及系统
CN115171886A (zh) * 2022-07-25 2022-10-11 北京戴来科技有限公司 一种基于随机森林算法的冻结步态检测方法、装置、和存储介质
CN115171886B (zh) * 2022-07-25 2023-03-24 北京戴来科技有限公司 一种基于随机森林算法的冻结步态检测方法、装置、设备和存储介质

Also Published As

Publication number Publication date
CN112057080B (zh) 2021-07-02

Similar Documents

Publication Publication Date Title
CN112057080B (zh) 一种基于分阶段特征提取的冻结步态检测系统
Zhang et al. Heart sound classification based on scaled spectrogram and partial least squares regression
CN105550659A (zh) 基于随机投影的实时心电分类方法
CN106228200A (zh) 一种不依赖于动作信息采集设备的动作识别方法
CN113095302B (zh) 用于心律失常分类的深度模型、利用该模型的方法及装置
Wang et al. A novel multi-scale dilated 3D CNN for epileptic seizure prediction
Zhao et al. An IoT-based wearable system using accelerometers and machine learning for fetal movement monitoring
CN107582077A (zh) 一种基于手机触摸行为的人体精神状态分析方法
CN108309304A (zh) 一种生成冻结步态智能监测系统的方法
CN106175754A (zh) 睡眠状态分析中清醒状态检测装置
Li et al. Patient-specific seizure prediction from electroencephalogram signal via multichannel feedback capsule network
CN108420406A (zh) 基于脉搏波睡眠分期的方法
CN111544005B (zh) 基于支持向量机的帕金森病人运动障碍量化及识别方法
Ye et al. Monitoring deterioration of knee osteoarthritis using vibration arthrography in daily activities
Ramović et al. Wavelet and Teager energy operator (TEO) for heart sound processing and identification
JP7489729B2 (ja) 転倒危険予防方法およびこのような方法を遂行する装置
CN113342166B (zh) 基于耳机的保护隐私的颈椎运动识别方法及系统
CN110507299A (zh) 一种心率信号检测装置及方法
CN109620241A (zh) 一种可穿戴设备及基于其的动作监测方法
Neili et al. Gammatonegram based pulmonary pathologies classification using convolutional neural networks
CN109044338A (zh) 房颤检测装置及存储介质
CN110477916B (zh) 一种表面肌电信号的运动起始时刻识别方法与装置
Assam et al. Activity recognition from sensors using dyadic wavelets and hidden markov model
CN113095379A (zh) 一种基于穿戴式六轴传感数据的人体运动状态识别方法
Şengür et al. Classification of amyotrophic lateral sclerosis and healthy electromyography signals based on transfer learning

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant