Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
A Unified Framework for Street-View Panorama Stitching
Next Article in Special Issue
Estimation of Ground Reaction Forces and Moments During Gait Using Only Inertial Motion Capture
Previous Article in Journal
Ambulatory Assessment of Instantaneous Velocity during Walking Using Inertial Sensor Measurements
Previous Article in Special Issue
Estimation of Full-Body Poses Using Only Five Inertial Sensors: An Eager or Lazy Learning Approach?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Measurement and Geometric Modelling of Human Spine Posture for Medical Rehabilitation Purposes Using a Wearable Monitoring System Based on Inertial Sensors

by
Gheorghe-Daniel Voinea
*,
Silviu Butnariu
and
Gheorghe Mogan
Department of Automotive and Transport Engineering, Transilvania University of Brasov, 29 Eroilor Blvd, RO-500036 Brasov, Romania
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(1), 3; https://doi.org/10.3390/s17010003
Submission received: 27 September 2016 / Revised: 28 November 2016 / Accepted: 11 December 2016 / Published: 22 December 2016
(This article belongs to the Special Issue Body Worn Behavior Sensing)

Abstract

:
This paper presents a mathematical model that can be used to virtually reconstruct the posture of the human spine. By using orientation angles from a wearable monitoring system based on inertial sensors, the model calculates and represents the curvature of the spine. Several hypotheses are taken into consideration to increase the model precision. An estimation of the postures that can be calculated is also presented. A non-invasive solution to identify the human back shape can help reducing the time needed for medical rehabilitation sessions. Moreover, it prevents future problems caused by poor posture.

1. Introduction

Spinal disorders are usually caused by incorrect postures adopted during daily activities. Once they appear, medication is the first solution that helps to ease the pain. Then, a special rehabilitation plan, which consists of physical exercises and constant posture control, is undertaken with the purpose of eliminating the spinal disorder. In general, these physical therapies last for long periods and the patient has the tendency to go back the bad posture after a while.
The ongoing technological innovations in the area of sensing systems are creating new and better medical applications for various diseases or spinal disorders. Medical applications are not complete without a software or a method to process the data obtained from the hardware [1]. Physicians and patients need to receive meaningful parameters from the system. This data processing requires rather complex algorithms.
A short summary of the state of the art in spine monitoring systems and spine shape reconstruction is presented in the following paragraphs. The aim of this work is to develop a method to generate the shape of the spine, with data obtained from a wearable monitoring system.
There are several approaches in assessing the shape of the human spine, each with advantages and limitations. In general, the back shape can be reconstructed with the help of mathematical formulas from data obtained from various acquisition systems. The radiological method is commonly used and it implies the segmentation of the image and then measuring the angles [2,3,4]. Recently, 3D models were generated using DICOM images that are obtained from an X-ray computed tomography (or CT scans) [5,6,7,8,9,10,11,12].
Monitoring spine movement with optoelectronic systems can yield very accurate results using a relatively small set of surface markers. However, this type of measurements are constrained to special indoor facilities and usually do not offer information on the rotations of the spine [13,14]. Other methods used by researchers to reconstruct the shape of the spine are based on images that were obtained through invasive [15,16,17] or non-invasive techniques [18,19,20,21].
Portable solutions based on strain-gauge sensors that can be used without supervision have the advantage of monitoring the spinal motion in an unobtrusive manner [22,23]. Usually, an accelerometer is also needed to determine the orientation of the sensor while considering Earth’s gravitational field. The system developed in [24] has the capability of storing measured angular data, at up to 100 Hz. Measurements were taken while participants performed several flexion and extension exercises. To reconstruct the shape of the spine, a 6th order polynomial interpolation was applied. Then, the results were compared with the data from the Vicon [25] optical system, which is considered a gold standard. The authors demonstrated that the SpineDMS system can be used to assess the posture of patients with normal body mass index, but there are accuracy problems in the case of adipose patients. O’Sullivan [26] used the BodyGuard™ strain-gauge based device to monitor lumbo-pelvic sagittal plane movements and found it to be reliable, but not suitable for some postural tasks due to high errors. While most of the orientation devices express motion in degrees, the BodyGuard™ calculates the spinal flexion or extension according to the strain gauge elongation and gives a relative range of motion. The system is able to monitor posture or send real-time postural feedback.
High precision potentiometers are at the base of the CA 6000 Spine Motion Analyzer (CA SMA); a system used for studying the kinematics of the lumbar [27] and cervical [28] spine. However, the CA SMA system is not suitable for measuring the range of motion of the spine.
Another unobtrusive method to evaluate the alignment of the spine is presented in [29]. The aim of the study was to monitor the spine by measuring spring indentation present in a mattress and combine this manual data with an auto generated body shape model.
Inertial sensors have been used to measure spinal motion or the range of motion in many studies. One important aspect in measuring spinal motion with inertial sensors is making the data intuitive and user friendly for the clinicians and patients who will use the system. The data can be transformed into meaningful parameters such as rotation, flexion-extension and lateral bending [30,31]. Theobald [32] measured cervical range of motion with inertial sensors. It was proven that they are a viable and objective method for evaluating spine shapes.
Williams [33] investigated the lumbar spine curvature using fiber-optic sensors. These sensors are able to detect inclination according to the intensity of the transmitted light. A spline curve populated with positional data was used to provide the shape of the spine. The authors found this method appropriate for clinical or home use, with good results in measuring spinal curvature over a longer period. However, it is not a solution for measuring rotations or lateral bending due to the inflexibility of the fiber-optic wire.
A new 3-D mobile measurement system that can assess trunk inclination and spine curvature was developed and evaluated in [34]. The system combines an ultrasound-based device with an inclinometer, which also comes with a personalized data analyzing software.
The novelty of this study consists in the development of a mathematical model of the human spine that can be used in medical rehabilitation sessions. The model can be implemented in a wearable spine monitoring system. It relies on the orientation data from inertial sensors and a known distance between the sensors. A classification of the identified models that correspond to possible spinal disorders is also presented. Thus, by considering the angles measured in different areas of the spine, the mathematical model generates the associated curves of the spine shape. The validation of the model was made offline, without using human subjects. The main advantage for the offline measurement is the predictability of how the system should work, thus allowing us to efficiently evaluate the performance of our system through a statistical analysis.

2. Materials and Methods

2.1. Measurement System

The concept of the spine monitoring system is presented in Figure 1a. It has the following components: a shirt on which the sensors able to detect movement are placed, a controller to collect data from the sensors and transmit it to a smartphone for the user to visualize his/her posture.
The first prototype uses five sensors (Figure 1b) that are mounted at equal distances, a condition that is needed and sufficient for reconstructing the spine shape with the developed mathematical model. The optimal number of sensors (five) was determined based on the following reasons: the five sensors are distributed evenly on the entire spine in order to detect the curvatures (less than five yield an inaccurate model), the post-processing is not very complex and the system is low cost. The mathematical model calculates and offers meaningful data that can be easily understood by physicians or patients. A spine monitoring system has been developed for acquiring data that describes the posture of the spine. It is based on inertial sensors that are positioned on a flexible frame, which is then fitted on the patient with flexible straps and fabric fasteners. It can communicate wirelessly with a portable device or send the data on a cloud server using GPRS connection. The architecture of the spine monitoring system is presented in Figure 2.
The data that represents the movement of the spine is acquired by the inertial sensors with a development board, stored on a microSD card. The data can be sent periodically to a portable device (smartphone, tablet) or straight to the online server for storage and processing. The Web based interface is used by patients or physicians to access the data, which was already processed and transformed into relevant information (Figure 3).
Inertial sensors are commonly used for monitoring movement. They are also used in smartphones or other small devices, such as fitness bracelets.

2.2. Mathematical Model of the Spine

The spine acts as a pillar to support the human body and it protects the spinal cord. When viewed from the side it normally has an S-shape due to three natural curves, which is why the spine can be divided into three regions: cervical, thoracic and lumbar spine.
The known data for the mathematical model are five orientation angles (from the inertial sensors) and four distances (between the sensors). We studied three aproaches:
  • A mathematical model based on polynomial functions: a 7th grade polynomial is required that is not time and computational efficient.
  • A mathematical model based on spline functions: cubic functions for each segment. In this case, we ended with a nonlinear sistem of equations that we found to be difficult to use.
  • A mathematical model based on circle arcs (Figure 4): the equations needed to calculate the radius and the coordinates of the imu sensors are only first grade relations. Thus, we get to a linear system of equations that is easier and faster to resolve.
We consider the following hypothesis:
  • The wearable monitoring system provides orientation in 3-axis, but we are only using one axis (using all the axes is the subject of another research article).
  • The proposed mathematical model is based on the approximation of the real curves with circle arcs and geometrical functions.
  • The inputs of the model are the distances measured between the sensors on the flexible frame. The orientation angles are provided by the inertial sensors.
  • The segments formed between two sensors are approximated with arcs of length l1, l2, l3, l4 and radius R1, R2, R3, R4. The circle arcs can have different radius, which in the common point have the same tangent to the slope that is equal with the measured value from the inertial sensors.
  • The circle arcs will be tangent to each other and a smooth curve will be obtained.
Based on the inflection point theory from [35,36], 14 possible models were proposed. These models were classified according to two criteria: the number of maximum points, determined by the coordinates of the point which is furthest from the reference, and the relations between the angles of each inertial sensor. The 14 models are thus grouped into four categories: A—single max point positive curve, B—single max point negative curve, C—double max points positive curve, D—double max negative points curve.
The model can calculate with a very good precision the (X, Y) coordinates of the position of the sensors and reconstruct the natural curves of the spine.
The mathematical model has the following steps (for C-shape), which have been used for all the modelled use cases:
  • Determine the length of the circle arcs generated between two successive sensors: l1, l2, l3 and l4 based on radius (R1, R2, R3, R4) and measured angles (p1, p2, p3, p4, p5);
  • Calculate the distances on X and Y axis for each sensor, except for the first, which is considered a reference: (d1x, d1y), (d2x, d2y), (d3x, d3y), (d4x, d4y), (d5x, d5y);
  • Solve the linear system of 14 equations (Equations (1)–(7) from Table 1) with 14 unknowns (R1, R2, R3, R4, d1x, d1y, d2x, d2y, d3x, d3y, d4x, d4y, d5x, d5y);
  • Calculate the (X, Y) coordinates for each sensor (x2, y2), (x3, y3), (x4, y4) and (x5, y5) and for maximum point (Xmax, Ymax);
  • Calculate the coordinates of the circle center corresponding to the four arcs (see Table 1).
The algorithm for the S-shape is similar with the C-shape algorithm with the difference that it has one more maximum point and the number of equations increases. In the following paragraphs, the equations used for the two scenarios (see Table 2) are presented: a single maximum point curve (a C-shaped curve, which characterizes a kyphosis) and a double maximum point curve (an S-shaped curve, which represents a healthy/normal spine), while the rest are classified according to a set of conditions and presented in Table A1 from the Appendix A. The five sensors of the spine monitoring system create four arcs that will be analytically determined. The distances on the X and Y-axis are calculated if the first sensor is considered as a reference, with the following equations:
The equations above form a system of linear equations that can be solved automatically with a software such as MATLAB or Maple. Equation (1) represent the formulas for the length of an arc. The following six are used to calculate the distances on the X and Y axis of the sensors. Equations (8)–(13) describe the coordinates of the sensors and the maximum points of the curves. The last step is to calculate the coordinates of the center of the circle, for each arc. In order for the prototype to work with the mathematical model the following must be taken into consideration:
  • The distance between the inertial sensors must be equal (and known), depending on the length of the spine;
  • The position of the first sensor to be at the C7 vertebrae, which in general can be easily palpated;
  • The prototype must be in contact with the skin. In the lumbar area of the spine, where there is a more pronounced curve, a flexible strap with a soft material can be used to push the sensing array closer to the spine.

2.3. Simulation of Mathematical Model

In order to test the mathematical model for the shape of the spine, the equations presented in Section 2.2 were implemented in Maple [37], a mathematical and analytical software. The models were simulated using possible angles that follow the classification criteria and a distance between sensors of 150 mm. This value was chosen because the average length of the male human spine is about 710 mm and the five sensors are equidistant.

2.4. Hardware Equipment Analysis

The main components of the spine shape monitoring system that were tested and presented are the acquisition and control board and the inertial measurement units.

2.4.1. Testing the Acquisition Board

We used the following methodology: identify three acquisition boards that could be used in the spine monitoring system, develop a program to collect data from inertial sensors and analyze the results, taking in consideration the ease of use, technical performance and size.
We reached the following conclusions: the Intel Edison board has an advantage regarding the technical performances, Bluetooth capability and microSD adapter, but the Teensy and Arduino Due are better in terms of price, size, ease of use and energy efficiency.

2.4.2. Testing the Sensors

A special device was used to test the inertial sensors. One sensor can be tested at a time, by precisely positioning it at a certain predetermined angle (Figure 5a). In order to simultaneously test a larger number of sensors we used a prototyping board and an I2C multiplexer. The test equipment needed to test five inertial sensors is presented in Figure 5b.
Using the same methodology as for the acquisition boards, we identified five different inertial measurement units: Bosch BNO055, MPU 9150, Flora LSM9DS0, AltIMU 10 v.4, MinIMU 9 V3. A comparative analysis is presented in Table 3.
After an analysis of the five IMUs regarding datasheet characteristics, ease of use and data reliability we chose the Bosch BNO055 sensors. These sensors can provide absolute orientation due to the built in sensor fusion algorithms and auto calibration function. The high speed ARM Cortex-M0 processor processes the raw data from the accelerometer, gyroscope and magnetometer and outputs data in the following usable formats: quaternions, Euler angles or vectors. Regarding the supply voltage and supported communication protocols, all tested sensors have similar results.
The next step was to simultaneously test five BNO055 sensors that were connected to a Teensy 3.1 development board, through an I2C multiplexer (Figure 5b).
The Bosch BNO055 (Sety Robotics, Bucharest, Romania) 9-axis absolute orientation sensor is used for the spine monitoring system due to its small size of 20 mm × 27 mm × 4 mm, auto-calibration function, real-time measurement capability and because of the efficient fusion algorithm that is already implemented. Although the sensor is capable of sensing orientation on 3-axis, in this study we focus only on one axis.
Before the spine monitoring system is ready to be used, it needs to be calibrated. There are two steps in calibrating the system: the first step is the IMU calibration: the BNO055 have an auto-calibration function, this requires moving it for a few seconds, while making a circular “∞” pattern. The second step is the patient calibration: the patient is instructed to stay in a specific position and the data captured by the system will serve as a reference. Every update received from the sensors will be compared with this reference values. The patient calibration is a subject for a future research.

2.4.3. Testing the Wireless Communication

We considered two methods of data transmission from the acquisition board: with a Bluetooth adapter that will send the data to a mobile device or with a GSM-GPRS module that will upload the data on the server. The two solutions were tested using the Sparkfun Bluetooth Mate Silver [38] and A-GSM 2.064 module [39].The central component of the A-GSM 2.064 module is the Quectel modem that is capable of connecting to the Internet using TCP, UDP, HTTP or FTP. We chose FTP to upload data to a server using AT commands. Use of a GSM-GPRS module has the advantage of directly connecting the wearable monitoring system to the Internet, but that comes with a smaller autonomy due to the energy that the module needs in order to use GSM or GPRS.

2.4.4. Testing the Local Data Storage

A wearable monitoring system must be capable of storing data throughout the day. The solution is to use a microSD card adapter because of the small size and ease of use. Using Serial Peripheral Interface we connected the card reader to the development board and tested saving the orientation data in a text file, with the values separated with a comma for easier post processing.

2.5. Measurement Protocol

The testing stand is presented in Figure 6. Measurements are taken in static positions and several use cases were studied and modeled. Calibration is repeated every time the system is disconnected from the power source, although calibration data can be saved and used for future measurements.
In Figure 7, the flow of the measurement protocol is presented. A theoretical spine shape is plotted and used on a testing table, and then inertial sensors capture reference values. Lastly, the spine monitoring system is positioned in order to correspond with the drawing. The data from the sensors are saved and are used by the mathematical model to calculate the coordinates of the sensors, using a Maple application. All use cases have been tested using this method.

2.6. Testing Methodology and Results

Two types of posture were chosen for further testing, specifically A4 and C2 (Figure 8). Ten measurements were taken for each shape model (a total of 20 tests), based on the procedure presented in Figure 7. Thus, the flexible spine monitoring system was positioned on top of the plotted shape on a paper and the orientation data was saved. After reading the data from the five sensors, the system was removed and then repositioned after verifying the IMUs calibration status. The repeatability of measurements was taken into consideration, as well as the precision. The saved data from the C2 and A4 postures was then processed using the XLSTAT software [40].

2.6.1. C2 Posture

The quantitative data that resulted is presented in Table 4.
Before analyzing the data using advanced statistical methods, it is customary to check the box chart type (or whisker diagram), which is a simple and complete representation of the obtained results. Thus, we can identify the data trends, anomalies locations and also visualize the minimum, maximum and medium data values. The box plots for the C2 posture are presented in Figure 9.
In statistics, dispersion is a method used to describe how spread out a set of data is. The most common examples of statistical dispersion are standard deviation, variance and interquartile range. Variance represents the arithmetic mean of squared deviations of individual values of a statistical set of experiences, from the arithmetic mean of the whole set. Standard deviation is used to measure the amount of variation or dispersion from the arithmetic mean of a set of data. It is obtained by calculating the square root of variance and it is the most commonly used indicator to characterize a data set, after the average indicator.

2.6.2. A4 Posture

The quantitative data that resulted is presented in Table 5.
The box plots for A4 posture are presented in Figure 10.

2.7. Visual Representation of Obtained Measurements

After an in depth analysis of the measurement data and the statistical results, we determined that there is a cumulative error in the mathematical model that we eliminated by introducing compensation values in the equations (the equations are solved in cascade; due to this fact, the error propagates from the first equation to the last with an increasing effect).
A graphical method was chosen in order to visualize the results of the measurements (Figure 11). Thus, the 10 results for each individual shape model were used with the Maple application we developed. It outputs the coordinates of points where the inertial sensors were positioned, the coordinates of the circle center and corresponding radius. The resulted data was then used with AUTOCAD to draw 10 curves (black color), that are compared with a posture (magenta color) obtained using the average values from the inertial sensors. The resulted shapes were compared with the initial curve plotted on the paper (blue color).
The graphical results of the proposed mathematical model was compared using the following method: the coordinates of the points P1 to P5, calculated by the mathematical model using an average of the sensor data, were graphically represented in AUTOCAD; two types of curves were obtained, one using the SPLINE command and the second is based on the arc circle coordinates. The results are presented in Figure 12 for posture C2 and for posture A4 in Figure 13.

Error assessment

The procedure for assessing the displacement error of the mathematical model is the following:
  • The spine monitoring system is positioned in a C and S-shape corresponding to models A4 and C2;
  • The orientation angles from the inertial sensors are given to the simulation software in which the mathematical model is implemented;
  • The mathematical model calculates the coordinates of the sensors;
  • The real coordinates of the sensors are manually measured on the testing table;
  • We analyze the results from the simulation software with the manually measured distances.
The cumulative error from the mathematical model was reduced to less than 5 mm and the reproduced spine shape has a high degree of fidelity (see Figure 14).
Considering that the maximum error percentage is less than 5%, we can conclude that the developed mathematical model can successfully reproduce the spine curvatures and is suitable for use in medical rehabilitation applications or postural monitoring.

3. Conclusions

A mathematical model that can be used to reproduce the shape of the spine with data from a wearable monitoring system was presented. The aim of the model is to be used in medical rehabilitation sessions and to help patients to avoid incorrect postures. The simulations presented very good results and demonstrated that the mathematical model can be used to estimate the coordinates of the sensors, which leads to a realistic reconstruction of the spine shape using inertial sensors.
The mathematical model was tested using real orientation angles of the spine monitoring system described in the previous sections. The repeatability of measurements is a key component of precision in any measurement system. The test-retest reliability has been assessed by using two scenarios, one using a C-shape model and one using an S-shape model.
The objective data from the inertial sensors throughout the day also helps physicians to give personalized physical exercises. Based on the inflections point of the spine, 14 possible models were identified and simulated.
The current mathematical model only uses one axis from the inertial sensors, but we plan to expand the model in order to reconstruct movement in all the directions. This advanced model will be used for spinal disorders, such as kyphosis, lordosis and scoliosis. Advances in technology have a direct impact in medical applications and together with software innovations could help to create new and better solutions for various human diseases.

Acknowledgments

This paper was realized within the Partnership Programme in priority domains—PN-II, which runs with the financial support of MEN-UEFISCDI, Project No. 227/2014, System for Diagnosis and Therapy of Spine Diseases (SPINE).

Author Contributions

Gheorghe-Daniel Voinea carried out the main research. Silviu Butnariu and Gheorghe Mogan contributed equally in the design of the experiments, experimental work, data analysis and manuscript writing.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Estimation of the postures that can be reproduced by the mathematical model.
Table A1. Estimation of the postures that can be reproduced by the mathematical model.
NameCriteriaResults
A1Max point between P2 and P3;
p1 > p2,
p2 < p3,
p3 < p4,
p4 < p5.
Sensors 17 00003 i001
A2Max point between P3 and P4;
p1 > p2,
p2 > p3,
p3 < p4,
p4 < p5.
Sensors 17 00003 i002
A3Max point between P1 and P2;
p1 > p2,
p2 < p3,
p3 < p4,
p4 < p5.
Sensors 17 00003 i003
A4Max point between P4 and P5;
p1 > p2,
p2 > p3,
p3 > p4,
p4 < p5.
Sensors 17 00003 i004
B1Max point between P2 and P3;
p1 > p2,
p2 < p3,
p3 < p4,
p4 < p5.
Sensors 17 00003 i005
B2Max point between P3 and P4;
p1 > p2,
p2 > p3,
p3 < p4,
p4 < p5.
Sensors 17 00003 i006
B3Max point between P1 and P2;
p1 > p2,
p2 < p3,
p3 < p4,
p4 < p5.
Sensors 17 00003 i007
B4Max point between P4 and P5;
p1 > p2,
p2 > p3,
p3 > p4,
p4 < p5.
Sensors 17 00003 i008
C1Max point between P1-P2 and P3-P4;
p1 > p2,
p2 < p3,
p3 > p4,
p4 < p5.
Sensors 17 00003 i009
C2Max point between P2-P3 and P4-P5;
p1 > p2,
p2 < p3,
p3 > p4,
p4 < p5.
Sensors 17 00003 i010
C3Max point between P1-P2 and P3-P4;
p1 > p2,
p2 < p3,
p3 > p4,
p4 < p5.
Sensors 17 00003 i011
D1Max point between P1-P2 and P3-P4;
p1 > p2,
p2 < p3,
p3 > p4,
p4 < p5.
Sensors 17 00003 i012
D2Max point between P2-P3 and P4-P5;
p1 > p2,
p2 < p3,
p3 > p4,
p4 < p5.
Sensors 17 00003 i013
D3Max point between P1-P2 and P3-P4;
p1 > p2,
p2 < p3,
p3 > p4,
p4 < p5.
Sensors 17 00003 i014

References

  1. Devedžić, G.; Ćuković, S.; Luković, V.; Milošević, D.; Subburaj, K.; Luković, T. ScolioMedIS: Web-oriented information system for idiopathic scoliosis visualization and monitoring. Comput. Meth. Program. Biomed. 2012, 108, 736–749. [Google Scholar] [CrossRef] [PubMed]
  2. Humbert, L.; De Guise, J.A.; Aubert, B.; Godbout, B.; Parent, S.; Mitton, D.; Skalli, W. 3D Reconstruction of the Spine from Biplanar X-rays Using Longitudinal and Transversal Inferences. J. Biomech. 2007, 40, S160. [Google Scholar] [CrossRef]
  3. Lecron, F.; Benjelloun, M.; Mahmoudi, S. Cervical spine mobility analysis on radiographs: A fully automatic approach. Comput. Med. Imaging Grap. 2012, 36, 634–642. [Google Scholar] [CrossRef] [PubMed]
  4. Milenkovic, S.M.; Kocijancic, R.I.; Belojevic, G.A. Left handedness and spine deformities in early adolescence. Eur. J. Epidemiol. 2004, 19, 969–972. [Google Scholar] [CrossRef] [PubMed]
  5. Vrtovec, T.; Likar, B.; Pernus, F. Quantitative analysis of spinal curvature in 3D: Application to CT images of normal spine. Phys. Med. Biol. 2008, 53, 1895–1908. [Google Scholar] [CrossRef] [PubMed]
  6. Salem, W.; Lenders, C.; Mathieu, J.; Hermanus, N.; Klein, P. In vivo three-dimensional kinematics of the cervical spine during maximal axial rotation. Man. Ther. 2013, 18, 339–344. [Google Scholar] [CrossRef] [PubMed]
  7. Mannion, A.F.; Knecht, K.; Balaban, G.; Dvorak, J.; Grob, D. A new skin-surface device for measuring the curvature and global and segmental ranges of motion of the spine: Reliability of measurements and comparison with data reviewed from the literature. Eur. Spine J. 2004, 13, 122–136. [Google Scholar] [CrossRef] [PubMed]
  8. Plamondon, A.; Delisle, A.; Larue, C.; Brouillette, D.; Mcfadden, D.; Desjardins, P.; Larivière, C. Evaluation of a hybrid system for three-dimensional measurement of trunk posture in motion. Appl. Ergon. 2007, 38, 697–712. [Google Scholar] [CrossRef] [PubMed]
  9. Kadoury, S.; Labelle, H.; Paragios, N. Automatic inference of articulated spine models in CT images using high-order Markov Random Fields. Med. Image Anal. 2011, 15, 426–437. [Google Scholar] [CrossRef] [PubMed]
  10. Laurent, C.P.; Jolivet, E.; Hodel, J.; Decq, P.; Skalli, W. New method for 3D reconstruction of the human cranial vault from CT-scan data. Med. Eng. Phys. 2011, 33, 1270–1275. [Google Scholar] [CrossRef] [PubMed]
  11. Reutlinger, C.; Gedet, P.; Buchler, P.; Kowal, J.; Rudolph, T.; Burger, J.; Scheffler, K.; Hasler, C. Combining 3D tracking and surgical instrumentation to determine the stiffness of spinal motion segments: A validation study. Med. Eng. Phys. 2011, 33, 340–346. [Google Scholar] [CrossRef] [PubMed]
  12. Yamauchi, T.; Yamazaki, M.; Okawa, A.; Furuya, T.; Hayashi, K.; Sakuma, T.; Takahashi, H.; Yanagawa, N.; Koda, M. Efficacy and reliability of highly functional open source DICOM software (OsiriX) in spine surgery. J. Clin. Neurosci. 2010, 17, 756–759. [Google Scholar] [CrossRef] [PubMed]
  13. Ranavolo, A.; Don, R.; Draicchio, F.; Bartolo, M.; Serrao, M.; Padua, L.; Cipolla, G.; Pierelli, F.; Iavicoli, S.; Sandrini, G. Modelling the spine as a deformable body: Feasibility of reconstruction using an optoelectronic system. Appl. Ergon. 2013, 44, 192–199. [Google Scholar] [CrossRef] [PubMed]
  14. Schmid, S.; Studer, D.; Hasler, C.C.; Romkes, J.; Taylor, W.R.; Lorenzetti, S.; Brunner, R. Quantifying spinal gait kinematics using an enhanced optical motion capture approach in adolescent idiopathic scoliosis. Gait Posture 2016, 44, 231–237. [Google Scholar] [CrossRef] [PubMed]
  15. Berthonnaud, E.; Dimnet, J.; Roussouly, P.; Labelle, H. Analysis of the sagittal balance of the spine and pelvis using shape and orientation parameters. J. Spinal Disord. Tech. 2005, 18, 40–47. [Google Scholar] [CrossRef] [PubMed]
  16. Meakin, J.R.; Gregory, J.S.; Smith, F.W.; Gilbert, F.J.; Aspden, R.M. Characterizing the shape of the lumbar spine using an active shape model: Reliability and precision of the method. Spine 2008, 33, 807–813. [Google Scholar] [CrossRef] [PubMed]
  17. Tillotson, K.M.; Burton, A.K. Noninvasive measurement of lumbar sagittal mobility. An assessment of the flexicurve technique. Spine 1991, 16, 29–33. [Google Scholar] [CrossRef] [PubMed]
  18. Remondino, F. 3-D reconstruction of static human body shape from image sequence. Comput. Vis. Image Underst. 2004, 93, 65–85. [Google Scholar] [CrossRef]
  19. Jia, W.; Yi, W.J.; Saniie, J.; Oruklu, E. 3D image reconstruction and human body tracking using stereo vision and Kinect technology. In Proceedings of the 2012 IEEE International Conference on Electro/Information Technology (EIT), Indianapolis, IN, USA, 6–8 May 2012; pp. 1–4.
  20. Li, Z.; Jia, W.; Mao, Z.H.; Li, J.; Chen, H.C.; Zuo, W.; Wang, K.; Sun, M. Anthropometric body measurements based on multi-view stereo image reconstruction. In Proceedings of the 2013 35th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Osaka, Japan, 3–7 July 2013; pp. 366–369.
  21. Yang, H.D.; Lee, S.W. Reconstruction of 3D human body pose from stereo image sequences based on top-down learning. Pattern Recognit. 2007, 40, 3120–3131. [Google Scholar] [CrossRef]
  22. Consmüller, T.; Rohlmann, A.; Weinland, D.; Druschel, C.; Duda, G.N. Comparative evaluation of a novel measurement tool to assess lumbar spine posture and range of motion. Eur. Spine J. 2012, 21, 2170–2180. [Google Scholar] [CrossRef] [PubMed]
  23. Dreischarf, M.; Pries, E.; Bashkuev, M.; Putzier, M.; Schmidt, H. Differences between clinical “snap-shot” and “real-life” assessments of lumbar spine alignment and motion – What is the “real” lumbar lordosis of a human being? J. Biomech. 2016, 49, 638–644. [Google Scholar] [CrossRef] [PubMed]
  24. Taylor, W.R.; Consmüller, T.; Rohlmann, A. A novel system for the dynamic assessment of back shape. Med. Eng. Phys. 2010, 32, 1080–1083. [Google Scholar] [CrossRef] [PubMed]
  25. Motion Capture Systems Vicon. Available online: http://www.vicon.com/ (accessed on 18 August 2016).
  26. O’Sullivan, K.; O’Sullivan, L.; Campbell, A.; O’Sullivan, P.; Dankaerts, W. Towards monitoring lumbo-pelvic posture in real-life situations: Concurrent validity of a novel posture monitor and a traditional laboratory-based motion analysis system. Man. Ther. 2012, 17, 77–83. [Google Scholar] [CrossRef] [PubMed]
  27. Mieritz, R.M.; Bronfort, G.; Jakobsen, M.D.; Aagaard, P.; Hartvigsen, J. Reliability and measurement error of sagittal spinal motion parameters in 220 chronic low back pain patients using a 3D measurement device. Spine J. 2014, 14, 1835–1843. [Google Scholar] [CrossRef] [PubMed]
  28. Petersen, C.M.; Schuit, D.; Johnson, R.D.; Knecht, H.; Levine, P. Agreement of measures obtained radiographically and by the OSI CA-6000 Spine Motion Analyzer for cervical spinal motion. Man. Ther. 2008, 13, 200–205. [Google Scholar] [CrossRef] [PubMed]
  29. Verhaert, V.; Druyts, H.; Van Deun, D.; Exadaktylos, V.; Verbraecken, J. Estimating spine shape in lateral sleep positions using silhouette-derived body shape models. Int. J. Ind. Ergon. 2012, 42, 489–498. [Google Scholar] [CrossRef]
  30. Goodvin, C.; Park, E.J.; Huang, K.; Sakaki, K. Development of a real-time three-dimensional spinal motion measurement system for clinical practice. Med. Biol. Eng. Comput. 2006, 44, 1061–1075. [Google Scholar] [CrossRef] [PubMed]
  31. Lee, J.K.; Desmoulin, G.T.; Khan, A.H.; Park, E.J. Comparison of 3D spinal motions during stair-climbing between individuals with and without low back pain. Gait Posture 2011, 34, 222–226. [Google Scholar] [CrossRef] [PubMed]
  32. Theobald, P.S.; Jones, M.D.; Williams, J.M. Do inertial sensors represent a viable method to reliably measure cervical spine range of motion? Man. Ther. 2011, 17, 92–96. [Google Scholar] [CrossRef] [PubMed]
  33. Williams, J.M.; Haq, I.; Lee, R.Y. Dynamic measurement of lumbar curvature using fibre-optic sensors. Med. Eng. Phys. 2010, 32, 1043–1049. [Google Scholar] [CrossRef] [PubMed]
  34. Wunderlich, M.; Rüther, T.; Essfeld, D.; Erren, T.C.; Piekarski, C.; Leyk, D. A new approach to assess movements and isometric postures of spine and trunk at the workplace. Eur. Spine J. 2011, 20, 1393–1402. [Google Scholar] [CrossRef] [PubMed]
  35. Roussouly, P.; Nnadi, C. Sagittal Plane Deformity: An Overview of Interpretation and Management. Eur. Spine J. 2010, 19, 1824–1836. [Google Scholar] [CrossRef] [PubMed]
  36. Devedžić, G.; Stojanović, R.; Luković, V.; Ćuković, S.; Milošević, D. Identification of anatomical landmarks for intelligent postural sensing. In Proceedings of the 2012 Mediterranean Conference on Embedded Computing (MECO), Bar, Montenegro, 19–21 June 2012; pp. 70–73.
  37. Maplesoft Technical Computing Software. Available online: http://www.maplesoft.com/ (accessed on 18 June 2016).
  38. SparkFun Bluetooth Mate Silver. Available online: https://www.sparkfun.com/products/12576 (accessed on 10 April 2016).
  39. Arduino & Rasberry PI GSM/GPRS/SMS/DTMF Shield. Available online: http://itbrainpower.net/a-gsm (accessed on 5 April 2016).
  40. Statistical Software for Excel. Available online: https://www.xlstat.com/en/ (accessed on 23 June 2016).
Figure 1. (a) Wearable spine monitoring system concept; (b) Wearable spine monitoring with flexible straps.
Figure 1. (a) Wearable spine monitoring system concept; (b) Wearable spine monitoring with flexible straps.
Sensors 17 00003 g001
Figure 2. The architecture of the spine monitoring system.
Figure 2. The architecture of the spine monitoring system.
Sensors 17 00003 g002
Figure 3. Data flow diagram of the spine monitoring system.
Figure 3. Data flow diagram of the spine monitoring system.
Sensors 17 00003 g003
Figure 4. Geometric model: (a) C-shape, (b) S-shape.
Figure 4. Geometric model: (a) C-shape, (b) S-shape.
Sensors 17 00003 g004
Figure 5. (a) Device for testing IMUs; (b) Testing equipment for simultaneously testing five IMUs.
Figure 5. (a) Device for testing IMUs; (b) Testing equipment for simultaneously testing five IMUs.
Sensors 17 00003 g005
Figure 6. Testing stand.
Figure 6. Testing stand.
Sensors 17 00003 g006
Figure 7. Procedure for testing the mathematical model: (a) A theoretical spine shape; (b) Getting the reference position; (c) The system is fitted to correspond to the plotted shape; (d) The inertial sensors collect orientation data.
Figure 7. Procedure for testing the mathematical model: (a) A theoretical spine shape; (b) Getting the reference position; (c) The system is fitted to correspond to the plotted shape; (d) The inertial sensors collect orientation data.
Sensors 17 00003 g007
Figure 8. (a) C2 posture and (b) A4 posture used for measurements.
Figure 8. (a) C2 posture and (b) A4 posture used for measurements.
Sensors 17 00003 g008
Figure 9. Box plots—C2.
Figure 9. Box plots—C2.
Sensors 17 00003 g009
Figure 10. Box plots—A4.
Figure 10. Box plots—A4.
Sensors 17 00003 g010
Figure 11. Visual representation of obtained results.
Figure 11. Visual representation of obtained results.
Sensors 17 00003 g011
Figure 12. (a) Plotted curve; (b) Reconstructed curve, using SPLINE command from AutoCAD; (c) Reconstructed curve using circle arcs; (d) Superimposed curves from (b) and (c).
Figure 12. (a) Plotted curve; (b) Reconstructed curve, using SPLINE command from AutoCAD; (c) Reconstructed curve using circle arcs; (d) Superimposed curves from (b) and (c).
Sensors 17 00003 g012
Figure 13. (a) Plotted curve; (b) Reconstructed curve, using SPLINE command from AutoCAD; (c) Reconstructed curve using circle arcs; (d) Superimposed curves from (b) and (c).
Figure 13. (a) Plotted curve; (b) Reconstructed curve, using SPLINE command from AutoCAD; (c) Reconstructed curve using circle arcs; (d) Superimposed curves from (b) and (c).
Sensors 17 00003 g013
Figure 14. Measured and original posture for C2 posture (a) and A4 posture (b); Graphical representation of errors (c).
Figure 14. Measured and original posture for C2 posture (a) and A4 posture (b); Graphical representation of errors (c).
Sensors 17 00003 g014
Table 1. Equations used by the mathematical model.
Table 1. Equations used by the mathematical model.
Single Max Point Curve (“C” Shape)Double Max Point Curve (“S” Shape)Equation Number
l 1 = R 1 × ( p 1 p 2 )
l 2 = R 2 × ( p 2 + p 3 )
l 3 = R 3 × ( p 4 p 3 )
l 4 = R 4 × ( p 5 p 4 )
l 1 = R 1 × ( p 1 p 2 )
l 2 = R 2 × ( p 2 + p 3 )
l 3 = R 3 × ( p 3 p 4 )
l 4 = R 4 × ( p 4 + p 5 )
(1)
d 1 x = R 1 × ( sin ( p 1 ) sin ( p 2 ) )
d 1 y = R 1 × ( cos ( p 2 ) cos ( p 1 ) )
d 1 x = R 1 × ( sin ( p 1 ) sin ( p 2 ) )
d 1 y = R 1 × ( cos ( p 2 ) cos ( p 1 ) )
(2)
d 2 x = R 2 × sin ( p 2 )
d 2 y = R 2 × ( 1 cos ( p 2 ) )
d 2 x = R 2 × sin ( p 2 )
d 2 y = R 2 × ( 1 cos ( p 2 ) )
(3)
d 3 x = R 2 × sin ( p 3 )
d 3 y = R 2 × ( 1 cos ( p 3 ) )
d 3 x = R 2 × sin ( p 3 )
d 3 y = R 2 × ( 1 cos ( p 3 ) )
(4)
d 4 x = R 3 × ( sin ( p 4 ) sin ( p 3 ) )
d 4 y = R 3 × ( cos ( p 3 ) cos ( p 4 ) )
d 4 x = R 3 × ( sin ( p 3 ) sin ( p 4 ) )
d 4 y = R 3 × ( cos ( p 4 ) cos ( p 3 ) )
(5)
d 5 x = R 4 × ( sin ( p 5 ) sin ( p 4 ) )
d 5 y = R 4 × ( cos ( p 4 ) cos ( p 5 ) )
d 5 x = R 4 × ( sin ( p 4 ) )
d 5 y = R 4 × ( 1 cos ( p 4 ) )
(6)
d 6 x = R 4 × ( sin ( p 5 ) )
d 6 y = R 4 × ( 1 cos ( p 5 ) )
(7)
x 2 = d 1 x
y 2 = d 1 y
x 2 = d 1 x
y 2 = d 1 y
(8)
x m a x = x 2 + d 2 x
y m a x = y 2 + d 2 y
x m a x 1 = x 2 + d 2 x
y m a x 1 = y 2 + d 2 y
(9)
x 3 = x m a x + d 3 x
y 3 = y m a x d 3 y
x 3 = x m a x 1 + d 3 x
y 3 = y m a x 1 d 3 y
(10)
x 4 = x 3 + d 4 x
y 4 = y 3 d 4 y
x 4 = x 3 + d 4 x
y 4 = y 3 d 4 y
(11)
x 5 = x 4 + d 5 x
y 5 = y 4 d 5 y
x m a x 2 = x 4 + d 5 x
y m a x 2 = y 4 d 5 y
(12)
x 5 = x m a x 2 + d 6 x
y 5 = y m a x 2 + d 6 y
(13)
x c 1 = R 1 × sin ( p 1 )
y c 1 = ( R 1 × cos ( p 1 ) )
x c 1 = R 1 × sin ( p 1 )
y c 1 = ( R 1 × cos ( p 2 ) d 1 y )
(14)
x c 2 = d 1 x + R 2 × sin ( p 2 )
y c 2 = ( R 2 d 1 y d 2 y )
x c 2 = d 1 x + R 2 × sin ( p 2 )
y c 2 = ( R 2 d 1 y d 2 y )
(15)
x c 3 = d 1 x + d 2 x + d 3 x R 3 × sin ( p 3 )
y c 3 = ( R 3 × cos ( p 3 ) d 1 y d 2 y + d 3 y )
x c 3 = d 1 x + d 2 x + d 3 x + R 3 × sin ( p 3 )
y c 3 = R 3 × cos ( p 3 ) + d 1 y + d 2 y d 3 y
(16)
x c 4 = d 1 x + d 2 x + d 3 x + d 4 x + d 5 x R 4 × sin ( p 5 )
y c 4 = ( R 4 × cos ( p 4 ) d 1 y d 2 y + d 3 y + d 4 y )
x c 4 = d 1 x + d 2 x + d 3 x + d 4 x + R 4 × sin ( p 4 )
y c 4 = R 4 + d 1 y + d 2 y d 3 y d 4 y d 5 y
(17)
Table 2. Estimation of the postures that can be reproduced by the mathematical model.
Table 2. Estimation of the postures that can be reproduced by the mathematical model.
NameCriteriaResults
A4Max point between P4 and P5;
p1 > p2,
p2 > p3,
p3 > p4,
p4 < p5.
Sensors 17 00003 i015
C2Max point between P2-P3 and P4-P5;
p1 > p2,
p2 < p3,
p3 > p4,
p4 < p5.
Sensors 17 00003 i016
Table 3. Inertial sensors comparison.
Table 3. Inertial sensors comparison.
CharacteristicBosch BNO055MPU 9150Flora LSM9DS0AltIMU 10 v.4MinIMU 9 V3
Accelerometer
Measurement range±2 g to ±16 g±2 g to ±16 g±2 g to ±16 g±2 g to ±16 g±2 g to ±16 g
Sensitivity1 LSB/g2 LSB/g0.732 LSB/g0.732 LSB/g0.732 LSB/g
Gyroscope
Measurement range±125°/s to ±2000°/s±250°/s to ±2000°/s±245°/s to ±2000°/s±245°/s to ±2000°/s±245°/s to ±2000°/s
Sensitivity16 LSB/°/s16.4 LSB/°/s70 LSB/°/s70 LSB/°/s70 LSB/°/s
Magnetometer
Measurement range±1300 μT (x-, y-axis); ±2500 μT (z-axis)±1200 µT200 µT to 1200 µT200 µT to 1200 µT200 µT to 1200 µT
Resolution0.3 μT0.3 μT0.48 μT0.47 μT0.47 μT
Temperature sensorYesYesYesNo, includes a 24 bits barometerNo
Size (mm)20 × 27 × 415.5 × 29 × 4Diameter of 16 mm, Width 0.8 mm25.4 × 12.7 × 2.5420 × 13 × 3
CommunicationHID-I2C/I2C/UARTI2CSPI/I2CI2CI2C
Supply voltage2.4 V to 3.6 V2.4 V to 3.46 V2.4 V to 3.6 V2.5 V to 5.5 V2.5 V to 5.5 V
Power managementYes, with three modesNoYesNoNo
Table 4. Descriptive statistics for posture C2 (quantitative data).
Table 4. Descriptive statistics for posture C2 (quantitative data).
StatisticZ1Z2Z3Z4Z5
Nbr. of observations1010101010
Minimum23.4606.470−35.110−0.98025.140
Maximum24.1308.170−33.0900.34031.440
1st Quartile23.6207.253−34.423−0.87828.123
Median23.7557.660−34.285−0.50028.790
3rd Quartile23.9008.028−33.865−0.21530.638
Mean23.7627.550−34.134−0.47728.962
Variance (n − 1)0.0500.3710.3690.1923.796
Standard deviation (n − 1)0.2230.6090.6070.4381.948
Table 5. Descriptive statistics for posture A4 (quantitative data).
Table 5. Descriptive statistics for posture A4 (quantitative data).
StatisticZ1Z2Z3Z4Z5
Nbr. of observations1010101010
Minimum24.06023.12016.4002.680−14.160
Maximum25.55023.92017.2804.000−12.740
1st Quartile24.47523.26516.7502.870−13.908
Median24.93023.34516.9553.130−13.510
3rd Quartile25.43023.48817.2103.463−13.263
Mean24.88923.40116.9483.238−13.536
Variance (n − 1)0.3340.050.0840.2110.188
Standard deviation (n − 1)0.5780.2230.2900.4600.434

Share and Cite

MDPI and ACS Style

Voinea, G.-D.; Butnariu, S.; Mogan, G. Measurement and Geometric Modelling of Human Spine Posture for Medical Rehabilitation Purposes Using a Wearable Monitoring System Based on Inertial Sensors. Sensors 2017, 17, 3. https://doi.org/10.3390/s17010003

AMA Style

Voinea G-D, Butnariu S, Mogan G. Measurement and Geometric Modelling of Human Spine Posture for Medical Rehabilitation Purposes Using a Wearable Monitoring System Based on Inertial Sensors. Sensors. 2017; 17(1):3. https://doi.org/10.3390/s17010003

Chicago/Turabian Style

Voinea, Gheorghe-Daniel, Silviu Butnariu, and Gheorghe Mogan. 2017. "Measurement and Geometric Modelling of Human Spine Posture for Medical Rehabilitation Purposes Using a Wearable Monitoring System Based on Inertial Sensors" Sensors 17, no. 1: 3. https://doi.org/10.3390/s17010003

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