The first days of diagnostic ultrasound

The Reflections series takes a look back on historical articles from The Journal of the Acoustical Society of America that have had a significant impact on the science and practice of acoustics.

A Qualitative Analysis of Clinician Perspectives of Ultrasound Biofeedback for Speech Sound Disorders

PURPOSE

Ultrasound biofeedback therapy (UBT) is a relatively new type of technology-assisted speech-language therapy and has shown promise in remediating speech sound disorders. However, there is a current lack of understanding of the barriers and benefits that may influence the usage behavior and clinical decision making for the implementation of UBT from a clinician perspective. In this qualitative study, we explore the perspectives of speech-language pathologists (SLPs) who have used ultrasound biofeedback in programs of speech sound therapy using the unified theory of acceptance and use of technology (UTAUT) model.

METHOD

Seven SLPs who had clinical experience treating speech sound disorders with UBT participated. Semistructured in-depth interviews were conducted and video-recorded. Two coders coded and categorized the transcribed data, with consensus established with a third coder. Using thematic analysis, the data were exploratorily grouped into themes along components of the UTAUT model.

RESULTS

The highest number of codes was sorted into the "effort expectancy" theme, followed by "performance expectancy," "social influence," and "facilitating conditions" themes of the UTAUT model. Clinicians identified multiple perceived barriers and benefits to the use of ultrasound technology. The top identified barrier was limited accessibility, and the top benefit was the ability to visualize a client's articulatory response to cues on a display.

CONCLUSIONS

Clinicians prioritized "effort expectancy" and "performance expectancy" when reflecting on the use of ultrasound biofeedback for speech sound disorders. Clinicians spoke favorably about using UBT for speech sound disorder treatment but acknowledged institutional barriers and limitations at organizational and social levels.

Classification of accurate and misarticulated /r/ for ultrasound biofeedback using tongue part displacement trajectories

Ultrasound biofeedback therapy (UBT), which incorporates real-time imaging of tongue articulation, has demonstrated generally positive speech remediation outcomes for individuals with residual speech sound disorder (RSSD). However, UBT requires high attentional demands and may therefore benefit from a simplified display of articulation targets that are easily interpretable and can be compared to real-time articulation. Identifying such targets requires automatic quantification and analysis of movement features relevant to accurate speech production. Our image-analysis program TonguePART automatically quantifies tongue movement as tongue part displacement trajectories from midsagittal ultrasound videos of the tongue, with real-time capability. The present study uses such displacement trajectories to compare accurate and misarticulated American-English rhotic /ɑr/ productions from 40 children, with degree of accuracy determined by auditory perceptual ratings. To identify relevant features of accurate articulation, support vector machine (SVM) classifiers were trained and evaluated on several candidate data representations. Classification accuracy was up to 85%, indicating that quantification of tongue part displacement trajectories captured tongue articulation characteristics that distinguish accurate from misarticulated production of /ɑr/. Regression models for perceptual ratings were also compared. The simplest data representation that retained high predictive ability, demonstrated by high classification accuracy and strong correlation between observed and predicted ratings, was displacements at the midpoint of /r/ relative to /ɑ/ for the tongue dorsum and blade. This indicates that movements of the dorsum and blade are especially relevant to accurate production of /r/, suggesting that a predictive parameter and biofeedback target based on this data representation may be usable for simplified UBT.

Reconstruction of thermoacoustic emission sources induced by proton irradiation using numerical time reversal

Objective.Mapping of dose delivery in proton beam therapy can potentially be performed by analyzing thermoacoustic emissions measured by ultrasound arrays. Here, a method is derived and demonstrated for spatial mapping of thermoacoustic sources using numerical time reversal, simulating re-transmission of measured emissions into the medium.Approach.Spatial distributions of thermoacoustic emission sources are shown to be approximated by the analytic-signal form of the time-reversed acoustic field, evaluated at the time of the initial proton pulse. Given calibration of the array sensitivity and knowledge of tissue properties, this approach approximately reconstructs the acoustic source amplitude, equal to the product of the time derivative of the radiation dose rate, mass density, and Grüneisen parameter. This approach was implemented using two models for acoustic fields of the array elements, one modeling elements as line sources and the other as rectangular radiators. Thermoacoustic source reconstructions employed previously reported measurements of emissions from proton energy deposition in tissue-mimicking phantoms. For a phantom incorporating a bone layer, reconstructions accounted for the higher sound speed in bone. Dependence of reconstruction quality on array aperture size and signal-to-noise ratio was consistent with previous acoustic simulation studies.Main results.Thermoacoustic source distributions were successfully reconstructed from acoustic emissions measured by a linear ultrasound array. Spatial resolution of reconstructions was significantly improved in the azimuthal (array) direction by incorporation of array element diffraction. Source localization agreed well with Monte Carlo simulations of energy deposition, and was improved by incorporating effects of inhomogeneous sound speed.Significance.The presented numerical time reversal approach reconstructs thermoacoustic sources from proton beam radiation, based on straightforward processing of acoustic emissions measured by ultrasound arrays. This approach may be useful for ranging and dosimetry of clinical proton beams, if acoustic emissions of sufficient amplitude and bandwidth can be generated by therapeutic proton sources.

Three-dimensional echo decorrelation monitoring of radiofrequency ablation in ex vivo bovine liver

Three-dimensional (3D) echo decorrelation imaging was investigated for monitoring radiofrequency ablation (RFA) in ex vivo bovine liver. RFA experiments (N = 14) were imaged by 3D ultrasound using a matrix array, with in-phase and quadrature complex echo volumes acquired about every 11 s. Tissue specimens were then frozen at -80 °C, sectioned, and semi-automatically segmented. Receiver operating characteristic (ROC) curves were constructed for assessing ablation prediction performance of 3D echo decorrelation with three potential normalization approaches, as well as 3D integrated backscatter (IBS). ROC analysis indicated that 3D echo decorrelation imaging is potentially a good predictor of local RFA, with the best prediction performance observed for globally normalized decorrelation. Tissue temperatures, recorded by four thermocouples integrated into the RFA probe, showed good correspondence with spatially averaged decorrelation and statistically significant but weak correlation with measured echo decorrelation at the same spatial locations. In tests predicting ablation zones using a weighted K-means clustering approach, echo decorrelation performed better than IBS, with smaller root mean square volume errors and higher Dice coefficients relative to measured ablation zones. These results suggest that 3D echo decorrelation and IBS imaging are capable of real-time monitoring of thermal ablation, with potential application to clinical treatment of liver tumors.

Reconstruction of shear wave speed in tissue-mimicking phantoms from aliased pulse-echo imaging of high-frequency wavefields

Quantitative elasticity estimation in medical and industrial applications may benefit from advancements in reconstruction of shear wave speed with enhanced resolution. Here, shear wave speed is reconstructed from pulse-echo ultrasound imaging of elastic waves induced by high-frequency (>400 Hz), time-harmonic mechanical excitation. Particle displacement in shear wavefields is mapped from measured interframe phase differences with compensation for timing of multiple scan lines, then processed by spatial Fourier analysis to estimate the predominant wave speed and analyzed by algebraic wavefield inversion to reconstruct wave speed maps. Reconstructions of shear wave speed from simulated wavefields illustrate the accuracy and spatial resolution available with both methods, as functions of signal-to-noise ratio and sizes of windows used for Fourier analysis or wavefield smoothing. The methods are applied to shear wavefields with frequencies up to six times the Nyquist rate, thus extending the frequency range measurable by a given imaging system. Wave speed measurements in tissue-mimicking phantoms are compared with supersonic shear imaging and mechanical tensile testing, demonstrating feasibility of the wavefield measurement and wave speed reconstruction methods employed.

Effect of Overpressure on Acoustic Emissions and Treated Tissue Histology in ex Vivo Bulk Ultrasound Ablation

Bulk ultrasound ablation is a thermal therapy approach in which tissue is heated by unfocused or weakly focused sonication (average intensities on the order of 100 W/cm²) to achieve coagulative necrosis within a few minutes exposure time. Assessing the role of bubble activity, including acoustic cavitation and tissue vaporization, in bulk ultrasound ablation may help in making bulk ultrasound ablation safer and more effective for clinical applications. Here, two series of ex vivo ablation trials were conducted to investigate the role of bubble activity and tissue vaporization in bulk ultrasound ablation. Fresh bovine liver tissue was ablated with unfocused, continuous-wave ultrasound using ultrasound image-ablate arrays sonicating at 31 W/cm² (0.9 MPa amplitude) for either 20 min at a frequency of 3.1 MHz or 10 min at 4.8 MHz. Tissue specimens were maintained at a static overpressure of either 0.52 or 1.2 MPa to suppress bubble activity and tissue vaporization or at atmospheric pressure for control groups. A passive cavitation detector was used to record subharmonic (1.55 or 2.4 MHz), broadband (1.2-1.5 MHz) and low-frequency (5-20 kHz) acoustic emissions. Treated tissue was stained with 2% triphenyl tetrazolium chloride to evaluate thermal lesion dimensions. Subharmonic emissions were significantly reduced in overpressure groups compared with control groups. Correlations observed between acoustic emissions and lesion dimensions were significant and positive for the 3.1-MHz series, but significant and negative for the 4.8-MHz series. The results indicate that for bulk ultrasound ablation, where both acoustic cavitation and tissue vaporization are possible, bubble activity can enhance ablation in the absence of tissue vaporization, but can reduce thermal lesion dimensions in the presence of vaporization.

Tongue Part Movement Trajectories for /r/ Using Ultrasound

Purpose

Because it shows the movement of different parts of the tongue in real time, ultrasound biofeedback therapy is a promising technology for speech research and remediation. One limitation is the difficulty of interpreting real-time ultrasound images of tongue motion. Our image processing system, TonguePART, tracks the tongue surface and allows for the acquisition of quantitative tongue part trajectories.

Method

TonguePART automatically identifies the tongue contour based on ultrasound image brightness and tracks motion of the tongue root, dorsum, and blade in real time. We present tongue part trajectory data from 2 children with residual sound errors on /r/ and 2 children with typical speech, focusing on /r/ (International Phonetic Alphabet ɹ) in the phonetic context /ɑr/. We compared the tongue trajectories to magnetic resonance images of sustained vowel /ɑ/ and /r/.

Results

Measured trajectories show larger overall displacement and greater differentiation of tongue part movements for children with typical speech during the production of /ɑr/, compared to children with residual speech sound disorders.

Conclusion

TonguePART is a fast, reliable method of tracking articulatory movement of tongue parts for syllables such as /ɑr/. It is extensible to other sounds and phonetic contexts. By tracking tongue parts, clinical researchers can investigate lingual coordination. TonguePART is suitable for real-time data collection and biofeedback. Ultrasound biofeedback therapy users may make more progress using simplified biofeedback of tongue movement.

Numerical analysis of three-dimensional echo decorrelation imaging

A numerical model for three-dimensional echo decorrelation imaging, a pulse-echo ultrasound method applicable to thermal ablation monitoring, is presented. Beam patterns for steered transmit and receive array apertures are combined with a three-dimensional numerical tissue model to yield beamformed scan lines in a pyramidal configuration, volumetric B-mode images, and spatial maps of normalized decorrelation between sequential image volumes. Simulated three-dimensional echo decorrelation images of random media are analyzed as estimators of local tissue reflectivity decoherence, mimicking thermal ablation effects. The estimation error is analyzed as a function of correlation window size, scan line density, and ensemble averaging of decorrelation maps.

In vivo ultrasound thermal ablation control using echo decorrelation imaging in rabbit liver and VX2 tumor

The utility of echo decorrelation imaging feedback for real-time control of in vivo ultrasound thermal ablation was assessed in rabbit liver with VX2 tumor. High-intensity focused ultrasound (HIFU) and unfocused (bulk) ablation were performed using 5 MHz linear image-ablate arrays. Treatments comprised up to nine lower-power sonications, followed by up to nine higher-power sonications, ceasing when the average cumulative echo decorrelation within a control region of interest exceeded a predefined threshold (- 2.3, log10-scaled echo decorrelation per millisecond, corresponding to 90% specificity for tumor ablation prediction in previous in vivo experiments). This threshold was exceeded in all cases for both HIFU (N = 12) and bulk (N = 8) ablation. Controlled HIFU trials achieved a significantly higher average ablation rate compared to comparable ablation trials without image-based control, reported previously. Both controlled HIFU and bulk ablation trials required significantly less treatment time than these previous uncontrolled trials. Prediction of local liver and VX2 tumor ablation using echo decorrelation was tested using receiver operator characteristic curve analysis, showing prediction capability statistically equivalent to uncontrolled trials. Compared to uncontrolled trials, controlled trials resulted in smaller thermal ablation regions and higher contrast between echo decorrelation in treated vs. untreated regions. These results indicate that control using echo decorrelation imaging may reduce treatment duration and increase treatment reliability for in vivo thermal ablation.

The effect of 220 kHz insonation scheme on rt-PA thrombolytic efficacy in vitro

Ultrasound-enhanced recombinant tissue plasminogen activator (rt-PA) thrombolysis is under development as an adjuvant to ischemic stroke therapy. The goal of this study was to design a pulsed ultrasound (US) exposure scheme that reduced intracranial constructive interference and tissue heating, and maintained thrombolytic efficacy relative to continuous wave (CW) insonation. Three 220 kHz US schemes were evaluated, two pulsed insonation schemes (15 cycles, 68 µs pulse duration, 33% or 62.5% duty cycle) and an intermittent CW insonation scheme (50 s active, 30 s quiescent) over a 30-min treatment period. An in silico study using a finite-difference model of transcranial US propagation was performed to estimate the intracranial acoustic field and temperature rise in the skull for each insonation scheme. In vitro measurements with flow were performed to assess thrombolysis using time-lapse microscopy. Intracranial constructive interference was not reduced with pulsed US using a pulse length of 15 cycles compared to intermittent CW US. The 33.3% duty cycle pulsed US scheme reduced heating in the temporal bone as much as 60% relative to the intermittent CW scheme. All insonation schemes promoted sustained stable cavitation in vitro and augmented thrombolysis compared to rt-PA alone (p   <  0.05). Ultraharmonic (UH) and harmonic cumulative energy over a 30 min treatment period was significantly higher (p   <  0.05) for the intermittent CW US scheme compared to either pulsed US scheme. Despite the difference in cavitation emissions, no difference was observed in the clot lysis between the three US schemes. These findings demonstrate that a 33.3% duty cycle pulsed US scheme with a 15-cycle burst can reduce bone heating and achieve equivalent thrombolytic efficacy as an intermittent CW scheme.

Characterization of cavitation-radiated acoustic power using diffraction correction

A method is developed for compensating absolute pressure measurements made by a calibrated passive cavitation detector (PCD) to estimate the average acoustic power radiated from a region of interest (ROI) defined to encompass all cavitating bubbles. A diffraction correction factor for conversion of PCD-measured pressures to cavitation-radiated acoustic power per unit area or volume is derived as a simple analytic expression, accounting for position- and frequency-dependent PCD sensitivity. This approach can be applied to measurements made by any PCD without precise knowledge of the number, spatial, or temporal distribution of cavitating bubbles. The diffraction correction factor is validated in simulation for a wide range of ROI dimensions and frequencies. The correction factor is also applied to emission measurements obtained during in vitro ultrasound-enhanced sonophoresis experiments, allowing comparison of stable cavitation levels between therapeutic configurations with different source center frequencies. Results incorporating sonication at both 0.41 and 2.0 MHz indicate that increases in skin permeability correlate strongly with the acoustic power of subharmonic emissions radiated per unit skin area.

Optimized Echo Decorrelation Imaging Feedback for Bulk Ultrasound Ablation Control

Feasibility of controlling bulk ultrasound (US) thermal ablation using echo decorrelation imaging was investigated in ex vivo bovine liver. The first of two ablation and control procedures used a sequence of constant-intensity sonication cycles, ceased when the minimum echo decorrelation within a control region of interest (ROI) exceeded a predetermined threshold. The second procedure used a variable-intensity sonication sequence, with spatially averaged decorrelation as the stopping criterion. US exposures and echo decorrelation imaging were performed by a linear image-ablate array. Based on preliminary experiments, control ROIs and thresholds for the minimum-decorrelation and average-decorrelation criteria were specified. Controlled trials for the minimum-decorrelation and average-decorrelation criteria were compared with uncontrolled trials employing 9 or 18 cycles of matching sonication sequences. Lesion dimensions, treatment times, ablation rates, and areas under receiver operating characteristic curves were statistically compared. Successfully controlled trials using both criteria required significantly shorter treatment times than corresponding 18-cycle treatments, with better ablation prediction performance than uncontrolled 9-cycle and 18-cycle treatments. Either control approach resulted in greater ablation rate than corresponding 9-cycle or 18-cycle uncontrolled approaches. A post hoc analysis studied the effect of exchanging control criteria between the two series of controlled experiments. For either group, the average time needed to exceed the alternative decorrelation threshold approximately matched the average duration of successfully controlled experimental trials. These results indicate that either approach, using minimum-decorrelation or average-decorrelation criteria, is feasible for control of bulk US ablation. In addition, use of a variable-intensity sonication sequence for bulk US thermal ablation can result in larger ablated regions compared to constant-intensity sonication sequences.

Real-Time Spatiotemporal Control of High-Intensity Focused Ultrasound Thermal Ablation Using Echo Decorrelation Imaging in ex Vivo Bovine Liver

The ability to control high-intensity focused ultrasound (HIFU) thermal ablation using echo decorrelation imaging feedback was evaluated in ex vivo bovine liver. Sonications were automatically ceased when the minimum cumulative echo decorrelation within the region of interest exceeded an ablation control threshold, determined from preliminary experiments as -2.7 (log-scaled decorrelation per millisecond), corresponding to 90% specificity for local ablation prediction. Controlled HIFU thermal ablation experiments were compared with uncontrolled experiments employing two, five or nine sonication cycles. Means and standard errors of the lesion width, area and depth, as well as receiver operating characteristic curves testing ablation prediction performance, were computed for each group. Controlled trials exhibited significantly smaller average lesion area, width and treatment time than five-cycle or nine-cycle uncontrolled trials and also had significantly greater prediction capability than two-cycle uncontrolled trials. These results suggest echo decorrelation imaging is an effective approach to real-time HIFU ablation control.

Quantitative Frequency-Domain Passive Cavitation Imaging

Passive cavitation detection has been an instrumental technique for measuring cavitation dynamics, elucidating concomitant bioeffects, and guiding ultrasound therapies. Recently, techniques have been developed to create images of cavitation activity to provide investigators with a more complete set of information. These techniques use arrays to record and subsequently beamform received cavitation emissions, rather than processing emissions received on a single-element transducer. In this paper, the methods for performing frequency-domain delay, sum, and integrate passive imaging are outlined. The method can be applied to any passively acquired acoustic scattering or emissions, including cavitation emissions. To compare data across different systems, techniques for normalizing Fourier transformed data and converting the data to the acoustic energy received by the array are described. A discussion of hardware requirements and alternative imaging approaches is additionally outlined. Examples are provided in MATLAB.

Echo Decorrelation Imaging of Rabbit Liver and VX2 Tumor during In Vivo Ultrasound Ablation

In open surgical procedures, image-ablate ultrasound arrays performed thermal ablation and imaging on rabbit liver lobes with implanted VX2 tumor. Treatments included unfocused (bulk ultrasound ablation, N = 10) and focused (high-intensity focused ultrasound ablation, N = 13) exposure conditions. Echo decorrelation and integrated backscatter images were formed from pulse-echo data recorded during rest periods after each therapy pulse. Echo decorrelation images were corrected for artifacts using decorrelation measured prior to ablation. Ablation prediction performance was assessed using receiver operating characteristic curves. Results revealed significantly increased echo decorrelation and integrated backscatter in both ablated liver and ablated tumor relative to unablated tissue, with larger differences observed in liver than in tumor. For receiver operating characteristic curves computed from all ablation exposures, both echo decorrelation and integrated backscatter predicted liver and tumor ablation with statistically significant success, and echo decorrelation was significantly better as a predictor of liver ablation. These results indicate echo decorrelation imaging is a successful predictor of local thermal ablation in both normal liver and tumor tissue, with potential for real-time therapy monitoring.

Dependence of ultrasound echo decorrelation on local tissue temperature during ex vivo radiofrequency ablation

This study investigates echo decorrelation imaging, an ultrasound method for thermal ablation monitoring. The effect of tissue temperature on the mapped echo decorrelation parameter was assessed in radiofrequency ablation experiments performed on ex vivo bovine liver tissue. Echo decorrelation maps were compared with corresponding tissue temperatures simulated using the finite element method. For both echo decorrelation imaging and integrated backscatter imaging, the mapped tissue parameters correlated significantly but weakly with local tissue temperature. Receiver operating characteristic (ROC) curves were used to assess the ability of echo decorrelation and integrated backscatter to predict tissue temperature greater than 40, 60, and 80 °C. Significantly higher area under the ROC curve (AUROC) values were obtained for prediction of tissue temperatures greater than 40, 60, and 80 °C using echo decorrelation imaging (AUROC = 0.871, 0.948 and 0.966) compared to integrated backscatter imaging (AUROC = 0.865, 0.877 and 0.832).

Accuracy of a bistatic scattering substitution technique for calibration of focused receivers

A recent method for calibrating single-element, focused passive cavitation detectors (PCD) compares bistatic scattering measurements by the PCD and a reference hydrophone. Here, effects of scatterer properties and PCD size on frequency-dependent receive calibration accuracy are investigated. Simulated scattering from silica and polystyrene spheres was compared for small hydrophone and spherically focused PCD receivers to assess the achievable calibration accuracy as a function of frequency, scatterer size, and PCD size. Good agreement between measurements was found when the scatterer diameter was sufficiently smaller than the focal beamwidth of the PCD; this relationship was dependent on the scatterer material. For conditions that result in significant disagreement between measurements, the numerical methods described here can be used to correct experimental calibrations.

Optimization of tissue physical parameters for accurate temperature estimation from finite-element simulation of radiofrequency ablation

Computational finite element models are commonly used for the simulation of radiofrequency ablation (RFA) treatments. However, the accuracy of these simulations is limited by the lack of precise knowledge of tissue parameters. In this technical note, an inverse solver based on the unscented Kalman filter (UKF) is proposed to optimize values for specific heat, thermal conductivity, and electrical conductivity resulting in accurately simulated temperature elevations. A total of 15 RFA treatments were performed on ex vivo bovine liver tissue. For each RFA treatment, 15 finite-element simulations were performed using a set of deterministically chosen tissue parameters to estimate the mean and variance of the resulting tissue ablation. The UKF was implemented as an inverse solver to recover the specific heat, thermal conductivity, and electrical conductivity corresponding to the measured area of the ablated tissue region, as determined from gross tissue histology. These tissue parameters were then employed in the finite element model to simulate the position- and time-dependent tissue temperature. Results show good agreement between simulated and measured temperature.

Using passive cavitation images to classify high-intensity focused ultrasound lesions

Passive cavitation imaging provides spatially resolved monitoring of cavitation emissions. However, the diffraction limit of a linear imaging array results in relatively poor range resolution. Poor range resolution has limited prior analyses of the spatial specificity and sensitivity of passive cavitation imaging in predicting thermal lesion formation. In this study, this limitation is overcome by orienting a linear array orthogonal to the high-intensity focused ultrasound propagation direction and performing passive imaging. Fourteen lesions were formed in ex vivo bovine liver samples as a result of 1.1-MHz continuous-wave ultrasound exposure. The lesions were classified as focal, "tadpole" or pre-focal based on their shape and location. Passive cavitation images were beamformed from emissions at the fundamental, harmonic, ultraharmonic and inharmonic frequencies with an established algorithm. Using the area under a receiver operating characteristic curve (AUROC), fundamental, harmonic and ultraharmonic emissions were found to be significant predictors of lesion formation for all lesion types. For both harmonic and ultraharmonic emissions, pre-focal lesions were classified most successfully (AUROC values of 0.87 and 0.88, respectively), followed by tadpole lesions (AUROC values of 0.77 and 0.64, respectively) and focal lesions (AUROC values of 0.65 and 0.60, respectively).

Methods to calibrate the absolute receive sensitivity of single-element, focused transducers

Absolute pressure measurements of acoustic emissions by single-element, focused passive cavitation detectors would be facilitated by improved wideband receive calibration techniques. Here, calibration methods were developed to characterize the absolute, frequency-dependent receive sensitivity of a spherically focused, single-element transducer using pulse-echo and pitch-catch techniques. Validation of these calibration methods on a focused receiver were made by generating a pulse from a small diameter source at the focus of the transducer and comparing the absolute pressure measured by a calibrated hydrophone to that of the focused transducer using the receive sensitivities determined here.

Analysis of tissue changes, measurement system effects, and motion artifacts in echo decorrelation imaging

Echo decorrelation imaging, a method for mapping ablation-induced ultrasound echo changes, is analyzed. Local echo decorrelation is shown to approximate the decoherence spectrum of tissue reflectivity. Effects of the ultrasound measurement system, echo signal windowing, electronic noise, and tissue motion on echo decorrelation images are determined theoretically, leading to a method for reduction of motion and noise artifacts. Theoretical analysis is validated by simulations and experiments. Simulated decoherence of the scattering medium was recovered with root-mean-square error less than 10% with accuracy dependent on the correlation window size. Motion-induced decorrelation measured in an ex vivo pubovisceral muscle model showed similar trends to theoretical motion-induced decorrelation for a 2.1 MHz curvilinear array with decorrelation approaching unity for 3-4 mm elevational displacement or 1-1.6 mm range displacement. For in vivo imaging of porcine liver by a 7 MHz linear array, theoretical decorrelation computed using image-based motion estimates correlated significantly with measured decorrelation (r = 0.931, N = 10). Echo decorrelation artifacts incurred during in vivo radiofrequency ablation in the same porcine liver were effectively compensated based on the theoretical echo decorrelation model and measured pre-treatment decorrelation. These results demonstrate the potential of echo decorrelation imaging for quantification of heat-induced changes to the scattering tissue medium during thermal ablation.

Prediction and suppression of HIFU-induced vessel rupture using passive cavitation detection in an ex vivo model

Background

Occlusion of blood vessels using high-intensity focused ultrasound (HIFU) is a potential treatment for arteriovenous malformations and other neurovascular disorders. However, attempting HIFU-induced vessel occlusion can also cause vessel rupture, resulting in hemorrhage. Possible rupture mechanisms include mechanical effects of acoustic cavitation and heating of the vessel wall.

Methods

HIFU exposures were performed on 18 ex vivo porcine femoral arteries with simultaneous passive cavitation detection. Vessels were insonified by a 3.3-MHz focused source with spatial-peak, temporal-peak focal intensity of 15,690–24,430 W/cm² (peak negative-pressure range 10.92–12.52 MPa) and a 50% duty cycle for durations up to 5 min. Time-dependent acoustic emissions were recorded by an unfocused passive cavitation detector and quantified within low-frequency (10–30 kHz), broadband (0.3–1.1 MHz), and subharmonic (1.65 MHz) bands. Vessel rupture was detected by inline metering of saline flow, recorded throughout each treatment. Recorded emissions were grouped into ‘pre-rupture’ (0–10 s prior to measured point of vessel rupture) and ‘intact-vessel’ (>10 s prior to measured point of vessel rupture) emissions. Receiver operating characteristic curve analysis was used to assess the ability of emissions within each frequency band to predict vessel rupture.

Based on these measurements associating acoustic emissions with vessel rupture, a real-time feedback control module was implemented to monitor acoustic emissions during HIFU treatment and adjust the ultrasound intensity, with the goal of maximizing acoustic power delivered to the vessel while avoiding rupture. This feedback control approach was tested on 10 paired HIFU exposures of porcine femoral and subclavian arteries, in which the focal intensity was stepwise increased from 9,117 W/cm² spatial-peak temporal-peak (SPTP) to a maximum of 21,980 W/cm², with power modulated based on the measured subharmonic emission amplitude. Time to rupture was compared between these feedback-controlled trials and paired controller-inactive trials using a paired Wilcoxon signed-rank test.

Results

Subharmonic emissions were found to be the most predictive of vessel rupture (areas under the receiver operating characteristic curve (AUROC) = 0.757, p < 10^-16) compared to low-frequency (AUROC = 0.657, p < 10^-11) and broadband (AUROC = 0.729, p < 10^-16) emissions. An independent-sample t test comparing pre-rupture to intact-vessel emissions revealed a statistically significant difference between the two groups for broadband and subharmonic emissions (p < 10^-3), but not for low-frequency emissions (p = 0.058).

In a one-sided paired Wilcoxon signed-rank test, activation of the control module was shown to increase the time to vessel rupture (T_- = 8, p = 0.0244, N = 10). In one-sided paired t tests, activation of the control module was shown to cause no significant difference in time-averaged focal intensity (t = 0.362, p = 0.363, N = 10), but was shown to cause delivery of significantly greater total acoustic energy (t = 2.037, p = 0.0361, N = 10).

Conclusions

These results suggest that acoustic cavitation plays an important role in HIFU-induced vessel rupture. In HIFU treatments for vessel occlusion, passive monitoring of acoustic emissions may be useful in avoiding hemorrhage due to vessel rupture, as shown in the rupture suppression experiments.

Relations between acoustic cavitation and skin resistance during intermediate- and high-frequency sonophoresis

Enhanced skin permeability is known to be achieved during sonophoresis due to ultrasound-induced cavitation. However, the mechanistic role of cavitation during sonophoresis has been extensively investigated only for low-frequency (LFS, <100 kHz) applications. Here, mechanisms of permeability-enhancing stable and inertial cavitation were investigated by passively monitoring subharmonic and broadband emissions arising from cavitation isolated within or external to porcine skin in vitro during intermediate- (IFS, 100-700 kHz) and high-frequency sonophoresis (HFS, >1 MHz). The electrical resistance of skin, a surrogate measure of the permeability of skin to a variety of compounds, was measured to quantify the reduction and subsequent recovery of the skin barrier during and after exposure to pulsed (1 second pulse, 20% duty cycle) 0.41 and 2.0 MHz ultrasound over a range of acoustic powers (0-21.7 W) for 30 min. During IFS, significant skin resistance reductions and acoustic emissions from cavitation were measured exclusively when cavitation was isolated outside of the skin. Time-dependent skin resistance reductions measured during IFS correlated significantly with subharmonic and broadband emission levels. During HFS, significant skin resistance reductions were accompanied by significant acoustic emissions from cavitation measured during trials that isolated cavitation activity either outside of skin or within skin. Time-dependent skin resistance reductions measured during HFS correlated significantly greater with subharmonic than with broadband emission levels. The reduction of the skin barrier due to sonophoresis was reversible in all trials; however, effects incurred during IFS recovered more slowly and persisted over a longer period of time than HFS. These results quantitatively demonstrate the significance of cavitation during sonophoresis and suggest that the mechanisms and post-treatment longevity of permeability enhancement due to IFS and HFS treatments are different.

In vivo thermal ablation monitoring using ultrasound echo decorrelation imaging

Previous work indicated that ultrasound echo decorrelation imaging can track and quantify changes in echo signals to predict thermal damage during in vitro radiofrequency ablation (RFA). In the in vivo studies reported here, the feasibility of using echo decorrelation imaging as a treatment monitoring tool was assessed. RFA was performed on normal swine liver (N = 5), and ultrasound ablation using image-ablate arrays was performed on rabbit liver implanted with VX2 tumors (N = 2). Echo decorrelation and integrated backscatter were computed from Hilbert transformed pulse-echo data acquired during RFA and ultrasound ablation treatments. Receiver operating characteristic (ROC) curves were employed to assess the ability of echo decorrelation imaging and integrated backscatter to predict ablation. Area under the ROC curves (AUROC) was determined for RFA and ultrasound ablation using echo decorrelation imaging. Ablation was predicted more accurately using echo decorrelation imaging (AUROC = 0.832 and 0.776 for RFA and ultrasound ablation, respectively) than using integrated backscatter (AUROC = 0.734 and 0.494).

Passive imaging with pulsed ultrasound insonations

Previously, passive cavitation imaging has been described in the context of continuous-wave high-intensity focused ultrasound thermal ablation. However, the technique has potential use as a feedback mechanism for pulsed-wave therapies, such as ultrasound-mediated drug delivery. In this paper, results of experiments and simulations are reported to demonstrate the feasibility of passive cavitation imaging using pulsed ultrasound insonations and how the images depend on pulsed ultrasound parameters. The passive cavitation images were formed from channel data that was beamformed in the frequency domain. Experiments were performed in an invitro flow phantom with an experimental echo contrast agent, echogenic liposomes, as cavitation nuclei. It was found that the pulse duration and envelope have minimal impact on the image resolution achieved. The passive cavitation image amplitude scales linearly with the cavitation emission energy. Cavitation images for both stable and inertial cavitation can be obtained from the same received data set.

Acoustic characterization of echogenic liposomes: frequency-dependent attenuation and backscatter

Ultrasound contrast agents (UCAs) are used clinically to aid detection and diagnosis of abnormal blood flow or perfusion. Characterization of UCAs can aid in the optimization of ultrasound parameters for enhanced image contrast. In this study echogenic liposomes (ELIPs) were characterized acoustically by measuring the frequency-dependent attenuation and backscatter coefficients at frequencies between 3 and 30 MHz using a broadband pulse-echo technique. The experimental methods were initially validated by comparing the attenuation and backscatter coefficients measured from 50-μm and 100-μm polystyrene microspheres with theoretical values. The size distribution of the ELIPs was measured and found to be polydisperse, ranging in size from 40 nm to 6 μm in diameter, with the highest number observed at 65 nm. The ELIP attenuation coefficients ranged from 3.7  ±  1.0 to 8.0  ±  3.3 dB/cm between 3 and 25 MHz. The backscatter coefficients were 0.011  ±  0.006 (cm str)(-1) between 6 and 9 MHz and 0.023 ± 0.006 (cm str)(-1) between 13 and 30 MHz. The measured scattering-to-attenuation ratio ranged from 8% to 22% between 6 and 25 MHz. Thus ELIPs can provide enhanced contrast over a broad range of frequencies and the scattering properties are suitable for various ultrasound imaging applications including diagnostic and intravascular ultrasound.

Treatment of rabbit liver cancer in vivo using miniaturized image-ablate ultrasound arrays

In the preclinical studies reported here, VX2 cancer within rabbit liver has been treated by bulk ultrasound ablation employing miniaturized image-ablate arrays. Array probes were constructed with 32 elements in a 2.3 × 20 mm(2) aperture, packaged within a 3.1 mm stainless steel tube with a cooling and coupling balloon for in vivo use. The probes were measured capable of 50% fractional bandwidth for pulse-echo imaging (center frequency 4.4 MHz) with >110 W/cm(2) surface intensity available at sonication frequencies 3.5 and 4.8 MHz. B-scan imaging performance of the arrays was measured to be comparable to larger diagnostic linear arrays, although nearfield image quality was reduced by ringdown artifacts. A series of in vivo ablation procedures was performed using an unfocused 32-element aperture firing at 4.8 MHz with exposure durations 20-70.5 s and in situ spatial average, temporal average intensities 22.4-38.5 W/cm(2). Ablation of a complete tumor cross-section was confirmed by vital staining in seven of 12 exposures, with four exposures ablating an additional margin >1 mm beyond the tumor in all directions. Analysis suggests a threshold ablation effect, with complete ablation of tumor cross-sections for exposures with delivery of >838 J acoustic energy. The results show feasibility for in vivo liver cancer ablation using miniaturized image-ablate arrays suitable for interstitial deployment.

Convolutional modeling of diffraction effects in pulse-echo ultrasound imaging

A model is presented for pulse-echo imaging of three-dimensional, linear, weakly-scattering continuum media by ultrasound array transducers. The model accounts for the diffracted fields of focused array subapertures in both transmit and receive modes, multiple transmit and receive focal zones, frequency-dependent attenuation, and aberration caused by mismatched medium and beamformer sound speeds. For a given medium reflectivity function, computation of a B-scan requires evaluation of a depth-dependent transmit/receive beam product, followed by two one-dimensional convolutions and a one-dimensional summation. Numerical results obtained using analytic expressions for transmit and receive beams agree favorably with measured B-scan images and speckle statistics.

Passive cavitation imaging with ultrasound arrays

A method is presented for passive imaging of cavitational acoustic emissions using an ultrasound array, with potential application in real-time monitoring of ultrasound ablation. To create such images, microbubble emissions were passively sensed by an imaging array and dynamically focused at multiple depths. In this paper, an analytic expression for a passive image is obtained by solving the Rayleigh-Sommerfield integral, under the Fresnel approximation, and passive images were simulated. A 192-element array was used to create passive images, in real time, from 520-kHz ultrasound scattered by a 1-mm steel wire. Azimuthal positions of this target were accurately estimated from the passive images. Next, stable and inertial cavitation was passively imaged in saline solution sonicated at 520 kHz. Bubble clusters formed in the saline samples were consistently located on both passive images and B-scans. Passive images were also created using broadband emissions from bovine liver sonicated at 2.2 MHz. Agreement was found between the images and source beam shape, indicating an ability to map therapeutic ultrasound beams in situ. The relation between these broadband emissions, sonication amplitude, and exposure conditions are discussed.

Ultrasound monitoring of in vitro radio frequency ablation by echo decorrelation imaging

OBJECTIVE

The purpose of this study was to test ultrasound echo decorrelation imaging for mapping and characterization of tissue effects caused by radio frequency ablation (RFA).

METHODS

Radio frequency ablation procedures (6-minute duration, 20-W power) were performed on fresh ex vivo bovine liver tissue (n = 9) with continuous acquisition of beam-formed ultrasound echo data from a 7-MHz linear array. Echo data were processed to form B-scan images, echo decorrelation images (related to rapid random changes in echo waveforms), and integrated backscatter images (related to local changes in received echo energy). Echo decorrelation and integrated backscatter values at the location of a low-noise thermocouple were assessed as functions of temperature. Echo decorrelation and integrated backscatter images were directly compared with ablated tissue cross sections and quantitatively evaluated as predictors of tissue ablation and overtreatment.

RESULTS

Echo decorrelation maps corresponded with local tissue temperature and ablation effects. Consistent echo decorrelation increases were observed for temperatures above 75 degrees C, whereas integrated backscatter maps showed a nonmonotonic temperature dependence complicated by acoustic shadowing, with high variance at large temperature elevations. In receiver operating characteristic curve analysis of echo decorrelation and integrated backscatter maps as predictors of local tissue ablation, echo decorrelation performed well (area under the receiver operating characteristic curve [AUROC] = 0.855 for ablation and 0.913 for overtreatment), whereas integrated backscatter performed poorly (AUROC < 0.6).

CONCLUSIONS

Echo decorrelation imaging can map tissue changes due to RFA in vitro, with local echo decorrelation corresponding strongly to local tissue temperature elevations and ablation effects. With further development and in vivo validation, echo decorrelation imaging is potentially useful for improved image guidance of clinical RFA procedures.

Characterization of ultrasound propagation through ex-vivo human temporal bone

Adjuvant therapies that lower the thrombolytic dose or increase its efficacy would represent a significant breakthrough in the treatment of patients with ischemic stroke. The objective of this study was to perform intracranial measurements of the acoustic pressure field generated by 0.12, 1.03 and 2.00-MHz ultrasound transducers to identify optimal ultrasound parameters that would maximize penetration and minimize aberration of the beam. To achieve this goal, in vitro experiments were conducted on five human skull specimens. In a water-filled tank, two unfocused transducers (0.12 and 1.03 MHz) and one focused transducer (2.00 MHz) were consecutively placed near the right temporal bone of each skull. A hydrophone, mounted on a micropositioning system, was moved to an estimated location of the middle cerebral artery (MCA) origin, and measurements of the surrounding acoustic pressure field were performed. For each measurement, the distance from the position of maximum acoustic pressure to the estimated origin of the MCA inside the skulls was quantified. The -3 dB depth-of-field and beamwidth in the skull were also investigated as a function of the three frequencies. Results show that the transducer alignment relative to the skull is a significant determinant of the detailed behavior of the acoustic field inside the skull. For optimal penetration, insonation normal to the temporal bone was needed. The shape of the 0.12-MHz intracranial beam was more distorted than those at 1.03 and 2.00 MHz because of the large aperture and beamwidth. However, lower ultrasound pressure reduction was observed at 0.12 MHz (22.5%). At 1.03 and 2.00 MHz, two skulls had an insufficient temporal bone window and attenuated the beam severely (up to 96.6% pressure reduction). For all frequencies, constructive and destructive interference patterns were seen near the contralateral skull wall at various elevations. The 0.12-MHz ultrasound beam depth-of-field was affected the most when passing through the temporal bone and showed a decrease in size of more than 55% on average. The speed of sound in the temporal bone of each skull was estimated at 1.03 MHz and demonstrated a large range (1752.1 to 3285.3 m/s). Attenuation coefficients at 1.03 and 2.00 MHz were also derived for each of the five skull specimens. This work provides needed information on ultrasound beam shapes inside the human skull, which is a necessary first step for the development of an optimal transcranial ultrasound-enhanced thrombolysis device.

Acoustic emissions during 3.1 MHz ultrasound bulk ablation in vitro

Acoustic emissions associated with cavitation and other bubble activity have previously been observed during ultrasound (US) ablation experiments. Because detectable bubble activity may be related to temperature, tissue state and sonication characteristics, these acoustic emissions are potentially useful for monitoring and control of US ablation. To investigate these relationships, US ablation experiments were performed with simultaneous measurements of acoustic emissions, tissue echogenicity and tissue temperature on fresh bovine liver. Ex vivo tissue was exposed to 0.9-3.3-s bursts of unfocused, continuous-wave, 3.10-MHz US from a miniaturized 32-element array, which performed B-scan imaging with the same piezoelectric elements during brief quiescent periods. Exposures used pressure amplitudes of 0.8-1.4 MPa for exposure times of 6-20 min, sufficient to achieve significant thermal coagulation in all cases. Acoustic emissions received by a 1-MHz, unfocused passive cavitation detector, beamformed A-line signals acquired by the array, and tissue temperature detected by a needle thermocouple were sampled 0.3-1.1 times per second. Tissue echogenicity was quantified by the backscattered echo energy from a fixed region-of-interest within the treated zone. Acoustic emission levels were quantified from the spectra of signals measured by the passive cavitation detector, including subharmonic signal components at 1.55 MHz, broadband signal components within the band 0.3-1.1 MHz and low-frequency components within the band 10-30 kHz. Tissue ablation rates, defined as the thermally ablated volumes per unit time, were assessed by quantitative analysis of digitally imaged, macroscopic tissue sections. Correlation analysis was performed among the averaged and time-dependent acoustic emissions in each band considered, B-mode tissue echogenicity, tissue temperature and ablation rate. Ablation rate correlated significantly with broadband and low-frequency emissions, but was uncorrelated with subharmonic emissions. Subharmonic emissions were found to depend strongly on temperature in a nonlinear manner, with significant emissions occurring within different temperature ranges for each sonication amplitude. These results suggest potential roles for passive detection of acoustic emissions in guidance and control of bulk US ablation treatments.

Ultrasound-enhanced thrombolysis using Definity as a cavitation nucleation agent

Ultrasound has been shown previously to act synergistically with a thrombolytic agent, such as recombinant tissue plasminogen activator (rt-PA) to accelerate thrombolysis. In this in vitro study, a commercial contrast agent, Definity, was used to promote and sustain the nucleation of cavitation during pulsed ultrasound exposure at 120 kHz. Ultraharmonic signals, broadband emissions and harmonics of the fundamental were measured acoustically by using a focused hydrophone as a passive cavitation detector and used to quantify the level of cavitation activity. Human whole blood clots suspended in human plasma were exposed to a combination of rt-PA, Definity and ultrasound at a range of ultrasound peak-to-peak pressure amplitudes, which were selected to expose clots to various degrees of cavitation activity. Thrombolytic efficacy was determined by measuring clot mass loss before and after the treatment and correlated with the degree of cavitation activity. The penetration depth of rt-PA and plasminogen was also evaluated in the presence of cavitating microbubbles using a dual-antibody fluorescence imaging technique. The largest mass loss (26.2%) was observed for clots treated with 120-kHz ultrasound (0.32-MPa peak-to-peak pressure amplitude), rt-PA and stable cavitation nucleated by Definity. A significant correlation was observed between mass loss and ultraharmonic signals (r = 0.85, p < 0.0001, n = 24). The largest mean penetration depth of rt-PA (222 microm) and plasminogen (241 microm) was observed in the presence of stable cavitation activity. Stable cavitation activity plays an important role in enhancement of thrombolysis and can be monitored to evaluate the efficacy of thrombolytic treatment.

Ultrasound-induced thermal elevation in clotted blood and cranial bone

Ultrasound thermal effects have been hypothesized to contribute to ultrasound-assisted thrombolysis. To explore the thermal mechanism of ultrasound-enhanced thrombolysis with recombinant tissue plasminogen activator (rt-PA) for the treatment of ischemic stroke, a detailed investigation is needed of the heating produced in skull, brain and blood clots. A theoretical model is developed to provide an estimate for the worst-case scenario of the temperature increase in blood clots and on the surface of cranial bone exposed to 0.12- to 3.5-MHz ultrasound. Thermal elevation was also assessed experimentally in human temporal bone, human clots and porcine clots exposed to 0.12 to 3.5-MHz pulsed ultrasound in vitro with a peak-to-peak pressure of 0.25 MPa and 80% duty cycle. Blood clots exposed to 0.12-MHz pulsed ultrasound exhibited a small temperature increase (0.25 degrees C) and bone exposed to 1.0-MHz pulsed ultrasound exhibited the highest temperature increase (1.0 degrees C). These experimental results were compared with the predicted temperature elevations.

Fresnel approximations for acoustic fields of rectangularly symmetric sources

A general approach is presented for determining the acoustic fields of rectangularly symmetric, baffled, time-harmonic sources under the Fresnel approximation. This approach is applicable to a variety of separable source configurations, including uniform, exponential, Gaussian, sinusoidal, and error function surface velocity distributions, with and without focusing in either surface dimension. In each case, the radiated field is given by a formula similar to that for a uniform rectangular source, except for additional scaling of wave number and azimuthal distance parameters. The expressions presented are generalized to three different Fresnel approximations that correspond, respectively, to diffracted plane waves, diffracted spherical waves, or diffracted cylindrical waves. Numerical results, for several source geometries relevant to ultrasonic applications, show that these expressions accurately depict the radiated pressure fields, except for points very near the radiating aperture. Highest accuracy near the source is obtained by choice of the Fresnel approximation most suited to the source geometry, while the highest accuracy far from the source is obtained by the approximation corresponding to diffracted spherical waves. The methods are suitable for volumetric computations of acoustic fields including focusing, apodization, and attenuation effects.

Miniaturized ultrasound arrays for interstitial ablation and imaging

A potential alternative to extracorporeal, noninvasive HIFU therapy is minimally invasive intense ultrasound ablation that can be performed laparoscopically or percutaneously. An approach to minimally invasive ablation of soft tissue using miniaturized linear ultrasound arrays is presented here. Recently developed 32-element arrays with aperture 2.3 x 49 mm, therapy frequency 3.1 MHz, pulse-echo bandwidths >42% and surface acoustic energy density >80 W/cm2, are described. These arrays are integrated into a probe assembly, including a coupling balloon and piercing tip, suitable for interstitial ablation. An integrated electronic control system allows therapy planning and automated treatment guided by real-time interstitial B-scan imaging. Image quality, challenging because of limited probe dimensions and channel count, is aided by signal processing techniques that improve image definition and contrast, resulting in image quality comparable to typical transabdominal ultrasound imaging. Ablation results from ex vivo and in vivo experiments on mammalian liver tissue show that this approach is capable of ablation rates and volumes relevant to clinical applications of soft tissue ablation such as treatment of liver cancer.

Bulk ablation of soft tissue with intense ultrasound: modeling and experiments

Methods for the bulk ablation of soft tissue using intense ultrasound, with potential applications in the thermal treatment of focal tumors, are presented. An approximate analytic model for bulk ablation predicts the progress of ablation based on tissue properties, spatially averaged ultrasonic heat deposition, and perfusion. The approximate model allows the prediction of threshold acoustic powers required for ablation in vivo as well as the comparison of cases with different starting temperatures and perfusion characteristics, such as typical in vivo and ex vivo experiments. In a full three-dimensional numerical model, heat deposition from array transducers is computed using the Fresnel approximation and heat transfer in tissue is computed by finite differences, accounting for heating changes caused by boiling and thermal dose-dependent absorption. Similar ablation trends due to perfusion effects are predicted by both the simple analytic model and the full numerical model. Comparisons with experimental results show the efficacy of both models in predicting tissue ablation effects. Phenomena illustrated by the simulations and experiments include power thresholds for in vivo ablation, differences between in vivo and ex vivo lesioning for comparable source conditions, the effect of tissue boiling and absorption changes on ablation depth, and the performance of a continuous rotational scanning method suitable for interstitial bulk ablation of soft tissue.

Simulation of ultrasonic focus aberration and correction through human tissue

Ultrasonic focusing in two dimensions has been investigated by calculating the propagation of ultrasonic pulses through cross-sectional models of human abdominal wall and breast. Propagation calculations used a full-wave k-space method that accounts for spatial variations in density, sound speed, and frequency-dependent absorption and includes perfectly matched layer absorbing boundary conditions. To obtain a distorted receive wavefront, propagation from a point source through the tissue path was computed. Receive focusing used an angular spectrum method. Transmit focusing was accomplished by propagating a pressure wavefront from a virtual array through the tissue path. As well as uncompensated focusing, focusing that employed time-shift compensation and time-shift compensation after backpropagation was investigated in both transmit and receive and time reversal was investigated for transmit focusing in addition. The results indicate, consistent with measurements, that breast causes greater focus degradation than abdominal wall. The investigated compensation methods corrected the receive focus better than the transmit focus. Time-shift compensation after backpropagation improved the focus from that obtained using time-shift compensation alone but the improvement was less in transmit focusing than in receive focusing. Transmit focusing by time reversal resulted in lower sidelobes but larger mainlobes than the other investigated transmit focus compensation methods.

Aberration correction for time-domain ultrasound diffraction tomography

Extensions of a time-domain diffraction tomography method, which reconstructs spatially dependent sound speed variations from far-field time-domain acoustic scattering measurements, are presented and analyzed. The resulting reconstructions are quantitative images with applications including ultrasonic mammography, and can also be considered candidate solutions to the time-domain inverse scattering problem. Here, the linearized time-domain inverse scattering problem is shown to have no general solution for finite signal bandwidth. However, an approximate solution to the linearized problem is constructed using a simple delay-and-sum method analogous to "gold standard" ultrasonic beamforming. The form of this solution suggests that the full nonlinear inverse scattering problem can be approximated by applying appropriate angle- and space-dependent time shifts to the time-domain scattering data; this analogy leads to a general approach to aberration correction. Two related methods for aberration correction are presented: one in which delays are computed from estimates of the medium using an efficient straight-ray approximation, and one in which delays are applied directly to a time-dependent linearized reconstruction. Numerical results indicate that these correction methods achieve substantial quality improvements for imaging of large scatterers. The parametric range of applicability for the time-domain diffraction tomography method is increased by about a factor of 2 by aberration correction.

Frequency dependence of speckle in continuous-wave ultrasound with implications for blood perfusion measurements

Speckle in continuous wave (CW) Doppler has previously been found to cause large variations in detected Doppler power in blood perfusion measurements, where a large number of blood vessels are present in the sample volume. This artifact can be suppressed by using a number of simultaneously transmitted frequencies and averaging the detected signals. To optimize the strategy, statistical properties of speckle in CW ultrasound need to be known. This paper presents analysis of the frequency separation necessary to obtain independent values of the received power for CW ultrasound using a simplified mathematical model for insonation of a static, lossless, statistically homogeneous, weakly scattering medium. Specifically, the autocovariance function for received power is derived, which functionally is the square of the (deterministic) autocorrelation function of the effective sample volumes produced by the transducer pair for varying frequencies, at least if a delta correlated medium is assumed. A marginal broadening of the modeled autocovariance functions is expected for insonation of blood. The theory is applicable to any transducer aperture, but has been experimentally verified here with 5-MHz, 6.35-mm circular transducers using an agar phantom containing small, randomly dispersed glass particles. A similar experimental verification of a transducer used in multiple-frequency blood perfusion measurements shows that the model proposed in this paper is plausible for explaining the decorrelation between different channels in such a measurement.

A k-space method for coupled first-order acoustic propagation equations

A k-space method for large-scale simulation of ultrasonic pulse propagation is presented. The present method, which solves the coupled first-order differential equations for wave propagation in inhomogeneous media, is derived in a simple form analogous to previous finite-difference methods with staggered spatial and temporal grids. Like k-space methods based on second-order wave equations, the present method is exact for homogeneous media, unconditionally stable for "slow" [c(r) < or = c0] media, and highly accurate for general weakly scattering media. In addition, unlike previous k-space methods, the form of the method allows straightforward inclusion of relaxation absorption and perfectly matched layer (PML) nonreflecting boundary conditions. Numerical examples illustrate the capabilities of the present k-space method. For weakly inhomogeneous media, accurate results are obtained using coarser temporal and spatial steps than possible with comparable finite-difference and pseudospectral methods. The low dispersion of the k-space method allows accurate representation of frequency-dependent attenuation and phase velocity associated with relaxation absorption. A technique for reduction of Gibbs phenomenon artifacts, in which compressibility and exponentially scaled density functions are smoothed by half-band filtering, is introduced. When employed together with this smoothing technique, the k-space method provides high accuracy for media including discontinuities, high-contrast inhomogeneities, and scattering structures smaller than the spatial grid resolution.

A k-space method for large-scale models of wave propagation in tissue

Large-scale simulation of ultrasonic pulse propagation in inhomogeneous tissue is important for the study of ultrasound-tissue interaction as well as for development of new imaging methods. Typical scales of interest span hundreds of wavelengths; most current two-dimensional methods, such as finite-difference and finite-element methods, are unable to compute propagation on this scale with the efficiency needed for imaging studies. Furthermore, for most available methods of simulating ultrasonic propagation, large-scale, three-dimensional computations of ultrasonic scattering are infeasible. Some of these difficulties have been overcome by previous pseudospectral and k-space methods, which allow substantial portions of the necessary computations to be executed using fast Fourier transforms. This paper presents a simplified derivation of the k-space method for a medium of variable sound speed and density; the derivation clearly shows the relationship of this k-space method to both past k-space methods and pseudospectral methods. In the present method, the spatial differential equations are solved by a simple Fourier transform method, and temporal iteration is performed using a k-t space propagator. The temporal iteration procedure is shown to be exact for homogeneous media, unconditionally stable for "slow" (c(x) < or = c0) media, and highly accurate for general weakly scattering media. The applicability of the k-space method to large-scale soft tissue modeling is shown by simulating two-dimensional propagation of an incident plane wave through several tissue-mimicking cylinders as well as a model chest wall cross section. A three-dimensional implementation of the k-space method is also employed for the example problem of propagation through a tissue-mimicking sphere. Numerical results indicate that the k-space method is accurate for large-scale soft tissue computations with much greater efficiency than that of an analogous leapfrog pseudospectral method or a 2-4 finite difference time-domain method. However, numerical results also indicate that the k-space method is less accurate than the finite-difference method for a high contrast scatterer with bone-like properties, although qualitative results can still be obtained by the k-space method with high efficiency. Possible extensions to the method, including representation of absorption effects, absorbing boundary conditions, elastic-wave propagation, and acoustic nonlinearity, are discussed.

Quantitative flaw reconstruction from ultrasonic surface wavefields measured by electronic speckle pattern interferometry

A new method for imaging flaws in plate and shell structures is presented. The method employs two-dimensional ultrasonic surface wave data obtained by optical electronic speckle pattern interferometry (ESPI) techniques. In the imaging method, the measured out-of-plane displacement field associated with an externally excited ultrasonic Lamb wave is processed to obtain the spatial frequency domain spectrum of the wavefield. A free space Green's function is then deconvolved from the wavefield to obtain quantitative images of effective scattering sources. Because the strength of these effective sources is directly dependent on local variations in sample thickness and material properties, these images provide a direct map of internal inhomogeneities. Simulation results show that the method accurately images flaws for a wide range of sizes and material contrast ratios. These results also demonstrate that flaw features much smaller than an acoustic wavelength can be imaged, consistent with the theoretical capability of the imaging method to employ scattered evanescent waves. Reconstructions are also obtained from ultrasonic Lamb wave displacement fields recorded by ESPI in a flawed aluminum plate. These reconstructions indicate that the present method has potential for imaging flaws in complex structures for which ESPI wavefield measurements cannot be straightforwardly interpreted.

Simulation of ultrasonic pulse propagation, distortion, and attenuation in the human chest wall

A finite-difference time-domain model for ultrasonic pulse propagation through soft tissue has been extended to incorporate absorption effects as well as longitudinal-wave propagation in cartilage and bone. This extended model has been used to simulate ultrasonic propagation through anatomically detailed representations of chest wall structure. The inhomogeneous chest wall tissue is represented by two-dimensional maps determined by staining chest wall cross sections to distinguish between tissue types, digitally scanning the stained cross sections, and mapping each pixel of the scanned images to fat, muscle, connective tissue, cartilage, or bone. Each pixel of the tissue map is then assigned a sound speed, density, and absorption value determined from published measurements and assumed to be representative of the local tissue type. Computational results for energy level fluctuations and arrival time fluctuations show qualitative agreement with measurements performed on the same specimens, but show significantly less waveform distortion than measurements. Visualization of simulated tissue-ultrasound interactions in the chest wall shows possible mechanisms for image aberration in echocardiography, including effects associated with reflection and diffraction caused by rib structures. A comparison of distortion effects for varying pulse center frequencies shows that, for soft tissue paths through the chest wall, energy level and waveform distortion increase markedly with rising ultrasonic frequency and that arrival-time fluctuations increase to a lesser degree.

The effect of abdominal wall morphology on ultrasonic pulse distortion. Part II. Simulations

Wavefront propagation through the abdominal wall was simulated using a finite-difference time-domain implementation of the linearized wave propagation equations for a lossless, inhomogeneous, two-dimensional fluid as well as a simplified straight-ray model for a two-dimensional absorbing medium. Scanned images of six human abdominal wall cross sections provided the data for the propagation media in the simulations. The images were mapped into regions of fat, muscle, and connective tissue, each of which was assigned uniform sound speed, density, and absorption values. Propagation was simulated through each whole specimen as well as through each fat layer and muscle layer individually. Wavefronts computed by the finite-difference method contained arrival time, energy level, and wave shape distortion similar to that in measurements. Straight-ray simulations produced arrival time fluctuations similar to measurements but produced much smaller energy level fluctuations. These simulations confirm that both fat and muscle produce significant wavefront distortion and that distortion produced by fat sections differs from that produced by muscle sections. Spatial correlation of distortion with tissue composition suggests that most major arrival time fluctuations are caused by propagation through large-scale inhomogeneities such as fatty regions within muscle layers, while most amplitude and waveform variations are the result of scattering from smaller inhomogeneities such as septa within the subcutaneous fat. Additional finite-difference simulations performed using uniform-layer models of the abdominal wall indicate that wavefront distortion is primarily caused by tissue structures and inhomogeneities rather than by refraction at layer interfaces or by variations in layer thicknesses.

The effect of abdominal wall morphology on ultrasonic pulse distortion. Part I. Measurements

The relative importance of the fat and muscle layers of the human abdominal wall in producing ultrasonic wavefront distortion was assessed by means of direct measurements. Specimens employed included six whole abdominal wall specimens and twelve partial specimens obtained by dividing each whole specimen into a fat and a muscle layer. In the measurement technique employed, a hemispheric transducer transmitted a 3.75-MHz ultrasonic pulse through a tissue section. The received wavefront was measured by a linear array translated in the elevation direction to synthesize a two-dimensional aperture. Insertion loss was also measured at various locations on each specimen. Differences in arrival time and energy level between the measured waveforms and computed references that account for geometric delay and spreading were calculated. After correction for the effects of geometry, the received waveforms were synthetically focused. The characteristics of the distortion produced by each specimen and the quality of the resulting focus were analyzed and compared. The measurements show that muscle produces greater arrival time distortion than fat while fat produces greater energy level distortion than muscle, but that the distortion produced by the entire abdominal wall is not equivalent to a simple combination of distortion effects produced by the layers. The results also indicate that both fat and muscle layers contribute significantly to the distortion of ultrasonic beams by the abdominal wall. However, the spatial characteristics of the distortion produced by fat and muscle layers differ substantially. Distortion produced by muscle layers, as well as focal images aberrated by muscle layers, show considerable anisotropy associated with muscle fiber orientation. Distortion produced by fat layers shows smaller-scale, granular structure associated with scattering from the septa surrounding individual fat lobules. Thick layers of fat may be expected to cause poor image quality due to both scattering and bulk absorption effects, while thick muscle layers may be expected to cause focus aberration due to large arrival time fluctuations. Correction of aberrated focuses using time-shift compensation shows more complete correction for muscle sections than for fat sections, so that correction methods based on phase screen models may be more appropriate for muscle layers than for fat layers.