# Validation of 3D Knee Kinematics during Gait on Treadmill with an Instrumented Knee Brace

Nicolas Reneaud, Raphaël Zory, Olivier Guérin, Luc Thomas, Serge S Colson,

Pauline Gerus, Frédéric Chorin

► To cite this version:

Nicolas Reneaud, Raphaël Zory, Olivier Guérin, Luc Thomas, Serge S Colson, et al. Validation of 3D Knee Kinematics during Gait on Treadmill with an Instrumented Knee Brace. Sensors, 2023, 23, 10.3390/s23041812. [https://doi.org/10.3390/s23041812](https://doi.org/10.3390/s23041812)

HAL Id: hal-04978353

https://univ-cotedazur.hal.science/hal-04978353v1

Submitted on 5 Mar 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Validation of 3D Knee Kinematics during Gait on Treadmill with an Instrumented Knee Brace

Nicolas Reneaud<sup>1,2,3</sup> and Raphaël Zory<sup>1,4</sup>, Olivier Guérin<sup>5</sup>, Luc Thomas<sup>2</sup>, Serge S. Colson<sup>1,*</sup>, Pauline Gerus<sup>1</sup> and Frédéric Chorin<sup>1,3</sup>

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Reneaud, N.; Zory, R.; Guérin, O.; Thomas, L.; Colson, S.S.; Gerus, P.; Chorin, F. Validation of 3D Knee Kinematics during Gait on Treadmill with an Instrumented Knee Brace. *Sensors* **2023**, *23*, 1812. https://doi.org/10.3390/s23041812

Academic Editors: Miguel Correia and Leandro José Rodrigues Machado

Received: 31 October 2022

Revised: 27 January 2023

Accepted: 31 January 2023

Published: 6 February 2023

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** To test a novel instrumented knee brace intended for use as a rehabilitation system, based on inertial measurement units (IMU) to monitor home-based exercises, the device was compared to the gold standard of motion analysis. The purpose was to validate a new calibration method through functional tasks and assessed the value of adding magnetometers for motion analysis. Thirteen healthy young adults performed a 60-second gait test at a comfortable walking speed on a treadmill. Knee kinematics were captured simultaneously, using the instrumented knee brace and an optoelectronic camera system (OCS). The intraclass correlation coefficient (ICC) showed excellent reliability for the three axes of rotation with and without magnetometers, with values ranging between 0.900 and 0.972. Pearson's r coefficient showed good to excellent correlation for the three axes, with the root mean square error (RMSE) under 3° with the IMUs and slightly higher with the magnetometers. The instrumented knee brace obtained certain clinical parameters, as did the OCS. The instrumented knee brace seems to be a valid tool to assess ambulatory knee kinematics, with an RMSE of <3°, which is sufficient for clinical interpretations. Indeed, this portable system can obtain certain clinical parameters just as well as the gold standard of motion analysis. However, the addition of magnetometers showed no significant advantage in terms of enhancing accuracy.

**Keywords:** connected knee brace; inertial measurement unit; joint kinematics; gait; validity

## 1. Introduction

Knee joint function can be altered by several pathologies. To assess the severity of the effect of these pathologies on locomotion, gait analysis is the gold standard evaluation method for clinicians [1,2]. This type of analysis involves the quantification of the movement of the lower limbs, which characterizes human locomotion [2]. While motion capture with an optoelectronic camera system (OCS) can accurately quantify knee kinematics, such measurements are costly, complex for clinicians, and not ecologically sustainable [3,4]. Recent improvements in this field led to the use of markerless motion capture systems, making it easier to implement. However, this new type of motion capture system is still expensive [5]. Human motion can also be measured with a more affordable system. Placing inertial measurement units (IMUs) on different body segments can accurately quantify three-dimensional kinematics in natural movement situations [6,7]. The measurement system contains a tri-axial accelerometer and a tri-axial gyroscope. The accuracy of this sensor can be improved with the use of a tri-axial magnetometer to limit the drift issue in terms of orientation estimation. This hardware solution entails the use of a fusion algorithm to limit the magnetic field disturbance caused by the proximity of ferromagnetic objects [8].

The three-dimensional kinematics are obtained using the respective orientations in space of each of the magnetic-inertial measurement units (MIMU). The orientation estima-

<footnote n="1">¹ Université Côte d'Azur, LAMHESS, 06205 Nice, France</footnote>

<footnote n="2">² Ted Orthopedics, 37 Rue Guibal, 13003 Marseille, France</footnote>

<footnote n="3">³ Université Côte d'Azur, CHU, 06000 Nice, France</footnote>

<footnote n="4">⁴ Institut Universitaire de France, 75231 Paris, France</footnote>

<footnote n="5">⁵ Université Côte d'Azur, CNRS, INSERM, IRCAN, 06107 Nice, France</footnote>

<footnote>* Correspondence: serge.colson@univ-cotedazur.fr; Tel.: +33-4-89-15-39-05</footnote>
tion is defined using the Cardan angle with the pitch, roll, and yaw, or with quaternions. Firstly, the gyroscope is used to estimate the orientation of the sensor by integrating the angular velocity. However, this integration method suffers from drift effects. Therefore, the accelerometer is used to estimate the roll and correct the drift on this axis. Secondly, the magnetometer is used to correct the yaw axis by sensing the geomagnetic field. Both the accelerometer and the magnetometer are easily corrupted by external acceleration and magnetic disturbances. To deal with these limitations, sensor fusion has been introduced as a widely accepted method for accurate orientation estimation. Several studies have proposed different methods using a quaternion-based Kalman filter or a complementary filter [9,10]. The low computational cost of a complementary filter with equivalent performance to the Kalman filter methods has led to the popularization of the use of these sensor fusion algorithms for embedded systems [11].

The growing interest in telemedicine and patient monitoring has led to a rapid rise in the use of wearable technologies in healthcare systems [12]. However, a practical approach to monitoring physical rehabilitation or evaluating the physical capacity of a patient at home can be challenging because it is difficult to capture meaningful data in a low-cost and easy-to-use way [13]. The monitoring system must be minimally burdensome for the patient, in order to obtain natural movement data and to ensure high patient acceptance levels regarding the technology. Therefore, quick, easy-to-use calibration is increasingly becoming a benefit for remote patient monitoring [14,15].

In this vein, a novel instrumented knee brace was developed, with two MIMUs embedded in a compressive sleeve. These sensors measure the movement of the thigh relative to the shank. The integrated system enabled a quantified analysis of knee motion via a smartphone application, without the need for clinicians, due to its ease of use. Patients were able to put on the instrumented knee brace, carry out calibration movements using the smartphone application with complete autonomy, and thus undertake home-based rehabilitation through explanatory videos and biofeedback in real time. Data measured during the rehabilitation process allowed for a quantified rehabilitation program with clinicians and physiotherapists via success scores and pain-scale charts [16].

However, few studies have focused on the use of wearable sensors from the perspective of self-management [17,18]. Before such technology can be used routinely to measure 3D knee kinematics for gait evaluation at home, the device's reliability and validity need to be reviewed to compare its performance against the gold standard. In this context, a combination of functional tasks to estimate the axes of the anatomical knee-joint angles is an important aspect [19]. The majority of wearable sensors use a combination of functional tasks or postures, with different levels of accuracy [1]. Quick and easy-to-use calibration tends to be beneficial for remote patient monitoring. Moreover, the improvement in accuracy with better calibration allows for the extraction of the clinical parameters used for monitoring rehabilitation and may be relevant for clinicians [1].

Therefore, the purpose of this study was to evaluate 3D knee kinematics from an instrumented knee brace, calibrated with a new combination of functional tasks, and then compare the performance to an OCS that is calibrated anatomically. The first hypothesis was that the instrumented knee brace could measure 3D kinematics during a gait test with an acceptable level of error (under 5 degrees) in three axes. The second objective was to assess the value of adding magnetometer data in gait-test evaluations. The second hypothesis was that an instrumented knee brace with IMUs and magnetometers would be more accurate.

## 2. Materials and Methods

### 2.1. *Subjects*

Thirteen healthy young adults (six females, seven males, 25.4 ± 2.2 years, height: 173.4 ± 7.0 cm, body mass: 71.5 ± 17.8 kg) volunteered to participate in the study after providing informed written consent and GDPR consent. They were free of any known
disorders that would affect gait and functional mobility. The study was approved by the
local ethics committee and was conducted in accordance with the Declaration of Helsinki.

## 2.2. *Materials and Instrumentation*

Knee kinematics were captured simultaneously using the Ted K-Ortho instrumented knee brace (Ted Orthopedics, Marseille, France) on the dominant lower limb in all participants, working with a network of 9 Optitrack Prime 13 optoelectronic cameras (Natural-Point, Inc., Corvallis, OR, USA). The dominant leg was defined by establishing the preferred leg for dynamic tasks (kicking a ball) [20]. Knee kinematics with the OCS were obtained via two clusters of three reflective markers, held against the lateral side of the thigh and the anterior side of the shank with stretch bands. For knee calibration, ten anatomical landmarks were identified, using reflective markers fixed onto the skin (Figure 1). Knee kinematics with the instrumented knee brace were obtained with the orientation of two ICM-20948 IMUs positioned in a soft compression sleeve, with no mechanical properties except for fixing the sensors to the thigh and shank. The IMU offers 16-bit resolution for accelerometers, gyroscopes, and magnetometers. However, due to Bluetooth low-energy communication requirements, these resolutions were downgraded to 12 bits, with a sensitivity of 1.95 mg/LSB (milli acceleration/least significant bit) for the ±4 g range, a sensitivity of 0.24 dps/LSB (degrees per second/least significant bit) for a range of ±500 dps, + and a sensitivity of 2.39 µT/LSB (micro Tesla/least significant bit) for a range of ±4900 µT, respectively. Kinematic data obtained with the instrumented knee brace and the OCS had a sampling frequency of 100 Hz.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** The instrumented knee brace is highlighted in blue, the anatomical reflective markers used for bone calibration are in green, and clusters with reflective markers used for knee kinematics are in red. Anatomical reflective markers were fixed onto the medial malleolus (MM), the lateral malleolus (LM), the medial and lateral epicondyle of the tibia (MT, LT), the medial and lateral epicondyle of the femur (MF, LF), the greater trochanter, the right and left anterior superior iliac spine, and the right and left posterior superior iliac spine.</figcaption></figure>

## 2.3. *Protocol*

### 2.3.1. Gait Test

Before the gait test, all participants completed a familiarization period on the treadmill at their comfortable speed. Comfortable gait speeds were measured with a 10 m gait test on an OptoGAIT portable optometric system (Microgate, Bolzano, Italy) in barefoot conditions. The motion capture systems were calibrated after the familiarization period.
Each participant performed one 60-second gait test on the treadmill at their own comfortable speed, in barefoot conditions (1.39 ± 0.12 m·s⁻¹).

### 2.3.2. Calibration of Clusters

To evaluate the position of the clusters relative to the bones, the participants completed both a static and a functional calibration procedure. Static anatomical calibration was performed according to standard laboratory protocol. Participants performed a movement with the lower limb fully extended, making a circular hip movement. The hip joint can be modeled as a ball-and-socket joint, and the hip joint center (HJC) location can be characterized by a point that is invariant in any position of the joint [21]. The coordinates of the femoral head location were obtained by an optimization method to minimize the HJC movement during the circular movements, compared with the hip reference [22]. The distal femur location was determined at the middle of the MF and LF markers. The proximal and distal tibia locations were calculated between the MT and LT markers and the MM and LM markers, respectively.

The mechanical axis of the femur was defined by the line containing the femoral head location and the distal femur location. The femoral condylar axis was defined by the MF and LF markers. The femoral frontal plane was defined as the plane containing the mechanical axis and the femoral condylar axis. The femoral sagittal plane was defined using the femoral axis and the cross-product of the femoral axis and the femoral condylar axis [21]. The mechanical axis of the tibia was defined by the line containing the proximal and distal locations of the tibia. The tibial condylar axis was defined by the MT and LT markers. The tibial frontal plane was defined by the plane containing the mechanical axis and the tibial condylar axis. The tibial sagittal plane was defined using the tibial axis and the cross-product of the tibial axis and the tibial condylar axis [21].

### 2.3.3. Instrumented Knee Brace Calibration

Prior to calibration, the participant must put the brace with the thigh sensor on the lateral face of the thigh and the shank sensor on the anterior-medial surface of the tibia [17]. Two marks indicate the placement of the sensors in the instrumented knee brace. To evaluate the orientation of the IMUs relative to the bones, the participants conducted a static and also a functional calibration procedure. The system was calibrated by having the participant remain stationary in a neutral reference posture, with the lower limb fully extended during calibration. The mechanical axes of the femur and the tibia were defined by the gravity vector. Then, participants performed a functional movement comprising small knee flexions, considering the joint as a hinge joint, when extended, to define the extension/flexion axis [21]. The sagittal plane of the knee was defined by the movement of the IMU of the thigh, compared to the IMU of the shank, during the small flexions. To synchronize the two systems, a cross-correlation between flexion and extension was used to identify the time difference between the OCS and IMU systems.

### 2.3.4. Sensor Fusion Algorithm

The sensor fusion algorithm used by the instrumented knee brace is a two-step complementary filter [11]. This method decouples the pitch and roll estimation (attitude) from the yaw estimation (heading). The accelerometer and the gyroscope are used to estimate the attitude in quaternion form. This quaternion is used as a correction for the roll and pitch components. Then, the magnetometer data are used for the second step of the complementary filter, to correct the heading of the previous estimation by performing a rotation about the global z-axis, in order to align the current frame with the magnetic field. The remainder of this section closely follows the work of Valenti et al. [11].

The orientation in the 3D space of the sensor frame with respect to the earth frame can be represented by a unit quaternion, $S_E^q$, defined as follows:

$$S_E^q = [q_0 \ q_1 \ q_2 \ q_3]^T = [\cos \frac{\alpha}{2} \ e_x \sin \frac{\alpha}{2} \ e_y \sin \frac{\alpha}{2} \ e_z \sin \frac{\alpha}{2}]^T \quad (1)$$
where $\alpha$ is the rotation angle and $e$ is the unit vector that represents the rotation axis.

Unit quaternions can be applied to operate the rotations of 3D vectors. For example, the vector $^S v_q$, expressed with respect to the sensor frame, can also be expressed with respect to the earth frame using the following operation:

$$ ^E v_q = ^E S q \otimes ^S v_q \otimes ^E S q^* = ^E S q^* \otimes ^S v_q \otimes ^E S q \quad (2) $$

where the vector $v$ is expressed as a quaternion:

$$ v_q = [0 \quad v]^{\mathrm{T}} = [0 \quad v_x \quad v_y \quad v_z]^{\mathrm{T}}. \quad (3) $$

We defined the quaternion obtained with the magnetometer $q_{mag}$ to have only a single degree of freedom, to correct the heading around the z-axis only, by setting it to:

$$ q_{mag} = [q_{0mag} \quad 0 \quad 0 \quad q_{3mag}]^{\mathrm{T}}. \quad (4) $$

The rotation defined in Equation (2) can be expressed in a matrix, as in Equation (5):

$$ ^E v = R(^E S q)^S v \quad (5) $$

where $R(^E S q)$ is the direct cosine matrix, given in terms of the orientation quaternion $^E S q$, as shown below:

$$ R(^E S q) = \begin{bmatrix} q_0^2 + q_1^2 - q_2^2 - q_3^2 & 2(q_1q_2 - q_0q_3) & 2(q_1q_3 + q_0q_2) \\ 2(q_1q_2 + q_0q_3) & q_0^2 - q_1^2 + q_2^2 - q_3^2 & 2(q_2q_3 - q_0q_1) \\ 2(q_1q_3 - q_0q_2) & 2(q_2q_3 + q_0q_1) & q_0^2 - q_1^2 - q_2^2 + q_3^2 \end{bmatrix}. \quad (6) $$

This two-step complementary filter is based on the work of Valenti et al. [11]. The block diagram of the method is shown in Figure 2. Before the first step for the complementary filter, the gyroscope bias is removed with the calibration file that was obtained through the static phases during the instrumented knee-brace calibration procedure. Then, the complementary filter is initialized with a first estimation of the sensor orientation in quaternion form, with the accelerometer and the magnetometer data. Assuming this initial condition, the filter calculates the quaternion derivative, $^E \dot{q}_{\omega,t}$, describing the rate of change of the orientation, $\omega$, obtained through the gyroscope and the previous frame of the sensor [11].

$$ ^E \dot{q}_{\omega,t} = -\frac{1}{2} S \omega_{q,t} \otimes ^E S q_{t-1} \quad (7) $$

Then, the orientation of the sensor, $^E S q_{\omega,t}$, can be obtained by integrating the quaternion derivative. This first orientation estimation is used in the first step of the complementary filter to predict the gravity vector, $^E g_p$, and calculate the deviation from the real gravity vector, $^E g$, defined as $\Delta q_{acc}$. The delta quaternion, $\Delta q_{acc}$, is then filtered to reduce the noise with a linear interpolation (LERP) or a spherical linear interpolation (SLERP), depending on the deviation of the predicted gravity vector from the real gravity vector, according to the work of Valenti et al. [11]. The LERP is used in the case of small deviation, when $|| \sin(^E g_p, ^E g) || < 0.01$, using the equation below:

$$ \overline{\Delta q}_{acc} = (1 - K_a)q_I + K_a \Delta q_{acc} \quad (8) $$

where $q_I$ is the identity quaternion and $K_a \in [0, 1]$ is the gain that characterizes the level of confidence in the accelerometer, compared to the gyroscope.
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Block diagram of the implemented quaternion-based two-step sensor fusion algorithm, based on the work of Valenti et al. [11]. A canonical vector replaces the magnetometer for the IMU alone. $K_a = 0.2$ and $K_m = 0.1$; $\omega$: measurement from the gyroscope; $a$: measurement from the accelerometer; $m$: measurement from the magnetometer; $\hat{\Delta q}_{acc}$: orientation estimated from the accelerometer; $\hat{\Delta q}_{mag}$: the correction quaternion, calculated from the magnetometer; $\hat{S}_{E}q$: the orientation of the sensor frame, relative to the earth frame; $\otimes$: quaternion multiplication operator.</figcaption></figure>

Then, the delta quaternion $\overline{\Delta q}_{acc}$ is normalized:

$$ \hat{\Delta q}_{acc} = \frac{\overline{\Delta q}_{acc}}{\|\overline{\Delta q}_{acc}\|} \qquad (9) $$

If $\|\sin(Eg_p, Eg)\| \ge 0.01$, the SLERP is used to filter the $\Delta q_{acc}$ with the following equation:

$$ \hat{\Delta q}_{acc} = \frac{\sin((1 - K_a) * (Eg_p, Eg))}{\sin(Eg_p, Eg)} q_l + \frac{\sin(K_a * (Eg_p, Eg))}{\sin(Eg_p, Eg)} \Delta q_{acc}. \qquad (10) $$

Finally, the quaternion estimated from the gyroscope data is multiplied by the filtered delta quaternion, to correct the orientation estimation in the attitude component (Figure 3b).

$$ \hat{S}_{E}q'_{t} = \hat{S}_{E}q_{\omega,t} \otimes \hat{\Delta q}_{acc} \qquad (11) $$

A second step is performed to apply a correction on the heading component. In the same way as the previous step, the magnetometer data is used to compute a delta quaternion between the predicted geomagnetic field, $E l_p$, and the real one, $E l$. This delta quaternion, $\Delta q_{mag}$, performs a rotation that is only about the global z-axis, by aligning the global x-axis into the positive direction of magnetic North [11] (Figure 3c). If the magnetic field data is not provided, a canonical vector on the x-axis (Figure 2) replaces the magnetometer vector. The filtered delta quaternion of the magnetometer, $\hat{\Delta q}_{mag}$, is calculated with the same method of LERP or SLERP presented previously, with a constant gain, $K_m \in [0, 1]$ (Figure 2). The second correction of the orientation estimation is obtained by the multiplication of the filtered delta quaternion, $\hat{\Delta q}_{mag}$, with the orientation estimation in quaternion form, obtained via step 1 of the complementary filter.

$$ \hat{S}_{E}q_{t} = \hat{S}_{E}q'_{t} \otimes \hat{\Delta q}_{mag} = \hat{S}_{E}q_{\omega,t} \otimes \hat{\Delta q}_{acc} \otimes \hat{\Delta q}_{mag} \qquad (12) $$
<figure><img src="image_6.png" /><figcaption>**Figure 3.** The schematic of orientation estimation-solving with the two-step complementary filter: (a) is the initial orientation of the sensor, with the sensor frame $R_{MIMU} = (\vec{X}', \vec{Y}', \vec{Z}')$ in the Earth frame $R_{earth} = (\vec{X}, \vec{Y}, \vec{Z})$; (b) shows the first correction on the attitude, with the delta quaternion, $\Delta q_{acc}$; (c) shows the second correction on the heading, with the delta quaternion, $\Delta q_{mag}$.</figcaption></figure>

### 2.3.5. Knee Kinematics

The kinematics obtained with the instrumented knee brace and the OCS were calculated by the quaternion, $^{Th}_{Sh}q$, which defines the orientation of the thigh frame, $^{Th}_{E}q$, relative to the shank frame, $^{Sh}_{E}q$.

$$^{Th}_{Sh}q = ^{Th}_{E}q \otimes ^{Sh}_{E}q^* = ^{Th}_{E}q \otimes ^{E}_{Sh}q \quad (13)$$

Then, to obtain the flexion/extension, internal/external rotation, and abduction/adduction, the quaternion, $^{Th}_{Sh}q$, was converted to Cardan angles (pitch, roll, and yaw) in the local frame, with the following equations:

$$flex = \tan^{-1} \left( \frac{2(q_0 q_1 + q_2 q_3)}{q_0^2 - q_1^2 - q_2^2 + q_3^2} \right) \quad (14)$$

$$abd = \sin^{-1}(2(q_0 q_2 + q_1 q_3)) \quad (15)$$

$$rot = \tan^{-1} \left( \frac{2(q_0 q_3 + q_1 q_2)}{q_0^2 + q_1^2 - q_2^2 - q_3^2} \right). \quad (16)$$

### 2.4. Data Analysis

The intra-reliability of the instrumented knee brace was determined by the intraclass correlation coefficient (ICC) type (2,1) between participants and gait cycles, based on the range of motion (RoM) from the first 40 strides [23]. ICC was classified as poor (<0.50), moderate (0.50–0.75), good (between 0.75 and 0.90), or excellent (>0.90) [3]. The validity between the instrumented knee brace and the OCS was calculated, based on the complete acquisition of the three knee angles, using the root mean square error (RMSE) computed for the full curve, without offset difference and Pearson's correlation coefficient. An error of RMSE of <5° was considered excellent, and between 5° and 10° was considered good [6]. Pearson's coefficient was evaluated as weak (<0.65), moderate (0.65–0.75), good (0.75–0.85), very good (0.85–0.95), or excellent (>0.95) [3].

From the study by Robert-Lachaine et al., twelve clinical parameters were extracted from the three knee rotations, based on normalized gait cycles, to characterize gait (Table 1) [1]. The assumption of normality-distributed data was checked using the Shapiro-Wilk test.
The differences between the IMU system and the OCS, and between the IMUs and MIMUs, were calculated by paired *t*-tests.

<table><caption>**Table 1.** Definitions of the calculation of the clinical parameters.</caption><thead><tr><th>Clinical Parameter</th><th>Definition</th></tr></thead><tbody><tr><td>Initial flexion</td><td>Flexion at foot strike</td></tr><tr><td>Flexion during loading</td><td>Difference between the maximal flexion, from 0% to 20% of the gait cycle, and the initial flexion</td></tr><tr><td>Flexion during stance</td><td>Difference between the maximal flexion, from 0% to 20% of the gait cycle, and the minimal flexion from 20% to 68% of the gait cycle</td></tr><tr><td>Maximal flexion</td><td>Maximal flexion of the complete gait cycle</td></tr><tr><td>Flexion RoM</td><td>Flexion range of motion over the complete gait cycle</td></tr><tr><td>Initial adduction</td><td>Adduction at foot strike</td></tr><tr><td>Varus thrust</td><td>Difference between the maximal adduction, from 0% to 20% of the gait cycle, and the initial adduction</td></tr><tr><td>Valgus thrust</td><td>Difference between the initial adduction and the minimum adduction from 0% to 20% of the gait cycle</td></tr><tr><td>Adduction during stance</td><td>Mean adduction from 20% to 54% of the gait cycle</td></tr><tr><td>Adduction RoM</td><td>Adduction range of motion over the complete gait cycle</td></tr><tr><td>Initial tibial rotation</td><td>Tibial rotation at foot strike</td></tr><tr><td>Tibial rotation during loading</td><td>Mean tibial rotation from 0% to 20% of the gait cycle</td></tr><tr><td>Tibial rotation RoM</td><td>Tibial rotation range of motion over the complete gait cycle</td></tr></tbody></table>

## 3. Results

### 3.1. Reproducibility of the Kinematic Variables Recorded with the Instrumented Knee Brace

The ICC (95%) for the kinematic range of motion recorded with the instrumented knee brace showed excellent reliability for the three axes of rotation, both with and without magnetometers, with values ranging between 0.900 and 0.972 (Table 2).

<table><caption>**Table 2.** The intraclass reliability and validation of the IMU system. ICC values were assessed on the RoM of the gait cycle. The difference between the IMUs and MIMUs was highlighted by paired *t*-tests (* *p*-value < 0.05).</caption><thead><tr><th>Data Analysis</th><th>Flex/Ext (CI 95%)</th><th>Rot Int/Ext (CI 95%)</th><th>Abd/Add (CI 95%)</th></tr></thead><tbody><tr><td>ICC</td><td>0.922 (0.850–0.974)</td><td>0.923 (0.852–0.974)</td><td>0.941 (0.885–0.980)</td></tr><tr><td>ICC mag</td><td>0.900 (0.821–0.961)</td><td>0.901 (0.822–0.962)</td><td>0.972 (0.946–0.990)</td></tr><tr><td>Pearson's coefficient</td><td>0.997 (0.997–0.998)</td><td>0.755 (0.759–0.888)</td><td>0.895 (0.653–0.985)</td></tr><tr><td>Pearson's coefficient mag.</td><td>0.997 (0.996–0.998)</td><td>0.789 (0.747–0.893)</td><td>0.873 (0.614–0.972)</td></tr><tr><td>RMSE (°)</td><td>1.79 (1.61–1.96)</td><td>2.16 (1.82–2.50)</td><td>2.38 (1.59–3.17)</td></tr><tr><td>RMSE mag (°)</td><td>1.84 (1.61–2.07)</td><td>3.58 (2.00–5.15) *</td><td>2.76 (1.64–3.88) *</td></tr></tbody></table>

### 3.2. Concurrent Validity of the Instrumented Knee Brace during a Gait Test

In terms of concurrent validity, Pearson's *r* coefficient between systems showed an excellent correlation for the F/E axis (σ > 0.95), a good correlation for the I/E rotation axis (σ > 0.75), and a very good correlation for the A/A axis (σ > 0.85) (Table 2). The RMSE was less than 3° between the instrumented knee brace and the OCS. The addition of the magnetometer data to compute the knee kinematics showed a larger measurement error, especially for the I/E rotation axis, with an RMSE of 3.58°. Paired *t*-tests performed between the IMUs and MIMUs modalities showed a significant difference in the RMSE measured between the instrumented knee brace and the OCS for the tibial rotation and abduction axes (*p*-value < 0.05) (Table 2). The 3D kinematics of the knee are depicted in Figure 4, which compares the measurements with the IMUs, MIMUs, and OCS.
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Comparison of mean knee kinematics, normalized to the gait cycle, between the OCS (orange), IMUs (blue), and MIMUs (green), involving flexion, external rotation, and abduction for spontaneous gait speed.</figcaption></figure>

### 3.3. Concurrent Validity of the Gait Parameters Obtained with the Instrumented Knee Brace

Table 3 shows the difference in clinical parameters between the two systems and the two modalities (IMU and MIMU vs. OCS). The instrumented knee brace obtained a similar RoM, with a significantly positive 6° offset for flexion/extension for both modalities (*p*-value = 0.002). The RoM of adduction, the varus thrust, and the tibial rotation were overestimated. Other parameters showed no significant differences between the instrumented knee brace and the OCS (Table 3).

<table><caption>**Table 3.** Comparison of the clinical parameters. Paired *t*-tests were used and significant *p*-values are identified in bold (*p*-value < 0.05).</caption><thead><tr><th>Clinical Parameter</th><th>OCS</th><th>IMU</th><th>Δ</th><th>p</th><th>MIMU</th><th>Δ</th><th>p</th></tr></thead><tbody><tr><td>Initial flexion (°)</td><td>-0.4 (7.1)</td><td>5.6 (5.9)</td><td>-6.0 (1.5)</td><td><strong>0.002</strong></td><td>5.7 (5.9)</td><td>-6.0 (1.6)</td><td><strong>0.002</strong></td></tr><tr><td>Flexion during loading (°)</td><td>18.0 (6.0)</td><td>21.3 (6.6)</td><td>-3.3 (0.4)</td><td><strong>&lt;0.001</strong></td><td>21.3 (6.6)</td><td>-3.3 (0.4)</td><td><strong>&lt;0.001</strong></td></tr><tr><td>Flexion during stance (°)</td><td>14.0 (4.0)</td><td>15.8 (4.1)</td><td>-1.8 (0.3)</td><td><strong>&lt;0.001</strong></td><td>15.7 (4.1)</td><td>-1.7 (0.3)</td><td><strong>&lt;0.001</strong></td></tr><tr><td>Maximal flexion (°)</td><td>67.7 (7.1)</td><td>73.7 (5.8)</td><td>-5.9 (1.4)</td><td><strong>0.001</strong></td><td>73.5 (6.0)</td><td>-5.8 (1.4)</td><td><strong>0.001</strong></td></tr><tr><td>Flexion RoM (°)</td><td>69.1 (4.9)</td><td>68.8 (5.6)</td><td>0.2 (0.6)</td><td>0.676</td><td>68.7 (5.9)</td><td>0.4 (0.6)</td><td>0.562</td></tr><tr><td>Initial adduction (°)</td><td>-4.5 (4.9)</td><td>-3.6 (3.9)</td><td>-0.8 (1.9)</td><td>0.676</td><td>-3.7 (4.6)</td><td>-0.7 (2.0)</td><td>0.715</td></tr><tr><td>Varus thrust (°)</td><td>1.7 (1.2)</td><td>2.5 (1.1)</td><td>-0.9 (0.3)</td><td><strong>0.016</strong></td><td>2.8 (1.1)</td><td>-1.1 (0.3)</td><td><strong>0.017</strong></td></tr><tr><td>Valgus thrust (°)</td><td>0.5 (0.6)</td><td>0.5 (0.4)</td><td>0.0 (0.2)</td><td>0.961</td><td>0.7 (1.4)</td><td>-0.3 (0.4)</td><td>0.562</td></tr><tr><td>Adduction during stance (°)</td><td>-2.6 (4.8)</td><td>-0.7 (3.4)</td><td>-1.9 (1.5)</td><td>0.217</td><td>-0.8 (3.6)</td><td>-1.8 (1.5)</td><td>0.259</td></tr><tr><td>Adduction RoM (°)</td><td>10.1 (4.0)</td><td>14.3 (5.1)</td><td>-4.2 (1.4)</td><td><strong>0.011</strong></td><td>15.0 (6.9)</td><td>-4.9 (1.8)</td><td><strong>0.020</strong></td></tr><tr><td>Initial tibial rotation (°)</td><td>0.0 (4.6)</td><td>-3.0 (4.2)</td><td>3.1 (1.6)</td><td>0.079</td><td>-2.9 (5.6)</td><td>2.9 (1.7)</td><td>0.106</td></tr><tr><td>Tibial rotation during loading (°)</td><td>-1.6 (4.3)</td><td>-3.1 (2.0)</td><td>1.4 (1.5)</td><td>0.340</td><td>-3.0 (3.5)</td><td>1.3 (1.5)</td><td>0.392</td></tr><tr><td>Tibial rotation RoM (°)</td><td>11.6 (3.0)</td><td>16.2 (5.1)</td><td>-4.6 (1.1)</td><td><strong>&lt;0.001</strong></td><td>15.8 (4.8)</td><td>-4.3 (1.0)</td><td><strong>&lt;0.001</strong></td></tr></tbody></table>
## 4. Discussion

The purpose of this study was to evaluate the accuracy and reliability of an instrumented knee brace during a gait test, using a new calibration process and with the addition of magnetometers. The results showed excellent reliability and a strong correlation with the OCS. Kinematic data showed a positive offset of 6° for the flexion/extension computed by the instrumented knee brace with the new calibration. However, the first hypothesis was validated, with an acceptable level of error below 5° on the three axes. The second hypothesis was refuted because the study showed a lower accuracy on the abduction and rotation axes with the use of magnetometers, due to the ferromagnetic disturbance caused by the treadmill motors and the metallic structure supporting the harness.

Regarding the validation of the instrumented knee brace, several studies have been conducted to reliably detect knee kinematics with a wearable system. The reliability reported in our study is in line with previous studies, which found fair to excellent intrarater reliability through a systematic review and, more specifically, ICCs between 0.609 and 0.989, depending on the knee axes [4,6]. This level of reliability is associated with a strong correlation between the gold standard and the instrumented knee brace. The major axis (flexion/extension) showed an excellent correlation, with a lower Pearson's correlation coefficient, ranging from good to very good, for the minor axes. These results can be explained by a crosstalk effect, which projects the rotations of the sagittal axis onto the frontal and transverse axes, due to small errors in bone orientation estimation during the calibration phases [1,24]. The correlations between systems in the sagittal plane were also similar to those reported in the literature [25,26]. These correlation results were in line with the RMSE values that were obtained with the modalities of the IMUs and the MIMUs, compared to the OCS. The major axis in this study (flexion/extension) showed an RMSE of <2°, while the two minor axes showed a slightly higher RMSE compared with the 5° [17,26] and 3° [4] results reported in previously published studies. Therefore, both systems validate our first hypothesis; the new calibration process provides an acceptable error level under 5°.

Despite this low level of error, the comparison of clinical parameters showed some limitations to the calibration procedure. The results highlighted an offset of 6° on the sagittal axis and led to significant differences for all clinical parameters except the RoM, which is in line with the same overestimation of the initial flexion reported in the study by Kayaalp et al. [17]. Conversely, our study showed a significantly higher RoM in the frontal and transverse planes and no difference for the other parameters, as reported in the study by Robert-Lachaine et al. However, the results of that study are limited because the absolute differences between systems were not reported [1]. These differences in measurements that are highlighted in the current study imply that care should be taken when comparing the gait parameters obtained with the instrumented knee brace to those obtained with another motion capture system.

Based on the RMSE values and *t*-tests between the IMUs and MIMUs, the second hypothesis was rejected, as there was no significant improvement in accuracy in the different axes when incorporating the magnetometers. The results revealed a higher RMSE for abduction and rotation with the addition of magnetometers in the calculation of the kinematics. These findings were in line with those of Robert-Lachaine et al., who suggested that magnetic disturbances might affect the RMSE, specifically the axial axis [8]. The studies of Valenti et al. and Fan et al. explained this issue with a stress test of ferromagnetic disturbances, achieved by placing a magnet near a MIMU sensor during different periods. The results show the immunity of the attitude estimation but a large error in the heading estimation when using the two-step complementary filter [9,11]. This hypothesis was suggested because the steel structure supporting the harness and the motors of the treadmill for the gait test may have induced an inhomogeneous magnetic field during data acquisition [27].

A potential improvement of the instrumented knee brace would have been to use an adaptive gain for the magnetometer. The filter presented by Fan et al. showed an
adaptative gain, $K_m$, for the filter, depending on the magnetic disturbance level. The results demonstrated good performance by using a finite state machine, which switches between different gains depending on the magnetic disturbance level. The measurement error of heading estimation was reduced from more than nineteen degrees to less than one degree [9]. The implementation of this adaptive gain filter for the magnetometer would reduce the error on the heading and may outperform the IMU accuracy of the instrumented knee brace on the axial axis.

A few limitations of this study should be acknowledged when interpreting the results. The OCS was based on clusters and simple kinematics without constraints to obtain the knee angles. This method for calculating kinematics is impacted by soft tissue and skin artifacts (STA) that could introduce measurement errors in the femur and tibia orientations [28]. However, this method was used to compare similar kinematic data and validate the instrumented knee-brace calibration. Therefore, kinematics from the instrumented knee brace showed good axes estimations, but these could be subject to STA, as in our previous study [29].

To conclude, a comparison of the instrumented knee brace against OCS measurements during a gait test demonstrated excellent reliability and strong agreement. An IMU-based system seems to be a valid tool by which to assess ambulatory knee kinematics, with an RMSE of <3°, which is sufficient for clinical interpretations. Indeed, this portable system can obtain certain clinical parameters just as well as the gold standard OCS. However, the addition of magnetometers with constant gain showed no significant advantage in terms of enhancing the accuracy with a two-step complementary filter, due to ferromagnetic disturbances in the laboratory. This finding could be different in environments free of magnetic disturbances.

**Author Contributions:** Conceptualization, N.R., R.Z. and F.C.; methodology, N.R., R.Z., L.T., S.S.C., P.G. and F.C.; software, N.R. and L.T.; validation, N.R., R.Z., P.G. and F.C.; formal analysis, N.R., L.T., P.G. and F.C.; investigation, O.G., data curation, N.R. and L.T.; writing—original draft preparation, N.R.; writing—review and editing, N.R., R.Z., S.S.C., P.G. and F.C.; supervision, R.Z., P.G. and F.C.; project administration, O.G. and F.C. All authors have read and agreed to the published version of the manuscript.

**Funding:** The institutes of the authors have received funding from the company marketing the instrumented knee brace, Ted Orthopedics, to cover the cost involved in the execution of this study.

**Institutional Review Board Statement:** The study was conducted according to the guidelines of the Declaration of Helsinki and was approved by the national ethics committee (C.P.P. EST I-RCB 19.10.30.42054).

**Informed Consent Statement:** Informed consent was obtained from all participants involved in the study. Written informed consent has been obtained from the participant (s) to publish this paper.

**Data Availability Statement:** Not applicable.

**Acknowledgments:** The Authors acknowledge Ted Orthopedics for their technical and financial contribution to this study. The study sponsors had no role in the study design and interpretation of data, the writing of the manuscript, or in the decision to submit the manuscript for publication.

**Conflicts of Interest:** I, Nicolas RENEAUD, declare that Raphaël ZORY, Frédéric CHORIN, Pauline GERUS, Serge S. COLSON, and Olivier GUERIN, have no proprietary, financial, professional, or other personal interest of any nature or of any kind in any product, service, and/or company that could be construed as influencing the position presented in the manuscript.

## References

1. Robert-Lachaine, X.; Parent, G.; Fuentes, A.; Hagemeister, N.; Aissaoui, R. Inertial Motion Capture Validation of 3D Knee Kinematics at Various Gait Speed on the Treadmill with a Double-Pose Calibration. *Gait Posture* **2020**, *77*, 132–137. [CrossRef] [PubMed]

2. Yi, C.; Jiang, F.; Bhuiyan, M.Z.A.; Yang, C.; Gao, X.; Guo, H.; Ma, J.; Su, S. Smart Healthcare-Oriented Online Prediction of Lower-Limb Kinematics and Kinetics Based on Data-Driven Neural Signal Decoding. *Future Gener. Comput. Syst.* **2021**, *114*, 96–105. [CrossRef]
3. Cuesta-Vargas, A.I.; Galán-Mercant, A.; Williams, J.M. The Use of Inertial Sensors System for Human Motion Analysis. *Phys. Ther. Rev.* **2010**, *15*, 462–473. [CrossRef]

4. Cho, Y.S.; Jang, S.H.; Cho, J.S.; Kim, M.J.; Lee, H.D.; Lee, S.Y.; Moon, S.B. Evaluation of Validity and Reliability of Inertial Measurement Unit-Based Gait Analysis Systems. *Ann. Rehabil. Med.* **2018**, *42*, 872–883. [CrossRef] [PubMed]

5. Kanko, R.M.; Laende, E.K.; Davis, E.M.; Selbie, W.S.; Deluzio, K.J. Concurrent Assessment of Gait Kinematics Using Marker-Based and Markerless Motion Capture. *J. Biomech.* **2021**, *127*, 110665. [CrossRef]

6. Poitras, I.; Dupuis, F.; Bielmann, M.; Campeau-Lecours, A.; Mercier, C.; Bouyer, L.J.; Roy, J.S. Validity and Reliability Of wearable Sensors for Joint Angle Estimation: A Systematic Review. *Sensors* **2019**, *19*, 1555. [CrossRef] [PubMed]

7. Picerno, P. 25 Years of Lower Limb Joint Kinematics by Using Inertial and Magnetic Sensors: A Review of Methodological Approaches. *Gait Posture* **2017**, *51*, 239–246. [CrossRef] [PubMed]

8. Robert-Lachaine, X.; Mecheri, H.; Larue, C.; Plamondon, A. Effect of Local Magnetic Field Disturbances on Inertial Measurement Units Accuracy. *Appl. Ergon.* **2017**, *63*, 123–132. [CrossRef]

9. Fan, B.; Li, Q.; Liu, T. Improving the Accuracy of Wearable Sensor Orientation Using a Two-Step Complementary Filter with State Machine-Based Adaptive Strategy. *Meas Sci. Technol.* **2018**, *29*, 115104. [CrossRef]

10. del Rosario, M.B.; Lovell, N.H.; Redmond, S.J. Quaternion-Based Complementary Filter for Attitude Determination of a Smartphone. *IEEE Sens. J.* **2016**, *16*, 6008–6017. [CrossRef]

11. Valenti, R.G.; Dryanovski, I.; Xiao, J. Keeping a Good Attitude: A Quaternion-Based Orientation Filter for IMUs and MARGs. *Sensors* **2015**, *15*, 19302–19330. [CrossRef] [PubMed]

12. Lunardini, F.; Malavolti, M.; Pedrocchi, A.L.G.; Borghese, N.A.; Ferrante, S. A Mobile App to Transparently Distinguish Single- from Dual-Task Walking for the Ecological Monitoring of Age-Related Changes in Daily-Life Gait. *Gait Posture* **2021**, *86*, 27–32. [CrossRef] [PubMed]

13. Sarker, A.; Emenonye, D.R.; Kelliher, A.; Rikakis, T.; Buehrer, R.M.; Asbeck, A.T. Capturing Upper Body Kinematics and Localization with Low-Cost Sensors for Rehabilitation Applications. *Sensors* **2022**, *22*, 2300. [CrossRef]

14. Hamilton, R.I.; Williams, J.; Holt, C. Biomechanics beyond the Lab: Remote Technology for Osteoarthritis Patient Data—A Scoping Review. *Front. Rehabil. Sci.* **2022**, *3*, 235. [CrossRef]

15. Turner, A.J.; Carroll, W.; Kodithuwakku Arachchige, S.N.K.; Saucier, D.; Burch V, R.F.; Ball, J.E.; Smith, B.K.; Freeman, C.E.; Knight, A.C.; Chander, H. Closing the Wearable Gap—Part VIII: A Validation Study for a Smart Knee Brace to Capture Knee Joint Kinematics. *Biomechanics* **2021**, *1*, 152–162. [CrossRef]

16. Reneaud, N.; Gerus, P.; Chorin, F.; Zory, R.; Guérin, O.; Roux, C. Effects of Home Based-Progressive Rehabilitation Program on Assessment of Functional and Pain Outcomes in People with Knee Osteoarthritis: A Pilot Study. *Osteoarthr. Cartil.* **2021**, *29*, S10–S11. [CrossRef]

17. Kayaalp, M.E.; Agres, A.N.; Reichmann, J.; Bashkuev, M.; Duda, G.N.; Becker, R. Validation of a Novel Device for the Knee Monitoring of Orthopaedic Patients. *Sensors* **2019**, *19*, 5193. [CrossRef]

18. Bolam, S.M.; Batinica, B.; Yeung, T.C.; Weaver, S.; Cantamessa, A.; Vanderboor, T.C.; Yeung, S.; Munro, J.T.; Fernandez, J.W.; Besier, T.F.; et al. Remote Patient Monitoring with Wearable Sensors Following Knee Arthroplasty. *Sensors* **2021**, *21*, 5143. [CrossRef]

19. Nazarahari, M.; Noamani, A.; Ahmadian, N.; Rouhani, H. Sensor-to-Body Calibration Procedure for Clinical Motion Analysis of Lower Limb Using Magnetic and Inertial Measurement Units. *J. Biomech.* **2019**, *85*, 224–229. [CrossRef] [PubMed]

20. Promsri, A.; Haid, T.; Federolf, P. How Does Lower Limb Dominance Influence Postural Control Movements during Single Leg Stance? *Hum. Mov. Sci.* **2018**, *58*, 165–174. [CrossRef] [PubMed]

21. Jenny, J.; Boeri, C.; Picard, F.; Leitner, F. Reproducibility of Intra-Operative Measurement of the Mechanical Axes of the Lower Limb during Total Knee Replacement with a Non-Image-Based Navigation System. *Comput. Aided Surg.* **2004**, *9*, 161–165. [CrossRef]

22. Fiorentino, N.M.; Atkins, P.R.; Kutschke, M.J.; Foreman, K.B.; Anderson, A.E. In-Vivo Quantification of Dynamic Hip Joint Center Errors and Soft Tissue Artifact. *Gait Posture* **2016**, *50*, 246–251. [CrossRef]

23. König, N.; Singh, N.B.; von Beckerath, J.; Janke, L.; Taylor, W.R. Is Gait Variability Reliable? An Assessment of Spatio-Temporal Parameters of Gait Variability during Continuous Overground Walking. *Gait Posture* **2014**, *39*, 615–617. [CrossRef] [PubMed]

24. Baudet, A.; Morisset, C.; D'athis, P.; Maillefert, J.-F.; Casillas, J.-M.; Ornetti, P.; Laroche, D. Cross-Talk Correction Method for Knee Kinematics in Gait Analysis Using Principal Component Analysis (PCA): A New Proposal. *PLoS ONE* **2014**, *9*, e102098. [CrossRef] [PubMed]

25. Bolink, S.A.A.N.; Naisas, H.; Senden, R.; Essers, H.; Heyligers, I.C.; Meijer, K.; Grimm, B. Validity of an Inertial Measurement Unit to Assess Pelvic Orientation Angles during Gait, Sit-Stand Transfers and Step-up Transfers: Comparison with an Optoelectronic Motion Capture System. *Med. Eng. Phys.* **2016**, *38*, 225–231. [CrossRef]

26. Zügner, R.; Tranberg, R.; Timperley, J.; Hodgins, D.; Mohaddes, M.; Kärrholm, J. Validation of Inertial Measurement Units with Optical Tracking System in Patients Operated with Total Hip Arthroplasty. *BMC Musculoskelet. Disord.* **2019**, *20*, 52. [CrossRef] [PubMed]

27. de Vries, W.H.K.; Veeger, H.E.J.; Baten, C.T.M.; van der Helm, F.C.T. Magnetic Distortion in Motion Labs, Implications for Validating Inertial Magnetic Sensors. *Gait Posture* **2009**, *29*, 535–541. [CrossRef]
28. Barré, A.; Jolles, B.M.; Theumann, N.; Aminian, K. Soft Tissue Artifact Distribution on Lower Limbs during Treadmill Gait: In Fluctuation of Skin Markers' Location on Cluster Design. *J. Biomech.* **2015**, *48*, 1965–1971. [CrossRef]

29. Reneaud, N.; Pierre, E.; Zory, R.; Chorin, F.; Thomas, L.; Chavet, P.; Coyle, T.; Truchet, E.; Puech, S.; Ollivier, M.; et al. Validity of an Instrumented Knee Brace Compared to 3D Motion Navigation: A Cadaveric Investigation. *Measurement* **2021**, *173*, 108590. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.