CN111544006B - 帕金森病人运动障碍量化及识别的可穿戴设备 - Google Patents
帕金森病人运动障碍量化及识别的可穿戴设备 Download PDFInfo
- Publication number
- CN111544006B CN111544006B CN202010414484.XA CN202010414484A CN111544006B CN 111544006 B CN111544006 B CN 111544006B CN 202010414484 A CN202010414484 A CN 202010414484A CN 111544006 B CN111544006 B CN 111544006B
- Authority
- CN
- China
- Prior art keywords
- ankle
- parkinson
- wrist
- measurement unit
- signal
- 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.)
- Active
Links
- 208000012661 Dyskinesia Diseases 0.000 title claims abstract description 35
- 208000018737 Parkinson disease Diseases 0.000 title claims abstract description 35
- 210000003423 ankle Anatomy 0.000 claims abstract description 72
- 230000033001 locomotion Effects 0.000 claims abstract description 69
- 210000000707 wrist Anatomy 0.000 claims abstract description 64
- 238000000034 method Methods 0.000 claims abstract description 26
- 238000012545 processing Methods 0.000 claims abstract description 21
- 238000011002 quantification Methods 0.000 claims abstract description 19
- 238000005259 measurement Methods 0.000 claims description 52
- 230000005021 gait Effects 0.000 claims description 35
- 238000004891 communication Methods 0.000 claims description 13
- 230000006870 function Effects 0.000 claims description 13
- 239000013598 vector Substances 0.000 claims description 12
- 238000012706 support-vector machine Methods 0.000 claims description 11
- 230000001133 acceleration Effects 0.000 claims description 7
- 238000000546 chi-square test Methods 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 4
- 238000012952 Resampling Methods 0.000 claims description 3
- 230000002194 synthesizing effect Effects 0.000 claims description 3
- 239000000284 extract Substances 0.000 claims description 2
- 230000008447 perception Effects 0.000 abstract description 2
- 208000016285 Movement disease Diseases 0.000 abstract 2
- 206010044565 Tremor Diseases 0.000 description 7
- 208000024891 symptom Diseases 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 4
- 238000003745 diagnosis Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 4
- 238000000605 extraction Methods 0.000 description 4
- 210000001364 upper extremity Anatomy 0.000 description 4
- 210000003141 lower extremity Anatomy 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 206010006100 Bradykinesia Diseases 0.000 description 2
- 208000006083 Hypokinesia Diseases 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 210000002683 foot Anatomy 0.000 description 2
- 210000004247 hand Anatomy 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 208000017667 Chronic Disease Diseases 0.000 description 1
- 206010017577 Gait disturbance Diseases 0.000 description 1
- 208000002740 Muscle Rigidity Diseases 0.000 description 1
- 238000009098 adjuvant therapy Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 238000002790 cross-validation Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000008451 emotion Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000037230 mobility Effects 0.000 description 1
- 230000004770 neurodegeneration Effects 0.000 description 1
- 208000015122 neurodegenerative disease Diseases 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000284 resting effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000007958 sleep Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
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/1118—Determining activity level
-
- 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/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/1121—Determining geometric values, e.g. centre of rotation or angular range of movement
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/40—Detecting, measuring or recording for evaluating the nervous system
- A61B5/4076—Diagnosing or monitoring particular conditions of the nervous system
- A61B5/4082—Diagnosing or monitoring movement diseases, e.g. Parkinson, Huntington or Tourette
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
- A61B5/6801—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
- A61B5/6813—Specially adapted to be attached to a specific body part
- A61B5/6824—Arm or wrist
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
- A61B5/6801—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
- A61B5/6813—Specially adapted to be attached to a specific body part
- A61B5/6829—Foot or ankle
-
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Physiology (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Dentistry (AREA)
- Neurology (AREA)
- Neurosurgery (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- Geometry (AREA)
- Developmental Disabilities (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
本发明属于惯性信号感知识别和帕金森智能医疗技术领域,具体涉及了一种帕金森病人运动障碍量化及识别的可穿戴设备,旨在解决缺少可靠无创的可穿戴设备来准确地实现帕金森病人的运动障碍量化及识别的问题。本发明包括:通过穿戴于被测对象和健康者身体不同侧的手腕惯性测量单元、脚踝惯性测量单元采集手腕运动信号、脚踝运动信号,并提取信号的特征;将提取的特征发送到输出处理终端,通过预设的帕金森病人运动障碍量化及识别方法,获取被测对象的病情程度并输出。本发明结构简单、容易操作、尺寸合适方便佩戴,功耗低,并且结合手腕与脚踝的运动信号,判断结果准确可靠。
Description
技术领域
本发明属于惯性信号感知识别和帕金森智能医疗技术领域,具体涉及了一种帕金森病人运动障碍量化及识别的可穿戴设备。
背景技术
帕金森病(PD)是多发于老年人的神经退行性疾病,作为一种典型的老年慢性疾病,帕金森病在老年人群中患病率成倍增加,所以关于帕金森的辅助治疗不可忽视。帕金森病的主要临床特征为运动障碍包括动作迟缓、手部震颤、肢体僵直和姿势步态异常等。在临床工作中,临床医师对PD患者的判断主要依赖于主观的临床经验,难以提供客观可靠的定量评价。
随着惯性传感器的可穿戴计算技术的发展,基于惯性传感器以及智能终端的可穿戴设备已经成为研究热点之一,其中医疗设备是可穿戴设备最具前景的应用领域之一。理想的可穿戴设备必须满足操作简单,成本低,适合日常临床试验的要求。此外,要求可穿戴设备提供的数据是可靠,准确和易于使用的。PD运动障碍定量评估方法有足迹分析,震颤分析,智能脚环等。足迹分析和震颤分析简单易行,不受试验场地限制,但其测量精度低,远不能满足现代运动评估要求。智能脚环结构简单,但没有考虑到身体其他部分可能存在的运动症状。由于人体运动监测的复杂性和多样性,基于惯性传感器的帕金森检测***依旧面临着许多技术问题,如何减小测量偏差是值得探究的问题。
总的来说,目前尚无可靠无创的针对帕金森病人运动障碍量化的可穿戴设备,能够准确地实现帕金森病人的运动障碍量化及识别。
发明内容
为了解决现有技术中的上述问题,即缺少可靠无创的可穿戴设备来准确地实现帕金森病人的运动障碍量化及识别的问题,本发明提供了一种帕金森病人运动障碍量化及识别的可穿戴设备,该可穿戴设备包括运动信号测量模块、数据处理终端;
所述运动信号测量模块包括手腕惯性测量单元、脚踝惯性测量单元,用于获取被测对象及健康者行走、站立、转身时的手腕运动信号、脚踝运动信号,进行信号的特征提取并将提取的特征发送至数据处理终端;
所述数据处理终端基于所述提取的特征,通过预设的帕金森病人运动障碍量化及识别方法,获取被测对象的病情程度并输出。
在一些优选的实施例中,所述手腕惯性测量单元的X轴指向手指方向;所述脚踝惯性测量单元的X轴指向脚趾方向;所述手腕惯性测量单元、脚踝惯性测量单元的X轴分别为其内置的惯性传感器的X轴。
在一些优选的实施例中,所述手腕惯性测量单元、脚踝惯性测量单元分别设置于所述被测对象及健康者的身体不同侧。
在一些优选的实施例中,所述手腕惯性测量单元基于手腕运动信号提取的特征包括标准差、能量绝对值、运动信号最大值、分位数、惯性传感器中陀螺仪Y轴和Z轴的交叉过零点;所述脚踝惯性测量单元基于脚踝运动信号提取的特征包括分位数、运动信号最大值。
在一些优选的实施例中,所述脚踝惯性测量单元设置有一个低功耗蓝牙通讯模块,用于将提取的脚踝特征发送至所述手腕惯性测量单元。
在一些优选的实施例中,所述手腕惯性测量单元设置有两个低功耗蓝牙通讯模块,分别用于接收脚踝惯性测量单元发送的特征以及将手腕和脚踝的特征发送至所述数据处理终端。
在一些优选的实施例中,所述预设的帕金森病人运动障碍量化及识别方法为:
步骤S10,对手腕运动信号、脚踝运动信号进行重采样,并通过预设的第一算法将重采样后的信号合成加速度器信号,并通过五阶椭圆带通滤波器对所述加速度器信号、手腕运动信号、脚踝运动信号滤波;
步骤S20,基于滤波后的加速度器信号Z轴的信号,通过预设的第一滑动窗口获取行走间隔;
步骤S30,基于滤波后的手腕运动信号及脚踝运动信号,通过预设的第二滑动窗口在所述行走间隔内分别提取被测对象及健康者的步态特征,并通过卡方检验方法分别选择被测对象及健康者的设定数量的步态特征作为最佳步态特征;
步骤S40,将所述最佳步态特征归一化为步态特征向量,通过训练好的支持向量机进行所述步态特征向量的分类,并通过预设的第二算法计算被测对象相对于健康者的阳性率Pr;
步骤S50,若Pr>0.9,则被测对象为重度帕金森病人;若0.9>Pr≥0.6,则被测对象为中度帕金森病人;若0.6>Pr≥0.5,则被测对象为轻度帕金森病人;若Pr<0.5,则被测对象为非帕金森病人。
在一些优选的实施例中,所述预设的第一算法为:
其中,As为合成的加速度器信号,Ax、Ay、Az分别代表惯性传感器中三轴加速度计X、Y、Z轴分量的加速度信号。
在一些优选的实施例中,步骤S20包括:
步骤S21,通过预设的第一滑动窗口计算所述加速度器信号Z轴的信号功率,并以大于设定上限的信号功率对应的点作为起始位置,小于设定下限的信号功率对应的点作为结束位置;
步骤S22,计算所述起始位置与结束位置之间的距离,以所述距离大于设定阈值的作为一个时间间隔,得到行走间隔。
在一些优选的实施例中,步骤S30中“基于滤波后的手腕运动信号及脚踝运动信号,通过预设的第二滑动窗口在所述行走间隔内分别提取被测对象及健康者的步态特征”,其方法为:
步骤S31,通过预设的第二滑动窗口在所述行走间隔内分别获取被测对象及健康者的重采样的手腕运动信号、脚踝运动信号的窗口数据;
步骤S32,将所述窗口数据作为数据帧,并分别提取所述数据帧内被测对象及健康者的步态特征。
在一些优选的实施例中,所述支持向量机的核函数为径向基函数,所述径向基函数为:
在一些优选的实施例中,所述预设的第二算法为:
其中,Dp为支持向量机分类器的阳性结果输出的数量,D代表在设定时间段内第二滑动窗口获取的数据帧的总数量。
本发明的有益效果:
本发明帕金森病人运动障碍量化及识别的可穿戴设备,结构简单、容易操作、尺寸合适方便佩戴,功耗低,单次充电可使用24小时以上。并且,不只考虑单一部分的运动情况,而是通过手腕和脚踝的测量单元反映上肢和下肢的运动情况。通过SVM-RBF算法区分正常人和PD患者,并设置pr值(小于1)来表明帕金森病的严重程度,一共为四个程度,分别是健康、轻度、中度、重度,值越大说明越严重。通过实验验证本发明的准确率可达到87.01%,并且与专家的诊断具有较高的一致性(二次加权Kappa系数k=0.961),为日常生活中帕金森病症状的持续监测提供了一种客观可靠的辅助办法。
附图说明
通过阅读参照以下附图所作的对非限制性实施例所作的详细描述,本申请的其它特征、目的和优点将会变得更明显:
图1是本发明帕金森病人运动障碍量化及识别的可穿戴设备的框架示意图;
图2是本发明帕金森病人运动障碍量化及识别的可穿戴设备一种实施例的手腕、脚踝惯性测量单元佩戴位置示意图;
图3是本发明帕金森病人运动障碍量化及识别的可穿戴设备一种实施例的手腕惯性测量单元工作流程示意图;
图4是本发明帕金森病人运动障碍量化及识别的可穿戴设备一种实施例的脚踝惯性测量单元工作流程示意图。
具体实施方式
下面结合附图和实施例对本申请作进一步的详细说明。可以理解的是,此处所描述的具体实施例仅用于解释相关发明,而非对该发明的限定。另外还需要说明的是,为了便于描述,附图中仅示出了与有关发明相关的部分。
需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本申请。
本发明的一种帕金森病人运动障碍量化及识别的可穿戴设备,该可穿戴设备包括运动信号测量模块、数据处理终端;
所述运动信号测量模块包括手腕惯性测量单元、脚踝惯性测量单元,用于获取被测对象及健康者行走、站立、转身时的手腕运动信号、脚踝运动信号,进行信号的特征提取并将提取的特征发送至数据处理终端;
所述数据处理终端基于所述提取的特征,通过预设的帕金森病人运动障碍量化及识别方法,获取被测对象的病情程度并输出。
为了更清晰地对本发明帕金森病人运动障碍量化及识别的可穿戴设备进行说明,下面结合图1对本发明实施例中各模块展开详述。
本发明一种实施例的帕金森病人运动障碍量化及识别的可穿戴设备,包括运动信号测量模块、数据处理终端,各模块详细描述如下:
运动信号测量模块包括手腕惯性测量单元、脚踝惯性测量单元,用于获取被测对象及健康者的行走、站立、转身时的手腕运动信号、脚踝运动信号,进行手腕运动信号、脚踝运动信号的特征提取并将提取的特征发送至数据处理终端。
脚踝惯性测量单元主要包括惯性传感器、微程序控制器、数据缓存模块、低功耗蓝牙通讯模块:
惯性传感器:六轴惯性传感器,包含3D陀螺仪和3D加速度计,每隔十五分钟测量30秒的惯性传感器数据。
微程序控制器:基于UCOSⅢ操作***进行数据采集,特征提取,数据存储和无线通信。
数据缓存模块:微程序控制器对采集的数据进行特征提取,提取后的特征值存储于数据缓存模块。
低功耗蓝牙通讯模块:相比ZigBee和WIFI,蓝牙功耗更低、传输更加稳定。在脚踝惯性测量单元与手腕惯性测量单元之间,采用低功耗蓝牙传输特征数据。
手腕惯性测量单元比脚踝惯性测量单元多一个低功耗蓝牙通讯模块,而手腕惯性测量单元有两个低功耗蓝牙通讯模块。除了与脚踝惯性测量单元进行数据传输的低功耗蓝牙通讯模块,还有一个将手腕和脚踝的特征发送至数据处理终端的低功耗蓝牙通讯模块。
如图2所示,为本发明帕金森病人运动障碍量化及识别的可穿戴设备一种实施例的手腕、脚踝惯性测量单元佩戴位置示意图,脚踝惯性测量单元和手腕惯性测量单元的佩戴位置在身体的不同侧,并且需要调整惯性测量单元的位置,使其X轴分别指向脚趾和手指(惯性测量单元的X轴为其内置的惯性传感器的X轴,惯性传感器的说明书中指明了X轴、Y轴、Z轴)。由于70%PD病患的初始症状是手部震颤,而且多数开始于一侧上肢末端。当休息或停止动作时症状发作,情绪紧张时加重,睡觉时消失。在行走过程中,手部的静态震颤往往会加重。另一种常见症状是以动作迟缓,难以开始运动以及逐渐失去活动能力为特征的运动障碍。在临床诊断中,患者的表现是手部摆动幅度逐渐减小甚至消失,缓慢行走、步距变小且经常拖行。手腕惯性测量单元主要检测手部震颤和上肢摆动情况,脚踝惯性测量单元主要检测下肢运动情况。
惯性传感器每隔十五分钟选取30秒时间采集数据,并对采集到的数据进行特征提取,得到一组特征值。脚踝惯性测量单元的特征值为:分位数(0.9)和运动信号最大值,其代表了下肢的运动能力。手腕惯性测量单元的特征值分别为:标准差、能量绝对值、运动信号最大值、分位数(0.9)、惯性传感器中陀螺仪Y轴和Z轴的交叉过零点。其标准差、能量绝对值、分位数(0.9)代表行走时上肢的摆动情况,交叉过零点象征手部震颤情况。每隔十五分钟,脚踝惯性测量单元的特征值传输至手腕惯性测量单元,手腕惯性测量单元重新打包两个测量单元的特征值并存储于数据缓存模块。因此,脚踝惯性测量单元只有一个低功耗蓝牙通讯模块,用以传输脚踝特征值;而手腕惯性测量单元有两个低功耗蓝牙通讯模块,用以接收脚踝惯性测量单元的特征值和发送两个测量单元的特征值到数据处理终端。
数据处理终端基于所述提取的特征,通过预设的帕金森病人运动障碍量化及识别方法,获取被测对象的病情程度并输出。
本发明一个实施例中,选用智能手机作为数据处理终端,在其他的实施例中,还可以根据需要其他智能终端,本发明在此不做详述。
如图3所示,为本发明帕金森病人运动障碍量化及识别的可穿戴设备一种实施例的手腕惯性测量单元工作流程示意图,上电复位后,***外设初始化、操作***初始化;然后,创建蓝牙连接状态检查任务、创建LED指示任务、创建数据采集任务、创建***15分钟定时唤醒服务;若不到唤醒时,持续等待,若***唤醒延时,连接唤醒脚踝从机(即脚踝惯性测量单元);若连接不成功则持续尝试连接,若连接成功,则判断是否有数据上报;若有则接受数据后存储于FLASH(即数据缓存模块),向脚踝从机(即脚踝惯性测量单元)发送开始采集指令,若没有数据上报则直接跳至向脚踝从机(即脚踝惯性测量单元)发送开始采集指令;开始采集数据30s,进行数据特征提取,接受从机(即脚踝惯性测量单元)数据,进行数据特征提取后存储于FLASH(即数据缓存模块);判断是否连接移动设备(即数据处理终端),若没有,***进入低功耗,若有,则向移动终端(即数据处理终端)上传特征数据后***进入低功耗;返回等待***唤醒延时状态。
如图4所示,为本发明帕金森病人运动障碍量化及识别的可穿戴设备一种实施例的脚踝惯性测量单元工作流程示意图,上电复位后,***外设初始化、操作***初始化;然后,创建蓝牙连接状态检查任务、创建LED指示任务、创建数据采集任务;判断是否有主机蓝牙(主机即手腕惯性测量单元,主机蓝牙即手腕惯性测量单元与数据处理终端之间的蓝牙)连接,若没有,持续等待,若有,唤醒***并进入工作模式;上传更新数据,判断是否有数据采集指令,若没有则等待,若有,开始采集数据30s;进行数据特征提取后存储于FLASH(即数据缓存模块);判断蓝牙是否在连(即脚踝惯性测量单元与手腕惯性测量单元处于蓝牙连接状态),若没有则判断蓝牙连接是否超时,若不是,则返回等待蓝牙连接,若蓝牙连接超时,则***进入低功耗,若蓝牙在连,则向手腕主机(手腕惯性测量单元)上传更新数据后***进入低功耗;返回等待主机蓝牙连接状态。
预设的帕金森病人运动障碍量化及识别方法为:
步骤S10,对手腕运动信号、脚踝运动信号进行重采样,并通过预设的第一算法将重采样后的信号合成加速度器信号,并通过五阶椭圆带通滤波器对所述加速度器信号、手腕运动信号、脚踝运动信号滤波。
本发明一个实施例中,分别将采集的帕金森病人、健康者行走、站立、转身时的手腕运动信号、脚踝运动信号重采样到50Hz以此减少计算复杂程度,并通过预设的第一算法将重采样后的信号合成加速度器信号,如式(1)所示:
其中,As为合成的加速度器信号,Ax、Ay、Az分别代表惯性传感器中三轴加速度计X、Y、Z轴分量的加速度信号。
在进行步骤S30的行走间隔获取之前,还可以通过五阶椭圆带通滤波器对加速度器信号、重采样后的手腕运动信号、脚踝运动信号进行滤波,滤波是为了消除重力分量和高频噪声的影响。
步骤S20,基于滤波后的加速度器信号Z轴的信号,通过预设的第一滑动窗口获取行走间隔。
步骤S21,通过预设的第一滑动窗口计算所述加速度器信号Z轴的信号功率,并以大于设定上限的信号功率对应的点作为起始位置,小于设定下限的信号功率对应的点作为结束位置。本发明一个实施例中,预设的第一滑动窗口的长度设置为20个采样点数。
步骤S22,计算所述起始位置与结束位置之间的距离,以所述距离大于设定阈值的作为一个时间间隔,得到行走间隔。本发明一个实施例中,如果起始位置和结束位置之间的距离大于256个采样点数(即行走时间大于5s),则分割出一个时间间隔。
步骤S30,基于滤波后的手腕运动信号及脚踝运动信号,通过预设的第二滑动窗口在所述行走间隔内分别提取被测对象及健康者的步态特征,并通过卡方检验方法分别选择被测对象及健康者的设定数量的步态特征作为最佳步态特征。
“基于滤波后的手腕运动信号及脚踝运动信号,通过预设的第二滑动窗口在所述行走间隔内分别提取被测对象及健康者的步态特征”,其方法为:
步骤S31,通过预设的第二滑动窗口在所述行走间隔内分别获取被测对象及健康者的重采样的手腕运动信号、脚踝运动信号的窗口数据。
本发明一个实施例中,预设的第二滑动窗口大小为256,滑动步长为16个采样点数,将这个滑动窗口捕获的运动信号的窗口数据处理为数据帧。
步骤S32,将所述窗口数据作为数据帧,并分别提取所述数据帧内被测对象及健康者的步态特征。
为了提取最佳的步态特征,在进行特征归一化之前,本发明还通过卡方检验方法分别进行被测对象及健康者的步态特征筛选,最终得到不同被测对象及健康者的8个步态特征,分别为腕部惯性传感器上陀螺仪Z轴的数据标准差、陀螺仪Z轴的能量绝对值、陀螺仪Z轴过零点数量、陀螺仪Y轴过零点数量、合成加速度信号极值、合成加速度信号0.9分位数以及踝部合成加速度信号极值、合成加速度信号0.9分位数。
步骤S40,将所述最佳步态特征归一化为步态特征向量,通过训练好的支持向量机进行所述步态特征向量的分类,并通过预设的第二算法计算被测对象相对于健康者的阳性率Pr。
将步态特征归一化为步态特征向量,如式(2)所示:
其中,Xscaled代表归一化后的步态特征向量,X代表待归一化的步态特征,Xmax、Xmin分别为待归一化步态特征中的最大值和最小值。
这里使用支持向量机对帕金森病人和健康者的数据片段进行分类,支持向量机适用于解决小样本的非线性分类问题,这里类标签和所提取的特征之间呈现非线性关系,因此使用径向基函数(RFB,Radial Basis Function)作为支持向量机分类器的核函数,RBF核函数具有逼近非线性函数的能力,并能处理***中的不溶性规律,具有良好的泛化能力和较快的学习速度,如式(3)所示:
通过训练好的支持向量机对步态特征向量进行分类,并计算被测对象相对于健康者的阳性率Pr,如式(4)所示
其中,Dp为支持向量机分类器的阳性结果输出的数量,D代表在设定时间段内第二滑动窗口获取的数据帧的总数量。
Pr表示支持向量机分类器在SVM输出结果的30秒(1500个数据点)数据长度内,认为扫描数据段属于帕金森患者的比率。
步骤S50,若Pr>0.9,则被测对象为重度帕金森病人;若0.9>Pr≥0.6,则被测对象为中度帕金森病人;若0.6>Pr≥0.5,则被测对象为轻度帕金森病人;若Pr<0.5,则被测对象为非帕金森病人。
本发明还通过十倍交叉验证方法评估支持向量机分类器的性能,采用了常用的性能指标:准确度、召回率、F1分数。设置帕金森患者为阳性组,健康者为阴性组,其中准确度为将帕金森患者和健康者正确分类的比率,召回率为真阳性组的比率,F1分数为模型精度和召回率的加权平均值。本发明计算出帕金森患者与健康受试者的分类准确性大约是92.9%,敏感性(分类器将正例分类成正例的比率,也就是将帕金森病人分为病人的比例)为88.9%,特异性(分类器将反例分类成反例的比率,也就是将健康者分为健康者的比例)为100%,显示了本发明方法在识别帕金森患者与健康者方面的良好能力。
最后,通过二次加权Kappa系数对该方法进行定量的可信度评价,如式(5)所示:
其中,∑ωf0s是序数标度中的观测到的加权频率之和,∑ωfc是序数标度中的随机期望的加权频率之和。
二次加权ω的计算方法如式(6)所示:
其中,i-j是人工分类和本发明方法不一致的数量,k是序数标度上的等级数。
本发明一个实施例中,k取值为4,4个等级上的ω分别为1、0.89、0.56、0,表明,评估偏差越大,权重越大。经过实验测试,本发明与专家诊断相比,二次加权Kappa系数k=0.961,证明本发明方法与专家人工诊断具有较高的一致性,本发明诊断结果准确度高。
术语“第一”、“第二”等是用于区别类似的对象,而不是用于描述或表示特定的顺序或先后次序。
术语“包括”或者任何其它类似用语旨在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备/装置不仅包括那些要素,而且还包括没有明确列出的其它要素,或者还包括这些过程、方法、物品或者设备/装置所固有的要素。
至此,已经结合附图所示的优选实施方式描述了本发明的技术方案,但是,本领域技术人员容易理解的是,本发明的保护范围显然不局限于这些具体实施方式。在不偏离本发明的原理的前提下,本领域技术人员可以对相关技术特征作出等同的更改或替换,这些更改或替换之后的技术方案都将落入本发明的保护范围之内。
Claims (6)
1.一种帕金森病人运动障碍量化及识别的可穿戴设备,其特征在于,该可穿戴设备包括运动信号测量模块、数据处理终端;
所述运动信号测量模块包括手腕惯性测量单元、脚踝惯性测量单元,用于获取被测对象及健康者行走、站立、转身时的手腕运动信号、脚踝运动信号,进行信号的特征提取并将提取的特征发送至数据处理终端;
所述脚踝惯性测量单元设置有一个低功耗蓝牙通讯模块,用于将提取的脚踝特征发送至所述手腕惯性测量单元;所述脚踝惯性测量单元基于脚踝运动信号提取的特征包括分位数、运动信号最大值;
所述手腕惯性测量单元设置有两个低功耗蓝牙通讯模块,分别用于接收脚踝惯性测量单元发送的特征以及将手腕和脚踝的特征发送至所述数据处理终端;所述手腕惯性测量单元基于手腕运动信号提取的特征包括标准差、能量绝对值、运动信号最大值、分位数、惯性传感器中陀螺仪Y轴和Z轴的交叉过零点;所述手腕惯性测量单元的X轴指向手指方向,为其内置的惯性传感器的X轴;
所述数据处理终端基于所述提取的特征,通过预设的帕金森病人运动障碍量化及识别方法,获取被测对象的病情程度并输出;
其中,所述预设的帕金森病人运动障碍量化及识别方法为:
步骤S10,对手腕运动信号、脚踝运动信号进行重采样,并通过预设的第一算法将重采样后的信号合成加速度器信号,并通过五阶椭圆带通滤波器对所述加速度器信号、手腕运动信号、脚踝运动信号滤波;
步骤S20,通过预设的第一滑动窗口计算所述加速度器信号Z轴的信号功率,并以大于设定上限的信号功率对应的点作为起始位置,小于设定下限的信号功率对应的点作为结束位置,计算所述起始位置与结束位置之间的距离,以所述距离大于设定阈值的作为一个时间间隔,得到行走间隔;
步骤S30,通过预设的第二滑动窗口在所述行走间隔内分别获取被测对象及健康者的重采样的手腕运动信号、脚踝运动信号的窗口数据,将所述窗口数据作为数据帧,分别提取所述数据帧内被测对象及健康者的步态特征,并通过卡方检验方法分别选择被测对象及健康者的设定数量的步态特征作为最佳步态特征;
步骤S40,将所述最佳步态特征归一化为步态特征向量,通过训练好的支持向量机进行所述步态特征向量的分类,并通过预设的第二算法计算被测对象相对于健康者的阳性率Pr;所述支持向量机的核函数为径向基函数;
步骤S50,若Pr>0.9,则被测对象为重度帕金森病人;若0.9>Pr≥0.6,则被测对象为中度帕金森病人;若0.6>Pr≥0.5,则被测对象为轻度帕金森病人;若Pr<0.5,则被测对象为非帕金森病人。
2.根据权利要求1所述的帕金森病人运动障碍量化及识别的可穿戴设备,其特征在于,所述脚踝惯性测量单元的X轴指向脚趾方向;所述脚踝惯性测量单元的X轴为其内置的惯性传感器的X轴。
3.根据权利要求1或2所述的帕金森病人运动障碍量化及识别的可穿戴设备,其特征在于,所述手腕惯性测量单元、脚踝惯性测量单元分别设置于所述被测对象及健康者的身体不同侧。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010414484.XA CN111544006B (zh) | 2020-05-15 | 2020-05-15 | 帕金森病人运动障碍量化及识别的可穿戴设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010414484.XA CN111544006B (zh) | 2020-05-15 | 2020-05-15 | 帕金森病人运动障碍量化及识别的可穿戴设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111544006A CN111544006A (zh) | 2020-08-18 |
CN111544006B true CN111544006B (zh) | 2021-10-26 |
Family
ID=72004777
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010414484.XA Active CN111544006B (zh) | 2020-05-15 | 2020-05-15 | 帕金森病人运动障碍量化及识别的可穿戴设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111544006B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112674762A (zh) * | 2020-12-28 | 2021-04-20 | 江苏省省级机关医院 | 一种基于可穿戴惯性传感器的帕金森震颤评估装置 |
CN114171194B (zh) * | 2021-10-20 | 2022-09-06 | 中国科学院自动化研究所 | 帕金森多症状量化评估方法、装置、电子设备、及介质 |
CN117423452A (zh) * | 2022-07-11 | 2024-01-19 | 华为技术有限公司 | 一种用于评估帕金森病情的电子设备 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104398263A (zh) * | 2014-12-25 | 2015-03-11 | 中国科学院合肥物质科学研究院 | 一种基于近似熵和互近似熵的帕金森患者震颤症状量化评测方法 |
CN104729507A (zh) * | 2015-04-13 | 2015-06-24 | 大连理工大学 | 一种基于惯性传感器的步态识别方法 |
CN109069066A (zh) * | 2015-09-28 | 2018-12-21 | 卡斯西部储备大学 | 可穿戴和连接的步态分析*** |
CN109480858A (zh) * | 2018-12-29 | 2019-03-19 | 中国科学院合肥物质科学研究院 | 一种用于量化检测帕金森患者运动迟缓症状的可穿戴智能***及方法 |
CN110151189A (zh) * | 2019-04-30 | 2019-08-23 | 杭州电子科技大学 | 用于帕金森步态风险评估的非线性步态动力学判别方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7827011B2 (en) * | 2005-05-03 | 2010-11-02 | Aware, Inc. | Method and system for real-time signal classification |
CN104840191A (zh) * | 2015-04-30 | 2015-08-19 | 吴健康 | 测试心脏运动功能的设备、***和方法 |
US10226204B2 (en) * | 2016-06-17 | 2019-03-12 | Philips North America Llc | Method for detecting and responding to falls by residents within a facility |
CN106971059B (zh) * | 2017-03-01 | 2020-08-11 | 福州云开智能科技有限公司 | 一种基于神经网络自适应健康监测的可穿戴设备 |
US10182746B1 (en) * | 2017-07-25 | 2019-01-22 | Verily Life Sciences Llc | Decoupling body movement features from sensor location |
CN108720841A (zh) * | 2018-05-22 | 2018-11-02 | 上海交通大学 | 基于云检测的可穿戴下肢运动矫正*** |
CN108968918A (zh) * | 2018-06-28 | 2018-12-11 | 北京航空航天大学 | 早期帕金森的可穿戴辅助筛查设备 |
-
2020
- 2020-05-15 CN CN202010414484.XA patent/CN111544006B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104398263A (zh) * | 2014-12-25 | 2015-03-11 | 中国科学院合肥物质科学研究院 | 一种基于近似熵和互近似熵的帕金森患者震颤症状量化评测方法 |
CN104729507A (zh) * | 2015-04-13 | 2015-06-24 | 大连理工大学 | 一种基于惯性传感器的步态识别方法 |
CN109069066A (zh) * | 2015-09-28 | 2018-12-21 | 卡斯西部储备大学 | 可穿戴和连接的步态分析*** |
CN109480858A (zh) * | 2018-12-29 | 2019-03-19 | 中国科学院合肥物质科学研究院 | 一种用于量化检测帕金森患者运动迟缓症状的可穿戴智能***及方法 |
CN110151189A (zh) * | 2019-04-30 | 2019-08-23 | 杭州电子科技大学 | 用于帕金森步态风险评估的非线性步态动力学判别方法 |
Non-Patent Citations (1)
Title |
---|
基于运动传感器的帕金森自动分级研究;杨越等;《中国医疗设备》;20180910;第33卷(第09期);第37-41页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111544006A (zh) | 2020-08-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111544006B (zh) | 帕金森病人运动障碍量化及识别的可穿戴设备 | |
US10986465B2 (en) | Automated detection and configuration of wearable devices based on on-body status, location, and/or orientation | |
CN109480858B (zh) | 一种用于量化检测帕金森患者运动迟缓症状的可穿戴智能***及方法 | |
Bhat et al. | Online human activity recognition using low-power wearable devices | |
Curone et al. | A real-time and self-calibrating algorithm based on triaxial accelerometer signals for the detection of human posture and activity | |
JP5674766B2 (ja) | 着用位置を検出するためのセンシングデバイス | |
CN107580471B (zh) | 可穿戴脉搏感测设备信号质量估计 | |
JP5874130B2 (ja) | 人の歩行を検出するためのシステムおよび方法 | |
CN109069066A (zh) | 可穿戴和连接的步态分析*** | |
WO2015063520A1 (en) | Biomechanical activity monitoring | |
AU2006242134A1 (en) | Method and system for wearable vital signs and physiology, activity, and environmental monitoring | |
KR101307250B1 (ko) | 피검자의 트래머데이터 측정장치 | |
US8552871B2 (en) | Apparatus and method for sensing photoplethysmogram and fall | |
US20220287651A1 (en) | Method and system for analyzing biomechanical activity and exposure to a biomechanical risk factor on a human subject in a context of physical activity | |
CN103006195A (zh) | 一种基于图像处理的非接触式生命体征数据监测***和监测方法 | |
KR102376904B1 (ko) | 가속도 및 각속도 신호를 이용한 파킨슨병 서동증 진단 평가 장치 및 이를 이용 파킨슨병 서동증 결정 방법 | |
Rasheed et al. | Evaluation of human activity recognition and fall detection using android phone | |
Jatesiktat et al. | An elderly fall detection using a wrist-worn accelerometer and barometer | |
Lo et al. | A wireless sensor system for motion analysis of Parkinson's disease patients | |
WO2024012330A1 (zh) | 一种用于评估帕金森病情的电子设备 | |
EP4247262A1 (en) | Method and system for measuring and displaying biosignal data to a wearer of a wearable article | |
CN111544005B (zh) | 基于支持向量机的帕金森病人运动障碍量化及识别方法 | |
Dimitrakopoulos et al. | Tremor quantification through event-based movement trajectory modeling | |
KR101343506B1 (ko) | 피검자의 운동장애 진단을 위한 생체역학적 파라미터 추출장치 | |
CN115530774A (zh) | 癫痫检测方法和装置 |
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 |