Method and device for fall detection using inertial capture

Abstract

The invention relates to a device for fall detection, comprising at least one accelerometer (2) and means (8) which are programmed to: receive data of at least one acceleration signal originating from the accelerometer (2); identify, from the data from the accelerometer which exceed a certain threshold, data of at least one signal liable to constitute a fall event; identify, from the data of a signal liable to constitute a fall event and in a time window around this event, at least the longest duration for which the signal is uninterruptedly below the mean value of the detected signal; and classify, on the basis of the longest duration, the event as a fall or as a non-fall by means of a classifier based on at least one decision tree and / or based on a gradient descent (XBG or "gradient boosting tree") and / or based on an additional tree (or ETC, or "extremely randomised trees").

Description

[0001] DESCRIPTION

[0002] TITLE: Method and device for fall detection by inertial capture

[0003] TECHNICAL FIELD AND PRIOR ART

[0004] The invention relates to a method and a device for automatic fall detection based on one or more sensors worn by a physical person.

[0005] The accelerometer is the most commonly used sensor for fall detection, including in recent studies. The gyrometer is less common, as its performance contribution does not always justify the resulting increase in power consumption. The barometer is very rarely used due to its high noise and possible redundancy with the accelerometer.

[0006] The first detection models were based on one or more thresholds on characteristics extracted from the signals such as the acceleration norm, the speed at the point of impact, or the angular change. Despite correct announced detection rates, the very high number of false alarms (several per day) did not make them suitable for daily use. In institutions for people requiring care (for example: EHPAD), these devices were abandoned by the care teams because of these false positives which generated additional work for teams already very busy elsewhere. Raising the detection thresholds solves these problems of false alarms (better specificity), but deteriorates the performance in terms of good detection (less good sensitivity).

[0007] To solve the problem of false alarms, approaches based on supervised learning have emerged, as well as larger databases (several dozen subjects). The first models used (Support Vector Machines, k-Nearest Neighbors, Random Forests,) have gradually given way to ensemble methods such as Gradient Boosting, and to deep learning.

[0008] These typically use recurrent networks such as Long Short Term Memory to classify falls from raw accelerometer signals. They offer significant performance gains over traditional machine learning methods, but at the expense of increased complexity. The article by N. Zurbuchen et al. “A comparison of machine learning algorithms for fall detection using wearable sensors”, Conference Paper of the IEEE International Conference on Artificial Intelligence in Information and Communication, 19 / 02 / 20 - 21 / 02 / 20, DOI: 10.1109 / ICAIIC48513.2020.9065205, describes a fall detection system based on the analysis of accelerometric and gyrometric signals, and on the extraction of certain relatively complex and expensive temporal and frequency predictors in terms of computation time (such as the Power Spectral Density) and on the implementation of an ensemble classifier (“Gradient Boosting Trees”).

[0009] The article by TR Mauldin et al. "Smartfall: a smartwatch-based fall detection system using deep learning" published in Sensors, 2018, 18, 3363, NDPI, doi:10.3390 / sl8103363 proposes a system based on deep learning via Recurrent Neural Networks (RNN) which is very costly in terms of system resources, energy and computation time.

[0010] The problem arises of finding a new method and a new device for fall detection, faster, more relevant and preferably less energy-consuming than known devices.

[0011] STATEMENT OF THE INVENTION

[0012] The invention firstly relates to a fall detection device comprising at least one accelerometer and means, or digital means, programmed to:

[0013] - receiving data from said accelerometer relating to at least one acceleration signal or receiving data from at least one acceleration signal from said accelerometer;

[0014] - identify, for example by a filter, among the data from said accelerometer which exceed a certain threshold, data from at least one signal likely to constitute a fall event; - identify, among said data likely to constitute a fall event and in a time window around, or including, this event, at least the greatest duration, or 1 ère duration during which the signal is, uninterruptedly, below the average value of the detected signal;

[0015] - classify, preferably according to said greatest duration, said data or the event, as falling or not falling, by regression of a probability or using a classifier based for example on at least one decision tree and / or on a gradient descent (XGB, or “Gradient Boosting tree”), and / or an additional tree (or

[0016] “Extremely Randomized Tree”, or ETC, or even “Extra Trees”).

[0017] The invention implements a fall predictor (the longest duration during which the signal is, uninterruptedly, below the average value, also called

[0018] "longest_strike_below_mean") which, in combination with a classification, based for example on at least one decision tree and / or on a gradient descent (XGB, or a "Gradient Boosting tree"), and / or on an additional tree (or "Extremely Randomized Tree", or ETC, or "Extra Trees"), proves to be very relevant for the detection of falls, in particular because it is very differentiating between, on the one hand, falls and, on the other hand, activities of daily life, while being very fast and energy-efficient.

[0019] This predictor corresponds to the duration of the longest "strike" (continuous segment) located below the average of a signal, when this is likely to be a signal reflecting a fall. According to one implementation, it is calculated on a 1 ère time window of, for example, 1 second before the main event and on a 2 èmetime window of, for example, 2 seconds, after the main event. The signal used can be that of the accelerometer standard, which may have been filtered by a high-pass filter, for example again with a cut-off frequency fc=0.2Hz.

[0020] This predictor is similar to the duration of the fall. It was identified, among approximately 1000 time and frequency predictors, as being more relevant for fall detection than other commonly used predictors (speed at the point of impact in particular). It is even more effective when the sensor is worn on the wrist, which is the preferred solution for seniors.

[0021] The means, or digital means, may be programmed to calculate, from the data received from the accelerometer, the dynamic component of the acceleration norm, which is then used to identify data from at least one signal likely to constitute a fall event.

[0022] The classifier used may be:

[0023] - at least Random Forest type, with for example between 2 and 20 trees, each with a depth that can be between 2 and 12;

[0024] - or implement at least one Limited Random Forest (RFT), for example limited to 6 trees;

[0025] Alternatively, the event can be classified as a fall or a non-fall by random forest regression.

[0026] Applied to the results of the fall predictor mentioned above,

[0027] (“longest_strike_below_mean”), these classifiers give very good results in terms of accuracy and specificity of fall detection.

[0028] A device according to the invention may further comprise means, or digital means, for calculating, for the accelerometer signal data, or accelerometer data, one or more of the other values ​​among the maximum value of these accelerometer data, the minimum value of these accelerometer data, the median value of these accelerometer data, the standard deviation of these accelerometer data, the variance of these accelerometer data, the root mean square value of these accelerometer data, the sum of the values ​​of these accelerometer data, and for entering at least one of said other calculated values ​​into said classifier.

[0029] A device according to the invention may further comprise at least one other sensor chosen from a gyrometer and / or a barometer. The means, or digital means, may be programmed to calculate, from the data received respectively from a gyrometer and / or a pressure sensor, the angular velocity norm and / or the short-term variation of the local atmospheric pressure, which is / are then used to identify data from at least one signal likely to constitute a fall event.

[0030] A device according to the invention may further comprise means, or digital means, programmed to calculate:

[0031] - for the gyrometer signal data, or gyrometer data, one or more of the largest period during which the signal is continuously below the mean value of the detected signal, the maximum value of such gyrometer data, the median value of such gyrometer data, the minimum value of such gyrometer data, the root mean square value of such gyrometer data, the standard deviation of such gyrometer data, the variance of such gyrometer data, and to input at least one of these calculated values ​​for the gyrometer signal data into said classifier;

[0032] - and / or for the barometer signal data, or barometer data, one or more of the greatest period during which the signal is continuously below the mean value of the detected signal, the maximum value of such barometer data, the mean value of such barometer data, the median value of such barometer data, the minimum value of such barometer data, the root mean square value of such barometer data, the standard deviation of such barometer data, the sum of the values ​​of such barometer data, the variance of such barometer data and to input at least one of these calculated values ​​for the barometer signal data into said classifier.

[0033] In a device according to the invention the time window around, or including, the event can be at most 1.5 seconds before the event and / or at most 2.5 seconds after the event.

[0034] The means, or digital means, for identifying whether an event is likely to constitute a fall may be programmed to identify whether an event is being tracked, for a certain duration, or 2 ème duration, for example between 2 and 4 seconds, of another event, or is not followed by another event: repetitive events can thus be eliminated from events likely to constitute a fall. The means, or digital means, for identifying whether an event is likely to constitute a fall can implement a state machine.

[0035] A device according to the invention may further comprise means, or digital means, for, or programmed to, send an alert based on the results generated by the means, based on the classification result.

[0036] The invention also relates to a method for detecting a fall using at least one accelerometer, this method comprising at least the following steps: a) - collecting at least acceleration data from said accelerometer; b) - identifying, among the data from said accelerometer, data likely to constitute a fall; c) - identifying, using said data likely to constitute a fall event and in a time window around, or including, this event, at least the greatest duration, or 1 èreduration, during which the signal is, uninterruptedly, below the average value of the detected signal; d) - classify, preferably according to said greatest duration, said data or the event, as falling or not falling, by regression of a probability or using a classifier based for example on at least one decision tree and / or on a gradient descent (XGB, or “Gradient Boosting tree”), and / or on an additional tree (or

[0037] “Extremely Randomized Tree”, or ETC, or “Extra Trees”). This method can implement a device according to the invention, as presented above and / or in the remainder of this application.

[0038] The classifier is for example of the Random Forest type, with for example between 2 and 20 trees, each with a depth which can be between 2 and 12.

[0039] Said classifier can implement at least one limited Random Forest (RFT), for example limited to 6 trees (rfjight, or "Random Forest limited"). Before step b), it is possible to calculate, from the data received from the accelerometer, the dynamic component of the acceleration norm, which is then used during steps b) and following.

[0040] Alternatively, the event may be classified as a fall or a non-fall by random forest regression. Such a method may further comprise a step of calculating, for the accelerometer signal data, or accelerometer data, one or more of the other values ​​among the maximum value, the minimum value of these accelerometer data, the median value of these accelerometer data, the standard deviation of these accelerometer data, the variance of these accelerometer data, the root mean square value of these accelerometer data, the sum of the values ​​of these accelerometer data, and inputting at least one of these other values, or calculated values, into said classifier.

[0041] Such a method may further comprise measuring data from said movement using at least one other sensor chosen from a gyrometer and / or a barometer and / or collecting at least data from at least one other sensor chosen from a gyrometer and / or a barometer.

[0042] Before step b), it is possible to calculate, from the data received respectively from a gyrometer and / or a pressure sensor, the norm of the angular velocity and / or the short-term variation of the local atmospheric pressure, which is / are then used during steps b) and following.

[0043] Such a method may further include a calculation step:

[0044] - for the gyrometer signal data, or gyrometer data, one or more of the largest period during which the signal is continuously below the mean value of the detected signal, the maximum value of such gyrometer data, the median value of such gyrometer data, the minimum value of such gyrometer data, the root mean square value of such gyrometer data, the standard deviation of such gyrometer data, the variance of such gyrometer data, and inputting at least one of these values ​​calculated for the gyrometer signal data into said classifier;

[0045] - and / or for the barometer signal data, one or more of the largest period during which the signal is continuously below the mean value of the detected signal, the maximum value of such barometer data, the mean value of such barometer data, the median value of such barometer data, the minimum value of such barometer data, the root mean square value of such barometer data, the standard deviation of such barometer data, the sum of the values ​​of such barometer data, the variance of such barometer data, and inputting at least one of these values ​​calculated for the barometer signal data into said classifier.

[0046] In a method according to the invention, said time window around the event is for example a maximum of 1.5 seconds before the event and / or a maximum of 2.5 seconds after the event.

[0047] Step b) can be implemented by identifying whether an event is tracked, for a certain duration, or 2 ème duration, for example between 2 and 4 seconds, of another event, or is not followed by another event: repetitive events can thus be eliminated from events likely to constitute a fall.

[0048] According to one embodiment, step b) implements a state machine.

[0049] Finally, a method according to the invention may include a step of sending an alert based on the results of step d) of classification: if an event is classified as a fall, it is possible to send an alert to another person, for example to a caregiver.

[0050] A method according to the invention may be implemented using digital means, for example a microcontroller or a processor or a microprocessor, programmed for this purpose.

[0051] More generally, a device or method according to the invention can be implemented in a bracelet or in a pendant.

[0052] BRIEF DESCRIPTION OF THE DRAWINGS

[0053] [FIG 1] schematically illustrates an embodiment of a device according to the present invention.

[0054] [FIG 2] represents an example of a state machine implementation of event predetection. [FIG 3] represents an example of time windowing.

[0055] [FIG 4] represents the time interval corresponding to the predictor of greatest duration below the mean;

[0056] [FIG 5A] and [FIG 5B] represent the classification of the various predictors studied, carried out by a Random Forest type classifier, for a bracelet (figure 5A) and a pendant (figure 5B).

[0057] [FIG 6A] and [FIG 6B] represent performances obtained, with a Random Forest type classifier, for a bracelet (figure 6A) and a pendant (figure 6B);

[0058] [FIG 7] illustrates the contribution of different sensors to overall performance;

[0059] [FIG 8] compares the respective performances of different classifiers;

[0060] [FIG 9] schematically illustrates one embodiment of a method according to the present invention.

[0061] [FIG 10A] and [FIG 10B] represent the influence of the hyper-parameters of a random forest on the performance (measured in average recall).

[0062] DETAILED DESCRIPTION OF SPECIFIC EMBODIMENTS

[0063] A portable fall detection system 10 according to the invention is shown schematically in Figure 1 and comprises at least one sensor for detecting movements of a person, for example at least one accelerometer 2 (preferably 3-axis), and possibly at least one gyrometer 4 (preferably 3-axis) and / or at least one pressure sensor, or barometer, 6. Figure 1 represents a system comprising an accelerometer 2, a gyrometer 4 and a barometer 6, but other configurations can be produced, as explained later.

[0064] Such a system preferably comprises digital means 8, for example a microcontroller or a processor or a microprocessor, making it possible to receive the data emitted or measured by the sensor(s) of the device and programmed to implement an analysis and / or processing of the data as described below, preferably in real time. Optionally, these means 8 can trigger the emission of an alert signal 12. These means 8 can further comprise storage means for storing data or instructions implementing an analysis and / or processing of the data as described below and possibly storage means for storing all or part of the data measured by at least one of said sensor(s) and / or data which result from one or more of the processing operations described below.

[0065] Such a system can be incorporated into an everyday object worn by a physical person, for example a bracelet or a pendant.

[0066] In the case of the example of the sensors of Figure 1 and mentioned above, an accelerometer 2, a gyrometer 4 and a barometer 6, the signals that they produce during the movements, and in particular during a fall, of a person to whom they are linked, are respectively an acceleration (of components (accN_x), (accN_y), (accN_z)), an angular velocity and a measurement of the local atmospheric pressure. These signals are digitized by the means 8 and can be preconditioned, that is to say subjected to one or more signal preprocessing steps aimed at extracting one or more physical quantity(ies) of interest.

[0067] Concerning the accelerometer 2, the means 8 can for example be programmed to calculate the dynamic component of the acceleration norm, the norm being given by the formula:

[0068] The accelerometer standard (e.g. sampled at Fe = 100 Hz) makes it possible to obtain a quantity independent of the person's attitude (and therefore of the sensor's orientation). The dynamic component of the acceleration (as opposed to a component due to gravitation), denoted accNdyn_mg, is for example obtained by passing the acceleration data through a high-pass digital filter, for example of the Butterworth type of order 2 and with a cut-off frequency which is for example of the order of fc = 0.2 Hz. This eliminates the component due to gravitation.

[0069] For the 3-axis gyrometer, which can be sampled at, for example, 100 Hz, the means 8 can be programmed to calculate the norm gyrN_mg of the angular velocity, by a formula identical or similar to that given above. For the barometer, the means 8 are preferably programmed to retain the short-term variation of the local atmospheric pressure (and eliminate slow variations due to meteorology). The signals from the barometer arrive for example at 50 Hz; they can be resampled at, for example, 100 Hz to simplify the time windowing (in order to have a frequency common to the three sensors). The variations due to atmospheric pressure can be compensated by applying a high-pass filter, for example with a cut-off frequency of 0.2 Hz. The corrected pressure is denoted pressuredyn_mbar.

[0070] Preferably, the same filter is used for both the accelerometer and the barometer, which is a resource-efficient solution for processing resources (e.g. RAM and CPU resources).

[0071] Using the means 8, it is possible to achieve a 1 erprocessing, or pre-processing, called "event detection" of at least part of the data produced by the sensor(s). Indeed, certain events, which have nothing to do with a fall, for example a hand movement, can nevertheless cause at least one sensor to produce a signal of intensity greater than a certain threshold value, a signal which must therefore a priori be taken into account. In order to reduce energy consumption, pre-processing is preferably carried out to eliminate signals of this type which cannot result from a fall. This is particularly the case for accelerometer signals generated by repetitive events; an event, identified by a signal which exceeds a certain threshold (set by the operator in the means 8), is repetitive if it is repeated within a certain time window, for example of a duration of the order of 2 to 4 s, for example 3.5 s.On the other hand, an event not followed by another event during this period is likely to constitute a fall.

[0072] For example, using means 8:

[0073] - we identify cases of exceeding a threshold on the dynamic component of the acceleration;

[0074] - if this threshold is exceeded, a more in-depth analysis is triggered, as described in this application and in particular below; - in the case of repetitive movements, the repetitive activity is waited for to end before triggering the more in-depth analysis.

[0075] In other words, an analysis is triggered in the case of repetitive movements, but preferably only once, i.e. at the end of the succession of movements: if there was no pre-processing, for example by a state machine as explained below, an analysis would be triggered at each of the movements, thus increasing the energy consumption of the microcontroller.

[0076] This pre-processing implements for example a state machine using as input the signals of the dynamic component of the acceleration accNdyn_mg (obtained as explained above). It can be programmed in the means 8 and is illustrated in figure 2. The first state (2-1) "WAITING 0 » corresponds to the expectation of an event.

[0077] The dynamic component of the acceleration accNdyn_mg is calculated and compared to a threshold H_THR (which is previously defined by an operator in the means 8; for example it was defined by systematic search according to a regular grid, the criterion being the average recall (“recall_mean”). If this threshold is exceeded by the dynamic component, the machine goes to state (2-2) IMPACT_DETECTE and then initializes a counter (“impact_detected_timer”).

[0078] In state 2-2 (“IMPACT_DETECTED”), if a new threshold exceedance by accNdyn is detected, the counter is reset to zero. Otherwise, it is incremented by +1 for each new sample (for example at 104 Hz). When the counter exceeds a predetermined duration (“IMPACT_DETECTED_TIMOUT”, which can, for example, again have been defined by systematic search according to a regular grid, the criterion being the average recall (“recall_mean”)), the machine goes into state 2-3 “FALL_LIKE_EVENT” (“FALL_LIKE_EVENT”) and triggers the rest of the analysis process (windowing, predictor extraction and classification, see below). It then returns to state 2-1 “WAITING”.

[0079] The transition to the "FALL EVENT" state occurs after a number of samples defined by the IMPACT_DETECTED_TIMOUT duration, after the last threshold exceedance. To know the index of the sample corresponding to the main event (last threshold exceedance), we apply the following formula: Sample index (Event_index) = Current index (current_index) - IMPACT_DETECTED_TIMOUT

[0080] Table 1 below specifies examples of H_THR and IMPACT_DETECTED_TIMOUT threshold values, for a bracelet and for a pendant. [Table 1]

[0081] Time windowing involves selecting a range of samples on either side of the last threshold violation recorded by the state machine.

[0082] This last threshold exceedance is therefore located at the time defined by the duration “IMPACT_DETECTED_TIMOUT” samples before the transition to the “FALL_LIKE_EVENT” state (which is the literal translation of the equation above).

[0083] Another step, called time windowing, can then be implemented by the means 8 to select before and / or after the event (which may have been identified in the manner explained above), the data of interest to be analyzed more precisely in order to confirm, or not, that this event is a fall.To this end, as illustrated in Figure 3 (which represents the dynamic component of the acceleration, expressed in milligrams (mg)), a time window F is identified with a predetermined total duration (by an operator in the means 8), for example a few seconds, for example between 3 and 4 seconds, around the last transition to the IMPACT-DETECTED state: a first part of the window is located before the last transition to the IMPACT-DETECTED state, a second part of the window is located after it; according to an exemplary embodiment with a window with a total duration of 3 seconds, this is divided into seconds (corresponding for example to 100 samples) before the last transition to the IMPACT-DETECTED state and into 2 seconds (corresponding for example to 200 samples) after the last transition to the IMPACT-DETECTED state.Indeed, in the case of a real fall, only the data in such a limited time window around the last transition to the IMPACT-DETECTED state can contain relevant information. It is therefore this data that will make it possible to decide whether the detected signals really reflect a fall.

[0084] This windowing step can be applied to the signals of at least one of the sensors mentioned above, preferably to the accelerometer signals; more particularly, it can be applied to one or more of the 3 quantities previously calculated during the conditioning phase: to the dynamic component of the accelerometer norm accNdyn_mg, to the gyrometer norm gyroN_mdps, to the dynamic component of the pressure pressuredyn_mbar.

[0085] As explained further on, a device 1 according to the invention may also implement only an accelerometer 2, or an accelerometer 2 and a gyrometer 4, or an accelerometer 2 and a barometer 6, in which case the above explanations are transposed to the signals of the sensor(s) implemented.

[0086] The data identified within the time window described above can then be subjected to statistical analysis as explained below.

[0087] To this end, the inventors identified a statistical indicator, or predictor, of falling. A systematic search for relevant predictors of falling was carried out using the "tsfresh" tool described in the article by M. Christ et al. "Time Series FeatuRe Extraction on the basis of Scalable Hypothesis tests," published in Neurocomputing, vol. 307, pp. 72-77, Sept. 2018, doi: 10.1016 / j.neucom.2018.03.067.

[0088] This tool allows the extraction of approximately 1000 predictors most commonly used in time series analysis, and the selection of those that are relevant for separating the two classes of events (falls and non-falls). This selection is based on a calculation of p-value, or probability, which can be seen as the smallest significance threshold for which the null hypothesis is accepted, as explained in the article by J.D. Gibbons et al., “P-Values: Interpretation and Methodology,” Am. Stat., vol. 29, no. 1, p. 20-25, 1975, doi: 10.2307 / 2683674, with control of the false discovery rate by the Benjamini Yekutieli method (Yoav Benjamini, Daniel Yekutieli. "The control of the false discovery rate in multiple testing under dependency." Ann. Statist. 29 (4) 1165 - 1188, August 2001. https: / / doi.org / 10.1214 / aos / 1013699998). This method consists of controlling the false discovery rate by adjusting the p-value threshold according to the number of predictors tested.Among the tested predictors, 8 predictors per sensor can be selected (i.e. a maximum of 24 in total) according to a compromise, for the application to fall detection, between performance and simplicity of implementation.

[0089] These predictors are listed in Table 2 below:

[0090] [Table 2] >

[0091] In this table 2:

[0092] - the sorting order is simply alphabetical. The units appear in the predictor name (for the gyrometer, the unit is millidegrees per second (mdps));

[0093] - the “typical value” corresponds to that obtained for a given subject. It is therefore an example.

[0094] These predictors are therefore the "longest strike below mean", the maximum value, the mean value (but not for gyroscope data), the median value, the minimum value, the root mean square value, the sum of the values ​​(but not for gyroscope data), the variance. The "longest_strike_below_mean" can for example be realized (in means 8) using the algorithm below:

[0095] 1. For each value of x, test the condition: x < mean(x): we then obtain a vector of booleans;

[0096] 2. Group this vector of Booleans into groups of the same consecutive values. For example, the sequence (0,0,0, 1, 1,0,0) will give three groups: (0,0,0) (1,1) and (0,0)

[0097] 3. Measure the number of elements in each group (respectively 3, 2 and 2)

[0098] 4. Return the largest number of elements in its groups. In the example given above, this number is 3.

[0099] Figure 4 (the dynamic component of the acceleration, expressed in milliG (mg)) represents a time window, around a transition to the IMPACT-DETECTED state, which may possibly be that of a fall, which represents the greatest duration below the average. In practice, this time window can be partially before the event (the event = the peak in Figure 4) and partially after the event, for example 1 s before the event and 2 s after the event. The relevance of these predictors to separate falls from non-falls was evaluated using the method of Benjamini and Yekutieli (see the article by these authors already cited above). Here we are mainly interested in the importance attributed to them by the classifier used (random forest as explained below). The degree of relevance is represented in Figures 5A and 5B, respectively for a bracelet and a pendant.

[0100] In the case of the bracelet (Figure 5A) (containing the 3 sensors), the “longest strike below mean” predictor calculated on the acceleration comes in second position, after the average of the atmospheric pressure. In the case of the pendant (Figure 5B) (also containing the 3 sensors), it is considered as the most important predictor by the classifier.

[0101] The data can then be classified or merged, for example by a "random forest" type classifier, selected after comparing the performance of different classifiers. A classification algorithm can be implemented or programmed in the digital means 8. The input data for this classifier are one or more of the predictors mentioned above and / or in Table 2 above.

[0102] In this example, this random forest has 5 decision trees, each making a decision independently of the others. The final decision is obtained by voting.

[0103] The principle of Random Forest is as follows:

[0104] 1. Creation of B new training sets by a double sampling process:

[0105] On the observations (an observation is a time window containing the data of an event (fall or non-fall), using a draw with replacement of a number N of observations identical to that of the original data, a. And on the p predictors, by retaining only a sample of cardinal m < sqrt (p)

[0106] 2. On each sample, a decision tree is trained using one of the known techniques, limiting its growth by cross-validation. 3. The B predictions of the variable of interest are stored for each original observation.

[0107] The prediction of the random forest can then be a simple majority vote (Ensemble learning).

[0108] The parameters of this classifier can be optimized by constant step search grid.

[0109] For the implementation of this process, it has been shown that a number of trees less than 20 or 10 or even 6 was sufficient, with, for each tree, a depth of between 2 and 12, for example 10, or even between 5 and 7.

[0110] It is possible to implement only one decision tree, which is less efficient than a random forest, but may be of interest in the case of an implementation of the algorithm on a microcontroller with limited resources (CPU and / or RAM resources), even if this may be to the detriment of performance. We can estimate the influence of the 2 main parameters of the random forest:

[0111] - the number of trees (n_estimators);

[0112] - the depth (or maximum number of nodes) of each tree (“max_depth”).

[0113] To do this, we vary these parameters according to a regular grid of step 2. The results obtained are illustrated in figures 10A (case of a bracelet) and 10B (case of a pendant), on which the light colors correspond to the best performances, the performance being measured in terms of “average recall” (or “mean_recall”).

[0114] In these figures 10A and 10B we see a similar profile for the pendant and for the bracelet, that is to say an increase in performance with the increase in the number of trees and / or their depth.

[0115] The performance of this fall detection system was estimated on a separate training and validation database according to Table 3 below (ADL meaning “Activity of Daily Life”); this database consisted of fall data from 30 people (“subjects”).

[0116] [Table 3]

[0117] The performances obtained are gathered in figures 6A and 6B, given respectively for a bracelet and a pendant.

[0118] These figures show that:

[0119] - For the pendant (Figure 6A): 374 cases of activities of daily living and 230 cases of falls were correctly classified; only 7 cases of activities of daily living and 10 cases of falls were misclassified;

[0120] - For the bracelet (Figure 6B): 360 cases of activities of daily living and 211 cases of falls were correctly classified; only 21 cases of activities of daily living and 29 cases of falls were misclassified.

[0121] Tables 4 (for a bracelet) and 5 (for a pendant) below gather precision (column 1) and specificity (column 2) data; the 3rd column (fl score) is the harmonic mean between precision and recall. The "support" column corresponds to the number of observations on which the calculation of the indicators is based. The macro average corresponds to the average of the two upper lines. The weighted average corresponds to the average of the two lines weighted by the support.

[0122] [Table 4]

[0123] [Table 5]

[0124] Figure 7 illustrates the contribution of the different sensors to the overall performance. The

[0125] "mean_recall" (or average recall), which is on the abscissa and in Table 6 below, corresponds to the average between sensitivity and specificity; the designation of the sensors (accelerometer = acc, gyrometer = gyr, barometer = pressure) is followed by the letter N when the standard is extracted and by "dyn" when the dynamic component is used. In addition, for each of the combinations considered of the accelerometer, the gyrometer and the barometer, the upper horizontal bar in Figure 7 corresponds to the case of the pendant, the lower horizontal bar corresponds to the case of the bracelet. [Table 6]

[0126] From this figure and these data, it is understood that, within the framework of the present invention: - the accelerometer can be used alone;

[0127] - either the accelerometer can be used in combination with the gyroscope or with the barometer; - or the accelerometer can, as already explained above, be used with both the gyroscope and the barometer.

[0128] The best performance is obtained when all 3 quantities (acceleration, angular velocity and pressure) are used ("average recall" = 92% for the bracelet and 97% for the pendant). When the accelerometer is used alone, the performance is lower ("average recall" = 84% for the bracelet and 93% for the pendant), although still satisfactory.

[0129] Other variants are possible, by changing the classifier used. As illustrated in Figure 8, the respective performances of different classifiers were compared, these different classifiers being applied to fall data such as those presented above. Here again, the "average recall" is on the abscissa. The different classifiers are used with their default parameters. The abbreviations shown in Figure 8 correspond to: Etc: ExtraTrees Classifier (Extremely Randomized Trees) Trees, RF: Random Forest, xgb: Gradient Boosting Trees, rfjight: Random Forest limited to 5 trees, dt: Decision Trees, svm: Support Vector Machines knn: k-Nearest Neighbors, gnb: Gaussian Naive Bayes, mlp: MultiLayer Perception (Neural Network). For each of them, the upper horizontal bar in Figure 8 corresponds to the pendant case, the lower horizontal bar corresponds to the bracelet case.In this figure, we see that the “neural network” type classifier is the least efficient.

[0130] We see from this figure that, for the application to fall detection, the 4 best classifiers (RF, ETC, XGB, RF_LIGHT), notably the random forest, are of the ensemble type, that is to say that they aggregate the decisions of several classifiers. This makes it possible to make the decision more robust for the application targeted in the present application and to optimize the bias / variance trade-off.

[0131] It is possible to implement only one decision tree, which is less efficient than a random forest, but may be of interest in the case of porting the algorithm to a microcontroller with limited resources.

[0132] It is best to implement at least the 4 classifiers (RF, ETC, XGB, RF_LIGHT) mentioned above. Alternatively, one can use a regression of a probability, for example a Random Forest regression instead of a Random Forest classification. Unlike a classifier that provides a binary value (1 for fall or 0 for non-fall, or vice versa), regression provides a continuous value between 0 and 1 similar to a probability of fall. The decision to classify as a fall or not is then made by setting a threshold on this probability. If the probability is greater than the threshold, then an alert is triggered, otherwise, the event is ignored.

[0133] The advantage is that this allows the user to set the sensitivity threshold for fall detection himself, depending on the intended application. The threshold can be programmed in the means 8. It can be adjusted by an operator.

[0134] For example, in a context of high risk of falling (highly dependent people), the trigger threshold can be set very low to detect all falls. In the case of low risk of falling (independent and healthy people), the threshold can be set higher to reduce the number of false alarms.

[0135] The invention finds applications in particular:

[0136] To detect falls in elderly people (at home or in a specialized establishment (EHPAD or CHU);

[0137] To detect falls of workers working on difficult terrain (construction sites, roofs, scaffolding, snow, ice)...

[0138] Figure 9 represents an embodiment of a method according to the invention, comprising 5 steps executed chronologically:

[0139] A step S1 of conditioning the signals collected from the sensor(s); A step S2 of event detection;

[0140] A time windowing step S3;

[0141] A step S4 of calculating at least one fall predictor;

[0142] An S5 classification step.

[0143] Each of these steps has been described above. In this embodiment, comprising steps S1 and S2, steps S3 - S5 are preferably activated only in the event of a positive response from the event detection in order to save system energy. In certain cases, it is possible not to implement one or the other of steps S1 and S2:

[0144] - step S1 (which calculates the signal norm) allows working with orientation-independent data; in certain applications, it is preferable to implement this step, particularly for senior citizens, who do not always wear an object, for example a pendant or a bracelet, which implements the invention in the same direction or on the same hand. Calculating the norm overcomes this problem;

[0145] - step S2 contributes to savings in terms of power supply (battery in particular), because without it the calculation of predictors and classification must operate continuously.

[0146] Optionally, the result of step S5 may trigger the emission of an alert signal 12 (figure 1).

[0147] The "longest_strike_below_mean" predictor can be implemented using the following code:

[0148] 1. For each value of x, test the condition: x < mean(x). We obtain a vector of booleans

[0149] 2. Group this vector of Booleans into groups of similar consecutive values. For example, the sequence 0,0,0,1,1,0,0 will give three groups: (0,0,0) (1,1) and (0,0)

[0150] 3. Measure the number of elements in each group (respectively 3, 2 and 2)

[0151] 4. Return the largest number of elements in its groups; in the given example, this number is 3.

Claims

CLAIMS 1. Fall detection device comprising at least one accelerometer (2) and digital means (8) programmed to: - receiving data from at least one acceleration signal from said accelerometer (2); - identify, among the data from said accelerometer which exceed a certain threshold, data from at least one signal likely to constitute a fall event; - identify, among said data of a signal likely to constitute a fall event and in a time window including this event, at least the greatest duration, or 1 ère duration during which the signal is, uninterruptedly, below the average value of the detected signal; - classify, based on said greatest duration, the event as a fall or non-fall by regression of a probability or using a classifier based at least on a decision tree and / or on a gradient descent (XGB, or “Gradient Boosting tree”), and / or on an additional tree (or ETC, or “Extremely Randomized Trees”) 2. Device according to claim 1, the event being classified as a fall or a non-fall by random forest regression or said classifier being at least of the Random Forest type, with for example between 2 and 20 trees, each with a depth of between 2 and 12.

3. Device according to claim 1 or 2, said classifier implementing at least one limited Random Forest (RFT), for example limited to 6 trees 4. Device according to one of claims 1 to 3, further comprising means (8) programmed to calculate, for the data of the accelerometer signal, or accelerometer data, one or more of the other values ​​among the maximum value of these accelerometer data, the minimum value of these accelerometer data, the median value of these accelerometer data, the standard deviation of these accelerometer data, the variance of these accelerometer data, the average value quadratic of such accelerometer data, the sum of the values ​​of such accelerometer data, and inputting at least one of said other values ​​into said classifier.

5. Device according to one of claims 1 to 4, further comprising at least one other sensor (4, 6) chosen from a gyrometer and / or a barometer.

6. Device according to claim 5, further comprising means (8) programmed to calculate: - for the gyrometer signal data, or gyrometer data, one or more of the largest period during which the signal is continuously below the mean value of the detected signal, the maximum value of such gyrometer data, the median value of such gyrometer data, the minimum value of such gyrometer data, the root mean square value of such gyrometer data, the standard deviation of such gyrometer data, the variance of such gyrometer data, and to input at least one of these calculated values ​​for the gyrometer data into said classifier; - and / or for the barometer signal data, or barometer data, one or more of the greatest period during which the signal is continuously below the mean value of the detected signal, the maximum value of such barometer data, the mean value of such barometer data, the median value of such barometer data, the minimum value of such data, the root mean square value of such data, the standard deviation of such barometer data, the sum of the values ​​of such barometer data, the variance of such data and to input at least one of these calculated values ​​for the barometer data into said classifier.

7. Device according to one of claims 1 to 6, said time window comprising the event being at most 1.5 seconds before the event and 2.5 seconds after the event.

8. Device according to one of claims 1 to 7, the means (8) for identifying whether an event is likely to constitute a fall being programmed to identify whether an event is followed, for a certain duration, or 2 ème duration, of another event, or is not followed by another event.

9. Device according to claim 8, said 2 ème duration being between 2 and 4 seconds.

10. Device according to one of claims 1 to 9, the means (8) for identifying whether an event is likely to constitute a fall implementing a state machine.

11. Device according to one of claims 1 to 10, further comprising means for sending an alert on the basis of the results generated by the means (8), on the basis of the classification result.

12. Method for detecting a fall using at least one accelerometer (2), this method comprising at least the following steps: a) - collecting at least acceleration data from said accelerometer (2); b) - identifying, among the data from said accelerometer, data likely to constitute a fall; c) - identifying, using said data likely to constitute a fall event and in a time window including this event, at least the greatest duration, or 1 ère duration, during which the signal is, uninterruptedly, below the average value of the detected signal; d) - classify, according to said greatest duration, the event as falling or not falling, by regression of a probability or using a classifier based on at least one decision tree and / or on a gradient descent (XGB, or “Gradient Boosting tree”) and / or on an additional tree (or ETC, or “Extremely Randomized Trees”).

13. Method according to claim 12, the event being classified as a fall or a non-fall by random forest regression or said classifier being of the Random Forest type, with for example between 2 and 20 trees, each with a depth of between 2 and 12.

14. Method according to claim 12 or 13, said classifier implementing a limited Random Forest (RTF), for example limited to 6 trees 15. Method according to one of claims 12 to 14, further comprising a step of calculating, for the data of the accelerometer signal, or accelerometer data, one or more of the other values ​​among the maximum value of these data of accelerometer, the minimum value of such accelerometer data, the median value of such accelerometer data, the standard deviation of such accelerometer data, the variance of such accelerometer data, the root mean square value of such accelerometer data, the sum of the values ​​of such accelerometer data, and input at least one of said other values ​​into said classifier.

16. Method according to one of claims 12 to 15, further comprising measuring data from said movement using at least one other sensor chosen from a gyrometer and / or a barometer.

17. Method according to claim 16, further comprising a calculation step: - for the gyrometer signal data, or gyrometer data, one or more of the largest period during which the signal is continuously below the mean value of the detected signal, the maximum value of such gyrometer data, the median value of such gyrometer data, the minimum value of such gyrometer data, the root mean square value of such gyrometer data, the standard deviation of such gyrometer data, the variance of such gyrometer data, and inputting at least one of these values ​​calculated for the gyrometer data into said classifier; - and / or for the barometer signal data, or barometer data, one or more of the greatest period during which the signal is continuously below the mean value of the detected signal, the maximum value of such barometer data, the mean value of such barometer data, the median value of such barometer data, the minimum value of such barometer data, the root mean square value of such barometer data, the standard deviation of such barometer data, the sum of the values ​​of such barometer data, the variance of such barometer data, and inputting at least one of such values ​​calculated for the barometer data into said classifier.

18. Method according to one of claims 12 to 17, wherein said time window around the event is at most 1.5 seconds before the event and 2.5 seconds after the event.

19. Method according to one of claims 12 to 18, step b) being implemented by identifying whether an event is tracked, for a certain duration, or 2 ème duration, of another event, or is not followed by another event.

20. Method according to claim 19, said 2 ème duration being between 2 and 4 seconds.

21. Method according to one of claims 12 to 20, step b) implementing a state machine.

22. Method according to one of claims 12 to 21, further comprising a step of sending an alert based on the results of step d) of classification.

Images