Ultrasound parametric imaging of the calcaneus: in vivo results with a new device

Mechanisms of Interaction of Ultrasound With Cancellous Bone: A Review

Ultrasound is now a clinically accepted modality in the management of osteoporosis. The most common commercial clinical devices assess fracture risk from measurements of attenuation and sound speed in cancellous bone. This review discusses fundamental mechanisms underlying the interaction between ultrasound and cancellous bone. Because of its two-phase structure (mineralized trabecular network embedded in soft tissue-marrow), its anisotropy, and its inhomogeneity, cancellous bone is more difficult to characterize than most soft tissues. Experimental data for the dependencies of attenuation, sound speed, dispersion, and scattering on ultrasound frequency, bone mineral density, composition, microstructure, and mechanical properties are presented. The relative roles of absorption, scattering, and phase cancellation in determining attenuation measurements in vitro and in vivo are delineated. Common speed of sound metrics, which entail measurements of transit times of pulse leading edges (to avoid multipath interference), are greatly influenced by attenuation, dispersion, and system properties, including center frequency and bandwidth. However, a theoretical model has been shown to be effective for correction for these confounding factors in vitro and in vivo. Theoretical and phantom models are presented to elucidate why cancellous bone exhibits negative dispersion, unlike soft tissue, which exhibits positive dispersion. Signal processing methods are presented for separating "fast" and "slow" waves (predicted by poroelasticity theory and supported in cancellous bone) even when the two waves overlap in time and frequency domains. Models to explain dependencies of scattering on frequency and mean trabecular thickness are presented and compared with measurements. Anisotropy, the effect of the fluid filler medium (marrow in vivo or water in vitro), phantoms, computational modeling of ultrasound propagation, acoustic microscopy, and nonlinear properties in cancellous bone are also discussed.

Quantitative ultrasound imaging monitoring progressive disuse osteopenia and mechanical stimulation mitigation in calcaneus region through a 90-day bed rest human study

Background

Osteoporosis parallels aging and functional mechanical unloading (e.g., space flight and bed rest), jeopardizing mineral density, microstructure, and integrity of bone and leading to an increased risk of fracture. A way to combat this deterioration is to harness the sensitivity of bone to mechanical signals.

Objective

This study evaluates the longitudinal effect of a dynamic mechanical loading through the heel on human bone in vivo during 90-day bed rest, monitored by quantitative ultrasound (QUS) imaging and dual-energy X-ray absorptiometry (DXA) in localized regions of interests, i.e., calcaneus.

Methods

A total of 29 bed rest individuals were evaluated (11 control and 18 treatment) with a brief (10-minute) daily low-intensity (0.3g), high-frequency (30Hz) dynamic mechanical stimulation countermeasure through vibrational inhibition bone erosion (VIBE). Both QUS and DXA detected longitudinal bone density and quality changes.

Results

Ultrasound velocity (UV) decreased in the control group and increased in the group treated with low-intensity loading. The UV increased by 1.9% and 1.6% at 60- and 90-day bed rest (p=0.01) in VIBE over control groups. A trend was found in broadband ultrasound attenuation (BUA), with a VIBE benefit of 1.8% at day 60 and 0.5% at day 90 in comparison with control (p=0.5). Bone mineral density (BMD) assessed by DXA decreased -4.50% for control individuals and -2.18% for VIBE individuals, showing a moderate effect of the mechanical intervention (p=0.19). Significant correlations between QUS and DXA were observed, with a combined BUA and UV vs. BMD: r²=0.70.

Conclusion

These results indicated that low-intensity, high-frequency loading has the potential to mitigate regional bone loss induced by long-term bed rest and that QUS imaging may be able to assess the subtle changes in bone alteration.

Translational potential of this article

Quantitative ultrasound has shown the efficacy of noninvasively assessing bone mass and structural properties in cadaver and isolated trabecular bone samples. While its ability in measuring in vivo bone quality and density is still unclear, a scanning confocal ultrasound imaging is developed and can perform an instant assessment for the subtle changes of such bone loss. This ultrasound imaging modality can potentially be used in the clinical assessment of bone mass. Moreover, physical stimulation has shown the ability to prevent bone loss induced by functional disuse and estrogen deficiency in animal models. However, its treatment capability is unclear. This study has shown that low-magnitude mechanical signals, introduced using low-intensity vibration (LIV), can mitigate regional bone loss caused by functional disuse. Thus localized mechanical treatment, and the quantitative ultrasound imaging have shown translational potential to noninvasively attenuate bone loss, and assess bone mass in the clinic, e.g., in an extreme condition such as long-term space mission, and long-term bedrest such as in case of spinal cord injury.

The Ability of Ultrasonic Backscatter Parametric Imaging to Characterize Bovine Trabecular Bone

The ultrasonic backscatter technique holds the promise of characterizing bone density and microstructure. This paper conducts ultrasonic backscatter parametric imaging based on measurements of apparent integrated backscatter (AIB), spectral centroid shift (SCS), frequency slope of apparent backscatter (FSAB), and frequency intercept of apparent backscatter (FIAB) for representing trabecular bone mass and microstructure. We scanned 33 bovine trabecular bone samples using a 7.5 MHz focused transducer in a 20 mm × 20 mm region of interest (ROI) with a step interval of 0.05 mm. Images based on the ultrasonic backscatter parameters (i.e., AIB, SCS, FSAB, and FIAB) were constructed to compare with photographic images of the specimens as well as two-dimensional (2D) μ-CT images from approximately the same depth and location of the specimen. Similar structures and trabecular alignments can be observed among these images. Statistical analyses demonstrated that the means and standard deviations of the ultrasonic backscatter parameters exhibited significant correlations with bone density (|R| = 0.45-0.78, p < 0.01) and bone microstructure (|R| = 0.44-0.87, p < 0.001). Some bovine trabecular bone microstructure parameters were independently associated with the ultrasonic backscatter parameters (ΔR² = 4.18%-44.45%, p < 0.05) after adjustment for bone apparent density (BAD). The results show that ultrasonic backscatter parametric imaging can provide a direct view of the trabecular microstructure and can reflect information about the density and microstructure of trabecular bone.

Ultrasonic scanner for in vivo measurement of cancellous bone properties from backscattered data

A dedicated ultrasonic scanner for acquiring RF echoes backscattered from the trabecular bone was developed. The design of device is based on the goal of minimizing of custom electronics and computations executed solely on the main computer processor and the graphics card. The electronic encoder-digitizer module executing all of the transmission and reception functions is based on a single low-cost field programmable gate array (FPGA). The scanner is equipped with a mechanical sector-scan probe with a concave transducer with 50 mm focal length, center frequency of 1.5 MHz and 60% bandwidth at -6 dB. The example of femoral neck bone examination shows that the scanner can provide ultrasonic data from deeply located bones with the ultrasound penetrating the trabecular bone up to a depth of 20 mm. It is also shown that the RF echo data acquired with the scanner allow for the estimation of attenuation coefficient and frequency dependence of backscattering coefficient of trabecular bone. The values of the calculated parameters are in the range of corresponding in vitro data from the literature but their variation is relatively high.

Numerical analysis of variability in ultrasound propagation properties induced by trabecular microstructure in cancellous bone

The manner by which the trabecular microstructure affects the propagation of ultrasound waves through cancellous bone is numerically investigated by finite difference time-domain (FDTD) simulation. Sixteen 3-D numerical models of 6.45 x 6.45 x 6.45 mm with a voxel size of 64.5 microm are reconstructed using a 3-D microcomputed tomographic (microCT) image taken from a bovine cancellous bone specimen of approximately 20 x 20 x 9 mm. All cancellous bone models have an oriented trabecular structure, and their trabecular elements are gradually eroded to increase the porosity using an image processing technique. Three erosion procedures are presented to realize various changes in the trabecular microstructure with increasing porosity. FDTD simulations of the ultrasound pulse waves propagating through the cancellous bone models at each eroded step are performed in 2 cases of the propagations parallel and perpendicular to the major trabecular orientation. The propagation properties of the wave amplitudes and propagation speeds are derived as a function of the porosity, and their variability due to the trabecular microstructure is revealed. To elucidate an effect of the microstructure, the mean intercept length (MIL), which is a microstructural parameter, is introduced, and the correlations of the propagation properties with the MILs of the trabecular elements and pore spaces are investigated.

The application of quantitative ultrasound (QUS) on study of aging effects of achilles tendons

Recently, the quantitative ultrasound (QUS) technique has been widely applied for the characterization of tissues. For example, it can be used for the quantization of Achilles tendons based on the broadband ultrasound attenuation (BUA) and speed of sound (SOS) when the ultrasound wave passes through the tissues. The main purpose of this study is to develop an integrated system to investigate the correlations between the measurements of Achilles tendons using QUS (UBIS 5000) and 2D ultrasound instrument (HDI 5000). Subjects including young (32 females and 17 males; mean age: 23.7+/-2.0) and middle-age groups (8 female and 8 males; mean age: 47.3+/-8.5 s) were recruited and tested for this study. All the subjects were ensured to have no habit of regular exercise and no record of tendon injury. The results show that BUA is significantly higher for the young group (45.2+/-1.6 dB/MHz) than the middle-age group (40.5+/-1.9 dB/MHz), while SOS is significantly lower for the young (1601.9+/-11.2m/s) than the middle-age (1624.1+/-8.7m/s). The thickness of Achilles tendons for both groups (young: 4.31+/-0.23mm; middle age: 4.24+/-0.23mm) are very similar. The noninvasive ultrasonic assessment of Achilles tendons might be helpful in clinical diagnosis and in evaluation of a therapeutic regimen. In the future, further animal studies should be made by measuring the BUA and SOS invasively to evaluate the accuracy of the parameters measured using QUS parametric images.

Ultrasound emitter localization in heterogeneous media

A novel algorithm to accurately determine the location of an ultrasound source within heterogeneous media is presented. The method obtains a small spacial error of 748 microm+/-310 microm for 100 different measurements inside a circular area with 140 mm diameter. The new algorithm can be used in targeted drug delivery for cancer therapies as well as to accurately locate any kind of ultrasound sources in heterogeneous media, such as ultrasonically marked medical devices or tumors.

Broadband ultrasound attenuation in the calcaneal region: a comparative study of single-position versus scanning systems

This work describes a system developed to measure the broadband ultrasound attenuation (BUA) in the calcaneal region. The patient's calcanei were inspected using a microcomputer-controlled X-Y axis displacement unit with two 500-kHz, central-frequency, ultrasound transducers. The transducers facing each other are submerged in a small water tank with a support for the patient's foot between them. The system allows data to be collected from a single position or by scanning the calcaneal region to obtain a BUA map. Tests were carried out on 201 patients (110 using the single-position method, and 91 using the scanning method). The results were compared with those of densitometry tests performed using the dual energy X-ray absorptiometry (DEXA) technique (single position: r=0.50; P<0.0001; scanner: r=0.75; P<0.0001). It was concluded that the single position method is more susceptible to errors due to the difficulty in positioning the transducers relative to the calcaneus. The scanning method provides better results and can be used to screen patients before referring them for DEXA.

Instrumentation for in vivo ultrasonic characterization of bone strength

Although it has been more than 20 years since the first recorded use of a quantitative ultrasound (QUS) technology to predict bone fragility, the field has not yet reached its maturity. QUS has the potential to predict fracture risk in several clinical circumstances and has the advantages of being nonionizing, inexpensive, portable, highly acceptable to patients, and repeatable. However, the wide dissemination of QUS in clinical practice is still limited and suffering from the absence of clinical consensus on how to integrate QUS technologies in bone densitometry armamentarium. Several critical issues need to be addressed to develop the role of QUS within rheumatology. These include issues of technologies adapted to measure the central skeleton, data acquisition, and signal processing procedures to reveal bone properties beyond bone mineral quantity and elucidation of the complex interaction between ultrasound and bone structure. This article reviews the state-of-the art in technological developments applied to assess bone strength in vivo. We describe generic measurement and signal processing methods implemented in clinical ultrasound devices, the devices and their practical use, and performance measures. The article also points out the present limitations, especially those related to the absence of standardization, and the lack of comprehensive theoretical models. We conclude with suggestions of future lines and trends in technology challenges and research areas such as new acquisition modes, advanced signal processing techniques, and modelization.

Bone surface topology mapping and its role in trabecular bone quality assessment using scanning confocal ultrasound

INTRODUCTION

Quantitative ultrasound (QUS) has been used to assess non-invasively bone quality, in which ultrasound velocity (UV) is a primary acoustic property.

METHODS

While UV calculation requires the tissue thickness in the ultrasound path, a bone surface topology mapping (STM) method was developed in this study for enhancing the accuracy of the UV measurement. STM accuracy was verified by both aluminum and a QUS heel phantom, indicating that the STM can determine the phantom thickness within 0.02 mm thickness error and the aluminum calibration step within 0.1 mm thickness error. STM performance was further evaluated using 25 cadaveric human calcanei samples.

RESULTS

The UV calculations using STM had a significant better correlation to bone mineral density (BMD) (r = 0.75, p < 0.05), volume fraction (r = 0.72, p < 0.05) and modulus (r = 0.69, p < 0.05) than the UV with fixed thickness. The later correlation coefficients were r = 0.64 for BMD, r = 0.65 for volume fraction, and r = 0.58 for modulus. The nBUA value determined using STM was also highly correlated to BMD (r(2) = 0.74) and modulus (r(2) = 0.62). This was comparable to the correlation result for BUA (BMD: r(2) = 0.76; Modulus: r(2) = 0.64).

CONCLUSION

These results suggested that STM technique in scanning ultrasound is capable of determining calcaneus bone thickness and hence enhancing the accuracy of UV measurement.

A new anisotropy index on trabecular bone radiographic images using the fast Fourier transform

BACKGROUND

The degree of anisotropy (DA) on radiographs is related to bone structure, we present a new index to assess DA.

METHODS

In a region of interest from calcaneus radiographs, we applied a Fast Fourier Transform (FFT). All the FFT spectra involve the horizontal and vertical components corresponding respectively to longitudinal and transversal trabeculae. By visual inspection, we measured the spreading angles: Dispersion Longitudinal Index (DLI) and Dispersion Transverse Index (DTI) and calculated DA = 180/(DLI+DTI). To test the reliability of DA assessment, we synthesized images simulating radiological projections of periodic structures with elements more or less disoriented.

RESULTS

Firstly, we tested synthetic images which comprised a large variety of structures from highly anisotropic structure to the almost isotropic, DA was ranging from 1.3 to 3.8 respectively. The analysis of the FFT spectra was performed by two observers, the Coefficients of Variation were 1.5% and 3.1 % for intra-and inter-observer reproducibility, respectively. In 22 post-menopausal women with osteoporotic fracture cases and 44 age-matched controls, DA values were respectively 1.87 +/- 0.15 versus 1.72 +/- 0.18 (p = 0.001). From the ROC analysis, the Area Under Curve (AUC) were respectively 0.65, 0.62, 0.64, 0.77 for lumbar spine, femoral neck, total femoral BMD and DA.

CONCLUSION

The highest DA values in fracture cases suggest that the structure is more anisotropic in osteoporosis due to preferential deletion of trabeculae in some directions.

Comparison of imaging-guided and non-imaging-guided quantitative sonography of the calcaneus with dual X-ray absorptiometry of the spine and femur

OBJECTIVE

The purpose of our study was to evaluate the diagnostic agreement between imaging-guided and non-imaging-guided quantitative sonography of the calcaneus and dual X-ray absorptiometry of the spine and femur to show osteoporosis.

SUBJECTS AND METHODS

In 113 patients (73 women, 59 +/- 14 years old; 40 men, 48 +/- 16 years old), dual X-ray absorptiometry of the lumbar spine and the proximal femur, imaging-guided quantitative sonography, and non-imaging-guided quantitative sonography of the calcaneus were performed. The percentage of patients having a T-score equal to or less than a threshold of -2.5 SDs (prevalence of osteoporosis) was calculated for each imaging technique. The diagnostic agreement of the three techniques in identifying individuals with osteoporosis was assessed.

RESULTS

Eleven percent of the women and 8% of the men were classified as osteoporotic by imaging-guided quantitative sonography, and 38% of the women and 25% of the men were so classified by non-imaging-guided quantitative sonography. At dual X-ray absorptiometry of the spine, 44% of the women and 38% of the men were classified as osteoporotic, and, at different femoral regions, 19-60% of the women and 8-38% of the men were so classified. Kappa analysis for both quantitative sonography techniques was not significant. Kappa analysis for both quantitative sonography techniques and dual X-ray absorptiometry showed diagnostic agreement to be generally poor.

CONCLUSION

No advantage in diagnostic accuracy could be found for imaging-guided quantitative sonography. The considerable diagnostic disagreement between both quantitative sonography techniques and dual X-ray absorptiometry could be confusing in daily clinical practice.

Diagnostic disagreement of imaging quantitative sonography of the calcaneus with dual X-ray absorptiometry of the spine and femur

OBJECTIVE

This study evaluates the diagnostic agreement between imaging quantitative sonography of the calcaneus and dual X-ray absorptiometry of the spine and femur for diagnosing osteoporosis.

MATERIALS AND METHODS

In 498 female patients (56 +/- 18 years old), bone mineral density measurements by dual X-ray absorptiometry of the lumbar spine (posteroanterior, L1--L4) and the proximal femur and imaging quantitative sonography of the calcaneus were performed. The percentage of patients having T-scores less than or equal to a threshold of -2.5 standard deviations below a young normal reference was used to compare quantitative sonography with dual X-ray absorptiometry. The diagnostic agreement was assessed using kappa scores.

RESULTS

Approximately 30% of the patients had a T-score less than or equal to -2.5 standard deviations as assessed by imaging quantitative sonography (broadband ultrasound attenuation), 26.5% as assessed by dual X-ray absorptiometry of the spine, and 16.7--56.4% as assessed by dual X-ray absorptiometry of the different regions of interest at the femur. Kappa analysis showed that severe diagnostic disagreement exists among broadband ultrasound attenuation and dual X-ray absorptiometry (kappa = 0.28-0.42).

CONCLUSION

Considerable diagnostic disagreement exists between imaging quantitative sonography and dual X-ray absorptiometry of the spine and femur. The disagreement is in the same range as that reported recently in comparisons of dual X-ray absorptiometry and nonimaging quantitative sonography. In general, no distinct advantage for imaging quantitative sonography could be found when compared with other techniques.

Determination of a standard site for the measurement of bone mineral density of the human calcaneus

Ultrasound of the calcaneus may be used as a cheap, ionising radiation-free and easy to use indicator of skeletal status, and hence of osteoporotic fracture risk. At present ultrasound is not widely used as it suffers from high precision errors. As ultrasound parameters are determined in part by bone mineral density (BMD), an increase in the accuracy and precision of BMD measurements should reduce the precision error associated with ultrasound measurements. The aim of this study was to define an anatomical site on the calcaneus at which accurate and precise measurements of BMD can be made. Ten dry calcanei and 10 cadaveric feet were scanned using a DXA scanner; 9 anatomically defined regions (1 cm2) were selected in the posterior part of the calcaneus for analysis. The centre of region 1 was positioned halfway along the line joining the anterior border of the calcaneal tubercle and the peak of the posterior superior tubercle, and the remaining 8 regions were placed around this central area. The BMD in these 9 regions was compared with the whole bone BMD and the variability of BMD within each of the 9 regions was measured. The reproducibility of the technique was assessed by taking 10 repeated measurements of 2 bone and 2 cadaveric specimens, each specimen being removed and repositioned between measurements. Region 1 was found to be the most representative of total BMD in cadaveric feet. This region also showed the least variability of BMD and consistently gave the lowest coefficients of variation in the reproducibility study both in the bone and the cadaveric specimens. This region is hence the most suitable site on the calcaneus for measuring absolute values of and changes in BMD. The surface position of region 1 was found to be consistently 5/9 along the line at 45 degrees to the vertical, from the lateral malleolus to the heel. The identification of the surface location of region 1 relative to anatomical landmarks of the foot has enabled the same anatomical site to be measured in all subjects. This allows meaningful intersubject comparisons to be made. Preliminary data suggest that precision errors using ultrasound are also reduced when measurements are taken at this region of the calcaneus. The reduction in the precision error of ultrasound assessment of skeletal status may provide a cheap and safe way to identify individuals at risk from osteoporotic fracture.

Automatic detection of the boundary of the calcaneus from ultrasound parametric images using an active contour model; clinical assessment

This paper presents a computerized method for automated detection of the boundary of the os calcis on in vivo ultrasound parametric images, using an active dynamic contour model. The initial contour, defined without user interaction, is an iso-contour extracted from the textural feature space. The contour is deformed through the action of internal and external forces, until stability is reached. The external forces, which characterize image features, are a combination of gray-level information and second-order textural features arising from local cooccurrence matrices. The broadband ultrasound attenuation (BUA) value is then averaged within the contour obtained. The method was applied to 381 clinical images. The contour was correctly detected in the great majority of the cases. For the short-term reproducibility study, the mean coefficient of variation was equal to 1.81% for BUA values and 4.95% for areas in the detected region. Women with osteoporosis had a lower BUA than age-matched controls (p = 0.0005). In healthy women, the age-related decline was -0.45 dB/MHz/yr. In the group of healthy post-menopausal women, years since menopause, weight and age were significant predictors of BUA. These results are comparable to those obtained when averaging BUA values in a small region of interest.

Assessment of the relationship between broadband ultrasound attenuation and bone mineral density at the calcaneus using BUA imaging and DXA

The purposes of this study was to determine the relationship between broadband ultrasound attenuation (BUA) and bone mineral density (BMD) measured at different regions of the calcaneus with identical site-matched regions of interest (ROIs). Dual-energy X-ray absorptiometry (DXA) measurements of the calcaneus and BUA imaging were performed in 30 women (15 premenopausal and 15 postmenopausal). Four square ROIs were located in the great tuberosity and one square ROI in the foramen calcaneus. A ROI adapted to the shape and size of the whole calcaneus was also considered. All ROIs were analyzed three times with both techniques to minimize intra-observer variability. The correlation coefficient between attenuation and frequency was used as an index of BUA measurement error. Before accepting a measurement of BUA in inhomogeneous material, it could be useful to map the spatial variations of the measurement error. In all ROIs we found the BUA and BMD were strongly related (r = 0.78-0.91, p < 0.001). The correlation between BUA and BMD was slightly higher in the inferior part of the posterior tuberosity than in the superior part and in the foramen calcaneus. The very high correlation between attenuation and frequency found in all ROIs (r = 0.99) suggests that measurement errors of propagation were probably not significant. Ultrasound imaging yields the opportunity for studying the spatial acoustic properties in the calcaneus and their relation to bone mass or structural parameters provided by independent imaging techniques. BUA measured with current transmission techniques reflects mainly bone mass, and microarchitecture to a smaller extent.