PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
J Biomech. Author manuscript; available in PMC 2010 August 25.
Published in final edited form as:
PMCID: PMC2727693
NIHMSID: NIHMS119208

Using Two Palpable Measurements Improves the Subject Specific Femoral Modeling

Abstract

Subject-specific musculoskeletal models are essential to biomedical research and clinical applications, such as customized joint replacement, computer-aided surgical planning, gait analysis and automated segmentation. Generating these models from CT or MRI is time and resource intensive, requiring special skills. Therefore, in many studies individual bone models are approximated by scaling a generic template. Thus, the primary goal of this study was to determine a set of clinically available parameters (palpable measures and demographic data) that could improve the prediction of femoral dimensions, as compared to predicting these variables using uniform scaling based on palpable length. Similar to previous non-homogenous anthropometric scaling methods, the non-homogenous scaling method proposed in this study improved the prediction over uniform scaling of five key femoral measures. Homogenous scaling forces all dimensions of an object to be scaled equally, whereas non-homogenous scaling allows the dimensions to be scaled independently. The largest improvement was in femoral depth, where the coefficient of determination (r2) improved from 0.22 (homogenous) to 0.60 (non-homogeneous). In general, the major advantage of this non-homogenous scaling method is its ability to support the accurate and rapid generation of subject specific femoral models since all parameters can be collected clinically, without imaging or invasive methods.

Keywords: scaling, magnetic resonance imaging, femur

Introduction

Subject-specific musculoskeletal models are essential to biomedical research and clinical applications, such as customized joint replacement (Kessler et al., 2008), computer-aided surgery (Radermacher et al., 1998), gait analysis (Kepple et al., 1998) and automated image segmentation (Yao and Russell, 2003). For example, reliability and wear of implants are dependent upon the stress distribution in the bone (Kerner et al., 1999), which is primarily determined by implant properties (e.g. material, shape and component dimensions). The today’s longer life expectancies requires improving the durability and the reliability of joint implants, especially for younger patients with greater physical activity (Zanetti et al., 2005). Computer-aided orthopaedic surgery requires a geometrical model of the patient’s bone that can reproduce at least its basic morphological aspects and dimensions. In gait analysis, a realistic geometrical model is the first step towards creating an accurate musculoskeletal model (Kepple et al., 1998).

Accurate subject-specific models can be created from high-resolution in vivo imaging modalities, such as CT or MRI, but the creation of these models is expensive, time consuming and requires special skills. Thus, many studies have approximated an individual’s bone shape by scaling a generic template (Adam et al., 2002; Audenino et al., 1996; Bert, 1996; Brand et al., 1982; Lew and Lewis, 1977). The femur has been a particular focus of scaling studies due to its involvement in both the hip and knee joints. Scaling has usually been uniform and isotropic, requiring only a single scale factor, which is likely less accurate than multiple scaling factors. In more advanced studies the bone shape has been deformed by scaling with multiple factors (non-homogenous scaling). Lew et al. (1980) and Brand (1982) used three scale factors to represent femoral size along three orthogonal directions. Similarly, Sommer III et al. (1982) used an affine transformation, which allowed independent scaling along three orthogonal directions. This method was used by Kepple et al. (1998) to generate a 3D database of lower extremity musculotendon insertions and origins. More recent studies have used both x-ray (Gunay et al., 2007) and ultrasound (Rajamani et al., 2007) to generate a sparse subject-specific model. A higher resolution model was then generated by warping a generic high-resolution model to the low-resolution subject-specific model. The shortcomings of these non-homogenous methods are their dependence on internal landmark coordinates, requiring in vivo imaging.

Thus, the primary goal of this study was to determine a set of clinically available parameters (palpable measures and demographic data) that could accurately predict five key anatomic distances on the femur. “Accurately” was define as a coefficient of determination (r2) greater than 0.5.

Materials and method

The study included 76 femurs (52 skeletal femur specimens (cadavers) and 24 volunteers). All femurs were skeletally mature and free from known orthopaedic problems. Demographic data (height, weight, gender and age) were documented for each subject (172.8±8.8cm, 59.7±15.16kg, 42F/34M, 42.2±16.0 years).

Landmark coordinates from the 52 specimens were obtained from an existing database (Kepple et al, 1998). These coordinates were acquired using the Ascension Technologies tracking device (Burlington, VT), having a root-mean-square (RMS) calibration precision of 2 mm. Landmark coordinates from the 24 healthy volunteers’ femurs were obtained from geometric models generated using Geomagic Studio (Research Triangle Park, NC), based on T1-weighted spin echo MR images, segmented in MIPAV (National Institutes of Health, Bethesda, MD). All images had a 0.94 mm in-plane resolution. A slice thickness of 4 mm was used to acquire femoral volumes at the knee and hip where the anatomical landmarks were defined (Figure 1). A slice thickness of 15 mm was used for imaging the femoral shaft. Prior to extracting the anatomical landmarks from geometric models, the models were aligned based on their principal axes.

Figure 1
Key anatomical landmarks (red stars with red lower case labels) that were used to define key anatomic distances and palpable measures, which are listed in all caps. The palpable measures were palpable length, defined as the distance from the most lateral ...

Seven key measurements were calculated for each femoral dataset. The first two were palpable measures: palpable length and distal femoral width (Figure 1). The next five were key anatomic distances: hip length, shaft length, distal length, femoral depth, and total length (Figure 1). These five measures required in vivo imaging and were not attainable through palpation.

A first analysis was conducted to determine if uniform scaling using only palpable length was sufficient for creating subject-specific models from a normative femoral template (a model generated from one or more femurs that were free from any known defects or abnormalities). Four regression equations (linear with zero intercept, linear with y-intercept, non-linear semilog and non-linear second order polynomial) were calculated for each femoral measure using the palpable length as the independent variable. The methods were compared using the RMS error between the original value and the predicted value for each equation.

A second analysis was used to determine the improvement (over the first analysis) in predicting five of the seven key femoral measures using distal femoral width and subject demographic data (height, weight, and age) in conjunction with palpable length. These measures were added because they are clinical measures (i.e., measures that could be ascertained in a physical exam or demographic information). Stepwise linear regression was used to automatically select a set of parameters that constructed the most accurate predicting equations. The coefficient of determination (r2) determined the quality of fit.

Results

The standard deviation for the femoral measurements ranged from 0.5 to 2.9 cm (Table 1). The standard deviation for each measurement was weakly correlated with its average value (r=0.45, r= regression coefficient). All seven measurements were determined to be normally distributed with Lilliefors significance levels below 5%.

Table 1
Femoral measurements (n = 76)

Correlations between each of the six dependent measures and the palpable length ranged from 0.47 to 0.97 (Table 2). The quality of the fit was the same among the four regression equations, except for the hip length, which had a slightly higher RMS error for the polynomial function (Table 3).

Table 2
Correlation coefficients (r) among seven key lengths
Table 3
Comparison of the four fitting techniques for the first analysis (palpable length only)

Uniform scaling based on palpable length could not predict four of the six femoral measures (r2 < 0.5, Table 4). This, along with the large fluctuations in distal femoral width, depth, and hip length relative to palpable length (Figure 2), suggests that at least two independent variables are needed for accurate femoral scaling. Using additional clinical measures, along with the palpable length, improves the regression equations for each dependent variable and increases the coefficient of determination to above 0.5 in all cases (percent increase ranged from 2.1% to 75.0%, Table 4). Of all additional measures, distal femoral width has the greatest impact on improving prediction error and was the largest contributor to the prediction of distal length and depth (i.e., has the associated r2 value, Table 5). Step-wise regression analysis demonstrated r2 values, ranging from 0.41 to 1.00, with significance predictability (P<0.05, Table 5). The regression for the female and male populations did not improve the r2 values.

Figure 2
The data distribution is tight for total length and shaft length. The fit for each measurement versus palpable length is shown in a separate graph. The equations listed on each graph are a linear fit with zero-intercept, linear fit with y-intercept, the ...
Table 4
Comparison of the coefficient of determination (r2) between single regression with palpable length and stepwise multiple regression.
Table 5
Stepwise Linear Regression analysis

Discussion

This study advances previous work by using multiple scaling factors without requiring invasive measures or advanced high-resolution in vivo imaging. Similar to the non-homogenous anthropometric scaling method (Lewis et al., 1980), the non-homogenous scaling method proposed in this study improved the predictability over homogenous scaling using only palpable length and was able to predict 5 key femoral measures with a coefficient of determination (r2) greater than 0.5. Regression analysis demonstrated that using additional clinically available measurements (distal femoral width, subject height, weight, and age) improved the coefficient of determination from 2.1% to 75.0%. In general, the major advantage of this non-homogenous scaling method is its ability to be implemented quickly and accurately, since all parameters can be collected without imaging or invasive methods. Although, the homogenous method did not perform as well as the non-homogenous, the largest error was still under 12% (depth, average RMS error = 9.3mm). In studies where this is an acceptable level of error, uniform scaling based on palpable length could be implemented, providing a simpler alternative to the non-homogenous method.

Numerous scaling methodologies (Lew and Lewis, 1977; Lewis et al., 1980; Sommer, III et al., 1982) have been developed in an attempt to create subject-specific musculoskeletal models from generic femoral models. These studies focused on mathematical development of scaling techniques and used a minimal set of specimens (<3). Kepple et al (1998), using a larger specimen base (52 cadaver sets), demonstrated that uniform scaling based on subject height provided a RMS error of 6.6 mm for predicting femoral landmark locations. Yet, Lewis et al (1980) demonstrated that non-homogenous scaling reduced RMS errors in identifying anatomical landmark locations, particularly when a large number of points were used in the scaling. Since only a minimal set of anatomical landmarks are available through external palpation, the current work demonstrated that the most efficient application of non-homogeneous affine scaling was attainable through the use of three points, with RMS error of 5.0 mm (averaged over the six distance measures), which is less than the RMS error for individual points (6.5mm) reported by Kepple et al (1998).

The current study demonstrates a simple method to improve scaling by adding a single factor (distal femoral width) that can be easily measured through palpation. The strong correlations between palpable length and two measures of length (shaft length and total length) imply that the superior-inferior direction of the individual femoral sections change size proportionally to the palpable length. The lower correlation between distal femoral width and depth indicates that the two orthogonal dimensions do not. Thus, the use of uniform scaling based on generic models provides less accurate information than the use of multiple scaling factors.

The primary limitations of this study are two-fold. This study was limited to the femur. Yet, by simple extension, the same analysis could be applied to the entire lower extremity. The accuracy of predicting femoral measures is based on visualizing anatomical landmarks directly from a cadaver bone or a MR image. In practice, these landmarks would typically be found through palpation, which would likely introduce higher imprecision.

Acknowledgements

The Authors wish to thank Tom Kepple for providing data for this study. We also thank Bonnie Damaska, Susan O’Connell, Jere McLucas, and the Diagnostic Radiology Department at the National Institutes of Health for their support and research time. Any opinions, findings, and conclusions or recommendations expressed in this material are those of the author and do not necessarily reflect the views of the National Institutes of Health or the US Public Health Service. This research was supported in part by the Intramural Research Program of the NIH, CC and NICHD. The NIH did not have any involvement in the study design, in the collection, analysis and interpretation of data; in the writing of the manuscript; nor in the decision to submit the manuscript for publication.

Footnotes

Publisher's Disclaimer: This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

Conflict of interest statement No authors have any conflict of interest to disclose

Reference List

  • Adam F, Hammer DS, Pape D, Kohn D. Femoral anatomy, computed tomography and computer-aided design of prosthetic implants. Arch.Orthop.Trauma Surg. 2002;122:262–268. [PubMed]
  • Audenino AL, Zanetti EM, Calderale PM. Radiograph processing for quantitative assessment of bone remodeling. Med.Eng Phys. 1996;18:382–389. [PubMed]
  • Bert JM. Custom total hip arthroplasty. J.Arthroplasty. 1996;11:905–915. [PubMed]
  • Brand RA, Crowninshield RD, Wittstock CE, Pedersen DR, Clark CR, van Krieken FM. A model of lower extremity muscular anatomy. J.Biomech.Eng. 1982;104:304–310. [PubMed]
  • Gunay M, Shim MB, Shimada K. Cost- and time-effective three-dimensional bone-shape reconstruction from X-ray images. Int.J.Med.Robot. 2007;3:323–335. [PubMed]
  • Kepple TM, Sommer HJ, III, Lohmann SK, Stanhope SJ. A three-dimensional musculoskeletal database for the lower extremities. J.Biomech. 1998;31:77–80. [PubMed]
  • Kerner J, Huiskes R, van Lenthe GH, Weinans H, van RB, Engh CA, Amis AA. Correlation between pre-operative periprosthetic bone density and post-operative bone loss in THA can be explained by strain-adaptive remodeling. J.Biomech. 1999;32:695–703. [PubMed]
  • Kessler O, Patil S, Wirth S, Mayr E, Colwell CW, Jr., D’Lima DD. Bony impingement affects range of motion after total hip arthroplasty: A subject-specific approach. J.Orthop.Res. 2008;26:443–452. [PubMed]
  • Lew WD, Lewis JL. An anthropometric scaling method with application to the knee joint. J.Biomech. 1977;10:171–181. [PubMed]
  • Lewis JL, Lew WD, Zimmerman JR. A nonhomogeneous anthropometric scaling method based on finite element principles. J.Biomech. 1980;13:815–824. [PubMed]
  • Radermacher K, Portheine F, Anton M, Zimolong A, Kaspers G, Rau G, Staudte HW. Computer assisted orthopaedic surgery with image based individual templates. Clin.Orthop.Relat Res. 1998:28–38. [PubMed]
  • Rajamani KT, Styner MA, Talib H, Zheng G, Nolte LP, Ballester M. A. Gonzalez. Statistical deformable bone models for robust 3D surface extrapolation from sparse data. Med.Image Anal. 2007;11:99–109. [PubMed]
  • Sommer HJ, III, Miller NR, Pijanowski GJ. Three-dimensional osteometric scaling and normative modeling of skeletal segments. J.Biomech. 1982;15:171–180. [PubMed]
  • Yao J, Russell T. Non-rigid registration and correspondence finding in medical image analysis using multiple-layer flexible mesh template matching. International journal of pattern recognition and artificial intelligence. 2003;17:1145–1165.
  • Zanetti EM, Crupi V, Bignardi C, Calderale PM. Radiograph-based femur morphing method. Med.Biol.Eng Comput. 2005;43:181–188. [PubMed]