Next Article in Journal
Reliability of Vibroarthrography to Assess Knee Joint Sounds in Motion
Next Article in Special Issue
Inertial Sensor-Based Motion Tracking in Football with Movement Intensity Quantification
Previous Article in Journal
Accuracy Improvement of Binocular Vision Measurement System for Slope Deformation Monitoring
Previous Article in Special Issue
Timing and Modulation of Activity in the Lower Limb Muscles During Indoor Rowing: What Are the Key Muscles to Target in FES-Rowing Protocols?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ski Jumping Trajectory Reconstruction Using Wearable Sensors via Extended Rauch-Tung-Striebel Smoother with State Constraints

1
Institute of Flight System Dynamics, Technical University of Munich, 85748 Garching, Germany
2
Department of Biomechanics in Sports, Faculty of Sport and Health Sciences, Technical University of Munich, 80992 Munich, Germany
3
Olympic Training Center of Bavaria, 80809 Munich, Germany
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(7), 1995; https://doi.org/10.3390/s20071995
Submission received: 13 January 2020 / Revised: 28 March 2020 / Accepted: 30 March 2020 / Published: 2 April 2020

Abstract

:
To satisfy an increasing demand to reconstruct an athlete’s motion for performance analysis, this paper proposes a new method for reconstructing the position and velocity in the context of ski jumping trajectories. Therefore, state-of-the-art wearable sensors, including an inertial measurement unit, a magnetometer, and a GPS logger are used. The method employs an extended Rauch-Tung-Striebel smoother with state constraints to estimate state information offline from recorded raw measurements. In comparison to the classic inertial navigation system and GPS integration solution, the proposed method includes additional geometric shape information of the ski jumping hill, which are modeled as soft constraints and embedded into the estimation framework to improve the position and velocity estimation accuracy. Results for both simulated measurement data and real measurement data demonstrate the effectiveness of the proposed method. Moreover, a comparison between jump lengths obtained from the proposed method and video recordings shows the relative root-mean-square error of the reconstructed jump length is below 1.5 m depicting the accuracy of the algorithm.

1. Introduction

Acquiring an accurate estimation of the position and velocity of the athlete during a ski jump has always been of great interest to athletes, coaches, and sport researchers to analyze the movement for improving the jumping performance. Traditional solutions to meet this demand are video analysis techniques based on camera recorded videos [1]. However, such a system has the disadvantage that it generally only covers limited parts of the jumping area [2]. To record the entire jump, a large number of cameras would be required making the system expensive and hard to calibrate. Furthermore, camera systems merely provide information about the athlete’s position, while velocity and acceleration have to be estimated via numeric differentiation, which may be subject to large error.
In recent years, wearable sensors, such as inertial measurement units (IMUs) and global navigation satellite system (GNSS) receivers, have been widely employed for motion analysis in sports including ski jumping [3]. Chardonnens et al. [4] presented an IMU-based system measuring ski jumping dynamics including the position and velocity of the center of mass perpendicular to the table during take-off. Groh et al. [2] proposed a measurement system with IMUs and a light barrier to estimate the ski velocity and jump length. Blumenbach [5] designed and tested a helmet with a high-precision Global Positioning System (GPS) receiver for the positioning of ski jumpers. Fasel et al. [6] presented a solution to estimate the center of mass and its velocity using differential GNSS and IMUs for alpine skiing. In comparison to video analysis techniques, trajectory estimation using wearable IMUs and GNSS sensors has the advantages of being easy to use and maintain and having full coverage over all phases in a ski jump. For coaches and athletes, such a system provides the possibility to give feedback after each jump regarding the jumper’s performance, such as jump length and take-off velocity data.
Trajectory reconstruction techniques are a class of methods to achieve an accurate state estimation of a moving object by properly combining information from the kinematic model with measurement data from sensors, including but not limited to IMUs, magnetometers, and GNSS receivers. Trajectory reconstruction techniques are also termed flight path reconstruction in the aerospace field. A systematic overview of such methods for an aircraft can be found in [7]. Göttlicher and Holzapfel [8] reconstructed an aircraft’s flight path using low-cost sensors with the extended Rauch-Tung-Striebel (RTS) smoother. Similar techniques can be also applied to ski jumping applications.
From a methodological point of view, a trajectory reconstruction problem can be interpreted as a state estimation problem: One powerful tool to solve this is the classic RTS smoother, which was first proposed for linear systems in [9] using a linear Kalman filter [10] as the forward-time estimator. The extended RTS smoother is a modification of the classical RTS smoother capable of considering nonlinear systems. It employs the extended Kalman filter (EKF) as the forward-time filter to deal with the nonlinearity of the system. In this paper, the extended RTS smoother is adopted as an estimator. Moreover, we can formulate state equality constraints in a state estimation problem. Different methods to solve state estimation problems with state constraints, which are termed as constrained filtering (or smoothing) problems, are discussed in the survey [11].
In this paper, we present an offline method based on the extended RTS smoother to estimate the trajectory, including the position and velocity, of a ski jumper by the inertial navigation system and GPS integration. Additionally, to fully utilize all available information to improve the estimation quality, the geometric shape of the ski jumping hill is modeled as a set of additional soft constraints and included in the smoothing framework.

2. Data Acquisition

The measurement data were collected during a summer session in June 2018 on the jumping hill Schattenbergschanze (Hill Size = 106 m) in Oberstdorf, Germany. During this summer session, the hill slope was covered with water-soaked plastic. The measurement setup consists of the following sensors: a Qstarz BT Q1000eX GPS data logger and an iPhone 7, which contains an InvenSense ICM-20600 IMU, and an Alps HSCDTD008A magnetometer. Table 1 summarizes the performance characteristics of the sensors provided by the respective manufacturers. The GPS data logger provides position and velocity measurements with a 10 Hz sample rate. The IMU measures triaxial rotational rates as well as specific forces (i.e., non-gravitational accelerations) at 100 Hz frequency. The magnetometer provides triaxial local magnetic field flux density measurements also at 100 Hz. To the best of our knowledge, the IMU and magnetometer sensor chips in the iPhone 7 provide raw measurement data with a similar accuracy level as previous studies, e.g., [2,4,6,8]. The raw data measured by the smartphone was logged via the application Phyphox [12] and was attached to the jumper’s right upper arm using a running armband with the orientation shown in Figure 1. The GPS logger was fixed to the top of the helmet by adhesive tapes for better signal reception.

3. Methods

In this section, we first introduce the extended RTS smoother framework for a general system in Section 3.1. Next, modifications made for the joint state and parameter estimation as well as the state constraints in Section 3.2. Finally, the system model for ski jumping, which is implemented in the extended RTS smoother framework, is presented in Section 3.3.

3.1. Extended Rauch-Tung-Striebel Smoother

The extended RTS smoother is an offline estimation method utilizing the entire batch of measurements over a fixed time interval [16]. It combines a forward filtering pass using the EKF with a backward smoothing pass with the RTS smoother.
A state-space representation for a general nonlinear system described by ordinary differential equations (ODEs) can be written as
(1) x ˙ = f ( x , u , w ) , x ( t 0 ) = x 0 , (2) y = g ( x ) , (3) z = y + v ,
where x R n x is the state vector, u R n u is the known input vector, y R n y is the model outputs, and z R n y represents the measurements. Symbols w R n w and v R n y denote process and measurement noise vectors. This system can be linearized and discretized into the following form [16]:
(4) x k = Φ k 1 x k 1 + Γ u , k 1 u k 1 + Γ k 1 w k 1 , (5) y k = H k x k , (6) z k = y k + v k ,
where H k = g ( x k ) / x . The transformation matrices Φ k 1 , Γ u , k 1 , and Γ k 1 are computed as
Φ k 1 = exp ( A k 1 Δ t k 1 ) = i = 0 ( A k 1 Δ t k 1 ) i i ! ,
Γ u , k 1 = 0 Δ t k 1 exp ( A k 1 τ ) B k 1 d τ = i = 1 A k 1 i 1 Δ t k 1 i i ! B k 1 ,
Γ k 1 = 0 Δ t k 1 exp ( A k 1 τ ) F k 1 d τ = i = 1 A k 1 i 1 Δ t k 1 i i ! F k 1 ,
where A k 1 = f ( x k 1 , u k 1 , w k 1 ) x , B k 1 = f ( x k 1 , u k 1 , w k 1 ) u , F k 1 = f ( x k 1 , u k 1 , w k 1 ) w , and Δ t k 1 = t k t k 1 . The noise processes { w k } and { v k } are assumed to be zero-mean, uncorrelated, white Gaussians, i.e., w k N ( 0 , Q k ) , v k N ( 0 , R k ) , E [ w j w k ] = Q k δ j k , E [ v j v k ] = R k δ j k , and E [ w j v k ] = 0 , where δ j k is the Kronecker delta.

3.1.1. Extended Kalman Filter

The EKF is a sequential state estimation method for nonlinear systems [11,16]. First, the initial state x ^ 0 | 0 and covariance P 0 | 0 for the EKF are given as
x ^ 0 | 0 = x ¯ 0 ,
P 0 | 0 = P ¯ 0 ,
where x ¯ 0 is an a priori estimate of the initial states x 0 and P ¯ 0 is the covariance of the estimation error of x ¯ 0 . Then, for each time point t k with k { 1 , 2 , , N } , the EKF consist of a prediction and an update step as follows:
1) Prediction step:
x ^ k | k 1 = x ^ k 1 | k 1 + t k 1 t k f ( x ^ k 1 | k 1 , u ˜ k 1 , 0 ) d τ ,
P k | k 1 = Φ k 1 P k 1 | k 1 Φ k 1 + Γ k 1 Q k 1 Γ k 1 ,
where x ^ k | k 1 represents the a priori state estimate at time t k based on measurements up to time t k 1 , i.e., x ^ k | k 1 = E ( x k | z 1 , , z k 1 ) , and P k | k 1 denotes the covariance of the estimation error of x ^ k | k 1 . To evaluate Φ k 1 and Γ k 1 , matrices A k 1 and F k 1 in Equation (7) and Equation (9) are linearized as A k 1 = f ( x ^ k 1 | k 1 , u k 1 , 0 ) x and F k 1 = f ( x ^ k 1 | k 1 , u k 1 , 0 ) w . The integral in Equation (12) is numerically solved via the classic 4-th order Runge–Kutta method with u ˜ k 1 = t k τ t k t k 1 u k 1 + τ t k 1 t k t k 1 u k .
2) Update step:
K k = P k | k 1 H k [ H k P k | k 1 H k + R k ] 1 ,
x ^ k | k = x ^ k | k 1 + K k [ z k g ( x ^ k | k 1 ) ] ,
P k | k = [ I n x K k H k ] P k | k 1 ,
where x ^ k | k denotes the a posteriori state estimate at time t k based on measurements up to time t k , i.e., x ^ k | k = E ( x k | z 1 , , z k ) , and P k | k denotes the covariance of the estimation error of x ^ k | k . The symbol I n x denotes an n x by n x identity matrix, and H k is approximated by H k = g ( x ^ k | k 1 ) / x .

3.1.2. Rauch-Tung-Striebel Smoother

After the forward-time filtering pass, the EKF estimation of the states x ^ N | N and its error covariance P N | N is used to initialize the corresponding values in the RTS smoother for the final time point t N :
x ^ N s = x ^ N | N ,
P N s = P N | N ,
where the superscripts “s” indicate the smoothed estimates. Then, the RTS smoother runs backward in time for t k with k { N 1 , N 2 , , 0 } as
K k s = P k | k Φ k ( P k + 1 | k ) 1 ,
x ^ k s = x ^ k | k + K k s [ x ^ k + 1 s x ^ k + 1 | k ] ,
P k s = P k | k + K k s [ P k + 1 s P k + 1 | k ] ( K k s ) ,
where the smoothed estimate x ^ k s can be regarded as the state estimate of time point t k utilizing all measurements, i.e., x ^ k s = E ( x k | z 1 , , z N ) [17], and P k s is the covariance of the estimation error of x ^ k s .

3.2. Adaption for Parameter Estimation and Soft Constraints

To account for unknown parameters and known constraints for the estimation problem, the extended RTS smoother for state estimation is modified by the following adaptations.

3.2.1. Joint State and Parameter Estimation

When the model also includes parameters p R n p to be estimated, these parameters are treated as additional states with p ˙ = 0 so that they are compatible with the extended RTS smoother framework. Therefore, with the augmented state vector defined as x a = x p R n x + n p , the augmented system ODEs become
x ˙ a = f a ( x a , u , w ) = f ( x , u , p , w ) 0 .

3.2.2. Constrained Filtering

Suppose that a system satisfies the nonlinear state equality constraints as
c k = h ( x k ) .
In order to include these into the extended RTS smoother framework, the constraints can be treated as additional pseudo measurements in the form of
z k c k = y k h ( x k ) + v k v c , k .
For the case that equality constraints are satisfied exactly, which is termed as hard constraints in [11,18], the pseudo measurements can be treated as perfect with pseudo measurement noises v c , k = 0 . Contrary, soft constraints are only required to be fulfilled approximately, i.e., c k h ( x k ) . In this case, we assume v c , k are zero-mean Gaussian noises with small covariance Q c , i.e., v c , k N ( 0 , Q c ) .
We further write Equation (24) as
y a = g a ( x a ) ,
z a = y a + v a ,
where z a = z k c k , y a = y k h ( x k ) , and v a = v k v c , k . By Equations (22), (25), and (26), we formulated the state-space presentation of an augmented system for state and parameter estimation with state constraints in the form of Equations (1)–(3).

3.3. System Model

After introducing the estimation framework for a general system in the previous two subsections, the system model for the ski jumping application is presented in this subsection.

3.3.1. Coordinate Frame Definitions

The first necessary frame is the hill reference frame O x N y N z N . As shown in Figure 2, its origin O is located at the end of the in-run table. The x N axis and z N axis are both within the symmetric plane of the jumping hill and point horizontally forward and vertically down, respectively. The y N axis is perpendicular to x N and z N axes to form a right-handed coordinate system. r N = x N y N z N denotes the coordinate matrix in the hill reference frame of the jumper’s (or, to be exact, the IMU sensor’s) relative position vector r with respect to the origin O. Another necessary frame is the local north-east-down (NED) frame O x O y O z O . The origin of the local NED frame coincides with point O, and the three axes x O , y O and z O point to north, east, and down directions respectively. Finally, the sensor’s body frame O B x B y B z B is defined as following: the origin O B is located at the center of the cellphone, and the three axes are fixed to its three axes of symmetry as shown in Figure 3. The coordinates in the NED frame of the relative position vector r are denote by r O = x O y O z O .

3.3.2. Conversion of GPS Position Measurements to the Hill Reference Frame

The position measurements from GPS receivers are usually provided in geodetic coordinates, such as the World Geodetic System 1984 (WGS84) frame [19], as a time sequence of latitude μ , longitude λ , and height h. For ski jumping applications, it is more useful to calculate relative positions of the ski jumper with respect to the jumping hill.
With the frame definitions stated before, the GPS-measured positions r WGS , GPS = μ λ h in the WGS84 frame are first transformed into the local NED frame r O , GPS = x O y O z O by the following equations:
(27) x O = ( M μ + h ) ( μ μ 0 ) , (28) y O = ( N μ + h ) cos μ 0 ( λ λ 0 ) , (29) z O = ( h h 0 ) ,
where μ 0 , λ 0 , and h 0 are the latitude, longitude, as well as height of the hill reference frame’s origin O, and parameters for the WGS84 reference ellipsoid [19] can be found in Table 2.
Then, the GPS-measured positions r N , GPS in the hill reference frame are obtained by
r N , GPS = M NO r O , GPS ,
where M NO is the transformation matrix from local NED to hill reference frame as
M NO = cos ψ N sin ψ N 0 sin ψ N cos ψ N 0 0 0 1 ,
where ψ N denotes the hill azimuth angle, which is defined as the rotation angle from the x O axis (the north) to the x N axis, as shown in Figure 2.

3.3.3. Constraint Modeling of the Geometric Shape of a Ski Jumping Hill

The geometric shape of the ski jumping hill is used as additional a priori information and included in the trajectory reconstruction problem as a set of constraints for further improvement of the estimation quality. The mathematical model of a ski jumping hill is described in [20], which was published by the International Ski Federation (Fédération Internationale de Ski, FIS) as the guideline to construct ski jumping hills. Additionally, all necessary model parameters can be found in the jumping hill certificate, which is also issued by the FIS. In our case, we have the certificate of jumping hill - No. 215/GER 29 [21] for the Schattenbergschanze.
Figure 4 shows the geometric profile of a FIS-certified ski jumping hill in the longitudinal plane. The curve from the point A to the origin O is called the in-run where the jumper skis down following the track before take-off. For any point P 1 on the in-run curve, we assume its coordinates to be ( x P 1 , 0 , z IR ( x P 1 ) ) where z P 1 = z IR ( x P 1 ) denotes the mapping relationship for the in-run in the vertical plane. Since the trajectory of the jumper before take-off must coincide with the in-run curve, we define the constraints for the in-run curve in the longitudinal and lateral plane as
c V , IR = z N z IR ( x N ) ,
c H , IR = y N .
When the jumper is in the in-run area before take-off, i.e., t k t T O , the constraints should fulfill that c V , IR 0 and c H , IR 0 . The curve BC in Figure 4 is the landing area (including the landing slope and the out-run). For a point in the landing area P 2 , its position coordinate in the hill reference frame is ( x P 2 , y P 2 , z LA ( x P 2 ) ) where z P 2 = z LA ( x P 2 ) is the mapping relationship of the landing area in vertical plane. The jumper would land somewhere in between and afterward skiing along this surface. Therefore, we can formulate another constraint as follows:
c V , LA = z N z LA ( x N ) .
After the jumper’s landing, i.e., t k t T D , the constraints should be active such that c V , LA 0 . The method for detecting the take-off time point t T O and the touch-down time point t T D using raw measurements is introduced in Appendix A.

3.3.4. Measurement Error Models

To model the measurement error of the inertial sensor, constant bias terms Δ ω and Δ a for the gyroscope and accelerometer as well as zero-mean, white, Gaussian noise terms w gyro and w acc are explicitly considered as
ω B , gyro = ω B + Δ ω + w gyro ,
a B , acc = a B + Δ a + w acc ,
where ω B , gyro represents the rotational rates measurement from the gyroscope, a B , acc denotes the measurement from the accelerometer, and ω B as well as a B are the true rotational rates and specific forces respectively.
The measurement error model of the magnetometer is considered to be
m B , mag = ( 1 + diag ( Δ S m ) ) m B + Δ m + v mag ,
where m B , mag is the magnetometer measurement of the local magnetic field strength, Δ S m represents the vector of the scaling factor error and diag ( Δ S m ) denotes the corresponding square matrix with S m being its diagonal elements, Δ m denotes the constant bias vector, m B is the true local magnetic field strength, and v mag represents the zero-mean, white, Gaussian measurement noise.
The GPS measurements are modeled as
r N , GPS = r N + v pos ,
v O , GPS = v O + v vel ,
where r N is the true position in the hill reference frame, v O is the true velocity in the local NED frame, r N , GPS and v O , GPS are the GPS measurements accordingly, and v pos and v vel are the measurement noise terms.

3.3.5. State, Output, and Constraint Equations

The system states vector x is given by
x = r N v B q BO ,
where v B denotes the triaxial velocity in the body-fixed frame, and q BO represents the attitude quaternions of the body-fixed frame with respect to the local NED frame. The input vector u contains the IMU measurements ω B , gyro and a B , acc as
u = ω B , gyro a B , acc .
The parameter vector p to be estimated includes the biases of the gyroscope Δ ω , the accelerometer Δ a , the magnetometer Δ m , and the scaling factor error of the magnetometer Δ S m :
p = Δ ω Δ a Δ m Δ S m .
The process noise vector w consists of influences from both gyroscope and accelerometer:
w = w B , gyro w B , acc .
With the definitions above, the state equations x ˙ = f ( x , u , p , w ) stated in Equation (22) can be assembled by the equations of inertial navigation mechanism and kinematics as follows:
r ˙ N = M NO M BO v B ,
v ˙ B = a B + M BO g O + ω B × v B ,
q ˙ BO = 1 2 q 1 q 2 q 3 q 0 q 3 q 2 q 3 q 0 q 1 q 2 q 1 q 0 ω B ,
where M BO is the transformation matrix from the local NED frame to the body-fixed frame as
M BO = q 0 2 + q 1 2 q 2 2 q 3 2 2 q 1 q 2 + q 0 q 3 2 q 1 q 3 q 0 q 2 2 q 1 q 2 q 0 q 3 q 0 2 q 1 2 + q 2 2 q 3 2 2 q 2 q 3 + q 0 q 1 2 q 0 q 2 + q 1 q 3 2 q 2 q 3 q 0 q 1 q 0 2 q 1 2 q 2 2 + q 3 2 ,
and g O = 0 0 g 0 represents the gravitational vector in the local NED frame. The term a B in Equation (45) can be substituted by a B = a B , acc Δ a w acc using the accelerometer measurement model Equation (36). Similarly, ω B can be written as ω B = ω B , gyro Δ ω w gyro via the gyroscope measurement model Equation (35) and be further substituted into Equation (46).
The output vector y is written as
y = r N v O m B , meas ,
where r N can be directly obtained from the state vector, the velocity v O is calculated as
v O = M BO v B ,
and m B , meas is the estimated magnetometer measurement. By taking the magnetometer model Equation (37) into account, m B , meas is estimated by
m B , meas = ( 1 + diag ( Δ S m ) ) M BO m O , WMM + Δ m ,
where m O , WMM is the reference value of the local magnetic field strength in the local NED frame from the World Magnetic Model 2015 [22].
Accordingly, the measurement vector is z = r N , GPS v O , GPS m B , mag with measurements from the GPS logger ( r N , GPS and v O , GPS ) and the magnetometer ( m B , mag ). The measurement noise vector is therefore v = v pos v vel v mag .
The vector of constraints c is defined to be
c = c q c V , IR c H , IR c V , LA ,
where c q = q BO 2 1 is the constraint for attitude quaternions. Here, c q ensures a unit quaternion, i.e., q BO 2 = 1 . The constraints c V , IR , c H , IR , and c V , LA are geometric constraints given by the shape of the jumping hill, which are introduced in Section 3.3.3.

4. Results and Discussion

In this Section, trajectory reconstruction results from both simulated and real measurement data are presented. The configuration for the extended RTS smoother is first introduced in Section 4.1. The simulated measurements are used to theoretically validate the purposed method in Section 4.2. Results based on real measurements are presented in Section 4.3. The jump length obtained from both the proposed method and video recordings are compared in Section 4.4 to further validate the results based on real measurement data.

4.1. Setting

We first introduce necessary configurations for the extended RTS smoother algorithm. For the initial states x ¯ 0 = r ¯ N , 0 v ¯ O , 0 q ¯ 0 , the GPS measurement at initial time t 0 is adopted as an estimation for the initial position and velocity as r ¯ N , 0 = r N , GPS , 0 and v ¯ O , 0 = v O , GPS , 0 . The estimated value for initial attitude q ¯ 0 is calculated via the factored quaternion algorithm proposed in [23] using the accelerometer and magnetometer measurements at time t 0 . The initial guess for all measurement error parameters is chosen to be zero as p ¯ 0 = Δ ω 0 Δ a 0 Δ m 0 Δ S m , 0 = 0 . For the augmented initial states x ¯ a , 0 = x ¯ 0 p ¯ 0 , the estimated covariance matrix for the initial states error P ¯ 0 is set as shown in Table 3. In addition, the covariance matrix of noises Q and R is considered constant and the estimated values of Q and R are listed in Table 4 and Table 5 respectively. Within the model, the hill azimuth angle ψ N is set as 312.6° for the jumping hill Schattenbergschanze. The gravitational constant at the site is calculated to be g 0 = 9.8053 m / s 2 according to the Earth Gravitational Model 1996 [19].

4.2. Validation by Simulated Measurement Data

To validate the accuracy of the proposed algorithm, an artificial measurement data set is generated by assuming a true value and adding sensor error and noise. This data set is generated based on the result from real measurements as a reference trajectory, ensuring similar observability of the estimation problem as in the real world. Trajectory optimization techniques with least-squares costs are used to obtain a trajectory with similar acceleration and angular velocities as the measured data. By adding path constraints to the trajectory optimization, the true trajectory also fulfills the constraints in Equations (32)–(34). The true value for the generated trajectory also complies with the kinematic models in Equations (44)–(46). The sensor errors and noises are added to the true values according to Equations (35)–(39), with the noise covariance values shown in Table 4 and Table 5 to generate artificial measurements.
The trajectory is reconstructed with only the noisy artificial measurements known to the extended RTS smoother. The reconstruction result is presented in Figure 5. Here, it can be observed that the smoothed trajectory agrees closely with the true reference. To quantify the error of the estimated position and velocity, we calculated the error by
Δ r N = r ^ N s r N , ref ,
Δ v O = v ^ O s v O , ref ,
where variables with subscript “ref” are the true values from the generated data. The error of these estimated positions and velocities are presented in Figure 6a,b respectively. To intuitively display the attitude, the quaternions q BO are converted to the corresponding Euler angles (in z-y-x order): the bank angle ϕ , the pitch angle θ , and the azimuth angle ψ by
(54) ϕ = arctan 2 q 2 q 3 + q 0 q 1 q 0 2 q 1 2 q 2 2 + q 3 2 , (55) θ = arcsin 2 q 1 q 3 q 0 q 2 , (56) ψ = arctan 2 q 1 q 2 + q 0 q 3 q 0 2 + q 1 2 q 2 2 q 3 2 .
Then, we use the error in Euler angles to represent the attitude estimation error as
Δ ϕ = ϕ ^ s ϕ ref , Δ θ = θ ^ s θ ref , and Δ ψ = ψ ^ s ψ ref .
The Euler angle estimation error is presented in Figure 6c. The root-mean-square (RMS) error for position, velocities, and Euler angles are summarized in Table 6. This result shows that the proposed algorithm can successfully reconstruct the trajectory and that the position accuracy for the assumed noise covariance is in decimeter magnitudes.

4.3. Validation by Real Measurement Data

Now, we apply the extended RTS smoother on real measurement data to the ski jumping trajectory reconstruction problem. The results for a reconstructed ski jump are presented in the following part.
In Figure 7, the three-dimensional trajectories of the jumper are plotted with a digital ski jumping hill model based on [20,21]. The purple solid line with dots shows the converted raw GPS position measurements r N , GPS where each dot on the line represents a measurement sample. Furthermore, the line consisting of color-coded dots represents position reconstructed by the extended RTS smoother, where the color of one dot represent the velocity V s = v O 2 = v O , x 2 + v O , y 2 + v O , z 2 . Note that we use the relative local height h N here instead of z N in the figure for convenience. The relationship between h N and z N is directly h N = z N . The subplot on the top-right corner of Figure 7 presents the projection of the jumping trajectories on the x N - z N plane (focusing mainly on the flight part). The comparison between GPS-measured and extended RTS smoother reconstructed position as well as velocity time histories is presented in Figure 8 and Figure 9 respectively.
Both GPS measurements and the smoother estimation in Figure 7 show the trajectory of an entire jump from start to stop. By comparing the GPS-measured trajectory to the jumping hill model for the in-run and out-run parts, it shows that the GPS height channel contains relatively large errors. After the jumper starts his run along the track during the in-run phase, the GPS-measured trajectory indicates that the jumper starts to move horizontally, but remains at the same height vertically. For part after the touch-down in the landing area, the GPS-measured trajectory shows large error leading to the position being above the ground. Therefore, in this case, the GPS-measured trajectory is not adequate for the demand for analyzing the ski jump limited to its accuracy. It is worth mentioning that the GPS measured trajectory does not necessarily encounter the same type of error in the height channel as in our case. Nevertheless, due to the trilateration positioning principle of the GPS technology employed, the GPS height measurement is far less accurate than the latitude and longitude measurements. The GPS service standard [24] reported that the GPS has a global average positioning accuracy of 9 m ( 95 % ) horizontal error and 15 m ( 95 % ) vertical error.
On the other hand, besides the GPS measurements as a data source, the extended RTS smoother uses IMU and magnetic measurements as well as geometric information of the jumping hill to perform data fusion. The application of the RTS smoother reconstructs the trajectory despite a lower GPS measurement frequency and the existence of some GPS height errors as shown in Figure 7, Figure 8 and Figure 9. In Figure 7, the reconstructed trajectory shows the entire process of the jump. The jumper first slides down along the in-run track with increasing speed V S . For the flight part, this can be clearly observed from both the three-dimensional plot and the subplot in Figure 7. It can be observed that the height above the ground is decreasing until the jumper’s touch-down with the speed further increasing. After the touch-down, the jumper skis down on the landing area of the hill with his speed still accelerating. After reaching the flat part at the hill bottom, the jumper decelerates due to the friction and his break action and then comes to a full stop. In comparison to the GPS measurement, the estimated position and velocity are improved especially in the height channel, which reveals more information for the jump analysis, and the sampling frequency on the trajectory is also increased.
In Figure 8, we can observe that two lines fit well for x N . For y N in the extended RTS smoother result (the green line), the constraint c H , IR correct the position y N to 0 before take-off. Afterward, the reconstructed trajectory matches the GPS measurements with some differences. For the h N channel, as also mentioned before, the raw GPS measurement has a relatively large error which can be clearly observed in Figure 7. The smoother reconstructed height by considering the vertical constraints c V , I R and c V , L A is moving in the in-run curve before the take-off and on the surface of the landing area after the touch down (see Figure 7). The height history is closer to a real physical movement. In Figure 9, for horizontal velocity v O , x and v O , y , the reconstruction and measurement fit accurately. The extended RTS smoother result gives a bit different estimation of v O , z compared to the GPS measured one. Considering the fact that after the jumper starts gliding down around 610.5 s, the vertical velocity afterward should be positive (moving downwards). So, the vertical speed measured by GPS also contains an error for the in-run part which is associated with the height error. Since the extended RTS smoother considers kinetics h ˙ N = v O , z , the smoother estimated vertical speed is linked to the height estimation. Combining the result in Figure 7 and Figure 8 for the in-run part (around 610.5 s to 616.8 s), the positive v O , z estimates the smoother giving out is more close to the fact that jumper is moving vertically down in the in-run. This indicates that the extended RTS smoother improves the velocity comparing to the GPS with the assistance of information from other sources.
The improvement of the estimation accuracy is achieved by fusing data from different sources including measurement sensors, kinematics, and geometric information of the hill. The gyroscope and accelerometer data, after removing the estimated bias error, are used for propagating the states including attitude, velocity, and position together with a priori knowledge of the uncertainties. The GPS measurements and the magnetometer are included as additional data sources for the correction step. Furthermore, the introduction of the geometric constraints from the jumping hill model provides an additional information source, which contributes to the estimation. The constraints c V , IR and c V , LA , as shown in Figure 10, establish a strong relationship between states x N and z N on the vertically plane for each part of the trajectory (before take-off and after touchdown).
Although the major aim of applying the extended RTS smoother is to obtain a better position and velocity estimation, it is worth to mention that the RTS smoother also provides an attitude estimation of the IMU sensor. Furthermore, the estimation of the IMU attitude is essential to attribute the accelerometer measured specific forces into the correct direction (see Equations (45) and (49) ). By applying Equations (54)–(56), the converted Euler angles reconstructed by the extended RTS smoother are shown in Figure 11. From Equation (37), we know that the magnetic field measurements change with the attitude. Therefore, the similar tendency of extended RTS smoother estimated m B , meas and magnetometer measurements m B , mag in Figure 12 indicates good attitude estimation.
To conclude the result presentation, measurements from the gyroscope and the accelerometer are shown in Figure 13 and Figure 14.
From the measurement campaign, we collected valid data sets for five different jumps in total. The reconstruction results for the rest of the four jumps are shown in Figure 15, Figure 16, Figure 17 and Figure 18. Comparing to the result from Jump No. 1 (Figure 8), the results show similar characteristics as before which demonstrated the repeatability of the results.

4.4. Validation by Jump Length

To validate the trajectory reconstruction results, we calculate the jump length results for the five collected trials obtained from the trajectory reconstruction L TR for comparison with the jump length measurements from video recordings L VR . The process of determining the jump length from the video recordings is explained in Appendix B. The validation result of five different jumps is shown in Table 7.
The error Δ L in Table 7 is defined as Δ L = L TR L VR , and the RMS error for the jump length estimation is therefore Δ L ¯ = 1 5 i = 1 5 Δ L 2 = 1.0 m . Considering that the estimate of the L VR value could contain an error of 0.5 m as worst case, we calculate the upper bound of the RMS error as Δ L ¯ ub = 1 5 i = 1 5 ( | Δ L | + 0.5 ) 2 = 1.5 m .

5. Conclusions and Discussion

This paper presents a novel method for trajectory reconstruction in ski jumping using wearable sensors by applying the extended RTS smoother with state constraints. The result of simulated measurement data validates the proposed method and shows that it has a theoretical accuracy of decimeter magnitude in positions. The results based on real measurements show the method can successfully reconstruct the trajectories of ski jumps. By formulating the geometric information of the jumping hill as soft equality constraints, the estimation accuracy of positions and velocities is largely improved comparing to the GPS measurements. Jump length data for the five collected trails obtained from video recordings demonstrate good accordance with the reconstruction results with the root-mean-square error upper bound of 1.5 m . This provides strong supporting evidence for the validity of the proposed method.
For future research, a quantitative assessment of the accuracy of the proposed method by further real measurement data needs to be carried out. Therefore, another measuring system providing higher accuracy of position and velocity, e.g., a well-calibrated camera system, needs to be employed in order to compare the results of the two systems. For future device developments, wearable sensors providing better position and velocity measurement accuracy, such as differential GPS receivers, which could provide centimeter-level accuracy [6], are recommended to gain better estimation of the jumper’s trajectories. In addition, GPS receivers with higher sampling frequency can be tested which could potentially improve the results. On the other hand, due to the importance of size and weight in the performance of ski jumping, the wearability of such sensors is still a major aspect to be considered. Future measurements should include more trials performed by different athletes in order to increase the statistical confidence. After further development, this technique is promising to be implemented to provide quantified motion feedback in ski jumping training.

Author Contributions

Conceptualization, X.F. and F.H.; methodology, X.F.; software, X.F.; validation, X.F. and B.G.; formal analysis, X.F.; investigation, X.F., B.G., P.P., V.B. and J.P.; resources, F.H.; data curation, X.F.; writing–original draft preparation, X.F.; writing–review and editing, B.G., P.P., V.B., J.P. and F.H.; visualization, X.F.; supervision, F.H.; project administration, B.G., funding acquisition, F.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Deutsche Forschungsgemeinschaft (DFG) through the TUM International Graduate School of Science and Engineering (IGSSE), GSC 81. The first author was funded by the China Scholarship Council(CSC) for his doctoral study in TUM.

Acknowledgments

The authors are thankful to the athletes who participated in the study and to the availability of the coach. The authors would also like to thank Ansgar Schwirtz and Matthias Bittner for their support in the TUM IGSSE research project for ski jumping.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the result.

Appendix A. Take-Off and Touch-Down Time Point Detection

The method implemented for detecting the time points for take-off t T O and touch-down t T D is introduced in this section: Due to the change of the ground reaction forces, both take-off and touch-down phase show abrupt changes in acceleration, which is measured by the accelerometer. Previous studies on ski jumping take-off [25] and landing [26,27] also presented similar acceleration changing behavior in their measurement-based results. Therefore, we propose to determine t T O and t T D by detecting these abrupt changes in the accelerometer measurement.
First of all, since the take-off and touchdown would happen when the jumper is moving with a high speed, we limit the analyzing time interval to the part when the GPS measured speed V GPS ( t k ) = | | v GPS ( t k ) | | 2 = u GPS 2 ( t k ) + v GPS 2 ( t k ) + w GPS 2 ( t k ) is greater than 10 m / s . Next, within this time interval, the L2 norm of the triaxial specific force measured by accelerometer for discrete time points t k are calculated by a norm ( t k ) = | | a B , acc ( t k ) | | 2 = a x , acc 2 ( t k ) + a y , acc 2 ( t k ) + a z , acc 2 ( t k ) as shown in Figure A1. Then, two peaks (local maximums) with very large a norm and time difference more than 1.5 s (based on experience) are selected by the following rule as shown in Figure A1 as circles:
t p e a k 1 = argmax t k ( a norm ( t k ) | V G P S ( t k ) > 10 m / s ) ,
t p e a k 2 = argmax t k ( a norm ( t k ) | V G P S ( t k ) > 10 m / s , 1.5 s t k t p e a k 1 1.5 s ) .
These two peaks indicate large ground reaction forces for take-off and landing. Finally, the take-off time points t T O is determined as the earlier peak. As presented in [26,27], the actual touch-down point in the landing process is a short moment before the maximal ground reaction peak, approximately 0.05 s as we estimated. Therefore, the touch-down time points t T D is determined as 0.05 s before the later peak as:
t T O = min { t p e a k 1 , t p e a k 2 } ,
t T D = max { t p e a k 1 , t p e a k 2 } 0.05 s .
Figure A1. An illustrative example on detecting the take-off and touch-down time points by the raw measurement data.
Figure A1. An illustrative example on detecting the take-off and touch-down time points by the raw measurement data.
Sensors 20 01995 g0a1

Appendix B. Jump Length Detected by the Video Recordings

In this section, the process for detecting the jump length detected using the video recordings is explained. Cameras with a frame rate of 50 fps are set up at the side of the landing area to record the jumper’s landing movement. Figure A2 shows a picture of overlaid snapshots from one of the video recordings. The time difference between each pair of neighboring snapshots is 0.1 s (with a distance of 5 frames). First, the frame for the jumper’s touchdown determined indicated by both skis basically touching the ground. Then, the jump length is detected by the middle point of the ski with the help of the distance markers on the other side of the hill, as shown by the red dash lines in Figure A2. Overall, we consider the detected jump length is within an accuracy of 0.5 m .
Figure A2. A picture of overlaid snapshots from video recordings to illustrate the jump length detection from video recordings.
Figure A2. A picture of overlaid snapshots from video recordings to illustrate the jump length detection from video recordings.
Sensors 20 01995 g0a2

References

  1. Virmavirta, M.; Isolehto, J.; Komi, P.; Brüggemann, G.P.; Müller, E.; Schwameder, H. Characteristics of the early flight phase in the Olympic ski jumping competition. J. Biomech. 2005, 38, 2157–2163. [Google Scholar] [CrossRef] [PubMed]
  2. Groh, B.H.; Warschun, F.; Deininger, M.; Kautz, T.; Martindale, C.; Eskofier, B.M. Automated ski velocity and jump length determination in ski jumping based on unobtrusive and wearable sensors. Proc. ACM Interact. Mob. Wearable Ubiquitous Technol. 2017, 1, 53. [Google Scholar] [CrossRef]
  3. Adesida, Y.; Papi, E.; McGregor, A.H. Exploring the role of wearable technology in sport kinematics and kinetics: A systematic review. Sensors 2019, 19, 1597. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Chardonnens, J.; Favre, J.; Cuendet, F.; Gremion, G.; Aminian, K. Measurement of the dynamics in ski jumping using a wearable inertial sensor-based system. J. Sport. Sci. 2014, 32, 591–600. [Google Scholar] [CrossRef] [PubMed]
  5. Blumenbach, T. High precision kinematic GPS positioning of ski jumpers. In Proceedings of the 17th International Technical Meeting of the Satellite Division of The Institute of Navigation (ION GNSS 2004), Long Beach, CA, USA, 21–24 September 2004; pp. 761–765. [Google Scholar]
  6. Fasel, B.; Spörri, J.; Gilgien, M.; Boffi, G.; Chardonnens, J.; Müller, E.; Aminian, K. Three-dimensional body and centre of mass kinematics in alpine ski racing using differential GNSS and inertial sensors. Remote. Sens. 2016, 8, 671. [Google Scholar] [CrossRef] [Green Version]
  7. Mulder, J.; Chu, Q.; Sridhar, J.; Breeman, J.; Laban, M. Non-linear aircraft flight path reconstruction review and new advances. Prog. Aerosp. Sci. 1999, 35, 673–726. [Google Scholar] [CrossRef]
  8. Göttlicher, C.; Holzapfel, F. Flight path reconstruction for an unmanned aerial vehicle using low-cost sensors. In Proceedings of the ICAS 30th International Congress of the International Council of the Aeronautical Sciences, Daejeon, Korea, 25–30 September 2016; pp. 25–30. [Google Scholar]
  9. Rauch, H.E.; Striebel, C.; Tung, F. Maximum likelihood estimates of linear dynamic systems. AIAA J. 1965, 3, 1445–1450. [Google Scholar] [CrossRef]
  10. Kalman, R.E. A new approach to linear filtering and prediction problems. J. Basic Eng. 1960, 82, 35–45. [Google Scholar] [CrossRef] [Green Version]
  11. Simon, D. Kalman filtering with state constraints: A survey of linear and nonlinear algorithms. IET Control. Theory Appl. 2010, 4, 1303–1318. [Google Scholar] [CrossRef] [Green Version]
  12. Staacks, S.; Hütz, S.; Heinke, H.; Stampfer, C. Advanced tools for smartphone-based experiments: Phyphox. Phys. Educ. 2018, 53, 045009. [Google Scholar] [CrossRef] [Green Version]
  13. Qstarz BT Q1000eX GPS Data Logger Specifications. Available online: http://www.qstarz.com/Products/GPS%20Products/BT-Q1000EX-10HZ-S.html (accessed on 30 March 2020).
  14. InvenSense ICM-20600 6-Axis IMU Datasheet. Available online: https://www.invensense.com/download-pdf/icm-20600-datasheet/ (accessed on 30 March 2020).
  15. Datasheet Alps HSCDTD008A. Available online: https://tech.alpsalpine.com/prod/c/pdf/sensor/geomagnetic/hscd/hscdtd008a_data.pdf (accessed on 30 March 2020).
  16. Crassidis, J.L.; Junkins, J.L. Optimal Estimation of Dynamic Systems, 2nd ed.; Chapman and Hall/CRC: New York, NY, USA, 2011; ISBN 978-143-983-985-0. [Google Scholar]
  17. Särkkä, S. Bayesian Filtering and Smoothing; Cambridge University Press: Cambridge, UK, 2013; ISBN 978-110-703-065-7. [Google Scholar]
  18. Wang, L.S.; Chiang, Y.T.; Chang, F.R. Filtering method for nonlinear systems with constraints. IEE Proc. Control. Theory Appl. 2002, 149, 525–531. [Google Scholar] [CrossRef]
  19. Mularie, W. World Geodetic System 1984–Its Definition and Relationships with Local Geodetic Systems, 3rd ed.; U.S. Department of Defense: Washington, DC, USA, 2000. Available online: https://earth-info.nga.mil/GandG/publications/tr8350.2/wgs84fin.pdf (accessed on 30 March 2020).
  20. Gasser, H.H. Jumping Hills Construction Norm 2018; International Ski Federation: Oberhofen, Switzerland, 2018; Available online: https://assets.fis-ski.com/image/upload/v1542377868/fis-prod/assets/Construction-Norm_2018-2.pdf (accessed on 30 March 2020).
  21. International Ski Federation. Certificate of Jumping Hill—No. 215/ GER 29. Available online: http://www.skisprungschanzen.com/photos/ger/oberstdorf_arena/HS106_2011.jpg (accessed on 30 March 2020).
  22. Chulliat, A.; Macmillan, S.; Alken, P.; Beggan, C.; Nair, M.; Hamilton, B.; Wood, A.; Ridley, V.; Maus, S.; Thomson, A. The US/UK World Magnetic Model for 2015–2020: Technical Report; National Geophysical Data Center, NOAA: Boulder, CO, USA, 2015. [Google Scholar] [CrossRef]
  23. Yun, X.; Bachmann, E.R.; McGhee, R.B. A simplified quaternion-based algorithm for orientation estimation from earth gravity and magnetic field measurements. IEEE Trans. Instrum. Meas. 2008, 57, 638–650. [Google Scholar] [CrossRef]
  24. Grimes, J.G. Global Positioning System Standard Positioning Service Performance Standard, 4th ed.; U.S. Department of Defense: Washington, DC, USA, 2008. Available online: https://www.gps.gov/technical/ps/2008-SPS-performance-standard.pdf (accessed on 30 March 2020).
  25. Logar, G.; Munih, M. Estimation of joint forces and moments for the in-run and take-off in ski jumping based on measurements with wearable inertial sensors. Sensors 2015, 15, 11258–11276. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Bessone, V.; Petrat, J.; Schwirtz, A. Ground reaction forces and kinematics of ski jump landing using wearable sensors. Sensors 2019, 19, 2011. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Groh, B.H.; Fritz, J.; Deininger, M.; Schwameder, H.; Eskofier, B.M. Unobtrusive and wearable landing momentum estimation in Ski jumping with inertial-magnetic sensors. In Proceedings of the 2018 IEEE 15th International Conference on Wearable and Implantable Body Sensor Networks (BSN), Las Vegas, NV, USA, 4–7 March 2018; pp. 102–105. [Google Scholar] [CrossRef]
Figure 1. Illustrative figure on the placement of the sensors (based on a body scanning figure).
Figure 1. Illustrative figure on the placement of the sensors (based on a body scanning figure).
Sensors 20 01995 g001
Figure 2. Definitions of the local north-east-down frame and the hill reference frame.
Figure 2. Definitions of the local north-east-down frame and the hill reference frame.
Sensors 20 01995 g002
Figure 3. Definition of the sensor body frame.
Figure 3. Definition of the sensor body frame.
Sensors 20 01995 g003
Figure 4. The Longitudinal Profile of a Ski Jumping Hill.
Figure 4. The Longitudinal Profile of a Ski Jumping Hill.
Sensors 20 01995 g004
Figure 5. Comparison of the extended RTS smoother reconstructed trajectory (color-coded line), the generated reference trajectory (black line), and the GPS measured trajectory (purple dots) plotted on the ski jumping hill model.
Figure 5. Comparison of the extended RTS smoother reconstructed trajectory (color-coded line), the generated reference trajectory (black line), and the GPS measured trajectory (purple dots) plotted on the ski jumping hill model.
Sensors 20 01995 g005
Figure 6. The RTS smoother estimation error for the generated reference data.
Figure 6. The RTS smoother estimation error for the generated reference data.
Sensors 20 01995 g006
Figure 7. Comparison of the extended RTS smoother reconstructed trajectory (color-coded line) and the GPS measured trajectory (purple solid line with dots) plotted on the ski jumping hill model (orange lines and the gray surface).
Figure 7. Comparison of the extended RTS smoother reconstructed trajectory (color-coded line) and the GPS measured trajectory (purple solid line with dots) plotted on the ski jumping hill model (orange lines and the gray surface).
Sensors 20 01995 g007
Figure 8. Converted GPS position measurements r N , GPS (blue dots) and the extended RTS smoother reconstructed relative positions r N (green line).
Figure 8. Converted GPS position measurements r N , GPS (blue dots) and the extended RTS smoother reconstructed relative positions r N (green line).
Sensors 20 01995 g008
Figure 9. GPS velocity measurements v O , GPS (blue dots) and the extended RTS smoother reconstructed velocity v O (green line).
Figure 9. GPS velocity measurements v O , GPS (blue dots) and the extended RTS smoother reconstructed velocity v O (green line).
Sensors 20 01995 g009
Figure 10. The extended RTS smoother estimated geometric constraints c (green line) and pseudo measurements (blue dots).
Figure 10. The extended RTS smoother estimated geometric constraints c (green line) and pseudo measurements (blue dots).
Sensors 20 01995 g010
Figure 11. The extended RTS smoother estimated attitude angles: ϕ , θ , and ψ .
Figure 11. The extended RTS smoother estimated attitude angles: ϕ , θ , and ψ .
Sensors 20 01995 g011
Figure 12. Real magnetometer measurements m B , mag (blue dots) and the extended RTS smoother estimated magnetometer measurements m B , meas (green line).
Figure 12. Real magnetometer measurements m B , mag (blue dots) and the extended RTS smoother estimated magnetometer measurements m B , meas (green line).
Sensors 20 01995 g012
Figure 13. The gyroscope measurements ω B , gyro (red dots).
Figure 13. The gyroscope measurements ω B , gyro (red dots).
Sensors 20 01995 g013
Figure 14. The accelerometer measurements a B , acc (red dots).
Figure 14. The accelerometer measurements a B , acc (red dots).
Sensors 20 01995 g014
Figure 15. Trajectory Reconstruction result for Jump No. 2.
Figure 15. Trajectory Reconstruction result for Jump No. 2.
Sensors 20 01995 g015
Figure 16. Trajectory Reconstruction result for Jump No. 3.
Figure 16. Trajectory Reconstruction result for Jump No. 3.
Sensors 20 01995 g016
Figure 17. Trajectory Reconstruction result for Jump No. 4.
Figure 17. Trajectory Reconstruction result for Jump No. 4.
Sensors 20 01995 g017
Figure 18. Trajectory Reconstruction result for Jump No. 5.
Figure 18. Trajectory Reconstruction result for Jump No. 5.
Sensors 20 01995 g018
Table 1. Sensors and their key performance information.
Table 1. Sensors and their key performance information.
SensorTypeFrequencyPerformance Characteristics
GPS loggerQstarz BT Q1000eX10 HzPosition accuracy: 3 m circular error probable (50%), velocity accuracy: 0.1 m/s. [13]
IMUInvenSense ICM-20600100 HzGyroscope: measurement range: ± 2000 / s, rate noise spectral density: ± 0.004 ( / s ) / Hz ;
Accelerometer: measurement range: ±16 g, noise spectral density: 100 μ g / Hz . [14]
MagnetometerAlps HSCDTD008A100 HzMeasurement range: ±2.4 mT. [15]
Table 2. Parameters for the WGS84 reference ellipsoid.
Table 2. Parameters for the WGS84 reference ellipsoid.
NameSymbolValue
The semi-major axis a 6378137.0 m
The first eccentricity e 0.0818191908426
The meridian radius of curvature M μ a 1 e 2 ( 1 e 2 sin 2 μ 0 ) 3 2
The radius of curvature in the prime vertical N μ a 1 e 2 sin 2 μ 0
Table 3. Diagonal elements of the estimated initial states error covariance matrix P ¯ 0 .
Table 3. Diagonal elements of the estimated initial states error covariance matrix P ¯ 0 .
System StatesSymbolEstimated Variance
Position r N , 0 diag ( [ 5 5 15 ] ) 2 [ m 2 ]
Velocity v B , 0 diag ( [ 1 1 3 ] ) 2 [ ( m / s ) 2 ]
Attitude quaternions q BO , 0 ( 0.1 I 4 ) 2 [ ]
Gyroscope bias Δ ω 0 ( 0.0873 I 3 ) 2 [ ( rad / s ) 2 ]
Accelerometer bias Δ a 0 ( 0.2 I 3 ) 2 [ ( m / s 2 ) 2 ]
Magnetometer bias Δ m 0 ( 5 I 3 ) 2 [ ( μ T ) 2 ]
Magnetometer scaling error Δ S m , 0 ( 0.1 I 3 ) 2 [ ]
Table 4. Diagonal elements of the estimated process noise covariance matrix Q .
Table 4. Diagonal elements of the estimated process noise covariance matrix Q .
System InputsSymbolEstimated Variance
Gyroscope ω B , gyro ( 0.0175 I 3 ) 2 [ ( rad / s ) 2 ]
Accelerometer a B , acc ( 0.1 I 3 ) 2 [ ( m / s 2 ) 2 ]
Table 5. Diagonal elements of the estimated measurement noise covariance matrix R .
Table 5. Diagonal elements of the estimated measurement noise covariance matrix R .
System MeasurementsSymbolEstimated Variance
GPS position r N , GPS diag ( [ 3 3 5 ] ) 2 [ m 2 ]
GPS velocity v O , GPS diag ( [ 0.1 0.1 0.3 ] ) 2 [ ( m / s ) 2 ]
Magnetometer m B , mag ( 5 I 3 ) 2 [ ( μ T ) 2 ]
Quaternions constraint c q 0.01 2 [ ]
In-run vertical constraint c V , IR 0.1 2 [ m 2 ]
In-run horizontal constraint c H , IR 0.1 2 [ m 2 ]
Land area vertical constraint c V , LA 0.1 2 [ m 2 ]
Table 6. The root-mean-square error for estimated variables.
Table 6. The root-mean-square error for estimated variables.
Variables Δ x N ¯ Δ y N ¯ Δ z N ¯ Δ v O , x ¯ Δ v O , y ¯ Δ v O , z ¯ Δ ϕ ¯ Δ θ ¯ Δ ψ ¯
RMS Error:0.159 m0.123 m0.109 m0.021 m/s0.026 m/s0.023 m/s0.584 0.319 0.759
Table 7. Comparison between jump length obtained from the trajectory reconstruction results (proposed method) and video recordings (reference).
Table 7. Comparison between jump length obtained from the trajectory reconstruction results (proposed method) and video recordings (reference).
Jump Length:From Trajectory Reconstruction L TR From Video Recording L VR Error Δ L
Jump No. 191.6 m90.0 m 1.6 m
Jump No. 285.8 m85.0 m 0.8 m
Jump No. 398.6 m97.5 m 1.1 m
Jump No. 469.5 m70.0 m-0.5 m
Jump No. 586.0 m85.5 m 0.5 m

Share and Cite

MDPI and ACS Style

Fang, X.; Grüter, B.; Piprek, P.; Bessone, V.; Petrat, J.; Holzapfel, F. Ski Jumping Trajectory Reconstruction Using Wearable Sensors via Extended Rauch-Tung-Striebel Smoother with State Constraints. Sensors 2020, 20, 1995. https://doi.org/10.3390/s20071995

AMA Style

Fang X, Grüter B, Piprek P, Bessone V, Petrat J, Holzapfel F. Ski Jumping Trajectory Reconstruction Using Wearable Sensors via Extended Rauch-Tung-Striebel Smoother with State Constraints. Sensors. 2020; 20(7):1995. https://doi.org/10.3390/s20071995

Chicago/Turabian Style

Fang, Xiang, Benedikt Grüter, Patrick Piprek, Veronica Bessone, Johannes Petrat, and Florian Holzapfel. 2020. "Ski Jumping Trajectory Reconstruction Using Wearable Sensors via Extended Rauch-Tung-Striebel Smoother with State Constraints" Sensors 20, no. 7: 1995. https://doi.org/10.3390/s20071995

APA Style

Fang, X., Grüter, B., Piprek, P., Bessone, V., Petrat, J., & Holzapfel, F. (2020). Ski Jumping Trajectory Reconstruction Using Wearable Sensors via Extended Rauch-Tung-Striebel Smoother with State Constraints. Sensors, 20(7), 1995. https://doi.org/10.3390/s20071995

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop