Imaging Internal Structure of Long Bones Using Wave Scattering Theory

Femoral neck phantom imaging using time-domain topological energy method

Ultrasonic bone imaging is a complex task, primarily because of the low energy contained in the signals reflected from the internal bone structures. In this study, the reconstruction of a bone-mimicking phantom echographic image using time-domain topological energy (TDTE) is proposed. A TDTE image results from a combination of forward and adjoint fields. The first is a solution of a numerical model that reproduces the setup of the experimental data acquisition to the best extent possible. The second has similar characteristics, but the source term is the time-reversed residue between the forward field and signals obtained from the experiment. The acquisition-reconstruction system used a linear phased-array transducer with a 5 MHz center frequency to acquire the signals and was coupled with a k-wave toolbox to implement the numerical models and perform the image reconstruction. The results showed good agreement between the geometry of the real phantom and the ultrasonic images. However, thickness evaluation errors were observed, which may be due to incorrect assumptions about the velocity models throughout the medium, a priori assumed to be known. Thus, this method has shown promising results and should be applied to the real femoral neck as a long-term objective.

Estimation of Thickness and Speed of Sound for Transverse Cortical Bone Imaging Using Phase Aberration Correction Methods: An In Silico and Ex Vivo Validation Study

Delay-and-sum (DAS) beamforming of backscattered echoes is used for conventional ultrasound imaging. Although DAS beamforming is well suited for imaging in soft tissues, refraction, scattering, and absorption, porous mineralized tissues cause phase aberrations of reflected echoes and subsequent image degradation. The recently developed refraction corrected multi-focus technique uses subsequent focusing of waves at variable depths, the tracking of travel times of waves reflected from outer and inner cortical bone interfaces, the estimation of the shift needed to focus from one interface to another to determine cortical thickness (Ct.Th), and the speed of sound propagating in a radial bone direction (Ct.ν11). The method was validated previously in silico and ex vivo on plate shaped samples. The aim of this study was to correct phase aberration caused by bone geometry (i.e., curvature and tilt with respect to the transducer array) and intracortical pores for the multi-focus approach. The phase aberration correction methods are based on time delay estimation via bone geometry differences to flat bone plates and via the autocorrelation and cross correlation of the reflected ultrasound waves from the endosteal bone interface. We evaluate the multi-focus approach by incorporating the phase aberration correction methods by numerical simulation and one experiment on a human tibia bone, and analyze the precision and accuracy of measuring Ct.Th and Ct.ν11. Site-matched reference values of the cortical thickness of the human tibia bone were obtained from high-resolution peripheral computed tomography. The phase aberration correction methods resulted in a more precise (coefficient of variation of 5.7%) and accurate (root mean square error of 6.3%) estimation of Ct.Th, and a more precise (9.8%) and accurate (3.4%) Ct.ν11 estimation, than without any phase aberration correction. The developed multi-focus method including phase aberration corrections provides local estimations of both cortical thickness and sound velocity and is proposed as a biomarker of cortical bone quality with high clinical potential for the prevention of osteoporotic fractures.

Index-Rotated Fast Ultrasound Imaging of Cortical Bone Based on Predicted Velocity Model

Due to the significant acoustic impedance contrast at cortical boundaries, highly inside attenuation, and the unknown sound velocity distribution, accurate ultrasound cortical bone imaging remains a challenge, especially for the traditional pulse-echo modalities using unique sound velocity. Moreover, the large amounts of data recorded by multielement probe results in a relatively time-consuming reconstruction process. To overcome these limitations, this article proposed an index-rotated fast ultrasound imaging method based on predicted velocity model (IR-FUI-VP) for cortical cross section ultrasound tomography (UST) imaging, utilizing ray-tracing synthetic aperture (RTSA). In virtue of ring probe, the sound velocity model was predicted in advance using bent-ray inversion (BRI). With the predicted velocity model, index-rotated fast ultrasound imaging (IR-FUI) was further applied to image the cortical cross sections in the sectors corresponding to the dynamic apertures (DAs) and ring center. The final result was merged by all sector images. One cortical bone phantom and two ex vivo bovine femurs were utilized to demonstrate the performance of the proposed method. Compared to the conventional synthetic aperture (SA) imaging, the method can not only accurately image the outer cortical boundary but also precisely reconstruct the inner cortical surface. The mean relative errors of the predicted sound velocity in the region of interest (ROI) were all smaller than 7%, and the mean errors of cortical thickness are all less than 0.31 mm. The reconstructed images of bovine femurs were in good agreement with the reference images scanned by micro-computed tomography ( μ CT) with respect to the morphology and thickness. The speed of IR-FUI is about 3.73 times faster than the traditional SA. It is proved that the proposed IR-FUI-VP-based UST is an effective way for fast and accurate cortical bone imaging.

Estimation of Thickness and Speed of Sound in Cortical Bone Using Multifocus Pulse-Echo Ultrasound

Most bone loss during the development of osteoporosis occurs in cortical bone at the peripheral skeleton. Decreased cortical thickness (Ct.Th) and the prevalence of large pores at the tibia are associated with reduced bone strength at the hip. Ct.Th and cortical sound velocity, i.e., a surrogate marker for changes of cortical porosity (Ct.Po), are key biomarkers for the identification of patients at high fracture risk. In this study, we have developed a method using a conventional ultrasound array transducer to determine thickness (Ct.Th) and the compressional sound velocity propagating in the radial bone direction (Ct. ν₁₁ ) using a refraction-corrected multifocus imaging approach. The method was validated in-silico on porous bone plate models using a 2-D finite-difference time-domain method and ex vivo on plate-shaped plastic reference materials and on plate-shaped cortical bovine tibia samples. Plane-wave pulse-echo measurements provided reference values to assess precision and accuracy of our method. In-silico results revealed the necessity to account for inclination-dependent transmission losses at the bone surface. Moreover, the dependence of Ct. ν₁₁ on both porosity and pore density was observed. Ct.Th and Ct. ν₁₁ obtained ex vivo showed a high correlation ) with reference values. The ex-vivo accuracy and precision for Ct. ν₁₁ were 29.9 m/s and 0.94%, respectively, and those for Ct.Th were 0.04 mm and 1.09%, respectively. In conclusion, this numerical and experimental study demonstrates an accurate and precise estimation of Ct.Th and Ct. ν₁₁ . The developed multifocus technique may have high clinical potential to improve fracture risk prediction using noninvasive and nonionizing conventional ultrasound technology with image guidance.

Nonlinear Inversion of Ultrasonic Dispersion Curves for Cortical Bone Thickness and Elastic Velocities

In this study, a nonlinear grid-search inversion has been developed to estimate the thickness and elastic velocities of long cortical bones, which are important determinants of bone strength, from axially-transmitted ultrasonic data. The inversion scheme is formulated in the dispersive frequency-phase velocity domain to recover bone properties. The method uses ultrasonic guided waves to retrieve overlying soft tissue thickness, cortical thickness, compressional, and shear-wave velocities of the cortex. The inversion strategy requires systematic examination of a large set of trial dispersion-curve solutions within a pre-defined model space to match the data with minimum cost in a least-squares sense. The theoretical dispersion curves required to solve the inverse problem are computed for bilayered bone models using a semi-analytical finite-element method. The feasibility of the proposed approach was demonstrated by the numerically simulated data for a 1 mm soft tissue-5 mm bone bilayer and ex-vivo data from a bovine femur plate with an overlying 2 mm-thick soft-tissue mimic. The bootstrap method was employed to evaluate the inversion uncertainty and stability. Our results have shown that the cortical thickness and wave speeds could be recovered with fair accuracy.