CN102333483B - 用于检测人的步行的系统和方法 - Google Patents
用于检测人的步行的系统和方法 Download PDFInfo
- Publication number
- CN102333483B CN102333483B CN201080009404.7A CN201080009404A CN102333483B CN 102333483 B CN102333483 B CN 102333483B CN 201080009404 A CN201080009404 A CN 201080009404A CN 102333483 B CN102333483 B CN 102333483B
- Authority
- CN
- China
- Prior art keywords
- dominant frequency
- axis
- signal
- measurement axis
- measurement
- 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
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/1036—Measuring load distribution, e.g. podologic studies
- A61B5/1038—Measuring plantar pressure during gait
-
- 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/1113—Local tracking of patients, e.g. in a hospital or private home
- A61B5/1114—Tracking parts of the body
-
- 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/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/1123—Discriminating type of movement, e.g. walking or running
-
- 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/683—Means for maintaining contact with the body
- A61B5/6831—Straps, bands or harnesses
-
- 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/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/10—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
- G01C21/12—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
- G01C21/16—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
- G01C21/165—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation combined with non-inertial navigation instruments
- G01C21/1654—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation combined with non-inertial navigation instruments with electromagnetic compass
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C22/00—Measuring distance traversed on the ground by vehicles, persons, animals or other moving solid bodies, e.g. using odometers, using pedometers
- G01C22/006—Pedometers
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- 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
Abstract
一种用于检测人的步行的系统,布置有包括至少双轴运动传感器(CM)的壳体(BT),该壳体适合于附着到所述人的身体的上部,从而所述传感器(CM)的第一测量轴适于提供代表所述身体的前后轴(AP)或垂直轴(VT)的度量,并且,从而所述传感器(CM)的第二测量轴适于提供代表所述身体的中侧轴(ML)的度量,所述系统还布置有用于分析由该传感器(CM)递送的测量结果的分析装置(MA),其中,该分析装置(MA)包括:处理装置(MT),其用于在时间窗口上处理由所述传感器(CM)递送的测量信号,所述处理装置包括用于搜索所述信号中的主频(MRFD)的装置,和检测装置(MD),其用于在第一测量轴的信号的主频和第二测量轴的主频之间的比或由所述传感器(CM)传输的测量结果的向量的欧几里得范数的主频和第二测量轴的信号的主频之间的比大体上等于二时,检测人的步行。
Description
技术领域
本发明涉及用于检测人的步行的系统和方法,或者,换句话说,用于通过包括一系列步伐的运动的模式而检测人的运动的系统和方法。
背景技术
用于分析人的运动的系统在生物医学领域中日趋广泛,尤其是用于分析人的身体活动。
对人的行走活动的检测所提供的信息使得例如估计人的能量消耗、评价人的久坐水平或估计在外科手术或药物处置之后的身体机能的质量或丢失成为可能。
Najafi,B.、Aminian,K.、Paraschiv-Ionescu A.、Loew,F.、Bula C.J.以及Robert,P.的标题为“Ambulatory system for human motion analysis using akinematic sensor:monitoring of daily physical activity in the elderly,”Biomedical Engineering,IEEE Transactions on,2003年6月、第50卷、no.6、第711-723页的文献描述2A1G运动传感器(双轴加速度计和单轴陀螺测试仪),该运动传感器佩戴在人的躯干上,并且,由0.62-5.00Hz带通滤波器对垂直加速度信号进行滤波。在经滤波的信号上寻找到至少三个高于阈值的等间隔发生的峰值振幅。难以先验地设置通用阈值,这显著地涉及这样的系统的可靠性的缺乏。
Bernard Auvinet、Gilles Berrut、Claude Touzard、Laurent Moutel、NadineCollet、Denis Chaleil以及Eric Barrey的标题为“Reference data for normalsubjects obtained with an accelerometric device”,Gait & Posture、2002年10月、第16卷、第2期、第124-134页的文献描述对被认为是大体上周期性活动的行走活动的频率分析,其以取决于行走速度的频率产生功率峰值。对垂直加速度信号的偶次谐波和奇次谐波之间的比的研究使研究步行的稳定性成为可能。这并不涉及检测行走活动,而涉及在已知人在行走时对人的行走活动进行分析或描绘特性。
M.Sekine、T.Tamura、T.Togawa以及Y.Fukui的标题为“Classificationof waist-acceleration signals in a continuous walking record”,MedicalEngineering & Physics 22(4)(2000),第285-291页的文献描述使用小波变换而在代表人的步行的信号中区分人是否在水平表面上行走、爬楼梯或下楼。该文献的内容并不使检测行走活动成为可能。
发明内容
本发明的目的是,在人的走动信号的记录中检测该人的行走活动。
根据本发明的一个方面,提出了一种用于检测人的步行的系统,该系统布置有包括双轴或三轴运动传感器的壳体。该壳体适合于附着到所述人的身体的上部,从而所述传感器的第一测量轴适于提供代表所述身体的前后轴或垂直轴的度量,并且,从而所述传感器的第二测量轴适于提供代表所述身体的中侧轴的度量,所述系统还布置有用于分析由所述传感器递送的测量结果的分析装置。所述分析装置包括:
-处理装置,其用于在时间窗口上处理由所述传感器递送的测量信号,所述处理装置包括用于搜索所述信号中的主频的装置,和
-检测装置,其用于在第一测量轴的信号的主频和第二测量轴的主频之间或由所述传感器传输的测量结果的向量的欧几里得范数的主频和第二测量轴的信号的主频之间的比大体上等于二时,检测人的步行。
这样的系统使以降低的成本且在用户佩戴该系统时几乎没有不适感的情况下,以鲁棒且自动的方式检测人的步行成为可能。
例如,时间窗口是滑动窗口。
因此,该系统甚至在相当大的处理期期间都是极其精确的。
在一个实施例中,所述运动传感器是三轴的,所述传感器的第一测量轴与所述身体的前后轴重合,所述传感器的第二测量轴与所述身体的中侧轴重合,并且,所述传感器的第三测量轴与所述身体的垂直轴重合,所述检测装置适合于检测第一测量轴的信号的主频和第二测量轴的主频之间或第三测量轴的信号的主频和第二测量轴的主频之间或由所述传感器传输的测量结果的向量的欧几里得范数的主频和第二测量轴的信号的主频之间的大体上等于二的比。
因此,改进检测的精度。
根据一个实施例,该系统还包括高通滤波器。
因此,去除由运动传感器传输的信号的相应的连续分量,以便能够以高精度检测主频。
在一个实施例中,该系统还包括例如具有0.5和10Hz的频带的带通滤波器。
因此,大大地制约与步行无关的信噪或频率的影响。
根据一个实施例,所述分析装置在所述壳体的内部或外部,并且,所述运动传感器包括用于将其测量结果传输至该分析装置的有线或无线的传输装置。
分析装置可以并入壳体中或安装在远程库上,并且,来自壳体的输出信号无论是否被分析,都能够利用线或不利用线传输。
所述运动传感器可以是双轴或三轴加速度计、双轴或三轴磁强计或双轴或三轴陀螺测试仪。
本发明利用所有这些类型的运动传感器来工作。
例如,滑动时间窗口持续五秒,在偏移一秒的两个连续的窗口之间存在着四秒的部分重叠。
这些值特别适合于人的步行。
根据一个实施方式,用于搜索由运动传感器传输的信号的主频的所述装置适合于在每个时间窗口中通过谱分析而执行对主频的搜索。例如,该谱分析能够是谱图类型的。
使用与切趾窗口卷积的信号的傅里叶变换模的平方的谱图是一种搜索主频的简单、可靠且低成本的方法,所述主频即与最大信号功率相对应的频率。
用于搜索主频的所述装置可以适合于将对沿着第二轴vML的主频的搜索限制在0.25Hz和1Hz之间的频率。
用于搜索主频的所述装置可以适合于在第一轴与前后轴重合时,将对沿着第一轴的主频的搜索限制在预定的频率范围中的频率。该频率范围能够以沿着第二轴的主频fML Hz增加0.2Hz和3Hz为边界([fML+0.2;3])。
优选,该范围可以以沿着第二轴的主频fML增加0.25Hz和2Hz为边界([fML+0.25;2])。
用于搜索主频的所述装置可以适合于在第一轴与垂直轴重合时,将对沿着第一轴的主频的搜索限制在预定的频率范围中所包括的频率。该频率范围可以以沿着第二轴的主频fML Hz增加0.2Hz和3Hz为边界([fML+0.2;3])。
用于搜索主频的所述装置可以适合于将对针对由所述运动传感器传输的测量结果的向量的欧几里得范数的主频的搜索限制在预定的频率范围中所包括的频率。该频率范围可以以沿着第二轴的主频fML Hz增加0.2Hz和3Hz为边界([fML+0.2;3])。
所有这些值特别适合于步行活动。
根据一个实施例,所述检测装置适合于当在至少一个轴上的至少一个频率下的信号功率也高于阈值时,检测方差内的大体上等于二的所述主频的比。
换句话说,仅针对在至少一个轴上的至少一个频率下具有的功率高于所确定的阈值的窗口而确定或使用这样的比,该阈值能够被称为功率阈值。例如在测试阶段期间先验地确定或实验地调整该功率阈值。
不但可能将至少一个频率的功率的条件应用于至少一个轴,而且还应用于所有轴。应当指出,当该功率条件应用于各种轴时,每个功率阈值可以互不相同。
根据特定的实施例,所述确定装置适合于当主频在至少一个轴上具有高于阈值功率的功率时,确定与给定的时间窗口相对应的主频的比。
可能将该阈值条件仅应用于主频或者同样地应用于其他定义的频率或频带。
应当指出,该功率阈值标准还可以应用于测量向量的欧几里得范数。对于每个时间窗口,例如主频的(多个)频率被验证为具有高于功率阈值的功率。
功率阈值可以先验地或通过在例如解剖校准的未发生任何事情的时隙中预先确定而设置。
根据一个实施例,所述壳体适合于附着到所述人的躯干或骶骨。
当壳体附着到躯干时,躯干的振动幅度较高,这改进系统的精度。
对骶骨的附着特别容易且谨慎,例如借助于带。
根据本发明的另一方面,也提出一种用于基于由双轴或三轴运动传感器对沿着所述传感器的第一测量轴和沿着所述传感器的第二测量轴的运动进行的测量而检测人的步行的方法,第一测量轴适于提供代表所述人的身体的前后轴或垂直轴的度量,第二测量轴适于提供代表所述身体的中侧轴的度量,其中:
-在时间窗口上处理由所述运动传感器递送的测量信号,所述处理包括搜索所述信号中的主频,以及
-在第一测量轴的信号的主频和第二测量轴的主频之间或由所述传感器传输的测量结果的向量的欧几里得范数的主频和第二测量轴的信号的主频之间的比大体上等于二时,检测所述人的步行。
例如,在滑动时间窗口上执行处理。
附图说明
在研究作为非限制性的示例而描述并由附图图解说明的许多实施例的基础上,将更好地理解本发明,其中:
-图1示意地图解说明根据本发明的一个方面的系统的一个实施例;
-图2示意地表示人及其前后、垂直以及中侧解剖轴;
-图3图解说明根据图1的由系统进行的测量结果的示例,其中,运动传感器是双轴加速度计;
-图4和5图解说明分析装置的运行;
-图6a和6b图解说明根据本发明的一个方面的系统的第一实施例;
-图7a、7b以及7c图解说明根据本发明的一个方面的系统的第二实施例;
-图8a和8b图解说明根据本发明的一个方面的系统的第三实施例;
-图9a和9b图解说明根据本发明的一个方面的系统的第四实施例;以及
-图10a和10b图解说明根据本发明的一个方面的系统的第五实施例。
具体实施方式
在所有附图中,具有相同的标记的元件是相似的。
如在图1中所图解说明的,用于检测人的步行的系统包括壳体BT,该壳体BT包括双轴或三轴运动传感器CM。壳体BT适合于附着到所述人的身体的上部,在该实例中,借助于弹性附着带CEF,从而所述运动传感器的第一测量轴适于提供代表所述身体的前后轴AP或垂直轴VT的度量,并且,所述运动传感器的第二测量轴适于提供代表所述身体的中侧轴ML的度量。作为变型,任何其他附着装置都可以是合适的。
例如,运动传感器的第一测量轴可以与身体的前后轴AP或垂直轴VT重合,并且,运动传感器的第二测量轴可以与身体的中侧轴ML重合。
该重合可以例如通过解剖校准来实现,例如通过要求附着有壳体BT的人靠墙尽可能地站直几秒。该系统以已知的方式确定应用于测量的旋转矩阵,以便给予简化为中侧轴ML、前后轴AP或垂直轴VT的测量结果。运动传感器CM还提供有传输模块MTR,以便将测量结果通过在该示例中为无线传输的传输而传输至外部站SE,在这种情况下,外部站SE为膝上型计算机。
作为变型,传输能够通过线而进行。运动传感器可以是例如双轴或三轴加速度计、双轴或三轴磁强计或者双轴或三轴陀螺测试仪。
然而,在剩余的描述中,以非限制性的方式,运动传感器CM是双轴加速度计,其第一测量轴与人的身体的前后轴AP重合,并且,第二测量轴与人的身体的中侧轴ML重合。
作为变型,加速度计的第二测量轴可以与人的身体的中侧轴ML重合,并且,第一测量轴可以与人的身体的垂直轴VT重合。
膝上型计算机SE包括用于分析由加速度计CM传输的数据的分析模块MA。作为变型,分析模块可以并入壳体BT中。
分析模块适合于在低于或等于1kHz且典型地为10至200Hz的数量级的采样频率下对从加速度计CM接收的信号进行采样。
分析模块MA包括用于处理由加速度计CM递送的测量信号的处理模块MT。
作为变型,在诸如三轴加速度计的三轴运动传感器CM的情况下,有可能执行解剖校准,从而加速度计的第一测量轴与人的身体的前后轴AP重合,加速度计的第二测量轴与人的身体的中侧轴ML重合,并且,加速度计的第三测量轴与人的身体的垂直轴VT重合。在这种情况下,检测模块适合于检测第一测量轴的信号的主频和第二测量轴的主频之间的比或者第三测量轴的信号的主频和第二测量轴的主频之间的比或者由所述传感器传输的测量结果的向量的欧几里得范数和第二测量轴的信号的主频之间的比大体上等于二。于是,这给予改进的检测精度。
大体上等于2的比意指例如1.7和2.3之间的比,并且,优选为1.9和2.1之间。
可以预先确定该比,但同样可以试验地调整该比,尤其是在测试阶段期间。然后,当人在行走时,对该比所取的各种值精确地进行分析,然后,确定在算法中使用的临界值。考虑到与假阳性(当人不在行走时,设备发出他在行走的信号)或假阴性(当人在行走时,设备指示他不在行走)相关联的风险,该确定尤其可以统计的方式进行。
处理模块MT可以包括高通滤波器FPH,这使去除由加速度计CM传输的信号的相应的连续分量成为可能,以便以高精度检测主频。
处理模块MT还可以包括带通滤波器,以便大大地限制与行走无关的信噪或频率的影响。
此外,处理模块MT包括用于通过谱分析而搜索由运动传感器传输的信号的主频MRFD的模块。包含根据频率而估计信号的功率的谱分析是一种用于搜索信号的主频的在计算方面简单且低成本的已知的方法。“主频”被理解为与信号的功率密度的最大值相对应的频率。当然,作为变型,可以设想任何其他搜索主频MRFD的方法。例如,谱分析不但可以通过使用傅里叶变换来执行,而且也可以通过例如小波变换、更适合于非平稳信号的技术的本领域技术人员已知的其他技术来执行。
分析模块MA还包括检测模块MD,该检测模块MD用于在第一测量轴的信号的主频和第二测量轴的主频之间,或者由所述传感器传输的测量结果的向量的欧几里得范数的主频和第二测量轴的信号的主频之间的比大体上等于二时,检测人的步行。
本发明在不需要运动传感器CM的物理校准的情况下工作,或者,换句话说,根据本发明的系统基于以数值单元或伏特表示的原始数据工作,并且,运动传感器CM的增益和位移的知识并非必要。如果决定不将伏特转换为物理单位(例如,对于加速度计而言,m/s2),最小功率阈值的概念可以基于对休息状态下的人的测量而不基于运动学上的数据而决定。
图2示意地图解说明人及其三个解剖轴,即中侧轴ML、前后轴AP以及垂直轴VT,其取向为使得三面体(ML、VT、AP)是直接三面体。身体的中侧轴ML取向为从身体的左部至身体的右部,前后轴AP取向为从身体的后部至身体的前部,并且,垂直轴取向为从身体的上部至身体的下部。
最佳地,壳体可以放在躯干上或骶骨上。
带通滤波器FPB可以例如是在特别适合于行走的0.5和10Hz之间的频带中进行滤波的4阶巴特沃斯滤波器。
信号被分成用于分析的时间窗口,时间窗口优选为滑动时间窗口,例如,被分成五秒的窗口,在偏移一秒的两个连续的窗口之间存在着四秒的部分重叠。然后,对于每个时间窗口,寻找信号的主频。
用于通过对谱图类型的谱分析而搜索主频MRFD的模块可以适于将对沿着第一轴的主频fML的搜索限制在0.25Hz和1Hz之间的频率。
用于通过谱分析而搜索主频MRFD的模块也可以适于在第一轴与前后轴AP重合时将对沿着第一轴的主频的搜索限制在沿着第二轴的主频fML Hz增加0.2Hz和3Hz之间的频率或沿着第二轴的主频fML Hz增加0.25Hz和2Hz之间的频率([fML+0.25;2])。
用于通过谱分析而搜索主频MRFD的模块还可以适于在第一轴与垂直轴VT重合时将对沿着第一轴的主频的搜索限制在沿着第二轴的主频fML Hz增加0.2Hz和3Hz之间的频率。
用于通过谱分析而搜索主频MRFD的模块还可以适于将对针对由所述运动传感器CM传输的测量结果的向量的欧几里得范数的主频的搜索限制在沿着第二轴的主频fML Hz增加0.2Hz和3Hz之间的频率。
所有这些对主频的搜索的限制特别适合于行走,并且,使限制所使用的计算时间和内存大小成为可能。
此外,测试表明,对其中制约主频搜索的这样的频率范围的选择通过显著地降低假阳性或假阴性的风险而使增加可靠性且甚至稳定性成为可能。
图3示出根据时间的由根据本发明的一个方面的壳体BT传输的信号SML和SAP的示例,该壳体BT布置有双轴加速度计CM,其第一测量轴与前后轴AP重合,其第二测量轴与中侧轴ML重合。壳体BT例如放在其活动被监测的人的骶骨上。
图4针对图3的数据的情况而表示与信号SML和SAP相对应的根据时间的主频fML和fAP、由用于对主频MRFD进行搜索的模块通过滑动窗口计算主频fML和fAP。
图5针对图3和4的情况而表示由行走检测模块MD基于主频fML和fAP的比而进行的计算。然后,该系统检测与初始时刻0s相对应的时刻和与初始时刻之后182s相对应的时刻之间的行走活动,以及在与初始时刻之后221s相对应的时刻重新开始的行走活动。
在用户佩戴根据本发明的一个方面的行走检测系统的下列示范性的实施例中,用户在各种测试期间采取不动或移动的各种姿态。
在该实例中,在200Hz的采样率下进行测量。在持续10s的滑动窗口上,将这些数据聚集在一起,其中,两个连续的窗口之间重叠90%。由4阶IIR巴特沃斯滤波器将带通滤波[0.1Hz;10Hz]应用于这些测量结果。通过计算所测量的信号乘以切趾窗口的乘积的傅里叶变换模的平方而执行每个窗口的谱分析,从而使达到谱图成为可能。在每个窗口上确定第二轴上的主频fML。
在这些示例中,对于每个测量轴,在频率值的相应的优选范围中确定主频:fML在0.25Hz和1Hz之间,并且,fAP在fML+0.25和2Hz之间。
图6a和6b图解说明第一示例。图6a图解说明利用在带上佩戴的根据本发明的一个方面的系统来进行的记录。该检测不受沿着中侧轴ML测量的信号的功率的影响而起动。图6a表示根据窗口指数的主频fML和fAP。它们是相应地沿着中侧轴ML和前后轴AP针对每个时间窗口而确定的主频。
图6b表示针对每个时间窗口的主频fAP和fML之间的比。在该示例中,人仅在时间窗口160和210之间行走。
与这些时间窗口相对应的比大约处于值2。因此,该测试使确定阈值成为可能,该阈值大体上等于2,在该示例中,例如,等于1.9,高于该阈值,则认为人正在行走。可以根据该类型的测试或通过已知的统计分析技术而手动地确定该阈值,从而使估计假阳性和假阴性的风险成为可能。
现在,图解说明阈值大体上等于2,就是说,接近2,但并不严格等于2,在实验测试期间能够进行调整。该调整可以是手动的,或者例如通过确定fAP和fML之间的比的统计分布,并且,通过估计该分布的某些参数,例如在假设分布为正态的情况下的平均值和标准偏差,而自动地生成该调整。
图7a、7b以及图7c图解说明在带上携带行走监测系统的第二示例。例如在沿着中侧轴ML测量的信号上强加所测量的信号的功率阈值。在此,功率阈值设置在0.01g2(g=9.81m·s-2),低于该功率阈值,则不试图检测比fAP/fML。在根据窗口指数的功率的图7a中表示该功率。在该实例中,功率单位是重力常数的平方。因此,该曲线表示针对每个时间窗口沿着中侧轴ML确定的主频的功率。
图7b仅包含仅针对具有沿着中侧轴ML测量的信号的时间窗口的沿着前后轴AP的主频fAP,该信号的功率高于所提到的阈值。在该示例中,用户在时间窗口155和210之间以及时间窗口102和107之间行走。
在图7c中,针对每个与用户的步行相对应的窗口,获得大体上等于2的比,就是说,在1.8和2.2之间。在该实例中,阈值能够设置在1.8或1.9。将功率标准应用于沿着中侧轴ML测量的信号或沿着前后轴AP测量的信号,或者同时应用于这些信号,而且,阈值能够不同。
图8a和图8b图解说明第三示例,对于第三示例,图8b表示由所述传感器CM传输的测量结果的向量的欧几里得范数的主频和第二测量轴的信号的主频fML的比。
在图9a、9b和10a、10b中相应地图解说明两个其他示范性的实施例,系统相应地附着到带和躯干。
这些示例示出用户在时间窗口150和200之间行走。在这两个示例中特别清楚地图解说明,无论系统放在身体的上部的什么地方,系统都极其可靠。
本发明使以低成本且以高精度检测人的行走活动成为可能。
尽管本发明特别地描述了行走阶段的检测,但是该检测也能够应用于跑步阶段,则调整用于搜索主频的频率范围。
此外,本发明并不强制物理地校准运动传感器。
Claims (18)
1.一种用于检测人的步行的系统,布置有包括至少双轴运动传感器(CM)的壳体(BT),所述壳体适合于附着到所述人的身体的上部,从而所述传感器(CM)的第一测量轴适于提供代表所述身体的前后轴(AP)或垂直轴(VT)的度量,并且,从而所述传感器(CM)的第二测量轴适于提供代表所述身体的中侧轴(ML)的度量,所述系统还布置有用于分析由所述传感器(CM)递送的测量结果的分析装置(MA),其特征在于,所述分析装置(MA)包括:
-处理装置(MT),其用于在时间窗口上处理由所述传感器(CM)递送的测量信号,所述处理装置包括用于搜索所述信号中的主频(MRFD)的装置,以及
-检测装置(MD),其用于在所述第一测量轴的信号的主频和所述第二测量轴的信号的主频之间或由所述传感器(CM)传输的测量结果的向量的欧几里得范数的主频和所述第二测量轴的信号的主频之间的比大体上等于二时,检测所述人的步行。
2.如权利要求1所述的系统,其中,所述时间窗口是滑动窗口。
3.如权利要求1所述的系统,其中,所述运动传感器(CM)是三轴的,所述传感器(CM)的所述第一测量轴与所述身体的前后轴(AP)重合,所述传感器(CM)的所述第二测量轴与所述身体的中侧轴(ML)重合,并且,所述传感器(CM)的第三测量轴与所述身体的垂直轴重合,所述检测装置(MD)适合于检测所述第一测量轴的信号的主频和所述第二测量轴的信号的主频之间或所述第三测量轴的信号的主频和所述第二测量轴的信号的主频之间或由所述传感器传输的测量结果的向量的欧几里得范数的主频和所述第二测量轴的信号的主频之间的大体上等于二的比。
4.如权利要求1所述的系统,还包括高通滤波器(FPH)。
5.如权利要求1所述的系统,还包括带通滤波器(FPB)。
6.如权利要求1所述的系统,其中,所述分析装置(MA)在所述壳体(BT)的内部或外部,并且其中,所述运动传感器(CM)包括用于将其测量结果传输至所述分析装置(MA)的有线或无线的传输装置(MTR)。
7.如权利要求1所述的系统,其中,所述运动传感器(CM)是双轴或三轴加速度计或双轴或三轴磁强计或双轴或三轴陀螺测试仪。
8.如权利要求2所述的系统,其中,所述滑动窗口为五秒长,并且在偏移一秒的两个连续的窗口之间存在着四秒的部分重叠。
9.如权利要求1所述的系统,其中,用于搜索由所述运动传感器(CM)传输的信号的主频(MRFD)的所述装置适合于通过对谱图类型的谱分析而执行对主频的搜索。
10.如权利要求1所述的系统,其中,用于搜索主频(MRFD)的所述装置适合于将对沿着所述第二测量轴的主频的搜索限制在0.25Hz和1Hz之间的频率。
11.如权利要求1所述的系统,其中,用于搜索主频(MRFD)的所述装置适合于在所述第一测量轴与前后轴重合时,将对沿着所述第一测量轴的主频的搜索限制在沿着所述第二测量轴的主频Hz增加0.2Hz和3Hz之间的频率。
12.如权利要求1所述的系统,其中,用于搜索主频(MRFD)的所述装置适合于在所述第一测量轴与垂直轴重合时,将对沿着所述第一测量轴的主频的搜索限制在沿着所述第二测量轴的主频Hz增加0.2Hz和3Hz之间的频率。
13.如权利要求1所述的系统,其中,用于搜索主频(MEFD)的所述装置适合于将对针对由所述运动传感器传输的测量结果的向量的欧几里得范数的主频的搜索限制在沿着所述第二测量轴的主频Hz增加0.2Hz和3Hz之间的频率。
14.如权利要求1所述的系统,其中,所述检测装置(MD)适合于在相应的信号的功率也高于阈值时,检测在所述第一测量轴的信号的主频和所述第二测量轴的信号的主频之间或由所述传感器(CM)传输的测量结果的向量的欧几里得范数的主频和所述第二测量轴的信号的主频之间的比,所述比的范围在1.7和2.3之间。
15.如权利要求1所述的系统,其中,所述检测装置(MD)适合于在相应的信号的功率也高于阈值时,检测在所述第一测量轴的信号的主频和所述第二测量轴的信号的主频之间或由所述传感器(CM)传输的测量结果的向量的欧几里得范数的主频和所述第二测量轴的信号的主频之间的大体上等于二的比,其中,所述阈值是通过估计检测到假阳性和假阴性的风险而手动地确定的。
16.如权利要求1所述的系统,其中,所述壳体(BT)适合于附着到所述人的躯干或骶骨。
17.如权利要求14所述的系统,其中,所述比的范围在1.9和2.1之间。
18.一种用于基于由双轴或三轴运动传感器对沿着所述传感器(CM)的第一测量轴以及沿着所述传感器(CM)的第二测量轴的运动进行的测量而检测人的步行的方法,所述第一测量轴适于提供代表所述人的身体的前后轴(AP)或垂直轴(VT)的度量,所述第二测量轴适于提供代表所述身体的中侧轴(ML)的度量,其特征在于:
-在时间窗口上处理由所述运动传感器(CM)递送的测量信号,所述处理包括搜索所述信号中的主频,以及
-在所述第一测量轴的信号的主频和所述第二测量轴的信号的主频之间或由所述传感器传输的测量结果的向量的欧几里得范数的主频和所述第二测量轴的信号的主频之间的比大体上等于二时,检测所述人的步行。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR0951212 | 2009-02-26 | ||
FR0951212A FR2942388B1 (fr) | 2009-02-26 | 2009-02-26 | Systeme et procede de detection de marche d'une personne |
PCT/EP2010/052368 WO2010097422A1 (fr) | 2009-02-26 | 2010-02-25 | Systeme et procede de detection de marche d'une personne |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102333483A CN102333483A (zh) | 2012-01-25 |
CN102333483B true CN102333483B (zh) | 2015-07-22 |
Family
ID=41082361
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201080009404.7A Active CN102333483B (zh) | 2009-02-26 | 2010-02-25 | 用于检测人的步行的系统和方法 |
Country Status (7)
Country | Link |
---|---|
US (1) | US9265448B2 (zh) |
EP (1) | EP2400889A1 (zh) |
JP (1) | JP5874130B2 (zh) |
KR (1) | KR101718555B1 (zh) |
CN (1) | CN102333483B (zh) |
FR (1) | FR2942388B1 (zh) |
WO (1) | WO2010097422A1 (zh) |
Families Citing this family (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2972344B1 (fr) | 2011-03-07 | 2014-01-31 | Lape Medical | Dispositif de surveillance d'une prothese medicale et du corps humain |
RU2650586C2 (ru) * | 2011-11-28 | 2018-04-16 | Конинклейке Филипс Н.В. | Система мониторинга состояния здоровья для вычисления общей оценки риска |
FR2984511B1 (fr) * | 2011-12-19 | 2014-12-05 | Commissariat Energie Atomique | Systeme et procede de detection d'au moins une phase transitoire dans une activite stationnaire d'un etre anime |
FR3015072B1 (fr) | 2013-12-18 | 2017-03-17 | Movea | Procede de determination de l'orientation d'un repere capteur lie a un terminal mobile muni d'un ensemble capteur, porte par un utilisateur et comprenant au moins un capteur de mouvement lie en mouvement |
JP6567658B2 (ja) * | 2014-05-30 | 2019-08-28 | 日東電工株式会社 | ユーザーの活動を分類し及び/又はユーザーの歩数をカウントするデバイス及び方法 |
US9877668B1 (en) * | 2014-11-21 | 2018-01-30 | University Of South Florida | Orientation invariant gait matching |
EP3224575B1 (en) * | 2014-11-27 | 2023-08-23 | Razer (Asia-Pacific) Pte. Ltd. | Step counter devices and step counting methods |
EP3032455A1 (en) | 2014-12-09 | 2016-06-15 | Movea | Device and method for the classification and the reclassification of a user activity |
WO2016168610A1 (en) | 2015-04-15 | 2016-10-20 | Nike, Inc. | Activity monitoring device with assessment of exercise intensity |
WO2016196289A1 (en) * | 2015-05-29 | 2016-12-08 | Nike Innovate C.V. | Activity monitoring device with assessment of exercise intensity |
JP2017023689A (ja) * | 2015-07-24 | 2017-02-02 | 株式会社東芝 | モニタリングシステム、モニタ方法およびプログラム |
KR102556924B1 (ko) * | 2016-09-05 | 2023-07-18 | 삼성전자주식회사 | 보행 보조 방법 및 이를 수행하는 장치 |
WO2018081795A1 (en) | 2016-10-31 | 2018-05-03 | Zipline Medical, Inc. | Systems and methods for monitoring physical therapy of the knee and other joints |
GB2574074B (en) | 2018-07-27 | 2020-05-20 | Mclaren Applied Tech Ltd | Time synchronisation |
KR102399672B1 (ko) * | 2019-06-11 | 2022-05-20 | 한국과학기술연구원 | 보행 시간-주파수 분석에 기초한 개인 식별 방법 및 시스템 |
KR102395937B1 (ko) * | 2019-06-11 | 2022-05-11 | 한국과학기술연구원 | 보행 시간-주파수 분석에 기초한 건강 상태 예측 방법 및 시스템 |
GB2588236B (en) | 2019-10-18 | 2024-03-20 | Mclaren Applied Ltd | Gyroscope bias estimation |
KR102336580B1 (ko) * | 2019-10-30 | 2021-12-10 | 한국생산기술연구원 | 좌우 걸음의 보행 균형도 분석방법 |
KR102357770B1 (ko) * | 2020-03-30 | 2022-02-04 | 인제대학교 산학협력단 | 웨어러블 가속도계를 이용한 신경퇴행성 질환 약효 소진 증상 측정 장치 및 방법 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2875320Y (zh) * | 2006-02-28 | 2007-03-07 | 深圳市万机创意电子科技有限公司 | 人体步行传感器 |
CN101327125A (zh) * | 2007-06-23 | 2008-12-24 | 株式会社百利达 | 步行评价系统、步行计、步行评价程序和记录介质 |
Family Cites Families (30)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE19538925C2 (de) * | 1995-10-19 | 2000-07-27 | Wieland Friedmund | Vorrichtung zur Auswertung eines Narkose- oder Intensiv-EEG |
US6776766B2 (en) * | 1996-04-03 | 2004-08-17 | Rush-Presbyterian-St. Luke's Medical Center | Method and apparatus for characterizing gastrointestinal sounds |
US6522266B1 (en) * | 2000-05-17 | 2003-02-18 | Honeywell, Inc. | Navigation system, method and software for foot travel |
EP1195139A1 (en) * | 2000-10-05 | 2002-04-10 | Ecole Polytechnique Féderale de Lausanne (EPFL) | Body movement monitoring system and method |
IL147502A0 (en) * | 2002-01-07 | 2002-08-14 | Widemed Ltd | Self-adaptive system, for the analysis of biomedical signals of a patient |
US6886010B2 (en) * | 2002-09-30 | 2005-04-26 | The United States Of America As Represented By The Secretary Of The Navy | Method for data and text mining and literature-based discovery |
JP2004121539A (ja) * | 2002-10-02 | 2004-04-22 | Seiko Epson Corp | 体動検出装置 |
US7387611B2 (en) * | 2003-04-10 | 2008-06-17 | Matsushita Electric Industrial Co., Ltd. | Physical movement analyzer and physical movement analyzing method |
JP2005267152A (ja) * | 2004-03-18 | 2005-09-29 | Seiko Instruments Inc | 電子歩数計 |
JP2006026092A (ja) * | 2004-07-15 | 2006-02-02 | Matsushita Electric Ind Co Ltd | 加速度情報送信装置、身体運動解析装置および身体運動解析方法 |
US9820658B2 (en) * | 2006-06-30 | 2017-11-21 | Bao Q. Tran | Systems and methods for providing interoperability among healthcare devices |
US7117030B2 (en) * | 2004-12-02 | 2006-10-03 | The Research Foundation Of State University Of New York | Method and algorithm for spatially identifying sources of cardiac fibrillation |
JP2006187469A (ja) * | 2005-01-06 | 2006-07-20 | Seiko Instruments Inc | 運動強度評価装置 |
KR100601981B1 (ko) * | 2005-01-14 | 2006-07-18 | 삼성전자주식회사 | 활동패턴 감시 방법 및 장치 |
EP1867951B1 (en) * | 2005-03-28 | 2014-03-19 | Asahi Kasei EMD Corporation | Traveling direction measuring apparatus and traveling direction measuring method |
FR2886532B1 (fr) * | 2005-06-07 | 2008-03-28 | Commissariat Energie Atomique | Procede et systeme de detection de chute d'une personne |
JP5028751B2 (ja) * | 2005-06-09 | 2012-09-19 | ソニー株式会社 | 行動認識装置 |
JP5120795B2 (ja) * | 2005-11-15 | 2013-01-16 | 学校法人日本大学 | 人の姿勢動作判別装置およびエネルギー消費量算出装置 |
GB0602127D0 (en) * | 2006-02-02 | 2006-03-15 | Imp Innovations Ltd | Gait analysis |
US8055469B2 (en) * | 2006-03-03 | 2011-11-08 | Garmin Switzerland Gmbh | Method and apparatus for determining the attachment position of a motion sensing apparatus |
JP4785640B2 (ja) * | 2006-06-20 | 2011-10-05 | セイコーインスツル株式会社 | 歩数計 |
JP4885637B2 (ja) * | 2006-07-27 | 2012-02-29 | セイコーインスツル株式会社 | 腕装着型電子歩数計 |
JP4885664B2 (ja) * | 2006-09-21 | 2012-02-29 | セイコーインスツル株式会社 | 歩数計 |
JP4894500B2 (ja) | 2006-12-22 | 2012-03-14 | ソニー株式会社 | 歩行波形処理方法及び歩行波形処理装置 |
US7884727B2 (en) * | 2007-05-24 | 2011-02-08 | Bao Tran | Wireless occupancy and day-light sensing |
US8024024B2 (en) * | 2007-06-27 | 2011-09-20 | Stereotaxis, Inc. | Remote control of medical devices using real time location data |
JP4271711B2 (ja) * | 2007-10-02 | 2009-06-03 | 本田技研工業株式会社 | 運動補助装置 |
US8152734B2 (en) * | 2007-11-28 | 2012-04-10 | Pierson Precision Auscultation | System and method for diagnosis of bovine diseases using auscultation analysis |
US8023928B2 (en) * | 2008-01-16 | 2011-09-20 | Intuitive Research And Technology | System and method for monitoring an analog data signal |
KR20090082711A (ko) * | 2008-01-28 | 2009-07-31 | 삼성전자주식회사 | 보행자 항법 시스템에서의 보폭 추정 방법 및 시스템 |
-
2009
- 2009-02-26 FR FR0951212A patent/FR2942388B1/fr active Active
-
2010
- 2010-02-25 CN CN201080009404.7A patent/CN102333483B/zh active Active
- 2010-02-25 US US13/203,461 patent/US9265448B2/en active Active
- 2010-02-25 KR KR1020117022475A patent/KR101718555B1/ko active IP Right Grant
- 2010-02-25 JP JP2011551487A patent/JP5874130B2/ja active Active
- 2010-02-25 EP EP10705367A patent/EP2400889A1/fr not_active Ceased
- 2010-02-25 WO PCT/EP2010/052368 patent/WO2010097422A1/fr active Application Filing
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2875320Y (zh) * | 2006-02-28 | 2007-03-07 | 深圳市万机创意电子科技有限公司 | 人体步行传感器 |
CN101327125A (zh) * | 2007-06-23 | 2008-12-24 | 株式会社百利达 | 步行评价系统、步行计、步行评价程序和记录介质 |
Non-Patent Citations (1)
Title |
---|
Coordination of head and trunk accelerations during walking;J. J. Kavanagh et al;《European Journal of Applied Physiology》;20050701;第94卷(第4期);第468-475页 * |
Also Published As
Publication number | Publication date |
---|---|
KR20110125656A (ko) | 2011-11-21 |
KR101718555B1 (ko) | 2017-03-21 |
EP2400889A1 (fr) | 2012-01-04 |
FR2942388B1 (fr) | 2012-10-12 |
JP2012518506A (ja) | 2012-08-16 |
US9265448B2 (en) | 2016-02-23 |
US20120041713A1 (en) | 2012-02-16 |
JP5874130B2 (ja) | 2016-03-02 |
CN102333483A (zh) | 2012-01-25 |
WO2010097422A1 (fr) | 2010-09-02 |
FR2942388A1 (fr) | 2010-08-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102333483B (zh) | 用于检测人的步行的系统和方法 | |
Bourke et al. | Evaluation of waist-mounted tri-axial accelerometer based fall-detection algorithms during scripted and continuous unscripted activities | |
Shan et al. | A wearable pre-impact fall detector using feature selection and support vector machine | |
EP3563764B1 (en) | Determining body postures and activities | |
EP2925213B1 (en) | Method and apparatus for identifying transitions between sitting and standing postures | |
Bennett et al. | Inertial measurement unit-based wearable computers for assisted living applications: A signal processing perspective | |
Godfrey et al. | Activity classification using a single chest mounted tri-axial accelerometer | |
Nyan et al. | Classification of gait patterns in the time–frequency domain | |
US10258257B2 (en) | Quantitative falls risk assessment through inertial sensors and pressure sensitive platform | |
CN102753095B (zh) | 用于确定呼吸信号的方法和设备 | |
Zhang et al. | Continuous functional activity monitoring based on wearable tri-axial accelerometer and gyroscope | |
Wagenaar et al. | Continuous monitoring of functional activities using wearable, wireless gyroscope and accelerometer technology | |
Popp et al. | A novel algorithm for detecting active propulsion in wheelchair users following spinal cord injury | |
Lee et al. | A study on the activity classification using a triaxial accelerometer | |
Guenterberg et al. | An automatic segmentation technique in body sensor networks based on signal energy | |
Mezghani et al. | Context aware adaptable approach for fall detection bases on smart textile | |
Yang et al. | Assessing the utility of smart mobile phones in gait pattern analysis | |
Iervolino et al. | A wearable device for sport performance analysis and monitoring | |
Chen et al. | Gait pattern classification with integrated shoes | |
Rodríguez-Martín et al. | Posture transition analysis with barometers: contribution to accelerometer-based algorithms | |
Zhang et al. | A theoretic algorithm for fall and motionless detection | |
Madej et al. | Optimal placement of IMU sensor for the detection of children activity | |
Biswas et al. | On fall detection using smartphone sensors | |
Yang et al. | Algorithm design for real-time physical activity identification with accelerometry measurement | |
López-Nava et al. | Towards ubiquitous acquisition and processing of gait parameters |
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 |