Embedded inertial measurement unit reveals pole lean angle for cross-country skiing

This study introduces an inertial measurement unit-based measurement system for resolving the dynamic lean angle of a ski pole during double poling while cross-country skiing. The measurement system estimates both the pole lean angle and pole–terrain contact events. Reported are results from 20 trials providing validated estimates of ski pole lean angle and the timing of pole plant and pole lift events. The pole lean angle is estimated from a complementary filter that fuses estimates of orientation from the embedded accelerometer and angular rate gyro. Validation follows from comparison with video capture measurements. Bland–Altman analysis showed agreement between the two measurement modalities with less than 5% bias in the mean differences (relative to the lean angle range of motion). Companion correlation analysis confirms strong correlation (r=0.99\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r = 0.99$$\end{document}) between the inertial measurement unit and video-estimated lean angles and with mean root-mean-square errors below 4.5∘\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{\circ }$$\end{document}.


Introduction
Poling serves as a major source of propulsion in the classic (diagonal-stride) style of cross-country skiing and becomes the main source when double poling (poling synchronously with both arms) is employed [1]. The pole forces developed during double poling (DP) and the optimum timing to maximize those forces remain of keen interest to coaches and athletes, as they seek to increase poling power and hence skiing speed. This interest spans the entire range of skiing events from sprint courses to long-distance marathons where DP remains prevalent [2]. Prior studies highlight the overall importance and influence of the poling forces arising from DP [1,[3][4][5][6][7][8][9][10][11][12].
Simultaneous knowledge of (1) the total force applied by the pole to the terrain, (2) the angle of the pole with respect to the terrain, and (3) the timing of pole-terrain contact enables one to deduce the component of the pole force applied parallel to the terrain, the time interval over which it is applied, and consequently the forward impulse responsible for propulsion. The ability to measure these three factors while actively skiing outdoors on a course remains a clear challenge. These studies [1,[3][4][5][6][7][8][9][10][11][12] have measured these factors only in laboratory settings or in highly controlled conditions that preclude characterizing propulsion on outdoor courses.
To succeed in studying propulsion on outdoor courses, one could combine pole-embedded force transducers with inertial measurement units (IMUs) to estimate pole orientation with respect to the terrain. Doing so would provide the essential information to estimate the component of the poling force in the plane of the terrain that is responsible for propulsion (rather than weight bearing). However, simply measuring the pole orientation, such as the pole lean angle with respect to the vertical, remains a challenge in outdoor skiing because of the lack of portability of the measurement devices for this task, their high sensitivity to weather conditions, and the high number of accessories needed for their correct operation. Camera-based systems for measuring This  pole orientation based on video [8,[13][14][15] or infrared [1, 5, 10-12, 16, 17] motion capture are ill suited for the purpose of outdoor testing and over the large distances required. A potential solution proposed in this study is a pole-embedded IMU to estimate the dynamic pole lean angle. Moreover, IMUs represent an attractive solution due to their low cost and miniaturization.
While IMUs have been deployed to study skiing technique [18][19][20][21], they have not been used to study the motion of the pole beyond detecting pole-terrain contact events, cycle duration, cycle speed, and cycle length [19,22]. Nevertheless, a pole-embedded IMU also provides the requisite data (acceleration and angular velocity) needed to estimate the lean angle of the pole as well. Thus, the objective of this study was to introduce an IMU-based system to estimate the dynamic pole lean angle of skiers performing the DP technique and, importantly, to determine the accuracy of the estimated pole lean angles by comparison with those measured using video recordings.

Methods
A skilled skier, with experience in cross-country skiing competitions in the range of 30 km in national racing, acted as the test participant following requirements of The Finnish Advisory Board for Research Integrity. Data collecting and handling complied with the requirements of the Finnish Data Protection Ombudsman. The participant provided informed consent for the study.

Experimental facility and protocol
The test facility was a flat and level outdoor ski track equipped for video recording within a short (9 m) field of view; refer to schematic of ski track (top view) in Fig. 1. At the beginning of a trial, the participant warmed up according to individual routine. The participant then performed ten trials utilizing the DP technique with self-selected skiing speeds using classic-style skis. Following a 3-minute break, the participant performed an additional 10 trials (total 20 trials). The field of view of the video camera was sufficient to resolve two poling cycles per trial yielding a total of 40 poling cycles for the study. An embedded IMU provided simultaneous data (linear acceleration and angular velocity) sampled from the pole as described next.

Selection of the skiing technique
The DP technique induces two-dimensional motion of the pole (motion largely in the saggital plane) as seen in previous studies [1,27]. This motion can be measured using a camera system with its optical axes aligned perpendicular to this plane.

Equipment specifications and data analysis
The study employed a ski pole (1.450 m length, R2 Competition Pole AH03-4R2S145, Rex, Finland) that was modified with two markers to support video imaging and custom instrumentation to support inertial motion sensing; refer to Fig. 2a.
For video imaging, one (top) marker was attached near the handgrip (250 mm below the top of the pole) and a second (bottom) marker was attached near the tip of the pole (300 mm above the tip) (Fig. 2a). A stationary camera (Cannon EOS 700, Japan), having a 50 Hz frame rate, filmed the motion of the pole (Fig. 1). The (x, y) coordinates of the pole-attached markers were extracted automatically from the video recordings using the open-source program Tracker 4.11.0 (Open Source Physics, USA) with a position measurement resolution of 1 mm. Given the known distance between the two markers (950 mm), the expected resolution of the lean angle measured in the vertical plane from the video images is arctan(2/950)=0.12 • . The procedure utilized here was similar to that described in [19].
The custom instrumentation for inertial motion sensing consisted of a six-degree of freedom IMU (MPU-6050, InvenSense, USA) capable of resolving both acceleration (acceleration range of ±16 g and a noise floor of 400 μ g / associated data acquisition (DAC) system (prototype, LUT, Finland), both mounted to a miniature pc board measuring 13 mm × 100 mm (Fig. 2b). This slender design enabled the IMU and DAC system to fit within the hollow of the handgrip and with an overall added mass of 53 g (DAC board: 8 g/Battery: 45 g). The DAC sampled the acceleration and the angular velocity at a rate of 300 Hz and wrote the data to a 8 Gb micro SD card. This design enabled continuous sampling over 6 h. Figure 2b illustrates the custom instrumentation and its installation within the grip of the pole (Fig. 2c). A small cap was added to the handgrip that covers the instrumentation (not shown). The IMU was calibrated to reduce bias error and crossaxis sensitivity errors according to the procedure detailed in [23]. The raw IMU data was filtered using a 2 Hz low-pass second-order Butterworth digital filter to yield smoothed data similar to the procedure in [20]. Note that the dominant (pendulum-like) motion of the pole does not possess frequencies significantly above 2 Hz as confirmed by the power spectrum of the acceleration data. All data processing was performed using standard MATLAB functions (The MathWorks, Natick, MA, USA).
To compare the estimates of the pole motion measured by video and IMU data, the two data streams were time synchronized. The time synchronization was achieved using the following three-stage procedure. Each test session, consisting of ten trials, was filmed in one continuous video recording. At the start of each test session, three pole "movement events" were introduced, which were readily discernible in both the video recordings and the IMU data to establish synchronization times. First, the trial started with a static pose of the pole. The static pose period is readily identifiable in the IMU data (period of sustained zero angular velocity) and the start and end of this static pose period are also clearly identifiable in the video recordings. Second, the pole was rotated 90 • in the vertical plane with static poses between two successive rotations. These two static poses yield additional synchronization times. Finally, the pole was rotated (at a priori unknown angular velocity) starting and ending with a static pose. Again, the resulting static poses and angular motions are readily identifiable in both the video recordings and IMU readings. Simultaneous use of the resulting three sets of synchronization times enabled minimization of the synchronization error across these times. It is assumed that this time synchronization was maintained during the length of each of the test sessions attributed to any possible inaccuracy generated from the real-time clock utilized in the custom instrumentation.

Estimation of ski pole orientation in DP
The orientation of the ski pole in DP is largely determined by the lean angle (in the sagittal plane relative to the vertical and the direction of skiing) as presented in [20].
The lean angle is calculated by fusing independent estimates of the vertical direction using the accelerometer data and the angular velocity data by using a modification of the complementary filter presented in [24]. The modified filter leads to an estimate of the "true" angular velocity over each time step by correcting for the differing estimates of vertical (direction of gravity) found from the acceleration data versus the original angular velocity data. Doing so leverages the measured linear acceleration to estimate the bias error and noise present in the angular velocity data and to then remove those error sources in arriving at the "true angular velocity" as follows. Figure 3 illustrates a block diagram of the complementary filter that returns estimates of the lean angle at each time step following five major steps summarized as follows. In the first step, the direction of gravity at the next time step is estimated from the acceleration data assuming that the accelerometer functions as an inclinometer. The resulting direction of gravity is denoted a . Next, the orientation of the IMU at the next time step is estimated by integrating the angular velocity forward one time step. Doing so yields a second estimate of the direction of gravity (i.e., the vertical direction) and this estimate is denoted . In the third step, the (angular) difference of the two estimates yields the misalignment denoted e . The misalignment is then used to correct the original angular velocity over the same time step in arriving at the "true" angular velocity for that time step denoted . Finally, the corrected orientation of the IMU is found by integrating the true angular velocity forward one  Fig. 3, is repeated for each time step. Refer to [24] for further details.

Ski pole event detection
The magnitude of the acceleration measured by the IMU embedded in the grip reveals the major phases of ski pole motion. The motion phases include the plant phase that begins the instant the pole tip is first planted on the terrain, the instant that the tip is lifted from the snow, and the follow-on swing phase that precedes the next plant. Figure 4 illustrates these pole contact events over two complete poling cycles for a Fig. 3 Inputs, steps, and output of the complementary filter implemented in [24] ω IMU a IMU  Fig. 4 Magnitude of the acceleration measured by the IMU embedded in the handgrip reveals the pole plants and lifts over two poling cycles for diagonal stride skiing typical trial. The instants of pole plants are readily observable in this plot of acceleration magnitude versus time, since the highest acceleration peaks arise when the pole first contacts the terrain [25,26]. Pole lifts arise within the lower acceleration peaks between two successive pole plants. During the plant phase, the pole is under compression and that force (directed along the axis of the pole) generates a vertical ground reaction and, importantly, also a horizontal ground reaction that propels the skier forward [1,4,27]. Note that while the acceleration magnitude reveals a time region (± 0.13 s as observed in the plots) where pole lift occurs, the precise time of that lift is undetectable without further analysis.
Our proposed analysis exploits the kinematic constraint that, during the plant phase, the tip of the pole remains planted. Consequently, the acceleration measured at the grip must also depend on the angular velocity (and angular acceleration) of the pole as measured by the grip-embedded angular rate gyro. To reveal that, we introduce the three position vectors illustrated in Fig. 5 that locate the position of the tip of the pole (point T), the position of the IMU within the grip (point I), and the position of T relative to I.
During the plant phase, the angular velocity of the pole, denoted , can be used to estimate the acceleration measured at the location of embedded IMU, denoted I,est , through where is the angular acceleration of the pole (obtained from by central finite differencing). During the plant phase, the magnitude of I,est must remain within a small threshold of that measured directly by the embedded accelerometer, denoted I , per where the superimposed component of gravity in I is duly subtracted. To estimate the quantity threshold a threshold , four poling cycles representing 10% of total cycles in the experiment were randomly selected. A video analysis of these four poling cycles showed the instants of pole lift, and from this a threshold value of a threshold = 0.05 g succeeds in identifying the lifts.

Assessment of agreement
The agreement between the IMU-estimated lean angles and those measured by the video recording was assessed using Bland-Altman analysis [28]. Prior to applying the Bland-Altman analysis, a test for normality was carried out on the distribution of the differences using the Kolmogorov-Smirnov (KS) test [29]. In cases where the differences failed the normality test, a logarithmic transformation was applied. The range of motion (ROM) of the pole lean angle was used to define limits of 5% of the maximum bias for the mean difference and 10% of the maximum limits for the 95% confidence interval. These limits determine if the bias and confidence intervals from the Bland-Altman analysis yield acceptable agreement. The measurement error for the trials was qualitatively categorized (random or systematic) according to the variability of the differences observed in the Bland-Altman plots as detailed in [30].

Results and discussion
The pole lean angle estimated using the IMU data is compared to that measured from the video images. Subsequently, the IMU-derived pole-terrain contact events are compared with those deduced from the video images.

Ski pole lean angle
Across all 20 trials considered in this study, the pole lean angles estimated from IMU data agree with those measured from the video images. This agreement is first observed qualitatively in Fig. 6 which compares the angles estimated by both measurement modalities for poling cycles from two sample trials. These trials, which are representative of all trials, also show variation across poling cycles. In particular, the poling cycle time and range of motion for the trial on the left (Fig. 6a) differ from those shown on the right (Fig. 6b) by 0.21 s and 12 • , respectively. While the experiments were not intended to study cycle-to-cycle variations and possible differences in poling techniques, these variations are easy to observe and may well be useful in studying poling technique in future studies.
(2) �‖ I,est ‖ − ‖ I ‖� ≤ a threshold ,  Tables 1 and 2, respectively, across all 20 trials. Additionally, a sample set of Bland-Altman and correlation plots for a representative trial (that appearing in Fig. 6a) are presented in Fig. 7.
Inspection of the sample Bland-Altman plot, Fig. 7a, reveals that the two measurement modalities yield highly comparable estimates of the pole lean angle. The differences in the angular measurements remain well within the  The table reports the number of samples (n) in each trial, the difference mean ( d ), the standard deviation (SD) of the differences, the 95% confidence interval of the agreement limits (− 1.96d , + 1.96d ), the p values from the KS test to check for normality of the differences, the qualitative type of error observed on the scatter, the bias error, and range of motion (ROM) of the pole lean angle 95% confidence interval (with few exceptions) and possess a mean difference of only 0.99 • . In keeping, the correlation plot for this same example trial, Fig. 7b, reveals tight clustering ( r = 0.994 ) about the best fit line with nearly all data within the 95% confidence interval. The findings from this example trial are consistent with the results reported in Tables 1 and 2 across all trials. For the Bland-Altman results of Table 1, note that the mean differences and the bias errors for each trial lie within the range of (0.17 • < bias error < 5.71 • ), and remain less than 5% of the corresponding ROM of the pole lean angles. It is also noted that the normality test failed in two trials; refer to p values for the KS test for trials 5 and 7. After a visual inspection, no apparent qualitative difference was found between these two trials and those that passed the normality test and a logarithmic transformation was applied to these two trials, as recommended in [28]. For the correlation results of Table 2, note that the two methods remain correlated (0.982 < r < 0.997) with small root-mean-square differences (1.53 • < rms < 4.51 • ).
The above results clearly demonstrate agreement between the lean angles estimated by both measurement modalities. The remaining differences are likely unimportant for the stated goal of using the embedded IMU as the means to understand poling dynamics during skiing on outdoor courses. The pole may bend during the pole plant phase, and thus the angle near the top of the pole (where the IMU is mounted) might differ from the angle at the bottom of the pole. Previous attempts to quantify ski pole bending are based on direct distance measurement systems (such as wires) [31], or indirect measurements using laser or infrared video motion capture systems [10,31]. These methods are confined to laboratory conditions and not suitable for outdoor (e.g., competition-like) conditions which is the focus of this study.

Pole-terrain contact events
The IMU embedded within the grip yields data sufficient to determine the plant and swing phases of the poling cycle. These phases are separated by the instant when the tip first contacts the terrain (end of swing phase and beginning of plant phase) and the instant when the tip subsequently loses contact with the terrain (end of plant phase and start of swing phase). These critical pole-terrain contact events are detected as described in the Methods, and the results from a sample trial are reported next. Figure 8 illustrates three important experimental quantities plotted over a representative poling cycle. The black curve denotes the vertical position (y coordinate) of the lower marker on the pole as measured from the video images. The dotted portion of this curve denotes the plant phase when the tip of the pole remains planted on the terrain as confirmed by video images. The blue curve denotes the magnitude of the pole acceleration sampled at the grip (similar to that illustrated in Fig. 4). The instant when the tip first contacts the terrain occurs at the peak acceleration magnitude (marked by black asterisks). Shortly following this instant, the pole rotates about the planted tip of the pole and the acceleration at the grip satisfies the constraint condition given by Eq. (2). That constraint condition is plotted in green and it remains satisfied, hence the constraint magnitude remains close to zero, until the tip lifts from the terrain when the constraint is no longer satisfied, and the constraint magnitude sharply increases. As illustrated, this sharp increase also arises when the vertical position achieves its minimum and then increases thereafter as the tip lifts from the terrain. Thus, the duration of the pole plant phase, which is responsible for active skier propulsion, lies between these two contact events (i.e., between the pairs of asterisks) and these important events are detectable from the grip-embedded IMU data.
The times of pole plant and lift were estimated in this manner from the IMU data across all 40 poling cycles included in this study. The same times were estimated from the video images. The differences in these times from the two measurement methods were computed across all poling cycles and descriptive statistics are summarized in Table 3, which reports the mean, maximum and minimum differences in these times for both the pole plant and pole lift events. From these distributions, the mean and standard deviation of the differences in pole plant times are − 4 ± 15 ms, while those for the pole lift times are − 7 ± 11 ms. Thus, the two measurement methods yield similar times of pole plant and pole lift (within − 7 ms on average). The remaining differences likely arise from the challenge of identifying the pole-terrain contact events in the video images, some (small) sensitivity to the threshold used in the constraint Eq. (2), and any possible residual sliding of the pole tip while contacting the terrain. (2). The constraint condition is evaluated when the threshold is exceeded, and the value of the constraint is multiplied by a factor of 50 to aid visualization (color figure online)

Further remarks and future work
A future research effort could combine this method with a pole-embedded load cell to simultaneously estimate pole force. Combining the force and angle estimates provides the added capability of estimating the horizontal pole ground reaction that propels the skier forward. The resulting technology may further reveal poling performance across athletic populations.

Conclusions
This paper contributes a new and non-invasive technique to measure the pole lean angle while freely skiing on outdoor courses, thereby surpassing the constraints imposed in laboratory studies. The technique employs an IMU embedded in the grip of the pole and companion algorithms that estimate the instantaneous lean angle of the pole while also identifying the instants of pole plant and pole lift. An experiment (using 20 trials with a total of 40 poling cycles) reveals that the IMU-based estimates of the pole lean angle falls within 95% confidence intervals and that they remain correlated ( r = 0.99 ) with those estimated from video images. Similarly, the IMU-based estimates of the instants of pole plant and pole lift exhibit little difference (− 4 ± 15 ms for plant and − 7 ± 11 ms for lift) with those estimated from the video images. Consequently, this technology holds considerable promise for measuring and evaluating poling performance during free skiing on outdoor courses.