The present invention relates to a method for determining toe off. More precisely, the present invention relates to a method for determining toe off with an insole comprising a plurality of sensing cells.
The present invention also relates to a computer program product to carry out the steps of the method.
In various fields, such as medical fields or sport training, it is often desirable to determine, during the gait cycle, the time at which a foot of a user does not contact the ground anymore allowing the determination of the end of the stance phase and the beginning of the swing phase (respectively P1 and P2 in
To do so, sensing cells are often used. In order to obtain the toe off as precisely as possible, the signal measured by each sensing cell must reflect exactly the pressure applied to said sensing cell. However, dielectric materials usually used in sensing cells are subject to latency, also called elastic recovery. This latency is a phenomenon occurring during the decay of the pressure applied to a sensing cell and manifesting itself as a signal decay from the cell which is delayed compared to the variation of the pressure. Therefore, the sensing cell still measures a signal when the pressure applied to this cell is null. This phenomenon thus leads to a time delay in the order of several milliseconds when measuring the variation of pressure from sensing cells. In the case of the gait, latency occurs during the pre-swing phase (P15 in
Known methods for the determination of the toe off using sensing cells are not satisfactory.
Indeed, several known methods use thousands of sensing cells to capture the geometry and relative arrangement of each footfall as a function of time. Moreover, several known methods require combining data captured from various sensors such as pressure sensors, camera, gyroscope, angular velocity sensors . . . making the method complex, time consuming and tedious. Moreover, these systems are expensive and can only be used within clinical and controlled laboratory environments. In addition, since the analysis is performed in a confined and controlled environment, patients may not exhibit their natural movement patterns because a limited number of strides/gestures can be observed.
A purpose of the invention is to provide a toe off detection method that is more precise, more reliable or that may be performed at lower cost outside controlled laboratory environments.
To this end, the present invention relates to a method for determining toe off with an insole comprising a controller and a plurality of sensing cells, the method comprising:
Indeed, the method being performed only during the discharge period, this allows to reduce the energy consumption of the insole compared to known methods performed during the whole gait cycle.
Moreover, outputting, for each sensing cell, a null output if a signal is lower than a predetermined cell calibration threshold allows to suppress most of the noisy signal produced, among other things, by the latency of the sensing cells allowing a method which is more precise.
Furthermore, triggering output of the toe off when the total signal is lower than a predetermined toe off threshold allows to correct from any possible remaining noisy signal.
Finally, the method is based only on sensing cells comprised in an insole. This leads to a method which is cheaper than known methods allowing performing the method outside controlled laboratory environments.
According to another advantageous aspect of the invention, the predetermined toe off threshold corresponds to a force ranging between 5 and 40N.
According to another advantageous aspect of the invention, identifying a discharge period comprises:
According to another advantageous aspect of the invention, triggering the determination of the beginning of the discharge period is performed when the total signal corresponds to a force lower than 300N.
This is advantageous because this avoids any false triggering due to the noise of each sensing cell which may imply, out of the actual discharge period, a total signal computed at the preceding instant which is larger than the total current signal.
According to another advantageous aspect of the invention, the method further comprises, before identifying a discharge period, determining heel strike comprising:
This determination allows to increase the precision of the method. Indeed, by including this determination, the discharge period is identified only during the stance phase (i.e., when the foot is on the ground). This allows to avoid an identification of the discharge period during the swing phase due to a noise measurement from the sensing cells.
Moreover, since this allows to perform the computations for the identification of the discharge period only during the stance, this leads to a method which is less energy consuming compared to known methods performed during the whole gait cycle.
Finally, the determination of the heel strike also allows to determine the duration of the contact between the foot and the ground (i.e., the duration between the heel strike and the toe off) with a higher precision.
According to another advantageous aspect of the invention, the predetermined heel strike threshold corresponds to a force ranging between 5 and 40N.
According to another advantageous aspect of the invention, the total signal is computed as the sum of the signals outputted by the controller.
In an advantageous aspect of the invention, sensing cells measure a pressure and the total signal is the integration of a field of pressure, said field of pressure being a spatial interpolation of pressure signals outputted by the controller.
According to another advantageous aspect of the invention, at least one sensing cell is a capacitive pressure cell or a force sensing cell.
According to another advantageous aspect of the invention, the calibration phase is performed by measuring a dynamic response of sensing cell with a load decreasing from its maximum value to zero in 0.12 to 0.7 s, corresponding to usual gait behaviour.
This allows to increase the precision of the method. Indeed, measuring the dynamic response of the sensing cells allows to take the rate of load (pressure or force) variations into account since it modifies the response of the sensing cell.
According to another advantageous aspect of the invention, the sensing cells are distributed along the surface of the insole.
This allows to compute a distribution of the load along the surface of the insole with a higher precision. Moreover, the higher the number of sensing cells, the lower is the value error on the computed total signal.
According to another advantageous aspect of the invention, at least one of the sensing cells is located in the toe part of the insole.
This characteristic is advantageous in normal gait, where the toes are the last contact with the ground during the pre-swing phase. Indeed, since the toe off is the moment when the foot does not contact the ground anymore, in normal gait, locating at least one sensing cell in the toe part of the insole allows to compare the predetermined toe off threshold to a total signal to which at least the toe part of the insole contributes.
According to another advantageous aspect of the invention, at least one of the sensing cells is located in the heel part of the insole.
This characteristic is advantageous in abnormal gait, where the toes are not the last contact with the ground during the pre-swing phase. Indeed, in abnormal gait, locating at least one sensing cell in the heel part of the insole allows to compare the predetermined toe off threshold to a total signal to which at least the part of the insole other than the toe part contributes.
This invention also relates to a computer program product comprising instructions which, when the program is executed by a controller, cause the controller to carry out the steps of the method.
In the present invention, the following terms have the following meanings:
A method 100 according to the present invention allows to determine the toe off. The toe off is a moment separating two consecutive phases of the gait cycle (or stride) represented in
Advantageously, the method 100 determines the toe off with an insole 1 comprising a plurality of sensing cells 10 (an example of insole which may be used with the method 100 is shown in
For example, the insole 1 may be the insole disclosed in WO2017033036. Advantageously, the sensing cells 10 are distributed along the surface of the insole 1. This is advantageous because the distribution of the load along the surface of the insole 1 can be computed with a higher precision leading to a higher precision of the method 100. Moreover, the higher the number of sensing cells, the lower the error on the values computed from the overall sensing cells 10. Preferably, the insole 1 comprises a number of sensing cells 10 comprised in a range between 5 and 30, more preferably between 15 and 20.
In one embodiment, at least one of the sensing cells 10 is located in the toe part 2 of the insole 1 (
In one embodiment, at least one of the sensing cells 10 is located in the heel part 3 of the insole 1 (
The sensing cells 10 may be, for example, a pressure sensing cell or a force sensing cell. If the sensing cell is a pressure sensing cell, the sensing cell may be, for example, a capacitive pressure cell, a resistive pressure cell or a piezoelectric pressure cell. Capacitive pressure cells may comprise a thin sheet of deformable dielectric material placed between two electrodes. Indeed, the value of the capacitance C of a capacitive pressure cell can be determined as a function of the thickness L of the dielectric sheet of the capacitive pressure cells, the surface S of the upper electrode and the lower electrode of the capacitive pressure cells and the dielectric constant ε of the material between the electrodes by the following equation: C=εS/L. Consequently, when the thickness L of the dielectric sheet of the capacitive pressure cells is changed, the capacitance C varies. Change of thickness may be induced by the pressure applied by the user on the insole 1 while walking or exercising. If the sensing cell is a force sensing cell, the sensing cell may be, for example, a resistive force sensing cell—also known as force-sensing resistors FSR. A resistive force sensing cell is a sensor whose resistance changes when a force, pressure or mechanical stress is applied. It may be composed of superimposed layers arranged as: two substrate layers surrounding a conductive film and a plastic spacer, which includes an opening aligned with the conductive film. Above the plastic spacer there may be a conductive substrate. When an external force is applied to the sensing cell, the conductive film is deformed towards the substrate. Air from the opening of the plastic spacer is pushed through an air vent, and the conductive film comes into contact with the conductive substrate. The more of the conductive substrate gets in touch with the conductive film, the lower the resistance. Therefore, the more pressure applied on the sensor, the more the conductive substrate touches the conductive film making the resistance go down.
The insole 1 used for the method 100 also comprises a controller 15 to collect the measurements performed by the sensing cells 10 and/or to perform the computations associated to the method 100. The controller 15 may include a processor. The expression “processor” should not be construed to be restricted to hardware capable of executing software, and refers in a general way to a processing device, which can for example include a computer, a microprocessor, an integrated circuit, or a programmable logic device (PLD). The processor may also encompass one or more Graphics Processing Units (GPU), whether exploited for computer graphics and image processing or other functions. Additionally, the instructions and/or data enabling to perform associated and/or resulting functionalities may be stored on any processor-readable medium such as, e.g., an integrated circuit, a hard disk, a CD (Compact Disc), an optical disc such as a DVD (Digital Versatile Disc), a RAM (Random-Access Memory) or a ROM (Read-Only Memory). Instructions may be notably stored in hardware, software, firmware or in any combination thereof.
Preferably, the insole 1 is worn and used (by example, during a gait) by the user during the overall application of the method 100. In an alternative embodiment, the measurements from the sensing cells 10 are recorded by the controller 15 while the insole 1 is worn and used by the user and the method 100 is performed after the measurement session.
The signal measured by each sensing cell 10 may be pre-processed. The pre-processing may include a filtering, typically a smoothing filtering, so that high frequency variations in signals during stance are removed.
Preferably, the signal measured by each sensing cell 10 is a series of timestamped data, i.e., each signal is a succession of discrete data characterized by temporal information (measurement time) identifying the time when the data is measured. In another embodiment, the measurement session is divided in several measurement periods, each measurement period defining one measurement time. The overall data measured during a measurement period are then timestamped by the measurement time associated to said measurement period.
A first step 120 of the method 100 (
In one embodiment, the step 120 may comprise several steps represented in
The total signal may be then computed from the signals collected from the overall sensing cells 10 for said measurement time t′ (step 124). The total signal may be a simple sum of the signals output for the overall sensing cells 10. Alternatively, the total signal may be the result of the integration of a spatial interpolation of the signals collected from the overall sensing cells 10 as the method explained in the European patent application EP 20 306 521.4.
The total signal of the measurement time t′ may be then compared to a total signal computed at a preceding instant t′-Nτ (step 126), where 1/τ represents the frequency of signal acquisition. Preferably, the total signal of the measurement time t′ is compared to a total signal computed at the time just before, i.e., at the measurement time t′-τ.
From this comparison, the determination of the beginning of the discharge period may be triggered when the total signal computed at the preceding instant t′-Nτ is larger than the total signal of the current measurement time t′ (step 128).
To avoid any false triggering due to the noise of each sensing cell which may imply a total signal computed at the preceding instant t′-Nτ which is larger than the total signal of the measurement time t′ while being out of the actual discharge period (this effect can occur mainly during P13 and P14 in
Another way to avoid any false triggering is to trigger the determination of the beginning of the discharge period when the total signal corresponds to a force lower than 300N. Indeed, the minimal weight of an adult is about 40 kilograms corresponding to 400N. Setting the force threshold to 300N allows to be sure that the weight of the user has begun to move on the other feet and thus that the decay measured between two successive instants actually corresponds to the discharge period.
Yet another way to avoid any false triggering is to determine heel strike (step 110) before the step 120 of identification of the discharge period. In one embodiment, step 110 comprises several steps represented in
In one embodiment, the predetermined heel strike threshold is a signal value corresponding to a force ranging between 5 and 40N, preferably between 10 and 30N, even more preferably between 15 and 25N. The conversion between a value of the force and a signal measured by a sensing cell 10 may be performed thanks to the transfer function, determined during the calibration phase, between at least one signal measured by the sensing cell 10 and corresponding force applied. Therefore, the force values of each sensing cell 10 may be computed in order to determine the total force which will be compared to the predetermined heel strike threshold.
Once the discharge period is identified, the controller 15 of the insole 1 may collect, during a step 130, the signal measured by each sensing cell 10. The collected signals may then be analyzed by the controller 15 during a step 140 by comparing the signal of each sensing cell 10 to a predetermined cell calibration threshold.
The predetermined cell calibration threshold is determined during the calibration phase for each sensing cell 10. In other words, each sensing cell 10 is associated with its own predetermined cell calibration threshold which may be different from those of the other sensing cells. During the calibration phase, a load may be applied to a sensing cell 10 and then decreased (unloading phase) until the load becomes null. Preferably, as shown in
In one embodiment, to determine the calibration curve, at least one data from the signal measured by said sensing cell 10 is associated to the value of the corresponding load—a pressure or a force. The predetermined cell calibration threshold is then determined as the signal measured by said sensing cell 10 when the load becomes null. In some embodiments, the predetermined cell calibration threshold may be determined by interpolation or as the signal corresponding to the first signal measured when the load becomes null.
In a specific embodiment, the calibration phase is performed by measuring a dynamic response of a sensing cell 10 with a load decreasing from its maximum value to zero in 0.12 to 0.7 s, corresponding to actual gait behavior, from low pace (long time for load decrease) to rapid pace (short time for load decrease). Said sensing cell 10 is then characterized by several calibration curves depending on the variation rate of the applied load. To determine each calibration curve, at least one data from the signal measured by said sensing cell 10 for each variation rate of the load is associated to the value of the corresponding load. The method of dynamic calibration may be, for example, the method disclosed in the European patent application EP 21 306 275.5. The predetermined cell calibration threshold may be then determined for each calibration curve as the signal measured by said sensing cell 10 when the load becomes null. In some embodiments, the predetermined cell calibration threshold associated to a variation rate of the load is determined, for said variation rate of the load, by interpolation of the corresponding calibration curve or as the signal corresponding to the first signal measured when the load becomes null.
In this specific embodiment in combination with pressure sensing cells, the applied load is a controlled pressure varying with a rate ranging from 5 to 30 kg·cm−2·s−1, preferably from 10 to 20 kg·cm−2.
In this specific embodiment in combination with force sensing cells, the applied load is a controlled force varying with a rate ranging from 600 to 10000 N·s−1, preferably from 1400 to 4200 N·s−1.
The predetermined cell calibration threshold may be a signal value corresponding to a pressure ranging between 0.05 and 0.5 kg/cm2, preferably between 0.15 and 0.3 kg/cm2. The conversion between a value of the pressure and a signal measured by a sensing cell 10 may be performed thanks to the transfer function, determined during the calibration phase, between the signal measured by the sensing cell and the pressure applied.
At a step 150, for each sensing cell 10 for which the signal is larger than or equal to the predetermined cell calibration threshold, the controller outputs the signal measured by said sensing cell(s) 10. Advantageously, for all other sensing cell(s) 10 (i.e., those for which the measured signal is lower than the predetermined cell calibration threshold), the controller outputs a null signal (step 160). This is advantageous because this allows to suppress most of the noisy signal produced, among other things, by the latency of the sensing cells allowing a method which is more precise.
This is illustrated in
After the controller 15 has output the measured signals from the overall sensing cells, the total signal of the sensing cells is computed (step 180). The total signal may be a simple sum of the signals output for the overall sensing cells and for the same measurement time. Alternatively, the total signal may be the result of the integration of a spatial interpolation of the signals output for the overall sensing cells 10 at the same measurement time, for example using the method explained in the European patent application EP 20 306 521.4.
Finally, the toe off is triggered when the total signal is lower than a predetermined toe off threshold (step 190).
In one embodiment, the predetermined toe off threshold is a signal value corresponding to a force ranging between 5 and 40N, preferably between 10 and 30N, even more preferably between 15 and 25N. The conversion between signals measured by all sensing cells 10 and the value of the total force may be performed in various ways. In an embodiment, thanks to the transfer function between the signal measured by the sensing cell and the force applied, determined during the calibration phase, the force values of each sensing cell 10 are computed in order to determine the total force which will be compared to the predetermined toe off threshold. In another embodiment, the signals measured by all sensing cells 10 are pressure signals and may be spatially interpolated in a continuous pressure field, then integrated over the surface of the insole to yield the total force, as taught in European patent application EP 20 306 521.4.
In the embodiment where the method 100 is performed after the measurement session, the method 100 may be performed at once after the loading of the timestamped signal measured by the overall sensing cells 10. For example, the method 100 optionally starts with the step 110 for the determination of the heel strike by computing the total signal for each measurement time. The heel strike may be triggered as the first measurement time of a stride when the total signal exceeds the predetermined heel strike threshold. The discharge period is then identified (step 120) by comparing the total signal for each measurement time to the total signal at a preceding measurement time. The beginning of the discharge period may be triggered as the first measurement time of a stride (optionally occurring after the heel strike) when the total signal is lower than those at the preceding measurement time and optionally corresponding to a force lower than 300N. Considering only the signals measured during the discharge period, the step 140 is performed by comparing the signal of each sensing cell 10 to a predetermined cell calibration threshold. The step 150 or 160 is then performed according to the result of this comparison. The step 180 is then performed, by computing the total signal for each measurement time during the discharge period. The toe off is then triggered (step 190) as the first measurement time when the total signal is lower than a predetermined toe off threshold.
In the embodiment where the insole 1 is worn and used by the user during the overall application of the method 100, the method 100 may be performed at each measurement time after the signal measurement by the overall sensing cells 10. For example, considering the measurement time t, the method 100 optionally starts with the step 110 for the determination of the heel strike by computing the total signal measured at t. The heel strike t1 is then triggered when the total signal exceeds the predetermined heel strike threshold. The beginning of the discharge period t2 is then triggered when the total signal measured at a time t (optionally greater than t1) is lower than those at the preceding measurement time (t-Nτ) and optionally corresponding to a force lower than 300N. At t>t2, the step 140 is performed by comparing the signal of each sensing cell 10 to the corresponding predetermined cell calibration threshold. The step 150 or 160 is then performed according the result of this comparison. Still on said time t>t2, the step 180 is performed by computing the total signal and the toe off is triggered (step 190) when the total signal is lower than a predetermined toe off threshold.
The invention may also relate to a computer program product comprising instructions which, when executed by a computer, cause the computer to carry out the method 100 as disclosed above.
The expression “computer” should not be construed to be restricted to hardware capable of executing software, and refers in a general way to a processing device, which can for example include a computer, a microprocessor, an integrated circuit, or a programmable logic device (PLD). The computer may also encompass one or more Graphics Processing Units (GPU), whether exploited for computer graphics and image processing or other functions. Additionally, the instructions and/or data enabling to perform associated and/or resulting functionalities may be stored on any processor-readable medium such as, e.g., an integrated circuit, a hard disk, a CD (Compact Disc), an optical disc such as a DVD (Digital Versatile Disc), a RAM (Random-Access Memory) or a ROM (Read-Only Memory). Instructions may be notably stored in hardware, software, firmware or in any combination thereof.
The present invention is further illustrated by the following comparison performed between the stance duration measured using the method of the present invention (INV) and using a state of the art (SOA) method. Both methods were applied using an insole comprising 18 sensing cells, as shown in
The method of the present invention comprises the steps of determining heel strike (associated to the beginning of the stance) with a predetermined heel strike threshold of 20N, identifying a discharge period, and, during the discharge period, triggering output of the toe off (associated to the end of the stance) with a predetermined toe off threshold of 20N.
In the SOA method, the beginning of the stance is defined when the total load applied to the sensing cells of the insole becomes greater than a threshold of 20N and the end of the stance when the total load applied to the sensing cells of the insole becomes null. The duration of the stance is thus the difference between the end and the beginning of the stance.
In both methods (INV and SOA), the user is wearing shoes comprising the insole and walks on a BP400600 Biomechanics Force Platform supplied by AMTI. The force sensed by the platform is considered as the reference to determine the actual beginning and end of the stance.
Five users-without gait impairment-were asked to walk at low, normal and high speed. For each user and each walking pace, 60 stances were recorded, yielding a total of 900 stances representative of different gaits, and analyzed for each method. The experiment was reproduced with four insoles comprising different dielectric materials in their capacitive pressure cells. Table 1 below reports results obtained. For each insole, statistics over the 900 stances are reported for the method of the present invention and of the state of the art. The mean relative absolute error on the duration of the stance is defined as
where fpi is the measure obtained by the force platform and ini is the measure obtained with the insole. The bias of stance duration is the average difference between force platform and insole stance duration and is defined as
Last, the upper limit of the 95% confidence interval—i.e., bias plus 1.96 standard deviation—is reported for the difference between force platform and insole stance duration measures.
For the four types of dielectric materials, results show that the method according to the invention allows a much more precise determination of stance duration. The mean absolute relative error is significantly decreased, as well as the bias and the upper limit of the 95% confidence interval. In addition, measurements with insoles do not require the specialized laboratory equipment of a force platform, but can be performed in every day conditions, when a user is wearing insoles.
| Number | Date | Country | Kind |
|---|---|---|---|
| 21306522.0 | Oct 2021 | EP | regional |
| Filing Document | Filing Date | Country | Kind |
|---|---|---|---|
| PCT/EP2022/080243 | 10/28/2022 | WO |