[go: up one dir, main page]

CN107157450A - Quantitative estimation method and system are carried out for the hand exercise ability to patient Parkinson - Google Patents

Quantitative estimation method and system are carried out for the hand exercise ability to patient Parkinson Download PDF

Info

Publication number
CN107157450A
CN107157450A CN201710463520.XA CN201710463520A CN107157450A CN 107157450 A CN107157450 A CN 107157450A CN 201710463520 A CN201710463520 A CN 201710463520A CN 107157450 A CN107157450 A CN 107157450A
Authority
CN
China
Prior art keywords
specified action
data set
classifier
mrow
training data
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
CN201710463520.XA
Other languages
Chinese (zh)
Other versions
CN107157450B (en
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.)
Institute of Computing Technology of CAS
Original Assignee
Institute of Computing Technology of CAS
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 Institute of Computing Technology of CAS filed Critical Institute of Computing Technology of CAS
Priority to CN201710463520.XA priority Critical patent/CN107157450B/en
Publication of CN107157450A publication Critical patent/CN107157450A/en
Application granted granted Critical
Publication of CN107157450B publication Critical patent/CN107157450B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/48Other medical applications
    • A61B5/4848Monitoring or testing the effects of treatment, e.g. of medication
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Measuring devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/11Measuring movement of the entire body or parts thereof, e.g. head or hand tremor or mobility of a limb
    • A61B5/1124Determining motor skills
    • A61B5/1125Grasping motions of hands
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/389Electromyography [EMG]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/40Detecting, measuring or recording for evaluating the nervous system
    • A61B5/4076Diagnosing or monitoring particular conditions of the nervous system
    • A61B5/4082Diagnosing or monitoring movement diseases, e.g. Parkinson, Huntington or Tourette
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7264Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • A61B5/7267Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/21Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/214Generating training patterns; Bootstrap methods, e.g. bagging or boosting
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/24Classification techniques

Landscapes

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

Abstract

本发明提供了用于对帕金森病人的手部运动能力进行量化评估方法和系统,利用可穿戴的肌电传感器捕获被检测者在执行指定动作时的手部表面肌电信号,基于训练好的与所述指定动作对应的运动能力分类器以从表面肌电信号中提取时域特征和频域特征以及与指定动作的完成情况相关的特征来对被检测者的手部运动能力进行量化评估。通过该方法和系统可以更客观准确地对帕金森病人的手部运动能力进行等级评估。

The present invention provides a method and system for quantitatively evaluating the hand movement ability of Parkinson's patients, using a wearable myoelectric sensor to capture the hand surface electromyography signal of the detected person when performing a specified action, based on the trained The motion ability classifier corresponding to the specified action extracts time-domain features and frequency-domain features from the surface electromyographic signal, as well as features related to the completion of the specified action, to quantitatively evaluate the hand motion ability of the detected person. The method and system can more objectively and accurately evaluate the hand movement ability of Parkinson's patients.

Description

用于对帕金森病人的手部运动能力进行量化评估方法和系统Method and system for quantitative evaluation of hand motor ability of Parkinson's patients

技术领域technical field

本发明涉及运动能力的定量评估,尤其涉及帕金森病人的手部运动能力的量化评估方法和系统。The invention relates to the quantitative evaluation of motor ability, in particular to a quantitative evaluation method and system for Parkinson's patient's hand motor ability.

背景技术Background technique

帕金森病(Parkinson’s Disease,PD)患者的运动症状主要表现为运动迟缓、静止性震颤和肌肉僵直等。目前对于PD患者的运动能力的检测和评估主要是依赖于国际运动协会提供的PD评定量表(UnifiedParkinson's Disease Rating Scale,UPDRS)来进行等级评定,例如,通常帕金森病人的手部运动能力被分为0-4级,0级相当于正常人的手部运动能力,级数越高,手部运动能力越差。然而,这种评定方式通常受到评分医师操作经验和评估时被检测者的状态和情绪的影响,因此其评估的结果还是不够客观准确。The motor symptoms of patients with Parkinson's Disease (PD) mainly include slowness of movement, resting tremor, and muscle stiffness. At present, the detection and evaluation of the motor ability of PD patients mainly relies on the PD Rating Scale (Unified Parkinson's Disease Rating Scale, UPDRS) provided by the International Sports Association for grade evaluation. For example, the hand motor ability of Parkinson patients is usually divided into It is graded from 0 to 4. Grade 0 is equivalent to the hand motor ability of normal people. The higher the grade, the worse the hand motor ability. However, this evaluation method is usually affected by the scoring physician's operating experience and the state and emotion of the person being tested during the evaluation, so the evaluation results are still not objective and accurate enough.

发明内容Contents of the invention

因此,本发明的目的在于克服上述现有技术的缺陷,提供一种利用手部的表面肌电信号对帕金森病人的手部运动能力的量化评估方法和系统。Therefore, the object of the present invention is to overcome the defects of the above-mentioned prior art, and provide a method and system for quantitatively evaluating the hand motor ability of a Parkinson's patient by using the surface electromyographic signal of the hand.

本发明的目的是通过以下技术方案实现的:The purpose of the present invention is achieved through the following technical solutions:

一方面,本发明提供了一种用于对帕金森病人的手部运动能力进行量化评估方法,包括:On the one hand, the present invention provides a method for quantitatively evaluating the hand motor ability of Parkinson's patients, comprising:

步骤1,经由放置在被检测者手部的肌电传感器采集被检测者执行指定动作时的表面肌电信号;Step 1, collect the surface electromyography signal when the subject performs a specified action via the electromyographic sensor placed on the subject's hand;

步骤2,利用预先训练好的与所述指定动作对应的帕金森病人运动能力分类器来根据所采集的表面肌电信号判断被检测者的手部运动能力所属的等级;Step 2, using the pre-trained Parkinson's patient motor classifier corresponding to the specified action to judge the grade of the hand motor ability of the detected person according to the collected surface electromyographic signals;

其中用于训练与指定动作对应的帕金森病人运动能力分类器的特征包括从表面肌电信号中提取时域特征和频域特征以及基于表面肌电信号获取的与所述指定动作的完成情况相关的特征。The features used to train the Parkinson's patient's motor ability classifier corresponding to the specified action include extracting time-domain features and frequency-domain features from the surface electromyography signal and obtaining information related to the completion of the specified action based on the surface electromyography signal. Characteristics.

上述方法中,与所述指定动作的完成情况相关的特征可包括下列中的至少一个:完成指定动作时的最大肌电信号、完成指定动作所用的时间。In the above method, the feature related to the completion of the specified action may include at least one of the following: the maximum electromyographic signal when the specified action is completed, and the time taken to complete the specified action.

上述方法中,所述指定动作可包括下列中的至少一个:握拳、对指。In the above method, the specified action may include at least one of the following: making a fist and pointing fingers.

上述方法中,还可包括训练与指定动作对应的帕金森病人运动能力分类器的步骤,其包括:In the above method, it may also include the step of training a Parkinson's patient motor ability classifier corresponding to the specified action, which includes:

a)接收来自多个不同程度的帕金森病人和正常人执行所述指定动作时由其佩戴的肌电传感器采集的多个表面肌电信号作为训练数据集;a) receiving a plurality of surface myoelectric signals collected by a myoelectric sensor worn by a plurality of different degrees of Parkinson's patients and normal people when performing the specified action as a training data set;

b)从各表面肌电信号中提取时域特征和频域特征以及基于表面肌电信号获取与所述指定动作的完成情况相关的特征;b) extracting time-domain features and frequency-domain features from each surface electromyography signal and obtaining features related to the completion of the specified action based on the surface electromyography signal;

c)基于各个特征对于训练数据集的信息增益选择用于训练所述运动能力分类器的特征,其中每个特征对训练数据集的信息增益为训练数据集的经验熵与在给定该特征的条件下训练数据集的经验条件熵之间的差;C) select the feature for training the athletic ability classifier based on the information gain of each feature for the training data set, wherein the information gain of each feature for the training data set is the experience entropy of the training data set and given this feature The difference between the empirical conditional entropy of the training dataset under condition ;

d)基于所选择的特征来训练所述运动能力分类器。d) Training the athletic ability classifier based on the selected features.

上述方法中,当所述指定动作为对指时,可采用序列最小优化模型作为所述运动能力分类器;当所述指定动作为握拳时,可采用J48分类器作为所述运动能力分类器。In the above method, when the specified action is pointing fingers, a sequential minimal optimization model can be used as the athletic ability classifier; when the specified action is fist clenching, a J48 classifier can be used as the athletic ability classifier.

上述方法中,训练数据集的经验熵表示对训练数据集进行分类的不确定性,可以如下公式进行计算:In the above method, the empirical entropy of the training data set represents the uncertainty of classifying the training data set, which can be calculated by the following formula:

其中D表示训练数据集,H(D)表示训练数据集D的经验熵,n表示训练数据集D共分为n类,pi表示把数据分为第i类的概率。Among them, D represents the training data set, H(D) represents the experience entropy of the training data set D, n represents the training data set D is divided into n categories, and p i represents the probability of dividing the data into the i-th category.

上述方法中,在给定某个特征的条件下训练数据集的经验条件熵表示在给定该特征的条件下对训练数据集进行分类的不确定性,可以如下公式来进行计算:In the above method, the empirical conditional entropy of the training data set under the condition of a given feature represents the uncertainty of classifying the training data set under the condition of a given feature, which can be calculated by the following formula:

其中A表示某个特征,H(D|A)为在给定特征A的条件下训练数据集D的经验熵,pi(D|A)表示在已知特征A的情况下把数据分为第i类的概率。Where A represents a certain feature, H(D|A) is the experience entropy of the training data set D under the condition of given feature A, p i (D|A) means that the data is divided into The probability of class i.

又一方面,本发明提供了一种用于对帕金森病人的手部运动能力进行量化评估系统,包括:In yet another aspect, the present invention provides a system for quantitatively evaluating the hand motor ability of Parkinson's patients, comprising:

采集装置,用于经由放置在被检测者手部的肌电传感器采集被检测者执行指定动作时的表面肌电信号;The acquisition device is used to collect the surface electromyography signal of the subject when the subject performs a specified action via the myoelectric sensor placed on the subject's hand;

检测装置,用于利用预先训练好的与所述指定动作对应的帕金森病人运动能力分类器来根据所采集的表面肌电信号判断被检测者的手部运动能力所属的等级;The detection device is used to use the pre-trained Parkinson's patient movement ability classifier corresponding to the specified action to judge the grade of the hand movement ability of the detected person according to the collected surface electromyographic signal;

其中用于训练与指定动作对应的帕金森病人运动能力分类器的特征包括从表面肌电信号中提取时域特征和频域特征以及基于表面肌电信号获取的与所述指定动作的完成情况相关的特征。The features used to train the Parkinson's patient's motor ability classifier corresponding to the specified action include extracting time-domain features and frequency-domain features from the surface electromyography signal and obtaining information related to the completion of the specified action based on the surface electromyography signal. Characteristics.

上述系统中,与所述指定动作的完成情况相关的特征可包括下列中的至少一个:完成指定动作时的最大肌电信号、完成指定动作所用的时间。In the above system, the features related to the completion of the specified action may include at least one of the following: the maximum electromyographic signal when the specified action is completed, and the time taken to complete the specified action.

上述系统中,所述指定动作可包括下列中的至少一个:握拳、对指。In the above system, the specified action may include at least one of the following: making a fist and pointing fingers.

上述系统还可包括用于训练与指定动作对应的帕金森病人运动能力分类器的训练装置,其被配置为:The system described above may also include a training device for training a Parkinson's patient motor capacity classifier corresponding to a specified action, which is configured to:

a)接收来自多个不同程度的帕金森病人和正常人执行所述指定动作时由其佩戴的肌电传感器采集的多个表面肌电信号作为训练数据集;a) receiving a plurality of surface myoelectric signals collected by a myoelectric sensor worn by a plurality of different degrees of Parkinson's patients and normal people when performing the specified action as a training data set;

b)从各表面肌电信号中提取时域特征和频域特征以及基于表面肌电信号获取与所述指定动作的完成情况相关的特征;b) extracting time-domain features and frequency-domain features from each surface electromyography signal and obtaining features related to the completion of the specified action based on the surface electromyography signal;

c)基于各个特征对于训练数据集的信息增益选择用于训练所述运动能力分类器的特征,其中每个特征对训练数据集的信息增益为训练数据集的经验熵与在给定该特征的条件下训练数据集的经验条件熵之间的差;C) select the feature for training the athletic ability classifier based on the information gain of each feature for the training data set, wherein the information gain of each feature for the training data set is the experience entropy of the training data set and given this feature The difference between the empirical conditional entropy of the training dataset under condition ;

d)基于所选择的特征来训练所述运动能力分类器。d) Training the athletic ability classifier based on the selected features.

上述系统中,当所述指定动作为对指时,可采用序列最小优化模型作为所述运动能力分类器;当所述指定动作为握拳时,可采用J48分类器作为所述运动能力分类器。In the above system, when the specified action is pointing fingers, a sequential minimal optimization model can be used as the athletic ability classifier; when the specified action is fist clenching, a J48 classifier can be used as the athletic ability classifier.

与现有技术相比,本发明的优点在于:Compared with the prior art, the present invention has the advantages of:

利用可穿戴的肌电传感器捕获被检测者手部表面肌电信号,以从表面肌电信号中提取时域特征和频域特征以及与指定动作的完成情况相关的特征来对被检测者的手部运动能力进行量化评估,可以更客观准确地对帕金森病人的手部运动能力进行等级评估。Use the wearable electromyographic sensor to capture the surface electromyographic signal of the subject's hand to extract time-domain features and frequency-domain features from the surface electromyographic signal, as well as features related to the completion of the specified action to analyze the subject's hand. Quantitative assessment of the hand movement ability of Parkinson's patients can be performed more objectively and accurately.

附图说明Description of drawings

以下参照附图对本发明实施例作进一步说明,其中:Embodiments of the present invention will be further described below with reference to the accompanying drawings, wherein:

图1为根据本发明实施例的用于对帕金森病人的手部运动能力进行量化评估方法的流程示意图;FIG. 1 is a schematic flowchart of a method for quantitatively evaluating the hand motor ability of a Parkinson's patient according to an embodiment of the present invention;

图2(a)为与“对指”动作对应的各分类器的识别精度比较示意图;Figure 2(a) is a schematic diagram of the comparison of the recognition accuracy of each classifier corresponding to the "pointing finger" action;

图2(b)为与“握拳”动作对应的各分类器的识别精度比较示意图;Figure 2(b) is a schematic diagram of the comparison of the recognition accuracy of each classifier corresponding to the "clench fist" action;

图3(a)为对于“对指”动作利用表面肌电信号的信号识别特征训练的SMO分类器识别性能示意图;Fig. 3 (a) is the SMO classifier recognition performance schematic diagram that utilizes the signal recognition feature training of surface electromyographic signal for " pointing to " action;

图3(b)为对于“握拳”动作利用表面肌电信号的信号识别特征训练的J48分类器识别性能示意图;Figure 3 (b) is a schematic diagram of the J48 classifier recognition performance training using the signal recognition feature of the surface electromyography signal for the "clench fist" action;

图4(a)为采用根据本发明的方法训练的与“对指”动作对应的分类器识别性能示意图;Fig. 4 (a) is a schematic diagram of classifier recognition performance corresponding to the "finger" action trained by the method according to the present invention;

图4(b)为采用根据本发明的方法训练的与“握拳”动作对应的分类器识别性能示意图;Fig. 4 (b) is a schematic diagram of classifier recognition performance corresponding to the action of "clenching a fist" trained according to the method of the present invention;

图5(a)为采用不同特征训练的与“对指”动作对应的分类器的识别结果对比示意图;Figure 5(a) is a schematic diagram of the comparison of the recognition results of the classifiers corresponding to the "finger-to-finger" action trained with different features;

图5(b)为采用不同特征训练的与“握拳”动作对应的分类器的识别结果对比示意图。Figure 5(b) is a schematic diagram of the comparison of the recognition results of the classifier corresponding to the "fist" action trained with different features.

具体实施方式detailed description

为了使本发明的目的,技术方案及优点更加清楚明白,以下结合附图通过具体实施例对本发明进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。In order to make the object, technical solution and advantages of the present invention clearer, the present invention will be further described in detail below through specific embodiments in conjunction with the accompanying drawings. It should be understood that the specific embodiments described here are only used to explain the present invention, not to limit the present invention.

“对指”和“握拳”手势是用来评估帕金森病人的手部运动能力的两种常用的国际标准指标。在下文中,将以这两种指定动作为例来结合具体实施例说明如何评估帕金森病人的手部运动能力。但应理解,本发明的实施例对于具体的运动动作并不进行限制,本发明的原理也可以应用于其他动作或其他部位的运动能力评估。The "pointing" and "fist" gestures are two commonly used international standard indicators for assessing the hand motor ability of Parkinson's patients. In the following, the two specified actions will be taken as examples to describe how to evaluate the hand motor ability of Parkinson's patients in combination with specific embodiments. However, it should be understood that the embodiments of the present invention are not limited to specific sports actions, and the principles of the present invention may also be applied to other actions or the assessment of the athletic ability of other parts.

图1给出了根据本发明的一个实施例的用于对帕金森病人的手部运动能力进行量化评估方法,其包括经由放置在被检测者手部的相关部位的肌电传感器采集被检测者执行指定动作时的表面肌电信号(步骤1);利用预先训练好的与所述指定动作对应的帕金森病人运动能力分类器来根据所采集的表面肌电信号判断被检测者的手部运动能力所属的等级(步骤2)。其中与指定动作对应的帕金森病人运动能力分类器是在大量的不同程度的帕金森病人和一些正常人执行该指定动作时提取的样本数据的基础上进行训练的。Figure 1 shows a method for quantitatively evaluating the hand motor ability of a Parkinson's patient according to an embodiment of the present invention, which includes collecting the data of the subject via a myoelectric sensor placed on a relevant part of the subject's hand. The surface electromyographic signal (step 1) when performing the specified action; the Parkinson's patient's motor ability classifier corresponding to the specified action is used to judge the hand movement of the detected person according to the collected surface electromyographic signal The level to which the ability belongs (step 2). The motor capacity classifier of Parkinson's patients corresponding to the specified action is trained on the basis of sample data extracted from a large number of Parkinson's patients of different degrees and some normal people when performing the specified action.

更具体地,在步骤1,将肌电传感器放置在被检测者的指定位置,例如手臂等相关部位,然后让被检测者执行指定动作,例如握拳或对指。这里,可以采用能够采集到人体表面肌电信号的任何类型的肌电传感器。优选地,可以采用可穿戴肌电仪,例如在下文中采用由加拿大创业公司ThalmicLabs出品的臂环MYO来检测手部运动时表面肌电信号,其是一种手势控制臂环,可以通过蓝牙与其他电子产品进行无线连接和数据传输。被检测者在做握拳和对指动作时,主要是小臂上方的肌肉进行活动,产生相应的肌电信号,所以可以将MYO放置在被检测者的小臂上方。肌电传感器将采集到的表面肌电信号通过有线或无线的传输方式发送给相应的信号分析和处理装置,例如,诸如平板或手机之类的移动终端、台式电脑以及能够处理表面肌电信号的任何其他计算装置。More specifically, in step 1, the myoelectric sensor is placed on a designated position of the subject, such as an arm and other relevant parts, and then the subject is asked to perform a designated action, such as making a fist or pointing fingers. Here, any type of myoelectric sensor capable of collecting human surface myoelectric signals can be used. Preferably, a wearable electromyographic instrument can be used, for example, the armband MYO produced by the Canadian startup company ThalmicLabs is used below to detect the surface electromyographic signal during hand movement. It is a gesture control armband that can communicate with other Electronics make wireless connections and data transfers. When the subject is making a fist and pointing fingers, the muscles above the forearm are mainly active to generate corresponding myoelectric signals, so MYO can be placed on the forearm of the subject. The electromyographic sensor sends the collected surface electromyographic signals to corresponding signal analysis and processing devices through wired or wireless transmission, for example, mobile terminals such as tablets or mobile phones, desktop computers, and devices capable of processing surface electromyographic signals. any other computing device.

优选地,负责处理表面肌电信号的计算装置还可以对接收到的表面肌电信号进行一定程度的预处理,例如过滤掉噪声或者进行降维处理,例如对来自多个肌电传感器的信号进行合成,以尽量简化计算复杂度,减少对计算资源的消耗。以MYO为例,其由八个传感器芯片组成,可以采集到患者小臂一周共八维的肌电信号,分别用a1,a2,a3,a4,a5,a6,a7,a8来表示。当被检测者做“握拳”和“对指”时,小臂表面的肌肉会处于激活状态。例如当患者握拳时,小臂一周的肌肉同时处于紧绷状态,也就是激活状态,导致肌电仪上的八个肌电传感器的肌电信号同时增大。当患者手臂张开时,小臂上一周的肌肉同时处于放松状态,或称为非激活状态,导致八片肌电传感器的肌电信号同时减小或者消失。所采集到的八维肌电信号具有的相关性,由此可以采用合成加速度的方法来将所采集的肌电信号进行合成,例如公式(1)从而得到合成的肌电信号sample。使用合成肌电信号在减少计算复杂度的同时,还可以让佩戴者忽略佩戴肌电仪的方向,可以随意佩戴这八个肌电贴片组成的肌电仪。应理解,以上预处理进行是优选的方式而非对此进行任何限制,使用单个肌电传感器采集的表面肌电信号也可以实现手部运动能力的量化评估。在下文中提到的表面肌电信号时不再特指是单个肌电传感器采集的肌电信号还是合成后的肌电信号,因为二者的处理方式是类似的。Preferably, the computing device responsible for processing the surface electromyography signals can also perform a certain degree of preprocessing on the received surface electromyography signals, such as filtering out noise or performing dimensionality reduction processing, such as performing processing on signals from multiple electromyography sensors. Synthesis to simplify the computational complexity as much as possible and reduce the consumption of computing resources. Taking MYO as an example, it consists of eight sensor chips, which can collect eight - dimensional myoelectric signals of the patient's forearm in one week . , a 8 to represent. When the subject makes a "fist" and "fingers", the muscles on the surface of the forearm will be activated. For example, when the patient makes a fist, the muscles around the forearm are in a tense state at the same time, that is, the activation state, which causes the myoelectric signals of the eight myoelectric sensors on the electromyography instrument to increase simultaneously. When the patient's arm is opened, the muscles of the upper forearm are in a relaxed state at the same time, or called an inactive state, resulting in the reduction or disappearance of the myoelectric signals of the eight myoelectric sensors. The collected eight-dimensional myoelectric signal has a correlation, so the method of synthesizing acceleration can be used to synthesize the collected myoelectric signal, such as formula (1) Thus, the synthesized EMG signal sample is obtained. While using synthetic EMG signals to reduce the computational complexity, it also allows the wearer to ignore the direction in which the EMG is worn, and can wear the EMG composed of eight EMG patches at will. It should be understood that the above preprocessing is a preferred method without any limitation, and the quantitative evaluation of the hand movement ability can also be realized by using the surface electromyographic signal collected by a single electromyographic sensor. The surface electromyographic signal mentioned below no longer specifically refers to the electromyographic signal collected by a single electromyographic sensor or the synthesized electromyographic signal, because the processing methods of the two are similar.

继续参考图1,在步骤2,在获得了执行指定动作时的表面肌电信号之后,从其中提取与该指定动作对应的帕金森病人运动能力分类器所需要的特征,接着将所提取的特征提供给预先训练好的分类器来判断被检测者的手部运动能力所属的等级。对于每个指定动作,都有预先训练好的与其相对应的帕金森病人运动能力分类器,为下文描述方法,将与握拳动作对应的分类器简称为握拳分类器,以及将与对指动作对应的分类器简称为对指分类器。这两个分类器的训练方式类似,因此,下文仅以握拳分类器为例来对分类器的训练方法进行举例说明。Continuing to refer to Fig. 1, in step 2, after obtaining the surface electromyographic signal when performing the specified action, extract the features required by the Parkinson's patient's motor ability classifier corresponding to the specified action, and then extract the features It is provided to a pre-trained classifier to judge the class to which the detected person's hand motor ability belongs. For each specified action, there is a pre-trained Parkinson's patient's motor ability classifier corresponding to it. For the method described below, the classifier corresponding to the fisting action will be referred to as the fisting classifier for short, and the classifier corresponding to the fingering action will be The classifier of is referred to as the index classifier. The training methods of these two classifiers are similar, so the following only uses the fist classifier as an example to illustrate the training method of the classifier.

在训练时首先让大量的不同程度的帕金森病人和一些正常人佩戴肌电传感器并执行握拳动作,从而采集到大量的表面肌电信号作为基本的训练数据集。在训练分类器时关键是确定采用数据的哪些特征来进行训练,这些特征会直接影响分类结果的准确性。通常,可以对表面肌电信号进行时域和频域分析,从中提取多个时域特征和频域特征来训练分类器。这些时域和频域特征反映了肌电信号波形本身的特点。在本发明的优选实施例中,除了使用时域特征和频域特征之外,还采用了与指定动作的完成情况相关的特征。这是考虑到不同程度的PD患者在完成某个指定动作时所用的时间、力量等有较大差异,因此,与指定动作的完成情况相关的特征也可以对于手部运动能力有较好的区分性。为了提高分类器本身的准确性,在本发明的实施例中通过基于信息增益的方式来从上述的多个特征中选择最优的一些特征来训练分类器。训练分类器的过程主要包括:During training, a large number of Parkinson's patients and some normal people of different degrees are firstly allowed to wear myoelectric sensors and perform fisting movements, so as to collect a large number of surface electromyography signals as a basic training data set. When training a classifier, the key is to determine which features of the data are used for training, and these features will directly affect the accuracy of the classification results. Usually, time-domain and frequency-domain analysis can be performed on surface EMG signals, from which multiple time-domain features and frequency-domain features can be extracted to train a classifier. These time-domain and frequency-domain features reflect the characteristics of the EMG waveform itself. In a preferred embodiment of the invention, in addition to using time-domain and frequency-domain features, features related to the completion of specified actions are also used. This is because the time and strength used by PD patients with different degrees to complete a specified action are quite different. Therefore, the characteristics related to the completion of the specified action can also better distinguish the hand motor ability. sex. In order to improve the accuracy of the classifier itself, in the embodiment of the present invention, some optimal features are selected from the above multiple features based on information gain to train the classifier. The process of training a classifier mainly includes:

(1)从表面肌电信号中提取时域特征和频域特征(1) Extract time domain features and frequency domain features from surface EMG signals

表面肌电信号(surface electromyogram signal,SEMS)是从人体皮肤表面通过传感器或电极记录神经肌肉活动时发放的生物电信号,属于无创性。表面肌电信号的幅度具有随机性,通常可表示为准高斯分布函数。通常数字化分析是处理表面肌电信号的主要手段,包括时域分析和频域分析。时域分析是将表面肌电信号看作时间的函数,通过分析得到肌电信号的某些统计特征,例如,常见的时域特征包括肌电信号在时域上的均值、一个窗口内的最大值、方差、标准差、众数、过零次数、均方值、过均值率等。频域分析是通过傅里叶变换将时域信号转换为频域信号,常用频域特征包括峰值、频率分量、直流分量、能量、形状均值、形状标准差、形状偏度、形状峰度、幅度均值、幅度标准差、幅度偏度、幅度峰度、中值频率、均值频率、频率范围、最高波峰频率、最高波峰幅值等等。上述时域和频域特征也可以统称为信号识别特征。Surface electromyogram signal (SEMS) is a bioelectrical signal emitted when neuromuscular activity is recorded from the surface of the human skin through sensors or electrodes, and is noninvasive. The amplitude of the surface electromyography signal is random and can usually be expressed as a quasi-Gaussian distribution function. Usually digital analysis is the main means of processing surface electromyographic signals, including time domain analysis and frequency domain analysis. Time-domain analysis regards the surface EMG signal as a function of time, and obtains some statistical characteristics of the EMG signal through analysis. For example, common time-domain features include the average value of the EMG signal in the time domain, the maximum value, variance, standard deviation, mode, zero-crossing times, mean square value, mean-crossing rate, etc. Frequency domain analysis is to convert time domain signals into frequency domain signals through Fourier transform. Common frequency domain features include peak value, frequency component, DC component, energy, shape mean, shape standard deviation, shape skewness, shape kurtosis, amplitude Mean, amplitude standard deviation, amplitude skewness, amplitude kurtosis, median frequency, mean frequency, frequency range, highest peak frequency, highest peak amplitude, and more. The above time-domain and frequency-domain features may also be collectively referred to as signal identification features.

(2)基于表面肌电信号获取与指定动作的完成情况相关的特征(2) Obtain features related to the completion of the specified action based on the surface electromyographic signal

与指定动作的完成情况相关的特征(可简称为动作特征)包括完成指定动作时的最大肌电信号、完成指定动作所用的时间等。随着PD患者病情的加重,帕金森病运动症状会出现暂停和迟缓,这将导致肌电信号强度的下降。另外,当PD运动症状相对严重时,完成动作所用时间会随之变长。因此,与指定动作的完成情况相关的特征也能有效地区分不同程度的PD患者的手部运动能力。Features related to the completion of the specified action (may be referred to as action features for short) include the maximum electromyographic signal when the specified action is completed, the time taken to complete the specified action, and the like. With the aggravation of PD patients, motor symptoms of Parkinson's disease will pause and slow down, which will lead to a decrease in the strength of EMG. In addition, when the motor symptoms of PD are relatively severe, the time taken to complete the action will be longer. Therefore, features related to the completion of specified actions can also effectively distinguish the hand motor ability of different degrees of PD patients.

(3)基于信息增益从上述多个信号识别特征和动作特征中选择有效特征在来训练分类器。(3) Select effective features from the above-mentioned multiple signal recognition features and action features based on information gain to train a classifier.

首先,计算各个特征的信息增益。信息增益表示得知特征X的信息而使得类Y的信息的不确定性减少程度。特征A对训练数据集D的信息增益记为g(D,A),其为训练数据集合D的经验熵或信息熵H(D)与在特征A给定条件下训练数据集合D的经验条件熵H(D|A)之间的差,即:First, the information gain of each feature is calculated. Information gain represents the degree to which the uncertainty of class Y information is reduced by knowing the information of feature X. The information gain of feature A to training data set D is denoted as g(D,A), which is the empirical entropy or information entropy H(D) of training data set D and the empirical condition of training data set D under the given condition of feature A The difference between the entropy H(D|A), namely:

g(D,A)=H(D)-H(D|A)g(D,A)=H(D)-H(D|A)

其中经验熵H(D)表示对数据集D进行分类的不确定性,而经验条件熵H(D|A)表示在特征A给定的条件下对数据集D进行分类的不确定性,二者之间的差即为信息增益,表示由于特征A而使得对数据集D的分类的不确定性减少的程度,显然对于训练数据集D而言,信息增益依赖于特征,不同的特征往往具有不同的信息增益,信息增益大的特征具有更强的分类能力。经验熵H(D)例如可以通过下面的公式来计算:Among them, the empirical entropy H(D) represents the uncertainty of classifying the data set D, and the empirical conditional entropy H(D|A) represents the uncertainty of classifying the data set D under the given condition of feature A, two The difference between the two is the information gain, which indicates the degree of reduction of the uncertainty of the classification of the data set D due to the feature A. Obviously, for the training data set D, the information gain depends on the feature, and different features often have Different information gain, features with large information gain have stronger classification ability. The empirical entropy H(D) can be calculated, for example, by the following formula:

其中n表示训练数据集D共分为n类,pi表示把数据分为第i类的概率,由上式所知熵H越大,分类的不确定性就越大。经验条件熵H(D|A)例如可以通过下面的公式来计算:Among them, n indicates that the training data set D is divided into n categories, p i indicates the probability of dividing the data into the i-th category, and the greater the entropy H is known from the above formula, the greater the uncertainty of the classification. The empirical conditional entropy H(D|A) can be calculated by the following formula, for example:

其中n表示训练数据集D共分为n类,pi(D|A)表示在已知特征A的情况下把数据分为第i类的概率,H(D|A)表征在特征A给定的条件下对数据集D进行分类的不确定性程度。Among them, n indicates that the training data set D is divided into n categories, p i (D|A) indicates the probability of classifying the data into the i-th category when the feature A is known, and H(D|A) represents that the feature A is given The degree of uncertainty in classifying the data set D under certain conditions.

在根据上述公式可以计算PD运动症状每个特征的信息增益之后,比较各个特征的信息增益,选择那些信息增益较大的特征作为用于训练分类器的,从训练数据集中的表面肌电信号中提取这些特征作为样本数据来训练分类器。After the information gain of each feature of PD motor symptoms can be calculated according to the above formula, the information gain of each feature is compared, and those features with larger information gain are selected as the ones used to train the classifier, from the surface electromyographic signals in the training data set These features are extracted as sample data to train the classifier.

为了进一步说明根据上述实施例的方法的效果,发明人还进行了下列试验。In order to further illustrate the effects of the methods according to the above embodiments, the inventors also conducted the following experiments.

首先,仅使用从表面肌电信号中提取的时域和频域特征来训练分类器。例如分别使用支持向量机SVM(Support Vector Machine)、序列最小优化SMO(Sequential minimaloptimization)模型、诸如J48的决策树模型和随机森林RF(RandomForest)模型对“握拳”和“对指”动作的训练数据进行处理,实验结果如图2(a)和2(b)所示。对于“握拳”动作,使用J48分类器的识别精度最高,可达到86.2%;对于“对指”动作,利用SMO分类器的识别精度最高,可达到63.7%。其中,对于“对指”动作利用SMO分类器识别其所属等级(即0-4级)的精度(precision)、召回率(Recall)和F-score的结果如图3(a)所示,对于每个等级,这三个指标在图中从左到右依次排列。对“握拳”动作,利用J48分类器识别其所属等级(即0-4级)的精度、召回率和F-score的结果如图3(b)所示。通过观察肌电信号强度,发现对指的肌电信号较握拳的肌电信号弱,这是因为对指的肌肉激活程度较握拳的肌肉激活程度较弱,而肌电信号相对弱时较难进行分类,所以对指的肌电信号的区分较握拳时不明显,识别时对指的精度没有握拳的高。当分类等级为0和1时,错误率较高,也就是说正常人的对指等级和动作轻微变慢的病人对指等级不易区分,这是由于正常人在大量的工作劳动后,也会产生动作轻微变慢的症状。当分类等级为2和3时,精度较低,如图3(a)所示。这是因为患者在做对指测试时,等级为2和3的患者的肌肉激活程度较弱,造成肌电信号强度偏小,而强度不大的肌电信号发生细微变化时,容易造成分类混淆的结果。以上这些因素是造成对指分类错误率高的主要原因。First, a classifier is trained using only time- and frequency-domain features extracted from surface EMG signals. For example, use support vector machine SVM (Support Vector Machine), sequential minimal optimization SMO (Sequential minimal optimization) model, decision tree model such as J48 and random forest RF (RandomForest) model to train data of "clenching fist" and "fingering" actions After processing, the experimental results are shown in Figure 2(a) and 2(b). For the "clench fist" action, the recognition accuracy using the J48 classifier is the highest, which can reach 86.2%; for the "finger-to-finger" action, the recognition accuracy using the SMO classifier is the highest, which can reach 63.7%. Among them, the results of precision (precision), recall rate (Recall) and F-score of the level (ie, 0-4 level) identified by the SMO classifier for the "finger-to-finger" action are shown in Figure 3(a). For each level, the three indicators are arranged in order from left to right in the figure. For the "clench fist" action, the results of using the J48 classifier to identify its class (ie, 0-4 class) of precision, recall and F-score are shown in Figure 3(b). By observing the strength of the EMG signal, it is found that the EMG signal of the finger is weaker than that of making a fist. This is because the muscle activation of the finger is weaker than the muscle activation of the fist, and it is difficult to carry out when the EMG signal is relatively weak. Classification, so the distinction of the electromyographic signals of the fingers is less obvious than that of making a fist, and the accuracy of the recognition of the fingers is not as high as that of making a fist. When the classification level is 0 and 1, the error rate is higher, that is to say, it is not easy to distinguish the finger-to-finger grade of a normal person from that of a patient whose movements are slightly slowed down. A slight slowing of movement occurs. When the classification level is 2 and 3, the accuracy is lower, as shown in Fig. 3(a). This is because when patients do finger alignment tests, patients with grades 2 and 3 have weaker muscle activation, resulting in a small strength of the EMG signal, and subtle changes in the EMG signal with low strength are likely to cause classification confusion the result of. These factors are the main reasons for the high error rate of finger classification.

接着,对于“握拳”动作,使用J48分类器,对于“对指”动作,使用SMO分类器,但是仅使用完成指定动作时的最大肌电信号、完成指定动作所用的时间这些动作特征来训练分类器。当使用J48分类器时,握拳分类器的分类精度达到70.6%,当使用SMO分类器时,对指分类器的识别精度达到72.5%。Next, for the "clench fist" action, use the J48 classifier, and for the "pointing finger" action, use the SMO classifier, but only use the action characteristics of the maximum EMG signal when completing the specified action and the time it takes to complete the specified action to train the classification device. When using the J48 classifier, the classification accuracy of the fist classifier reaches 70.6%, and when using the SMO classifier, the recognition accuracy of the finger classifier reaches 72.5%.

最后,通过上述计算信息增益的方式从上述信号识别特征和动作特征中选择前N个有效特征来训练分类器,采用与上文相同的分类器模型。如当使用J48分类器时,握拳分类器的精度达到90.8%。当使用SMO分类器时,对指分类器的精度达到82.3%。这样的指定动作分类器对于手部能力0-4等级的识别精度、召回率和F-score如图4(a)-(b)所示,其中对于每个等级,这三个指标在图中从左到右依次排列。图5(a)-(b)所示的分别为用信号识别特征、动作特征和融合特征的实验结果对比图。如图5所示,与仅使用信号识别特征或者仅适用于动作特征相比,采用信息增益的方式选择有效特征(在图中称为融合特征)后的分类精度得到了明显提高。Finally, the first N effective features are selected from the above-mentioned signal recognition features and action features by the above-mentioned method of calculating information gain to train a classifier, and the same classifier model as above is used. For example, when using the J48 classifier, the accuracy of the fist classifier reaches 90.8%. When using the SMO classifier, the finger classifier achieves an accuracy of 82.3%. The recognition accuracy, recall rate and F-score of such a designated action classifier for grades 0-4 of hand ability are shown in Fig. 4(a)-(b), where for each grade, these three indicators are shown in Arranged from left to right. Figure 5(a)-(b) shows the comparison of experimental results using signal recognition features, action features and fusion features, respectively. As shown in Figure 5, compared with only using signal recognition features or only applying to action features, the classification accuracy after using information gain to select effective features (called fusion features in the figure) has been significantly improved.

虽然本发明已经通过优选实施例进行了描述,然而本发明并非局限于这里所描述的实施例,在不脱离本发明范围的情况下还包括所做出的各种改变以及变化。Although the present invention has been described in terms of preferred embodiments, the present invention is not limited to the embodiments described herein, and various changes and changes are included without departing from the scope of the present invention.

Claims (10)

1.一种用于对帕金森病人的手部运动能力进行量化评估方法,所述方法包括:1. A method for quantifying the hand motor ability of a Parkinson's patient, said method comprising: 步骤1,经由放置在被检测者手部的肌电传感器采集被检测者执行指定动作时的表面肌电信号;Step 1, collect the surface electromyography signal when the subject performs a specified action via the electromyographic sensor placed on the subject's hand; 步骤2,利用预先训练好的与所述指定动作对应的帕金森病人运动能力分类器来根据所采集的表面肌电信号判断被检测者的手部运动能力所属的等级;Step 2, using the pre-trained Parkinson's patient motor classifier corresponding to the specified action to judge the grade of the hand motor ability of the detected person according to the collected surface electromyographic signals; 其中用于训练与指定动作对应的帕金森病人运动能力分类器的特征包括从表面肌电信号中提取时域特征和频域特征以及基于表面肌电信号获取的与所述指定动作的完成情况相关的特征。The features used to train the Parkinson's patient's motor ability classifier corresponding to the specified action include extracting time-domain features and frequency-domain features from the surface electromyography signal and obtaining information related to the completion of the specified action based on the surface electromyography signal. Characteristics. 2.根据权利要求1所述的方法,其中与所述指定动作的完成情况相关的特征包括下列中的至少一个:完成指定动作时的最大肌电信号、完成指定动作所用的时间。2. The method according to claim 1, wherein the features related to the completion of the specified action include at least one of the following: the maximum electromyographic signal when the specified action is completed, and the time taken to complete the specified action. 3.根据权利要求1所述的方法,其中所述指定动作包括下列中的至少一个:握拳、对指。3. The method according to claim 1, wherein the specified action comprises at least one of the following: making a fist, pointing fingers. 4.根据权利要求1-3中任一项所述的方法,还包括训练与指定动作对应的帕金森病人运动能力分类器的步骤,其包括:4. The method according to any one of claims 1-3, further comprising the step of training a Parkinson's patient's motor capacity classifier corresponding to specified actions, comprising: a)接收来自多个不同程度的帕金森病人和正常人执行所述指定动作时由其佩戴的肌电传感器采集的多个表面肌电信号作为训练数据集;a) receiving a plurality of surface myoelectric signals collected by a myoelectric sensor worn by a plurality of different degrees of Parkinson's patients and normal people when performing the specified action as a training data set; b)从各表面肌电信号中提取时域特征和频域特征以及基于表面肌电信号获取与所述指定动作的完成情况相关的特征;b) extracting time-domain features and frequency-domain features from each surface electromyography signal and obtaining features related to the completion of the specified action based on the surface electromyography signal; c)基于各个特征对于训练数据集的信息增益选择用于训练所述运动能力分类器的特征,其中每个特征对训练数据集的信息增益为训练数据集的经验熵与在给定该特征的条件下训练数据集的经验条件熵之间的差;C) select the feature for training the athletic ability classifier based on the information gain of each feature for the training data set, wherein the information gain of each feature for the training data set is the experience entropy of the training data set and given this feature The difference between the empirical conditional entropy of the training dataset under condition ; d)基于所选择的特征来训练所述运动能力分类器。d) Training the athletic ability classifier based on the selected features. 5.根据权利要求4所述的方法,其中当所述指定动作为对指时,采用序列最小优化模型作为所述运动能力分类器;当所述指定动作为握拳时,采用J48分类器作为所述运动能力分类器。5. The method according to claim 4, wherein when the specified action is pointing fingers, adopting the sequence minimum optimization model as the athletic ability classifier; when the specified action is making a fist, adopting the J48 classifier as the Ability classifier described above. 6.根据权利要求4所述的方法,其中,训练数据集的经验熵表示对训练数据集进行分类的不确定性,以如下公式进行计算:6. The method according to claim 4, wherein the empirical entropy of the training data set represents the uncertainty of classifying the training data set, and is calculated with the following formula: <mrow> <mi>H</mi> <mrow> <mo>(</mo> <mi>D</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msub> <mi>p</mi> <mi>i</mi> </msub> <mi>l</mi> <mi>o</mi> <mi>g</mi> <mi> </mi> <msub> <mi>p</mi> <mi>i</mi> </msub> </mrow> <mrow> <mi>H</mi> <mrow> <mo>(</mo> <mi>D</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msub> <mi>p</mi> <mi>i</mi> </msub> <mi>l</mi> <mi>o</mi> <mi>g</mi> <mi> </mi> <msub> <mi>p</mi> <mi>i</mi> </msub> </mrow> 其中D表示训练数据集,H(D)表示训练数据集D的经验熵,n表示训练数据集D共分为n类,pi表示把数据分为第i类的概率。Among them, D represents the training data set, H(D) represents the experience entropy of the training data set D, n represents the training data set D is divided into n categories, and p i represents the probability of dividing the data into the i-th category. 7.根据权利要求6所述的方法,其中在给定某个特征的条件下训练数据集的经验条件熵表示在给定该特征的条件下对训练数据集进行分类的不确定性,以如下公式来进行计算:7. The method according to claim 6, wherein the empirical conditional entropy of the training data set given a certain feature represents the uncertainty of classifying the training data set given the feature, as follows formula to calculate: <mrow> <mi>H</mi> <mrow> <mo>(</mo> <mi>D</mi> <mo>|</mo> <mi>A</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <msubsup> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </msubsup> <msub> <mi>p</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>D</mi> <mo>|</mo> <mi>A</mi> <mo>)</mo> </mrow> <msub> <mi>logp</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>D</mi> <mo>|</mo> <mi>A</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>H</mi> <mrow> <mo>(</mo> <mi>D</mi> <mo>|</mo> <mi>A</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <msubsup> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </msubsup> <msub> <mi>p</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>D</mi> <mo>|</mo> <mi>A</mi> <mo>)</mo> </mrow> <msub> <mi>logp</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>D</mi> <mo>|</mo> <mi>A</mi> <mo>)</mo> </mrow> </mrow> 其中A表示某个特征,H(D|A)为在给定特征A的条件下训练数据集D的经验熵,pi(D|A)表示在已知特征A的情况下把数据分为第i类的概率。Where A represents a certain feature, H(D|A) is the experience entropy of the training data set D under the condition of given feature A, p i (D|A) means that the data is divided into The probability of class i. 8.一种用于对帕金森病人的手部运动能力进行量化评估系统,所述系统包括:8. A system for quantitative evaluation of the hand motor ability of Parkinson's patients, said system comprising: 采集装置,用于经由放置在被检测者手部的肌电传感器采集被检测者执行指定动作时的表面肌电信号;The acquisition device is used to collect the surface electromyography signal of the subject when the subject performs a specified action via the myoelectric sensor placed on the subject's hand; 检测装置,用于利用预先训练好的与所述指定动作对应的帕金森病人运动能力分类器来根据所采集的表面肌电信号判断被检测者的手部运动能力所属的等级;The detection device is used to use the pre-trained Parkinson's patient movement ability classifier corresponding to the specified action to judge the grade of the hand movement ability of the detected person according to the collected surface electromyographic signal; 其中用于训练与指定动作对应的帕金森病人运动能力分类器的特征包括从表面肌电信号中提取时域特征和频域特征以及基于表面肌电信号获取的与所述指定动作的完成情况相关的特征。The features used to train the Parkinson's patient's motor ability classifier corresponding to the specified action include extracting time-domain features and frequency-domain features from the surface electromyography signal and obtaining information related to the completion of the specified action based on the surface electromyography signal. Characteristics. 9.根据权利要求1所述的系统,其中与所述指定动作的完成情况相关的特征包括下列中的至少一个:完成指定动作时的最大肌电信号、完成指定动作所用的时间。9. The system according to claim 1, wherein the features related to the completion of the specified action include at least one of the following: the maximum electromyographic signal when the specified action is completed, and the time taken to complete the specified action. 10.根据权利要求1或2所述的方法,还包括用于训练与指定动作对应的帕金森病人运动能力分类器的训练装置,其被配置为:10. The method according to claim 1 or 2, further comprising a training device for training a Parkinson's patient motor capacity classifier corresponding to a specified action, configured to: a)来自多个不同程度的帕金森病人和正常人执行所述指定动作时由其佩戴的肌电传感器采集的多个表面肌电信号作为训练数据集;a) multiple surface electromyographic signals collected by the electromyographic sensors worn by the Parkinson's patients and normal people of different degrees when performing the specified actions as training data sets; b)从各表面肌电信号中提取时域特征和频域特征以及基于表面肌电信号获取与所述指定动作的完成情况相关的特征;b) extracting time-domain features and frequency-domain features from each surface electromyography signal and obtaining features related to the completion of the specified action based on the surface electromyography signal; c)基于各个特征对于训练数据集的信息增益选择用于训练所述运动能力分类器的特征,其中每个特征对训练数据集的信息增益为训练数据集的经验熵与在给定该特征的条件下训练数据集的经验条件熵之间的差;C) select the feature for training the athletic ability classifier based on the information gain of each feature for the training data set, wherein the information gain of each feature for the training data set is the experience entropy of the training data set and given this feature The difference between the empirical conditional entropy of the training dataset under condition ; d)基于所选择的特征来训练所述运动能力分类器。d) Training the athletic ability classifier based on the selected features.
CN201710463520.XA 2017-06-19 2017-06-19 Quantitative assessment method and system for hand motion ability of Parkinson's disease people Active CN107157450B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710463520.XA CN107157450B (en) 2017-06-19 2017-06-19 Quantitative assessment method and system for hand motion ability of Parkinson's disease people

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710463520.XA CN107157450B (en) 2017-06-19 2017-06-19 Quantitative assessment method and system for hand motion ability of Parkinson's disease people

Publications (2)

Publication Number Publication Date
CN107157450A true CN107157450A (en) 2017-09-15
CN107157450B CN107157450B (en) 2020-03-31

Family

ID=59818815

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710463520.XA Active CN107157450B (en) 2017-06-19 2017-06-19 Quantitative assessment method and system for hand motion ability of Parkinson's disease people

Country Status (1)

Country Link
CN (1) CN107157450B (en)

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108309295A (en) * 2018-02-11 2018-07-24 宁波工程学院 A kind of arm muscular strength assessment method
CN108670713A (en) * 2018-04-23 2018-10-19 中国科学院深圳先进技术研究院 A kind of exoskeleton robot
CN108742612A (en) * 2018-05-24 2018-11-06 杭州航弈生物科技有限责任公司 DBS validation checking equipment based on myoelectricity marker
CN109009148A (en) * 2018-08-24 2018-12-18 广东工业大学 A kind of gait function appraisal procedure
CN110516762A (en) * 2019-10-10 2019-11-29 深圳大学 A method, device, storage medium and intelligent terminal for quantitative assessment of muscle state
CN110547807A (en) * 2019-09-17 2019-12-10 深圳市赛梅斯凯科技有限公司 driving behavior analysis method, device, equipment and computer readable storage medium
CN111210912A (en) * 2020-01-14 2020-05-29 上海恩睦信息科技有限公司 Parkinson prediction method and device
CN111415746A (en) * 2020-04-22 2020-07-14 上海邦邦机器人有限公司 Physical function evaluation model generation method, physical function evaluation method, and physical function evaluation apparatus
CN111528842A (en) * 2020-05-26 2020-08-14 复嶂环洲生物科技(上海)有限公司 Quantitative assessment method for Parkinson disease symptoms based on physiological and behavioral indexes
CN112826504A (en) * 2021-01-07 2021-05-25 中新国际联合研究院 A gamified Parkinson's symptom level assessment method and device
CN116467647A (en) * 2023-03-08 2023-07-21 浙江大学 Parkinson's disease artificial intelligence classification method and system based on three-dimensional finger movement
CN117643456A (en) * 2024-01-29 2024-03-05 北京航空航天大学 Auxiliary evaluation system, method and storage medium for parkinsonism
CN118430798A (en) * 2023-11-20 2024-08-02 河南千唯电子科技有限公司 Parkinson's disease intelligent screening and assessment system based on multi-data drive

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050234309A1 (en) * 2004-01-07 2005-10-20 David Klapper Method and apparatus for classification of movement states in Parkinson's disease
US20100176952A1 (en) * 2008-12-04 2010-07-15 The Regents Of The University Of California System for detection of body motion
CN104207782A (en) * 2014-08-22 2014-12-17 南昌大学 Hand motion function rating device for Parkinson's disease patient
CN104398263A (en) * 2014-12-25 2015-03-11 中国科学院合肥物质科学研究院 Method for quantitatively evaluating symptoms of tremor of patient with Parkinson's disease according to approximate entropy and cross approximate entropy
CN104434129A (en) * 2014-12-25 2015-03-25 中国科学院合肥物质科学研究院 Quantization evaluating device and method for dyskinesia symptoms of Parkinson and related extrapyramidal diseases
CN104522949A (en) * 2015-01-15 2015-04-22 中国科学院苏州生物医学工程技术研究所 Smart wristband for quantitatively evaluating motion function of Parkinson patient
CN105426696A (en) * 2015-12-24 2016-03-23 中国科学院苏州生物医学工程技术研究所 Multi-node quantitative assessment system and method for symptoms of Parkinson's disease
CN105426842A (en) * 2015-11-19 2016-03-23 浙江大学 Support vector machine based surface electromyogram signal multi-hand action identification method
CN105550583A (en) * 2015-12-22 2016-05-04 电子科技大学 Random forest classification method based detection method for malicious application in Android platform
CN105930663A (en) * 2016-04-26 2016-09-07 北京科技大学 Parkinson's disease early diagnosis method

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050234309A1 (en) * 2004-01-07 2005-10-20 David Klapper Method and apparatus for classification of movement states in Parkinson's disease
US20100176952A1 (en) * 2008-12-04 2010-07-15 The Regents Of The University Of California System for detection of body motion
CN104207782A (en) * 2014-08-22 2014-12-17 南昌大学 Hand motion function rating device for Parkinson's disease patient
CN104398263A (en) * 2014-12-25 2015-03-11 中国科学院合肥物质科学研究院 Method for quantitatively evaluating symptoms of tremor of patient with Parkinson's disease according to approximate entropy and cross approximate entropy
CN104434129A (en) * 2014-12-25 2015-03-25 中国科学院合肥物质科学研究院 Quantization evaluating device and method for dyskinesia symptoms of Parkinson and related extrapyramidal diseases
CN104522949A (en) * 2015-01-15 2015-04-22 中国科学院苏州生物医学工程技术研究所 Smart wristband for quantitatively evaluating motion function of Parkinson patient
CN105426842A (en) * 2015-11-19 2016-03-23 浙江大学 Support vector machine based surface electromyogram signal multi-hand action identification method
CN105550583A (en) * 2015-12-22 2016-05-04 电子科技大学 Random forest classification method based detection method for malicious application in Android platform
CN105426696A (en) * 2015-12-24 2016-03-23 中国科学院苏州生物医学工程技术研究所 Multi-node quantitative assessment system and method for symptoms of Parkinson's disease
CN105930663A (en) * 2016-04-26 2016-09-07 北京科技大学 Parkinson's disease early diagnosis method

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108309295A (en) * 2018-02-11 2018-07-24 宁波工程学院 A kind of arm muscular strength assessment method
CN108670713A (en) * 2018-04-23 2018-10-19 中国科学院深圳先进技术研究院 A kind of exoskeleton robot
CN108742612A (en) * 2018-05-24 2018-11-06 杭州航弈生物科技有限责任公司 DBS validation checking equipment based on myoelectricity marker
CN108742612B (en) * 2018-05-24 2023-12-12 王开亮 DBS effectiveness detection equipment based on myoelectric marker
CN109009148A (en) * 2018-08-24 2018-12-18 广东工业大学 A kind of gait function appraisal procedure
CN110547807A (en) * 2019-09-17 2019-12-10 深圳市赛梅斯凯科技有限公司 driving behavior analysis method, device, equipment and computer readable storage medium
CN110516762B (en) * 2019-10-10 2022-11-15 深圳大学 A method, device, storage medium and intelligent terminal for quantitative evaluation of muscle state
CN110516762A (en) * 2019-10-10 2019-11-29 深圳大学 A method, device, storage medium and intelligent terminal for quantitative assessment of muscle state
CN111210912A (en) * 2020-01-14 2020-05-29 上海恩睦信息科技有限公司 Parkinson prediction method and device
CN111415746A (en) * 2020-04-22 2020-07-14 上海邦邦机器人有限公司 Physical function evaluation model generation method, physical function evaluation method, and physical function evaluation apparatus
CN111528842A (en) * 2020-05-26 2020-08-14 复嶂环洲生物科技(上海)有限公司 Quantitative assessment method for Parkinson disease symptoms based on physiological and behavioral indexes
CN111528842B (en) * 2020-05-26 2023-01-03 复嶂环洲生物科技(上海)有限公司 Quantitative assessment method for Parkinson disease symptoms based on physiological and behavioral indexes
CN112826504A (en) * 2021-01-07 2021-05-25 中新国际联合研究院 A gamified Parkinson's symptom level assessment method and device
CN112826504B (en) * 2021-01-07 2024-03-26 中新国际联合研究院 Game parkinsonism grade assessment method and device
CN116467647A (en) * 2023-03-08 2023-07-21 浙江大学 Parkinson's disease artificial intelligence classification method and system based on three-dimensional finger movement
CN116467647B (en) * 2023-03-08 2025-11-25 浙江大学 An Artificial Intelligence Classification Method and System for Parkinson's Disease Based on Three-Dimensional Finger Movement
CN118430798A (en) * 2023-11-20 2024-08-02 河南千唯电子科技有限公司 Parkinson's disease intelligent screening and assessment system based on multi-data drive
CN117643456A (en) * 2024-01-29 2024-03-05 北京航空航天大学 Auxiliary evaluation system, method and storage medium for parkinsonism
CN117643456B (en) * 2024-01-29 2024-05-24 北京航空航天大学 Auxiliary evaluation system, method and storage medium for Parkinson's neuropathy

Also Published As

Publication number Publication date
CN107157450B (en) 2020-03-31

Similar Documents

Publication Publication Date Title
CN107157450A (en) Quantitative estimation method and system are carried out for the hand exercise ability to patient Parkinson
Zhang et al. Heart sound classification based on scaled spectrogram and partial least squares regression
Zhao et al. Towards low-cost sign language gesture recognition leveraging wearables
Arora et al. High accuracy discrimination of Parkinson's disease participants from healthy controls using smartphones
CN102961203B (en) Method for identifying surface electromyography (sEMG) on basis of empirical mode decomposition (EMD) sample entropy
WO2019100560A1 (en) Artificial intelligence self-learning-based automatic electrocardiography analysis method and apparatus
CN110598676B (en) Deep Learning Gesture EMG Recognition Method Based on Confidence Score Model
CN104127181A (en) Electromyographic signal tumble detection method based on WKFDA
CN108305680A (en) Intelligent parkinsonism aided diagnosis method based on multi-element biologic feature and device
Gallón et al. Comparison of machine learning algorithms and feature extraction techniques for the automatic detection of surface EMG activation timing
McCool et al. Lower arm electromyography (EMG) activity detection using local binary patterns
Pereira et al. Physiotherapy Exercises Evaluation using a Combined Approach based on sEMG and Wearable Inertial Sensors.
Li et al. Adaptive detection of Ahead-sEMG based on short-time energy of local-detail difference and recognition in advance of upper-limb movements
Choi et al. Earppg: Securing your identity with your ears
Li et al. Muscle fatigue identification and prediction in motion using wearable device with power and torque-based features
JP6103591B2 (en) Auscultation heart sound signal processing method, auscultation heart sound signal processing apparatus, and program for processing auscultation heart sound signal
Wang et al. OFGST-Swin: swin transformer utilizing overlap fusion-based generalized S-transform for respiratory cycle classification
Jiménez-Grande et al. A machine learning approach for the identification of kinematic biomarkers of chronic neck pain during single-and dual-task gait
CN112863515B (en) Method for recognizing human tongue and jaw movement using mobile phone speaker and microphone
Islam et al. Myoelectric pattern recognition performance enhancement using nonlinear features
Tomczyński et al. Hand gesture-based interface with multichannel sEMG band enabling unknown gesture discrimination
Kanoga et al. Subject-transfer framework with unlabeled data based on multiple distance measures for surface electromyogram pattern recognition
Song et al. A wearable system for monitoring neurological disorder events with multi-class classification model in daily life
Patel et al. Different Transfer Learning Approaches for Recognition of Lung Sounds
Siddiqi et al. mmJaw: Remote jaw gesture recognition with COTS mmWave radar

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