Three-dimensional image registration of MR proximal femur images for the analysis of trabecular bone parameters

Use of image registration and fusion algorithms and techniques in radiotherapy: Report of the AAPM Radiation Therapy Committee Task Group No. 132

Image registration and fusion algorithms exist in almost every software system that creates or uses images in radiotherapy. Most treatment planning systems support some form of image registration and fusion to allow the use of multimodality and time-series image data and even anatomical atlases to assist in target volume and normal tissue delineation. Treatment delivery systems perform registration and fusion between the planning images and the in-room images acquired during the treatment to assist patient positioning. Advanced applications are beginning to support daily dose assessment and enable adaptive radiotherapy using image registration and fusion to propagate contours and accumulate dose between image data taken over the course of therapy to provide up-to-date estimates of anatomical changes and delivered dose. This information aids in the detection of anatomical and functional changes that might elicit changes in the treatment plan or prescription. As the output of the image registration process is always used as the input of another process for planning or delivery, it is important to understand and communicate the uncertainty associated with the software in general and the result of a specific registration. Unfortunately, there is no standard mathematical formalism to perform this for real-world situations where noise, distortion, and complex anatomical variations can occur. Validation of the software systems performance is also complicated by the lack of documentation available from commercial systems leading to use of these systems in undesirable 'black-box' fashion. In view of this situation and the central role that image registration and fusion play in treatment planning and delivery, the Therapy Physics Committee of the American Association of Physicists in Medicine commissioned Task Group 132 to review current approaches and solutions for image registration (both rigid and deformable) in radiotherapy and to provide recommendations for quality assurance and quality control of these clinical processes.

Bone marrow lesions are associated with altered trabecular morphometry

OBJECTIVE

Bone marrow lesions (BMLs) are a common magnetic resonance (MR) feature in patients with osteoarthritis, however their pathological basis remains poorly understood and has not been evaluated in vivo. Our aim was to evaluate the trabecular structure associated with the presence and size of BMLs present in the same regions of interest (ROI) using quantitative MR-based trabecular morphometry.

DESIGN

158 participants in the Osteoarthritis Initiative (OAI) were imaged with a coronal 3D fast imaging with steady state precession (FISP) sequence for trabecular morphometry in the same session as the OAI 3 T MR knee evaluation. The proximal medial tibial subchondral bone in the central weight-bearing ROI on these knee 3D FISP images were quantitatively evaluated for apparent bone volume fraction, trabecular number, spacing, and thickness. BMLs were also evaluated in the subchondral bone immediately adjacent to the articular cartilage. BML volume was also evaluated within the same trabecular morphometry ROI and semi-quantitatively classified as none, small, or large. Kruskal-Wallis test was used to determine if mean apparent bone volume fraction, trabecular number, spacing, or thickness differed by BML score.

RESULTS

Compared to knees with ROIs containing no BMLs, knees with small or large BMLs had statistically higher apparent bone volume fraction (P < 0.01), trabecular number (P < 0.01), and thickness (P = 0.02), and lower trabecular spacing (P < 0.01).

CONCLUSIONS

Compared to knees with ROIs containing no BMLs, knees with ROIs containing small or large BMLs had higher apparent bone volume fraction, trabecular number and thickness, but lower trabecular spacing. These findings may represent areas of locally increased bone remodeling or compression.

Noncontact ultrasound imaging applied to cortical bone phantoms

PURPOSE

The purpose of this paper was to take the first steps toward applying noncontact ultrasound (NCU) to the tasks of monitoring osteoporosis and quantitative ultrasound imaging (QUS) of cortical bone. The authors also focused on the advantages of NCU, such as its lack of reliance on a technologist to apply transducers and a layer of acoustical coupling gel, the ability of the transducers to operate autonomously as specified by preprogrammed software, and the likely reduction in statistical and systematic errors associated with the variability in the pressure applied by the clinician to the transmitting transducer that NCU might provide. The authors also undertook this study in order to find additional applications of NCU beyond its past limited usage in assessing the severity of third degree burns.

METHODS

A noncontact ultrasound imaging system using a pair of specially designed broadband, 1.5 MHz noncontact piezoelectric transducers and cortical bone phantoms, were used to determine bone mineral density (BMD), speed of sound (SOS), integrated response (IR), and ultrasonic transmittance. Air gaps of greater than 3 cm, two transmission and two reflection paths, and a digital signal processor were also used in the collection of data from phantoms of nominal mass densities that varied from 1.17 to 2.25 g/cm(3) and in bone mineral density from 0 to 1.7 g/cm(3).

RESULTS

Good correlations between known BMD and measured SOS, IR, and transmittance were obtained for all 17 phantoms, and methods for quantifying and minimizing sources of systematic errors were outlined. The BMD of the phantom sets extended through most of the in vivo range found in cortical bone. A total of 16-20 repeated measurements of the SOS, thickness, and IR for the phantom set that were conducted over a period of several months showed a small variation in the range of measurements of ±1%-2%. These NCU data were shown to be in agreement with similar results using contact ultrasound to be within 1%-2%. Transmittance images of cortical bone phantoms showed differences in the nominal overall BMD values of the phantoms that were large enough to be distinguished by a visual examination. A list of possible sources of errors in quantitative NCU was also included in this study.

CONCLUSIONS

The results of this paper suggest that NCU might find additional applications in medical imaging, beyond its original and only previous usage in assessing third degree burns. The fact that the authors' phantom measurements using conventional, gel coupled ultrasound are in agreement with those obtained with NCU demonstrates that in spite of large additional levels of attenuation of up to 150 dB and new error sources, NCU could have comparable levels of accuracy to those of conventional quantitative ultrasound, while providing the medical and patient comfort-related advantages of not involving direct contact.

Image interpolation allows accurate quantitative bone morphometry in registered micro-computed tomography scans

Time-lapsed in vivo micro-computed tomography is a powerful tool to analyse longitudinal changes in the bone micro-architecture. Registration can overcome problems associated with spatial misalignment between scans; however, it requires image interpolation which might affect the outcome of a subsequent bone morphometric analysis. The impact of the interpolation error itself, though, has not been quantified to date. Therefore, the purpose of this ex vivo study was to elaborate the effect of different interpolator schemes [nearest neighbour, tri-linear and B-spline (BSP)] on bone morphometric indices. None of the interpolator schemes led to significant differences between interpolated and non-interpolated images, with the lowest interpolation error found for BSPs (1.4%). Furthermore, depending on the interpolator, the processing order of registration, Gaussian filtration and binarisation played a role. Independent from the interpolator, the present findings suggest that the evaluation of bone morphometry should be done with images registered using greyscale information.

Generation of an atlas of the proximal femur and its application to trabecular bone analysis

Automatic placement of anatomically corresponding volumes of interest and comparison of parameters against a standard of reference are essential components in studies of trabecular bone. Only recently, in vivo MR images of the proximal femur, an important fracture site, could be acquired with high-spatial resolution. The purpose of this MRI trabecular bone study was two-fold: (1) to generate an atlas of the proximal femur to automatically place anatomically corresponding volumes of interest in a population study and (2) to demonstrate how mean models of geodesic topological analysis parameters can be generated to be used as potential standard of reference. Ten females were used to generate the atlas and geodesic topological analysis models, and 10 females were used to demonstrate the atlas-based trabecular bone analysis. All alignments were based on three-dimensional (3D) multiresolution affine transformations followed by 3D multiresolution free-form deformations. Mean distances less than 1 mm between aligned femora, and sharp edges in the atlas and in fused gray-level images of registered femora indicated that the anatomical variability was well accommodated and explained by the free-form deformations.

Bone densitometry using x-ray spectra

In contrast to the two distinct energy regions that are involved in dual-energy x-ray absorptiometry for bone densitometry, the complete spectrum of a beam transmitted through two layers of different materials is utilized in this study to calculate the areal density of each material. Test objects constructed from aluminum and Plexiglas were used to simulate cortical bone and soft tissue, respectively. Solid-state HPGe (high-purity germanium) detectors provided high-resolution x-ray spectra over an energy range of approximately 20-80 keV. Areal densities were obtained from spectra using two methods: a system of equations for two spectral regions and a nonlinear fit of the entire spectrum. Good agreement with the known areal densities of aluminum was obtained over a wide range of PMMA thicknesses. The spectral method presented here can be used to decrease beam hardening at a small number of bodily points selected for examination.

High-resolution imaging techniques for the assessment of osteoporosis

The importance of assessing the bone's microarchitectural make-up in addition to its mineral density in the context of osteoporosis has been emphasized in several publications. The high spatial resolution required to resolve the bone's microstructure in a clinically feasible scan time is challenging. At present, the best suited modalities meeting these requirements in vivo are high-resolution peripheral quantitative imaging (HR-pQCT) and magnetic resonance imaging (MRI). Whereas HR-pQCT is limited to peripheral skeleton regions like the wrist and ankle, MRI can also image other sites like the proximal femur but usually with lower spatial resolution. In addition, multidetector computed tomography has been used for high-resolution imaging of trabecular bone structure; however, the radiation dose is a limiting factor. This article provides an overview of the different modalities, technical requirements, and recent developments in this emerging field. Details regarding imaging protocols as well as image postprocessing methods for bone structure quantification are discussed.

Quantifying the material and structural determinants of bone strength

The ability of a bone to resist fracture depends on the amount of bone present, the spatial distribution of the bone mass as cortical and trabecular bone and the intrinsic properties of the bone material. Whereas low areal bone mineral density (aBMD) predicts fractures, its sensitivity and specificity is low, as over 50% of fractures occur in persons without osteoporosis by BMD testing and most women with osteoporosis do not sustain a fracture. New non-invasive imaging techniques, including three-dimensional (3D) assessments of bone density and geometry, microarchitecture and integrated measurements of bone strength such as finite element analysis (FEA), provide estimates of bone strength that can be used to increase the sensitivity and specificity of fracture risk assessment. Initial observations have shown that these techniques provide information that will improve our understanding of the pathophysiology of skeletal fragility and suggest that these techniques are likely to have a role in the clinical management of individuals at risk for fracture.

Evaluation of correction methods for coil-induced intensity inhomogeneities and their influence on trabecular bone structure parameters from MR images

Magnetic resonance (MR) imaging-based quantitative trabecular bone structure analysis has gained increasing interest in osteoporotic fracture risk assessment and treatment evaluation related to osteoporosis. In vivo MR images of anatomic regions such as the proximal femur and distal tibia are generally acquired with a surface coil in order to obtain sufficient sensitivity and resolution for quantification of the trabeculae. However, these coils introduce intensity inhomogeneities which affect the trabecular bone structure analysis. This work evaluates the applicability of a fully automatic coil correction by nonparametric nonuniform intensity normalization (N3) in the analysis of trabecular bone parameters. The ability to correct for coil-induced intensity inhomogeneity was evaluated ex vivo on proximal femur specimens scanned with both a surface coil and a volume coil, which allowed for a direct evaluation of the performance of the coil correction methods without any major confounding factors. In addition, trabecular bone parameter values were correlated with values from high-resolution peripheral computed tomography (HR-pQCT) scans, and the reproducibility of trabecular bone parameters was evaluated in an in vivo study of repeat hip MR scans. The trabecular bone parameters determined from MR surface coil scans processed with the N3 coil correction method showed significant correlation (p < 0.05) with corresponding values from homogeneous intensity data in the ex vivo study. This can be compared to the correlation without coil correction (p < 0.5), and coil correction using low-pass filtering (LPF) (p < 0.53). The in vivo interscan variability was reduced from 8.9% to 12.8% using LPF-based to 3.6%-8.4% (CV) using N3 coil correction; hence the results showed that N3 is advantageous to LPF-based coil correction. No significant differences in correlation to HR-pQCT data were found for the coil correction methods. The significant correlations with volume coil data and high reproducibility of the N3 processed data imply that N3 coil correction preserve image information while accurately correcting for coil-induced intensity inhomogeneities, which makes it suitable for quantitative analysis of trabecular bone structure from MR images acquired with surface coils.