Direct assessment of 3D foot bone kinematics using biplanar X-ray fluoroscopy and an automatic model registration method
© Ito et al. 2015
Received: 9 February 2015
Accepted: 1 June 2015
Published: 10 June 2015
Quantifying detailed 3-dimensional (3D) kinematics of the foot in contact with the ground during locomotion is crucial for understanding the biomechanical functions of the complex musculoskeletal structure of the foot. Biplanar X-ray fluoroscopic systems and model-based registration techniques have recently been employed to capture and visualise 3D foot bone movements in vivo, but such techniques have generally been performed manually. In the present study, we developed an automatic model-registration method with biplanar fluoroscopy for accurate measurement of 3D movements of the skeletal foot.
Three-dimensional surface models of foot bones were generated prior to motion measurement based on computed tomography. The bone models generated were then registered to biplanar fluoroscopic images in a frame-by-frame manner using an optimisation technique, to maximise similarity measures between occluding contours of the bone surface models with edge-enhanced fluoroscopic images, while avoiding mutual penetration of bones. A template-matching method was also introduced to estimate the amount of bone translation and rotation prior to automatic registration.
We analysed 3D skeletal movements of a cadaver foot mobilized by a robotic gait simulator. The 3D kinematics of the calcaneus, talus, navicular and cuboid in the stance phase of the gait were successfully reconstructed and quantified using the proposed model-registration method. The accuracy of bone registration was evaluated as 0.27 ± 0.19 mm and 0.24 ± 0.19° (mean ± standard deviation) in translation and rotation, respectively, under static conditions, and 0.36 ± 0.19 mm and 0.42 ± 0.30° in translation and rotation, respectively, under dynamic conditions.
The measurement was confirmed to be sufficiently accurate for actual analysis of foot kinematics. The proposed method may serve as an effective tool for understanding the biomechanical function of the human foot during locomotion.
Human bipedal locomotion is a mechanical phenomenon that transfers the body centre of mass by generating appropriate ground reaction forces. Therefore, quantifying the detailed 3-dimensional (3D) kinematics of the foot in contact with the ground during human bipedal locomotion is crucial for improving our understanding of the biomechanical functions of the complex foot musculoskeletal structures in the generation of stable and efficient bipedal locomotion. Many studies have been conducted to capture 3D kinematics of the foot bones during human walking, mainly using marker-based motion capture systems [1, 2]. However, such quantification is less accurate because of skin marker artefacts [3–6]. The use of bone pins allows direct measurement of 3D bone kinematics, but is very invasive [4–8]. Accurate measurement of the 3D kinematics of foot skeletal movement during locomotion thus remains a challenging problem.
Radiographic techniques provide a more direct, but less-invasive way to measure bone kinematics. Kinematic analyses of skeletal motions using single or biplanar X-ray fluoroscopic systems have recently become widespread in the field of biomechanics, yielding significant results. For example, dynamic kinematic measurements of the human knee , shoulder  and spine  have been conducted using X-ray fluoroscopy to capture precise joint kinematics in vivo. Furthermore, such techniques have been applied to the analysis of animal locomotion, bringing substantial benefits to uncovering morphofunctional relationships between animal anatomy and locomotion [12, 13]. In the same vein, X-ray fluoroscopic systems have recently been applied to analyses of 3D foot and ankle kinematics in vivo. For example, Yamaguchi et al.  and Fukano et al.  have employed single-plane fluoroscopy and Wan et al. , de Asla et al. , and Kozanek et al.  have used biplanar fluoroscopy to successfully reconstruct 3D kinematics of the ankle and subtalar joints in vivo.
In many of these foot studies, however, 3D surface models of the bones are manually matched to fluoroscopic images, by manipulating each bone in a virtual space to be positioned and rotated using a computer mouse. Manual registration is indeed challenging and time-consuming work that requires extreme patience, possibly resulting in inaccuracy of the reconstructed positions and orientations of the bones. The potential for inaccuracy in model registration may be relatively larger in the case of foot analyses, since bone registration is more difficult because the foot skeleton consists of 27 small, overlapping bones. Although many automatic registration methods have been proposed for relatively simple skeletal structures such as the knee [19–21], hip , shoulder  and spinal vertebrae , no studies have attempted to develop an accurate registration method for foot bones. Very recently, Campbell et al.  reported automatic registration of the tibia, fibula and calcaneus in walking, but did not include foot bones other than the calcaneus.
In the present study, we developed an automatic model-registration method with biplanar fluoroscopy for accurate measurement of 3D foot skeletal movements. To evaluate the feasibility and validity of the proposed method, we analysed 3D skeletal movements of a cadaver foot mobilized by the robotic gait simulator using a biplanar X-ray fluoroscopy system and a proposed model-registration method.
To reconstruct 3D movements of the foot bones from the 2D images, spatial calibration of the biplanar fluoroscopic system is necessary. We used two calibration objects with spherical metal markers (a planar L-shaped structure and a linear rigid bar) to obtain the projection geometry of the biplanar fluoroscopic system based on the direct linear transformation method , assuming that a fluoroscopic system can be modelled as an ideal perspective projection of a point-source X-ray onto a flat panel. Correction of non-linear distortion due to the photocathode curvature of the image intensifier was not necessary in the present study, as we used flat panel detectors, but we corrected for linear trapezoidal distortion of fluoroscopic images due to imperfect orthogonality of the X-ray axis to the flat panel detector. Mean and standard deviation of the calibration error was 0.076 ± 0.058 mm in the present study.
3D bone surface models
To create bone surface models registered to the fluoroscopic images, a cadaver foot was scanned using a medical computed tomography (CT) scanner (Aquilion ONE; Toshiba Medical Systems, Japan), and cross-sectional images were reconstructed at 0.25-mm intervals, with a pixel size of 0.316 mm. Images were then transferred to medical imaging software (Analyze version 9.0; Mayo Clinic, USA) for segmentation of the foot bones. Surface models of the bones were generated as triangular polygonal models by the marching cube method. The polygonal mesh models were then regenerated after polygon reduction using reverse-engineering software (RapidForm 2006; INUS Technology, Korea). A local coordinate system was created to describe the relative positions and orientations of each of the four tarsal bones, i.e., calcaneus, talus, cuboid, and navicular. Specifically, three or four landmarks were defined on each of the bones and an orthogonal coordinate system was calculated . The x-, y- and z-axes in the present study approximately correspond to the anteroposterior, mediolateral, and dorsoplantar directions, respectively. Origins of the coordinate systems were defined at the centroid of each polygonal mesh bone model. We used y-x-z Euler angles to describe the absolute and relative orientations of a bone, with ϕ, θ and ψ as the rotational angles around the y (plantarflexion-dorsiflexion), x (inversion-eversion), and z (adduction-abduction) axes, respectively.
Automatic model-registration method
where subscript i represents the bone number, and subscripts L and R denote left and right images, respectively. This similarity measure represents a summation of correlation coefficients between the two edge-enhanced images of the left and right fluoroscopic systems.
where D ij is the penetration depth between the ith and jth bones, and μ is the coefficient (=1000). The penalty function was calculated by a collision-detection library (SmartCollision version 2.3; I-NET Corporation, Japan), which can efficiently detect collisions and calculate penetration depths between 3D mesh models. We used a quasi-Newtonian method for optimisation, and the calculation was conducted with custom-made software using Microsoft Visual Studio C++ 2010. Due to the complexity of the edge-enhanced images from foot fluoroscopy, the optimisation function involves numerous local minima. In the present study, therefore, partial derivatives were computed with multiple different step sizes for this gradient-based optimisation algorithm.
The initial position of the bone model on the first frame was defined manually. Automatic registration was then performed to search for the true position and orientation of the bones on the first frame. This result was used as the initial guess for optimisation on the next frame and registration proceeds one after another. However, if movement of the bone is relatively fast compared to the frame rate of the fluoroscopic system, bone translation and rotation in two consecutive frames might be quite large, possibly resulting in tracking failure. Prior to the automatic registration, therefore, we roughly estimated the amount of bone translation and rotation in two consecutive frames using a template-matching method . Specifically, small areas (30 × 30 pixels) corresponding to bony landmarks were defined manually on the first frame as template images, and corresponding matching areas were automatically identified by comparing template images against the subsequent image based on cross-correlation. The same set of bony landmarks was traced on both fluoroscopic images. The 3D positions of bony landmarks were calculated based on triangulation for estimation of initial guesses for the position and orientation of each bone. Template images were updated for template-matching of the next frame.
To evaluate the accuracy of the proposed registration method, we used four articulated dry tarsal bones with four spherical metal markers attached to each bone. We placed the articulated tarsal bones in the fluoroscopic system in five different positions and orientations, and corresponding fluoroscopic images were obtained. We solved for the positions and orientations of each bone in two ways, one by the proposed registration method and the other using the 3D positions of the four markers calculated based on calibrated projection geometry. The positions and orientations of the bones calculated by the latter were taken as true values and the registration errors were evaluated (static evaluation). We also used a single unarticulated talus and calcaneus to see how the overlap of articulated tarsal bones in fluoroscopic images affects registration performance. In addition, we conducted the same evaluation experiment while manually translating and rotating the articulated tarsal bones to evaluate how movements affect registration performance.
The dry tarsal bones used in this experimental evaluation were obtained from the Laboratory of Physical Anthropology at Kyoto University, and were scanned using a medical CT scanner. Cross-sectional images were reconstructed at 0.5-mm intervals, with a pixel size of 0.234 mm. The bone models were constructed as described previously.
The present study reconstructed the dynamic skeletal movements of a cadaver foot (male, 78 years old) mobilized by a robotic gait simulator to evaluate how well the proposed method can track and register foot bones in motion during the stance phase of bipedal walking. The cadaver used in our study was donated to the Clinical Anatomy Laboratory at Keio University School of Medicine with the consent of the family. The present study was approved by the ethics committee of the School of Medicine and the Faculty of Science and Technology at Keio University.
Briefly, the simulator has three legs (fore, middle and hind legs), with the cadaver foot fixed to the middle leg arranged radially in the sagittal plane. The middle leg replicates the movements of the cadaver foot from heel-contact to toe-off, while the hind and fore legs replicate toe-off and heel-contact of the contralateral limb, respectively, so the X-ray image is not interrupted by the swing leg. The simulator initially stands only on the hind leg, and is then released to fall forward so that the middle cadaver foot makes heel contact. The simulator continuously rotates forward so that the fore leg makes the foot-ground contact and the cadaver foot finally makes toe-off. Tendons of the tibialis anterior and soleus were connected to pneumatic actuators to apply forces at the appropriate moment to reproduce how the foot would function during walking. Development of the robotic gait simulator will be described elsewhere.
Using fluoroscopic images of this cadaveric walking, we quantified 3D skeletal movements of the calcaneus, talus, navicular and cuboid in the stance phase of the gait and evaluated system performance. Template images were selected at the following six positions: 1) most anterodorsal point of the navicular; 2) most anterodorsal point of the talus; 3) most posterior point of the talus; 4) most superior point; 5) most posterior point; and 6) most inferior point of the calcaneal tuberosity (Fig. 2) for this analysis.
Mean (standard deviation) errors of articulated bones and unarticulated bones in static evaluation
Mean (standard deviation) errors of articulated bones and unarticulated bones in dynamic evaluation
In the present study, we developed an automatic model-registration method with biplanar fluoroscopy for accurate measurement of 3D foot skeletal movements, and demonstrated that 3D kinematics of the calcaneus, talus, navicular and cuboid in the stance phase of gait were successfully reconstructed using the proposed method. In this study, we used the similarity measure (Equation 1) for the registration of bone models with fluoroscopic images. However, the similarity measure alone did not produce satisfactory results, since foot bones overlapped heavily with one another on the fluoroscopic images. Furthermore, registration failed when bone translation and rotation in two consecutive frames were large. We therefore introduced the penetration penalty function, landmark tracking based on template matching, and partial derivatives of the objective function with multiple step sizes for our registration process. With these efforts, automatic tracking of the positions and orientations of foot bones from the two fluoroscopic images become possible as presented in Fig. 3.
The accuracy of foot bone registration was confirmed to be 0.27 ± 0.19 mm in translation and 0.24 ± 0.19° in rotation in the static condition, and 0.36 ± 0.19 mm in translation and 0.42 ± 0.30° in rotation in the dynamic condition. The registration error of published automatic model-registration methods was reportedly 0.16-0.33 mm and 0.16-0.58° for translation and rotation, respectively [20, 22, 23], indicating that similar levels of accuracy were obtained using the proposed method. In those previous studies, however, registrations were performed on relatively simple skeletal structures such as knee, hip and shoulder joints, where there is very little overlap between bones. Furthermore, we tried herein to measure movements of small foot bones such as the navicular and cuboid in which relatively few characteristic morphological features could be detected. We therefore believe that the accuracy of the proposed method should be evaluated as relatively good. Based on kinematic measurements using bone pins, the range of subtalar joint motion is reported to be about 2.5-3.8 mm in translation and 2.8-9.7° in rotation in cadaveric studies [29–31] during human walking, and about 3.3-9.8° from in vivo studies [6–8]. Range of motion of tarsometatarsal joints has been reported as about 4.6-12.9° in cadaveric studies [30, 31], and 4.1-13.3° from in vivo studies [6–8]. The present registration method should therefore be considered as sufficiently accurate for quantitative description of the talar and tarsometatarsal joint kinematics. It must be noted, however, that evaluation of accuracy using dry bones has probably underestimated the actual error, as bone edges are less distinct when surrounded by soft tissues. This should be confirmed in future studies.
Manual registration requires a great deal of patience, particularly if there are many small bones to be registered as in the case of the foot. On the other hand, the present study allows automatic registration of multiple bones while possibly eliminating some of the errors associated with the tediousness of manual registration. The accuracy of the biplanar manual registration is reportedly about 0.3-0.4 mm in translation and 0.4-0.5° in rotation [10, 11]. Consequently, the present framework may provide useful support for clarifying 3D kinematics of the foot in contact with the ground during human bipedal locomotion.
Recently, a digitally reconstructed radiograph (DRR) method has frequently been applied to assess skeletal movements from fluoroscopic images, as the DRR method offers more information for matching CT volumes to fluoroscopic images, since it can use internal bone morphology such as cortical thickness [9, 22, 25, 32]. Although the DRR method probably requires more processing time than the present contour-matching method, the DRR certainly provides more information than contours and theoretically yields better performance. Miranda et al.  recently proposed DRR-based tracking software using GPU acceleration and computational time is decreasing. Incorporating the DRR method in the present study may be an option in future studies to achieve better registration results.
However, the proposed method has some limitations. One is that the computational cost could become very high if the number of bones to be registered is large. Ideally, the positions and orientations of all 27 foot bones should be registered simultaneously. However, it is probably more practical to perform calculations in a step-by-step manner for simplicity, such as calculating the positions and orientations of the four bones as the first step, then proceeding to the metatarsals, cuneiforms, tibia, and so on, to reconstruct 3D skeletal movements of the entire foot. Another limitation is that the present registration process is not fully automatic, because the initial guess for the positions and orientations of bones on the first frame and the bony landmarks to be tracked should be determined manually. A simple method to place bone models in a vicinity of the true positions and orientations on the first frame seems necessary. For example, Haase et al.  proposed automatic annotation of landmarks in biplanar fluoroscopy based on active appearance models. Applying such a method may be useful for placing bones in the vicinity of the true position and orientation on the first frame. Lastly, although the present method can be applied to quantify foot skeletal movements in vivo, we quantified 3D skeletal movements of a cadaver foot mobilized by a robotic gait simulator using the proposed method. However, as shown in Fig. 4, reconstructed subtalar joint kinematics are currently not in good agreement with those of actual human walking as reported in the literature, mainly due to difficulties associated with the robotic gait simulator. More work is clearly needed to mobilize the cadaver foot in a manner similar to the way the human foot actually moves and mechanically interacts with the ground, to further elucidate the biomechanical functions of the complex musculoskeletal structures of the human foot during bipedal locomotion.
In the present study, we developed an automatic model-registration method with biplanar fluoroscopy for accurate measurement of 3D movements of the skeletal foot. Specifically, 3D surface models of foot bones were generated prior to motion measurement based on computed tomography, and the bone models generated were then registered to biplanar fluoroscopic images in a frame-by-frame manner to maximise similarity measures between occluding contours of the bone surface models with edge-enhanced fluoroscopic images, while avoiding mutual penetration of bones. Measurement was confirmed to be sufficiently accurate for actual analysis of foot kinematics. The proposed method may serve as an effective tool for understanding the biomechanical function of the human foot during locomotion.
We wish to express our sincere gratitude to Shimadzu Corporation, Kyoto, Japan, for the development of the biplanar X-ray fluoroscopic system used in the present study, and Hiroyuki Seki, Masateru Kitashiro and Taiki Kobayashi for assistance with the experiment. This study was financially supported by a Grant-in-Aid for Scientific Research (#10252610) from the Japan Society for the Promotion of Science.
- Carson MC, Harrington ME, Thompson N, O’Connor JJ, Theologis TN. Kinematic analysis of a multi-segment foot model for research and clinical applications : a repeatability analysis. J Biomech. 2001;34:1299–307.View ArticlePubMedGoogle Scholar
- Dixon PC, Böhm H, Döderlein L. Ankle and midfoot kinetics during normal gait: a multi-segment approach. J Biomech. 2012;45:1011–6.View ArticlePubMedGoogle Scholar
- Tranberg R, Karlsson D. The relative skin movement of the foot: a 2-D roentgen photogrammetry study. Clin Biomech. 1998;13:71–6.View ArticleGoogle Scholar
- Reinschmidt C, van Den Bogert AJ, Murphy N, Lundberg A, Nigg BM. Tibiocalcaneal motion during running, measured with external and bone markers. Clin Biomech. 1997;12:8–16.View ArticleGoogle Scholar
- Westblad P, Hashimoto T, Winson I, Lundberg A, Arundt A. Differences in ankle-joint complex during the stance phase of walking as measured by superficial and bone-anchored markers. Foot Ankle Int. 2002;23:856–63.PubMedGoogle Scholar
- Nester C, Jones RK, Liu A, Howard D, Lundberg A, Arndt A, et al. Foot kinematics during walking measured using bone and surface mounted markers. J Biomech. 2007;40:3412–23.View ArticlePubMedGoogle Scholar
- Arndt A, Westblad P, Winson I, Hashimoto T, Lundberg A. Ankle and subtalar kinematics measured with intracortical pins during the stance phase of walking. Foot Ankle Int. 2004;25:357–63.PubMedGoogle Scholar
- Lundgren P, Nester C, Liu A, Arndt A, Jones R, Stacoff A, et al. Invasive in vivo measurement of rear-, mid- and forefoot motion during walking. Gait Posture. 2008;28:93–100.View ArticlePubMedGoogle Scholar
- You BM, Anderst W, Tashman S. In vivo measurement of 3-D skeletal kinematics from sequences of biplane radiographs: application to knee kinematics. IEEE Trans Med Imaging. 2001;20:514–25.View ArticlePubMedGoogle Scholar
- Massimini DF, Warner JJ, Li G. Non-invasive determination of coupled motion of the scapula and humerus – an in vitro validation. J Biomech. 2011;44:408–12.View ArticlePubMedGoogle Scholar
- Wang S, Passias P, Li G, Wood K. Measurement of vertebral kinematics using noninvasive image matching method-validation and application. Spine. 2008;33:E355–61.View ArticlePubMedGoogle Scholar
- Baier DB, Gatesy SM, Dial KP. Three-dimensional, high-resolution skeletal kinematics of the avian wing and shoulder during ascending flapping flight and uphill flap-running. PLoS One. 2013;8:e63982.View ArticlePubMedPubMed CentralGoogle Scholar
- Brainerd EL, Baier DB, Gatesy SM, Hedrick TL, Metzger KA, Gilbert SL, et al. X-ray reconstruction of moving morphology (XROMM): precision, accuracy and applications in comparative biomechanics research. J Exp Zool. 2010;311A:262–79.Google Scholar
- Yamaguchi S, Sasho T, Kato H, Kuroyanagi Y, Banks SA. Ankle and subtalar kinematics during dorsiflexion-plantarflexion activities. Foot Ankle Int. 2009;30:361–6.View ArticlePubMedGoogle Scholar
- Fukano M, Kuroyanagi Y, Fukubayashi T, Banks S. Three-dimensional kinematics of the talocrural and subtalar joints during drop landing. J Appl Biomech. 2014;30:160–5.View ArticlePubMedGoogle Scholar
- Wan L, de Asla RJ, Rubash HE, Li G. Determination of in-vivo articular cartilage contact areas of human talocrural joint under weightbearing conditions. OsteoArth Cartil. 2006;14:1294–301.View ArticleGoogle Scholar
- de Asla RJ, Wan J, Rubash HE, Li G. Six DOF in vivo kinematics of the ankle joint complex: application of a combined dual-orthogonal fluoroscopic and magnetic resonance imaging technique. J Orthop Res. 2006;24:1019–27.View ArticlePubMedGoogle Scholar
- Kozanek M, Rubash HE, Li G, de Asla RJ. Effect of post-traumatic tibiotalar osteoarthritis on kinematics of the ankle joint complex. Foot Ankle Int. 2009;30:734–40.View ArticlePubMedGoogle Scholar
- Baka N, Kaptein BL, Giphart JE, Staring M, de Bruijne M, Lelieveldt BP, et al. Evaluation of automated statistical shape model based knee kinematics from biplane fluoroscopy. J Biomech. 2014;47:122–9.View ArticlePubMedGoogle Scholar
- Li G, Van de Velde SK, Bingham JT. Validation of a non-invasive fluoroscopic imaging technique for the measurement of dynamic knee joint motion. J Biomech. 2008;41:1616–22.View ArticlePubMedGoogle Scholar
- Tsai TY, Lu TW, Chen CM, Kuo MY, Hsu HC. A volumetric model-based 2D to 3D registration method for measurement kinematics of natural knees with single-plane fluoroscopy. Med Phys. 2010;37:1274–84.Google Scholar
- Tsai TY, Li JS, Wang S, Lin H, Malchau H, Li G, et al. A novel dual fluoroscopic imaging method for determination of THA kinematics: in-vitro and in-vivo study. J Biomech. 2013;46:1300–4.View ArticlePubMedGoogle Scholar
- Zhu Z, Massimini DF, Wang G, Warner JJ, Li G. The accuracy and repeatability of an automatic 2D-3D fluoroscopic image-model registration technique for determining shoulder joint kinematics. Med Eng Phys. 2012;34:1303–9.View ArticlePubMedGoogle Scholar
- Lin CC, Lu TW, Shih TF, Tsai TY, Wang TM, Hsu SJ. Intervertebral anticollision constraints improve out-of-plane translation accuracy of a single-plane fluoroscopy-to-CT registration method for measuring spinal motion. Med Phys. 2013;40:031912.View ArticlePubMedGoogle Scholar
- Campbell KL, Wilson KJ, LaPrade RF, Clanton TO. Normative rearfoot motion during barefoot and shod walking using biplane fluoroscopy. Knee Surg Sports Traumatok Arthrosc. 2014. doi:https://doi.org/10.1007/s00167-014-3084-4.Google Scholar
- Nigg BM, Herzog W. Biomechanics of the Musculo-skeletal System. 2nd ed. Chichester: John Wiley & Sons; 2007.Google Scholar
- Gutekunst DJ, Liu L, Ju T, Prior FW, Sinacore DR. Reliability of clinically relevant 3D foot bone angles from quantitative computed tomography. J Foot Ankle Res. 2013;6:38.View ArticlePubMedPubMed CentralGoogle Scholar
- Brunelli R. Template matching techniques in computer vision: theory and practice. New York: Wiley; 2009.View ArticleGoogle Scholar
- Peeters K, Natsakis T, Burg J, Spaepen P, Jonkers I, Dereymaeker G, et al. An in vitro approach to the evaluation of foot-ankle kinematics: Performance evaluation of a custom-built gait simulator. J Engineering Medicine. 2013;227:955–67.View ArticleGoogle Scholar
- Whittaker EC, Aubin PM, Ledoux WR. Foot bone kinematics as measured in a cadaveric robotic gait simulator. Gait Posture. 2010;33:645–60.View ArticleGoogle Scholar
- Nester CJ, Liu AM, Ward E, Howard D, Cocheba J, Derrick T, et al. In vitro study of foot kinematics using a dynamic walking cadaver model. J Biomech. 2007;40:1927–37.View ArticlePubMedGoogle Scholar
- Miranda DL, Schwartz JB, Loomis AC, Brainerd EL, Fleming BC, Crisco JJ. Static and Dynamic Error of a Biplanar Videoradiography System Using Marker-Based and Markerless Tracking Techniques. J Biomech Engineering. 2011;133:1–8.View ArticleGoogle Scholar
- Haase D, Denzler J. 2D and 3D analysis of animal locomotion from biplanar X-ray videos using augmented active appearance models. EURASIP J Image Video Process. 2013;2013:45.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.