Inertial Sensor-to-Segment Calibration for Accurate 3D Joint Angle Calculation for Use in OpenSim

Inertial capture (InCap) systems combined with musculoskeletal (MSK) models are an attractive option for monitoring 3D joint kinematics in an ecological context. However, the primary limiting factor is the sensor-to-segment calibration, which is crucial to estimate the body segment orientations. Walking, running, and stair ascent and descent trials were measured in eleven healthy subjects with the Xsens InCap system and the Vicon 3D motion capture (MoCap) system at a self-selected speed. A novel integrated method that combines previous sensor-to-segment calibration approaches was developed for use in a MSK model with three degree of freedom (DOF) hip and knee joints. The following were compared: RMSE, range of motion (ROM), peaks, and R2 between InCap kinematics estimated with different calibration methods and gold standard MoCap kinematics. The integrated method reduced the RSME for both the hip and the knee joints below 5°, and no statistically significant differences were found between MoCap and InCap kinematics. This was consistent across all the different analyzed movements. The developed method was integrated on an MSK model workflow, and it increased the sensor-to-segment calibration accuracy for an accurate estimate of 3D joint kinematics compared to MoCap, guaranteeing a clinical easy-to-use approach.


Introduction
The assessment of locomotor function needs to be taken out of the lab and into the patient's ecological environment to evaluate the true impact of specific clinical treatments such as gait retraining or pain medication on joint function. To do so, joint kinematics and ultimately loading should be accurately measured in real-life conditions. An inertial capture (InCap) system combined with a musculoskeletal (MSK) model that represents the complex kinematics patterns and the mechanical interactions of the joint can be used for instance to evaluate patients' knee adaptations to instructed gait pattern modifications (toe-in/out, trunk leaning), aiming to decrease knee loading and eventually contributing to intermediate pain alleviation [1][2][3][4][5] and slowing down degenerative joint disease progression in the long term.
Laboratory-based 3D motion capture (MoCap) combined with in-ground embedded force plate systems has long been considered the gold standard method for capturing human movement and musculoskeletal function. Consequently, the use of 3D motion capture data as an input for MSK modelling workflows has become an attractive method to allow the estimation of different movement parameters (e.g., muscle forces and joint contact forces), however restricting the evaluation of locomotor function to highly controlled laboratory settings. In the context of degenerative joint diseases, the validity of MoCap based methods for the evaluation and the analysis of human movements may not reproduce the hip and knee sagittal plane movement and correlations from 0.36-0.76 [32]. However, frontal and transverse plane angle kinematics such as knee adduction, rotation, and ankle pronation were not evaluated. Furthermore, MoCap data is still required for defining the initial pose therefore limiting its use as a pure substitute for MoCap. A third calibration method is a functional calibration method that assumes a perfect alignment of the sensor and the joint rotation axes during specific movements (e.g., hip flexion-extension, knee flexion-extension, sit-to-stand, walking, and circling movement). This method has proven to be more appropriate for estimating frontal and transverse angles. For example, methods based on hip adduction or cycling motion show errors below 3.8 • for sagittal plane angles and below 6 • for non-sagittal plane angles during walking and cycling [33][34][35][36][37]. However, the assumption of pure rotational movement around the segment axes during calibration movement, imposes specific challenges in patients with a limited range of motion. The final method is anatomical calibration which uses external devices with integrated sensors to locate the joint axes with respect to the sensor axes system [38,39]. This anatomical calibration method shows errors <4 • for sagittal and non-sagittal angles except for ankle pronation (9.7 • ) during walking [40]. Although anatomical calibrations require supplementary tools and trained operators to properly locate the landmarks, this method can be applied to patients with limited joint mobility. Overall, only a few studies [35,41] combined static and functional calibration methods and reported RMSE < 5 • for sagittal plane and non-sagittal plane angle during walking only.
Advanced algorithms have been developed to overcome other current calibration accuracy limitations such as soft tissue motion artifacts caused by no-rigid connection between the body and the sensors [42,43] and magnetic distortions [9] caused by ferromagnetic objects that may affect the IMU sensor's orientation. For instance, closed-source commercially available models, such as the Xsens system can successfully compensate for magnetic distortions and drifts based on an in-house anatomical model and algorithms during any motion with an easy-to-use calibration strategy [27,44]. However, an open-source platform to compute 3D joint kinematics from IMU signals (independent from the InCap system used) for any type of motion, with no need of MoCap, would be highly beneficial, and it is currently unavailable.
This study aims to (1) develop an integrated functional sensor-to-segment calibration open-source workflow applicable in a clinical context and therefore independent of 3D MoCap for computing three-dimensional pelvis, torso, and lower limb joint kinematics (hip, knee, ankle) with IMU sensors based on an MSK model in the open-source software OpenSim; and (2) evaluate its accuracy compared to the gold standard MoCap system during walking, stair ascent and descent, and running. Ideally, the accuracy of such a workflow would allow estimation of kinematics within an accuracy of ±5 • , which allows discriminating joint kinematic of healthy and OA subjects, and it would allow for the measuring of differences in joint kinematics targeted during gait retraining. Furthermore, performance would be expected to be similar during a variety of different relevant, functional movements.

Data Collection
3D MoCap data was collected at the Movement & Posture Analysis Laboratory Leuven (MALL), at the Department of Movement Science, KU Leuven, Leuven, Belgium. This research was in accordance with the ethical guidelines provided by the ethical research committee KU Leuven (reference no. G-2021-3436). Eleven healthy adults with no musculoskeletal or neurological disorders volunteered for the study (6 males and 5 females; age: 26.0 ± 3.1 yr; height: 1.77 ± 0.9 m; weight 68.4 ± 9.9 kg). All participants provided written informed consent, prior to data collection.

Instrumentation
Full-body InCap data were collected using the Xsens MVN Awinda system (Xsens Technologies B.V., Enschede, The Netherlands), for which 17 IMUs were positioned based on the manufacturer guidelines (MVN User Manual). To record the motion at 60 Hz, Xsens MVN Analyze Pro 2021.0 was used. Simultaneously, a 10 infrared camera MoCap system (VICON, Oxford Metrics Group, UK) was used as the gold standard reference system to record the 3D position of the 54 reflective markers attached on anatomical landmarks of the different body segments (Plug-in Gait Reference Guide-Nexus 2.12 Documentation-Vicon Documentation) at 120 Hz ( Figure S1 and Table S1). Both systems were time synchronized based on the manufacturer's guidelines with a specific trigger at the start/stop recording time.

Protocol
3D InCap and MoCap data were collected, while each subject performed four barefoot movement tasks at self-selected speed: (a) 10 m of over ground walking, (b) 10 s of treadmill running, (c) ascending, and (d) descending stairs. For calibration, participants started each task with specific postures/movements: (i) N-pose, standing with arms to the sides, feet shoulder width apart, and facing forward for a period of five seconds, (ii) sitting in a wooden chair facing forward, (iii) knee extension followed by knee flexion for each leg (iv) standing in N-pose for three seconds, (v) hip abduction/adduction movement (three repetitions for each leg) with no movements in knee or ankle joints; the latter movement execution enforced similar angular velocities of the three segments-thigh, shank and foot, and (vi) participants performed the tasks (a-b-c-d). Raw data recorded by the InCap and the MoCap system were exported to .mvnx and .c3d files using Xsens MVN Analyzer Pro 2021.0 and Nexus 2.12, respectively, and they were imported into MATLAB (R2020a, MathWorks). For both systems, all movements were time-normalized (101 data points at time intervals of 1% of a gait cycle) and processed using custom-built MATLAB scripts.

Kinematics Analysis
3D MoCap kinematics were calculated using the inverse kinematics tool in OpenSim v4.2 using an MSK model with 21 segments and 41 degrees-of-freedom (Hamner 2010 [45]). The knee joint was modeled as a ball-and-socket joint (3 degrees-of-freedom). 3D InCap kinematics were calculated using the inverse kinematics tool in the OpenSense toolkit [32], with a newly developed and integrated sensor-to-segment calibration.

Sensor-to-Segment Calibration
This integrated sensor-to-segment calibration method (method 3-M3) can be used independent of 3D MoCap data, and it relies on a combination of previously published methods [33,[46][47][48] (Figures 1 and 2). Firstly, functional calibration of thigh, shank, and foot sensors was based on analytical hip abduction-adduction motion [33] to improve the estimation of sagittal and non-sagittal plane knee and ankle kinematics (method 1-M1). The reference frame misalignment of the thigh and the shank sensors was estimated by comparing the angular velocity vectors during a hip ab-adduction movement. Each sensor's angular rotation along the global vertical axis was estimated given that the angular velocity of the thigh and the shank sensors should be identical, as during the hip ab-adduction, the knee motion can be assumed to be approximately zero. Therefore, the thigh and the shank can be considered as a single segment with identical angular velocity. This procedure was then extended and applied to the foot segment. Secondly, functional calibration of the torso, pelvis, and thigh sensors was based on a walking and a sit-to-stand motion (STS) [48] to improve the estimation of sagittal and non-sagittal plane lumbar, pelvic, and hip kinematics (method 2-M2). Torso and pelvis sensors were aligned based on STS motion, whereas thigh sensors were aligned based on the walking motion under the assumption that these movements occur mainly around the global medio-lateral axis. To this end, principal component analysis (PCA) was used to determine the sensor rotation that maximized the movement around the medio-lateral axis. Finally, the initial pose of the model segments was estimated using the approach presented in Fasel et al. [47], with adaptations regarding the utilized OpenSim MSK model and Xsens coordinate reference systems [31,49].The segment global reference frame was defined with the vertical axis aligned with gravity vector, the anterior-posterior axis perpendicular to the gravity vector and pointing in the MSK model forward direction and the medio-lateral axis defined as the cross-product between the vertical and the anterior-posterior axes, pointing to the right of the MSK model. It was assumed that the torso, pelvis, and lower limb segments were orientated in the same direction while standing. The segments' inclinations were determined using the accelerometer data and the gravity vector. The initial joint angles of the MSK model were then computed based on the segment orientation.
Performance of the developed method M3 (combination of M1 and M2) was evaluated against the anatomical calibration that uses an anatomical position as an initial guess for the standard static calibration with the OpenSense toolkit (M0). Specifically, the anteriorposterior axis of the pelvis sensor was aligned with the anterior-posterior axis of the pelvis segment of the MSK model; then, torso, thigh, shank, and foot sensors were assumed to be aligned with that of the pelvis. For the calibration reference system, segment and joint coordinate systems were defined according to ISB recommendations [50,51]. motion (STS) [48] to improve the estimation of sagittal and non-sagittal plane lumbar, pelvic, and hip kinematics (method 2-M2). Torso and pelvis sensors were aligned based on STS motion, whereas thigh sensors were aligned based on the walking motion under the assumption that these movements occur mainly around the global medio-lateral axis. To this end, principal component analysis (PCA) was used to determine the sensor rotation that maximized the movement around the medio-lateral axis. Finally, the initial pose of the model segments was estimated using the approach presented in Fasel et al. [47], with adaptations regarding the utilized OpenSim MSK model and Xsens coordinate reference systems [31,49].The segment global reference frame was defined with the vertical axis aligned with gravity vector, the anterior-posterior axis perpendicular to the gravity vector and pointing in the MSK model forward direction and the medio-lateral axis defined as the cross-product between the vertical and the anterior-posterior axes, pointing to the right of the MSK model. It was assumed that the torso, pelvis, and lower limb segments were orientated in the same direction while standing. The segments' inclinations were determined using the accelerometer data and the gravity vector. The initial joint angles of the MSK model were then computed based on the segment orientation. Performance of the developed method M3 (combination of M1 and M2) was evaluated against the anatomical calibration that uses an anatomical position as an initial guess for the standard static calibration with the OpenSense toolkit (M0). Specifically, the anterior-posterior axis of the pelvis sensor was aligned with the anterior-posterior axis of the pelvis segment of the MSK model; then, torso, thigh, shank, and foot sensors were assumed to be aligned with that of the pelvis. For the calibration reference system, segment and joint coordinate systems were defined according to ISB recommendations [50,51]. and with all the other sensors-static calibration approach [46]; Method 1: hip abduction-adduction motion to align thigh-shank-foot sensors-functional calibration [34]; Method 2: PCA-method based on sit-to-stand and walking to align torso-pelvis-thigh sensors-functional calibration [46]; Method 3: combination of Methods 1 and 2 for the alignment of all sensors-functional calibration [34,46]. and with all the other sensors-static calibration approach [46]; Method 1: hip abduction-adduction motion to align thigh-shank-foot sensors-functional calibration [34]; Method 2: PCA-method based on sit-to-stand and walking to align torso-pelvis-thigh sensors-functional calibration [46]; Method 3: combination of Methods 1 and 2 for the alignment of all sensors-functional calibration [34,46].

Statistics
For multiple comparisons, the Wilcoxon rank test was implemented in MATLAB, and it was used to statistically assess significant differences in the root mean squared error (RMSE), the determination coefficient R 2 , range of motion (ROM) difference, and the maximum peak difference between each calibration method. Significance was set at p ≤ 0.05. A false discovery rate (FDR) statistical approach in multiple assumptions testing was used for multiple comparisons in order to correct for random events that falsely appear significant. A q-value threshold of 0.05 (FDR of 5%) among all significant variables was used. To evaluate the functional relevance of the observed differences between methods for the different joints, the error in kinematics between different InCap methods and MoCap was evaluated against reported differences in kinematics between control and OA subjects [52][53][54][55].

Statistics
For multiple comparisons, the Wilcoxon rank test was implemented in MATLAB, and it was used to statistically assess significant differences in the root mean squared error (RMSE), the determination coefficient R 2 , range of motion (ROM) difference, and the maximum peak difference between each calibration method. Significance was set at p ≤ 0.05. A false discovery rate (FDR) statistical approach in multiple assumptions testing was used for multiple comparisons in order to correct for random events that falsely appear significant. A q-value threshold of 0.05 (FDR of 5%) among all significant variables was used. To evaluate the functional relevance of the observed differences between methods for the different joints, the error in kinematics between different InCap methods and MoCap was evaluated against reported differences in kinematics between control and OA subjects [52][53][54][55].

Results
Complete tables and figures for comparison of all sagittal and non-sagittal plane kinematics are presented in the Supplementary Materials.

Walking:
Knee: Figure 3 shows a representative result of estimated knee angles with the integrated method M3 during the different movements. During walking, M3 significantly reduced errors in estimated frontal and transverse average plane offset compared to the other methods: M0 (8°), M1 (6°), M2 (4°), M3 (2°). These offsets with respect to MoCap

Results
Complete tables and figures for comparison of all sagittal and non-sagittal plane kinematics are presented in the Supplementary Materials.

Walking
Knee: Figure 3 shows a representative result of estimated knee angles with the integrated method M3 during the different movements. During walking, M3 significantly reduced errors in estimated frontal and transverse average plane offset compared to the other methods: . These offsets with respect to MoCap angular values can be noticed at the start of the curves in Figure 3. For M3, the overall average RMSE of the knee flexion, abduction, and internal rotation angle were 3.2 ± 1.2 • , 2.5 ± 2.2 • , and 4.5 ± 1.7 • , respectively, with average R 2 > 0.7 ± 0.2. ROM errors were 1.1 ± 0.6 • , 2.2 ± 1.8 • , and 2.9 ± 2.2 • , and peak magnitude errors were 1.6 ± 2.5 • , 3.2 ± 2.1 • and 4.3 ± 2.7 • for knee flexion, abduction, and internal rotation angle, respectively (see Tables 1 and 2 and 2). M3 showed the highest accuracy compared to the other InCap methods. Where average RMSE reduced significantly from 7.4 ± 2.4° to 2.1 ± 2.4° and from 8.4 ± 3.5° to 4.3 ± 2.0°, for knee adduction and rotation, respectively, compared to M0. Compared to MoCap, significant differences were shown in peak knee adduction and rotation in M0 and between M1 and M3 for knee rotation. RMSE averaged over all joint planes improved for M3: 3.2 ± 1.6° compared to M1: 5.0 ± 2.4° and M2: 5.1 ± 2.8°.    The RMSE between kinematics estimated from MoCap and the different methods (M0-M1-M2-M3) during walking are shown in Figure 4. M3 showed an average RMSE < 5 • for all joints and DOFs excluding ankle pro/supination and pelvis internal/external rotation. M0 showed an average RMSE above 5 • for all joints, M1 < 5 • for hip and knee adduction and M2 < 5 • for hip and knee flexion. Statistically significant differences were found between MoCap and InCap in M0 for hip adduction, knee adduction, knee rotation, pelvic tilt, and pelvic rotation during walking, with similar trends during running, and the ascent and the descent of stairs.

Running and Stair Ascent and Descent
For running (Table S3), accuracy relative to Mocap was lower compared to walking across all methods. In particular, in non-sagittal plane angles, M3 showed for hip, knee, and ankle frontal and transverse plane angles on average a RMSE 5.1 ± 2.9 • and 5.2 ± 2.1 • compared to 3.5 ± 1.9 • and 4.0 ± 1.5 • during walking. Similarly, stair ascent and descent movement (Tables S4 and S5) showed smaller accuracy compared to walking across all methods. In particular, lumbar and pelvic angles during stair ascent showed in sagittal, frontal and transverse planes on average a RMSE 5.5 ± 3.2 • , 5.3 ± 2.2 • and 5.2 ± 2.1 • compared to 3.5 ±1.8 • , 3.0 ± 0.9 • and 4.1 ± 1.8 • during walking (Figures S3-S5).

Summary and Main Findings
The aim of this study was to develop an integrated functional sensor-to-segment calibration applicable in a clinical context for computing three-dimensional pelvis, torso, and lower limb joint kinematics (hip, knee, ankle) during common daily life activities based on IMU sensors combined with an OpenSim MSK model. The developed method addresses previously documented challenges within literature, specifically sensor-to-segment calibration to estimate the body segment orientations from the sensor orientations without MoCap data and the integration with a complex MSK modeling workflow to estimate 3D joint kinematics [15,32,56,57]. In general, the integrated method based on the combination of PCA and functional calibration (M3) was for the sagittal, frontal, and transverse joint angles more in agreement with the gold standard MoCap (RMSE< 5 • -except transverse plane in ankle joint) than previously reported methods. These results confirm that calibration methods based on functional movements such as hip abduction/adduction, sit-to-stand, and walking improved the sensor-to-segment accuracy, particularly for non-sagittal plane kinematics.
In contrast to previous studies that only focused on specific activities [28,35,58], the integrated method (M3) can be used for evaluating different daily life movement tasks, such as over ground walking, running, and stair ascent/descent. This method does not require specific subject postures, the presence of expert operators, nor the use of external devices as some previous works proposed [38,40,59,60]. Moreover, the integrated method (M3) shows overall more accurate kinematics compared to previous calibration methods (i.e., M0, M1, M2 based on previous works) [33,36,46]. In fact, this method measures estimates of sagittal and non-sagittal plane kinematics that are only less than 5 • difference compared to MoCap during walking (where errors between 2 • -5 • are considered clinically acceptable). The M0 method based on anatomical/static calibration (M0) showed the lowest accuracy, particularly for non-sagittal plane angles, indicating that functional methods are required to improve the non-sagittal plane kinematic estimations. The method based only on hip ab-adduction calibration (M1) showed good accuracy for non-sagittal plane angles comparable with previous works [34,61] during walking but lower accuracy in sagittal plane angles, especially during stair ascent and descent. The walking based calibration method (M2) showed good accuracy for sagittal plane angles in agreement with previous work [48] but lower accuracy in non-sagittal plane angles across all motions. Therefore, method M3 was the most reliable method for accurate sensor-to-segment calibration, resulting in computing sagittal and non-sagittal angles joint kinematics with an error <5 • for out of the lab contexts compared to MoCap.
In addition, we evaluated the performance of the integrated method M3 by calculating the differences in ROM and in angle peaks compared to MoCap. The ∆ROM and ∆peak error of sagittal and non-sagittal plane angles varied between 2 • and 5 • . Such ROM and peak differences have been demonstrated to be important parameters for a variety of clinical applications [62,63]. For instance, ROM in knee osteoarthritis patients was reduced compared to healthy subjects during walking [19,55,64]. M3 accuracy would allow detecting previously reported differences between healthy and OA patients, which are in the order of 6-10 • in knee flexion ROM, 5-10 • in peak knee adduction, 5-10 • in hip adduction and rotation, and 5-8 • in pelvic list ROM, which is within the range of accuracy indicated in this preliminary investigation in healthy individuals. However, a detailed analysis and repetition of this research would be required for a dedicate knee OA cohort to confirm this, and it will form the basis of future work [52,53,55]. However, changes in pelvic list, tilt or knee adduction introduced during gait retraining such as toe-in or toe-out in the order of 1-3 • for pelvic list and tilt and 2-5 • for knee adduction may be difficult to detect, given the M3 average error being 3-5 • . Moreover, during stair ascent and descent the error is slightly over 5 • compared to MoCap, which in some cases can be considered not clinically acceptable.
Given that the InCap system is an easy-to-use system capable of accurately detecting the difference in ROM of 5 • , it is a valid alternative to conventional MoCap systems for the assessment of 3D joint kinematics in clinical and ecological settings. However, its accuracy decreases for frontal and transverse plane angles: in particular, for joints with small ROM such as knee adduction (2-5 • ), where a RMSE threshold of 5 • may represent an error of 100% of ROM. Here, correspondence in terms of correlation of the kinematics may yield a more accurate comparison for evaluating the accuracy of the InCap systems [38,53,65]. Consequently, we are aware that obtaining accurate non-sagittal plane angles with an InCap system during real-life conditions is complex, especially when the differences in ROM and peaks between for instance patient populations or conditions are <5 • . Nevertheless, due to the combination of functional calibration methodologies applied for developing the integrated M3, kinematic estimation accuracy improved across all angles and planes, with an ascent to levels comparable to MoCap. Furthermore, our developed method does not rely on manual offset correction as previously used during for instance stair ascent and descent [66,67].
Therefore, for applications where the specific clinical intervention, e.g., gait retraining in patients with knee OA, needs accurate measurement of knee non-sagittal plane angles the developed functional calibration approaches in combination with an MSK model with multiple degrees-of-freedom could be used if future studies confirm similar levels of accuracy in OA patients [68,69]. However, the selected calibration method depends on the degree-of-freedom of the joints in the MSK model. Indeed, an MSK model with one degree-of-freedom for the knee joint and the ankle joint as previously described for use with OpenSense, does not need functional calibration and, in this case, methods based only on anatomical/static calibration are sufficient to accurately estimate sagittal plane joint kinematics [32].
Our study showed that the accuracy of the developed calibration method is higher than previously developed methods-M0, M1, M2-and that the effect is consistent across multiple movement tasks. Nevertheless, there seems to be a negative effect on the calculated error due to the more dynamic movement tasks such as running. This may be due to the limited sample frequency used by the InCap system, which was only 60 Hz. Therefore, future studies should consider a higher sample frequency during such high dynamic movements. Another limitation of the presented method is that it relies on the ability of the subjects to perform specific joint movements-the hip abduction/adduction movement, and the sit-to-stand movement. Therefore, functional calibration may be challenging in subjects with limited joint mobility or impairments. Furthermore, the repeatability of sensor placement was not investigated in this study. Sensor placement repeatability is operator dependent and the introduced misalignment is known to account for up to a 55% error in the sagittal plane angles, while the sensor-to-segment calibration misalignment is known to account for up to 40-50% in non-sagittal plane angles [26,70]. By introducing the proposed sensor-to-segment calibration, sensor misalignment, including that introduced due to between operator variability in sensor placement, is expected to be reduced as this will automatically be accounted for during the functional movements [71]. Moreover, the integrated method was applied during movement tasks with a duration under two-minute motions, hence, it may be necessary to mitigate errors due to drift during movement tasks with a longer duration. However, method 3 may be applied with drift-reduced estimation as presented in the Fasel et al., 2018 study [36]. Drift is estimated and reduced based on the accelerations' vector in the global frame.
This study impacts the use of wearable sensor technology for documenting 3D knee joint kinematics (e.g., patients with knee OA). However, to evaluate the impact of specific clinical treatments out-of-the-lab such as gait retraining or pain medication, the knee loading should ideally be accurately measured in real-life conditions, e.g., to estimate the impact of the patient's gait pattern modifications (toe-in/out, trunk lean) on knee loading [2,4,[72][73][74][75]. Future work should therefore validate the method developed for the individual activities of daily living with more natural activity sequences, and it should focus on the estimation of kinetics, in particular the ground reaction forces and joint moments estimation to enable true real-world monitoring of the impact of feedback and gait interventions on an individual patient's locomotor function and joint loading. Such an approach would require valid and reliable methods for the 3D plane angle estimation [76] but also numerical methods to estimate ground reaction forces and moments without force platforms, i.e., solely based on kinematics data combined with a model of the footground contact [77,78]. In this context, machine learning based and probabilistic methods are now being introduced to estimate ground reaction forces and moments using InCap systems [79][80][81][82].

Conclusions
In conclusion, in this study we developed a novel integrated sensor-to-segment calibration method for IMU sensors that improves non-sagittal plane kinematic estimates in a range of movements. Moreover, keeping the focus on future clinical applications, this method aimed to calculate the 3D kinematics of the torso, pelvis, and lower limbs, allowing for a not-predefined location of the IMU-sensors thereby guaranteeing a clinical easy-to-use approach.