CN103177530B - 一种基于三轴加速度的人体摔倒检测方法 - Google Patents

一种基于三轴加速度的人体摔倒检测方法 Download PDF

Info

Publication number
CN103177530B
CN103177530B CN201310059523.9A CN201310059523A CN103177530B CN 103177530 B CN103177530 B CN 103177530B CN 201310059523 A CN201310059523 A CN 201310059523A CN 103177530 B CN103177530 B CN 103177530B
Authority
CN
China
Prior art keywords
accekeration
tstart
timpact
axis
tmax
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
Application number
CN201310059523.9A
Other languages
English (en)
Other versions
CN103177530A (zh
Inventor
陈辰
黄晁
王鸿龙
张从连
任梁
候开宇
丁华锋
杨洪
刘晓金
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
NINGBO ZHONGKE IC DESIGN CENTER CO Ltd
Original Assignee
NINGBO ZHONGKE IC DESIGN CENTER CO Ltd
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 NINGBO ZHONGKE IC DESIGN CENTER CO Ltd filed Critical NINGBO ZHONGKE IC DESIGN CENTER CO Ltd
Priority to CN201310059523.9A priority Critical patent/CN103177530B/zh
Publication of CN103177530A publication Critical patent/CN103177530A/zh
Application granted granted Critical
Publication of CN103177530B publication Critical patent/CN103177530B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

本发明公开了一种基于三轴加速度的人体摔倒检测方法,其首先获取由前2秒对应的所有加速度值和当前1秒的所有加速度值构成的加速度值集合,然后根据该加速度值集合确定用于表示疑似摔倒的加速度值序列,接着根据该加速度值集合获取用于表示疑似摔倒的加速度值序列的特征向量,再根据特征向量计算摔倒判定值,最后根据摔倒判定值确定人体是否摔倒,以及根据当前获取的所有加速度值随后的报警判定时间内的任意一个采集时间采集的加速度值,确定摔倒的人体是否失去自救能力,这种摔倒检测方法,不仅实现简单,而且计算复杂度低、检测精度高。

Description

一种基于三轴加速度的人体摔倒检测方法
技术领域
本发明涉及一种人体摔倒检测技术,尤其是涉及一种基于三轴加速度的人体摔倒检测方法。
背景技术
摔倒是日常生活中对人造成意外伤害的主要原因之一。当摔倒者是老年人时,情况往往会更为严重。据民政部门对老年人养老意愿的调查,九成以上的老年人仍然希望在家养老,而传统的儿女照料老年人的模式现在已经难以应对不断增长的养老需求,这就造成了大量独居老人,由于没有亲人的照料,独居老人在遇到意外摔倒时,往往得不到及时有效的救助。
日常生活中,导致摔倒的原因主要有两种:一是行动中意外跌倒或滑倒;二是心脑血管等疾病引起的摔倒。对于意外摔伤,如果没有及时地治疗,则往往会导致伤情加重。对于心肌梗死、脑梗塞突然发作摔倒时,如果几个小时内不能得到及时有效的治疗,则患者甚至会付出生命的代价,而此类心脑血管疾病造成的休克、昏迷等状况通常也会让患者失去主动求救的能力。
目前,人体摔倒判定的技术主要有三类:基于视频和声音的摔倒检测技术、基于压力的摔倒检测技术及基于穿戴式传感器的摔倒检测技术。其中,基于视频和声音的摔倒检测技术是利用摄像头以及其它声音采集工具对一定区域进行监控,再通过图像处理中人体识别和姿势识别来检测监控图像中的人体摔倒事件;基于压力的摔倒检测技术是利用铺设于地板上的压力传感器,通过物体对地面的压力变化来确定撞击等事件,从而检测摔倒。这两种摔倒检测技术都只能检测一定区域内的摔倒事件,无法对特定的对象进行持续的监控。而基于穿戴式传感器的摔倒检测技术是通过佩戴在人体上的传感器采集人体活动数据,从而检测目标的摔倒事件。
在基于穿戴式传感器的摔倒检测技术中,目前应用较多的摔倒检测方法主要有判断阈值的方法和利用模式识别的方法,前者简单有效,但识别正确率很难达到实际应用的标准;后者在正确率上比较不错,但检测过程过于复杂,不适用于单片机的运行,并且在正确率方面也尚有提升空间。
发明内容
本发明所要解决的技术问题是提供一种实现简单,且计算复杂度低、检测精度高的基于三轴加速度的人体摔倒检测方法。
本发明解决上述技术问题所采用的技术方案为:一种基于三轴加速度的人体摔倒检测方法,其特征在于包括以下步骤:
①采用佩戴于人体上的三轴加速度传感器,按设定的采集频率实时采集人体的三轴加速度数据,将采集到的人体的X轴加速度值、Y轴加速度值和Z轴加速度值对应记为Ax、Ay和Az,然后根据Ax、Ay和Az,计算人体的加速度值,记为A, A = A x 2 + A y 2 + A z 2 ;
②当获得第3秒对应的所有加速度值时,将第3秒对应的所有加速度值作为当前加速度值序列;
③将当前加速度值序列的前2秒对应的所有加速度值和当前加速度值序列中的所有加速度值按采集时间顺序构成的集合定义为当前加速度值集合,然后根据当前加速度值集合,确定用于表示疑似摔倒的加速度值序列;
④根据当前加速度值集合,获取用于表示疑似摔倒的加速度值序列的特征向量,记为X;
⑤根据X计算摔倒判定值,记为f(X),然后判断f(X)是否大于0,如果是,则判定人体已摔倒,再执行步骤⑥,否则,判定人体未摔倒,再执行步骤⑦,其中,exp()表示以自然对数e为底的指数函数,e=2.718281828459,“||||”为欧氏距离运算符号,w、b和γ均为判定参数;
⑥判断当前加速度值序列随后的报警判定时间内的任意一个采集时间采集的加速度值是否大于重力加速度G的0.75倍且小于重力加速度G的1.5倍,如果是,则确定人体摔倒严重,需报警,然后执行步骤⑦,否则,确定人体摔倒不严重,无需报警,然后执行步骤⑦;
⑦将下一秒获得的所有加速度值作为当前加速度值序列,然后返回步骤③继续执行。
所述的设定的采集频率为大于或等于100HZ。
所述的步骤③中用于表示疑似摔倒的加速度值序列的获取过程为:
③-1、从当前加速度值集合中找出值最大的加速度值,并记为A(Tmax)将A(Tmax)的采集时间记为Tmax,其中,Ax(Tmax)表示A(Tmax)中的X轴加速度值,Ay(Tmax)表示A(Tmax)中的Y轴加速度值,Az(Tmax)表示A(Tmax)中的Z轴加速度值;
③-2、在当前加速度值集合中寻找第一个同时满足下列条件的加速度值,条件一:要求该加速度值大于重力加速度值G的0.75倍;条件二:要求该加速度值的采集时间的下一个采集时间采集的加速度值小于或等于重力加速度值G的0.75倍;条件三:要求Tmax与该加速度值的采集时间的差小于1秒;如果存在,则将该加速度值确定为人体失重开始的加速度值,记为A(Tstart) A ( Tstart ) = A x ( Tstart ) 2 + A y ( Tstart ) 2 + A z ( Tstart ) 2 , 将A(Tstart)的采集时间记为Tstart,其中,A(Tstart)>0.75×G,A(Tstart+1)≤0.75×G,A(Tstart+1)表示Tstart的下一个采集时间采集的加速度值,Ax(Tstart)表示A(Tstart)中的X轴加速度值,Ay(Tstart)表示A(Tstart)中的Y轴加速度值,Az(Tsatrt)表示A(Tstart)中的Z轴加速度值,Tmax-Tstart<1秒,然后执行步骤③-3;如果不存在,则执行步骤⑦;
③-3、计算A(Tstart)对应的方向与A(Tmax)对应的方向之间的夹角,记为θ, &theta; = arccos A x ( T max ) &times; A x ( Tstart ) + A y ( T max ) &times; A y ( Tstart ) + A z ( T max ) &times; A z ( Tstart ) A ( T max ) &times; A ( Tstart ) , 然后判断θ是否大于或等于45°,如果是,则将Tstart到Tmax之间采集的所有加速度值确定为用于表示疑似摔倒的加速度值序列,然后执行步骤④;否则,执行步骤⑦。
所述的步骤④的具体过程为:
④-1、在当前加速度值集合中寻找唯一一个同时满足下列条件的加速度值,条件一:要求Tstart到该加速度值的采集时间之间的任意一个采集时间采集的加速度值小于或等于重力加速度值G的0.75倍;条件二:要求该加速度值大于重力加速度值G的0.75倍;然后将该加速度值确定为人体失重结束的加速度值,记为A(Tend) 将A(Tend)的采集时间记为Tend,其中,A(Tend)>0.75×G,Ax(Tend)表示A(Tend)中的X轴加速度值,Ay(Tend)表示A(Tend)中的Y轴加速度值,Az(Tend)表示A(Tend)中的Z轴加速度值,对于任何满足Tstart<T<Tend的采集时间T,有A(T)≤0.75×G,A(T)表示采集时间为T时采集的加速度值;
④-2、在当前加速度值集合中寻找唯一一个同时满足下列条件的加速度值,条件一:要求该加速度值的采集时间到Tmax之间的任意一个采集时间采集的加速度值大于或等于重力加速度值G的1.2倍;条件二:要求该加速度值小于重力加速度值G的1.2倍;然后将该加速度值确定为人体撞击的加速度值,记为A(Timpact) A ( Timpact ) = A x ( Timpact ) 2 + A y ( Timpact ) 2 + A z ( Timpact ) 2 , 将A(Timpact)的采集时间记为Timpact,其中,A(Timpact)<1.2×G,Ax(Timpact)表示A(Timpact)中的X轴加速度值,Ay(Timpact)表示A(Timpact)中的Y轴加速度值,Az(Timpact)表示A(Timpact)中的Z轴加速度值,对于任何满足Timpact<T'<Tmax的采集时间T',有A(T')≥1.2×G,A(T')表示采集时间为T'时采集的加速度值;
④-3、计算用于表示疑似摔倒的加速度值序列的特征向量,记为X,X=[t1,μ1,σ1,t2,μ2,σ2,t3,μ3,σ3],其中,t1=Tend-Tstart,μ1表示Tstart到Tend之间采集的所有加速度值的平均值,σ1表示Tstart到Tend之间采集的所有加速度值的方差,t2=Timpact-Tend,μ2表示Tend到Timpact之间采集的所有加速度值的平均值,σ2表示Tend到Timpact之间采集的所有加速度值的方差,t3=T max-Timpact,μ3表示Timpact到Tmax之间采集的所有加速度值的平均值,σ3表示Timpact到Tmax之间采集的所有加速度值的方差。
所述的步骤⑤中w为向量(-0.975,0.0106383,-1,-0.978022,0.824114,0.206675,0.35,-0.53481,-0.564189),b=1.5,γ=1。
所述的步骤⑥中的报警判定时间为3~10秒。
所述的步骤⑥中的报警判定时间为5秒。
与现有技术相比,本发明的优点在于:
1)本发明方法首先获取由前2秒对应的所有加速度值和当前1秒获取的所有加速度值构成的加速度值集合,然后根据该加速度值集合确定用于表示疑似摔倒的加速度值序列,接着根据该加速度值集合获取用于表示疑似摔倒的加速度值序列的特征向量,再根据特征向量计算摔倒判定值,最后根据摔倒判定值确定人体是否摔倒,以及根据当前获取的所有加速度值随后的报警判定时间内的任意一个采集时间采集的加速度值,确定摔倒的人体是否失去自救能力,这种摔倒检测方法,不仅实现简单,而且计算复杂度低、检测精度高。
2)通过计算人体摔倒前与人体摔倒后的角度差即通过计算一段时间内人体失重开始的加速度值的方向与最大的加速度值的方向之间的角度差,联合失重与超重作为可疑摔倒数据的筛选门限,有效地减少了传统的利用模式识别的摔倒检测方法中多余的计算,使得本发明方法能够方便地在单片机实现,同时极大地减少了误报率。
3)通过对用于表示疑似摔倒的加速度值序列的特征向量的优化,将摔倒过程分为三段,这样能够更合理的提取出摔倒数据的特征,有利于通过计算摔倒判定值更好的区分摔倒和非摔倒,从而有效地提高了摔倒判断的准确率。
4)通过筛选疑似摔倒加速度序列以及利用摔倒判定公式计算摔倒判定值进一步区分摔倒和非摔倒,相比现有的基于判断阈值的摔倒检测方法,极大的提高了摔倒检测的准确率。
5)通过设置报警判定时间和可疑摔倒数据的筛选,最大限度地减少了误报,因此在计算摔倒判定值的参数选择上,更倾向于将一次疑似摔倒数据判断为摔倒,这样在减少误报的同时又减少了漏报。
具体实施方式
以下结合实施例对本发明作进一步详细描述。
本发明提出的一种基于三轴加速度的人体摔倒检测方法,其具体包括以下步骤:
①采用佩戴于人体上身的三轴加速度传感器,按设定的采集频率实时采集人体的三轴加速度数据,将采集到的人体的X轴加速度值、Y轴加速度值和Z轴加速度值对应记为Ax、Ay和Az,然后根据Ax、Ay和Az,计算人体的加速度值,记为A, A = A x 2 + A y 2 + A z 2 .
在本实施例中,设定的采集频率为大于或等于100HZ,即每隔小于或等于10ms对人体的三轴加速度数据进行一次采集,在实际实施过程中,可根据实际需求设置三轴加速度传感器的采集频率,一般情况下,如果人体的三轴加速度数据采集的时间间隔较短,则会有效提高最终的检测精度,但同时会增加计算复杂度,通过实验验证每隔10ms对人体的三轴加速度数据进行一次采集,能够很好地提高检测的综合性能。
②当获得第3秒对应的所有加速度值时,将第3秒对应的所有加速度值作为当前加速度值序列。
③将当前加速度值序列的前2秒对应的所有加速度值和当前加速度值序列中的所有加速度值按采集时间顺序构成的集合定义为当前加速度值集合,然后根据当前加速度值集合,确定用于表示疑似摔倒的加速度值序列。
在此具体实施例中,步骤③中用于表示疑似摔倒的加速度值序列的获取过程为:
③-1、从当前加速度值集合中找出值最大的加速度值,并记为A(Tmax)将A(Tmax)的采集时间记为Tmax,其中,Ax(Tmax)表示A(Tmax)中的X轴加速度值,Ay(Tmax)表示A(Tmax)中的Y轴加速度值,Az(Tmax)表示A(Tmax)中的Z轴加速度值。
③-2、在当前加速度值集合中寻找第一个同时满足下列条件的加速度值,条件一:要求该加速度值大于重力加速度值G的0.75倍;条件二:要求该加速度值的采集时间的下一个采集时间采集的加速度值小于或等于重力加速度值G的0.75倍;条件三:要求Tmax与该加速度值的采集时间的差小于1秒;如果存在,则将该加速度值确定为人体失重开始的加速度值,记为A(Tstart) A ( Tstart ) = A x ( Tstart ) 2 + A y ( Tstart ) 2 + A z ( Tstart ) 2 , 将A(Tstart)的采集时间记为Tstart,其中,A(Tstart)>0.75×G,A(Tstart+1)≤0.75×G,A(Tstart+1)表示Tstart的下一个采集时间采集的加速度值,Ax(Tstart)表示A(Tstart)中的X轴加速度值,Ay(Tstart)表示A(Tstart)中的Y轴加速度值,Az(Tstart)表示A(Tstart)中的Z轴加速度值,Tmax-Tstart<1秒,然后执行步骤③-3;如果不存在,则执行步骤⑦。
在此,如果三轴加速度传感器的采集频率为100HZ,则采集时间为10ms。
③-3、计算A(Tstart)对应的方向与A(Tmax)对应的方向之间的夹角,记为θ, &theta; = arccos A x ( T max ) &times; A x ( Tstart ) + A y ( T max ) &times; A y ( Tstart ) + A z ( T max ) &times; A z ( Tstart ) A ( T max ) &times; A ( Tstart ) , 然后判断θ是否大于或等于45°,如果是,则将Tstart到Tmax之间采集的所有加速度值确定为用于表示疑似摔倒的加速度值序列,然后执行步骤④;否则,执行步骤⑦。
④根据当前加速度值集合,获取用于表示疑似摔倒的加速度值序列的特征向量,记为X。
在此具体实施例中,步骤④的具体过程为:
④-1、在当前加速度值集合中寻找唯一一个同时满足下列条件的加速度值,条件一:要求Tstart到该加速度值的采集时间之间的任意一个采集时间采集的加速度值小于或等于重力加速度值G的0.75倍;条件二:要求该加速度值大于重力加速度值G的0.75倍;然后将该加速度值确定为人体失重结束的加速度值,记为A(Tend)将A(Tend)的采集时间记为Tend,其中,A(Tend)>0.75×G,Ax(Tend)表示A(Tend)中的X轴加速度值,Ay(Tend)表示A(Tend)中的Y轴加速度值,Az(Tend)表示A(Tend)中的Z轴加速度值,对于任何满足Tstart<T<Tend的采集时间T,有A(T)≤0.75×G,A(T)表示采集时间为T时采集的加速度值。
④-2、在当前加速度值集合中寻找唯一一个同时满足下列条件的加速度值,条件一:要求该加速度值的采集时间到Tmax之间的任意一个采集时间采集的加速度值大于或等于重力加速度值G的1.2倍;条件二:要求该加速度值小于重力加速度值G的1.2倍;然后将该加速度值确定为人体撞击的加速度值,记为A(Timpact) A ( Timpact ) = A x ( Timpact ) 2 + A y ( Timpact ) 2 + A z ( Timpact ) 2 , 将A(Timpact)的采集时间记为Timpact,其中,A(Timpact)<1.2×G,Ax(Timpact)表示A(Timpact)中的X轴加速度值,Ay(Timpact)表示A(Timpact)中的Y轴加速度值,Az(Timpact)表示A(Timpact)中的Z轴加速度值,对于任何满足Timpact<T'<Tmax的采集时间T',有A(T)≥1.2×G,A(T')表示采集时间为T'时采集的加速度值。
④-3、计算用于表示疑似摔倒的加速度值序列的特征向量,记为X,X=[t1,μ1,σ1,t2,μ2,σ2,t3,μ3,σ3],其中,t1=Tend-Tstart,μ1表示Tstart到Tend之间采集的所有加速度值的平均值,σ1表示Tstart到Tend之间采集的所有加速度值的方差,t2=Timpact-Tend,μ2表示Tend到Timpact之间采集的所有加速度值的平均值,σ2表示Tend到Timpact之间采集的所有加速度值的方差,t3=T max-Timpact,μ3表示Timpact到Tmax之间采集的所有加速度值的平均值,σ3表示Timpact到Tmax之间采集的所有加速度值的方差。
⑤根据X计算摔倒判定值,记为f(X),然后判断f(X)是否大于0,如果是,则判定人体已摔倒,再执行步骤⑥,否则,判定人体未摔倒,再执行步骤⑦,其中,exp()表示以自然对数e为底的指数函数,e=2.718281828459,“||||”为欧氏距离运算符号,w、b和γ均为判定参数,在本实施例中,w是一个向量,取值为向量(-0.975,0.0106383,-1,-0.978022,0.824114,0.206675,0.35,-0.53481,-0.564189),b=1.5,γ=1。
⑥判断当前加速度值序列随后的报警判定时间内的任意一个采集时间采集的加速度值是否大于重力加速度G的0.75倍且小于重力加速度G的1.5倍,如果是,则确定人体摔倒严重,需报警,然后执行步骤⑦,否则,确定人体摔倒不严重,无需报警,然后执行步骤⑦。
在本实施例中,报警判定时间可设为3~10秒,在实际操作过程中,可根据实际情况设定,如可将报警判定时间设为5秒,通常人体摔倒后5秒能够确定人体摔倒后是否失去自救能力。
⑦将下一秒获得的所有加速度值作为当前加速度值序列,然后返回步骤③继续执行。
为更有力地说明本发明方法的可行性和有效性,对本发明方法的误报率和漏报率性能进行测试。
误报率是指目标没有摔倒,却将其动作误判成摔倒;漏报率是指目标摔倒后没有成功检测到的比率,漏报率要求尽可能的低,这是因为对老年人来说,误报最多是虚惊一场,而漏报却可能造成严重的后果。
对于漏报率的测试,要求20位测试人员每人模拟30次摔倒,并且记录本发明方法的检测结果。原则上并不规定测试人员摔倒的方向和姿势,但是要求测试人员每一次摔倒的方向或姿势各不相同。在测试过程中通过冲击力来判断这一次摔倒是易受伤摔倒还是不易受伤摔倒,并对两种情况下的漏报做出统计,统计结果如表1所列。
表1不易受伤摔倒和易受伤摔倒两种情况下的漏报率统计结果
漏报次数 测试次数 漏报率
不易受伤摔倒 19 236 8.05%
易受伤摔倒 0 364 0%
所有摔倒 19 600 3.17%
从表1中可以看出,对于易受伤的摔倒,本发明方法的漏报率为0,而对于不易受伤的摔倒,本发明方法的漏报率为8.05%。
由于本发明方法在步骤⑤计算摔倒判定值的过程中采用了高维度的向量运算,因此比普通标量运算更能够精确的检测摔倒。本发明方法利用计算摔倒判定值f(X),更倾向于将一次疑似摔倒判定为摔倒,因而漏报率比传统的利用模式识别的摔倒检测技术有明显的降低。
表1中,易受伤摔倒包括站立时突然摔倒、运动中摔倒等对人体伤害较大的摔倒情况,本发明方法中通过w、b和γ的合理取值,可以完全消除对于易受伤摔倒的漏报。
表1中,不易受伤摔倒包括重心较低时(如蹲下的状态或坐在比较低的凳子上)的摔倒,这一类摔倒通常不会对身体造成伤害。此类摔倒后,人体目标基本无需外部救助。本发明方法中对此类摔倒的漏报率为8.05%。
现有的基于判断阈值的摔倒检测技术对不同方向的摔倒漏报率往往不同,本发明方法中,由于只需计算摔倒前与摔倒后的角度差,对摔倒的方向并不敏感,消除了不同方向摔倒的差异。
对于误报率的测试,要求20位测试人员佩戴应用了本发明方法的报警器,佩戴方式是将报警器固定在测试人员腰部,每人佩戴一周,每天5小时以上。佩戴期间正常进行工作生活,对误报原因及次数进行统计,结果如表2所列。
表2不同误报原因下的误报次数统计结果
误报原因 误报次数
随意抛掷或掉落报警器 17
用力坐下 2
户外运动及其他 0
从表2中可以看出,平均每人每天发生误报0.14次。其中大多数是因为测试人员将报警器随意抛掷在桌子上或是掉在地上,少部分是因为测试人员比较猛烈的坐倒在沙发上。由于所选的测试人员的年龄都在20岁到35岁之间,而本发明方法的应用对象都是老年人,老年人正常情况下不会像年轻人这样动作剧烈。因此,测试中的绝大多数误报可以避免。
本发明方法侧重于消除运动时的误报,在现有的利用模式识别的摔倒检测技术中,由于运动时的一些身体特征和摔倒类似,因此在运动中的误报较多。
本发明方法通过合理选取特征向量,因此能够很好的区分运动(如跑、跳、上楼梯等)和摔倒,误报率较现有的利用模式识别的摔倒检测技术有明显提高。

Claims (7)

1.一种基于三轴加速度的人体摔倒检测方法,其特征在于包括以下步骤:
①采用佩戴于人体上的三轴加速度传感器,按设定的采集频率实时采集人体的三轴加速度数据,将采集到的人体的X轴加速度值、Y轴加速度值和Z轴加速度值对应记为Ax、Ay和Az,然后根据Ax、Ay和Az,计算人体的加速度值,记为A, A = A x 2 + A y 2 + A z 2 ;
②当获得第3秒对应的所有加速度值时,将第3秒对应的所有加速度值作为当前加速度值序列;
③将当前加速度值序列的前2秒对应的所有加速度值和当前加速度值序列中的所有加速度值按采集时间顺序构成的集合定义为当前加速度值集合,然后根据当前加速度值集合,确定用于表示疑似摔倒的加速度值序列;
④根据当前加速度值集合,获取用于表示疑似摔倒的加速度值序列的特征向量,记为X;
⑤根据X计算摔倒判定值,记为f(X),然后判断f(X)是否大于0,如果是,则判定人体已摔倒,再执行步骤⑥,否则,判定人体未摔倒,再执行步骤⑦,其中,exp()表示以自然对数e为底的指数函数,e=2.718281828459,“||||”为欧氏距离运算符号,w、b和γ均为判定参数;
⑥判断当前加速度值序列随后的报警判定时间内的任意一个采集时间采集的加速度值是否大于重力加速度G的0.75倍且小于重力加速度G的1.5倍,如果是,则确定人体摔倒严重,需报警,然后执行步骤⑦,否则,确定人体摔倒不严重,无需报警,然后执行步骤⑦;
⑦将下一秒获得的所有加速度值作为当前加速度值序列,然后返回步骤③继续执行。
2.根据权利要求1所述的一种基于三轴加速度的人体摔倒检测方法,其特征在于所述的设定的采集频率为大于或等于100HZ。
3.根据权利要求1或2所述的一种基于三轴加速度的人体摔倒检测方法,其特征在于所述的步骤③中用于表示疑似摔倒的加速度值序列的获取过程为:
③-1、从当前加速度值集合中找出值最大的加速度值,并记为A(Tmax)将A(Tmax)的采集时间记为Tmax,其中,Ax(Tmax)表示A(Tmax)中的X轴加速度值,Ay(Tmax)表示A(Tmax)中的Y轴加速度值,Az(Tmax)表示A(Tmax)中的Z轴加速度值;
③-2、在当前加速度值集合中寻找第一个同时满足下列条件的加速度值,条件一:要求该加速度值大于重力加速度值G的0.75倍;条件二:要求该加速度值的采集时间的下一个采集时间采集的加速度值小于或等于重力加速度值G的0.75倍;条件三:要求Tmax与该加速度值的采集时间的差小于1秒;如果存在,则将该加速度值确定为人体失重开始的加速度值,记为A(Tstart)将A(Tstart)的采集时间记为Tstart,其中,A(Tstart)>0.75×G,A(Tstart+1)≤0.75×G,A(Tstart+1)表示Tstart的下一个采集时间采集的加速度值,Ax(Tstart)表示A(Tstart)中的X轴加速度值,Ay(Tstart)表示A(Tstart)中的Y轴加速度值,Az(Tstart)表示A(Tstart)中的Z轴加速度值,Tmax-Tstart<1秒,然后执行步骤③-3;如果不存在,则执行步骤⑦;
③-3、计算A(Tstart)对应的方向与A(Tmax)对应的方向之间的夹角,记为θ,
&theta; = arccos A x ( T max ) &times; A x ( Tstart ) + A y ( T max ) &times; A y ( Tstart ) + A z ( T max ) &times; A z ( Tstart ) A ( T max ) &times; A ( Tstart ) , 然后判断θ是否大于或等于45°,如果是,则将Tstart到Tmax之间采集的所有加速度值确定为用于表示疑似摔倒的加速度值序列,然后执行步骤④;否则,执行步骤⑦。
4.根据权利要求3所述的一种基于三轴加速度的人体摔倒检测方法,其特征在于所述的步骤④的具体过程为:
④-1、在当前加速度值集合中寻找唯一一个同时满足下列条件的加速度值,条件一:要求Tstart到该加速度值的采集时间之间的任意一个采集时间采集的加速度值小于或等于重力加速度值G的0.75倍;条件二:要求该加速度值大于重力加速度值G的0.75倍;然后将该加速度值确定为人体失重结束的加速度值,记为A(Tend)将A(Tend)的采集时间记为Tend,其中,A(Tend)>0.75×G,Ax(Tend)表示A(Tend)中的X轴加速度值,Ay(Tend)表示A(Tend)中的Y轴加速度值,Az(Tend)表示A(Tend)中的Z轴加速度值,对于任何满足Tstart<T<Tend的采集时间T,有A(T)≤0.75×G,A(T)表示采集时间为T时采集的加速度值;
④-2、在当前加速度值集合中寻找唯一一个同时满足下列条件的加速度值,条件一:要求该加速度值的采集时间到Tmax之间的任意一个采集时间采集的加速度值大于或等于重力加速度值G的1.2倍;条件二:要求该加速度值小于重力加速度值G的1.2倍;然后将该加速度值确定为人体撞击的加速度值,记为A(Timpact) A ( Timpact ) = A x ( Timpact ) 2 + A y ( Timpact ) 2 + A z ( Timpact ) 2 , 将A(Timpact)的采集时间记为Timpact,其中,A(Timpact)<1.2×G,Ax(Timpact)表示A(Timpact)中的X轴加速度值,Ay(Timpact)表示A(Timpact)中的Y轴加速度值,Az(Timpact)表示A(Timpact)中的Z轴加速度值,对于任何满足Timpact<T'<Tmax的采集时间T',有A(T')≥1.2×G,A(T')表示采集时间为T'时采集的加速度值;
④-3、计算用于表示疑似摔倒的加速度值序列的特征向量,记为X,X=[t1,μ1,σ1,t2,μ2,σ2,t3,μ3,σ3],其中,t1=Tend-Tstart,μ1表示Tstart到Tend之间采集的所有加速度值的平均值,σ1表示Tstart到Tend之间采集的所有加速度值的方差,t2=Timpact-Tend,μ2表示Tend到Timpact之间采集的所有加速度值的平均值,σ2表示Tend到Timpact之间采集的所有加速度值的方差,t3=Tmax-Timpact,μ3表示Timpact到Tmax之间采集的所有加速度值的平均值,σ3表示Timpact到Tmax之间采集的所有加速度值的方差。
5.根据权利要求4所述的一种基于三轴加速度的人体摔倒检测方法,其特征在于所述的步骤⑤中w为向量(-0.975,0.0106383,-1,-0.978022,0.824114,0.206675,0.35,-0.53481,-0.564189),b=1.5,γ=1。
6.根据权利要求5所述的一种基于三轴加速度的人体摔倒检测方法,其特征在于所述的步骤⑥中的报警判定时间为3~10秒。
7.根据权利要求6所述的一种基于三轴加速度的人体摔倒检测方法,其特征在于所述的步骤⑥中的报警判定时间为5秒。
CN201310059523.9A 2013-02-25 2013-02-25 一种基于三轴加速度的人体摔倒检测方法 Active CN103177530B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310059523.9A CN103177530B (zh) 2013-02-25 2013-02-25 一种基于三轴加速度的人体摔倒检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310059523.9A CN103177530B (zh) 2013-02-25 2013-02-25 一种基于三轴加速度的人体摔倒检测方法

Publications (2)

Publication Number Publication Date
CN103177530A CN103177530A (zh) 2013-06-26
CN103177530B true CN103177530B (zh) 2015-07-08

Family

ID=48637352

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310059523.9A Active CN103177530B (zh) 2013-02-25 2013-02-25 一种基于三轴加速度的人体摔倒检测方法

Country Status (1)

Country Link
CN (1) CN103177530B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103581852B (zh) * 2013-09-30 2018-03-06 吴家宝 人体摔倒检测的方法、装置及移动终端系统
CN103577836B (zh) 2013-09-30 2018-01-23 吴家宝 人体摔倒检测模型建立方法及模型系统
CN105894730A (zh) * 2014-12-26 2016-08-24 陆婷 一种便携式gsm&gps跌倒检测器的研究与设计
CN108257352B (zh) * 2017-12-30 2019-12-20 广州柏颐信息科技有限公司 一种基于智能穿戴设备的跌倒检测预警方法
CN110021144B (zh) * 2018-01-10 2022-08-19 厦门雅迅网络股份有限公司 一种摔倒检测方法、终端设备及存储介质
CN110638459B (zh) * 2019-09-03 2023-05-05 宁波路晟电器科技有限公司 一种基于加速度传感器的人体运动跌倒检测装置及方法
CN112258785B (zh) * 2019-11-27 2022-06-07 北京京东尚科信息技术有限公司 一种跌倒检测方法、装置及存储介质

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2323196A (en) * 1997-03-11 1998-09-16 Keith Henderson Cameron Automatic fall alarm
CN101950464A (zh) * 2010-09-17 2011-01-19 中国科学院深圳先进技术研究院 跌倒监测与报警的方法和系统

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7145461B2 (en) * 2001-01-31 2006-12-05 Ilife Solutions, Inc. System and method for analyzing activity of a body
JP2008032521A (ja) * 2006-07-28 2008-02-14 Icom Inc 転倒検出装置、転倒検出方法及びコンピュータプログラム

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2323196A (en) * 1997-03-11 1998-09-16 Keith Henderson Cameron Automatic fall alarm
CN101950464A (zh) * 2010-09-17 2011-01-19 中国科学院深圳先进技术研究院 跌倒监测与报警的方法和系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于加速度传感器的老年人跌倒检测装置设计;李冬等;《传感器与微系统》;20081231;第27卷(第9期);第85-88页 *

Also Published As

Publication number Publication date
CN103177530A (zh) 2013-06-26

Similar Documents

Publication Publication Date Title
CN103177530B (zh) 一种基于三轴加速度的人体摔倒检测方法
Hildebrand et al. Evaluation of raw acceleration sedentary thresholds in children and adults
Bianchi et al. Barometric pressure and triaxial accelerometry-based falls event detection
CN105125221B (zh) 云服务实时摔倒检测系统及方法
JP5674766B2 (ja) 着用位置を検出するためのセンシングデバイス
Kangas et al. Sensitivity and false alarm rate of a fall sensor in long-term fall detection in the elderly
US9877667B2 (en) Method for quantifying the risk of falling of an elderly adult using an instrumented version of the FTSS test
Bourke et al. A threshold-based fall-detection algorithm using a bi-axial gyroscope sensor
Prosperini et al. The clinical relevance of force platform measures in multiple sclerosis: a review
Majumder et al. iPrevention: Towards a novel real-time smartphone-based fall prevention system
Hashim et al. Accurate fall detection for patients with Parkinson's disease based on a data event algorithm and wireless sensor nodes
AT507941B1 (de) Verfahren zur detektion einer aussergewöhnlichen situation
CN103927851B (zh) 一种个人化多阈值跌倒检测方法及系统
CN110047248A (zh) 一种基于可穿戴的老人跌倒检测方法及设备
CN106456024A (zh) 静息心率监测器系统
TW201136569A (en) Apparatus for identifying falls and activities of daily living
CN106618487A (zh) 一种老年人平衡能力评定方法
Zhao et al. Recognition of human fall events based on single tri-axial gyroscope
CN108354610A (zh) 一种基于三轴传感器和心电传感器的老人跌倒检测方法及检测系统
CN106971502B (zh) 一种基于多传感器联合的快速跌倒检测系统及方法
CN110313917A (zh) 以日常生活动作为判断基准的跌倒传感系统及其方法
CN107146378A (zh) 一种人体跌倒判定方法及装置
Liu et al. Validation of an accelerometer-based fall prediction model
Locatelli et al. Differentiating essential tremor and Parkinson's disease using a wearable sensor—A pilot study
US11410521B2 (en) Fall detector incorporating physiological sensing

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant