 Research
 Open Access
 Published:
A biarticular model for scapularhumeral rhythm reconstruction through data from wearable sensors
Journal of NeuroEngineering and Rehabilitation volume 13, Article number: 40 (2016)
Abstract
Background
Patientspecific performance assessment of arm movements in daily life activities is fundamental for neurological rehabilitation therapy. In most applications, the shoulder movement is simplified through a socketball joint, neglecting the movement of the scapularthoracic complex. This may lead to significant errors. We propose an innovative biarticular model of the human shoulder for estimating the position of the hand in relation to the sternum. The model takes into account both the scapulartoracic and glenohumeral movements and their ratio governed by the scapularhumeral rhythm, fusing the information of inertial and textilebased strain sensors.
Method
To feed the reconstruction algorithm based on the biarticular model, an adhoc sensing shirt was developed. The shirt was equipped with two inertial measurement units (IMUs) and an integrated textile strain sensor. We built the biarticular model starting from the data obtained in two planar movements (arm abduction and flexion in the sagittal plane) and analysing the error between the reference data  measured through an optical reference system  and the socketball approximation of the shoulder. The 3D model was developed by extending the behaviour of the kinematic chain revealed in the planar trajectories through a parameter identification that takes into account the body structure of the subject.
Result
The biarticular model was evaluated in five subjects in comparison with the optical reference system. The errors were computed in terms of distance between the reference position of the trochlea (endeffector) and the correspondent model estimation. The introduced method remarkably improved the estimation of the position of the trochlea (and consequently the estimation of the hand position during reaching activities) reducing position errors from 11.5 cm to 1.8 cm.
Conclusion
Thanks to the developed biarticular model, we demonstrated a reliable estimation of the upper arm kinematics with a minimal sensing system suitable for daily life monitoring of recovery.
Background
Dailylife evaluation of patient’s functional performance is essential for optimal guidance in neurological rehabilitation therapy and requires unobtrusive and ambulatory monitoring systems. An accurate estimation of the hand position with respect to the sternum is important to assess the recovery induced by the treatment and to prevent compensatory movements (e.g. spine flexion or scapular elevation in reaching tasks to compensate the functional deficiency). From a technological perspective, the challenge is to obtain a reliable and quantitative evaluation with a minimal set of unobtrusive sensors.
As widely described in the kinesiology literature [1, 2], the human shoulder can be considered as a biarticular six degrees of freedom kinematic chain. Flexion on the sagittal plane and abduction of the arm can be described as the composition of the movements of the scapularthoracic joint or scapular girdle (i.e. the scapula sliding and spinning on the rib cage), the glenohumeral joint (i.e. the humerus spinning, rolling and sliding into the glenoid cavity), and the spine. In addition, the two joints move together in a precise variable ratio, according to the scapularhumeral rhythm described by Crosbie et al. in [3]. In the last few decades several studies on shoulder mechanics have been performed and several models have been developed with different levels of complexity. In 1987, Högfors et al. [4] provided a significant description of the kinematics of the shoulder, based on coordinate frames fixed with bones of the scapular girdle, starting from human dissection. Two years later, A. E. Engin et al. [5] presented their “ThreeDimensional Kinematic Modelling of the Human Shoulder Complex” based on sinus cone joints and trajectories. In 1994, Van der Helm [6] proposed his finite elements model for the shoulder mechanism, introducing the concept of rolling on differentiable manifolds for joint surfaces. The following year, Happee et al. [7] presented a study on shoulder modelling including 95 muscles modelled as a thirdorder non linear system.
Despite the complex mechanical structure of shoulder, most research works performed the kinematic reconstruction of the upper limb by considering the shoulder as a socketball joint and neglecting the scapular contribution. In 2001, Prokopenko et al. [8] assessed a 7DOF model of the upper limb where the shoulder is modelled as a socketball joint by using a laboratory setup based on an electromagnetic motion tracker. Similar approaches, setups and results were obtained by Biryukova [9], Hingtgen [10] and Rab [11]. Few works considered the multiarticular approach. Among the others, we report the work of Retting [12] that assessed, through an optical system and skin applied markers, an advanced model of the shoulder on a large number of subjects. To the best of our knowledge, very few works combine the biarticular approach with an ambulatory measurement setup. A relevant example is the work of Cutti et al. [13] that estimated the scapulothoracic and humerothoracic movements by using three inertial measurements units (placed on the arm, the sternum and the scapula).
To promote a step forward in unobtrusive human motion analysis, De Rossi and Veltink proposed an hybrid approach based on the fusion of data derived from inertial and etextile sensors [14]. This strategy have since been implemented in studies and research projects to obtain patient’s performance assessments in daily life through textile integrated sensing [15, 16].
Inertial measurement units (IMUs) combine the information of accelerometers, gyroscopes and magnetometers, and are now widely used in wearable motion tracking [17–19]. The use of different IMUs, placed on connected body segments, and the additional information on the kinematic constraints enable most joint angles to be measured [20, 21].
In the last decade, several works have reported about textilebased or etextile sensors for unobtrusive human motion detection. Notable contributes include [22–27]. Textilebased sensors have several advantages: low cost, lightweight, low thickness, flexibility, and the possibility of adapting them to different body structures. The main drawbacks are the reduced accuracy, the nonnegligible transient time and the hysteresis. In our recent works [28–30] we employed and characterized textilebased sensors based on knitted piezoresistive fabrics (KPF) that shown reliable performances as strain and angular transducers.
In this work, we developed a new method to estimate upper arm kinematics, based on a biarticular model of the shoulder, combining the widely used socketball model  which mimics the glenohumeral joint  with an additional joint capable of describing the movement of the scapularthoracic complex and taking into account the constraint given by the scapularhumeral rhythm. According to the model requirements, we developed a dedicated sensing shirt provided with two IMUs  to measure the armtosternum orientation  and a KPF strain sensor able to to detect the scapular sliding. Note that the armtosternum orientation (IMUs output) is a combination of the scapularthoracic and glenohumeral movements (i.e. governed by the scapularhumeral rhythm), while the strain sensor output depends only on the scapularthoracic movement. The shoulder model we developed in this paper is a general mathematical structure based on differential rototranslational geometry whose parameters are identified in subjects using nonlinear constraint optimizations. We built the biarticular model starting from the data obtained in two planar movements (arm abduction and flexion in the sagittal plane) and analysing the error between the reference data  measured through an optical system  and the socketball approximation of the shoulder. The introduced method remarkably improves the estimation of the position of the trochlea (and consequently the estimation of the hand position during reaching activities) reducing position errors from 11.5 cm to 1.8 cm. The novelty of our method consisted in the embodiment of the scapularhumeral rhythm  expressed as a subject dependent relation between glenohumeral and scapulothoracic movements  in the biarticular model. This additional constraint allowed us to perform the reconstruction with a reduced sensor set (i.e. two IMUs and a monodimensional textileintegrated sensor) thus increasing the potential usability in dailylife.
Methods
The scapularhumeral rhythm
The scapula and humerus move simultaneously in a precise variable ratio: the scapularhumeral rhythm (Crosbie et al., [3]). The proportion of such movements varies according to the instantaneous position of the kinematic chain. The scapularhumeral rhythm is highly subject dependant through the particular morphology of the joint surfaces and the mechanical properties of the periarticular tissue. External loads may also influence the scapularhumeral rhythm.
A widely accepted approximation, described in [1], summarises the arm abduction and flexion on a sagittal plane as follows:

In the initial phase of the arm abduction, from 0° to 100°−110°, the movement of the scapula is overall about 10°−20°, while the humerus joint moves 90° with respect to the scapula. At the threshold of 100°−110°, the glenohumeral joint goes into the closed pack position. The remaining abduction of the arm, up to 150°, is entirely due to a scapular movement. In the final step, the arm reaches the vertical position through a hyperextension of the spine.

The arm flexion is composed of a spinning movement of the glenohumeral joint of about 60° along with a 60° sliding of the scapula on the rib cage through the compliance of the subscapularis muscle. Again, the missing 60° necessary to reach the vertical position are due to spine movements.
The ratio between scapular and glenohumeral movements can be expressed in a mathematical form as:
and
where ϕ _{ arm },ψ _{ arm } represent the flexion on a sagittal plane and abduction angles of the upper arm and ϕ _{ a },ψ _{ a } are the abductionadduction and the intraextrarotation angles of the scapula (i.e. under the hypothesis that the coordinate systems for sternum, scapula and humerus return ϕ _{ arm }=0 and ψ _{ arm }=0 in anatomical position).
If the relations 1 and 2 were exact and not dependent on the particular subject, the measurement of ϕ _{ arm },ψ _{ arm } would be sufficient to describe the whole shoulder movement. Since the body structure of a particular subject influences the mechanics of his/her shoulder, relationships 1 and 2 are approximations which may lead to significant errors. However, we used the described relations as the starting point for the personalised refinement of the shoulder mechanics of individual subjects.
Instrumentation
We collected the data derived from two IMUs (MTw produced by XSens [31]), placed on the sternum and arm respectively, and a textilebased strain sensor positioned on the back, from the spine to scapula, as shown in Fig. 1. The wearable sensors were integrated in a sensing shirt specifically produced and tested for clinical applications in stroke recovery, as we have previously described in [15, 16]. Simultaneously, we measured the 3D position of repere points (sternum, acromion, medial epicondyle of the trochlea of humerus) through an optical system (Smart DX 100 produced by BTS Bioengineering [32]) according to the marker setup shown in Fig. 1.
The MTw used in our experiments are miniature IMUs containing triaxial accelerometers, gyroscopes and magnetometers. The MTw applies proprietary fusion algorithms to perform real time estimation of the sensor orientation. The orientation and the raw sensor data (accelerometer, gyroscope, magnetometer) are sent to a remote receiver through a wireless link. We performed the IMU to segment alignment by asking the subject to stand upright in the anatomical position with the arms neutral beside the body. The transformation from the sensor to the body segment was then calculated by matching the sensor orientations in the global frame with the know orientations of the upper arm and sternum in the considered anatomical pose, as indicated by the sensor producer in [33]. We evaluated the arm flexion on a horizontal plane (θ _{ arm }) and abduction (ψ _{ arm }) from the relative orientation of the arm IMU with respect to the sternum one.
The strain sensor was made of a knitted piezoresistive fabric (KPF) produced by Smartex [34–36]. As we demonstrated in [27], the electrical resistance of a single layer KPF sample (R) can be approximated, for small elongations, as a linear function of the sensor length (l) and curvature (Δ α):
where k _{1},k _{2} are positive constants that depend on the properties of the sensing textile. In this work, we used the KPF strain sensor to detect the scapula sliding and, given the sensor configuration shown in Fig. 1, we considered negligible the contribute of the curvature (Δ α) in Eq. 3.
Data collection and initial analysis
Five volunteers without any pathology in their shoulders were asked to repeat (five times) two different movements involving the scapular girdle:

flexions of the arm on a sagittal plane in the range 0°−120° (starting from the anatomical position),

abductions of the arm on a frontal plane in the range 0°−150° (starting from the anatomical position).
We fed the socketball model of the shoulder with the IMU data (i.e. the arm horizontal flexion and abduction angles) to estimate the trochlea position thanks to the knowledge of the length of the anatomical segments. We then compared the estimated trochlea trajectories with the reference trajectories obtained by the optical system through the trochlea marker.
In the following subsections we will analyse the errors introduced by the socketball approximation in the described arm movements. This analysis will represent the starting point for the development of our refined biarticular model of the shoulder.
Flexions on a sagittal plane
The trochlea trajectory obtained by the socketball model was plotted in projection on a frontal plane (Fig. 2 a). The socketball model returned circular trajectories with the center in a point that undergoes small oscillations. These oscillations can be explained as a reaction of the trunk to the arm movements induced to maintain the balance of the center of mass and are more evident in Fig. 2 c where the estimated trochlea trajectory is projected on a horizontal plane. After having mathematically removed the trunk oscillations, by subtracting the sternum movement, we obtained a perfect circular trajectory for the trochlea (Fig. 3 a). Figure 2 b and d shows the acromion and the trochlea reference trajectories detected by the optical system. Figure 3 b, again obtained after removing the sternal oscillation, highlights the difference with respect to the socketball case (i.e. the reference trajectory of the trochlea is not circular). The discrepancy is above all due to the scapular movements that are neglected in the socketball model. We then interpolated the acromion trajectory with an elliptic arc, by committing errors less than 1 cm. Subsequently, we interpolated the trochlea trajectory as a generalized cycloid with a constant radius (i.e. the length of the humerus) and the center moving on the ellipse spanned by the humeral head (note that the humeral head is in a fixed position with respect to the acromion). The differences between the socketball model reconstruction and the reference data, computed for fixed angles on the repetitions of the flexion movements, are reported in Table 1.
Abductions on a frontal plane
The same procedure was performed to gather data on the frontal plane during an abduction. Figure 4 a and c shows the trochlea trajectory obtained by the socketball model in the arm abduction experiments, again projected in the frontal and horizontal planes. The reference trajectories for the same movements are plotted in Fig. 4 b and d. After having subtracted the sternum oscillations, Fig. 5 a and b highlight the difference between the reference data and the socketball model (i.e. perfect circular trajectory of the trochlea). Also in this case we represented the trochlea trajectory as a generalized cycloid whose center moves on a curve which can be approximated by an elliptic arc by committing errors in the range of 1 cm. Table 2 reports the errors committed by modelling the shoulder as a socketball joint computed for fixed angles on the repetitions of the abduction movements. Following these considerations, we can improve the description of the double movement of the scapula and humerus according to the scapular rhythm. Figure 6 shows a frame of the estimated output, where the red line represents the trajectory of the acromion and the green line is the trochlea reference path that clearly differs from the one predicted by the socketball model (large circle, in red). The shoulder joint movements during the arm flexion, previously described, can be considered in analogous way. These considerations led us to the development of a 3D model which fits the collected data.
Model creation and data fusion
Leveraging on the findings obtained for the two planar movements, we built the function that estimates the threedimensional trajectory of the trochlea by fusing the IMUs and strain sensor outputs, formally expressed as:
where θ _{ arm },ψ _{ arm } represent the horizontal flexion and abduction angles of the arm detected through the arm and sternum IMUs, R indicates the output of the strain sensor and X _{ Tr },Y _{ Tr }, Z _{ Tr } are the coordinates of the trochlea evaluated with respect to the sternum (see Fig. 7).
Figures 3 b and 5 b suggest that the acromion moves on elliptic trajectories. Thus, the 3D acromion state manifold can be described, with respect to the frame fixed with the sternum, by equation:
where A _{1,0}, A _{2,0}, A _{3,0} are the coordinates of the centre of the ellipsoid containing the acromion trajectories and A _{1,1}, A _{2,1}, A _{3,1} are the ellipsoid semiaxes. θ _{ a } and ψ _{ a } are the two variables that describe the scapula orientation with respect to the sternum and are dependent on θ _{ arm } and ψ _{ arm } through the scapularhumeral rhythm mechanism.
Assuming that the rotation center of the humeral head is in a fixed position with respect to the acromion and considering that the mutual position of the trochlea and the rotation center of the humerus remain unchanged, the trochlea coordinates can be expressed by:
where B _{1},B _{2},B _{3} represent the coordinates of the rotation centre of the humeral head with respect to a reference frame fixed with the acromion and D _{ ht } represents the distance between the center of the humeral head and the trochlea. Assuming that B _{1},B _{2},B _{3} are negligible with respect to A _{1,0},A _{2,0},A _{3,0} (which are comparable with the scapula dimension), by combining Eqs. 5 and 6 we obtain:
Note that, at this stage, A _{ i,j } with i=1,2,3;j=0,1 and D _{ ht } are unknown parameters for the model that varies among the different subjects.
To take into account the subjectdependent relation between scapular (θ _{ a }, ψ _{ a }) and glenohumeral (θ _{ arm }, ψ _{ arm }) movements, the mechanical model proposed by Eq. 7 can be combined with the humeralscapular rhythm relations (Eqs. 1 and 2). In addition, it is possible to introduce the information R derived from the strain sensor placed on the scapula. Assuming that the behaviour of the sensor is linear, its value can be related to the acromion position with respect to the frame fixed with the sternum. Under this hypothesis, we can introduce an additional constraint to the F determination, expressed by equations:
where c _{0,1}, c _{1,1}, c _{0,2} c _{1,2},c _{0,3} and c _{1,3} are six coefficients depending on the structure of the body of the person wearing the sensing garment. Equation 8 introduces new parameters and increases the computational cost of the identification of F, however it reinforces the relation between the shoulder movement and the scapular sliding. With this last step, the whole sensing system becomes redundant. This redundancy produces an overdetermined problem. Hence the solution can be found in the leastsquare sense by solving the variable surplus that determines the system parameters.
In summary, the construction of F can be obtained by considering Eq. 7 under the constraints given by Relations 1, 2 and 8. The unknown parameters A _{ i,j }, D _{ ht } and c _{ i,j } can be identified by the minimization of the following functional:
on the integration domain:
where \(\hat {X}_{Tr}\), \(\hat {Y}_{Tr}\) and \(\hat {Z}_{Tr}\) are the real coordinates of the trochlea derived by the optical system. The minimization of the functional 9 is executed through Lagrange’s multiplier method under the constraints, 1, 2 and 8. Symbolic computation was executed by MapleV to obtain formal derivatives, while numerical computation for solving Lagrange’s equations was performed using in MATLAB. After identifying A _{ i,j } s, c _{ i,j } s and D _{ ht }, the function F is completely determined as an iterative process which corrects the initial socketball evaluation through the strain sensor and the double joint movement. The status manifold of the coordinates of the trochlea is shaped on the particular subject wearing the shirt, and the creation of F through its parameter identification needs to be executed for each different subject using the described measurement setup. Figure 8 shows the action of F. Starting from the sphere (which corresponds to the socketball status manifold), the biarticular ellipsoid model described by Eq. 7 is obtained through F.
Results
We assessed the performance of F on the five subjects in a series of upper limb vertical movements lying on different planes. The subjects repeated the movements three times.
We compared the estimated trochlea position  obtained by fusing IMUs and strain data through F  with the reference trajectory obtained by the optical system (trochlea marker). Table 3 reports the results of this comparison grouped in mean and standard deviation of the distance between the estimated and reference position of the trochlea, computed over the entire set of trials for each arm position. The maximal value of the average error is 1.8 cm and is obtained for 40° horizontal flexion and 120° abduction.
In addition, to highlight the improvement of our biarticular model, we evaluated the error of the socketball approach within the same set of trials. Table 4 reports the socketball error expressed as the mean and standard deviation of the distance between the socketball estimation and the reference position of the trochlea, over the entire set of trials for each arm position. As shown in Table 4, the socketball maximum error is greater than 11.5 cm which is 6 times higher than the biarticular model error reported in Table 3.
Figure 9 shows a visual comparison between the biarticular and socketball errors, again highlighting the better performance of the estimation through F.
Discussion
Our biarticular model remarkably improved the reconstruction of the position of the endeffector (trochlea) with respect to the simple socket ball approximation. The measurement of the elbow kinematics, out of the scope of this research, would allow to extend these results to the estimation of the hand position in reaching activities. Considering the elbow completely extended and an average armlength of about 80 cm, the error in hand position can be roughly extrapolated: 3.5 cm for the biarticular model vs. 20 cm for the socket ball approximation. The bidimensional error map (Fig. 9) is roughly symmetrical with respect to 40° of horizontal flexion where the error is maximum (corresponding to the 40°column of Table 3). This is due to the estimation of the model parameters described in the Methods section. Indeed, the integral to be minimized is defined on the domain I (Eq. 10) and the points of the shoulder workspace corresponding to the 40°column are the farthest from the ones used for the identification process.
The obtained results can be considered satisfactory for a reliable monitoring of reaching activity in dailylife context. For a comparison with the recent literature, we considered the contribute of Van Meulen et al. reported in [37]. Here, the authors validated the upper limb reconstruction accuracy of the IMUbased MVN Biomech motion capturing system [33] against an optical reference system. They obtained a maximum error around 6 cm and they were able, through the described setup, to correlate the patient’s reaching performance with the upper extremity part of the FuglMeyer scale [38], in a simulated inhome task. Note that MVN Biomech estimates the shoulder movement through a biarticular approach with three IMUs, following a configuration similar to the one of Cutti et al. [13], already discussed in the Background.
Our methodology has the potential to be applied in monitoring the quality of the movement and evaluating the recovery obtained by adequate physical training. In stroke patients, nerve deafferentation influences the muscle control in different ways. Proximal muscles, with a cranial innervation are easier to control than distal muscles. This phenomenon can activate pathological synergies when performing activities which need to be prevented and corrected. A classical case is the replacement of the use of the deltoid muscle in antigravitational actions by a torque exerted by trapezius and serratus anterioris muscles. This pathological motor scheme produces an excessive rotation of the scapula to compensate for the scarce elevation of the humerus, thus violating the rhythm. Other pathological conditions that affects the scapularhumeral rhythm are the subacromial and internal shoulder impingement [39] that may represent a future objective of our analysis.
Several limitations of our work should be acknowledged. First, the upper limb reconstruction through our biarticular model and setup can be affected by measurement errors introduced by the textilebased strain sensor, which may present a nonnegligible hysteresis [27]. A more intensive testing phase on an higher number of subjects would be needed to better quantify these aspects. However, our method is general and highlights the importance of measuring the scapular sliding independently on the particular strain sensor used (i.e. the more accurate the strain, the better the reconstruction). Secondly, the scapularhumeral rhythm  embodied in our biarticular model  depends on the external load and this may affect our reconstruction performance. Indeed, Yoshiaki et al. demonstrated that a 3kg load considerably modified the scapulothoracic vs. glenohumeral ratio [40]. However, we are confident that for reduced loads  typical of our target application  our approximation may still be valid. Finally, the calibration procedure for the identification of the model parameters requires the acquisition of the endeffector position in planar movements through an external referenced system. This aspect may limit the potential use of our system, especially in daily life contexts where there is the need of simple and inexpensive motion tracking systems. In future works, a simplified calibration based on a reduced number of predefined positions will be investigated.
Conclusions
This paper presented a new method for the reconstruction of the upper limb kinematics starting from minimal set of inertial and textile sensors aimed at ambulatory monitoring of the patient’s reaching activity in daily life. The idea consisted in using a biarticular approach, describing both the glenohumeral and scapularthoracic joints, and in implementing a constraint given by the subject dependent ratio between the two joint’s movements, governed by the scapularhumeral rhythm. We demonstrated a reliable estimation of the handtosternum position with a minimal sensing system, much less obtrusive than the existing ambulatory and multiarticular approaches.
Ethical statement
The study protocol is a subset of a larger protocol developed within the INTERACTION (FP7ICT20117287351) project and was carried out in accordance with the recommendations of Swiss medic, medical and ethical committee. All participants were healthy subjects and gave written informed consent in accordance with the declaration of Helsinki.
References
 1
Kapanji IA. the physiology of the joints: Upper limb. London (UK): Churchill Livingstone; 1983.
 2
Freivalds A. Biomechanics of the upper limb. Boca Raton FL (US): CRC Press; 2004.
 3
Crosbie J, Kilbreath SL, Hollmann L, York S. Scapulohumeral rhythm and associated spinal motion. Clin Biomech. 2008; 23(2):184–92.
 4
Högfors C, Sigholm G, Herberts P. Biomechanical model of the human shoulderI, Elements. J Biomech. 1987; 20(2):157–66.
 5
Engin AE, Tümer ST. Threedimensional kinematic modelling of the human shoulder complex  part I: Physical model and determination of joint sinus cones. J Biomech. 1989; 11(2):107–12.
 6
Van der Helm FCT. A finite element musculoskeletal model of the shoulder mechanism. J Biomech. 1994; 27(5):551–69.
 7
Happee R, Van der Helm FCT. The control of shoulder muscles during goal directed movements, an inverse dynamic analysis. Med Eng Phys. 1995; 8(10):1779–991.
 8
Prokopenko RA, Frolov AA, Biryukova EV, RobyBrami A. Assessment of the accuracy of a human arm model with seven degrees of freedom. J Biomech. 2001; 34(2):177–85.
 9
Biryukova EV, RobyBrami A, Frolov AA, Mokhtari M. Kinematics of human arm reconstructed from spatial tracking system recordings. J Biomech. 2000; 33(8):985–95.
 10
Hingtgen B, McGuire JR, Wang M, Harris GF. An upper extremity kinematic model for evaluation of hemiparetic stroke. J Biomech. 2006; 39(4):681–8.
 11
Rab G, Petuskey K, Bagley A. A method for determination of upper extremity kinematics. Gait Posture. 2002; 15(2):113–9.
 12
Retting O, Fradet L, Kasten P, Raiss P, Wolf S. A new kinematic model of the upper extremity based on functional joint parameter determination for shoulder and elbow. Gait Posture. 2009; 30(4):469–76.
 13
Cutti A, Giovanardi A, Rocchi L, Davalli A, Sacchetti R. Ambulatory measurement of shoulder and elbow kinematics through inertial and magnetic sensors. Med Biol Eng Comput. 2008; 46(2):169–78.
 14
Veltink PH, Rossi DD. Wearable technology for biomechanics: etextile or micromechanical sensors?[conversations in bme]. IEEE Eng Med Biol Mag. 2010; 29(3):37–43.
 15
Tognetti A, Lorussi F, Carbonaro N, De Rossi D, De Toma G, Mancuso C, Paradiso R, Luinge H, Reenalda J, Droog E, Veltink PH. Dailylife monitoring of stroke survivors motor performance: The interaction sensing system. In: 36th Annual International Conference of the IEEE Engineering in Medicine and Biology Society. US: IEEE: 2014. p. 4099–102.
 16
Klaassen B, van Beijnum BJ, Weusthof M, Hofs D, van Meulen F, Droog E, Luinge H, Slot L, Tognetti A, Lorussi F, Paradiso R, Held J, Luft A, Reenalda J, Nikamp C, Buurke J, Hermens H, Veltink P. A full body sensing system for monitoring stroke patients in a home environment. Commun Comput Inf Sci. 2015; 511:378–93.
 17
Roetenberg D, Luinge HJ, Baten C, Veltink PH. Compensation of magnetic disturbances improves inertial and magnetic sensing of human body segment orientation. Neural Syst Rehabil Eng IEEE Trans. 2005; 13(3):395–405.
 18
Sabatini AM. Estimating threedimensional orientation of human body parts by inertial/magnetic sensing. Sensors. 2011; 11(2):1489–525.
 19
Liu K, Liu T, Shibata K, Inoue Y, Zheng R. Novel approach to ambulatory assessment of human segmental orientation on a wearable sensor system. J Biomech. 2009; 42(16):2747–752.
 20
Luinge HJ, Veltink PH, Baten CT. Ambulatory measurement of arm orientation. J Biomech. 2007; 40(1):78–85.
 21
Roetenberg D, Slycke PJ, Veltink PH. Ambulatory position and orientation tracking fusing magnetic and inertial sensing. IEEE Trans Biomed Eng. 2007; 54(5):883–90.
 22
De Rossi D, Carpi F, Lorussi F, Paradiso R, Scilingo E, Tognetti A. Electroactive fabrics and wearable manmachine interfaces. Wearable Electronics Photonics. Cambridge (UK): Woodhead Publishing; 2005, pp. 59–80.
 23
Gibbs PT, Asada HH. Wearable conductive fiber sensors for multiaxis human joint angle measurements. J NeuroEng Rehabil. 2005; 2(1):7.
 24
Mattmann C, Clemens F, Tröster G. Sensor for measuring strain in textile. Sensors. 2008; 8(6):3719–32.
 25
Gioberto G, Dunne L. Theory and characterization of a topthread coverstitched stretch sensor. In: IEEE International Conference on Systems, Man, and Cybernetics. US: IEEE: 2012. p. 3275–80.
 26
Vanello N, Hartwig V, Tesconi M, Ricciardi E, Tognetti A, Zupone G, Gassert R, Chapuis D, Sgambelluri N, Scilingo EP, et al.Sensing glove for brain studies: design and assessment of its compatibility for fmri with a robust test. IEEE/ASME Trans Mechatron. 2008; 13(3):345–54.
 27
Tognetti A, Lorussi F, Dalle Mura G, Carbonaro N, Pacelli M, Paradiso R, Rossi DD. New generation of wearable goniometers for motion capture systems. J NeuroEng Rehabil. 2014; 11(1):56.
 28
Dalle Mura G, Lorussi F, Tognetti A, Anania G, Carbonaro N, Pacelli M, Paradiso R, De Rossi D. Piezoresistive goniometer network for sensing gloves. In: IFMBE Proceedings. Zurich (Switzerland): Springer International Publishing: 2014. p. 1547–50.
 29
Carbonaro N, Dalle Mura G, Lorussi F, Paradiso R, De Rossi D, Tognetti A. Exploiting wearable goniometer technology for motion sensing gloves. IEEE J Biomed Health Inf. 2014; 18(6):1788–1795.
 30
Tognetti A, Lorussi F, Carbonaro N, de Rossi D. Wearable goniometer and accelerometer sensory fusion for knee joint angle measurement in daily life. Sensors (Switzerland). 2015; 15(11):28435–55.
 31
Xsens MTw. Available online: https://www.xsens.com/products/mtwdevelopmentkitlite/mtw/. Accessed on 1 March 2016.
 32
BTS SMARTDX. Available online: www.btsbioengineering.com/products/kinematics/btssmartdx. Accessed on 1 March 2016.
 33
Roetenberg D, Luinge H, Slycke P. Xsens MVN: full 6DOF human motion tracking using miniature inertial sensors. 2009. Available online: www.xsens.com/wpcontent/uploads/2013/12/MVN\_white\_paper1.pdf. Accessed on 1 March 2016.
 34
Smartex srl. Available online: www.smartex.it. Accessed on 1 March 2016.
 35
Pacelli M, Caldani L, Paradiso R. Textile piezoresistive sensors for biomechanical variables monitoring. In: 28th Annual International Conference of the IEEE Engineering in Medicine and Biology Society. US: IEEE: 2006. p. 5358–61.
 36
Pacelli M, Caldani L, Paradiso R. Performances evaluation of piezoresistive fabric sensors as function of yarn structure. In: 35th Annual International Conference of the IEEE Engineering in Medicine and Biology Society. US: IEEE: 2013. p. 6502–5.
 37
van Meulen FB, Reenalda J, Buurke JH, Veltink PH. Assessment of dailylife reaching performance after stroke. Ann Biomed Eng. 2015; 43(2):478–86.
 38
Fugl Meyer AR, Jaasko L, Leyman I. The post stroke hemiplegic patient. i. a method for evaluation of physical performance. Scand J Rehabil Med. 1975; 7(1):13–31.
 39
Kibler WB. Scapular involvement in impingement: signs and symptoms. Instr Course Lect. 2006; 55:35–43.
 40
Kon Y, Nishinaka N, Gamada K, Tsutsui H, Banks SA. The influence of handheld weight on the scapulohumeral rhythm. J Should Elb Surg. 2008; 17(6):943–6.
Acknowledgements
This research was supported by the EU 7th framework project INTERACTION (FP7ICT20117287351).
Author information
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
FL and AT drafted the paper, AT, FL and NC performed the data acquisition, FL formulated the mathematical model, AT and FL performed the data elaboration, NC and DDR revised the draft. All authors read and approved the final manuscript.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License(http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver(http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Lorussi, F., Carbonaro, N., Rossi, D.D. et al. A biarticular model for scapularhumeral rhythm reconstruction through data from wearable sensors. J NeuroEngineering Rehabil 13, 40 (2016) doi:10.1186/s1298401601492
Received
Accepted
Published
DOI
Keywords
 Reaching activity
 Scapular girdle movement
 Hand posture estimation
 Scapularhumeral rhythm
 Wearable sensing