Next Article in Journal
Efficient High-Dimensional Quantum Key Distribution with Hybrid Encoding
Next Article in Special Issue
New Construction of Maximum Distance Separable (MDS) Self-Dual Codes over Finite Fields
Previous Article in Journal
Uncertainty Assessment of Hyperspectral Image Classification: Deep Learning vs. Random Forest
Previous Article in Special Issue
A New Efficient Expression for the Conditional Expectation of the Blind Adaptive Deconvolution Problem Valid for the Entire Range ofSignal-to-Noise Ratio
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quaternion Entropy for Analysis of Gait Data

by
Agnieszka Szczęsna
Institute of Informatics, Silesian University of Technology, Akademicka 16, 44-100 Gliwice, Poland
Entropy 2019, 21(1), 79; https://doi.org/10.3390/e21010079
Submission received: 9 December 2018 / Revised: 9 January 2019 / Accepted: 15 January 2019 / Published: 17 January 2019
(This article belongs to the Special Issue Information Theory Applications in Signal Processing)

Abstract

:
Nonlinear dynamical analysis is a powerful approach to understanding biological systems. One of the most used metrics of system complexities is the Kolmogorov entropy. Long input signals without noise are required for the calculation, which are very hard to obtain in real situations. Techniques allowing the estimation of entropy directly from time signals are statistics like approximate and sample entropy. Based on that, the new measurement for quaternion signal is introduced. This work presents an example of application of a nonlinear time series analysis by using the new quaternion, approximate entropy to analyse human gait kinematic data. The quaternion entropy was applied to analyse the quaternion signal which represents the segments orientations in time during the human gait. The research was aimed at the assessment of the influence of both walking speed and ground slope on the gait control during treadmill walking. Gait data was obtained by the optical motion capture system.

1. Introduction

The parameters associated with chaos are measures of dimension, rate of information (entropy) and the Lyapunov determinant. Kolmogorov entropy K is known as a chaos metrics and the value of entropy can be used for the classification of underling dynamic systems [1]. The connection of Kolmogorov entropy and Lyapunov determinants of the system is defined by Pessin’s theorem. Sensitive dependence on initial conditions is a distinguishing feature of chaotic behavior. Initially, arbitrarily close points in the phase space produce significantly different trajectories. Characteristically, trajectories in chaotic systems diverge exponentially and Lyapunov exponents ( LLE ) were proved to be a good quantitative measure for the average rate of exponential divergence of two trajectories. Hence, positive LLE is indicative of unpredictable behavior. K is basically equal to the sum of the positive LLE of the system.
Theoretical defined entropy K is not able to be obtained on a base on finite, noisy data (signals) obtained from real measurements. The amount of data typically required to archive convergence is between 10 d to 30 d points to fill out a d-dimensional strange attractor for the chaotic system [2]. To overcome this, the approximate entropy, as a family of statistics, has been proposed [3]. It was shown that it can potentially distinguish low dimensional deterministic systems, periods and multiply periodic, high dimensional chaotic, stochastic and mixed systems [4]. The approximate entropy statistics is based on an Eckmann–Ruelle entropy formula [5], for the physical invariant measure for use with time series data.
This work presents the new quaternion, approximate entropy and its application to analysis of human gait kinematic data represented as a quaternion signal. This signal represents the segments orientations in time during the human gait. The research was aimed at the assessment of the influence of both walking speed and ground slope on the proposed quaternion entropy during treadmill walking. Gait data was obtained by the optical motion capture system. The analysis was carried out for femur, tibia and foot segments.

2. Materials and Methods

2.1. Background

Letting q 1 , q 2 , , q N be a time series signal of measurements equally spaced in time, the sequence of vectors in R m can be defined as x ( 1 ) , x ( 2 ) , , x ( N m + 1 ) , where x ( i ) = q i , q i + 1 , , q i + m 1 . Next, for each i and j, 1 i N m + 1 , and 1 j N m + 1 , the following function is defined as
C i m ( r , N ) = # ( d x ( i ) , x ( j ) r ) ( N m + 1 ) ,
where # is the number of elements in set, r is vector comparison distance and m is the dimension of created and compared vectors. The distance d is defined as follows:
d x ( i ) , x ( j ) = max k = 1 , 2 , , m q i + k 1 , q j + k 1 .
The definitions above are based on correlation dimension derivation and are needed to define Φ m ( r ) as
Φ m ( r , N ) = N m + 1 1 i = 1 N m + 1 log C i m ( r , N ) .
After fixing m and r, we obtain a definition of approximate entropy [3]:
ApEn ( m , r ) = lim N Φ m ( r , N ) Φ m + 1 ( r , N ) .
Given N samples, we can define the following formula as statistics:
ApEn ( m , r , N ) = Φ m ( r , N ) Φ m + 1 ( r , N ) .
Informally, given N signal’s samples, the family of statistics ApEn ( m , r , N ) is a measure of probability that two sequences that are similar for m points remain similar, that is, within a tolerance r, in the next sample in a signal. Thus, a low value of ApEn reflects a high degree of regularity. Following family statistics, ApEn has been widely used in clinical cardiovascular studies [6,7,8,9] and neurology, among others, for analysis of electroencephalogram signals [10,11,12,13].
Similar derivation has related to complexity measure sample entropy ( SampEn ) [6]. Sample entropy in counting the similar vectors does not count self-matches. The ApEn entropy algorithm counts each sequence as matching itself, which is a way to avoid occurrence of log ( 0 ) in calculations. In practice, this causes some bias value in results. This was a very detailed discussion in [6,14]. However, ApEn and SampEn also indicate more self-similarity in the time series.
Entropy to gait analysis was used only for times series consisting of spatio-temporal parameters like step time, length and width, stride intervals [14,15] or segment trajectories [16]. Analysis of center of COP (center of pressure) trajectory by sample entropy can be found in [17]. In [18], multiscale entropy [19] was applied on trunk acceleration data collected during a gait of subjects of different ages: toddlers, pre-scholar and scholar children, adolescents, young adults, adults and elderlies. Control entropy was used to analyse the acceleration data a i = a x , i , a y , i , a z , i as three independent signals [20]. Each cited entropy used to analyse gait data is based on ApEnt formula.
Based on actual knowledge, there is no proposition to compute entropy to analyse gait data based on unit quaternions’ time series. This conception allows for processing correlated data and to obtain results taking into account rotations in 3D. The same preliminary studies were done based on quaternion energy and entropy to the classification of people based on gait data [21].
In [22], the authors describe an example of application of a nonlinear time series analysis directed at identifying the presence of deterministic chaos in human gait kinematic data by means of the largest Lyapunov exponent ( LLE ). A positive LLE value is interpreted as an indicator of local instability. However, the LLE was computed based on separate Euler angles (three rotations angles about the axes of a coordinate system) or only angles of rotation in quaternion representation. Such analysis in quaternion form did not consider axis of rotation. In this research experiment, data from this same trial was used so the results will be compared. Figure 1a,c present the axis of rotation (imaginary part of quaternion) during 300 samples of gait (about 3 strides). Figure 1b,d show the angle of rotation (scalar part of quaternion). To analyse all aspects of movement during gait, all of the information should be processed. The new propositions of quaternion motion analysis tools are widely discussed in [23,24,25,26].
The o orientation indicates the orientation of the rigid body defined in the reference system (also named as the reference frame). Rotation means a change in orientation, o 1 o 2 . The group SO ( 3 ) (special orthogonal group) contains all rotations around the origin of the Euclidean coordinate system R 3 . Rotations can be described by orthogonal matrices. Rotations can also be written as a combination of three Euler angles around the three coordinate system axes. The angles are usually referred to as roll, pitch and yaw. In this representation, the order in which the rotation is performed is important. There is also the phenomenon of blocking one degree of freedom of movement (gimbal lock).
In computer graphics, the most common description of rotation is unit quaternion [27,28]. Quaternions [27,28] are an extension of complex numbers, q H where H is a quaternion algebra. Quaternion q = q 0 + i · q 1 + j · q 2 + k · q 3 , consists of real q 0 and imaginary part i 2 + j 2 + k 2 = i · j · k = 1 . Multiplication of two quaternions is marked as ⊗. We can write quaternion as vector and scalar part: q = q 0 , u , u = q 1 , q 2 , q 3 . Using the Euler rule for complex numbers, we have a rotation around the axis u by an angle α :
q = cos α 2 + u · sin α 2 .
In order to represent the rotation, unit quaternions are used, which fulfill the condition q , q = 1 , where . is a scalar product. Unit quaternions lie on the hypersphere (denoted as H 1 ) embedded in Euclidean space 4D. The H 1 space is a SO ( 3 ) map with double mapping. This means that each rotation (or orientation) can be represented by two unit quaternions, called antipode. They represent a rotation around this same axis, but with a positive or negative angle. Conjunction quaternion represents inverse rotation q * , which for unit quaternions is equal to inverse quaternion q 1 .
In summary, unit quaternions are a suitable mathematical tool for describing orientation and rotation in 3D space. It is possible to describe the rotation composition as multiplying quaternions. Conjunction quaternion represents inverse rotation. Such a description is free from the phenomenon of blocking one degree of freedom. The distance metrics between quaternions [29,30] have also been defined.
Using quaternions and building analysis in the domain H 1 , we can process rotational data as correlated data with simultaneous rotation information around each axis of rotation. It also eliminates the problem of blocking one degree of freedom (gimbal lock) and allows a consistent mathematical record.

2.2. Quaternion Approximate Entropy

Let us assume that input a motion signal consists of quaternions: q 1 , q 2 , , q n where q i H 1 and n = 2 k for some k N .
Furthermore, the signal is processed by the selective negation (hemispherization), that is, every quaternion q i i > 1 is converted to q i if q i , q i 1 < 0 due to duality of unit quaternions which represent rotations. It satisfies the requirement according to which two adjacent quaternions are located on the same hemisphere.
The new proposed metric, the quaternion, approximate entropy ApQuatEn , is based on Equation (5). Quaternions are of unit length, which means they are located only on a hypersphere H 1 . Thus, to compare rotations, it is sufficient to calculate cosine distances between related quaternions, which is reflected by angles between vectors formed by quaternions’ components. The scalar product q 1 , q 2 can be used to accomplish the task:
d cos ine q i , q j = 1 q i , q j 2 .
The distance between two quaternions q i + k 1 , q j + k 1 in Equation (2) is defined as d cos ine , so the equation is:
d x ( i ) , x ( j ) = max k = 1 , 2 , , m d cos ine q i + k 1 , q j + k 1 .

2.3. Treadmill Experiments

Time series were extracted from treadmill gait sequences which were recorded in the Human Motion Laboratory (HML) of the Polish-Japanese Academy of Information Technology (Bytom, Poland) by the optical motion capture system. Data was recorded with the use of the Vicon Nexsus optical motion capture system (Vicon Motion Systems Ltd., Oxford, UK). The skeleton model of 22 segments is applied and positions are traced based on 39 markers in standard full body Plug-In Gait marker placement. For analysis, only Euler angles’ orientations of femur, tibia, foot (left and right) segments were converted to unit quaternions time series signal.
Application of the AC5000M treadmill (SCIFIT Corporate Headquarters, United States) allowed recordings in three variants: at the preferred walking speed (PWS) of each subject (denoted as Normal speed), at 80% of the PWS (denoted as Slower speed) and at 120% of the PWS (denoted as Faster speed). Additionally, the recordings with the PWS on an inclined treadmill at slope of + 7 deg (denoted as slope Up) and 3 deg (denoted as slope Down). Three sequences of continuous walking of lengths of several dozen seconds were recorded with a frequency of 100 Hz at a given walking speed for every person. In this way, the total number of sequences recorded for every subject was equal to 15. At any time, participants could take a rest upon request. The mean duration of a sequence was equal to 71.9 s with standard deviation of 10.7. A sequence consisted in a mean of 57.8 strides with standard deviation of 5.3. Seventeen healthy subjects (6 women and 11 men) participated in the experiments. Mean values of their age 27.18 (standard deviation 8.72), height 174.47 cm (standard deviation 6.46), weight 71.14 (standard deviation 10.06) and preferred walking speed 2.38 m/s (standard deviation 0.63). A stride interval (i.e., the time elapsed between consecutive ipsilateral heel strikes) varied across subjects and variants of the walking speed and treadmill slope. Mean values of a stride interval for every variant are 1.20 s for Normal speed (standard deviation 0.15), 1.40 m/s for Slower speed (standard deviation 0.21), 1.13 m/s for Faster speed (standard deviation 0.14), 1.27 m/s for Slope-up (standard deviation 0.18) and 1.24 m/s for Slope-down (standard deviation 0.16). This same dataset was also used in [22,31].

3. Results and Discussion

For all time series, data of 17 participants in each configuration (Normal, Slower, Faster, Up, Down) repeated 3 times gives 255 data streams. For all data streams, the ApQuatEn was calculated. The following length of vectors m are used in calculations: 2 , 3 , 4 . In addition, the value of threshold r has to be defined. The proposition is to use the mean value of distance d cosine between each following q i and q i + 1 quaternion in time series. The r value was calculated for each sequence (each segment and each experiment configuration). All computations were performed using Matlab software (2016a, MathWorks Inc., Natick, MA, USA).
The results for ApQuatEn ( m = 2 , r = mean ( d cosine ) ) are presented using box plots comprising all three segments for left and right side, speed and slope (Figure 2, Figure 3 and Figure 4). The other parameters’ configurations give results presenting similar dependencies. The influence of parameters to result in quaternion approximate entropy is presented in Figure 5. Movements of all considered body parts and experiments configurations are characterized by the positive value of entropy, which detect and quantify chaotic behavior. This is consistent with results obtained with use of LLE in [22], for which all times series are characterized by positive values of LLE . The median values for all segments are also presented in Table 1. In addition, the aggregated median values for segments without taking into account the side of the body are included.
The values of entropy are the smallest for Slower speed for segments femur (median 0.286 ), tibia (median 0.233 ) and foot (median 0.410 ), which can indicate more regularity and predictability movement in those segments [16]. The highest values for different speed configurations are for all segments in Faster speed configuration. The values of entropy for different treadmill inclinations (Up and Down), where participants maintain normal speed, are somewhat
Higher than for configuration with Normal speed and flat treadmill. Differences of entropy values between configurations with treadmill inclinations for different segments are somewhat higher for Up than Down configuration. Big differences between entropy values for left and right tibia segments for Slower speed are interesting. Generally, values for femur segments (median value: Normal 0.337 , Faster 0.388 , Up 0.387 , Down 0.371 ) are smaller than for other segments, which can indicate more regularity in movement in this segment. The exception is Slower speed configuration where the smallest value is for tibia segments (median 0.233 ).
In [22], the two variants of LLE , based on quaternion angle, were estimated: the short-term LLE for the first stride and the long-term LLE over a fixed interval between the fourth and the tenth stride. The Pearson correlation coefficient for values of ApQuatEn ( m = 2 , r = mean ( d cosine ) ) for left and right femurs throughout all experiment configurations is 0.735 for the long-term LLE and 0.709 for the short-term LLE .
Table 2, Table 3 and Table 4 present the Pearson correlation coefficient of ApQuatEn ( m = 2 , r = mean ( d cosine ) ) in each segment and configuration. We can observe that values of proposed quaternion approximate entropy are characteristics for two groups of experiment configurations—treadmill velocity (Normal and Faster) and slope (Up and Down). All Pearson correlation coefficients in those two groups for all segments are greater than 0.5 . The Pearson correlation coefficients are greater for Slower and Up, Slower and Down than for Slower and Normal, Slower and Faster for femur and foot segments. This dependency is not so significant for tibia segments. The gaits with Normal and Faster speed have similar properties in the context of chaotic behaviors. The same can be observed for gaits on an inclined treadmill at slope Up and Down.
The detailed influence of lengths of vector m and value of threshold r into result ApQuatEn for left femur segments is presented in Figure 5. All values of entropy have similar dependencies within experiment configurations and segments. Generally, the ApQuatEn for m = 4 is smaller than for m = 2 . This means that differences between Φ 4 ( r ) and Φ 5 ( r ) are smaller than Φ 2 ( r ) and Φ 3 ( r ) . Typically, it is suggested that, for clinical data, m is to be set at 2 when utilizing the ApEn algorithm [7].
In [14], ApEn values were quantified on generated data from 100 to 10,000 data points, in increments of 100. As can be seen, the entropy values stabilize around N of 2000. In order to check the stabilization of entropy ApQuatEn values as N increased in experiment data, the chaotic logistic map was subjected to entropy analysis up to an N of 3700 data points. For this particular analysis, ( m = 2 , r = mean ( d cosine ) ) was chosen. Entropy was quantified on data in increments of 100. As can be seen in Figure 6, the entropy values stabilize around an N of 2700 time samples. For gait with Normal velocity, this is after a twenty-second stride.

4. Conclusions

This publication concerns the new quaternion approximate entropy ApQuatEn for quaternion time series represented orientations of human skeleton segments. The presented approach was used to quantify chaotic behavior in the context of regularity in time series. The example of application for treadmill walking data, in five variants of walking speed and treadmill slope, was discussed. Data was obtained from the optical motion capture system. It was confirmed that all considered time series are characterized by positive values of entropy values which quantify a local instability. The use of entropy as a mathematical tool to quantify predictability of gait parameters is an entirely emerging task in human movement research. This is the first proposition to consider quaternion signal and has analysed the correlated data which describe rotation in 3D. A possible future work can be related to develop quaternion version of negentropies to develop hypothesis testing for normality [32].

Funding

This publication was supported by the Rector’s grant in the field of scientific research and development works, Silesian University of Technology (Gliwice, Poland), Grant No. 02/020/RGJ/18/0125, (2018–2019).

Acknowledgments

The data was captured in the Human Motion Lab (HML) of the Polish-Japanese Academy of Information Technology in Bytom, Poland (http://bytom.pja.edu.pl/).

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Schuster, H.G.; Just, W. Deterministic Chaos: An Introduction; John Wiley & Sons: New York, NY, USA, 2006. [Google Scholar]
  2. Wolf, A.; Swift, J.B.; Swinney, H.L.; Vastano, J.A. Determining Lyapunov exponents from a time series. Phys. D Nonlinear Phenom. 1985, 16, 285–317. [Google Scholar] [CrossRef] [Green Version]
  3. Pincus, S. Approximate entropy (ApEn) as a complexity measure. Chaos Interdiscip. J. Nonlinear Sci. 1995, 5, 110–117. [Google Scholar] [CrossRef]
  4. Pincus, S.M. Approximate entropy as a measure of system complexity. Proc. Natl. Acad. Sci. USA 1991, 88, 2297–2301. [Google Scholar] [CrossRef] [PubMed]
  5. Eckmann, J.P.; Ruelle, D. Ergodic theory of chaos and strange attractors. In The Theory of Chaotic Attractors; Springer: Berlin/Heidelberg, Germany, 1985; pp. 273–312. [Google Scholar]
  6. Richman, J.S.; Moorman, J.R. Physiological time-series analysis using approximate entropy and sample entropy. Am. J. Physiol.-Heart Circ. Physiol. 2000, 278, H2039–H2049. [Google Scholar] [CrossRef] [PubMed]
  7. Pincus, S.M.; Viscarello, R.R. Approximate entropy: A regularity measure for fetal heart rate analysis. Obstet Gynecol 1992, 79, 249–255. [Google Scholar] [PubMed]
  8. Kaplan, D.; Furman, M.; Pincus, S.; Ryan, S.; Lipsitz, L.; Goldberger, A. Aging and the complexity of cardiovascular dynamics. Biophys. J. 1991, 59, 945. [Google Scholar] [CrossRef]
  9. Liu, C.; Liu, C.; Shao, P.; Li, L.; Sun, X.; Wang, X.; Liu, F. Comparison of different threshold values r for approximate entropy: Application to investigate the heart rate variability between heart failure and healthy control groups. Physiol. Meas. 2010, 32, 167. [Google Scholar] [CrossRef] [PubMed]
  10. Abásolo, D.; Hornero, R.; Espino, P.; Poza, J.; Sánchez, C.I.; de la Rosa, R. Analysis of regularity in the EEG background activity of Alzheimer’s disease patients with Approximate Entropy. Clin. Neurophysiol. 2005, 116, 1826–1834. [Google Scholar] [CrossRef] [PubMed]
  11. Zhang, X.S.; Roy, R.J.; Jensen, E.W. EEG complexity as a measure of depth of anesthesia for patients. IEEE Trans. Biomed. Eng. 2001, 48, 1424–1433. [Google Scholar] [CrossRef] [PubMed]
  12. Kannathal, N.; Choo, M.L.; Acharya, U.R.; Sadasivan, P. Entropies for detection of epilepsy in EEG. Comput. Methods Progr. Biomed. 2005, 80, 187–194. [Google Scholar] [CrossRef] [PubMed]
  13. Acharya, U.R.; Molinari, F.; Sree, S.V.; Chattopadhyay, S.; Ng, K.H.; Suri, J.S. Automated diagnosis of epileptic EEG using entropies. Biomed. Signal Process. Control 2012, 7, 401–408. [Google Scholar] [CrossRef]
  14. Yentes, J.M.; Hunt, N.; Schmid, K.K.; Kaipust, J.P.; McGrath, D.; Stergiou, N. The appropriate use of approximate entropy and sample entropy with short data sets. Ann. Biomed. Eng. 2013, 41, 349–365. [Google Scholar] [CrossRef] [PubMed]
  15. Costa, M.; Peng, C.K.; Goldberger, A.L.; Hausdorff, J.M. Multiscale entropy analysis of human gait dynamics. Phys. A Stat. Mech. Its Appl. 2003, 330, 53–60. [Google Scholar] [CrossRef]
  16. Ahmadi, S.; Sepehri, N.; Wu, C.; Szturm, T. Sample Entropy of Human Gait Center of Pressure Displacement: A Systematic Methodological Analysis. Entropy 2018, 20, 579. [Google Scholar] [CrossRef]
  17. Ramdani, S.; Seigle, B.; Lagarde, J.; Bouchara, F.; Bernard, P.L. On the use of sample entropy to analyze human postural sway data. Med. Eng. Phys. 2009, 31, 1023–1031. [Google Scholar] [CrossRef] [PubMed]
  18. Bisi, M.; Stagni, R. Complexity of human gait pattern at different ages assessed using multiscale entropy: From development to decline. Gait Posture 2016, 47, 37–42. [Google Scholar] [CrossRef] [PubMed]
  19. Costa, M.; Goldberger, A.L.; Peng, C.K. Multiscale entropy analysis of biological signals. Phys. Rev. E 2005, 71, 021906. [Google Scholar] [CrossRef] [PubMed]
  20. McGregor, S.J.; Busa, M.A.; Skufca, J.; Yaggie, J.A.; Bollt, E.M. Control entropy identifies differential changes in complexity of walking and running gait patterns with increasing speed in highly trained runners. Chaos Interdiscip. J. Nonlinear Sci. 2009, 19, 026109. [Google Scholar] [CrossRef]
  21. Szczęsna, A.; Świtoński, A.; Słupik, J.; Josiński, H.; Wojciechowski, K. Wavelet features in motion data classification. In Proceedings of the International Conference on Numerical Analysis and Applied Mathematics 2015 (ICNAAM-2015), Rhodes, Greece, 23–29 September 2015; AIP Publishing: Melville, NY, USA, 2016; Volume 1738. [Google Scholar]
  22. Piórek, M.; Josiński, H.; Michalczuk, A.; Świtoński, A.; Szczęsna, A. Quaternions and joint angles in an analysis of local stability of gait for different variants of walking speed and treadmill slope. Inf. Sci. 2017, 384, 263–280. [Google Scholar] [CrossRef]
  23. Szczęsna, A. Verification of the blobby quaternion model of human joint limits. Biomed. Signal Process. Control 2018, 39, 130–138. [Google Scholar] [CrossRef]
  24. Szczęsna, A. Quaternion wavelet-based energy and entropy to analysis human gait data. In Proceedings of the International Conference on Numerical Analysis and Applied Mathematics 2017 (ICNAAM-2017), Thessaloniki, Greece, 25–30 September 2017; AIP Publishing: Melville, NY, USA, 2018; Volume 1978. [Google Scholar]
  25. Szczęsna, A.; Świtoński, A.; Słupik, J.; Zghidi, H.; Josiński, H.; Wojciechowski, K. Quaternion Lifting Scheme Applied to the Classification of Motion Data. Inf. Sci. 2018. [Google Scholar] [CrossRef]
  26. Switonski, A.; Josinski, H.; Wojciechowski, K. Dynamic time warping in classification and selection of motion capture data. Multidimens. Syst. Signal Process. 2018. [Google Scholar] [CrossRef]
  27. Hanson, A.J. Visualizing Quaternions; Morgan Kaufmann: Burlington, MA, USA, 2005. [Google Scholar]
  28. Goldman, R. Understanding quaternions. Gr. Models 2011, 73, 21–49. [Google Scholar] [CrossRef]
  29. Huynh, D.Q. Metrics for 3D rotations: Comparison and analysis. J. Math. Imaging Vis. 2009, 35, 155–164. [Google Scholar] [CrossRef]
  30. Gramkow, C. On averaging rotations. J. Math. Imaging Vis. 2001, 15, 7–16. [Google Scholar] [CrossRef]
  31. Josiński, H.; Michalczuk, A.; Świtoński, A.; Mucha, R.; Wojciechowski, K. Quantifying chaotic behavior in treadmill walking. In Proceedings of the Asian Conference on Intelligent Information and Database Systems, Bali, Indonesia, 23–25 March 2015; pp. 317–326. [Google Scholar]
  32. Arellano-Valle, R.B.; Contreras-Reyes, J.E.; Stehlík, M. Generalized skew-normal negentropy and its application to fish condition factor time series. Entropy 2017, 19, 528. [Google Scholar] [CrossRef]
Figure 1. Axis and angle of rotation of femur and foot segments during gait (three strides).
Figure 1. Axis and angle of rotation of femur and foot segments during gait (three strides).
Entropy 21 00079 g001
Figure 2. Results values of ApQuatEn ( m = 2 , r = mean ( d cosine ) ) for left and right femur segments.
Figure 2. Results values of ApQuatEn ( m = 2 , r = mean ( d cosine ) ) for left and right femur segments.
Entropy 21 00079 g002
Figure 3. Results values of ApQuatEn ( m = 2 , r = mean ( d cosine ) ) for left and right tibia segments.
Figure 3. Results values of ApQuatEn ( m = 2 , r = mean ( d cosine ) ) for left and right tibia segments.
Entropy 21 00079 g003
Figure 4. Results values of ApQuatEn ( m = 2 , r = mean ( d cosine ) ) for left and right foot segments.
Figure 4. Results values of ApQuatEn ( m = 2 , r = mean ( d cosine ) ) for left and right foot segments.
Entropy 21 00079 g004
Figure 5. The value of entropy ApQuatEn in relation to the length of vector (m) and threshold distance r value for left femur segments.
Figure 5. The value of entropy ApQuatEn in relation to the length of vector (m) and threshold distance r value for left femur segments.
Entropy 21 00079 g005
Figure 6. The value of entropy for left femur segments (Normal speed) ApQuatEn ( m = 2 , r = mean ( d cosine ) ) in relation to data length N.
Figure 6. The value of entropy for left femur segments (Normal speed) ApQuatEn ( m = 2 , r = mean ( d cosine ) ) in relation to data length N.
Entropy 21 00079 g006
Table 1. Median values of ApQuatEn ( m = 2 , r = mean ( d cosine ) ).
Table 1. Median values of ApQuatEn ( m = 2 , r = mean ( d cosine ) ).
NormalFasterSlowerUpDown
lfemur0.3340.3710.2740.4050.344
rfemur0.3450.4060.3000.3710.396
femur0.3370.3880.2860.3870.371
ltibia0.2520.4460.1380.5990.353
rtibia0.5540.5190.5700.4630.396
tibia0.3380.5190.2330.5370.385
lfoot0.4780.5250.4960.4670.529
rfoot0.4200.5520.3560.4470.396
foot0.4430.5440.4100.4470.447
Table 2. The Pearson correlation coefficient of ApQuatEn ( m = 2 , r = mean ( d cosine ) ) calculated for left and right femur segments.
Table 2. The Pearson correlation coefficient of ApQuatEn ( m = 2 , r = mean ( d cosine ) ) calculated for left and right femur segments.
Left Femur Right Femur
NormalFasterSlowerUpDown NormalFasterSlowerUpDown
Normal1.0000.6990.3780.1570.334Normal1.0000.8620.4120.3470.343
Faster0.6991.0000.4620.3030.476Faster0.8621.0000.6990.5370.604
Slower0.3780.4621.0000.9210.941Slower0.4120.6991.0000.8990.951
Up0.1570.3030.9211.0000.944Up0.3470.5370.8991.0000.834
Down0.3340.4760.9400.9441.000Down0.3430.6040.9510.8341.000
Table 3. The Pearson correlation coefficient of ApQuatEn ( m = 2 , r = mean ( d cosine ) ) calculated for left and right tibia segments.
Table 3. The Pearson correlation coefficient of ApQuatEn ( m = 2 , r = mean ( d cosine ) ) calculated for left and right tibia segments.
Left Tibia Right Tibia
NormalFasterSlowerUpDown NormalFasterSlowerUpDown
Normal1.0000.6950.7330.1430.564Normal1.0000.6110.7180.8740.458
Faster0.6951.0000.6320.3280.549Faster0.6111.0000.5050.4960.439
Slower0.7330.6321.0000.4090.744Slower0.7180.5051.0000.7970.685
Up0.1430.3280.4091.0000.542Up0.8740.4960.7971.0000.646
Down0.5640.5490.7440.5421.000Down0.4580.4390.6850.6461.000
Table 4. The Pearson correlation coefficient of ApQuatEn ( m = 2 , r = mean ( d cosine ) ) calculated for left and right foot segments.
Table 4. The Pearson correlation coefficient of ApQuatEn ( m = 2 , r = mean ( d cosine ) ) calculated for left and right foot segments.
Left Foot Right Foot
NormalFasterSlowerUpDown NormalFasterSlowerUpDown
Normal1.0000.7280.5100.4390.491Normal1.0000.7210.4910.1350.390
Faster0.7281.0000.7780.5730.773Faster0.7211.0000.4490.2410.237
Slower0.5100.7781.0000.8370.923Slower0.4910.4491.0000.9010.909
Up0.4390.5730.8371.0000.751Up0.1350.2410.9011.0000.841
Down0.4910.7730.9230.7511.000Down0.3900.2370.9090.8411.000

Share and Cite

MDPI and ACS Style

Szczęsna, A. Quaternion Entropy for Analysis of Gait Data. Entropy 2019, 21, 79. https://doi.org/10.3390/e21010079

AMA Style

Szczęsna A. Quaternion Entropy for Analysis of Gait Data. Entropy. 2019; 21(1):79. https://doi.org/10.3390/e21010079

Chicago/Turabian Style

Szczęsna, Agnieszka. 2019. "Quaternion Entropy for Analysis of Gait Data" Entropy 21, no. 1: 79. https://doi.org/10.3390/e21010079

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

Article Metrics

Back to TopTop