1. Introduction
Human posture and movement analysis is fundamental for a wide range of applications and disciplines, such as physical and neuro-rehabilitation, sports medicine, human performance assessment and virtual training. Although standard motion analysis instruments are widely used in these fields, the development and validation of ambulatory and unobtrusive sensing systems, which provide a reliable measurement of human motion and activity in out-of-lab daily life, are still open issues in the current literature. The development and validation of wearable technologies aimed at allowing physicians and therapists to remotely supervise and coach the patients during their rehabilitation exercise in their recovery phase are some of the current challenges.
Recently, as anticipated by De Rossi and Veltink in [
1], wearable textile-integrated sensing of human movements and electrophysiological acquisition devices have been developed. The authors of the current paper collaborated, within the EU project INTERACTION [
2], in the development of an e-textile platform endowed with wireless 3D inertial movement and textile-based stretch and goniometric sensing [
3]. The INTERACTION system was tested in post-stroke patients to assess the quality of their mobility and their reaching and grasping capacities.
Over the last decade, inertial measurement units (IMUs), based on the integration of accelerometers, gyroscopes and magnetometers, have gained the reputation of being at the cutting edge of wearable motion tracking [
4,
5,
6,
7,
8]. IMUs estimate the orientation of the body segments where they are attached by combining multi-sensor information through dedicated optimal sensor fusion algorithms mainly based on Kalman filtering [
4,
5,
7]. The general approach is to apply strap-down integration of the gyroscope signal [
9] and to correct the inclination and heading drifts through the accelerometer and magnetometer measurements. The combination of different IMUs, placed on connected body segments, and the additional information on the kinematic constraints enable joint angles to be measured [
10,
11]. IMU-based systems can also be employed to track the velocity of the human body, as demonstrated in the study performed by Yuan [
12].
Several wearable motion analysis systems based on IMU technology are now on the market. The well-known XSens MVN system [
13,
14] uses proprietary wireless motion trackers (MTw) applied to an adherent suit to perform full-body motion tracking (upper/lower limbs, trunk and neck). Each MTw includes a tri-axial accelerometer, a triaxial gyroscope and a tri-axial magnetometer and computes the fusion algorithms to calculate the global orientation of the body segment to which it is attached. The orientations of the single body segments are further fused by means of proprietary bio-mechanical models to estimate the full kinematics of the human body. XSens systems, which were used in the current study to benchmark the proposed technique, have been adopted and tested in many studies. For example, Zhang [
15] showed a mean error below two degrees by comparing the knee flexion-extension detected by two MTws with the same angle retrieved by a stereophotogrammetric system. Similar IMU-based systems focusing on gait analysis [
16,
17] have been marketed by INSENCO Co. (composed of nine inertial sensors placed on the lower limbs combined with wireless six-axis force sensors [
18]) and by Tec Gihan (the M3D system that combines inertial measurement units and wearable force plates [
19]).
Beside their widespread adoption in the field, IMU-based systems suffer from a loss of accuracy due to magnetic disturbance; they can be bulky, and their cost can be considerable. The contribution of inertial and translational accelerations can limit the accuracy of the accelerometer update [
20]. Most importantly, the presence of ferromagnetic materials or other magnetic disturbances due to environmental noise, which is inevitable during daily life applications, limits the accuracy of the heading compensation [
21,
22] and can compromise the full orientation estimation. Considering our context of wearable unobtrusive sensing in daily life, an ambulatory, yet reliable monitoring system that avoids the use of magnetometers would be beneficial, as already pointed out by Seel in [
23].
Thus, the new approach described in this work exploits only the acceleration components of IMUs in combination with other measurement sources taken from a particular class of e-textile sensors.
Textile-based or e-textile solutions have been developed, and their possible application for ambulatory and unobtrusive motion detection has been described in several research papers, including [
24,
25,
26,
27,
28]. Textile-based solutions have several advantages compared to solid-state sensors: low cost, lightweight, low thickness, flexibility and the possibility of adaptation to different body structures. By exploiting these kinds of sensors, it is possible to design sensing garments with sensor strips applied to specific locations on normal cloth. Despite these attractive characteristics, textile sensor adoption is still limited mainly due the low reliability, which thus limits their use to the reconstruction of large and slow movements. We recently developed a new generation of textile-based goniometers, obtained by coupling two layers of knitted piezoresistive fabrics (KPF) through an electrically-insulated layer, as described in [
29]. Compared to previously-developed e-textile sensors, textile goniometers can provide an accurate and reliable measurement of the angles between connected body segments [
29,
30,
31]. Although textile goniometers represent an important step forward in wearable human motion detection, they still have their limitations. In particular, the sensor calibration, which has to be performed in two known angular positions, as described in
Section 2.1, needs to be computed once the sensor has been integrated into the textile and worn by the user. Due to the intrinsic nature of the textile materials, which tend to modify their mechanical properties over time, the calibration coefficients can slowly drift, and the calibration procedure needs to be repeated periodically. This is highly undesirable in ambulatory and daily life applications where the subject has to use the system as is, without performing too many operations for the system to work properly. Moreover, we pointed out in [
29] that although KPF goniometers represent a consistent step forward in motion sensing through e-textiles, a slight hysteresis can still affect the sensor performance.
The aim of this paper is to show that by combining e-textiles and low cost inertial sensing, unobtrusive and reliable human motion monitoring is possible. This can be achieved by addressing the drawbacks of the two measurements, by fusing the information derived from KPF goniometers and tri-axial accelerometers. The idea behind this is that single triaxial accelerometers are less expensive and complex compared to the current complete IMU technology and do not suffer from magnetic interference. In addition, they can be easily attached or integrated into the extremities of the KPF goniometer in order to obtain a reliable and unobtrusive measurement system.
Without loss of generality, we demonstrated our concept on the movement of the femoral-tibial joint, which along with the hip and ankle is one of the three major articulations of the lower limb that allows ambulation.
Several research groups have investigated performing reliable knee joint kinematic assessments through simplified and reduced sensing systems based on accelerometers and gyroscopes and avoiding the use of the magnetometer information. Relevant studies are reported in [
23,
32,
33,
34,
35,
36,
37] and exploit the inertial information derived by two accelerometer/gyroscope pairs that are attached to the thigh and shank segments. In most of the cited publications, the reconstruction technique is based on the estimation of thigh and shank orientation to calculate the knee flexion-extension angle [
33,
35,
36]. The accuracy of these techniques relies on the precision of the orientation estimation, which is a challenge due to the absence of the magnetometer measurement, and the good alignment of the accelerometer/gyroscope reference frames with respect to the body segment reference frames. Significant results, root mean square error (RMSE) around two degrees in comparison with an optical system, have been obtained by Favre [
33] through an alignment procedure that exploited a static standing posture and predefined hip abduction/adduction movements. Another interesting work was recently presented by Seel in [
23], who proposed a new method for joint axis identification that does not rely on the accuracy with which the subject performs predefined postures or movements, as well as an innovative reconstruction technique. Their technique is based on the integration of the angular rate on the identified joint axis and the subsequent correction of the slowly-drifting angle by the accelerometer measurement (RMSE of three degrees in comparison with an optical reference system).
In this work, we propose an innovative solution for the implementation of a simplified and reduced sensing system for knee joint flexion-extension assessment. We exploit a hybrid system that combines the non-inertial angular measurement of the textile goniometer with the inertial information derived from the two accelerometers. More specifically, we developed and conducted a preliminary evaluation of a dedicated sensory fusion approach, based on a Kalman filter implementation. This method exploits the accelerometer measurement to correct the angular estimation performed by the KPF goniometer by continuously adjusting its calibration parameters. We demonstrated the capabilities of the new hybrid measurement system, which is made up of KPF goniometers and tri-axial accelerometers, through the experimental setup described in
Section 2. We built a sensing knee band prototype with pockets used to host a KPF goniometer, which detects the angle between femur and tibial segments, and two IMUs attached to the the thigh and shank segments. Raw accelerometer data and KPF outputs were processed through the fusion technique (see
Section 3) to estimate the knee flexion-extension angle. At the same time, the whole information set derived from the thigh and shank IMUs (
i.e., the full orientation matrix extracted by fusing accelerometer, gyroscope and magnetometer information) was employed as the standard reference for benchmarking our new approach. In the experimental session, the prototype was employed to acquire goniometer and IMU data while subjects performed different tasks involving moving their knee joints. To introduce the sensing fusion technique, the KPF device and inertial sensing system working principles are described in
Section 2.1 and
Section 2.2, together with the calibrations of the two systems needed to start the experimentation. The calibration procedure in
Section 2.1, which identifies the parameters of the KPF goniometer, is fundamental to determine the initial status of the Kalman filter procedure described in
Section 3. The main part of the paper describes the algorithm that improves the quality of the knee position estimation, by refining the goniometer parameters. Finally, in
Section 4, the knee flexion-extension estimate performed by the hybrid sensing system, detected in six different movement tasks, is compared to the reference measurement system to evaluate its performance. In the whole set of tests, the fusion technique showed a reliable performance and was effective at measuring dynamic knee movements (across trials, mean and standard deviation of the RMSE equal to
and
degrees, respectively). These results make our approach comparable with other ambulatory measurement systems. The results also demonstrated that the fusion method significantly improved the measurement performance of the individual sensor classes (accelerometers and goniometers considered individually
vs. the reference measurements). In addition, note that no hypotheses were made on the accelerometer position and orientation with respect to the joint reference frame. Note also that the system can be easily calibrated, through the dedicated procedure conceived of within this work, without requiring the user to perform complex tasks. These achievements and the promising measurement performance make the hybrid system a good candidate for ambulatory and unobtrusive wearable measurement systems.
2. Experimental Setup
In order to acquire data concerning the flexion-extension of the femoral-tibial joint, both from KPF goniometers and inertial sensors, a knee band prototype was created using a Lycra fabric. The band was specifically designed with pockets to contain the KPF and inertial sensing, as shown in
Figure 1a. As shown in in
Figure 1b, the knee was simply modeled as a hinge joint, and the flexion-extension angle (
θ) was defined as the angle between the two consecutive segments of the model (
i.e., the angle between the
x unit vectors of the
and
frames).
The KPF goniometer (length 30 cm, width 2.5 cm and thickness 1.5 mm) was integrated into the garment to entirely cover the considered joint. The working principle of the goniometer, produced by SMARTEX [
38] and based on our proprietary design, is reported in
Section 2.1. The goniometer was calibrated before being integrated into the knee band according to the procedure described in
Section 2.1.
In the same knee band, two inertial measurement units (MTw by XSens, [
13]) were fixed inside the pockets at the thigh (
) and the shank (
). IMUs are positioned on the knee band without any assumptions regarding their orientation and position with respect to the knee joint segments. The XSens MTws [
13] used in the experiments are compact IMUs containing 3D linear accelerometers, 3D rate gyroscopes and 3D magnetometers. Each MTw performs real-time signal elaboration and transmits, through a wireless link, the 3D orientation (Euler, quaternions or rotation matrices) and raw sensor data (acceleration, angular velocities, Earth magnetic field orientation). The raw accelerometer component available from the IMU signal was used in combination with the textile goniometer, as detailed in
Section 3, hereafter named the accelerometer sensor. The knee flexion-extension angle derived from the IMU outputs can be computed by taking into account the components of the full rotation matrix, which describes the orientation of the IMU frame on the shank with respect to the frame of the IMU placed on the thigh [
10,
14]. The IMU devices were calibrated according to the procedure provided by the producer [
14]. The resulting knee flexion extension angle (
) was considered as the reference measurement to evaluate the textile goniometer/accelerometer hybrid system that we propose in this work.
Figure 1.
Sensing prototype and sensor placement around the knee joint. (a) A double-layer knitted piezoresistive fabrics (KPF) goniometer and two inertial measurement units (IMUs) were applied to the knee band; (b) A geometrical scheme of the reference frames fixed with the body segments and the IMUs. The knee was simply modeled as a hinge joint, and the flexion-extension angle (θ) was defined as the angle between the two consecutive model segments (i.e., the angle between the x unit vectors of the and frames). IMU and accelerometer reference frames ( and ) are not aligned with the corresponding segment reference frame.
Figure 1.
Sensing prototype and sensor placement around the knee joint. (a) A double-layer knitted piezoresistive fabrics (KPF) goniometer and two inertial measurement units (IMUs) were applied to the knee band; (b) A geometrical scheme of the reference frames fixed with the body segments and the IMUs. The knee was simply modeled as a hinge joint, and the flexion-extension angle (θ) was defined as the angle between the two consecutive model segments (i.e., the angle between the x unit vectors of the and frames). IMU and accelerometer reference frames ( and ) are not aligned with the corresponding segment reference frame.
In order to collect data on the knee kinematics, in the experimental sessions, the band was worn by five healthy subjects, who were asked to perform several tasks, grouped into a set of six trials and repeated three times. In the first two trials, the subjects were asked to perform knee flexion-extension movements at different velocities in monopodalic contralateral standing. During these phases, the subjects were asked to flex/extend the knee approximately in the 0–
range. The last four trials concerned walking activities with increasing velocities. The subjects were asked to walk freely on a straight line for about 40 s, starting with slow-speed walking (first trial) up to very fast walking (fourth trial). The subjects were free to select the velocities without external constraints. Before starting to collect the full set of data, each subject was asked to perform the simple movements necessary to align the thigh and shank accelerometer frames following the calibration protocol described in
Section 2.2.
KPF goniometer output was gathered through a dedicated wireless electronic unit, and data were sampled at 100 Hz and transmitted to a remote PC for elaboration. The IMU sampling rate was set to 100 Hz, and the data (raw accelerometer outputs and orientation expressed as the rotation matrix) were transmitted via a wireless link to the same PC. IMU and goniometer data were synchronized through a dedicated XSens software tool. After recording the sensor outputs, our algorithm for the estimation of the knee flexion-extension angle, which fuses data from the KPF goniometer and the accelerometers (described in
Section 3), runs offline on a remote PC. The knee flexion-extension angle (
θ), defined in
Figure 1b, was extracted by the fusion method applied to goniometer and accelerometer outputs with the same update rate of the input signals (100 Hz).
2.1. KPF Goniometers
This section briefly reports the background and the working principle of the KPF goniometer used in our experiments. As described in our previous studies [
29,
30,
31] and reported in
Figure 2, our textile goniometers were developed by coupling two piezoresistive layers through an electrically-insulating layer. The sensing layers were made of knitted piezoresistive fabrics (KPFs), previously employed as single-layer strain sensors for monitoring biomechanical and cardio-pulmonary parameters [
39,
40].
Figure 2.
Schematic diagram of a double-layer KPF goniometer. The black stripes represent the two identical piezoresistive layers, while the gray stripe is the insulating layer. When the sensor is in the flat position, the resistance difference () between the two layers is zero. When the sensor is flexed, is proportional to the bending angle (θ), defined as the angle between the tangent planes to the sensor extremities (green dashed line in the picture).
Figure 2.
Schematic diagram of a double-layer KPF goniometer. The black stripes represent the two identical piezoresistive layers, while the gray stripe is the insulating layer. When the sensor is in the flat position, the resistance difference () between the two layers is zero. When the sensor is flexed, is proportional to the bending angle (θ), defined as the angle between the tangent planes to the sensor extremities (green dashed line in the picture).
Ideally, if the two KPF layers were geometrically and electrically equivalent, the sensor output, represented by the resistance difference between the two sensing layers (
), vanishes when the device is flat and is proportional to the flexion angle (
θ), except for a second order infinitesimal function [
29]:
In practice, Equation (
1) is not verified, due to the differences in the electrical properties between the two piezoresistive layers. In [
29], the
θ vs.
relation can be reasonably approximated by the following linear function:
where
and
represent the goniometer sensitivity and offset, respectively. The angle values
θ can be explained by Equation (
2) as:
In Equation (
3), parameters
and
remain unknown, and it is necessary to perform an initial calibration to determine them. Note that
and
are the process variables in the Kalman filtering described in
Section 3.1 and that their initial estimation is thus the zero status for the iterative process of data refinement. According to Equation (
3), by acquiring
and
in two different angular positions
and
, it is possible to compute
,
as:
Before starting to collect the data in our trials, the goniometer was calibrated in two angular positions ( and ) to obtain the and the values to feed the Kalman cycle.
2.2. Accelerometer Alignment
The data fusion algorithm, which is the core of the current research and is described in
Section 3, combines the information from the KPF goniometer and the triaxial accelerometers fixed on the thigh and shank (raw accelerometer components of the IMU signal). Considering the unknown and unmeasured orientation and the position of the accelerometers with respect to the knee joint segments (as shown in
Figure 1a,b, the accelerometer frames need to be aligned with the corresponding joint segments (
,
).
For this reason, we designed the sensor-to-segment calibration procedure outlined in
Figure 3. The required re-orientation of the accelerometer frames can be obtained as an accelerometer calibration protocol, which is performed on-body and requires knowledge of the goniometer output. The only assumption is that the position/orientation of the accelerometer frames with respect to the joint segments remains constant (
i.e., we neglect muscle artifacts that could otherwise be reduced with a strategic positioning of the sensors).
In the first step, using the accelerometer output in a static position, it is possible to determine the transformation that maps the
,
axes of the two frames into the corresponding axes of the segment frames. This transformation is performed by computing the
and
angles that can be obtained from the output vector of the uncalibrated accelerometer
, as:
Figure 3.
The calibration procedure for the accelerometers. In the first step, using the accelerometer output only, acquired from a subject in standing position, the axes of the accelerometer frame are aligned with the corresponding axis of the bone frames by computing the and angles. In the second step, using data collected by the goniometer and the accelerometers in a dynamic acquisition, the remaining axes of the inertial frames are aligned with the corresponding axes of the frames fixed with the joint segments.
Figure 3.
The calibration procedure for the accelerometers. In the first step, using the accelerometer output only, acquired from a subject in standing position, the axes of the accelerometer frame are aligned with the corresponding axis of the bone frames by computing the and angles. In the second step, using data collected by the goniometer and the accelerometers in a dynamic acquisition, the remaining axes of the inertial frames are aligned with the corresponding axes of the frames fixed with the joint segments.
To gather static information from the accelerometers, the subject was asked to stand upright in a natural position. In this position, the
’s are aligned with the gravity vector, except for the physiological angles of the leg bones in standing position, which are negligible with respect to the range of motion of the femoral-tibial joint. In the second step of the alignment procedure, starting from the upright position, the subject was asked to slowly flex the knee three times (in monopodalic standing). Let us suppose that the movement is spanned by a time variable
. Then, the angles
and
, necessary to complete the frame re-orientation, are computed as follows:
where:
describes a rotation around the
,
axis,
are the accelerometer outputs after the first calibration step,
is the flexion angle obtained by the nominal goniometer output and
is the rotation of the angle detected by the goniometer along the knee anatomical rotation axis (parallel to
in
Figure 1b). The calibrated accelerometer outputs (
) are then obtained by simply rotating
by the corresponding rotation matrix
.
To test the reliability of the accelerometer calibration procedure, let us introduce the following metric:
Where R, S are a special orthogonal matrix and I the identity in
. The misalignment function
is a real map defined on orthogonal frames (the columns of
R), and it is possible to prove, thanks to the triangle inequality, that its values lie in the range
. For
, matrices
R and
S are completely aligned, while for
, matrices
R and
S are “intuitively” completely misaligned. Function
provides an indication of the reliability of the presented calibration method for the accelerometers. In particular, we adopted the
metric to compare the accelerometer sensor-to-segment calibration matrices with the corresponding IMU calibration matrices obtained by the standard procedure provided by the producer. Below is the result of the alignment for the accelerometer and the IMU frame placed on the shank. The sensor-to-segment calibration matrix derived from the XSens MTw is given by:
While the accelerometer calibration matrix derived by the method described above and built using (
) holds:
Since IMUs and accelerometers refer to the same reference frame, the misalignment before calibration can be determined considering, in the
map, no transformations applied to the accelerometer signal (
i.e.,
) and the same IMU sensor-to-segment calibration matrix S. In this case, the misalignment function holds:
While using the matrix
R generated in the calibration procedure, we obtain:
Which ensures that and are, in practice, aligned. Similar results were obtained for the accelerometer alignment across the whole set of trials carried out in this work.
4. Results
In this section, the knee flexion extension angle estimated by the hybrid system, obtained by fusing goniometer and accelerometer data as described in
Section 3, is compared to the output of the reference measurement system (knee flexion extension angle obtained by the IMUs). In addition, to highlight the improvement of our fusion technique, a comparison between the reliability of our method (hybrid system
vs. IMUs) and the errors committed by the individual sensors is presented (
i.e., accelerometers
vs. IMUs and textile goniometer
vs. IMUs).
Figure 5 compares the results of our estimation technique and the reference measurement system in two representative plots taken from the slow and fast monopodalic flexion trials (
Figure 5a,b).
Figure 6 shows the same comparison for representative plots of the four walking conditions ordered from the slowest (
Figure 6a) to the fastest (
Figure 6d). As can be observed from
Figure 5 and
Figure 6, our estimation technique shows a reliable performance and a good capability to follow the knee movement, both in slow and fast tasks.
Figure 5.
Dynamic comparison between our estimation technique and the reference measurement during contralateral monopodalic standing tasks with different knee flexion-extension velocities. The blue line represents our estimation, while the red line is the reference measurement. (a) Slow knee flexion; (b) Fast knee flexion.
Figure 5.
Dynamic comparison between our estimation technique and the reference measurement during contralateral monopodalic standing tasks with different knee flexion-extension velocities. The blue line represents our estimation, while the red line is the reference measurement. (a) Slow knee flexion; (b) Fast knee flexion.
For the six trials considered (slow flexion, fast flexion and the four walking conditions with increasing speed), we calculated the deviation between the estimation performed by the proposed technique applied to the hybrid system (
) and the reference measurement (
extracted from the IMUs) in terms of root mean square error (RMSE). Likewise, we evaluated the RMSE of the KPF goniometer (
) and accelerometer (
) estimations against the reference measurement system. The relations describing the RMSE deviations are reported in the following equations:
where i represents the
i-th trial (
slow flexion,
fast flexion,
the slowest walking,
the fastest walking),
j denotes different subjects,
is the number of subjects,
is the number of samples of the
i-th trial for the
j-th subject and the suffix
m indicates the different measurement systems. In particular,
holds for the hybrid system,
for the goniometer, and
refers to the accelerometers. Data obtained from the statistic Equation (
29) are reported in
Table 1 with the average value
μ and the standard deviation
σ computed over the six trials (last two columns of
Table 1).
Figure 6.
Dynamic comparison between our estimation technique and the reference measurement during walking at different velocities. Velocities increase from (a) to (d). The blue line represents our estimation, while the red line is the reference measurement.
Figure 6.
Dynamic comparison between our estimation technique and the reference measurement during walking at different velocities. Velocities increase from (a) to (d). The blue line represents our estimation, while the red line is the reference measurement.
Table 1.
Root mean square errors (RMSEs) obtained in the six trials for the various estimation methods. The first row () reports the errors of the hybrid system by applying the fusion technique described in this paper. The second and third rows show the textile goniometer () and accelerometer () errors. The last two columns report the mean and standard deviation of the RMSE across the trials for the different measurement systems.
Table 1.
Root mean square errors (RMSEs) obtained in the six trials for the various estimation methods. The first row () reports the errors of the hybrid system by applying the fusion technique described in this paper. The second and third rows show the textile goniometer () and accelerometer () errors. The last two columns report the mean and standard deviation of the RMSE across the trials for the different measurement systems.
| Slow Flexion | Fast Flexion | Walking No. 1 | Walking No. 2 | Walking No. 3 | Walking No. 4 | Average | Standard |
---|
| (Slowest) | (Fastest) | Value μ | Deviation σ |
---|
| 0.97 | 3.50 | 1.07 | 1.6 | 2.1 | 2.5 | 1.96 | 0.96 |
| 5.12 | 4.60 | 5.40 | 4.6 | 5.5 | 5.7 | 5.15 | 0.47 |
| 1.48 | 10 | 5.80 | 6.7 | 7.1 | 8.2 | 6.55 | 2.87 |
5. Discussion
Our hybrid system’s reliable performance is confirmed by the results reported in the first row of
Table 1. The mean and standard deviations of the RMSE are
and
, respectively. The minimum deviation was detected in the slow flexion trial (
), while the maximum error was in the fast flexion experiment (
). The deviations in the walking tasks range from
for the slowest task to
for the fastest. The first row of
Table 1 and
Figure 5 and
Figure 6 highlight that in the hybrid system, RMSE increases as the speed of movement increases. In addition, fast walking trials show a lower error than fast flexion trials in monopodalic standing (as can be seen by comparing
Figure 5b and
Figure 6d and from the related RMSEs in
Table 1). This is reinforced by the fact that the maximum error was in the fast knee flexion-extension trial reported in
Figure 5b. This issue can be well explained by considering the estimation technique described in
Section 3.1. Indeed, in the fast flexion trial, performed at maximum velocity in the range
–
, the rotational acceleration measured by the shank accelerometer is far from being small with respect to gravity and translational contributions. The smaller error in the walking trials may be explained considering that, during walking, even at high speeds, the rotational acceleration is sufficiently small. A comparison with the current literature reveals that joint angle measurements with commercial solid state electrogoniometers, widely used for the ambulatory evaluation of the range of motion and movement frequency/velocity/acceleration of the joints for both clinical and occupational evaluations [
47,
48], shows errors greater than two degrees, with a strong dependence on the sensor positioning and on the cross-talk between joints [
49]. The reduced sensing systems, described in
Section 1 [
23,
32,
33,
34,
35,
36,
37] and based on two accelerometer/gyroscope pairs, achieved average RMSEs in the range
–
.
We performed a further analysis to quantify the improvement introduced by the sensor fusion technique with respect to the single sensing subsystems (
i.e., accelerometers and goniometer applied individually without fusing the information). The RMSEs related to the goniometer estimation of the knee flexion-extension angle are reported in the second row of
Table 1. The mean and standard deviations of the RMSE are
and
(in accordance with the results we obtained in a previous work on textile goniometer angular measurement [
29]). In addition, given how the KPF goniometer works, the RMSE variability across the trials, performed at different velocities, is quite low. The accelerometer RMSEs are reported in the third row of
Table 1. The mean and standard deviation of the error hold
and
, and as expected, there is a strong dependence on the execution velocity of the trial (from
for the slow flexion, to
for fast flexion). The last two columns of
Table 1 highlight that the error introduced by the hybrid system (
,
), evaluated in terms of RMSE, is considerably smaller than both the average errors of the goniometer (
,
) and the mean of the accelerometers (
,
). In addition, the error of the hybrid system is smaller than the goniometer and accelerometer systems in each of the tasks considered.
Figure 7 highlights the angle estimation
vs. the reference for the fusion method made on the hybrid system and the corresponding estimation made with the goniometer and accelerometers individually (representative plots of monopodalic flexion and walking tasks).
In the accelerometer signal, the errors are due to the acceleration spikes that are present when changing movement direction and when interacting with the environment (e.g., when the foot hits the floor in walking) and, thus, justify the large mean RMSE and the related standard deviation. Conversely, the goniometer signal is affected by errors due to the limitation introduced by the use of the nominal calibration parameters and slight hysteresis, which produced the error reported in
Table 1. In conclusion, our fusion method compensates for these errors and consistently improves the accuracy of the estimation with respect to the single sensing subsystems taken individually.
Figure 7.
Signal comparison between the angle reconstruction by the accelerometers (green dotted line), the goniometer (black dotted line), the hybrid system (accelerometer + goniometer, blue solid line) and the reference measurement. (a) Walking; (b) Fast flexion.
Figure 7.
Signal comparison between the angle reconstruction by the accelerometers (green dotted line), the goniometer (black dotted line), the hybrid system (accelerometer + goniometer, blue solid line) and the reference measurement. (a) Walking; (b) Fast flexion.
One limitation of the current study is that the hybrid system and the fusion technique were evaluated against inertial measurement units. It is well known that the golden standard in biomechanics is marker-based optical tracking. Optical motion capture systems based on external markers have been extensively used in gait analysis and have shown, for the knee flexion-extension angle, RMSE deviations in the order of
in comparison with bone fixed markers [
50]. However, we decided to evaluate our system with respect to the IMU reference, so that we would be able to track the subject’s activity in an unconstrained environment (free walking tasks whose duration was about 40 seconds). In any case, we have now begun evaluating the hybrid system during the daily life activity of subjects where the IMU reference represents the only valid option. On the other hand, as demonstrated in [
15], the estimation of the knee joint angle through the IMU-based system we used has shown reliable performance (errors in the order of two degrees).
A second limitation, which we are currently addressing, is that only one degree of freedom can be monitored per joint. We did not consider the rotation of the knee on the horizontal plane, since KPF goniometers are not influenced by twisting. From [
29], where the complete form of relationship Equation (
1) is reported in Taylor expansion form, it is possible to verify that there is no second order derivative term and, consequently, that the goniometer output
does not depend on its torsion. To detect knee torsion, another goniometer is needed, which is placed in a crossed position, with one extremity on the thigh and the other one on the shank. In our solution, the global curvature of the second goniometer is related to knee torsion. We are currently working on the modeling of this sensing configuration and on the relative calibration strategy.