An investigation of pulse-timing techniques for broadband ultrasonic velocity determination in cancellous bone: a simulation study

Ultrasound Imaging of Hip Displacement in Children With Cerebral Palsy

OBJECTIVE

An approach to estimation of hip displacement on ultrasound (US) images is described. Its accuracy is validated through numerical simulation, an in vitro study with 3-D-printed hip phantoms and pilot in vivo data.

METHODS

A diagnostic index, migration percentage (MP), is defined by the ratio of acetabulum-femoral head distance to femoral head width. The acetabulum-femoral head distance could be measured directly on hip US images, while the femoral head width was estimated from the diameter of a best-fit circle. Simulation was performed to evaluate the accuracy of circle fitting with noiseless and noisy data. Surface roughness was also considered. Nine hip phantoms (three different sizes of femur head × three MP values) and 10 US hip images were used in this study.

RESULTS

The maximum diameter error was 16.1 ± 8.5% when the roughness and noise were 20% of the original radius and 20% of the wavelet peak, respectively. In the phantom study, the percentage errors of MPs between the 3-D-design US and X-ray US were 0.3%-6.6% and 0.0%-5.7%, respectively. From the pilot clinical trial, the mean absolute difference between the X-ray-US MPs was 3.5 ± 2.8% (1%-9%).

CONCLUSION

This study indicates that the US method can be used to evaluate hip displacement in children.

Investigation of ultrasonic soft tissue-bone reflection coefficients correlating with curve severity in children with adolescent idiopathic scoliosis

Adolescent idiopathic scoliosis (AIS) is a three-dimensional curvature of spine. Children with AIS and low bone quality have higher chance to get curve progression leading to bigger spinal curvature. In addition, bone quality affects acoustic impedance of bone, thus influencing the reflection coefficient of ultrasound signal from the soft tissue–bone interface. This study aimed to estimate the bone quality of AIS patients based on the reflection coefficients to determine the correlation of the bone quality with curve severity. A simple bone model was used to develop an equation to calculate the reflection coefficient value. Experiments were conducted on five different phantoms. Acrylic was used to design a vertebral shape to study the effect of surface roughness and inclination, including: smooth flat surface (SFS), smooth curved surface (SCS), rough curved surface (RCS), and the rough curved inclined surface (RCIS). A clinical study with 37 AIS patients were recruited. The estimated reflection coefficient values of plate phantoms agreed well with the predicted values and the maximum error was 6.7%. The reflection coefficients measured from the acrylic-water interface for the SFS, SCS, RCS, RCIS (3° and 5°) were 0.37, 0.33, 0.28, (0.23 and 0.12), respectively. The surface roughness and inclination increased the reflection loss. From the clinical data, the average reflection coefficients for children with AIS were 0.11 and 0.07 for the mild curve group and the moderate curve group, respectively. A moderate linear correlation was found between the reflection coefficients and curve severity (r2 = 0.3). Patients with lower bone quality have observed to have larger spinal curvature.

Mussel-Inspired Adhesive Double-Network Hydrogel for Intraoral Ultrasound Imaging

Periodontal diseases could be diagnosed through intraoral ultrasound imaging with the advantages of simple operation procedures, low cost, and low safety risks. A couplant is normally placed between transducers and tissues for better ultrasound image quality. If applied intraorally, the couplants should possess good stability in water and robust mechanical properties, as well as strong adhesiveness to transducers and tissues. However, commercial couplants, such as Aquaflex (AF) cannot fulfill these requirements. In this work, inspired by the mussel adhesion mechanism, we reported a poly(vinyl alcohol)-polyacrylamide-polydopamine (PVA-PAM-PDA) hydrogel synthesized by incorporating PDA into the PAM-PVA double-network for intraoral ultrasound imaging. The hydrogel maintains good stability in water as well as exceptional mechanical properties and can adhere to different substrates (i.e., metal, glass, and porcine skin) without losing the original adhesion strength after multiple adhesion-strip cycles. Besides, when applied to porcine mandibular incisor imaging, the PVA-PAM-PDA hydrogel possesses good image quality for diagnosis as AF does. This work provides practical insights into the fabrication of multifunctional hydrogel-based interfaces between human tissues and medical devices for disease diagnosis applications.

Development and characterization of a tissue mimicking psyllium husk gelatin phantom for ultrasound and magnetic resonance imaging

Purpose: To develop and characterize a tissue-mimicking phantom that enables the direct comparison of magnetic resonance (MR) and ultrasound (US) imaging techniques useful for monitoring high-intensity focused ultrasound (HIFU) treatments. With no additions, gelatin phantoms produce little if any scattering required for US imaging. This study characterizes the MR and US image characteristics as a function of psyllium husk concentration, which was added to increase US scattering.Methods: Gelatin phantoms were constructed with varying concentrations of psyllium husk. The effects of psyllium husk concentration on US B-mode and MR imaging were evaluated at nine different concentrations. T1, T2, and T2* MR maps were acquired. Acoustic properties (attenuation and speed of sound) were measured at frequencies of 0.6, 1.0, 1.8, and 3.0 MHz using a through-transmission technique. Phantom elastic properties were evaluated for both time and temperature dependence.Results: Ultrasound image echogenicity increased with increasing psyllium husk concentration while quality of gradient-recalled echo MR images decreased with increasing concentration. For all phantoms, the measured speed of sound ranged between 1567-1569 m/s and the attenuation ranged between 0.42-0.44 dB/(cm·MHz). Measured T1 ranged from 974-1051 ms. The T2 and T2* values ranged from 97-108 ms and 48-88 ms, respectively, with both showing a decreasing trend with increased psyllium husk concentration. Phantom stiffness, measured using US shear-wave speed measurements, increased with age and decreased with increasing temperature.Conclusions: The presented dual-use tissue-mimicking phantom is easy to manufacture and can be used to compare and evaluate US-guided and MR-guided HIFU imaging protocols.

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.

China's First Demonstration of Cobalt-rich Manganese Crust Thickness Measurement in the Western Pacific with a Parametric Acoustic Probe

Cobalt-rich manganese crusts (CRCs) are important as a potential mineral source that could occur throughout the Pacific on seamounts, ridges, and plateaus. We built a prototype parametric acoustic probe to complete the task of in-situ thickness measurements to estimate the volumetric distribution of deep-sea mineral. The prototype is designed with dual-channels for receiving the primary and secondary signal, which lays a foundation for improving the thickness extraction algorithm. Considering that the signal quality is degraded by the system interference and ambient noise, some improvements to the algorithm are proposed by including the wavelet-based envelope extraction method and the adaptive estimation strategy based on the dual-channel information. Additionally, wavelet regression is applied to reduce the measuring noise assuming that the CRCs have local thickness invariability. The algorithm is suitable for the CRCs with the structure of the multilayers at the top surface and one single layer at the bottom surface. A laboratory experiment is performed to validate the effectiveness of the algorithm. The experiments carried out on the China Ocean 51th voyage in the Western Pacific Ocean on Aug 30, 2018, are described and the data obtained by using the sit-on-bottom stationary measurement are processed to validate the design of the prototype.

A tissue preparation to characterize uterine fibroid tissue properties for thermal therapies

PURPOSE

Treating uterine fibroids with less invasive therapies such as magnetic resonance-guided focused ultrasound (MRgFUS) is an attractive alternative to surgery. Treatment planning can improve MRgFUS procedures and reduce treatment times, but the tissue properties that currently inform treatment planning tools are not adequate. This study aims to develop an ex vivo uterine fibroid model that can emulate the in vivo environment allowing for characterization of the uterus and fibroid MR, acoustic, and thermal tissue properties while maintaining viability for the necessary postsurgical histopathological assessments.

METHODS

Women undergoing a hysterectomy due to fibroid-related symptoms were invited to undergo a preoperative pelvic MRI and to permit postoperative testing of their uterine specimen. Patients that declined or could not be scheduled for a pre-operative MRI were still able to allow post-operative testing of their excised tissue. Following surgical removal of the uterus, nonmorcellated tissues were reperfused with a Krebs-Henseleit buffer solution. An MR-compatible perfusion system was designed to maintain tissue viability inside the MR suite during scanning. MR imaging protocols utilized preoperatively were repeated on whole sample, reperfused ex vivo uterus specimens. Thermal properties including thermal diffusivity and thermal conductivity of the uterus and fibroids were determined using an invasive needle sensor device in 50% of the specimens. Acoustic property measurements (density, speed of sound and attenuation) were obtained for approximately 20% of the tissue samples using both through-transmission and radiation force balance techniques. Differences between fibroid and uterus and in vivo and ex vivo measurements were evaluated with a two-tailed Student t test.

RESULTS

Fourteen patients participated in the study and measurements were obtained from 22 unique fibroids. Of the 16 fibroids available for preoperative MRI testing, 69% demonstrated classic hypo-intensity relative to the myometrium, with the remainder presenting with iso- (25%) or hyper-intensity (6%). While thermal diffusivity was not significantly different between fibroid and myometrium tissues (0.217 ± 0.047 and 0.204 ± 0.039 mm² /s, respectively), the acoustic attenuation in fibroid tissue was significantly higher than myometrium (0.092 ± 0.021 and 0.052 ± 0.023 Np/cm/MHz, respectively). When comparing in vivo with ex vivo MRI T1 and T2 measurements in fibroids and myometrium tissue, the only difference was found in the fibroid T2 property (P < 0.05). Finally, the developed perfusion protocol successfully maintained tissue viability in ex vivo tissues as evaluated through histological analysis.

CONCLUSIONS

This study developed an MR-compatible extracorporeal perfusion technique that effectively maintains tissue viability, allowing for the direct measurement of patient-specific MR, thermal, and acoustic property values for both fibroid and myometrium tissues. These measured tissue property values will enable further development and validation of treatment planning models that can be utilized during MRgFUS uterine fibroid treatments.

Experimental assessment of phase aberration correction for breast MRgFUS therapy

PURPOSE

This study validates that phase aberrations in breast magnetic resonance-guided focussed ultrasound (MRgFUS) therapies can be corrected in a clinically relevant time frame to generate more intense, smaller and more spatially accurate foci.

MATERIALS AND METHODS

Hybrid angular spectrum (HAS) ultrasound calculations in an magnetic resonance imaging (MRI)-based tissue model, were used to compute phase aberration corrections for improved experimental MRgFUS heating in four heterogeneous breast-mimicking phantoms (n = 18 total locations). Magnetic resonance(MR) temperature imaging was used to evaluate the maximum temperature rise, focus volume and focus accuracy for uncorrected and phase aberration-corrected sonications. Thermal simulations assessed the effectiveness of the phase aberration correction implementation.

RESULTS

In 13 of 18 locations, the maximum temperature rise increased by an average of 30%, focus volume was reduced by 40% and focus accuracy improved from 4.6 to 3.6 mm. Mixed results were observed in five of the 18 locations, with focus accuracy improving from 6.1 to 2.5 mm and the maximum temperature rise decreasing by 8% and focus volume increasing by 10%. Overall, the study demonstrated significant improvements (p < 0.005) in maximum temperature rise, focus volume and focus accuracy. Simulations predicted greater improvements than observed experimentally, suggesting potential for improvement in implementing the technique. The complete phase aberration correction procedure, including model generation, segmentation and phase aberration computations, required less than 45 min per sonication location.

CONCLUSION

The significant improvements demonstrated in this study i.e., focus intensity, size and accuracy from phase aberration correction have the potential to improve the efficacy, time-efficiency and safety of breast MRgFUS therapies.

Relationships of the group velocity of the time-reversed Lamb wave with bone properties in cortical bone in vitro

The present study aims to investigate the feasibility of using the time-reversed Lamb wave as a new method for noninvasive characterization of long cortical bones. The group velocity of the time-reversed Lamb wave launched by using the modified time reversal method was measured in 15 bovine tibiae, and their correlations with the bone properties of the tibia were examined. The group velocity of the time-reversed Lamb wave showed significant positive correlations with the bone properties (r=0.55-0.81). The best univariate predictor of the group velocity of the time-reversed Lamb wave was the cortical thickness, yielding an adjusted squared correlation coefficient (r²) of 0.64. These results imply that the group velocity of the time-reversed Lamb wave, in addition to the velocities of the first arriving signal and the slow guided wave, could potentially be used as a discriminator for osteoporosis.

Development and validation of a MRgHIFU non-invasive tissue acoustic property estimation technique

MR-guided high-intensity focussed ultrasound (MRgHIFU) non-invasive ablative surgeries have advanced into clinical trials for treating many pathologies and cancers. A remaining challenge of these surgeries is accurately planning and monitoring tissue heating in the face of patient-specific and dynamic acoustic properties of tissues. Currently, non-invasive measurements of acoustic properties have not been implemented in MRgHIFU treatment planning and monitoring procedures. This methods-driven study presents a technique using MR temperature imaging (MRTI) during low-temperature HIFU sonications to non-invasively estimate sample-specific acoustic absorption and speed of sound values in tissue-mimicking phantoms. Using measured thermal properties, specific absorption rate (SAR) patterns are calculated from the MRTI data and compared to simulated SAR patterns iteratively generated via the Hybrid Angular Spectrum (HAS) method. Once the error between the simulated and measured patterns is minimised, the estimated acoustic property values are compared to the true phantom values obtained via an independent technique. The estimated values are then used to simulate temperature profiles in the phantoms, and compared to experimental temperature profiles. This study demonstrates that trends in acoustic absorption and speed of sound can be non-invasively estimated with average errors of 21% and 1%, respectively. Additionally, temperature predictions using the estimated properties on average match within 1.2 °C of the experimental peak temperature rises in the phantoms. The positive results achieved in tissue-mimicking phantoms presented in this study indicate that this technique may be extended to in vivo applications, improving HIFU sonication temperature rise predictions and treatment assessment.

Micro-scale finite element modeling of ultrasound propagation in aluminum trabecular bone-mimicking phantoms: A comparison between numerical simulation and experimental results

The present study investigated the accuracy of micro-scale finite element modeling for simulating broadband ultrasound propagation in water-saturated trabecular bone-mimicking phantoms. To this end, five commercially manufactured aluminum foam samples as trabecular bone-mimicking phantoms were utilized for ultrasonic immersion through-transmission experiments. Based on micro-computed tomography images of the same physical samples, three-dimensional high-resolution computational samples were generated to be implemented in the micro-scale finite element models. The finite element models employed the standard Galerkin finite element method (FEM) in time domain to simulate the ultrasonic experiments. The numerical simulations did not include energy dissipative mechanisms of ultrasonic attenuation; however, they expectedly simulated reflection, refraction, scattering, and wave mode conversion. The accuracy of the finite element simulations were evaluated by comparing the simulated ultrasonic attenuation and velocity with the experimental data. The maximum and the average relative errors between the experimental and simulated attenuation coefficients in the frequency range of 0.6-1.4 MHz were 17% and 6% respectively. Moreover, the simulations closely predicted the time-of-flight based velocities and the phase velocities of ultrasound with maximum relative errors of 20 m/s and 11 m/s respectively. The results of this study strongly suggest that micro-scale finite element modeling can effectively simulate broadband ultrasound propagation in water-saturated trabecular bone-mimicking structures.

Ultrasound measurements of segmental temperature distribution in solids: Method and its high-temperature validation

A novel approach that uses noninvasive ultrasound to measure the temperature distribution in solid materials is described and validated in high-temperature laboratory experiments. The approach utilizes an ultrasound propagation path with naturally occurring or purposefully introduced echogenic features that partially redirect the energy of an ultrasound excitation pulse back to the transducer, resulting in a train of echoes. Their time of flight depends on the velocity of ultrasound propagation, which changes with temperature distribution in different segments of the propagation path. We reconstruct segmental temperature distributions under different parameterizations. Several parameterizations are discussed, including piecewise constant and piecewise linear, and the parametrization that requires that the estimated temperature profile satisfies an appropriate heat conduction model. The experimental validation of the proposed approach with an alumina sample shows that even with simple parameterizations, the temperature profile is correctly captured with an accuracy that may be comparable to that of the traditional pointwise sensors. The advantages of the approach are discussed, including its suitability for real time and non-destructive temperature measurements in extreme environments and locations inaccessible to the traditional insertion sensors.

Phase aberration simulation study of MRgFUS breast treatments

PURPOSE

This simulation study evaluates the effects of phase aberration in breast MR-guided focused ultrasound (MRgFUS) ablation treatments performed with a phased-array transducer positioned laterally to the breast. A quantification of these effects in terms of thermal dose delivery and the potential benefits of phase correction is demonstrated in four heterogeneous breast numerical models.

METHODS

To evaluate the effects of varying breast tissue properties on the quality of the focus, four female volunteers with confirmed benign fibroadenomas were imaged using 3T MRI. These images were segmented into numerical models with six tissue types, with each tissue type assigned standard acoustic properties from the literature. Simulations for a single-plane 16-point raster-scan treatment trajectory centered in a fibroadenoma in each modeled breast were performed for a breast-specific MRgFUS system. At each of the 16 points, pressure patterns both with and without applying a phase correction technique were determined with the hybrid-angular spectrum method. Corrected phase patterns were obtained using a simulation-based phase aberration correction technique to adjust each element's transmit phase to obtain maximized constructive interference at the desired focus. Thermal simulations were performed for both the corrected and uncorrected pressure patterns using a finite-difference implementation of the Pennes bioheat equation. The effect of phase correction was evaluated through comparison of thermal dose accumulation both within and outside a defined treatment volume. Treatment results using corrected and uncorrected phase aberration simulations were compared by evaluating the power required to achieve a 20 °C temperature rise at the first treatment location. The extent of the volumes that received a minimum thermal dose of 240 CEM at 43 °C inside the intended treatment volume as well as the volume in the remaining breast tissues was also evaluated in the form of a dose volume ratio (DVR), a DVR percent change between corrected and uncorrected phases, and an additional metric that measured phase spread.

RESULTS

With phase aberration correction applied, there was an improvement in the focus for all breast anatomies as quantified by a reduction in power required (13%-102%) to reach 20 °C when compared to uncorrected simulations. Also, the DVR percent change increased by 5%-77% in seven out of eight cases, indicating an improvement to the treatment as measured by a reduction in thermal dose deposited to the nontreatment tissues. Breast compositions with a higher degree of heterogeneity along the ultrasound beam path showed greater reductions in thermal dose delivered outside of the treatment volume with correction applied than beam trajectories that propagated through more homogeneous breast compositions. An increasing linear trend was observed between the DVR percent change and the phase-spread metric (R(2) = 0.68).

CONCLUSIONS

These results indicate that performing phase aberration correction for breast MRgFUS treatments is beneficial for the small-aperture transducer (14.4 × 9.8 cm) evaluated in this work. While all breast anatomies could benefit from phase aberration correction, greater benefits are observed in more heterogeneous anatomies.

A simulation technique for 3D MR-guided acoustic radiation force imaging

PURPOSE

In magnetic resonance-guided focused ultrasound (MRgFUS) therapies, the in situ characterization of the focal spot location and quality is critical. MR acoustic radiation force imaging (MR-ARFI) is a technique that measures the tissue displacement caused by the radiation force exerted by the ultrasound beam. This work presents a new technique to model the displacements caused by the radiation force of an ultrasound beam in a homogeneous tissue model.

METHODS

When a steady-state point-source force acts internally in an infinite homogeneous medium, the displacement of the material in all directions is given by the Somigliana elastostatic tensor. The radiation force field, which is caused by absorption and reflection of the incident ultrasound intensity pattern, will be spatially distributed, and the tensor formulation takes the form of a convolution of a 3D Green's function with the force field. The dynamic accumulation of MR phase during the ultrasound pulse can be theoretically accounted for through a time-of-arrival weighting of the Green's function. This theoretical model was evaluated experimentally in gelatin phantoms of varied stiffness (125-, 175-, and 250-bloom). The acoustic and mechanical properties of the phantoms used as parameters of the model were measured using independent techniques. Displacements at focal depths of 30- and 45-mm in the phantoms were measured by a 3D spin echo MR-ARFI segmented-EPI sequence.

RESULTS

The simulated displacements agreed with the MR-ARFI measured displacements for all bloom values and focal depths with a normalized RMS difference of 0.055 (range 0.028-0.12). The displacement magnitude decreased and the displacement pattern broadened with increased bloom value for both focal depths, as predicted by the theory.

CONCLUSIONS

A new technique that models the displacements caused by the radiation force of an ultrasound beam in a homogeneous tissue model theory has been rigorously validated through comparison with experimentally obtained 3D displacement data in homogeneous gelatin phantoms using a 3D MR-ARFI sequence. The agreement of the experimentally measured and simulated results demonstrates the potential to use MR-ARFI displacement data in MRgFUS therapies.

Characterization and evaluation of tissue-mimicking gelatin phantoms for use with MRgFUS

Background

A tissue-mimicking phantom that accurately represents human-tissue properties is important for safety testing and for validating new imaging techniques. To achieve a variety of desired human-tissue properties, we have fabricated and tested several variations of gelatin phantoms. These phantoms are simple to manufacture and have properties in the same order of magnitude as those of soft tissues. This is important for quality-assurance verification as well as validation of magnetic resonance-guided focused ultrasound (MRgFUS) treatment techniques.

Methods

The phantoms presented in this work were constructed from gelatin powders with three different bloom values (125, 175, and 250), each one allowing for a different mechanical stiffness of the phantom. Evaporated milk was used to replace half of the water in the recipe for the gelatin phantoms in order to achieve attenuation and speed of sound values in soft tissue ranges. These acoustic properties, along with MR (T₁ and T₂*), mechanical (density and Young’s modulus), and thermal properties (thermal diffusivity and specific heat capacity), were obtained through independent measurements for all three bloom types to characterize the gelatin phantoms. Thermal repeatability of the phantoms was also assessed using MRgFUS and MR thermometry.

Results

All the measured values fell within the literature-reported ranges of soft tissues. In heating tests using low-power (6.6 W) sonications, interleaved with high-power (up to 22.0 W) sonications, each of the three different bloom phantoms demonstrated repeatable temperature increases (10.4 ± 0.3 °C for 125-bloom, 10.2 ± 0.3 °C for 175-bloom, and 10.8 ± 0.2 °C for 250-bloom for all 6.6-W sonications) for heating durations of 18.1 s.

Conclusion

These evaporated milk-modified gelatin phantoms should serve as reliable, general soft tissue-mimicking MRgFUS phantoms.

Propagation of time-reversed Lamb waves in bovine cortical bone in vitro

The present study aims to investigate the propagation of time-reversed Lamb waves in bovine cortical bone in vitro. The time-reversed Lamb waves were successfully launched at 200 kHz in 18 bovine tibiae through a time reversal process of Lamb waves. The group velocities of the time-reversed Lamb waves in the bovine tibiae were measured using the axial transmission technique. They showed a significant correlation with the cortical thickness and tended to follow the theoretical group velocity of the lowest order antisymmetrical Lamb wave fairly well, consistent with the behavior of the slow guided wave in long cortical bones.

Multichannel filtering and reconstruction of ultrasonic guided wave fields using time intercept-slowness transform

Multichannel ultrasonic axial-transmission data are multimodal by nature. As guided waves are commonly used in nondestructive material testing, wave field filtering becomes important because the analysis is usually limited to a few lower-order modes and requires their extraction. An application of the Radon transform to enhance signal-to-noise ratio and separate wave fields in ultrasonic records is presented. The method considers guided wave fields as superpositions of plane waves defined by ray parameters (p) and time intercepts (τ) and stacks the amplitudes along linear trajectories, mapping time-offset (t - x) data to a τ - p or Radon panel. The transform is implemented using a least-squares strategy with Cauchy-norm regularization that serves to enhance the focusing power. The method was verified using simulated data and applied to an uneven spatially sampled bovine-bone-plate data set. The results demonstrate the Radon panels show isolated amplitude clusters and the Cauchy-norm constraint provides a more focused Radon image than the damped least-squares regularization. Wave field separation can be achieved by selectively windowing the τ - p signals and inverse transformation, which is illustrated by the successful extraction of the A0 mode in bone plate. In addition, the method effectively attenuates noise, enhances the coherency of the guided wave modes, and reconstructs the missing records. The proposed transform presents a powerful signal-enhancement tool to process guided waves for further analysis and inversion.

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.

Measurements of ultrasonic phase velocities and attenuation of slow waves in cellular aluminum foams as cancellous bone-mimicking phantoms

The water-saturated aluminum foams with an open network of interconnected ligaments were investigated by ultrasonic transmission technique for the suitability as cancellous bone-mimicking phantoms. The phase velocities and attenuation of nine samples covering three pores per inch (5, 10, and 20 PPI) and three aluminum volume fractions (5, 8, and 12% AVF) were measured over a frequency range of 0.7-1.3 MHz. The ligament thickness and pore sizes of the phantoms and low-density human cancellous bones are similar. A strong slow wave and a weak fast wave are observed for all samples while the latter is not visible without significant amplification (100x). This study reports the characteristics of slow wave, whose speeds are less than the sound speed of the saturating water and decrease mildly with AVF and PPI with an average 1469 m/s. Seven out of nine samples show positive dispersion and the rest show minor negative dispersion. Attenuation increases with AVF, PPI, and frequency except for the 20 PPI samples, which exhibit non-increasing attenuation level with fluctuations due to scattering. The phase velocities agree with Biot's porous medium theory. The RMSE is 16.0 m/s (1%) at n = 1.5. Below and above this value, the RMSE decreases mildly and rises sharply, respectively.

Measurement of tortuosity in aluminum foams using airborne ultrasound

The slow compressional wave in air-saturated aluminum foams was studied by means of ultrasonic transverse transmission method over a frequency range from 0.2 MHz to 0.8 MHz. The samples investigated have three different cell sizes or pores per inch (5, 10 and 20 ppi) and each size has three aluminum volume fractions (5%, 8% and 12% AVF). Phase velocities show minor dispersion at low frequencies but remain constant after 0.7 MHz. Pulse broadening and amplitude attenuation are obvious and increase with increasing ppi. Attenuation increases considerably with AVF for 20 ppi foams. Tortuosity ranges from 1.003 to 1.032 and increases with AVF and ppi. However, the increase of tortuosity with AVF is very small for 10 and 20 ppi samples.

Broadband ultrasound attenuation measurement of long bone using peak frequency of the echoes

The traditional peak frequency formulation requires the knowledge of the signal wavelet. We improved and applied the method to estimate attenuation for homogeneous silicon rubber and bovine cortical bone without recourse to the wavelet assumption. The estimated values for rubber and bone samples are 6.59 +/- 0.28 dB/MHz/cm and 4.59 +/- 1.09 dB/MHz/cm, respectively, as compared with 6.33 +/- 0.19 dB/ MHz/cm and 5.00 +/- 1.10 dB/MHz/cm, respectively, by the spectral ratio method.

A novel approach to acoustic liquid density measurements using a buffer rod based measuring cell

A new method for measuring the pressure reflection coefficient in a buffer rod configuration is presented, together with experimental results for acoustic measurements of the liquid density, based on the measurement of the liquid's acoustic impedance. The method consists of using 2 buffers enclosing the liquid in a symmetrical arrangement with a transducer fixed to each buffer. One of the transducers is used in a pulse-echo mode while the other transducer operates as a receiver. The echo amplitudes leading to the pressure reflection coefficient as found by this method possess advantages such as reduced attenuation due to a shorter liquid transmission path and reduced interference, as compared with the ABC method. Measurements with distilled water and with special density calibration oil qualities have been performed using both the new method and the ABC method and are shown for the new method to give a density span within +/- 0.15% of the reference values. A comparison of the measured densities based on both a time-domain and a l(2)-norm frequency domain integration signal processing approach is given, along with a recommendation as to how the signal processing should be performed.

Model-based estimation of quantitative ultrasound variables at the proximal femur

To improve the prediction of the osteoporotic fracture risk at the proximal femur we are developing a scanner for quantitative ultrasound (QUS) measurements at this site. Due to multipath transmission in this complex shaped bone, conventional signal processing techniques developed for QUS measurements at peripheral sites frequently fail. Therefore, we propose a model-based estimation of the QUS variables and analyze the performance of the new algorithm. Applying the proposed method to QUS scans of excised proximal femurs increased the fraction of evaluable signals from approx. 60% (using conventional algorithms) to 97%. The correlation of the standard QUS variables broadband ultrasound attenuation (BUA) and speed of sound (SOS) with the established variable bone mineral density (BMD) reported in previous studies is maintained (BUA/BMD: r(2) = 0.69; SOS/BMD: r(2) = 0.71; SOS+BUA/BMD: r(2) = 0.88). Additionally, different wave types could be clearly detected and characterized in the trochanteric region. The ability to separate superimposed signals with this approach opens up further diagnostic potential for evaluating waves of different sound paths and wave types through bone tissue.

Spectral ratio method to estimate broadband ultrasound attenuation of cortical bonesin vitrousing multiple reflections

Broadband ultrasound attenuation (BUA) is commonly measured by the spectral ratio method. Conventionally BUA is measured in transverse transmission mode where ultrasound signal is recorded with and without the sample. The spectral ratio method was extended to estimate nBUA (BUA normalized by thickness) in axial transmission mode using spectral amplitudes of the primary reflection and multiple reflection, which echoes more than once between the material interfaces within a layer. We performed three experiments. First, reflections were numerically simulated to verify the accuracy of the method. We then applied the method to estimate attenuation of silicon rubber and the cortex of a bovine femur. The center frequency of the transducers is 2.25 MHz. We obtained 93% accuracy for a simulated data set with 10% random noise after bandpass filtering. For the silicon rubber, 15 measurements were collected and the mean attenuation was 6.33 +/- 0.19 dB MHz(-1) cm(-1). For the bovine bone, eight measurements were performed in the middle portion of the femur. The mean attenuation was 4.91 +/- 0.65 dB MHz(-1) cm(-1) and compared well with those reported in the literature. The results demonstrate that the proposed method has the potential to provide a quick, reliable and robust cortical attenuation assessment in vivo.