Linear measurements of cortical bone and dental enamel by computed tomography: Applications and problems

This paper explores the potential of high-resolution computed tomography (CT) as a morphometric tool in paleoanthropology. The accuracy of linear measurements of enamel thickness and cortical bone thickness taken from CT scans is evaluated by making comparison with measurements taken directly from physical sections. The measurements of cortical bone are taken on extant and fossil specimens with and without attached matrix, and the dental specimens studied include a sample of 12 extant human molars. Local CT numbers (representing X-ray, attenuation) are used to determine the exact position of the boundaries of a structure. Using this technique most studied dimensions, including four of human molar enamel thickness, could be obtained from CT scans with a maximum error range of ±0.1 mm. The limitations of the method are discussed with special reference to problems associated with highly mineralized fossils. © 1993 Wiley-Liss, Inc.     

An accurate estimation of bone density improves the accuracy of subject-specific finite element models

An experimental-numerical study was performed to investigate the relationships between computed tomography (CT)-density and ash density, and between ash density and apparent density for bone tissue, to evaluate their influence on the accuracy of subject-specific FE models of human bones. Sixty cylindrical bone specimens were examined. CT-densities were computed from CT images while apparent and ash densities were measured experimentally. The CT/ash-density and ash/apparent-density relationships were calculated. Finite element models of eight human femurs were generated considering these relationships to assess their effect on strain prediction accuracy. CT and ash density were linearly correlated (R(2)=0.997) over the whole density range but not equivalent (intercep t <0, slope >1). A constant ash/apparent-density ratio (0.598+/-0.004) was found for cortical bone. A lower ratio, with a larger dispersion, was found for trabecular bone (0.459+/-0.100), but it became less dispersed, and equal to that of cortical tissue, when testing smaller trabecular specimens (0.598+/-0.036). This suggests that an experimental error occurred in apparent-density measurements for large trabecular specimens and a constant ratio can be assumed valid for the whole density range. Introducing the obtained relationships in the FE modelling procedure improved strain prediction accuracy (R(2)=0.95, RMSE=7%). The results suggest that: (i) a correction of the densitometric calibration should be used when evaluating bone ash-density from clinical CT scans, to avoid ash-density underestimation and overestimation for low- and high-density bone tissue, respectively; (ii) the ash/apparent-density ratio can be assumed constant in human femurs and (iii) the correction improves significantly the model accuracy and should be considered in subject-specific bone modelling.     

A semi-automated method for hexahedral mesh construction of human vertebrae from CT scans

Generation of finite element (FE) meshes of vertebrae from computed tomography (CT) scans is labour intensive due to their geometric complexity. As such, techniques that simplify creation of meshes of vertebrae are needed to make FE analysis feasible for large studies and clinical applications. Techniques to obtain a geometric representation of bone contours from CT scans of vertebrae and construct a hexahedral mesh from the contours were developed. An automated edge detection technique was developed to identify surface contours of the vertebrae, followed by atlas based B-spline curve fitting to construct curves from the edge points. The method was automatic and robust to missing data, with a controllable degree of smoothing and interpolation. Parametric mapping was then used to generate nodes for each CT slice, which were connected between slices to obtain a hexahedral mesh. This method could be adapted for modelling a variety of orthopaedic structures.     

Comparative analysis of bone remodelling models with respect to computerised tomography-based finite element models of bone

Subject-specific finite element models are an extensively used tool for the numerical analysis of the biomechanical behaviour of human bones. However, bone modelling is not an easy task due to the complex behaviour of bone tissue, involving non-homogeneous and anisotropic mechanical properties. Moreover, bone is a living tissue and therefore its microstructure and mechanical properties evolve with time in a known process called bone remodelling. This phenomenon has been widely studied, many being the numerical models that have been formulated to predict density distribution and its evolution in several bones. The aim of the present study is to assess the capability of a bone remodelling model to predict the bone density distribution of different types of human bone (femur, tibia and mandible) comparing the obtained results with the bone density estimated by means of computerised tomography. Good accuracy was observed for the bone remodelling predictions including the thickness of the cortical layer.     

About the inevitable compromise between spatial resolution and accuracy of strain measurement for bone tissue: A 3D zero-strain study

The accurate measurement of local strain is necessary to study bone mechanics and to validate micro computed tomography (µCT) based finite element (FE) models at the tissue scale. Digital volume correlation (DVC) has been used to provide a volumetric estimation of local strain in trabecular bone sample with a reasonable accuracy. However, nothing has been reported so far for µCT based analysis of cortical bone. The goal of this study was to evaluate accuracy and precision of a deformable registration method for prediction of local zero-strains in bovine cortical and trabecular bone samples. The accuracy and precision were analyzed by comparing scans virtually displaced, repeated scans without any repositioning of the sample in the scanner and repeated scans with repositioning of the samples.     

Cortical thickness measurements and photon absorptiometry for determination of bone quantity.

Reproducibility of estimates of cortical bone quantity in the appendicular skeleton by two methods was studied in healthy individuals and patients undergoing long-term hemodialysis. Repeated measurements of cortical thickness (CT) at the midpoint of the second metacarpal were taken from single radiographs of both hands by two independent observers. Repeated measurements by the same observer were more reproducible and the degree of reproducibility was far greater in healthy subjects than in dialysis patients. Repeated measurements were made of bone mineral content (BMC) and bone width (W) of the distal radius by photon absorptiometry. Repeated BMC/W determinations were highly reproducible in both healthy subjects and dialysis patients. High correlation was found between BMC and cross-sectional cortical area and between both simple cortical thickness and cortical area/width. Thus the photon absorptiometric technique is superior for the serial monitoring of bone quantity, particularly in patients with uremic osteodystrophy, but results obtained by the two methods in group studies should be comparable.     

Effects of the Neurogranin Variant rs12807809 on Thalamocortical Morphology in Schizophrenia

Although the genome wide supported psychosis susceptibility neurogranin (NRGN) gene is expressed in human brains, it is unclear how it impacts brain morphology in schizophrenia. We investigated the influence of NRGN rs12807809 on cortical thickness, subcortical volumes and shapes in patients with schizophrenia. One hundred and fifty six subjects (91 patients with schizophrenia and 65 healthy controls) underwent structural MRI scans and their blood samples were genotyped. A brain mapping algorithm, large deformation diffeomorphic metric mapping, was used to perform group analysis of subcortical shapes and cortical thickness. Patients with risk TT genotype were associated with widespread cortical thinning involving frontal, parietal and temporal cortices compared with controls with TT genotype. No volumetric difference in subcortical structures (hippocampus, thalamus, amygdala, basal ganglia) was observed between risk TT genotype in patients and controls. However, patients with risk TT genotype were associated with thalamic shape abnormalities involving regions related to pulvinar and medial dorsal nuclei. Our results revealed the influence of the NRGN gene on thalamocortical morphology in schizophrenia involving widespread cortical thinning and thalamic shape abnormalities. These findings help to clarify underlying NRGN mediated pathophysiological mechanisms involving cortical-subcortical brain networks in schizophrenia.     

Relationships between bone mass and micro-architecture at the mandible and iliac bone in edentulous subjects: a dual X-ray absorptiometry, computerised tomography and microcomputed tomography study

doi: 10.1111/j.1741‐2358.2011.00527.x Relationships between bone mass and micro‐architecture at the mandible and iliac bone in edentulous subjects: a dual X‐ray absorptiometry, computerised tomography and microcomputed tomography study Objectives: To compare bone volume, bone mineral density, cortical thickness and bone micro‐architecture in a series of paired mandibular and iliac bone samples analysed by various imagery techniques to see whether relationships exist between the various techniques and between mandibular and iliac bone. Materials and methods: Bone samples from the mandible and ilium were harvested in 20 cadavers and analysed by dual energy X‐ray absorptiometry (DXA), computerised tomography (CT) on a conventional hospital machine and microCT. Results: Significant correlations were found between Hounsfield density obtained by CT, and bone mass determined by microCT but not with DXA values. Cortical thickness measurements were well correlated between CT and microCT. No relationships were found between mandibular and iliac bone, when considering mineral density, cortical thickness, bone volume or micro‐architecture. Conclusion: In clinical practice, CT remains the most appropriate routine means for bone qualitative and quantitative evaluation at the mandible. In this ex vivo study, these results confirm that mandibular bone status does not reflect the axial skeletal one and assist in the placement of implants with dental prostheses in old or osteoporotic patients.     

The stiffness of different sections of the cartilaginous surface of the femur head

As a result of investigations on microstiffness of the cartilagenous covering of the human femoral bone head, normal (20), different values of microstiffness have been defined in various zones, depending on their functional load and thickness of the cartilage. The data obtained can serve for construction of subtotal endoprostheses of the coxal joint.     

A Comparison of Micro-CT and Dental CT in Assessing Cortical Bone Morphology and Trabecular Bone Microarchitecture

Objective

The objective of this study was to evaluate the relationship between the trabecular bone microarchitecture and cortical bone morphology by using micro-computed tomography (micro-CT) and dental cone-beam computed tomography (dental CT).

Materials and Methods

Sixteen femurs and eight fifth lumbar vertebrae were collected from eight male Sprague Dawley rats. Four trabecular bone microarchitecture parameters related to the fifth lumbar vertebral body (percent bone volume [BV/TV], trabecular thickness [TbTh], trabecular separation [TbSp], and trabecular number [TbN]) were calculated using micro-CT. In addition, the volumetric cancellous bone grayscale value (vCanGrayscale) of the fifth lumbar vertebral body was measured using dental CT. Furthermore, four cortical bone morphology parameters of the femoral diaphysis (total cross-sectional area [TtAr], cortical area [CtAr], cortical bone area fraction [CtAr/TtAr], and cortical thickness [CtTh]) were calculated using both micro-CT and dental CT. Pearson analysis was conducted to calculate the correlation coefficients (r) of the micro-CT and dental CT measurements. Paired-sample t tests were used to compare the differences between the measurements of the four cortical bone morphology parameters obtained using micro-CT and dental CT.

Results

High correlations between the vCanGrayscale measured using dental CT and the trabecular bone microarchitecture parameters (BV/TV [r = 0.84] and TbTh [r = 0.84]) measured using micro-CT were observed. The absolute value of the four cortical bone morphology parameters may be different between the dental CT and micro-CT approaches. However, high correlations (r ranged from 0.71 to 0.90) among these four cortical bone morphology parameters measured using the two approaches were obtained.

Conclusion

We observed high correlations between the vCanGrayscale measured using dental CT and the trabecular bone microarchitecture parameters (BV/TV and TbTh) measured using micro-CT, in addition to high correlations between the cortical bone morphology measured using micro-CT and dental CT. Further experiments are necessary to validate the use of dental CT on human bone.     

Semiautomatic extraction of cortical thickness and diaphyseal curvature from CT scans

The understanding of locomotor patterns, activity schemes, and biological variations has been enhanced by the study of the geometrical properties and cortical bone thickness of the long bones measured using CT scan cross‐sections. With the development of scanning procedures, the internal architecture of the long bones can be explored along the entire diaphysis. Recently, several methods that map cortical thickness along the whole femoral diaphysis have been developed. Precise homology is vital for statistical examination of the data; however, the repeatability of these methods is unknown and some do not account for the curvature of the bones. We have designed a semiautomatic workflow that improves the morphometric analysis of cortical thickness, including robust data acquisition with minimal user interaction and considering the bone curvature. The proposed algorithm also performs automatic landmark refinement and rigid registration on the extracted morphometric maps of the cortical thickness. Because our algorithm automatically reslices the diaphysis into 100 cross‐sections along the medial axis and uses an adaptive thresholding method, it is usable on CT scans that contain soft tissues as well as on bones that have not been oriented specifically prior to scanning. Our approach exhibits considerable robustness to error in user‐supplied landmarks, suppresses distortion caused by the curvature of the bones, and calculates the curvature of the medial axis.     

A novel approach to estimate trabecular bone anisotropy using a database approach

Continuum finite element (FE) models of bones have become a standard pre-clinical tool to estimate bone strength. These models are usually based on clinical CT scans and material properties assigned are chosen as isotropic based only on the density distribution. It has been shown, however, that trabecular bone elastic behavior is best described as orthotropic. Unfortunately, the use of orthotropic models in FE analysis derived from CT scans is hampered by the fact that the measurement of a trabecular orientation (fabric) is not possible from clinical CT images due to the low resolution of such images. In this study, we explore the concept of using a database (DB) of high-resolution bone models to derive the fabric information that is missing in clinical images. The goal of this study was to investigate if models with fabric derived from a relatively small database can already produce more accurate results than isotropic models.     

Magnetic resonance imaging for in vivo assessment of three-dimensional patellar tracking

We have developed a non-invasive measurement technique which can ultimately be used to quantify three-dimensional patellar kinematics of human subjects for a range of static positions of loaded flexion and assessed its accuracy. Knee models obtained by segmenting and reconstructing one high-resolution scan of the knee were registered to bone outlines obtained by segmenting fast, low-resolution scans of the knee in static loaded flexion. We compared patellar tracking measurements made using the new method to measurements made using Roentgen stereophotogrammetric analysis in three cadaver knee specimens loaded through a range of flexion in a test rig. The error in patellar spin and tilt measurements was less than 1.02 degrees and the error in lateral patellar shift was 0.88 mm. Sagittal plane scans provided more accurate final measurements of patellar spin and tilt, whereas axial plane scans provided more accurate measurements of lateral translation and patellar flexion. Halving the number of slices did not increase measurement error significantly, which suggests that scan times can be reduced without reducing accuracy significantly. The method is particularly useful for multiple measurements on the same subject because the high-resolution bone-models need only be created once; thus, the potential variability in coordinate axes assignment and model segmentation during subsequent measurements is removed.     

Development and validation of a subject-specific finite element model of a human clavicle

This study aimed to develop and validate a finite element (FE) model of a human clavicle which can predict the structural response and bone fractures under both axial compression and anterior–posterior three-point bending loads. Quasi-static non-injurious axial compression and three-point bending tests were first conducted on a male clavicle followed by a dynamic three-point bending test to fracture. Then, two types of FE models of the clavicle were developed using bone material properties which were set to vary with the computed tomography image density of the bone. A volumetric solid FE model comprised solely of hexahedral elements was first developed. A solid-shell FE model was then created which modelled the trabecular bone as hexahedral elements and the cortical bone as quadrilateral shell elements. Finally, simulations were carried out using these models to evaluate the influence of variations in cortical thickness, mesh density, bone material properties and modelling approach on the biomechanical responses of the clavicle, compared with experimental data. The FE results indicate that the inclusion of density-based bone material properties can provide a more accurate reproduction of the force–displacement response and bone fracture timing than a model with uniform bone material properties. Inclusion of a variable cortical thickness distribution also slightly improves the ability of the model to predict the experimental response. The methods developed in this study will be useful for creating subject-specific FE models to better understand the biomechanics and injury mechanism of the clavicle.     

A technique for the quantification of the 3D connectivity of thin articulations in Bony sutures

The anatomy and development of cranial and facial sutures have been studied in detail using histological sections, 2D radiographs and more recently CT imaging. However, little attention has been paid to evaluating and quantifying the connectivity of these thin cortical bone articulations. More recent technological advances such as micro-CT imaging has the potential to be used to provide quantitative measurements of 3D connectivity in bony articulations. This study presents a new technique for quantifying the connectivity of bony projections inside cranial and facial sutures using a combination of skeletonization, thinning algorithms and 3D intensity mapping. The technique is demonstrated in five sutures through semi-automated analysis and image processing of μCT scans. In the sagittal, coronal and frontozygomatic sutures an average bone connectivity of 6.6–11.6% was found with multiple bony projections providing an interlocking structure between adjacent bones. Much higher bone connectivity was present in the zygomaticotemporal and zygomaticomaxillary sutures (22.7–37.4%) with few bony projections. This method combining μCT scanning and image processing techniques was successfully used to quantify the connectivity of thin bone articulations and allowed detailed assessment of sutural fusion in 3D. The wider application of this technique may allow quantification of connectivity in other structures, in particular fracture healing of long bones.     

Incidence and mechanical significance of pneumatization in the long bones of birds

The incidence of pneumatization in avian long bones was studied, by direct observation, in a large sample of species. Only proximal bones (humerus and femur) presented pneumatization in the sample studied. The incidence obtained was related to the variation of the maximum cortical thickness and mechanical properties, such as bending strength and flexural Young's modulus. Cortical thickness, bending strength and flexural Young's modulus were significantly lower in pneumatized bones than in marrow-filled bones. Furthermore, some congruence was found between pneumatization and systematic groups when compared. In this sense, Charadriformes was the only order studied with total absence of long bone pneumatization. Results on cortical thickness appear to be in agreement with modelling predictions previously made and with results obtained on other groups of flying vertebrates. The possible selective advantage of reduction in cortical thickness in relation to flying is suggested.     

Cortical bone mapping: An application to hand and foot bones in hominoids

Bone form reflects both the genetic profile and behavioural history of an individual. As cortical bone is able to remodel in response to mechanical stimuli, interspecific differences in cortical bone thickness may relate to loading during locomotion or manual behaviours during object manipulation. Here, we test the application of a novel method of cortical bone mapping to the third metacarpal (Mc3) and talus of Pan, Pongo, and Homo. This method of analysis allows measurement of cortical thickness throughout the bone, and as such is applicable to elements with complex morphology. In addition, it allows for registration of each specimen to a canonical surface, and identifies regions where cortical thickness differs significantly between groups. Cortical bone mapping has potential for application to palaeoanthropological studies; however, due to the complexity of correctly registering homologous regions across varied morphology, further methodological development would be advantageous.     

Subject-specific finite element models of long bones: An in vitro evaluation of the overall accuracy

The determination of the mechanical stresses induced in human bones is of great importance in both research and clinical practice. Since the stresses in bones cannot be measured non-invasively in vivo, the only way to estimate them is through subject-specific finite element modelling. Several methods exist for the automatic generation of these models from CT data, but before bringing them in the clinical practice it is necessary to assess their accuracy in the predictions of the bone stresses. Particular attention should be paid to those regions, like the epiphyseal and metaphyseal parts of long bones, where the automatic methods are typically less accurate. Aim of the present study was to implement a general procedure to automatically generate subject-specific finite element models of bones from CT data and estimate the accuracy of this general procedure by applying it to one real femur. This femur was tested in vitro under five different loading scenarios and the results of these tests were used to verify how the adoption of a simplified two-material homogeneous model would change the accuracy with respect to the density-based inhomogeneous one, with special attention paid to the epiphyseal and metaphyseal proximal regions of the bone. The results showed that the density-based inhomogeneous model predicts with a very good accuracy the measured stresses (R(2)=0.91, RMSE=8.6%, peak error=27%), while the two-material model was less accurate (R(2)=0.89, RMSE=9.6%, peak error=35%). The results showed that it is possible to automatically generate accurate finite element models of bones from CT data and that the strategy of material properties mapping has a significant influence on its accuracy.     

Intraoperative computer tomography control within the scope of maxillofacial traumatology using a mobile scanner]

Advances in intra-operative imaging and the development of new minimally invasive techniques are having an ever greater impact on modern surgery. Mobile CT scanners in the operating room is a new technique that permits image-guided surgery, and helps minimize postoperative complications. We report on our initial experience with intraoperative CT scanning during surgery on patients suffering lateral midface trauma. A mobile CT unit, the Tomoscan M (Philips, Utrecht, Netherlands) set up in the operating room, was evaluated in 6 patients with zygomatic bone fractures. The patients were placed on the CT scanner table, which is detachable from the gantry. The unit is powered by batteries charged from an ordinary ring mains supply via a conventional plug. The CT images obtained were of good quality in all cases. No technical problems were observed during surgery. Using repeat CT scans, the procedure also permits accurate intraoperative monitoring of the anatomical repositioning of the bone fragments, and accurate implantation. No intraoperative or early postoperative complications were observed. This new technical aid ensures highly accurate reduction of the bone fragments, and minimizes the need for reoperation. High-quality intraoperative imaging with surgical navigation increase surgical outcome, and expand the spectrum of minimally invasive surgery.     

Stress distribution and displacement analysis during an intermaxillary disjunction--a three-dimensional FEM study of a human skull

The goal of this study was to contribute to an understanding of how much expansion force is needed during a maxillary expansion (ME) and where bony reaction takes place. A finite element (FE) model of a dry human male skull was generated from CT scans. The FE model, which consists of cortical and cancellous bone and teeth, was loaded with the same force magnitudes, directions and working points as in rapid maxillary expansion (RME). A three-dimensional finite element stress analysis (FESA) of the forces and displacement was performed. The highest stress was observed in the maxilla in the region where the forces were applied, and spreads more or less throughout almost the whole frontal skull structures. The displacement distribution which causes stress in the skull is highly dependant on the thickness of the bone and its structure. All areas with high compressive and tensile stress are exactly the regions which determine the maximal amount of force to be used during the maxillary expansion and should be examined in case of any complication during a patient's treatment. Regions with significant compressive and tensile stress are the regions observed to have an increase in cellular activity. Further simulations with a given displacement (0.5mm) showed that displacement simulations need extra caution otherwise they will lead to very high forces which are not realistic in an orthodontic treatment.