Measuring anisotropy of elastic wave velocity with ultrasound imaging and an autofocus method: application to cortical bone

Assessing Transducer Parameters for Accurate Medium Sound Speed Estimation and Image Reconstruction

The influence of the transducer lens on image reconstruction is often overlooked. Lenses usually exhibit a lower sound speed than soft biological tissues. In academic research, the exact lens sound speed and thickness are typically unknown. Here, we present a simple and nondestructive method to characterize the lens sound speed and thickness as well as the time to peak of the round-trip ultrasound waveform, another key parameter for optimal image reconstruction. We applied our method to three transducers with center frequencies of 2.5, 7.5, and 15 MHz. We estimated the three parameters with an element-by-element transmission sequence that records internal reflections within the lens. We validated the retrieved parameters using an autofocusing approach that estimates sound speed in water. We show that the combination of our parameters estimation method with two-layer ray tracing outperforms standard image reconstruction. For all transducers, we successfully improved the accuracy of medium sound speed estimation, spatial resolution, and contrast. The proposed method is simple and robust and provides an accurate estimation of the transducer lens parameters and the time to peak of the ultrasound waveform, which leads to improved ultrasound image quality.

Bi-Directional Axial Transmission measurements applied in a clinical environment

Accurate measurement of cortical bone parameters may improve fracture risk assessment and help clinicians on the best treatment strategy. Patients at risk of fracture are currently detected using the current X-Ray gold standard DXA (Dual XRay Absorptiometry). Different alternatives, such as 3D X-Rays, Magnetic Resonance Imaging or Quantitative Ultrasound (QUS) devices, have been proposed, the latter having advantages of being portable and sensitive to mechanical and geometrical properties. The objective of this cross-sectional study was to evaluate the performance of a Bi-Directional Axial Transmission (BDAT) device used by trained operators in a clinical environment with older subjects. The device, positioned at one-third distal radius, provides two velocities: VFAS (first arriving signal) and VA0 (first anti-symmetrical guided mode). Moreover, two parameters are obtained from an inverse approach: Ct.Th (cortical thickness) and Ct.Po (cortical porosity), along with their ratio Ct.Po/Ct.Th. The areal bone mineral density (aBMD) was obtained using DXA at the femur and spine. One hundred and six patients (81 women, 25 men) from Marien Hospital and St. Anna Hospital (Herne, Germany) were included in this study. Age ranged from 41 to 95 years, while body mass index (BMI) ranged from 16 to 47 kg.m-2. Three groups were considered: 79 non-fractured patients (NF, 75±13years), 27 with non-traumatic fractures (F, 80±9years) including 14 patients with non-vertebral fractures (NVF, 84±7years). Weak to moderate significant Spearman correlations (R ranging from 0.23 to 0.53, p < 0.05) were found between ultrasound parameters and age, BMI. Using multivariate Partial Least Square discrimination analyses with Leave-One-Out Cross-Validation (PLS-LOOCV), we found the combination of VFAS and the ratio Ct.Po/Ct.Th to be predictive for all non traumatic fractures (F) with the odds ratio (OR) equals to 2.5 [1.6-3.4] and the area under the ROC curve (AUC) equal to 0.63 [0.62-0.65]. For the group NVF, combination of four parameters VA0. Ct.Th, Ct.Po and Ct.Po/Ct.Po, along with age provides a discrimination model with OR and AUC equals to 7.5 [6.0-9.1] and 0.75 [0.73-0.76]. When restricted to a smaller population (87 patients) common to both BDAT and DXA, BDAT ORs and AUCs are comparable or slightly higher to values obtained with DXA. The fracture risk assessment by BDAT method in older patients, in a clinical setting, suggests the benefit of the affordable and transportable device for the routine use.

Fabrication and mechanical properties of different types of carbon fiber reinforced polyetheretherketone: A comparative study

OBJECTIVES

To find alternative non-metallic materials as dental implants for clinical application, different types of carbon fiber reinforced polyetheretherketone were fabricated and investigated.

METHODS

Continuous carbon fiber reinforced polyetheretherketone fabrics were fabricated with polyetheretherketone fibers and carbon fibers. Different kinds of carbon fiber reinforced polyetheretherketone were synthesized by setting specific experiment parameters of injection or hot press molding. Various mechanical tests were performed to determine the mechanical properties of different carbon fiber reinforced polyetheretherketone, pure polyetheretherketone and pure titanium.

RESULTS

Polyetheretherketone composites presented outstanding mechanical and thermal properties after incorporating carbon fiber. The bending and tensile strength of short carbon fiber reinforced polyetheretherketone were close to human bone, and the bending strength of continuous carbon fiber reinforced polyetheretherketone reached 644 MPa, even higher than that of pure titanium.

CONCLUSIONS

The mechanical properties of polyetheretherketone composites are more similar to bone tissue than titanium, and the stress shielding phenomenon may be inhibited. They may become promising materials as substitutions for titanium and prospective materials in bone tissue engineering.

Accelerated 2-D Real-Time Refraction-Corrected Transcranial Ultrasound Imaging

In a recent study, we proposed a technique to correct aberration caused by the skull and reconstruct a transcranial B-mode image with a refraction-corrected synthetic aperture imaging (SAI) scheme. Given a sound speed map, the arrival times were calculated using a fast marching technique (FMT), which solves the Eikonal equation and, therefore, is computationally expensive for real-time imaging. In this article, we introduce a two-point ray tracing method, based on Fermat's principle, for fast calculation of the travel times in the presence of a layered aberrator in front of the ultrasound probe. The ray tracing method along with the reconstruction technique is implemented on a graphical processing unite (GPU). The point spread function (PSF) in a wire phantom image reconstructed with the FMT and the GPU implementation was studied with numerical synthetic data and experiments with a bone-mimicking plate and a sagittally cut human skull. The numerical analysis showed that the error on travel times is less than 10% of the ultrasound temporal period at 2.5 MHz. As a result, the lateral resolution was not significantly degraded compared with images reconstructed with FMT-calculated travel times. The results using the synthetic, bone-mimicking plate, and skull dataset showed that the GPU implementation causes a lateral/axial localization error of 0.10/0.20, 0.15/0.13, and 0.26/0.32 mm compared with a reference measurement (no aberrator in front of the ultrasound probe), respectively. For an imaging depth of 70 mm, the proposed GPU implementation allows reconstructing 19 frames/s with full synthetic aperture (96 transmission events) and 32 frames/s with multiangle plane wave imaging schemes (with 11 steering angles) for a pixel size of [Formula: see text]. Finally, refraction-corrected power Doppler imaging is demonstrated with a string phantom and a bone-mimicking plate placed between the probe and the moving string. The proposed approach achieves a suitable frame rate for clinical scanning while maintaining the image quality.

Improvement of Patient Classification Using Feature Selection Applied to Bidirectional Axial Transmission

Osteoporosis is still a worldwide problem, particularly due to associated fragility fractures. Patients at risk of fracture are currently detected using the X-Ray gold standard dual-energy X-ray absorptiometry (DXA), based on a calibrated 2-D image. Different alternatives, such as 3-D X-rays, magnetic resonance imaging (MRI) or ultrasound, have been proposed, the latter having advantages of being portable and sensitive to mechanical and geometrical properties. Bidirectional axial transmission (BDAT) has been used to classify between patients with or without nontraumatic fractures using "classical" ultrasonic parameters, such as velocities, as well as cortical thickness and porosity, obtained from an inverse problems. Recently, complementary parameters acquired with structural and textural analysis of guided wave spectrum images (GWSIs) have been introduced. These parameters are not limited by solution ambiguities, as for inverse problem. The aim of the study is to improve the patient classification using a feature selection strategy for all available ultrasound features completed by clinical parameters. To this end, three classical feature ranking methods were considered: analysis of variance (ANOVA), recursive feature elimination (RFE), and extreme gradient boosting importance feature (XGBI). In order to evaluate the performance of the feature selection techniques, three classical classification methods were used: logistic regression (LR), support vector machine (SVM), and extreme gradient boosting (XGB). The database was obtained from a previous clinical study [Minonzio et al., 2019]. Results indicate that the best accuracy of 71 [66-76]% was achieved by using RFE and SVM with 22 (out of 43) ultrasonic and clinical features. This value outperformed the accuracy of 68 [64-73]% reached with 2 (out of 6) DXA and clinical features. These values open promising perspectives toward improved and generalizable classification of patients at risk of fracture.

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.

Refraction-Corrected Transcranial Ultrasound Imaging Through the Human Temporal Window Using a Single Probe

Transcranial ultrasound imaging (TUI) is a diagnostic modality with numerous applications, but unfortunately, it is hindered by phase aberration caused by the skull. In this article, we propose to reconstruct a transcranial B-mode image with a refraction-corrected synthetic aperture imaging (SAI) scheme. First, the compressional sound velocity of the aberrator (i.e., the skull) is estimated using the bidirectional headwave technique. The medium is described with four layers (i.e., lens, water, skull, and water), and a fast marching method calculates the travel times between individual array elements and image pixels. Finally, a delay-and-sum algorithm is used for image reconstruction with coherent compounding. The point spread function (PSF) in a wire phantom image and reconstructed with the conventional technique (using a constant sound speed throughout the medium), and the proposed method was quantified with numerical synthetic data and experiments with a bone-mimicking plate and a human skull, compared with the PSF achieved in a ground truth image of the medium without the aberrator (i.e., the bone plate or skull). A phased-array transducer (P4-1, ATL/Philips, 2.5 MHz, 96 elements, pitch = 0.295 mm) was used for the experiments. The results with the synthetic signals, the bone-mimicking plate, and the skull indicated that the proposed method reconstructs the scatterers with an average lateral/axial localization error of 0.06/0.14 mm, 0.11/0.13 mm, and 1.0/0.32 mm, respectively. With the human skull, an average contrast ratio (CR) and full-width-half-maximum (FWHM) of 37.1 dB and 1.75 mm were obtained with the proposed approach, respectively. This corresponds to an improvement of CR and FWHM by 7.1 dB and 36% compared with the conventional method, respectively. These numbers were 12.7 dB and 41% with the bone-mimicking plate.

Revealing Intraosseous Blood Flow in the Human Tibia With Ultrasound

Intraosseous blood circulation is thought to have a critical role in bone growth and remodeling, fracture healing, and bone disorders. However, it is rarely considered in clinical practice because of the absence of a suitable noninvasive in vivo measurement technique. In this work, we assessed blood perfusion in tibial cortical bone simultaneously with blood flow in the superficial femoral artery with ultrasound imaging in five healthy volunteers. After suppression of stationary signal with singular‐value‐decomposition, pulsatile blood flow in cortical bone tissue is revealed, following the heart rate measured in the femoral artery. Using a method combining transverse oscillations and phase‐based motion estimation, 2D vector flow was obtained in the cortex of the tibia. After spatial averaging over the cortex, the peak blood velocity along the long axis of the tibia was measured at four times larger than the peak blood velocity across the bone cortex. This suggests that blood flow in central (Haversian) canals is larger than in perforating (Volkmann's) canals, as expected from the intracortical vascular organization in humans. The peak blood velocity indicates a flow from the endosteum to the periosteum and from the heart to the foot for all subjects. Because aging and the development of bone disorders are thought to modify the direction and velocity of intracortical blood flow, their quantification is crucial. This work reports for the first time an in vivo quantification of the direction and velocity of blood flow in human cortical bone. © 2021 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

In vivopulse-echo measurement of apparent broadband attenuation andQfactor in cortical bone: a preliminary study

Quantitative ultrasound (QUS) methods have been introduced to assess cortical bone health at the radius and tibia through the assessment of cortical thickness (Ct.Th), cortical porosity and bulk wave velocities. Ultrasonic attenuation is another QUS parameter which is not currently used. We assessed the feasibility ofin vivomeasurement of ultrasonic attenuation in cortical bone with a broadband transducer with 3.5 MHz center frequency. Echoes from the periosteal and endosteal interfaces were fitted with Gaussian pulses using sparse signal processing. Then, the slope of the broadband ultrasonic attenuation (Ct.nBUA) in cortical bone and quality factorQ11-1were calculated with a parametric approach based on the center-frequency shift. Five human subjects were measured at the one-third distal radius with pulse-echo ultrasound, and reference data was obtained with high-resolution x-ray peripheral computed tomography (Ct.Th and cortical volumetric bone mineral density (Ct.vBMD)). Ct.Th was used in the calculation of Ct.nBUA whileQ11-1is obtained solely from ultrasound data. The values of Ct.nBUA (6.7 ± 2.2 dB MHz^-1.cm^-1) andQ11-1(8.6 ± 3.1%) were consistent with the literature data and were correlated to Ct.vBMD (R2=0.92,p<0.01, RMSE = 0.56 dB.MHz^-1.cm^-1, andR2=0.93,p<0.01, RMSE = 0.76%). This preliminary study suggests that the attenuation of an ultrasound signal propagating in cortical bone can be measuredin vivoat the one-third distal radius and that it provides an information on bone quality as attenuation values were correlated to Ct.vBMD. It remains to ascertain that Ct.nBUA andQ11-1measured here exactly reflect the true (intrinsic) ultrasonic attenuation in cortical bone. Measurement of attenuation may be considered useful for assessing bone health combined with the measurement of Ct.Th, porosity and bulk wave velocities in multimodal cortical bone QUS methods.