We examine the usefulness of event-based sampling approaches for reducing communication in inertial-sensor-based analysis of human motion. To this end we consider realtime measurement of the knee joint angle during walking, employing a recently developed sensor fusion algorithm. We simulate the effects of different event-based sampling methods on a large set of experimental data with ground truth obtained from an external motion capture system. This results in a reduced wireless communication load at the cost of a slightly increased error in the calculated angles. The proposed methods are compared in terms of best balance of these two aspects. We show that the transmitted data can be reduced by 66% while maintaining the same level of accuracy.

Human motion analysis has a multitude of applications in sports, rehabilitation, orthopedics and similar fields. The gold standard for clinical usage, i.e. optical motion analysis, is expensive and restricted to laboratory environments. Wireless inertial sensor networks have recently become popular as an alternative, enabling ambulatory human motion analysis at much lower costs. However, in such networks large amounts of data must be transmitted wirelessly, which represents a limiting factor[1] in realtime applications, such as biofeedback, neuroprosthetics and feedback control of active orthoses. One possible solution is to use a more sophisticated sampling and transmission strategy instead of standard fixed-rate sampling.

In contrast to standard sampling, in which samples are transmitted at a fixed rate, the core idea of event-based sampling is to only transmit samples when certain events, which are derived from the measured signals, occur.[2] Different event-based sampling methods have been proposed; see [1], [2] for an overview. Send-on-delta [3] is a basic sampling scheme, where a threshold on the difference between the current sample and the last sent sample is used as the transmission criterion. When integrating this difference over time the send-on-area scheme is obtained, which was proposed in [4] and [5]. In the context of state estimation, innovation-based triggering [6], [7] and variance-based triggering [8] have been proposed, where an update is sent whenever the estimator innovation (i.e. the difference between a measurement and its prediction) or its variance exceeds a threshold. These algorithms have also been used for feedback control, such as stabilization of a balancing sculpture [6].

Two related approaches are called adaptive sampling and model-based active sampling [9]. In adaptive sampling, the sampling rate is adjusted based on the signal history. See [10] for a realization of adaptive sampling as an analog circuit and [11] for an example of adaptive sampling applied to general sensor networks. For model-based active sampling, a model is used to predict sensor readings. Those predictions are either used to decide if a measured sample should be transmitted or even to avoid activating the sensor at all. An example using autoregressive (AR) models can be found in [12].

In the present contribution, we examine the usefulness of event-based sampling approaches for inertial-sensor-based analysis of human motion. To this end we consider realtime measurement of the knee joint angle during walking. We employ a recently developed sensor fusion algorithm [13] that avoids the use of magnetometer readings. Applying this algorithm to a large set of experimental data yields an average root-mean-square error (RMSE) of about 3° if standard full-rate sampling is used. We propose four different event-based sampling algorithms, apply them to the recorded data sets in a simulated-online manner and feed the resulting signals to the joint angle calculation algorithm. This leads to data compression (i.e. a lower communication load) at the cost of an increased error in the calculated angles. The proposed methods are compared in terms of best balance of these two aspects.

First, consider the following standard sampling case: A sensor measures a vector-valued signal at a fixed sampling rate and transmits it with a fixed sampling rate to a receiver, which then processes the data. In contrast to this standard case, we introduce an event-based sampling protocol that defines under which conditions the sensor sends the current measurement sample to the receiver. At every sampling instant *t*, the sensor evaluates a transmission criterion based on current and old measurements. If the criterion is fulfilled, the current sample is sent. If a sample is sent, it is used directly for data processing. If no sample is sent, the receiver reuses the last received sample instead[3], such that the reconstructed data exhibits the same constant sampling rate as the original data. The following two approaches are considered:

**Send-on-delta (SOD):** In analogy to [3], the sample is sent if the Euclidean norm of the difference between the current measurement and the last sent measurement exceeds a certain threshold.

**Send-on-area (SOA):** In analogy to [4], the Euclidean norm of the aforementioned difference is added up. If this sum exceeds a certain threshold, the current sample is sent and the sum is reset to zero.[4]

**Decreasing threshold:** In their current form, both approaches can result in very long communication gaps. This might be undesirable from an application point of view. Therefore, we propose the following extension that can be combined with both send-on-delta and send-on-area: For each series of consecutively skipped samples, the threshold *c* is linearly decreased in
*n* consecutive samples are skipped. For example, with *n* = 3 the actual thresholds are *c*,
*c*,
*c* and 0.

For the present application, consider two or more Inertial Measurement Units (IMUs), which wirelessly transmit their acceleration **a**(*t*) ∈ *ℝ*^{3} and angular rate **g**(*t*) ∈ *ℝ*^{3} readings to a receiver, which then uses a realtime algorithm to estimate motion parameters from the reconstructed data. Each IMU applies either send-on-delta or send-on-area with either constant or decreasing threshold to the current measurements **a**(*t*) and **g**(*t*) separately. If one or both criteria are fulfilled, the IMU sends a combined data sample consisting of both **a**(*t*) and **g**(*t*) to the receiver.[5]

In our test setup, two IMUs per leg are employed to calculate the knee angle for each leg using the algorithm described in [13], which works roughly as follows: The angle *α*_{gyr}(*t*) is calculated from the angular rate **g**(*t*) and the angle *α*_{acc}(*t*) from the acceleration **a**(*t*). The angle *α*_{gyr}(*t*) is accurate but drifting, while *α*_{acc}(*t*) does not exhibit drift, but is less accurate, especially during fast movement. Sensor fusion of both angles is used to obtain the final angle *α*_{IMU}(*t*). For details, please refer to [13].

A large collection of walking data, consisting of 396 trials from 11 healthy subjects, along with reference angles *α*_{opt}(*t*) obtained by an optical 3D motion capture system, is available. If we use standard full-rate sampling, the RMSE[6] between IMU angles and reference angles is 3.18°. The proposed algorithms are implemented as follows: The recorded original IMU data is used as input for a simulated sender, which evaluates the described criteria. A simulated receiver feeds the reconstructed data to the knee angle calculation algorithm (cf. Figure 1).

We consider four different sampling methods as shown in Table 1. Send-on-area is only applied to the angular rate **g**(*t*). The thresholds *a*_{th} and *g*_{th} are varied from zero to a manually chosen maximum in a 16 × 16 grid. For comparison, we also test sending only every second and third sample, effectively downsampling the signal.

Table 1:

Method | a(t) |
g(t) |
Thresholds a_{th} and g_{th} |
---|---|---|---|

Method 1 | SOD | SOD | constant |

Method 2 | SOD | SOA | constant |

Method 3 | SOD | SOD | decreasing, n = 5 |

Method 4 | SOD | SOA | decreasing, n = 5 |

For evaluation, we analyze how the RMSE with respect to the optical reference angles changes when event-based sampling is applied. For each method, we calculate the RMSE and the compression ratio when varying the thresholds *a*_{th} and *g*_{th} independently, as depicted in Figure 2. The compression ratio is calculated as

To compare the methods, we create a plot of the (smallest achievable) RMSE over the compression ratio. Figure 3 compares all four tested methods. This plot is obtained in the following way: The 16 × 16 grids of the RMSE and the compression ratio are interpolated to a much higher resolution. We then split the range from zero to the maximum compression ratio in 20 intervals of equal size. For each interval, we find the pair of RMSE and compression ratio which has the smallest RMSE and a compression ratio within the given interval. Connecting those pairs results in the curves shown in Figure 3. For illustration, the connection line of the described pairs (for Method 2) is shown in the contour plot in Figure 2.

Send-on-delta event-based sampling surprisingly performs far worse than simple downsampling. Using send-on-area for the angular rate significantly reduces the error, which is less surprising, since the calculation of *α*_{gyr} involves integration of the angular rates. Compared to downsampling, however, the algorithm yields only slightly better results. The same is true for using send-on-delta in combination with decreasing thresholds. However, when combining send-on-area sampling with decreasing thresholds (Method 4), the error can be reduced further. For example, it is possible to reduce the transmitted data by 66% while increasing the total average RMSE by only 0.3°, as shown in Figure 4.

We presented four event-based sampling algorithms for human motion analysis via inertial sensor networks in realtime applications like biofeedback, neuroprosthetics and the control of active orthoses. To test the proposed methods, we applied them to a large collection of walking data in a simulated-online manner and analyzed how the error increases when less data samples are transmitted. The results show a potential for greatly reducing wireless communication load with a low impact on accuracy. For example, it was found that 64% of the transmitted samples can be dropped while increasing the RMSE of a knee joint angle by only 0.3°.

The presented methods were tested for one very specific scenario. Considering other setups and algorithms for human motion analysis is planned for future work. Employing model-based predictions [6], [8] instead of zero-order hold has the potential for reducing wireless communication even further [14]. Investigating the usefulness of such approaches for the considered application is subject to future research.

We would like to express our deep gratitude to Noelia Chia Bejarano and Simona Ferrante from Politecnico di Milano for providing the datasets that were analyzed in this study. Furthermore, we sincerely thank David Graurock for his skillful support in programming for data evaluation.

Research funding: Being conducted in the research project BeMobil, this work is funded by the German Federal Ministry of Research and Education (FKZ 16SV7069K). The work by S. Trimpe was funded by the Max Planck Society. Conflict of interest: Authors state no conflict of interest. Material and Methods: Informed consent: Informed consent has been obtained from all individuals included in this study. Ethical approval: The research related to human use complies with all the relevant national regulations, institutional policies and was performed in accordance with the tenets of the Helsinki Declaration, and has been approved by the authors’ institutional review board or equivalent committee.

[1] Liu Q, Wang Z, He X, Zhou DH. A survey of event-based strategies on control and estimation. Systems Sci Control Engineer. 2014;2:90–7. Search in Google Scholar

[2] Miskowicz M. Event-based control and signal processing. CRC Press; 2016. Search in Google Scholar

[3] Miskowicz M. Send-on-delta concept: an event-based data reporting strategy. Sensors. 2006;6:49–63. Search in Google Scholar

[4] Miskowicz M. Asymptotic effectiveness of the event-based sampling according to the integral criterion. Sensors. 2007;7:16–37. Search in Google Scholar

[5] Nguyen VH, Suh YS. Networked estimation with an area-triggered transmission method. Sensors. 2008;8:897–909. Search in Google Scholar

[6] Trimpe S, D’Andrea R. An experimental demonstration of a distributed and event-based state estimation algorithm. In: Proceeding of IFAC World Congress, 2011. p. 8811–8. Search in Google Scholar

[7] Wu J, Jia QS, Johansson KH, Shi L. Event-based sensor data scheduling: trade-off between communication rate and estimation quality. Automatic Control, IEEE Transactions on, 2013;58:1041–6. Search in Google Scholar

[8] Trimpe S, D’Andrea R. Event-based state estimation with variance-based triggering. IEEE Trans Autom Control 2014;59:3266–81. Search in Google Scholar

[9] Raghunathan V, Ganeriwal S, Srivastava M. Emerging techniques for long lived wireless sensor networks. IEEE Commun Mag. 2006;44:108–14. Search in Google Scholar

[10] Rieger R, Taylor JT. An adaptive sampling system for sensor nodes in body area networks. IEEE Trans Neural Syst Rehabil Eng. 2009;17:183–9. Search in Google Scholar

[11] Jain A, Chang EY. Adaptive sampling for sensor networks. In Proceeedings of the 1st international workshop on Data management for sensor networks: in conjunction with VLDB 2004, ACM; 2004, p. 10–6. Search in Google Scholar

[12] Tulone D, Madden S. PAQ: Time series forecasting for approximate query answering in sensor networks. In Wireless Sensor Networks, Springer; 2006. p. 21–37. Search in Google Scholar

[13] Seel T, Schauer T, Raisch J. IMU-based joint angle measurement for gait analysis. Sensors. 2014;14:6891–909. Search in Google Scholar

[14] Trimpe S, Campi M. On the choice of the event trigger in event-based estimation. In: International Conference on Event-based Control, Communication, and Signal Processing, p. 1–8, Krakow, Poland, June 2015. Search in Google Scholar