This article has Open Peer Review reports available.
Identification of hip fracture patients from radiographs using Fourier analysis of the trabecular structure: a cross-sectional study
© Gregory et al; licensee BioMed Central Ltd. 2004
Received: 30 April 2004
Accepted: 06 October 2004
Published: 06 October 2004
This study presents an analysis of trabecular bone structure in standard radiographs using Fourier transforms and principal components analysis (PCA) to identify contributions to hip fracture risk.
Radiographs were obtained from 26 hip fracture patients and 24 controls. They were digitised and five regions of interest (ROI) were identified from the femoral head and neck for analysis. The power spectrum was obtained from the Fourier transform of each region and three profiles were produced; a circular profile and profiles parallel and perpendicular to the preferred orientation of the trabeculae. PCA was used to generate a score from each profile, which we hypothesised could be used to discriminate between the fracture and control groups. The fractal dimension was also calculated for comparison. The area under the receiver operating characteristic curve (A z ) discriminating the hip fracture cases from controls was calculated for each analysis.
Texture analysis of standard radiographs using the fast Fourier transform yielded variables that were significantly associated with fracture and not significantly correlated with age, body mass index or femoral neck bone mineral density. The anisotropy of the trabecular structure was important; both the perpendicular and circular profiles were significantly better than the parallel-profile (P < 0.05). No significant differences resulted from using the various ROI within the proximal femur. For the best three groupings of profile (circular, parallel or perpendicular), method (PCA or fractal) and ROI (A z = 0.84 – 0.93), there were no significant correlations with femoral neck bone mineral density, age, or body mass index. PCA analysis was found to perform better than fractal analysis (P = 0.019).
Both PCA and fractal analysis of the FFT data could discriminate successfully between the fracture and control groups, although PCA was significantly stronger than fractal dimension. This method appears to provide a powerful tool for the assessment of bone structure in vivo with advantages over standard fractal methods.
The NIH Consensus Statement defines Osteoporosis as "a skeletal disorder characterised by compromised bone strength predisposing to an increased risk of fracture" . Bone strength was defined as "the integration of two main features: bone density and bone quality". Currently, clinical diagnosis is based solely on bone mineral density (BMD) in accordance with the World Health Organisation guidelines . Previous studies, however, have found that trabecular bone structure also plays a significant role in determining bone strength [3–5] with BMD explaining only 60 to 80 % of the variability in mechanical resistance .
Trabecular bone structure is visible on standard pelvic radiographs and many attempts have been made to quantify the quality of the structure and assess its relationship to osteoporosis and BMD. These range from visual scoring systems, such as the Singh index , through to sophisticated computerised methods based on fractals [8–10] and other image processing methods [11–13]. A review of the literature suggests that fractal analysis has been a method of choice in recent years for the analysis of trabecular bone structure in CT scans [14, 15], MRI , histology  and radiographs [18–21], although it has not been established categorically that it is preferable to other methods of texture analysis [22, 23]. By reducing all the information in the image to one descriptor, the fractal dimension , a large part of the information is lost. The Fourier transform of an image expresses the information in the image in terms of spatial frequencies rather than distances. Various methods can be applied to extract information from the Fourier transform , including the fractal dimension . However such methods have not been fully exploited for analysing bone structure [8, 26–30].
In this study we investigate the use of Fourier transforms and Principal Components Analysis to generate a mathematical model of the data which can be used to help classify individuals according to the presence or absence of a hip fracture. Principal component analysis (PCA)  is a data reduction technique that has been applied in many fields of study, including investigation of gene expression , development of an electronic nose  and tracing of the evolutionary changes in fish morphometry . It describes data in terms of a small number of orthogonal, linearly independent components which contain the majority of the information. PCA has no preconditions, such as relying on the data to fit a normal or fractal distribution, but builds a mathematical model based on the correlations present in the data. An eigenanalysis of the correlation or covariance matrix is used to perform PCA. The resulting components are then selected in order of the amount of variance they account for, enabling an efficient mapping of the data. As the first few components account for the vast majority of the variance in the original data, they can be selected for analysis whilst the remainder are discarded as 'noise'. In this way, the number of variables can be greatly reduced whilst maintaining the information present in the original data. In this pilot study we used these methods to investigate the similarities and differences between trabecular bone structure in fracture and control groups using standard radiographs of the proximal femur.
A set of digitised standard pelvic radiographs was available from a previous investigation into the morphology of the proximal femur . These radiographs were taken from an earlier study , that had examined three groups (osteoporotic, osteoarthritic and control) of age matched, postmenopausal women (30 subjects per group). Subjects with osteoarthritis were excluded from the present study. All patients had undergone a scan of the unfractured hip by dual-energy x-ray absorptiometry (DXA) using a Norland XR-26 scanner (CooperSurgical Inc, Trumbull, CT). The controls had had their left hip scanned. All patients and controls had had a pelvic antero-posterior radiograph recorded within a year of the DXA scan. We used those radiographs and the femoral neck BMD (Neck-BMD) data in the current study. A data set of 50 digitised radiographs was available comprising 26 hip fracture patients (HIP) and 24 controls (CNT). The radiographs were digitised, using a Howtek MultiRAD 850 scanner (Howtek, Hudson, New Hampshire) at a resolution of 584 dpi (44 μm per pixel) and a depth of 12 bits. The age, height and weight of each subject were also recorded.
Each ROI was 256 × 256 pixels (11.3 mm square), to enable use of the fast Fourier transform, and were selected as follows. The upper region of the head (UH) lies on the upper part of the principal compressive trabeculae, the central region of the head (CH) is at the intersection of the principal compressive and tensile trabeculae, the upper region of the neck (UN) lies on the principal tensile trabeculae, the lower region of the neck (LN) is at the base of the principal compressive trabeculae and finally Ward's triangle (WA) which lies between these structures. The points and regions were identified using a macro written for Image Pro Plus software (version 184.108.40.206, Media Cybernetics, Silver Spring, Maryland). The femoral head was described by a best-fit circle, calculated from a series of manually marked points around the outline of the femoral head. Between 15 and 20 evenly spaced points were used to describe the outline, depending on the size of the head. The radius and centre (marked as F in Figure 1) of the femoral head were then taken from this circle. The narrowest part of the neck (neck-width) was determined using two automatic edge traces, marking the upper and lower outlines of the femoral neck. The first point and the direction for each trace were marked manually; the edge of the neck could then be identified automatically by the software. The neck width (A – E in Figure 1) was calculated by finding the smallest Euclidean distance between the traces. The centre of the neck was located at the mid-point of this line (point C) and the axis of the femoral neck was taken to be a line perpendicular to this through the centre of the neck (dashed line). The top right corner of the WA region was located at the midpoint of the neck width (point C). Points B and D were placed 25% and 75% of the way along the neck width and used as the midpoints of the UN and LN regions respectively. Point F, the centre of the femoral head marked the centre of the CH region and point G, the centre of the base of the UH region. Point G was placed one half of the femoral head radius above point F, at a 45-degree angle to the neck width (A-E).
Analysis was performed using Matlab software (version 6.1.0, MathWorks Inc, Natick, Massachusetts). A fast Fourier transform was generated for each ROI and three profiles were generated using data from the power spectrum. Firstly a global or circular profile (CircP) was generated, composed of the magnitude at each spatial frequency averaged across all angles, resulting in a profile with 128 data points. To create this profile, each pixel in the Fourier transform was assigned to the integer spatial frequency that most closely matched its' distance from the zero'th component.
Principal component analysis
Fractal analysis was performed on each profile using a method similar to the Fourier transform technique described by Majumdar et al . The average power spectrum of the circular profile was plotted on a log-log scale, three approximately linear regions were defined and the gradient (slope) of a straight line fitted to each region was found; slopeA, a 'coarse' slope, where the log of the spatial frequency is less than or equal to 1.0, slopeB a 'medium' slope, where the log of the spatial frequency lies between 1.0 and 1.75 and slopeC, a 'fine' slope where the log of the spatial frequency is above 1.75. The fractal dimension was calculated for each slope using the formula suggested by Majumdar et al 
Stepwise discriminant analysis was used to select principal components that could be combined to build a linear classifier. If the stepwise procedure failed to select any components, the most accurate of the individual components was chosen. The same procedure was used to discriminate between the groups using the fractal dimension. Measurement of the area under the ROC curve was used to compare the classifiers built using the discriminant analysis . A three way ANOVA was applied in order to determine whether there were significant differences between the performance of classifiers depending on the type of analysis, the profile used or the region analysed. Pearson product moment correlation was applied to examine the relationship with age, BMI and Neck BMD for the strongest classifiers. A one-way ANOVA was used to test for significant differences in the performance of the slopes from each spatial frequency band used in the fractal analysis. T-tests, correlation and ANOVA were performed using SigmaStat (version 2.03, SPSS Science, Chicago). Principal component analysis, discriminant analysis, and measurement of the area under the ROC curve were calculated using SPSS (version 10 SPSS Science, Chicago).
Summary of anthropometric variables for the fracture and controls groups. Mean and standard deviation (SD) of the age, height, weight, BMI and BMD of the fracture and control groups. P values were obtained from a two-tailed t-test.
Control Group (n = 24)
Fracture Group (n = 26)
Body Mass Index, kg/m2
Femoral neck BMD (g cm -2 )
Classification accuracy for each region-profile combination. Area under the ROC curve for principal component analysis of each profile by region of the femoral neck. Analysis using three-way ANOVA found that the area under the ROC curve was significantly higher in the perpendicular profile than in the parallel profile. (P < 0.05)
Circular (95% CI)
Parallel (95% CI)
Perpendicular (95% CI)
0.76 (0.63 – 0.89)
0.57 (0.41 – 0.73)
0.84 (0.73 – 0.95)
0.59 (0.43 – 0.75)
0.56 (0.40 – 0.73)
0.84 (0.72 – 0.95)
0.72 (0.58 – 0.86)
0.72 (0.57 – 0.86)
0.67 (0.52 – 0.82)
0.74 (0.61 – 0.88)
0.61 (0.45 – 0.76)
0.71 (0.56 – 0.86)
0.71 (0.56 – 0.85)
0.55 (0.39 – 0.71)
0.93 (0.87 – 1.00)
The best five classifiers: Area under the curve and correlation with BMD, age and BMI. Area under the ROC curve (A z ) for each of the best 5 classifiers and the correlation with age Rage, femoral neck BMD (RBMD) and body mass index (RBMI) and associated significance values (P).
A z (95% CI)
0.93 (0.87 – 1.00)
0.84 (0.73 – 0.95)
0.84 (0.72 – 0.95)
0.76 (0.63 – 0.89)
0.75 (0.61 – 0.89)
Table 3 also shows the correlations between the top five classifiers with age, BMI and Neck-BMD. No significant correlations were found between any of these classifiers and either age or BMI and, for the top three, there was also no significant correlation with Neck-BMD (P > 0.05). The fifth placed classifier, fractal analysis of the parallel profile in the upper neck region, was the only one significantly associated with Neck-BMD (P = 0.034).
A three-way analysis of variance was used to examine differences in performance due to the region, profile or type of analysis used. It showed that overall PCA analysis performed significantly better than fractal analysis (P = 0.019) and that analysis of both the perpendicular and circular profiles performed significantly better than the parallel profile (P = 0.003 and 0.011 respectively). No significant differences were found between the different regions of the femoral neck (P = 0.241) (despite the apparently large differences in A z ). The power of this test was 0.69, 0.97 and 0.15 for the investigation of differences due to the method of analysis, type of profile used and the region analysed respectively.
Comparing slopeA, slopeB and slopeC. The average and standard deviation of the area under the ROC curve (A z ) are presented for each of the slopes used in the fractal analysis for all regions of interest. A significant difference was found between slopeA and slopeC in the circular profile, however when all the profiles were compared, no significant differences were found.
Discussion and conclusions
In these short series, this study found that texture analysis of standard radiographs using the fast Fourier transform can yield variables that are significantly associated with fracture but not significantly correlated with age, body mass index or Neck-BMD. Both PCA and fractal analysis of the FFT data could be used to discriminate successfully between the groups, although overall PCA was significantly stronger than fractal dimension. The best results from this study were not significantly correlated with femoral neck-BMD, age or BMI, indicating their potential for use as an independent predictor of fracture. The radiographic appearance of bone is known to be affected by factors including the size of the patient. As there was no significant difference in the BMI of the fracture and control groups, it is unlikely that this has influenced the results, however it is an issue that will need addressing in future studies.
Previous studies using non-fractal analysis of the Fourier power spectrum have focussed on images of the spine or wrist, where the alignment of trabeculae is generally orthogonal [28–30]. In such images, analysis of trabecular orientation can be performed by examining the vertical and horizontal sectors as the trabeculae lie predominantly in these directions. The trabecular structure of the femur is more complicated as the trabeculae are aligned in arcs, so the preferred orientation changes throughout the proximal femur. Analysis parallel to the preferred orientation of the trabeculae was significantly poorer than analysis using either the perpendicular or circular profiles (P < 0.05). Analysis in the perpendicular direction was strongest overall, although it was not significantly better than the circular profile. This accords with the increasingly anisotropic nature of trabecular bone with aging; bone loss is not evenly distributed but is lost primarily at angles perpendicular and oblique to the preferred orientation of the trabeculae . This loss heightens the risk of fracture, especially if the impact is from the side, as expected from a typical fall from standing height, as there are fewer trabeculae orientated in this direction to absorb the force of impact.
In summary, this paper presents a new method for analysing the structure of trabecular bone from standard radiographs. It demonstrates that the Fourier transform can be used to describe structural information in images which may be related to fracture, independently of BMD. This study is limited by the small size of the data set and further analysis is needed to validate these findings. This should be performed on a similar series of radiographs, consisting of fracture and control subjects scanned at the same resolution. The methods from this study could then be applied directly to this group (without recalculating the PCA) to evaluate whether they were generally applicable. However the success of both this and our previous study, using similar techniques to analyse histological sections, indicates that this may be an effective method with clinical utility for describing bone quality statistically in terms of structural parameters.
We thank The PPP Foundation and the Arthritis Research Campaign for funding this study and the MRC for a Senior Fellowship for RMA. We are grateful to Mr G. Turner for expert technical assistance. AS is an Arthritis Research Campaign Postdoctoral Research Fellow.
- Osteoporosis Prevention, Diagnosis and Therapy. NIM Consens Statement. 2000, 17: 1-45.Google Scholar
- Group WHO Study: Assessment of fracture risk and its application to screening for postmenopausal osteoporosis. 1994, Geneva, World Health Organisation, WHO Technical Report Series 843:Google Scholar
- Veenland JF, Link TM, Konermann W, Meier N, Grashuis JL, Gelsema ES: Unravelling the role of structure and density in determining vertebral bone strength. Calcif Tissue Int. 1997, 61: 474-479. 10.1007/s002239900370.View ArticlePubMedGoogle Scholar
- Heaney R: Is there a role for Bone Quality in Fragility Fractures. Calcif Tissue Int. 1993, 53: S3-S6.View ArticlePubMedGoogle Scholar
- Snyder BD, Piazza S, Edwards WT, Hayes WC: Role of Trabecular Morphology in the Etiology of Age-Related Vertebral Fractures. Calcif Tissue Int. 1993, 53: S14-S22.View ArticlePubMedGoogle Scholar
- Cortet B, Marchandise X: Bone microarchitecture and mechanical resistance. Joint Bone Spine. 2001, 68: 297-305. 10.1016/S1297-319X(01)00283-4.View ArticleGoogle Scholar
- Singh M, Nagrath AR, Mani PS: Changes in Trabecular Pattern of the Upper End of the Femur as an Index of Osteoporosis. J Bone Joint Surg Am. 1970, 52-A: 457-467.Google Scholar
- Chung H, Chu C, Underweiser M, Wehrli FW: On the fractal nature of trabecular structure. Medical Physics. 1994, 21: 1535-1540. 10.1118/1.597263.View ArticlePubMedGoogle Scholar
- Lin JC, Grampp S, Link T, Kothari M, Newitt DC, Felsenberg D, Majumdar S: Fractal Analysis of Proximal Femur Radiographs: Correlation with Biomechanical Properties and Bone Mineral Density. Osteoporos Int. 1999, 9: 516-524.View ArticlePubMedGoogle Scholar
- Pornprasertsuk S, Ludlow JB, Webber RL, Tyndall DA, Sanhueza AI, Yamauchi M: Fractal Dimension Analysis of Weight-bearing Bones of Rats During Skeletal Unloading. Bone. 2001, 29: 180-184. 10.1016/S8756-3282(01)00493-8.View ArticlePubMedGoogle Scholar
- Geraets WGM, Van der Stelt PF, Lips P, Van Ginkel FC: The radiographic trabecular pattern of hips in patients with hip fractures and in elderly control subjects. Bone. 1998, 22: 165-173. 10.1016/S8756-3282(97)00233-0.View ArticlePubMedGoogle Scholar
- Lee RL, Dacre JE, Hart DJ, Spector TD: Femoral neck trabecular patterns predict osteoporotic fractures. Medical Physics. 2002, 29: 1391-1396. 10.1118/1.1481514.View ArticlePubMedGoogle Scholar
- Smyth PP, Adams JE, Whitehouse RW, Taylor CJ: Application of computer texture analysis to the Singh index. Br J Radiol. 1997, 70: 242-247.View ArticlePubMedGoogle Scholar
- Haidekker MA, Andresen R, Werner HJ: Relationship between structural parameters, bone mineral density and fracture load in lumbar vertebrae, based on high- resolution computed tomography, quantitative computed tomography and compression tests. Osteoporos Int. 1999, 9: 433-440. 10.1007/s001980050168.View ArticlePubMedGoogle Scholar
- Dougherty G, Henebry GM: Lacunarity analysis of spatial pattern in CT images of vertebral trabecular bone for assessing osteoporosis. Medical Engineering & Physics. 2002, 24: 129-138. 10.1016/S1350-4533(01)00106-0.View ArticleGoogle Scholar
- Link TM, Majumdar S, Augat P, Lin JC, Newitt D, Lu Y, Lane NE, Genant HK: In vivo high resolution MRI of the calcaneus: Differences in trabecular structure in osteoporosis patients. J Bone Miner Res. 1998, 13: 1175-1182.View ArticlePubMedGoogle Scholar
- Parkinson IH, Fazzalari NL: Methodological principles for fractal analysis of trabecular bone. Journal of Microscopy. 2000, 198: 134-142. 10.1046/j.1365-2818.2000.00684.x.View ArticlePubMedGoogle Scholar
- Benhamou CL, Poupon S, Lespessailles E, Loiseau S, Jennane R, Siroux V, Ohley W, Pothuaud L: Fractal analysis of radiographic trabecular bone texture and bone mineral density: Two complementary parameters related to osteoporotic fractures. J Bone Miner Res. 2001, 16: 697-704.View ArticlePubMedGoogle Scholar
- Southard TE, Southard KA: Detection of simulated osteoporosis in maxillae using radiographic texture analysis. IEEE Trans Biomed Eng. 1996, 43: 123-132. 10.1109/10.481981.View ArticlePubMedGoogle Scholar
- Caligiuri P, Giger ML, Favus M: Multifractal radiographic analysis of osteoporosis. Medical Physics. 1994, 21: 503-508. 10.1118/1.597390.View ArticlePubMedGoogle Scholar
- Lynch JA, Hawkes DJ, Buckland-Wright JC: Analysis of texture in macroradiographs of osteoarthritic knees using the fractal signature. Phys Med Biol. 1991, 36: 709-722. 10.1088/0031-9155/36/6/001.View ArticlePubMedGoogle Scholar
- Lee RL, Dacre JE, James MF: Image processing assessment of femoral osteopenia. Journal of Digital Imaging. 1997, 10: 218-221.View ArticlePubMedPubMed CentralGoogle Scholar
- Cortet B, Dubois P, Boutry N, Bourel P, Cotten A, Marchandise X: Image Analysis of the Distal Radius Trabecular Network Using Computed Tomography. Osteoporos Int. 1999, 9: 410-419. 10.1007/s001980050165.View ArticlePubMedGoogle Scholar
- Pentland AP: Fractal-Based Description of Natural Scenes. IEEE Transactions On Pattern Analysis and Machine Intelligence. 1984, PAMI-6: 661-674.View ArticleGoogle Scholar
- Sonka M, Hlavac V, Boyle RD: Image Processing, Analysis and Machine Vision. 1998, Pacific Grove, CA, PWS Publishing, 2ndGoogle Scholar
- Gregory JS, Junold RM, Undrill PE, Aspden RM: Analysis of trabecular bone structure using Fourier transforms and neural networks. IEEE Transactions On Technology In Biomedicine. 1999, 3: 289-294. 10.1109/4233.809173.View ArticleGoogle Scholar
- Caligiuri P, Giger ML, Favus MJ, Jia H, Doi K, Dixon LB: Computerized radiographic analysis of osteoporosis: preliminary evaluation. Radiology. 1993, 186: 471-474.View ArticlePubMedGoogle Scholar
- Wigderowitz CA, Abel EW, Rowley DI: Evaluation of cancellous structure in the distal radius using spectral analysis. Clin Orthop. 1997, 152-161.Google Scholar
- Wigderowitz CA, Paterson CR, Dashti H, McGurty D, Rowley DI: Prediction of bone strength from cancellous structure of the distal radius: Can we improve on DXA?. Osteoporos Int. 2000, 11: 840-846. 10.1007/s001980070042.View ArticlePubMedGoogle Scholar
- Buck AM, Price RI, Sweetman IM, Oxnard CE: An investigation of thoracic and lumbar cancellous vertebral architecture using power-spectral analysis of plain radiographs. Journal of Anatomy. 2002, 200: 445-456. 10.1046/j.1469-7580.2002.00051.x.View ArticlePubMedPubMed CentralGoogle Scholar
- Chatfield C, Collins AJ: Principal component analysis. Introduction to Multivariate Analysis. 1989, Chapman and Hall, 4: 57-81. 4Google Scholar
- Yeung KY, Ruzzo WL: Principal component analysis for clustering gene expression data. Bioinformatics. 2001, 17: 763-774. 10.1093/bioinformatics/17.9.763.View ArticlePubMedGoogle Scholar
- Boothe DDH, Arnold JW: Electronic nose analysis of volatile compounds from poultry meat samples, fresh and after refrigerated storage. Journal of the science of food and agriculture. 2002, 82: 315-322. 10.1002/jsfa.1036.View ArticleGoogle Scholar
- Walker JA, Bell MA: Net evolutionary trajectories of body shape evolution within a microgeographic radiation of threespine sticklebacks (Gasterosteus aculeatus). Journal Of Zoology. 2000, 252: 293-302. 10.1017/S0952836900000030.View ArticleGoogle Scholar
- Gregory JS, Testi D, Stewart A, Undrill PE, Reid DM, Aspden RM: A method for assessment of the shape of the proximal femur and its relationship to osteoporotic hip fracture. Osteoporos Int. 2004, 15: 5-11. 10.1007/s00198-003-1451-y.View ArticlePubMedGoogle Scholar
- Stewart A, Black A, Robins SP, Reid DM: Bone Density and Bone Turnover in Patients with Osteoarthritis and Osteoporosis. J Rheumatol. 1999, 26: 622-626.PubMedGoogle Scholar
- Crabtree N, Lunt M, Holt G, Kroger H, Burger H, Grazio S, Khaw KT, Lorenc RS, Nijs J, Stepan J, Falch JA, Miazgowski T, Raptou P, Pols H, Dequeker J, Havelka S, Hoszowski K, Jajic I, Czekalski S, Lyritis G, Silman AJ, Reeve J: Hip Geometry, Bone Mineral Distribution, and Bone Strength in European Men and Women: The EPOS Study. Bone. 2000, 27: 151-159. 10.1016/S8756-3282(00)00300-8.View ArticlePubMedGoogle Scholar
- Défossez H, Hall RM, Walker PG, Wroblewski BM, Siney PD, Purbach B: Determination of the trabecular bone direction from digitised radiographs. Medical Engineering & Physics. 2003, 25: 719-729. 10.1016/S1350-4533(03)00123-1.View ArticleGoogle Scholar
- Cattell RB: The scree test for a number of factors. Multivariate Behavioural Research. 1966, 1: 245-276.View ArticleGoogle Scholar
- Majumdar S, Lin T, Link TM, Millard J, Augat P, Ouyang X, Newitt D, Gould R, Kothari M, Genant HK: Fractal analysis of radiographs: Assessment of trabecular bone structure and prediction of elastic modulus and strength. Medical Physics. 1999, 26: 1330-1340. 10.1118/1.598628.View ArticlePubMedGoogle Scholar
- Metz CE: ROC Analysis in Radiologic Imaging. Investigative Radiology. 1986, 21: 720-733.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2342/4/4/prepub
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/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.