Method for Recovering Lost Ultrasonic Information Using the Echo-intensity Mean.

The axial resolution of a B-mode (or intensity) image is limited by the bandwidth of the pulse envelope. In this report, we investigate the source of this limitation by examining the transfer of high-resolution information from the tissue impedance variance throughout the imaging process. For that purpose, we express the mean and variance of the echo-intensity signal as a linear system to track the flow of object information along the image-formation chain. The results reveal how demodulation influences the available information by discarding high spatial-frequency information. This analysis further points to a simple way to recover lost information with only a minor addition to signal processing. Software phantoms are used to show that under ideal conditions, information from small-scale high-contrast reflectors, such as microcalcifications, can be significantly enhanced with this simple change to echo processing.

Improvements to ultrasonic beamformer design and implementation derived from the task-based analytical framework.

The task-based framework, previously developed for beamformer comparison [Nguyen, Prager, and Insana, J. Acoust. Soc. Am. 140, 1048-1059 (2016)], is extended to design a new beamformer with potential applications in breast cancer diagnosis. The beamformer is based on a better approximation of the Bayesian strategy. It is a combination of the Wiener-filtered beamformer and an iterative process that adapts the generated image to specific features of the object. Through numerical studies, the new method is shown to outperform other beamformers drawn from the framework, but at an increase in computational cost. It requires a preprocessing step where the scattering field is segmented into regions with distinct statistical properties. Segmentation errors become a major limitation to the beamformer performance. All the beamformers under investigation are tested using data obtained from an instrumented ultrasound machine. They are implemented using a new time delay calculation, recently developed in the pixel-based beamforming studies presented here, which helps to overcome the challenge posed by the shift-variant nature of the imaging system. The efficacy of each beamformer is evaluated based on the quality of generated images in the context of the task-based framework. The in vitro results confirm the conclusions drawn from the simulations.

Limitations on estimation of effective scatterer diameters.

The influence of spatial diversity in acoustic scattering properties on estimates of the effective scatterer diameter (ESD) applied to soft biological tissues is investigated. This study is based on two-dimensional simulations of scattering media, beginning with random distributions of simple disk structures where all scattering features are known exactly. It concludes with an analysis of histology maps from healthy and fatty rabbit liver. Further, the liver histology is decomposed using an orthonormal basis to separate acoustic scattering at various spatial scales and observe how it influences ESD estimates. Overall, the goal is to quantitatively interpret ESD results for diagnostic assessments despite wide variations in tissue scatterer properties.

A task-based analytical framework for ultrasonic beamformer comparison.

A task-based approach is employed to develop an analytical framework for ultrasound beamformer design and evaluation. In this approach, a Bayesian ideal-observer provides an idealized starting point and a way to measure information loss in practical beamformer designs. Different approximations of this ideal strategy are shown to lead to popular beamformers in the literature, including the matched filter, minimum variance (MV), and Wiener filter (WF) beamformers. Analysis of the approximations indicates that the WF beamformer should outperform the MV approach, especially in low echo signal-to-noise conditions. The beamformers are applied to five typical tasks from the BIRADS lexicon. Their performance is evaluated based on ability to discriminate idealized malignant and benign features. The numerical results show the advantages of the WF over the MV technique in general; although performance varies predictably in some contrast-limited tasks because of the model modifications required for the MV algorithm to avoid ill-conditioning.

Indentation Measurements to Validate Dynamic Elasticity Imaging Methods.

We describe macro-indentation techniques for estimating the elastic modulus of soft hydrogels. Our study describes (a) conditions under which quasi-static indentation can validate dynamic shear-wave imaging estimates and (b) how each of these techniques uniquely biases modulus estimates as they couple to the sample geometry. Harmonic shear waves between 25 and 400 Hz were imaged using ultrasonic Doppler and optical coherence tomography methods to estimate shear dispersion. From the shear-wave speed of sound, average elastic moduli of homogeneous samples were estimated. These results are compared directly with macroscopic indentation measurements measured two ways. One set of measurements applied Hertzian theory to the loading phase of the force-displacement curves using samples treated to minimize surface adhesion forces. A second set of measurements applied Johnson-Kendall-Roberts theory to the unloading phase of the force-displacement curve when surface adhesions were significant. All measurements were made using gelatin hydrogel samples of different sizes and concentrations. Agreement within 5% among elastic modulus estimates was achieved for a range of experimental conditions. Consequently, a simple quasi-static indentation measurement using a common gel can provide elastic modulus measurements that help validate dynamic shear-wave imaging estimates.

Stromal-epithelial responses to fractionated radiotherapy in a breast cancer microenvironment.

BACKGROUND: The stromal-epithelial-cell interactions that are responsible for directing normal breast-tissue development and maintenance play a central role in the progression of breast cancer. In the present study, we continued our development of three-dimensional (3-D) cell co-cultures used to study cancerous mammary cell responses to fractionated radiotherapy. In particular, we focused on the role of the reactive stroma in determining the therapeutic ratio for post-surgical treatment. METHODS: Cancerous human mammary epithelial cells (MRC-7) were cultured in a 3-D collagen matrix with human fibroblasts (MRC-5) stimulated by various concentrations of transforming growth factor beta 1 (TGF-ß1). These culture samples were designed to model the post-lumpectomy mammary stroma in the presence of residual cancer cells. We tracked over time the changes in medium stiffness, fibroblast-cell activation (MRC-5 converted to cancer activated fibroblasts (CAFs)), and proliferation of both cell types under a variety of fractionated radiotherapy protocols. Samples were exposed to 6 MV X-rays from a linear accelerator in daily fraction sizes of 90, 180 and 360 cGy over five days in a manner consistent with irradiation exposure during radiotherapy. RESULTS: We found in fractionation studies with MRC-5 fibroblasts and CAFs that higher doses per fraction may be more effective early on in deactivating cancer-harboring cellular environments. Higher-dose fraction schemes inhibit contractility in CAFs and prevent differentiation of fibroblasts, thereby metabolically uncoupling tumor cells from their surrounding stroma. However, higher dose fraction appears to increase ECM stiffening. CONCLUSIONS: The findings suggest that dose escalation to the region with residual disease can deactivate the reactive stroma, thus minimizing the cancer promoting features of the cellular environment. Large-fraction irradiation may be used to sterilize residual tumor cells and inhibit activation of intracellular transduction pathways that are promoted during the post-surgical wound-healing period. The higher dose fractions may slow wound healing and increase ECM stiffening that could stimulate proliferation of surviving cancer cells.

Power-Doppler Ultrasonic Imaging of Peripheral Perfusion in Diabetic Mice.

OBJECTIVE: We explored the capabilities of power-Doppler ultrasonic (PD-US) imaging without contrast enhancement for monitoring changes in muscle perfusion over time. METHODS: Ischemic recovery was observed in healthy and type II diabetic male and female mice with and without exercise. In separate studies, perfusion was measured during and after 5-min ischemic periods and during four-week recovery periods following irreversible femoral ligation. A goal was to assess how well PD-US estimates tracked the diabetic-related changes in endothelial function that influenced perfusion. RESULTS: The average perfusion recovery time following femoral ligation increased 47% in diabetic males and 74% in diabetic females compared with non-diabetic mice. Flow-mediated dilation in conduit arteries and the reactive hyperemia index in resistive vessels each declined by one half in sedentary diabetic mice compared with sedentary non-diabetic mice. We found that exercise reduced the loss of endothelial function from diabetes in both sexes. The reproducibility of perfusion measurements was limited primarily by our ability to select the same region in muscle and to effectively filter tissue clutter. CONCLUSIONS/SIGNIFICANCE: PD-US measurements can precisely follow site-specific changes in skeletal muscle perfusion related to diabetes over time, which fills the need for techniques capable of regularly monitoring atherosclerotic changes leading to ischemic vascular pathologies.

Deep Q-learning to globally optimize a k-D parameter search for medical imaging.

Background: Estimation of the global optima of multiple model parameters is valuable for precisely extracting parameters that characterize a physical environment. This is especially useful for imaging purposes, to form reliable, meaningful physical images with good reproducibility. However, it is challenging to avoid different local minima when the objective function is nonconvex. The problem of global searching of multiple parameters was formulated to be a k-D move in the parameter space and the parameter updating scheme was converted to be a state-action decision-making problem. Methods: We proposed a novel Deep Q-learning of Model Parameters (DQMP) method for global optimization which updated the parameter configurations through actions that maximized the Q-value and employed a Deep Reward Network (DRN) designed to learn global reward values from both visible fitting errors and hidden parameter errors. The DRN was constructed with Long Short-Term Memory (LSTM) layers followed by fully connected layers and a rectified linear unit (ReLU) nonlinearity. The depth of the DRN depended on the number of parameters. Through DQMP, the k-D parameter search in each step resembled the decision-making of action selections from 3k configurations in a k-D board game. Results: The DQMP method was evaluated by widely used general functions that can express a variety of experimental data and further validated on imaging applications. The convergence of the proposed DRN was evaluated, which showed that the loss values of six general functions all converged after 12 epochs. The parameters estimated by the DQMP method had relative errors of less than 4% for all cases, whereas the relative errors achieved by Q-learning (QL) and the Least Squares Method (LSM) were 17% and 21%, respectively. Furthermore, the imaging experiments demonstrated that the imaging of the parameters estimated by the proposed DQMP method were the closest to the ground truth simulation images when compared to other methods. Conclusions: The proposed DQMP method was able to achieve global optima, thus yielding accurate model parameter estimates. DQMP is promising for estimating multiple high-dimensional parameters and can be generalized to global optimization for many other complex nonconvex functions and imaging of physical parameters.

Monitoring Muscle Perfusion in Rodents During Short-Term Ischemia Using Power Doppler Ultrasound.

OBJECTIVE: The aim of this work was to evaluate the reliability of power Doppler ultrasound (PD-US) measurements made without contrast enhancement to monitor temporal changes in peripheral blood perfusion. METHODS: On the basis of pre-clinical rodent studies, we found that combinations of spatial registration and clutter filtering techniques applied to PD-US signals reproducibly tracked blood perfusion in skeletal muscle. Perfusion is monitored while modulating hindlimb blood flow. First, in invasive studies, PD-US measurements in deep muscle with laser speckle contrast imaging (LSCI) of superficial tissues made before, during and after short-term arterial clamping were compared. Then, in non-invasive studies, a pressure cuff was employed to generate longer-duration hindlimb ischemia. Here, B-mode imaging was also applied to measure flow-mediated dilation of the femoral artery while, simultaneously, PD-US was used to monitor downstream muscle perfusion to quantify reactive hyperemia. Measurements in adult male and female mice and rats, some with exercise conditioning, were included to explore biological variables. RESULTS: PD-US methods are validated through comparisons with LSCI measurements. As expected, no significant differences were found between sexes or fitness levels in flow-mediated dilation or reactive hyperemia estimates, although post-ischemic perfusion was enhanced with exercise conditioning, suggesting there could be differences between the hyperemic responses of conduit and resistive vessels. CONCLUSION: Overall, we found non-contrast PD-US imaging can reliably monitor relative spatiotemporal changes in muscle perfusion. This study supports the development of PD-US methods for monitoring perfusion changes in patients at risk for peripheral artery disease.

A data-driven approach to characterizing nonlinear elastic behavior of soft materials.

The Autoprogressive (AutoP) method is a data-driven inverse method that leverages finite element analysis (FEA) and machine learning (ML) techniques to build constitutive relationships from measured force and displacement data. Previous applications of AutoP in tissue-like media have focused on linear elastic mechanical behavior as the target object is infinitesimally compressed. In this study, we extended the application of AutoP in characterizing nonlinear elastic mechanical behavior as the target object undergoes finite compressive deformation. Guided by the prior of nonlinear media, we modified the training data generated by AutoP to speed its ability to learn to model deformations. AutoP training was validated using both synthetic and experimental data recorded from 3D objects. Force-displacement measurements were obtained using ultrasonic imaging from heterogeneous agar-gelatin phantoms. Measurement on samples of phantom components were analyzed to obtain independent measurements of material properties. Comparisons validated the material properties found from neural network constitutive models (NNCMs) trained using AutoP. Results were found to be robust to measurement errors and spatial variations in material properties.

Combining Spatial Registration With Clutter Filtering for Power-Doppler Imaging in Peripheral Perfusion Applications.

Power-Doppler ultrasonic (PD-US) imaging is sensitive to echoes from blood cell motion in the microvasculature but generally nonspecific because of difficulties with filtering nonblood-echo sources. We are studying the potential for using PD-US imaging for routine assessments of peripheral blood perfusion without contrast media. The strategy developed is based on an experimentally verified computational model of tissue perfusion that simulates typical in vivo conditions. The model considers directed and diffuse blood perfusion states in a field of moving clutter and noise. A spatial registration method is applied to minimize tissue motion prior to clutter and noise filtering. The results show that in-plane clutter motion is effectively minimized. While out-of-plane motion remains a strong source of clutter-filter leakage, those registration errors are readily minimized by straightforward modification of scanning techniques and spatial averaging.

Physics-guided machine learning for 3-D quantitative quasi-static elasticity imaging.

We present a 3D extension of the Autoprogressive Method (AutoP) for quantitative quasi-static ultrasonic elastography (QUSE) based on sparse sampling of force-displacement measurements. Compared to current model-based inverse methods, our approach requires neither geometric nor constitutive model assumptions. We build upon our previous report for 2D QUSE and demonstrate the feasibility of recovering the 3D linear-elastic material property distribution of gelatin phantoms under compressive loads. Measurements of boundary geometry, applied surface forces, and axial displacements enter into AutoP where a Cartesian neural network constitutive model (CaNNCM) interacts with finite element analyses to learn physically consistent material properties with no prior constitutive model assumption. We introduce a new regularization term uniquely suited to AutoP that improves the ability of CaNNCMs to extract information about spatial stress distributions from measurement data. Results of our study demonstrate that acquiring multiple sets of force-displacement measurements by moving the US probe to different locations on the phantom surface not only provides AutoP with the necessary information for a CaNNCM to learn the 3D material property distribution, but may significantly improve the accuracy of the Young's modulus estimates. Furthermore, we investigate the trade-offs of decreasing the contact area between the US transducer and phantom surface in an effort to increase sensitivity to surface force variations without additional instrumentation. Each of these modifications improves the ability of CaNNCMs trained in AutoP to learn the spatial distribution of Young's modulus from force-displacement measurements.

Experimental Validation of Perfusion Imaging With HOSVD Clutter Filters.

Novel pulsed-Doppler methods for perfusion imaging are validated using dialysis cartridges as perfusion phantoms. Techniques that were demonstrated qualitatively at 24 MHz, in vivo, are here examined quantitatively at 5 and 12.5 MHz using phantoms with the blood-mimicking fluid flow within cellulose microfibers. One goal is to explore a variety of flow states to optimize measurement sensitivity and flow accuracy. The results show that 2-3-s echo acquisitions at roughly 10 frames/s yield the highest sensitivity to flows of 1-4 mL/min. A second goal is to examine methods for setting the parameters of higher order singular value decomposition (HOSVD) clutter filters. For stationary or moving clutter, the velocity of the blood-mimicking fluid in the microfibers is consistently estimated within measurement uncertainty (mean coefficient of variation = 0.26). Power Doppler signals were equivalent for stationary and moving clutter after clutter filtering, increasing approximately 3 dB/mL/min of blood-mimicking fluid flow for 0 ≤ q ≤ 4 mL/min. Comparisons between phantom and preclinical images show that peripheral perfusion imaging can be reliably achieved without contrast enhancement.

Fluidity and elasticity form a concise set of viscoelastic biomarkers for breast cancer diagnosis based on Kelvin-Voigt fractional derivative modeling.

Cancer progression involves biomechanical changes within transformed cells and the surrounding extracellular matrix (ECM). The viscoelastic features of fluidity and elasticity that are based on a novel Kelvin-Voigt fractional derivative (KVFD) model were found capable of discriminating normal, benign and malignant breast biopsy tissues on the cellular scale. The improved specificity of KVFD model parameters derives from greater accuracy of fitting the entire approaching force-indentation measurement curve ([Formula: see text] > 0.99) compared with traditional elastic models ([Formula: see text] < 0.86). Moreover, model parameters can be interpreted in terms of histopathological features. First, statistical comparisons reveal there are significant differences (p < 0.001) in elasticity E0, fluidity [Formula: see text], and viscosity [Formula: see text] among healthy, benign, and malignant groups. Malignant breast tissues show low-value, broad-distributions in E0 and with high fluidity [Formula: see text] as compared with healthy and benign tissues. Second, histograms of E0 and [Formula: see text] provide distinctive features by fitting to Gaussian mixture (GM) models. The histograms of E0 and [Formula: see text] are best fit by two kernels GM for malignant tissues, indicating that the cells are soft but with high fluidity and the ECM is stiff but with low fluidity. However, the data suggest one-kernel GM model for benign tissue and a patched uniform distribution for healthy tissue. Third, using fluidity [Formula: see text] as the test statistic, the area under the receiver operator characteristic curve (AUC) is 0.701 ± 0.012 (p < 0.0001) for control versus malignant and 0.706 ± 0.013 (p < 0.0001) for benign versus malignant group. Variations in tissue fluidity and elasticity offer a concise set of viscoelastic biomarkers that correlate well with histopathological features.

Multi-parameter Sub-Hertz Analysis of Viscoelasticity With a Quality Metric for Differentiation of Breast Masses.

We applied sub-Hertz analysis of viscoelasticity (SAVE) to differentiate breast masses in pre-biopsy patients. Tissue response during external ramp-and-hold stress was ultrasonically detected. Displacements were used to acquire tissue viscoelastic parameters. The fast instantaneous response and slow creep-like deformations were modeled as the response of a linear standard solid from which viscoelastic parameters were estimated. These parameters were used in a multi-variable classification framework to differentiate malignant from benign masses identified by pathology. When employing all viscoelasticity parameters, SAVE resulted in 71.43% accuracy in differentiating lesions. When combined with ultrasound features and lesion size, accuracy was 82.24%. Adding a quality metric based on uniaxial motion increased the accuracy to 81.25%. When all three were combined with SAVE, accuracy was 91.3%. These results confirm the utility of SAVE as a robust ultrasound-based diagnostic tool for non-invasive differentiation of breast masses when used as stand-alone biomarkers or in conjunction with ultrasonic features.

Imaging the Landmarks of Vascular Recovery.

Background: Peripheral arterial disease (PAD) is a major worldwide health concern. Since the late 1990s therapeutic angiogenesis has been investigated as an alternative to traditional PAD treatments. Although positive preclinical results abound in the literature, the outcomes of human clinical trials have been discouraging. Among the challenges the field has faced has been a lack of standardization of the timings and measures used to validate new treatment approaches. Methods: In order to study the spatiotemporal dynamics of both perfusion and neovascularization in mice subjected to surgically-induced hindlimb ischemia (n= 30), we employed three label-free imaging modalities (a novel high-sensitivity ultrasonic Power Doppler methodology, laser speckle contrast, and photoacoustic imaging), as well as a tandem of radio-labeled molecular probes, 99mTc-NC100692 and 99mTc-BRU-5921 respectively, designed to detect two key modulators of angiogenic activity, αVß3 and HIF-1α , via scintigraphic imaging. Results: The multimodal imaging strategy reveals a set of "landmarks"-key physiological and molecular events in the healing process-that can serve as a standardized framework for describing the impact of emerging PAD treatments. These landmarks span the entire process of neovascularization, beginning with the rapid decreases in perfusion and oxygenation associated with ligation surgery, extending through pro-angiogenic changes in gene expression driven by the master regulator HIF-1α , and ultimately leading to complete functional revascularization of the affected tissues. Conclusions: This study represents an important step in the development of multimodal non-invasive imaging strategies for vascular research; the combined results offer more insight than can be gleaned through any of the individual imaging methods alone. Researchers adopting similar imaging strategies and will be better able to describe changes in the onset, duration, and strength of each of the landmarks of vascular recovery, yielding greater biological insight, and enabling more comprehensive cross-study comparisons. Perhaps most important, this study paves the road for more efficient translation of PAD research; emerging experimental treatments can be more effectively assessed and refined at the preclinical stage, ultimately leading to better next-generation therapies.

Acoustomotive optical coherence elastography for measuring material mechanical properties.

Acoustomotive optical coherence elastography (AM-OCE), a dynamic and internal excitation optical coherence elastography technique, is reported. Acoustic radiation force was used for internal mechanical excitation, and spectral-domain optical coherence tomography was used for detection. Mechanical properties of gelatin tissue phantoms were measured by AM-OCE and verified using rheometry results. Measured mechanical properties including shear moduli and shear damping parameters of the gelatin samples double when their polymer concentration increases from 3% to 4%. Spectral analysis was also performed on the acquired data, which improved the processing speed by a factor of 5 compared with a least-square fitting approach. Quantitative measurement, microscale resolution, and remote excitation are the main features of AM-OCE, which make the technique promising for measuring biomechanical properties.

Material properties from acoustic radiation force step response.

An ultrasonic technique for estimating viscoelastic properties of hydrogels, including engineered biological tissues, is being developed. An acoustic radiation force is applied to deform the gel locally while Doppler pulses track the induced movement. The system efficiently couples radiation force to the medium through an embedded scattering sphere. A single-element, spherically-focused, circular piston element transmits a continuous-wave burst to suddenly apply and remove a radiation force to the sphere. Simultaneously, a linear array and spectral Doppler technique are applied to track the position of the sphere over time. The complex shear modulus of the gel was estimated by applying a harmonic oscillator model to measurements of time-varying sphere displacement. Assuming that the stress-strain response of the surrounding gel is linear, this model yields an impulse response function for the gel system that may be used to estimate material properties for other load functions. The method is designed to explore the force-frequency landscape of cell-matrix viscoelasticity. Reported measurements of the shear modulus of gelatin gels at two concentrations are in close agreement with independent rheometer measurements of the same gels. Accurate modulus measurements require that the rate of Doppler-pulse transmission be matched to a priori estimates of gel properties.

Data-Driven Elasticity Imaging Using Cartesian Neural Network Constitutive Models and the Autoprogressive Method.

Quasi-static elasticity imaging techniques rely on model-based mathematical inverse methods to estimate mechanical parameters from force-displacement measurements. These techniques introduce simplifying assumptions that preclude exploration of unknown mechanical properties with potential diagnostic value. We previously reported a data-driven approach to elasticity imaging using artificial neural networks (NNs) that circumvents limitations associated with model-based inverse methods. NN constitutive models can learn stress-strain behavior from force-displacement measurements using the autoprogressive (AutoP) method without prior assumptions of the underlying constitutive model. However, information about internal structure was required. We invented Cartesian NN constitutive models (CaNNCMs) that learn the spatial variations of material properties. We are presenting the first implementation of CaNNCMs trained with AutoP to develop data-driven models of 2-D linear-elastic materials. Both simulated and experimental force-displacement data were used as input to AutoP to show that CaNNCMs are able to model both continuous and discrete material property distributions with no prior information of internal object structure. Furthermore, we demonstrate that CaNNCMs are robust to measurement noise and can reconstruct reasonably accurate Young's modulus images from a sparse sampling of measurement data. CaNNCMs are an important step toward clinical use of data-driven elasticity imaging using AutoP.

Role of ultrasonic shear rate estimation errors in assessing inflammatory response and vascular risk.

Atherosclerotic lesions preferentially originate in arterial regions that experience low wall shear stress (WSS) and reversing flow patterns. Therefore, routinely monitoring arterial WSS may help to identify the potential sites of early atherosclerosis. A new noninvasive ultrasonic method implemented with coded excitation techniques was utilized to improve WSS estimation accuracy and precision by providing high spatial and temporal resolution. WSS measurement errors were quantified in a model system by scanning a linearly varying WSS field (0.3 to 1.9 Pa) within a flow chamber. A 13-bit optimal code (Opt) was found to be most effective in reducing bias and standard deviation in WSS estimates down to approximately 10% and approximately 8%. The measurement errors slowly increased with input WSS for all imaging pulses. The expression of endothelial cellular adhesion molecules vascular cell adhesion molecule-1 (VCAM-1) and endothelial-leukocyte adhesion molecule-1 (E-selectin) was investigated over a similar shear range (0 to 1.6 Pa) to study the impact of relating shear-mediated cellular adhesion molecule (CAM) expression to inaccuracies in WSS measurements. We quantified this influence as the prediction error, which accounts for the ultrasonic measurement errors and the sensitivity of CAM expression within certain shear ranges. The highest prediction errors were observed at WSS <0.8 Pa, where CAM expression is most responsive to WSS. The results emphasize the importance of minimizing estimation errors, especially within low shear regions. Preliminary two-dimensional in vivo shear imaging is also presented to provide information about the spatial heterogeneity in arterial WSS distribution.