Automatic construction of an anatomical coordinate system for three-dimensional bone models of the lower extremities – Pelvis,femur, and tibia

Automated methods for constructing patient-specific anatomical coordinate systems (ACSs) for the pelvis, femur and tibia were developed based on the bony geometry of each, derived from computed tomography (CT). The methods used principal axes of inertia, principal component analysis (PCA), cross-sectional area, and spherical and ellipsoidal surface fitting to eliminate the influence of rater's bias on reference landmark selection. Automatic ACSs for the pelvis, femur, and tibia were successfully constructed on each 3D bone model using the developed algorithm. All constructions were performed within 30 s; furthermore, between- and within- rater errors were zero for a given CT-based 3D bone model, owing to the automated nature of the algorithm. ACSs recommended by the International Society of Biomechanics (ISB) were compared with the automatically constructed ACS, to evaluate the potential differences caused by the selection of the coordinate system. The pelvis ACSs constructed using the ISB-recommended system were tilted significantly more anteriorly than those constructed automatically (range, 9.6–18.8°). There were no significant differences between the two methods for the femur. For the tibia, significant differences were found in the direction of the anteroposterior axis; the anteroposterior axes identified by ISB were more external than those in the automatic ACS (range, 17.5–25.0°).     

Pair-wise vs group-wise registration in statistical shape model construction: representation of physiological and pathological variability of bony surface morphology

Statistical shape models (SSM) of bony surfaces have been widely proposed in orthopedics, especially for anatomical bone modeling, joint kinematic analysis, staging of morphological abnormality, and pre- and intra-operative shape reconstruction. In the SSM computation, reference shape selection, shape registration and point correspondence computation are fundamental aspects determining the quality (generality, specificity and compactness) of the SSM. Such procedures can be made critical by the presence of large morphological dissimilarities within the surfaces, not only because of anthropometrical variability but also mainly due to pathological abnormalities. In this work, we proposed a SW pipeline for SSM construction based on pair-wise (PW) shape registration, which requires the a-priori selection of the reference shape, and on a custom iterative point correspondence algorithm. We addressed large morphological deformations in five different bony surface sets, namely proximal femur, distal femur, patella, proximal fibula and proximal tibia, extracted from a retrospective patient dataset. The technique was compared to a method from the literature, based on group-wise (GW) shape registration. As a main finding, the proposed technique provided generalization and specificity median errors, for all the five bony regions, lower than 2?mm. The comparative analysis provided basically similar results. Particularly, for the distal femur that was the shape affected by the largest pathological deformations, the differences in generalization, specificity and compactness were lower than 0.5?mm, 0.5?mm, and 1%, respectively. We can argue the proposed pipeline, along with the robust correspondence algorithm, is able to compute high-quality SSM of bony shapes, even affected by large morphological variability.     

Articular-surface-based automatic anatomical coordinate systems for the knee bones

Increasing use of patient-specific surgical procedures in orthopaedics means that patient-specific anatomical coordinate systems (ACSs) need to be determined. For knee bones, automatic algorithms constructing ACSs exist and are assumed to be more reliable than manual methods, although both approaches are based on non-unique numerical reconstructions of true bone geometries. Furthermore, determining the best algorithms is difficult, as algorithms are evaluated on different datasets. Thus, in this study, we developed 3 algorithms, each with 3 variants, and compared them with 5 from the literature on a dataset comprising 24 lower-limb CT-scans. To evaluate algorithms’ sensitivity to the operator-dependent reconstruction procedure, the tibia, patella and femur of each CT-scan were each reconstructed once by three different operators.Our algorithms use principal inertia axis (PIA), cross-sectional area, surface normal orientations and curvature data to identify the bone region underneath articular surfaces (ASs). Then geometric primitives are fitted to ASs, and the ACSs are constructed from the geometric primitive points and/or axes. For each bone type, the algorithm displaying the least inter-operator variability is identified. The best femur algorithm fits a cylinder to posterior condyle ASs and a sphere to the femoral head, average axis deviations: 0.12°, position differences: 0.20 mm. The best patella algorithm identifies the AS PIAs, average axis deviations: 0.91°, position differences: 0.19 mm. The best tibia algorithm finds the ankle AS center and the 1st PIA of a layer around a plane fitted to condyle ASs, average axis deviations: 0.38°, position differences: 0.27 mm.     

Comparative evaluation of scapular and humeral coordinate systems based on biomedical images of the glenohumeral joint

Bio-imaging techniques represent a powerful tool for shoulder joint biomechanical analysis. However, the restricted field of view may prevent the acquisition of complete scapula and humerus bone models and hence limiting the applicability of standardized anatomical coordinate system (ACS) definitions. The aim of this study was to propose ACS definitions for both scapula and humerus which can be implemented when limited portions of the relevant bones are available. Magnetic resonance (MR) images of twenty right humeri and scapulae were acquired. The proposed ACSs were assessed in terms of (1) sensitivity to bone morphological variation, (2) intra – and inter – operator repeatability and (3) consistency with the anatomical cardinal directions. A comparison with alternative ACS definitions was also performed. Overall, our ACS scapular proposal and that presented in Kedgley and Dunning (2010) were found to be the least sensitive to the morphometric variability (mean angular absolute deviation lower than 8.3 deg) and they were characterized by a high intra – and inter – operator repeatability (mean angular absolute deviation lower than 1.5 deg). The humeral ACS proposal showed a morphometric variability similar to Amadi et al. (2009b) (mean angular absolute deviation lower than 8.3 deg) but a higher reproducibility. The scapular and humeral ACS mean angular deviation from the reference anatomical cardinal directions were smaller than 15 deg and 8.6 deg, respectively. The proposed scapular and humeral ACS definitions are therefore suitable to be applied when a limited portion of the glenohumeral joint is available as it may occur in standard shoulder clinical exams.     

Accuracy and precision of digital volume correlation in quantifying displacements and strains in trabecular bone

Strain measurement is an essential tool in the study of trabecular bone structure-function relationships. Digital volume correlation (DVC) is a measurement technique that quantifies strains throughout the interior of a specimen, rather than simply those on the surface. DVC relies on tracking the movement of microstructural features, and as such, the accuracy and precision of this technique may depend on trabecular structure. This study quantified displacement and strain measurement errors in six types of trabecular bone that spanned a wide range of volume fraction and trabecular architecture. Accuracy and precision were compared across bone type and also across three DVC methods. Both simulated and real displacement fields were analyzed using micro-computed tomography images of specimens from the bovine distal femur, bovine proximal tibia, rabbit distal femur, rabbit proximal tibia, rabbit vertebra, and human vertebra. Differences as large as three-fold in accuracy and precision of the displacements and strains were found among DVC methods and among bone types. The displacement precision and the strain accuracy and precision were correlated with measures of trabecular structure such as structural model index. These results demonstrate that the performance of the DVC technique can depend on trabecular structure. Across all bone types, the displacement and strain errors ranged 1.86-3.39 microm and 345-794 microepsilon, respectively. For specimens from the human vertebra and bovine distal femur, the measurement errors were approximately 20 times smaller than the yield strain. In these cases, DVC is a viable technique for measuring pre- and post-yield strains throughout trabecular bone specimens and the trabecular compartment of whole bones.     

Scanner influence on the mechanical response of QCT-based finite element analysis of long bones

Patient-specific QCT-based finite element (QCTFE) analyses enable highly accurate quantification of bone strength. We evaluated CT scanner influence on QCTFE models of long bones.A femur, humerus, and proximal femur without the head were scanned with K₂HPO₄ phantoms by seven CT scanners (four models) using typical clinical protocols. QCTFE models were constructed. The geometrical dimensions, as well as the QCT-values expressed in Hounsfield unit (HU) distribution was compared. Principal strains at representative regions of interest (ROIs), and maximum principal strains (associated with fracture risk) were compared. Intraclass correlation coefficients (ICCs) were calculated to evaluate strain prediction reliability for different scanners. Repeatability was examined by scanning the femur twice and comparing resulting QCTFE models.Maximum difference in geometry was 2.3%. HU histograms before phantom calibration showed wide variation between QCT scans; however, bone density histogram variability was reduced after calibration and algorithmic manipulation. Relative standard deviation (RSD) in principal strains at ROIs was <10.7%. ICC estimates between scanners were >0.9. Fracture-associated strain had 6.7%, 8.1%, and 13.3% maximum RSD for the femur, humerus, and proximal femur, respectively. The difference in maximum strain location was <2 mm. The average difference with repeat scans was 2.7%.Quantification of strain differences showed mean RSD bounded by ∼6% in ROIs. Fracture-associated strains in “regular” bones showed a mean RSD bounded by ∼8%. Strains were obtained within a ±10% difference relative to the mean; thus, in a longitudinal study only changes larger than 20% in the principal strains may be significant. ICCs indicated high reliability of QCTFE models derived from different scanners.     

运动与17β-雌二醇对去卵巢大鼠血清17β-雌二醇和后肢骨骨密度的影响

目的：观察跑台运动和17β-雌二醇（E2）对去卵巢大鼠血清E2水平和后肢骨骨密度（BMD）的影响。方法：按体重将120只成年雌性SD大鼠随机分为假手术、假手术运动、去卵巢、去卵巢运动、雌激素、雌激素加运动6个组。假手术运动、去卵巢运动和雌激素加运动组每周进行5次60min、18m／min的平坡跑台运动训练,雌激素和雌激素加运动组每周按体重颈部皮下注射3次E2,每次50μg／kg体重。分别在运动和给药正式处理7和14周时,用放射免疫法检测血清E2水平;用双能X线骨密度仪检测右侧胫骨和股骨BMD的变化。结果：运动和给药正式处理7周时,去卵巢组胫骨近端、股骨近端和远端BMD以及血清E2水平均显著低于假手术组;去卵巢运动加E2组股骨近端和远端BMD显著高于去卵巢组,E2组和去卵巢运动加E2组大鼠血清E2水平显著高于去卵巢组。运动和给药正式处理14周时,去卵巢组大鼠胫骨近端、股骨近端和远端BMD以及血清E2水平均显著低于假手术组,假手术运动组股骨近端BMD显著高于假手术组;去卵巢运动组、E2组和去卵巢运动加E2组大鼠血清E2水平显著高于去卵巢组,去卵巢运动组股骨远端BMD显著高于去卵巢组,E2组和去卵巢运动加E2组胫骨近端、股骨近端和远端BMD均显著高于去卵巢组;去卵巢运动组大鼠胫骨近端BMD和血清E2水平显著低于去卵巢运动加E2组,E2组胫骨近端BMD显著高于去卵巢运动加E2组。结论：E2和较高中等强度跑台运动对去卵巢大鼠股骨和胫骨松质骨骨量减缓的效应是独立的。     

Automated image registration for assessing three-dimensional alignment of entire lower extremity and implant position using bi-plane radiography

An automated image-matching technique is presented to assess alignment of the entire lower extremity for normal and implanted knees and the positioning of implants with respect to bone. Sawbone femur and tibia and femoral and tibial components of a total knee arthroplasty system were used. Three spherical markers were attached to each sawbone and each component to define the local coordinate system. Outlines of the three-dimensional (3D) bone models and component computer-aided design (CAD) models were projected onto extracted contours of the femur, tibia, and implants in frontal and oblique X-ray images. Three-dimensional position of each model was recovered by minimizing the difference between the projected outline and the contour. Median values of the absolute error in estimating relative positions were within 0.5 mm and 0.6° for the femur with respect to the tibia, 0.5 mm and 0.5° for the femoral component with respect to the tibial component, 0.6 mm and 0.6° for the femoral component with respect to the femur, and 0.5 mm and 0.4° for the tibial component with respect to the tibia, indicating significant improvements when compared to manually obtained results.     

中等强度跑台运动预防去卵巢大鼠骨质疏松

目的观察中等强度跑台运动对去卵巢大鼠骨质疏松的预防作用。方法将30只3月龄未经产雌性SD大鼠随机分为假手术、去卵巢静止和去卵巢运动三个组。去卵巢运动组每周进行4次时间45min、速度18m/min、坡度5°的跑台训练。实验结束时,检测血清雌二醇（E2）、碱性磷酸酶（ALP）、抗酒石酸酸性磷酸酶（TRAP）和骨钙素（BGP）水平以及右侧游离股骨和胫骨的骨密度（BMD）和骨矿物含量（BMC）;同时观察左侧股骨远端和胫骨近端组织形态学变化。结果与假手术组比较,去卵巢静止组大鼠血清ALP活性和BGP含量显著升高,血清TRAP活性和E2含量显著下降,股骨近段和远端以及胫骨近端BMD和BMC显著下降,股骨远端和胫骨近端骨小梁断裂增加、数目减少;与去卵巢静止组比较,去卵巢运动组大鼠血清E2和BGP含量显著上升,股骨三个部位以及胫骨近端BMD和BMC显著增加,股骨远端和胫骨近端骨小梁断裂减少、数目增加。结论中等强度跑台运动能增加去卵巢大鼠血清E2和BGP含量,改善去卵巢大鼠骨组织学结构。     

Three-dimensional reconstruction of fracture callus morphogenesis.

Rat and mouse femur and tibia fracture calluses were collected over various time increments of healing. Serial sections were produced at spatial segments across the fracture callus. Standard histological methods and in situ hybridization to col1a1 and col2a1 mRNAs were used to define areas of cartilage and bone formation as well as tissue areas undergoing remodeling. Computer-assisted reconstructions of histological sections were used to generate three-dimensional images of the spatial morphogenesis of the fracture calluses. Endochondral bone formation occurred in an asymmetrical manner in both the femur and tibia, with cartilage tissues seen primarily proximal or distal to the fractures in the respective calluses of these bones. Remodeling of the calcified cartilage proceeded from the edges of the callus inward toward the fracture producing an inner-supporting trabecular structure over which a thin outer cortical shell forms. These data suggest that the specific developmental mechanisms that control the asymmetrical pattern of endochondral bone formation in fracture healing recapitulated the original asymmetry of development of a given bone because femur and tibia grow predominantly from their respective distal and proximal physis. These data further show that remodeling of the calcified cartilage produces a trabecular bone structure unique to fracture healing that provides the rapid regain in weight-bearing capacity to the injured bone.     

Accuracy of femur reconstruction from sparse geometric data using a statistical shape model

Sparse geometric information from limited field-of-view medical images is often used to reconstruct the femur in biomechanical models of the hip and knee. However, the full femur geometry is needed to establish boundary conditions such as muscle attachment sites and joint axes which define the orientation of joint loads. Statistical shape models have been used to estimate the geometry of the full femur from varying amounts of sparse geometric information. However, the effect that different amounts of sparse data have on reconstruction accuracy has not been systematically assessed. In this study, we compared shape model and linear scaling reconstruction of the full femur surface from varying proportions of proximal and distal partial femur geometry in combination with morphometric and landmark data. We quantified reconstruction error in terms of surface-to-surface error as well as deviations in the reconstructed femur’s anatomical coordinate system which is important for biomechanical models. Using a partial proximal femur surface, mean shape model-based reconstruction surface error was 1.8 mm with 0.15° or less anatomic axis error, compared to 19.1 mm and 2.7–5.6° for linear scaling. Similar results were found when using a partial distal surface. However, varying amounts of proximal or distal partial surface data had a negligible effect on reconstruction accuracy. Our results show that given an appropriate set of sparse geometric data, a shape model can reconstruct full femur geometry with far greater accuracy than simple scaling.     

Mapping ligament insertion sites onto bone surfaces in knee by co-registration of CT and digitization data

This paper describes a new methodology that enables mapping of the ligament insertion sites onto bone surfaces in the knee joint by co-registration of the data acquired using digitization and computed tomography (CT). Local coordinate systems on the distal femur and proximal tibia were established by three spherical fiducial markers rigidly affixed to each bone. The fiducial marker centroid locations were identified by a least-squares sphere-fitting algorithm. An optimization correction procedure was proposed to mitigate the effect of the target registration error (TRE) on the alignment of coordinate systems for co-registration. A test with four cadaveric specimens demonstrated successful mapping of insertion sites for five ligaments. Fiducial registration error (FRE) as measured by the differences in inter-marker distances between the two modalities was smaller than 2%. The optimization procedure corrected the insertion site invisibility or partial visibility problem and improved the overall mapping quality, as indicated by substantial reduction of the mean and dispersion of distances from digitized insertion site points to the bone surfaces.     

Protective actions of green tea polyphenols and alfacalcidol on bone microstructure in female rats with chronic inflammation

This study investigated the effects of green tea polyphenols (GTP) and alfacalcidol on bone microstructure and strength along with possible mechanisms in rats with chronic inflammation. A 12-week study using a 2 (no GTP vs. 0.5%, w/v GTP in drinking water)×2 (no alfacalcidol vs. 0.05 μg/kg alfacalcidol orally, 5×/week) factorial design was employed in lipopolysaccharide (LPS)-administered female rats. A group receiving placebo administration was used to compare with a group receiving LPS administration only to evaluate the effect of LPS. Changes in tibial and femoral microarchitecture and strength of femur were evaluated. Difference in expression of tumor necrosis factor-α (TNF-α) in proximal tibia using immunohistochemistry was examined. Compared to the placebo group, the LPS-administered-only group had significantly lower femoral mass, trabecular volume, thickness and number in proximal tibia and femur, and lower periosteal bone formation rate in tibial shafts but had significantly higher trabecular separation and osteoclast number in proximal tibia and eroded surface in endocortical tibial shafts. Both GTP and alfacalcidol reversed these LPS-induced detrimental changes in femur, proximal tibia and endocortical tibial shaft. Both GTP and alfacalcidol also significantly improved femoral strength, while significantly suppressed TNF-α expression in proximal tibia. There were significant interactions in femoral mass and strength, trabecular separation, osteoclast number and TNF-α expression in proximal tibia. A combination of both showed to sustain bone microarchitecture and strength. We conclude that a protective impact of GTP and alfacalcidol in bone microarchitecture during chronic inflammation may be due to a suppression of TNF-α.     

Distal skeletal tibia assessed by HR-pQCT is highly correlated with femoral and lumbar vertebra failure loads

Dual energy X-ray absorptiometry (DXA) is the standard for assessing fragility fracture risk using areal bone mineral density (aBMD), but only explains 60–70% of the variation in bone strength. High-resolution peripheral quantitative computed tomography (HR-pQCT) provides 3D-measures of bone microarchitecture and volumetric bone mineral density (vBMD), but only at the wrist and ankle. Finite element (FE) models can estimate bone strength with 86–95% precision. The purpose of this study is to determine how well vBMD and FE bone strength at the wrist and ankle relate to fracture strength at the hip and spine, and to compare these relationships with DXA measured directly at those axial sites. Cadaveric samples (radius, tibia, femur and L4 vertebra) were compared within the same body. The radius and tibia specimens were assessed using HR-pQCT to determine vBMD and FE failure load. aBMD from DXA was measured at the femur and L4 vertebra. The femur and L4 vertebra specimens were biomechanically tested to determine failure load. aBMD measures of the axial skeletal sites strongly correlated with the biomechanical strength for the L4 vertebra (r = 0.77) and proximal femur (r = 0.89). The radius correlated significantly with biomechanical strength of the L4 vertebra for vBMD (r = 0.85) and FE-derived strength (r = 0.72), but not with femur strength. vBMD at the tibia correlated significantly with femoral biomechanical strength (r = 0.74) and FE-estimated strength (r = 0.83), and vertebral biomechanical strength for vBMD (r = 0.97) and FE-estimated strength (r = 0.91). The higher correlations at the tibia compared to radius are likely due to the tibia’s weight-bearing function.     

A three-dimensional musculoskeletal model for gait analysis. Anatomical variability estimates

Three-dimensional coordinates defining the origin and insertion of 40 muscle units, and bony landmarks for osteometric scaling were identified on dry bone specimens. Interspecimen coordinate differences along the anterior-posterior axis of the pelvis and the long bone axes of the pelvis, femur and leg were reduced by scaling but landmark differences along the other axes were not. The coordinates were mapped to living subjects using close-range photogrammetry to locate superficial reference markers. The error of predicting the positions of internal coordinates was assessed by comparing joint centre locations calculated from local axes defining the orientation of segments superior and inferior to a joint. A difference was attributed to: anatomical variability not accounted for by scaling; errors in identifying and placing reference landmarks; the accuracy of locating markers using photogrammetry and error introduced by marker oscillation during movement. Anatomical differences between specimens are one source of error in defining a musculoskeletal model but larger errors are introduced when such models are mapped to living subjects.     

Relative variation in human proximal and distal limb segment lengths

The pattern of variation and covariation of proximal and distal limb segment lengths was examined within and between 20 geographically diverse skeletal samples of modern humans. Analyses of variance-covariance matrices (VCMs) of logarithmically transformed (ln) variates of humerus, radius, femur, and tibia length were performed to test the following hypotheses: first, within populations, the distal and proximal segments will have equal relative (i.e., size-independent) variability. However, between populations, the tibia is predicted to be more variable than the other segments. Tests of fit of computed VCMs to theoretical matrices by an iterative procedure (Anderson [1973] Ann. Stat. 1:135-141) reject the equal variance hypotheses, rather suggesting that the relative variances of the distal limb segments are greater than are those of the proximal. Males and females differ somewhat in that within females, the distal segments of both limbs have equal variance, while within males, the tibia has greater relative variance than the radius. The second hypothesis, regarding between-group variability, is somewhat supported in that between human populations, one cannot reject that the tibia has greater relative variance than the other limb segments. However, neither can one reject an alternative hypothesis that both distal limb segments (tibia and radius) are more variable than the proximal segments. Differential growth allometry is explored, and likely plays a major role in differences seen both within and between human populations.     

Resveratrol Treatment Delays Growth Plate Fusion and Improves Bone Growth in Female Rabbits

Trans-resveratrol (RES), naturally produced by many plants, has a structure similar to synthetic estrogen diethylstilbestrol, but any effect on bone growth has not yet been clarified. Pre-pubertal ovary-intact New Zealand white rabbits received daily oral administration of either vehicle (control) or RES (200 mg/kg) until growth plate fusion occurred. Bone growth and growth plate size were longitudinally monitored by X-ray imaging, while at the endpoint, bone length was assessed by a digital caliper. In addition, pubertal ovariectomized (OVX) rabbits were treated with vehicle, RES or estradiol cypionate (positive control) for 7 or 10 weeks and fetal rat metatarsal bones were cultured in vitro with RES (0.03 µM–50 µM) and followed for up to 19 days. In ovary-intact rabbits, sixteen-week treatment with RES increased tibiae and vertebrae bone growth and subsequently improved final length. In OVX rabbits, RES delayed fusion of the distal tibia, distal femur and proximal tibia epiphyses and femur length and vertebral bone growth increased when compared with controls. Histomorphometrical analysis showed that RES-treated OVX rabbits had a wider distal femur growth plate, enlarged resting zone, increased number/size of hypertrophic chondrocytes, increased height of the hypertrophic zone, and suppressed chondrocyte expression of VEGF and laminin. In cultured fetal rat metatarsal bones, RES stimulated growth at 0.3 µM while at higher concentrations (10 μM and 50 μM) growth was inhibited. We conclude that RES has the potential to improve longitudinal bone growth. The effect was associated with a delay of growth plate fusion resulting in increased final length. These effects were accompanied by a profound suppression of VEGF and laminin expression suggesting that impairment of growth plate vascularization might be an underlying mechanism.     

Evaluation of a morphing based method to estimate muscle attachment sites of the lower extremity

To generate subject-specific musculoskeletal models for clinical use, the location of muscle attachment sites needs to be estimated with accurate, fast and preferably automated tools. For this purpose, an automatic method was used to estimate the muscle attachment sites of the lower extremity, based on the assumption of a relation between the bone geometry and the location of muscle attachment sites. The aim of this study was to evaluate the accuracy of this morphing based method. Two cadaver dissections were performed to measure the contours of 72 muscle attachment sites on the pelvis, femur, tibia and calcaneus. The geometry of the bones including the muscle attachment sites was morphed from one cadaver to the other and vice versa. For 69% of the muscle attachment sites, the mean distance between the measured and morphed muscle attachment sites was smaller than 15 mm. Furthermore, the muscle attachment sites that had relatively large distances had shown low sensitivity to these deviations. Therefore, this morphing based method is a promising tool for estimating subject-specific muscle attachment sites in the lower extremity in a fast and automated manner.     

Elasticity–density and viscoelasticity–density relationships at the tibia mid-diaphysis assessed from resonant ultrasound spectroscopy measurements

Cortical bone tissue is an anisotropic material characterized by typically five independent elastic coefficients (for transverse isotropy) governing shear and longitudinal deformations in the different anatomical directions. It is well established that the Young’s modulus in the direction of the bone axis of long bones has a strong relationship with mass density. It is not clear, however, whether relationships of similar strength exist for the other elastic coefficients, for they have seldom been investigated, and the results available in the literature are contradictory. The objectives of the present work were to document the anisotropic elastic properties of cortical bone at the tibia mid-diaphysis and to elucidate their relationships with mass density. Resonant ultrasound spectroscopy (RUS) was used to measure the transverse isotropic stiffness tensor of 55 specimens from 19 donors. Except for Poisson’s ratios and the non-diagonal stiffness coefficient, strong linear correlations between the different elastic coefficients \((0.7 < {r^{2}} < 0.99)\) and between these coefficients and density \((0.79 < {r^{2}} < 0.89)\) were found. Comparison with previously published data from femur specimens suggested that the strong correlations evidenced in this study may not only be valid for the mid-tibia. RUS also measures the viscous part of the stiffness tensor. An anisotropy ratio close to two was found for damping coefficients. Damping increased as the mass density decreased. The data suggest that a relatively accurate estimation of all the mid-tibia elastic coefficients can be derived from mass density. This is of particular interest (1) to design organ-scale bone models in which elastic coefficients are mapped according to Hounsfield values from computed tomography scans as a surrogate for mass density and (2) to model ultrasound propagation at the mid-tibia, which is an important site for the in vivo assessment of bone status with axial transmission techniques.