CN111657918A - 一种结合心电和惯性传感数据的跌倒检测方法及系统 - Google Patents
一种结合心电和惯性传感数据的跌倒检测方法及系统 Download PDFInfo
- Publication number
- CN111657918A CN111657918A CN202010537555.5A CN202010537555A CN111657918A CN 111657918 A CN111657918 A CN 111657918A CN 202010537555 A CN202010537555 A CN 202010537555A CN 111657918 A CN111657918 A CN 111657918A
- Authority
- CN
- China
- Prior art keywords
- data
- inertial sensing
- fall
- ecg
- daily
- 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
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 59
- 230000000694 effects Effects 0.000 claims abstract description 56
- 230000009467 reduction Effects 0.000 claims abstract description 44
- 238000012706 support-vector machine Methods 0.000 claims abstract description 31
- 238000000034 method Methods 0.000 claims abstract description 29
- 238000013145 classification model Methods 0.000 claims abstract description 20
- 238000012545 processing Methods 0.000 claims abstract description 17
- 230000033001 locomotion Effects 0.000 claims abstract description 16
- 230000008569 process Effects 0.000 claims abstract description 16
- 238000000513 principal component analysis Methods 0.000 claims abstract description 13
- 238000010606 normalization Methods 0.000 claims abstract description 12
- 239000011159 matrix material Substances 0.000 claims description 33
- 238000000718 qrs complex Methods 0.000 claims description 21
- 230000009471 action Effects 0.000 claims description 16
- 210000000707 wrist Anatomy 0.000 claims description 13
- 238000000605 extraction Methods 0.000 claims description 9
- 230000001133 acceleration Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000010276 construction Methods 0.000 claims description 3
- 230000008859 change Effects 0.000 abstract description 8
- 230000007547 defect Effects 0.000 abstract description 4
- 230000001766 physiological effect Effects 0.000 abstract description 3
- 230000036544 posture Effects 0.000 description 8
- 230000000747 cardiac effect Effects 0.000 description 5
- 230000003321 amplification Effects 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 4
- 238000001914 filtration Methods 0.000 description 4
- 238000003199 nucleic acid amplification method Methods 0.000 description 4
- 238000001228 spectrum Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000034225 regulation of ventricular cardiomyocyte membrane depolarization Effects 0.000 description 2
- 230000013577 regulation of ventricular cardiomyocyte membrane repolarization Effects 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- 238000005452 bending Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013329 compounding Methods 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 210000002837 heart atrium Anatomy 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000028161 membrane depolarization Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000004660 morphological change Effects 0.000 description 1
- 210000004165 myocardium Anatomy 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000037081 physical activity Effects 0.000 description 1
- 230000035790 physiological processes and functions Effects 0.000 description 1
- 230000004213 regulation of atrial cardiomyocyte membrane depolarization Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction 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/1116—Determining posture transitions
-
- 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/1116—Determining posture transitions
- A61B5/1117—Fall detection
-
- 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/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
-
- 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
-
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2562/00—Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
- A61B2562/02—Details of sensors specially adapted for in-vivo measurements
- A61B2562/0219—Inertial sensors, e.g. accelerometers, gyroscopes, tilt switches
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Biophysics (AREA)
- Public Health (AREA)
- Pathology (AREA)
- Physiology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Veterinary Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Dentistry (AREA)
- Evolutionary Computation (AREA)
- Fuzzy Systems (AREA)
- Mathematical Physics (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
本发明公开了一种结合心电和惯性传感数据的跌倒检测方法及系统,包括:利用包括ECG和惯性传感单元的设备采集人体日常活动过程的跌倒数据和日常动作数据;将收集的所述跌倒数据和日常动作数据按相同时间段进行分割后,分别进行提取得到惯性传感数据特征和ECG数据特征,并进行归一化处理得到归一化样本数据特征;采用主成分分析算法对所述归一化样本数据特征进行降维处理得到降维数据特征,采用支持向量机算法对所述降维数据特征进行分类,得到支持向量机分类模型;将新采集的数据输入支持向量机分类模型得到跌倒检测结果。分别从人体生理活动和姿态两个不同维度反映人体跌倒过程的变化,避免单一类型数据不准确的缺陷,提高跌倒检测精度。
Description
技术领域
本发明涉及人体跌倒检测领域,具体涉及一种结合心电和惯性传感数据的跌倒检测方法 及系统。
背景技术
跌倒会造成老年人心理和身体伤害,是威胁老年人健康的重要因素之一。跌倒检测是指 采集人体活动数据,通过数据分析发现跌倒事件。用于人体活动数据采集的设备分为环境感 知设备和穿戴式设备。环境感知设备无需人体佩戴,但监测范围受限。穿戴设备主要是惯性 传感器,这类设备体积小,可随身佩戴,在跌倒检测领域应用较广,但其准确率受佩戴位置 影响。
近年来有研究将穿戴式心电图(ECG)检测仪用于检测跌倒。ECG信号代表心脏周期性 收缩和心肌松弛引起的生物电活动,可以反映人体在不同姿态及姿态转换过程中心脏活动的 变化情况。中国专利“CN201510697530.0一种基于心电图的跌倒检测的方法和系统”提供一 种基于心率的跌倒检测方法。通过检测终端获取气压和心率数据,首先通过气压判断人体高 度变化,若高度变化达到阈值,则启动心率检测,若心率变化值大于阈值,则判断为跌倒。 使用可穿戴设备观察人体在躺卧和站立姿势的ECG形态变化程度进行跌倒检测,准确率为 77.30%。
发明内容
本发明所要解决的技术问题是现有的检测技术中基于环境感知设备的跌倒检测存在隐私、 检测范围受限,基于惯性传感设备的跌倒检测精度受佩戴位置影响,而基于心电图的跌倒检 测主要通过心率是否达到阈值检测人体生理状态变换,未能充分利用心电图信息,检测精度 受影响的问题,目的在于提供一种结合心电和惯性传感数据的跌倒检测方法及系统,解决上 述问题。
本发明通过下述技术方案实现:
一种结合心电和惯性传感数据的跌倒检测方法,包括:
S1:利用包括ECG和惯性传感单元的腕部可穿戴设备采集人体日常活动过程的跌倒数据 和日常动作数据;
S2:将收集的所述跌倒数据和日常动作数据按相同时间段进行分割;
S3:分别提取分割后的所述跌倒数据和日常动作数据的特征得到惯性传感数据特征和 ECG数据特征;
S4:对所述惯性传感数据特征和ECG数据特征进行归一化处理得到归一化样本数据特征;
S5:采用主成分分析算法对所述归一化样本数据特征进行降维处理得到降维数据特征, 采用支持向量机算法对所述降维数据特征进行分类,得到支持向量机分类模型;
S6:将所述腕部可穿戴设备新采集的日常活动数据输入所述支持向量机分类模型,得到 跌倒检测结果。
本发明通过采集不同的日常活动数据即跌倒数据和正常活动数据,对包括向前滑倒、向 后滑倒、侧向滑倒、绊倒和晕倒的跌倒动作和包括走、跑、上下楼梯、姿态转换、弯腰的正 常日常活动进行模拟得到支持向量机分类模型,所述支持向量机分类模型得到跌倒或非跌倒 两种分类,因此,将输入的新的日常活动数据输入所述支持向量机分类模型即可得到跌倒检 测结果即跌倒或非跌倒。
本发明采用包括ECG即穿戴式心电图检测仪和惯性传感器的腕部可穿戴设备进行日常 活动中不同类型数据的采集,ECG数据在不同活动状态下的数据不同,可以弥补只有惯性传 感器进行检测时其检测范围受限以及跌倒检测精度受到佩戴位置影响而造成跌倒检测结果不 准确的不足,将两者进行结合,分别从人体生理活动和姿态两个不同维度反映人体跌倒过程 的变化,避免单一类型数据不准确的缺陷,从而得到更准确的跌倒数据,提高跌倒检测精度。
本发明的日常活动包含跌倒及正常日常活动。其中跌倒动作在软垫上进行模拟,包括向 前滑倒、向后滑倒、侧向滑倒、绊倒和晕倒;正常日常活动包括走、跑、上下楼梯、姿态转 换、弯腰等。
本发明中ECG数据与惯性传感数据的类型完全不同,且所述两种数据在正常活动和跌倒 得到的数据是不同的,取值范围不同以及量纲不同,因此需要将所有的数据通过相同时间段 进行分割并提取有效的惯性传感数据和ECG数据,并将惯性传感数据和ECG数据采用线性 归一化处理,将所有的数据映射到0~1;采用归一化函数将有量纲的表达式,经过变换,化 为无量纲的表达式,成为标量,可以简化计算过程从而减少检测过程的计算时间;为了在更 大程度上提高跌倒检测的精度,需要从惯性传感数据中提取多方面的特征,具体包括:算术 平均值、标准偏差、中位数绝对偏差、最大值、最小值、频率信号偏斜、频率信号峰度、最 大频率成分、平均能量;在ECG采集的数据中ECG信号的一个正常周期由QRS波、T波和 P波组成,其中P波对应左右心房去极化阶段,频谱范围10-15Hz。QRS复合波对应于左心 室和右心室去极化,频谱范围10-40Hz。T波反映心室复极化过程,在QRS复合波后300毫 秒左右的范围。因此需要对ECG信号的每一心动周期先进行QRS复合波确定,然后进行P 波和T波确定。
进一步的,所述S3包括:
S31:提取惯性传感单元中分割后的所述跌倒数据和日常动作数据的惯性传感数据得到惯 性传感数据特征;
所述惯性传感数据特征包括:算术平均值、标准偏差、中位数绝对偏差、最大值、最小 值、频率信号偏斜、频率信号峰度、最大频率成分、平均能量;
S32:提取分割后的所述跌倒数据和日常动作数据的ECG数据中的QRS复合波、P波和 T波数据特征得到所述ECG数据特征;
所述S4中的归一化处理函数:
其中max表示样本数据最大值,min表示样本数据最小值。x表示原始传感器数据,x* 表示归一化结果。
进一步的,所述QRS复合波的提取:
根据ECG数据的峰值提取QRS复合波的峰值点,在峰值点前后分别提取最低点,确定 QRS复合波的起始位置和偏移量,得到QRS复合波。
进一步的,使用以下时间差方程进行滤波处理得到ynT:
ynT=-2*χ(n-2)T-χ(n-1)T+χ(n+1)T+2*χ(n+2)T
其中,χnT表示时间T上大小为n的数据样本;
所述P波的提取公式:
θ=0.30*(μn-μHF)
其中,θ表示某一时间段ECG数据的动态阈值表达式,μn指的是在最近一段信号时间序 列内所有数据点的斜率平均值,μHF指一次跳动中的高频噪音;
其中,
μHF指一次跳动中的高频噪音,为了计算高频噪音,首先将信号传递到高通滤波器:
y′nT=χnT-2*χ(n-1)T+χ(n-2)T
其中,RA是R峰的振幅,K是通过实验获取的放大因子;
当θ达到阈值时,提取所述时间段ECG数据的峰值后,通过时间序列顺序,根据峰值前 后波形数据的斜率,提取斜率最低值点,确定P波起始位置和偏移量,得到P波。
进一步的,所述T波提取公式:
θ=0.30*(μn-μ′HF)
其中,θ表示某一时间段ECG数据的动态阈值表达式,μn指的是在最近一段信号时间序 列内所有数据点的斜率平均值,μ′HF指一次跳动中的低频噪音;
μ′HF指一次跳动中的低频噪音,为了计算低频噪音,首先通过2阶低通滤波器进行平滑处 理:
y″nT=2*y(n-1)T-y(n-2)T+χnT-2*χ(n-4)T+χ(n-8)T
其中,RA是R峰的振幅,K是通过实验获取的放大因子;
当θ达到阈值时,采用时间导数提取所述时间段ECG数据的峰值后,通过时间序列顺序, 根据峰值前后最小曲率半径位置确定T波起始位置和偏移量,得到T波。
进一步的,所述S5包括:
S51:对所述归一化样本数据特征进行标记,得到跌倒样本数据和日常活动样本数据;
S52:采用主成分分析算法对所述跌倒样本数据和日常活动样本数据进行降维得到降维数 据特征;
所述主成分分析算法:
Y=PX
其中,Y表示降维数据特征,P表示特征向量矩阵,X表示特征矩阵;
S53:采用支持向量机算法对所述降维数据特征进行分类,得到支持向量机分类模型;
给定一个训练数据集D={(xi,yi)|xi∈R,yi∈{-1,1}}
其中,xi是降维后的惯性传感单元和心电图特征向量,每个xi包含k维特征,yi为跌倒类 别,本方案取值为1或-1(对应跌倒与非跌倒)。两类线性分类的超平面可定义为:对于yi=1,wTxi+b≥1,对于yi=-1,wTxi+b≤-1,其中w是加权向量,b是偏差。目标 是通过最大限度减少两个超平面的距离||w||,可表示为一个二次优化问题:
判断函数可以将两个类别(跌倒和非跌倒)分开:
f(xi)=sign(wTxi+b)
其中,f(xi)表示数据特征的类别,xi表示降维数据特征向量,w表示加权向量,b表示偏 差;
对于所有xi,当f(xi)=1时,检测类别为跌倒;当f(xi)=-1时,检测类别为日常动作。
进一步的,所述主成分分析算法步骤:
将所述跌倒样本数据和日常活动样本数据按列组成n行m列的特征矩阵X;
将所述特征矩阵X的每一行元素进行零均值化操作;
计算所述特征矩阵X中不同特征之间的协方差得到协方差矩阵,所述协方差计算公式:
其中表示跌倒样本数据和日常活动样本中特征数据Xi,Xj在第k个跌倒样本数 据和日常活动样本数据中的取值;而则表示在所有跌倒样本数据和日常活动样本中特征数 据Xi的均值,则表示在所有跌倒样本数据和日常活动样本中特征数据Xj的均值,n表示每 个跌倒样本数据和日常活动样本数据的特征维度;
计算所述协方差矩阵的特征值及对应的特征向量;
将所述特征向量按对应特征值大小从上到下按行排列成矩阵,取前k行组成特征向量矩 阵P;
根据不同维度的跌倒样本数据和日常活动样本数据进行跌倒检测的分类,确定k的取值, 得到降维数据特征:
Y=PX
其中,Y表示降维数据特征,P表示特征向量矩阵,X表示特征矩阵。
进一步的,所述惯性传感单元包括三轴加速度计和三轴螺旋仪。
进一步的,所述跌倒数据和日常动作数据包括ECG数据、加速度数据和螺旋仪数据。
一种结合心电和惯性传感数据的跌倒检测系统,包括:
采集模块,包括ECG和惯性传感单元,用于采集人体日常活动过程的跌倒数据和日常动 作数据;
数据处理模块,用于提取分割后的所述跌倒数据和日常动作数据的特征得到惯性传感数 据特征和ECG数据特征,并对提取后的得到惯性传感数据特征和ECG数据特征进行归一化 处理;
模型构建模块,用于对所述归一化样本数据特征进行降维处理得到降维数据特征,并采 用支持向量机算法对所述降维数据特征进行分类,得到支持向量机分类模型;
跌倒检测模块,将所述腕部可穿戴设备新采集的日常活动数据输入所述支持向量机分类 模型,得到跌倒检测结果。
进一步的,所述惯性传感数据特征包括:算术平均值、标准偏差、中位数绝对偏差、最 大值、最小值、频率信号偏斜、频率信号峰度、最大频率成分、平均能量;
所述ECG数据包括QRS复合波、P波和T波数据特征。
本发明与现有技术相比,具有如下的优点和有益效果:
本发明一种结合心电和惯性传感数据的跌倒检测方法及系统,结合ECG和惯性传感数据, 分别从人体生理活动和姿态两个不同维度反映人体跌倒过程的变化,信息更加全面,避免单 一类型数据的缺陷,提高跌倒检测精确度。
附图说明
此处所说明的附图用来提供对本发明实施例的进一步理解,构成本申请的一部分,并不 构成对本发明实施例的限定。在附图中:
图1为本发明整体流程示意图;
图2为本发明另一检测流程图;
图3为本发明ECG数据示意图;
图4为本发明系统结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施例和附图,对本发明 作进一步的详细说明,本发明的示意性实施方式及其说明仅用于解释本发明,并不作为对本 发明的限定。
实施例1
如图1、图2所示,一种结合心电和惯性传感数据的跌倒检测方法,包括:
S1:利用包括ECG和惯性传感单元的腕部可穿戴设备采集人体日常活动过程的跌倒数据 和日常动作数据;
S2:将收集的所述跌倒数据和日常动作数据按相同时间段进行分割;
S3:分别提取分割后的所述跌倒数据和日常动作数据的特征得到惯性传感数据特征和 ECG数据特征;
S4:对所述惯性传感数据特征和ECG数据特征进行归一化处理得到归一化样本数据特征;
S5:采用主成分分析算法对所述归一化样本数据特征进行降维处理得到降维数据特征, 采用支持向量机算法对所述降维数据特征进行分类,得到支持向量机分类模型;
S6:将所述腕部可穿戴设备新采集的日常活动数据输入所述支持向量机分类模型,得到 跌倒检测结果。
本发明的日常活动包含跌倒及正常日常活动。其中跌倒动作在软垫上进行模拟,包括向 前滑倒、向后滑倒、侧向滑倒、绊倒和晕倒;正常日常活动包括走、跑、上下楼梯、姿态转 换、弯腰等。
进一步的,所述S3包括:
S31:提取惯性传感单元中分割后的所述跌倒数据和日常动作数据的惯性传感数据得到惯 性传感数据特征;
所述惯性传感数据特征包括:算术平均值、标准偏差、中位数绝对偏差、最大值、最小 值、频率信号偏斜、频率信号峰度、最大频率成分、平均能量;
惯性传感单元包含三轴加速度计和三轴螺旋仪两种传感单元,人体日常活动过程的跌倒 数据和日常动作数据包括加速度和螺旋仪数据两种,其中每种数据提取九种特征;
因此共提取2*3*9=54种惯性传感数据特征,其中3表示加速度和螺旋仪的角速度均有三 轴,所述惯性传感数据需要对加速度和螺旋仪的角速度每一轴的数据进行提取,故需要乘3。
S32:提取分割后的所述跌倒数据和日常动作数据的ECG数据中的QRS复合波、P波和 T波数据特征得到所述ECG数据特征;
ECG信号的一个正常周期由QRS波、T波和P波组成,如图3所示。P波对应左右心房去极化阶段,频谱范围10-15Hz。QRS复合波对应于左心室和右心室去极化,频谱范围10-40Hz。 T波反映心室复极化过程,在QRS复合波后300毫秒左右的范围。
所述S4中的归一化处理函数:
其中max表示样本数据最大值,min表示样本数据最小值。x表示原始传感器数据,x* 表示归一化结果。
进一步的,所述QRS复合波的提取:
根据ECG数据的峰值提取QRS复合波的峰值点,在峰值点前后分别提取最低点,确定 QRS复合波的起始位置和偏移量,得到QRS复合波。
进一步的,使用以下时间差方程进行滤波处理得到ynT:
ynT=-2*x(n-2)T-χ(n-1)T+χ(n+1)T+2*χ(n+2)T
其中,χnT表示时间T上大小为n的数据样本;
所述P波的提取公式:
θ=0.30*(μn-μHF)
其中,θ表示某一时间段ECG数据的动态阈值表达式,μn指的是在最近一段信号时间序 列内所有数据点的斜率平均值,μHF指一次跳动中的高频噪音;
μHF指一次跳动中的高频噪音,为了计算高频噪音,首先将信号传递到高通滤波器:
y′nT=χnT-2*χ(n-1)T+χ(n-2)T
其中,RA是R峰的振幅,K是通过实验获取的放大因子,本实施例设置为40;
当θ达到阈值时,提取所述时间段ECG数据的峰值后,通过时间序列顺序,根据峰值前 后波形数据的斜率,提取斜率最低值点,确定P波起始位置和偏移量,得到P波。
进一步的,所述T波提取公式:θ=0.30*(μn-μ′HF)
其中,θ表示某一时间段ECG数据的动态阈值表达式,μn指的是在最近一段信号时间序 列内所有数据点的斜率平均值,μ′HF指一次跳动中的低频噪音;
μ′HF指一次跳动中的低频噪音,为了计算低频噪音,首先通过2阶低通滤波器进行平滑处 理:y″nT=2*y(n-1)T-y(n-2)T+χnT-2*χ(n-4)T+χ(n-8)T
其中,RA是R峰的振幅,K是通过实验获取的放大因子;
当θ达到阈值时,采用时间导数提取所述时间段ECG数据的峰值后,通过时间序列顺序, 根据峰值前后最小曲率半径位置确定T波起始位置和偏移量,得到T波。
本实施例通过上述步骤从ECG数据中提取14种特征,如表1所示;
表1 ECG信号特征
序号 | 特征 |
1 | P波持续时间(毫秒) |
2 | QRS波持续时间(毫秒) |
3 | Q波持续时间(毫秒) |
4 | R波持续时间(毫秒) |
5 | S波持续时间(毫秒) |
6 | T波持续时间(毫秒) |
7 | QRS波偏转(毫秒) |
8 | P波波峰振幅(μV) |
9 | QRS波峰幅度(μV) |
10 | Q波振幅(μV) |
11 | R波振幅(μV) |
12 | S波振幅(μV) |
13 | T波波峰振幅(μV) |
14 | QRS波面积(μV) |
所述S5包括:
S51:对所述归一化样本数据特征进行标记,得到跌倒样本数据和日常活动样本数据;
S52:采用主成分分析算法对所述跌倒样本数据和日常活动样本数据进行降维得到降维数 据特征;
所述主成分分析算法:
Y=PX
其中,Y表示降维数据特征,P表示特征向量矩阵,X表示特征矩阵;
S53:采用支持向量机算法对所述降维数据特征进行分类,得到支持向量机分类模型;
给定一个训练数据集D={(xi,yi)|xi∈R,yi∈{-1,1}}
其中,xi是降维后的惯性传感单元和心电图特征向量,每个xi包含k维特征,yi为跌倒类 别,本方案取值为1或-1(对应跌倒与非跌倒)。两类线性分类的超平面可定义为:对于yi=1,wTxi+b≥1,对于yi=-1,wTxi+b≤-1,其中w是加权向量,b是偏差。目标 是通过最大限度减少两个超平面的距离||w||,可表示为一个二次优化问题:
判断函数可以将两个类别(跌倒和非跌倒)分开:
f(xi)=sign(wTxi+b)
其中,f(xi)表示数据特征的类别,xi表示降维数据特征向量,w表示加权向量,b表示偏 差;
对于所有xi,当f(xi)=1时,检测类别为跌倒;当f(xi)=-1时,检测类别为日常动作。
进一步的,所述主成分分析算法步骤:
设有m条68维数据,将所述跌倒样本数据和日常活动样本数据按列组成68行m列的特 征矩阵X;
将所述特征矩阵X的每一行元素进行零均值化操作;
计算所述特征矩阵X中不同特征之间的协方差得到协方差矩阵,所述协方差计算公式:
其中表示跌倒样本数据和日常活动样本中特征数据Xi,Xj在第k个跌倒样本数 据和日常活动样本数据中的取值;而则表示在所有跌倒样本数据和日常活动样本中特征数 据Xi的均值,则表示在所有跌倒样本数据和日常活动样本中特征数据Xj的均值,n表示每 个跌倒样本数据和日常活动样本数据的特征维度;
计算所述协方差矩阵的特征值及对应的特征向量;
将所述特征向量按对应特征值大小从上到下按行排列成矩阵,取前k行组成特征向量矩 阵P;
根据不同维度的跌倒样本数据和日常活动样本数据进行跌倒检测的分类,确定k的取值, 得到降维数据特征:
Y=PX
其中,Y表示降维数据特征,P表示特征向量矩阵,X表示特征矩阵。
进一步的,所述惯性传感单元包括三轴加速度计和三轴螺旋仪。
进一步的,所述跌倒数据和日常动作数据包括ECG数据、加速度数据和螺旋仪数据。
如图4所示,一种结合心电和惯性传感数据的跌倒检测系统,包括:
采集模块,包括ECG和惯性传感单元,用于采集人体日常活动过程的跌倒数据和日常动 作数据;
数据处理模块,用于提取分割后的所述跌倒数据和日常动作数据的特征得到惯性传感数 据特征和ECG数据特征,并对提取后的得到惯性传感数据特征和ECG数据特征进行归一化 处理;
模型构建模块,用于对所述归一化样本数据特征进行降维处理得到降维数据特征,并采 用支持向量机算法对所述降维数据特征进行分类,得到支持向量机分类模型;
跌倒检测模块,将所述腕部可穿戴设备新采集的日常活动数据输入所述支持向量机分类 模型,得到跌倒检测结果。
进一步的,所述惯性传感数据特征包括:算术平均值、标准偏差、中位数绝对偏差、最 大值、最小值、频率信号偏斜、频率信号峰度、最大频率成分、平均能量;
所述ECG数据包括QRS复合波、P波和T波数据特征。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说 明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护 范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本 发明的保护范围之内。
Claims (10)
1.一种结合心电和惯性传感数据的跌倒检测方法,其特征在于,包括:
S1:利用包括ECG和惯性传感单元的腕部可穿戴设备采集人体日常活动过程的跌倒数据和日常动作数据;
S2:将收集的所述跌倒数据和日常动作数据按相同时间段进行分割;
S3:分别提取分割后的所述跌倒数据和日常动作数据的特征得到惯性传感数据特征和ECG数据特征;
S4:对所述惯性传感数据特征和ECG数据特征进行归一化处理得到归一化样本数据特征;
S5:采用主成分分析算法对所述归一化样本数据特征进行降维处理得到降维数据特征,采用支持向量机算法对所述降维数据特征进行分类,得到支持向量机分类模型;
S6:将所述腕部可穿戴设备新采集的日常活动数据输入所述支持向量机分类模型,得到跌倒检测结果。
3.根据权利要求2所述的一种结合心电和惯性传感数据的跌倒检测方法,其特征在于,所述QRS复合波的提取:
根据ECG数据的峰值提取QRS复合波的峰值点,在峰值点前后分别提取最低点,确定QRS复合波的起始位置和偏移量,得到QRS复合波。
4.根据权利要求2所述的一种结合心电和惯性传感数据的跌倒检测方法,其特征在于,
所述P波的提取公式:
θ=0.30*(μn-μHF)
其中,θ表示某一时间段ECG数据的动态阈值表达式,μn指的是在最近一段信号时间序列内所有数据点的斜率平均值,μHF指一次跳动中的高频噪音;
当θ达到阈值时,提取所述时间段ECG数据的峰值后,通过时间序列顺序,根据峰值前后波形数据的斜率,提取斜率最低值点,确定P波起始位置和偏移量,得到P波。
5.根据权利要求2所述的一种结合心电和惯性传感数据的跌倒检测方法,其特征在于,
所述T波提取公式:
θ=0.30*(μn-μ′HF)
其中,θ表示某一时间段ECG数据的动态阈值表达式,μn指的是在最近一段信号时间序列内所有数据点的斜率平均值,μ′HF指一次跳动中的低频噪音;
当θ达到阈值时,采用时间导数提取所述时间段ECG数据的峰值后,通过时间序列顺序,根据峰值前后最小曲率半径位置确定T波起始位置和偏移量,得到T波。
6.根据权利要求1所述的一种结合心电和惯性传感数据的跌倒检测方法,其特征在于,所述S5包括:
S51:对所述归一化样本数据特征进行标记,得到跌倒样本数据和日常活动样本数据;
S52:采用主成分分析算法对所述跌倒样本数据和日常活动样本数据进行降维得到降维数据特征;
所述主成分分析算法:
Y=PX
其中,Y表示降维数据特征,P表示特征向量矩阵,X表示特征矩阵;
S53:采用支持向量机算法对所述降维数据特征进行分类,得到支持向量机分类模型;
所述支持向量机算法:
f(xi)=sign(wTxi+b)
其中,f(xi)表示数据特征的类别,xi表示降维数据特征向量,w表示加权向量,b表示偏差。
7.根据权利要求6所述的一种结合心电和惯性传感数据的跌倒检测方法,其特征在于,所述主成分分析算法步骤:
将所述跌倒样本数据和日常活动样本数据按列组成n行m列的特征矩阵X;
将所述特征矩阵X的每一行元素进行零均值化操作;
计算所述特征矩阵X中不同特征之间的协方差得到协方差矩阵,所述协方差计算公式:
其中Xi (k),Xj (k)表示跌倒样本数据和日常活动样本中特征数据Xi,Xj在第k个跌倒样本数据和日常活动样本数据中的取值,而则表示在所有跌倒样本数据和日常活动样本中特征数据Xi的均值,则表示在所有跌倒样本数据和日常活动样本中特征数据Xj的均值,n表示每个跌倒样本数据和日常活动样本数据的特征维度;
计算所述协方差矩阵的特征值及对应的特征向量;
将所述特征向量按对应特征值大小从上到下按行排列成矩阵,取前k行组成特征向量矩阵P;
根据不同维度的跌倒样本数据和日常活动样本数据进行跌倒检测的分类,确定k的取值,得到降维数据特征:
Y=PX
其中,Y表示降维数据特征,P表示特征向量矩阵,X表示特征矩阵。
8.根据权利要求1所述的一种结合心电和惯性传感数据的跌倒检测方法,其特征在于,所述惯性传感单元包括三轴加速度计和三轴螺旋仪。
9.根据权利要求1所述的一种结合心电和惯性传感数据的跌倒检测方法,其特征在于,所述跌倒数据和日常动作数据包括ECG数据、加速度数据和螺旋仪数据。
10.一种结合心电和惯性传感数据的跌倒检测系统,其特征在于,包括:
采集模块,包括ECG和惯性传感单元,用于采集人体日常活动过程的跌倒数据和日常动作数据;
数据处理模块,用于提取分割后的所述跌倒数据和日常动作数据的特征得到惯性传感数据特征和ECG数据特征,并对提取后的得到惯性传感数据特征和ECG数据特征进行归一化处理;
模型构建模块,用于对所述归一化样本数据特征进行降维处理得到降维数据特征,并采用支持向量机算法对所述降维数据特征进行分类,得到支持向量机分类模型;
跌倒检测模块,将所述腕部可穿戴设备新采集的日常活动数据输入所述支持向量机分类模型,得到跌倒检测结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010537555.5A CN111657918B (zh) | 2020-06-12 | 2020-06-12 | 一种结合心电和惯性传感数据的跌倒检测方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010537555.5A CN111657918B (zh) | 2020-06-12 | 2020-06-12 | 一种结合心电和惯性传感数据的跌倒检测方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111657918A true CN111657918A (zh) | 2020-09-15 |
CN111657918B CN111657918B (zh) | 2021-09-24 |
Family
ID=72387602
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010537555.5A Expired - Fee Related CN111657918B (zh) | 2020-06-12 | 2020-06-12 | 一种结合心电和惯性传感数据的跌倒检测方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111657918B (zh) |
Citations (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4184493A (en) * | 1975-09-30 | 1980-01-22 | Mieczyslaw Mirowski | Circuit for monitoring a heart and for effecting cardioversion of a needy heart |
US6070097A (en) * | 1998-12-30 | 2000-05-30 | General Electric Company | Method for generating a gating signal for cardiac MRI |
US7218960B1 (en) * | 2003-06-24 | 2007-05-15 | Pacesetter, Inc. | System and method for detecting cardiac ischemia based on T-waves using an implantable medical device |
CN104771178A (zh) * | 2015-04-13 | 2015-07-15 | 深圳市飞马与星月科技研究有限公司 | 身份识别的方法及装置 |
CN105796112A (zh) * | 2016-05-11 | 2016-07-27 | 江苏财经职业技术学院 | 老人跌倒检测监护系统 |
CN105869354A (zh) * | 2016-04-29 | 2016-08-17 | 华南理工大学 | 一种基于吸引子传播算法的老人跌倒检测方法 |
CN106407996A (zh) * | 2016-06-30 | 2017-02-15 | 华南理工大学 | 一种基于机器学习的老人跌倒检测方法及其检测系统 |
CN106725495A (zh) * | 2017-01-13 | 2017-05-31 | 深圳先进技术研究院 | 一种跌倒检测方法、装置及系统 |
US9943244B2 (en) * | 2007-09-14 | 2018-04-17 | Angel Medical Systems, Inc. | Waveform feature value averaging system and methods for the detection of cardiac events |
CN108245172A (zh) * | 2018-01-10 | 2018-07-06 | 山东大学 | 一种不受位置约束的人体姿态识别方法 |
CN108257352A (zh) * | 2017-12-30 | 2018-07-06 | 广州柏颐信息科技有限公司 | 一种基于智能穿戴设备的跌倒检测预警方法 |
CN108354610A (zh) * | 2017-08-29 | 2018-08-03 | 浙江好络维医疗技术有限公司 | 一种基于三轴传感器和心电传感器的老人跌倒检测方法及检测系统 |
CN108498091A (zh) * | 2018-05-08 | 2018-09-07 | 复旦大学附属中山医院 | 一种两电极心电图导联智能定位系统 |
CN109171738A (zh) * | 2018-07-13 | 2019-01-11 | 杭州电子科技大学 | 基于人体加速度多特征融合和knn的跌倒检测方法 |
CN109691994A (zh) * | 2019-01-31 | 2019-04-30 | 英菲泰克(天津)科技有限公司 | 一种基于心电图的心率监测分析方法 |
CN110353694A (zh) * | 2019-07-18 | 2019-10-22 | 哈尔滨工业大学 | 一种基于特征选取的运动识别方法 |
CN110659595A (zh) * | 2019-09-10 | 2020-01-07 | 电子科技大学 | 一种基于特征分类的跌倒类型及伤害部位检测方法 |
CN110742599A (zh) * | 2019-11-01 | 2020-02-04 | 广东工业大学 | 一种心电信号特征提取分类方法及系统 |
CN110946573A (zh) * | 2019-11-01 | 2020-04-03 | 东软集团股份有限公司 | 心脏停搏的检测装置、检测模型的训练装置、方法及设备 |
CN111199187A (zh) * | 2019-12-11 | 2020-05-26 | 中国科学院计算机网络信息中心 | 一种基于算法的动物行为识别方法及相应存储介质与电子装置 |
US20200196898A1 (en) * | 2018-12-20 | 2020-06-25 | Queen's University At Kingston | Long QT Syndrome Diagnosis and Classification |
CN111603170A (zh) * | 2020-06-01 | 2020-09-01 | 河南华南医电科技有限公司 | 一种基于心电向量图的人体体位检测方法 |
-
2020
- 2020-06-12 CN CN202010537555.5A patent/CN111657918B/zh not_active Expired - Fee Related
Patent Citations (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4184493A (en) * | 1975-09-30 | 1980-01-22 | Mieczyslaw Mirowski | Circuit for monitoring a heart and for effecting cardioversion of a needy heart |
US6070097A (en) * | 1998-12-30 | 2000-05-30 | General Electric Company | Method for generating a gating signal for cardiac MRI |
US7218960B1 (en) * | 2003-06-24 | 2007-05-15 | Pacesetter, Inc. | System and method for detecting cardiac ischemia based on T-waves using an implantable medical device |
US9943244B2 (en) * | 2007-09-14 | 2018-04-17 | Angel Medical Systems, Inc. | Waveform feature value averaging system and methods for the detection of cardiac events |
CN104771178A (zh) * | 2015-04-13 | 2015-07-15 | 深圳市飞马与星月科技研究有限公司 | 身份识别的方法及装置 |
CN105869354A (zh) * | 2016-04-29 | 2016-08-17 | 华南理工大学 | 一种基于吸引子传播算法的老人跌倒检测方法 |
CN105796112A (zh) * | 2016-05-11 | 2016-07-27 | 江苏财经职业技术学院 | 老人跌倒检测监护系统 |
CN106407996A (zh) * | 2016-06-30 | 2017-02-15 | 华南理工大学 | 一种基于机器学习的老人跌倒检测方法及其检测系统 |
CN106725495A (zh) * | 2017-01-13 | 2017-05-31 | 深圳先进技术研究院 | 一种跌倒检测方法、装置及系统 |
CN108354610A (zh) * | 2017-08-29 | 2018-08-03 | 浙江好络维医疗技术有限公司 | 一种基于三轴传感器和心电传感器的老人跌倒检测方法及检测系统 |
CN108257352A (zh) * | 2017-12-30 | 2018-07-06 | 广州柏颐信息科技有限公司 | 一种基于智能穿戴设备的跌倒检测预警方法 |
CN108245172A (zh) * | 2018-01-10 | 2018-07-06 | 山东大学 | 一种不受位置约束的人体姿态识别方法 |
CN108498091A (zh) * | 2018-05-08 | 2018-09-07 | 复旦大学附属中山医院 | 一种两电极心电图导联智能定位系统 |
CN109171738A (zh) * | 2018-07-13 | 2019-01-11 | 杭州电子科技大学 | 基于人体加速度多特征融合和knn的跌倒检测方法 |
US20200196898A1 (en) * | 2018-12-20 | 2020-06-25 | Queen's University At Kingston | Long QT Syndrome Diagnosis and Classification |
CN109691994A (zh) * | 2019-01-31 | 2019-04-30 | 英菲泰克(天津)科技有限公司 | 一种基于心电图的心率监测分析方法 |
CN110353694A (zh) * | 2019-07-18 | 2019-10-22 | 哈尔滨工业大学 | 一种基于特征选取的运动识别方法 |
CN110659595A (zh) * | 2019-09-10 | 2020-01-07 | 电子科技大学 | 一种基于特征分类的跌倒类型及伤害部位检测方法 |
CN110742599A (zh) * | 2019-11-01 | 2020-02-04 | 广东工业大学 | 一种心电信号特征提取分类方法及系统 |
CN110946573A (zh) * | 2019-11-01 | 2020-04-03 | 东软集团股份有限公司 | 心脏停搏的检测装置、检测模型的训练装置、方法及设备 |
CN111199187A (zh) * | 2019-12-11 | 2020-05-26 | 中国科学院计算机网络信息中心 | 一种基于算法的动物行为识别方法及相应存储介质与电子装置 |
CN111603170A (zh) * | 2020-06-01 | 2020-09-01 | 河南华南医电科技有限公司 | 一种基于心电向量图的人体体位检测方法 |
Non-Patent Citations (13)
Title |
---|
(GONZALEZ-FERNANDEZ, RENE);(RIVERO-VARONA, MARTHA);: "Detection of P Wave in Electrocardiogram", 《COMPUTING IN CARDIOLOGY CONFERENCE》 * |
CHATTERJEE, H.K.; MITRA, M.; GUPTA, R.: "Real-time detection of electrocardiogram wave features using template matching and implementation in FPGA", 《INTERNATIONAL JOURNAL OF BIOMEDICAL ENGINEERING AND TECHNOLOGY》 * |
LI, HQ (LI, HONGQIANG) ; WANG, XF (WANG, XIAOFEI);: "Denoising and R-Peak Detection of Electrocardiogram Signal Based on EMD and Improved Approximate Envelope", 《CIRCUITS SYSTEMS AND SIGNAL PROCESSING》 * |
R. CASTALDO,L. PECCHIA: "Preliminary Results from a Proof of Concept Study for Fall Detection via ECG Morphology", 《IFMBE PROCEEDINGS》 * |
TAN, KF (TAN, KF); CHAN, KL (CHAN, KL); CHOI, K (CHOI, K): "DETECTION OF THE QRS COMPLEX, P WAVE AND T WAVE IN", 《IEE CONFERENCE PUBLICATIONS》 * |
ZHU, JIAJING;LIU, YONGGUO;YANG, CHANGHONG;YANG, WEN;: "A degree-based block model and a local expansion optimization algorithm for anti-community detection in networks", 《PLOS ONE》 * |
尹登峰.: "动态ECG分析中QRS波检测算法的研究", 《中国优秀硕士学位论文全文数据库》 * |
李泽君: "基于心电与脉搏信息的心血管功能检测与监护系统的研究", 《中国优秀硕士学位论文全文数据库》 * |
杨尚明: "盲信号分离ICA理论与应用", 《中国博士学位论文全文数据库》 * |
王莉,郭晓东,惠延波,何宇宸: "基于小波变换的 QRS 波特征提取算法研究与实现", 《专题研究与综述》 * |
陈日清,吴剑,韩永贵.: "基于双侧累计面积的心电信号P波检测", 《中国医药导刊》 * |
陈毛毛: "远程心电监护系统P波自动检测算法的研究与实现", 《中国优秀硕士学位论文全文数据库》 * |
马金伟,刘盛平.: "心电信号识别分类算法综述", 《重庆理工大学学报(自然科学)》 * |
Also Published As
Publication number | Publication date |
---|---|
CN111657918B (zh) | 2021-09-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3065628B1 (en) | Biomechanical activity monitoring | |
Tapia et al. | Real-time recognition of physical activities and their intensities using wireless accelerometers and a heart rate monitor | |
Khan et al. | Accelerometer's position free human activity recognition using a hierarchical recognition model | |
Khan et al. | Accelerometer signal-based human activity recognition using augmented autoregressive model coefficients and artificial neural nets | |
CN102138789B (zh) | 一种动态心电和运动记录与分析系统 | |
CN107137071A (zh) | 一种分析心冲击信号用来计算短期心率值的方法 | |
CN111248879B (zh) | 一种基于多模态注意力融合的高血压老年人活动分析方法 | |
Anastasopoulou et al. | Classification of human physical activity and energy expenditure estimation by accelerometry and barometry | |
Krassnig et al. | User-friendly system for recognition of activities with an accelerometer | |
Jia et al. | Human daily activity recognition by fusing accelerometer and multi-lead ECG data | |
Gjoreski et al. | Telehealth using ECG sensor and accelerometer | |
US20190320944A1 (en) | Biomechanical activity monitoring | |
Fan et al. | Wearable motion attitude detection and data analysis based on Internet of Things | |
Zhang et al. | Upper limb motion capturing and classification for unsupervised stroke rehabilitation | |
Jiang et al. | Exploration of gait parameters affecting the accuracy of force myography-based gait phase detection | |
CN105796091B (zh) | 一种用于去除心电信号车辆运动噪声的智能终端 | |
CN112263244A (zh) | 基于步态的疲劳度评估系统及方法 | |
CN114469041A (zh) | 一种运动过程中心率变化数据特征分析方法 | |
Erin et al. | Spectral Analysis of Cardiogenic Vibrations to Distinguish Between Valvular Heart Diseases. | |
CN111657918B (zh) | 一种结合心电和惯性传感数据的跌倒检测方法及系统 | |
CN116999057A (zh) | 基于穿戴式传感器的偏瘫步态识别及偏瘫步态评估方法 | |
Pawar et al. | Impact of ambulation in wearable-ECG | |
CN116746913A (zh) | 一种基于压力传感器监测步态和身体健康的设备 | |
Nouredanesh et al. | Detection of compensatory balance responses using wearable electromyography sensors for fall-risk assessment | |
CN111297366B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210924 |
|
CF01 | Termination of patent right due to non-payment of annual fee |