Analytic Global Regularized Backscatter Quantitative Ultrasound

Inter- and Intra-Operator Variability of Regularized Backscatter Quantitative Ultrasound for the Characterization of Breast Masses

OBJECTIVES

Here we report on the intra- and inter-operator variability of the backscatter coefficient (BSC) estimated with a new low-variance quantitative ultrasound (QUS) approach applied to breast lesions in vivo.

METHODS

Radiofrequency (RF) echo signals were acquired from 29 BIRADS 4 and 5 breast lesions in 2 sequential cohorts following 2 imaging protocols: cohort 1) radial and antiradial views, and cohort 2) short- and long-axis views. Protocol 2 was implemented after retraining and discussion on how to improve reproducibility. Each patient was scanned by at least 2 of 3 radiologists; each performed 3 acquisitions with transducer and patient repositioning in between acquisitions. BSC was estimated using a low-variance QUS approach based on regularization. Intra- and inter-operator variability of the intra-lesion median BSC was evaluated with a multifactorial ANOVA test (P-values) and the intraclass correlation coefficient (ICC).

RESULTS

Inter-operator variability was only significant in the first protocol (P < .007); ICCinter  = .77 (95% CI .71-.82), indicating good inter-operator agreement. In the second protocol, the inter-operator variability was not significant (P > .05) and agreement was excellent (ICCinter  = .92 [.89-.94]). In both protocols, the intra-operator variability was not significant.

CONCLUSIONS

Our findings demonstrate the need for standardizing image acquisition protocols for backscatter-based QUS to reduce inter-operator variability and ensure its successful translation to the characterization of suspicious breast masses.

Physics-Inspired Regularized Pulse-Echo Quantitative Ultrasound: Efficient Optimization With ADMM

Pulse-echo quantitative ultrasound (PEQUS), which estimates the quantitative properties of tissue microstructure, entails estimating the average attenuation and the backscatter coefficient (BSC). Growing recent research has focused on the regularized estimation of these parameters. Herein, we make two contributions to this field: first, we consider the physics of the average attenuation and backscattering to devise regularization terms accordingly. More specifically, since the average attenuation gradually alters in different parts of the tissue, while BSC can vary markedly from tissue to tissue, we apply L2 and L1 norms for the average attenuation and the BSC, respectively. Second, we multiply different frequencies and depths of the power spectra with different weights according to their noise levels. Our rationale is that the high-frequency contents of the power spectra at deep regions have a low signal-to-noise ratio (SNR). We exploit the alternating direction method of multipliers (ADMM) for optimizing the cost function. The qualitative and quantitative evaluations of bias and variance exhibit that our proposed algorithm improves the estimations of the average attenuation and the BSC up to about 100%.

Simultaneous imaging of ultrasonic relative backscatter and attenuation coefficients for quantitative liver steatosis assessment

Prevalence of liver disease is continuously increasing and nonalcoholic fatty liver disease (NAFLD) is the most common etiology. We present an approach to detect the progression of liver steatosis based on quantitative ultrasound (QUS) imaging. This study was performed on a group of 55 rats that were subjected to a control or methionine and choline deficient (MCD) diet known to induce NAFLD. Ultrasound (US) measurements were performed at 2 and 6 weeks. Thereafter, animals were humanely euthanized and livers excised for histological analysis. Relative backscatter and attenuation coefficients were simultaneously estimated from the US data and envelope signal-to-noise ratio was calculated to train a regression model for: (1) fat fraction percentage estimation and (2) performing classification according to Brunt's criteria in grades (0 <5%; 1, 5-33%; 2, >33-66%; 3, >66%) of liver steatosis. The trained regression model achieved an [Formula: see text] of 0.97 (p-value < 0.01) and a RMSE of 3.64. Moreover, the classification task reached an accuracy of 94.55%. Our results suggest that in vivo QUS is a promising noninvasive imaging modality for the early assessment of NAFLD.

High-Resolution Ultrasound Characterization of Local Scattering in Cancer Tissue

Ultrasound (US) has afforded an approach to tissue characterization for more than a decade. The challenge is to reveal hidden patterns in the US data that describe tissue function and pathology that cannot be seen in conventional US images. Our group has developed a high-resolution analysis technique for tissue characterization termed H-scan US, an imaging method used to interpret the relative size of acoustic scatterers. In the present study, the objective was to compare local H-scan US image intensity with registered histologic measurements made directly at the cellular level. Human breast cancer cells (MDA-MB 231, American Type Culture Collection, Manassas, VA, USA) were orthotopically implanted into female mice (N = 5). Tumors were allowed to grow for approximately 4 wk before the study started. In vivo imaging of tumor tissue was performed using a US system (Vantage 256, Verasonics Inc., Kirkland, WA, USA) equipped with a broadband capacitive micromachined ultrasonic linear array transducer (Kolo Medical, San Jose, CA, USA). A 15-MHz center frequency was used for plane wave imaging with five angles for spatial compounding. H-scan US image reconstruction involved use of parallel convolution filters to measure the relative strength of backscattered US signals. Color codes were applied to filter outputs to form the final H-scan US image display. For histologic processing, US imaging cross-sections were carefully marked on the tumor surface, and tumors were excised and sliced along the same plane. By use of optical microscopy, whole tumor tissue sections were scanned and digitized after nuclear staining. US images were interpolated to have the same number of pixels as the histology images and then spatially aligned. Each nucleus from the histologic sections was automatically segmented using custom MATLAB software (The MathWorks Inc., Natick, MA, USA). Nuclear size and spacing from the histology images were then compared with local H-scan US image features. Overall, local H-scan US image intensity exhibited a significant correlation with both cancer cell nuclear size (R2 > 0.27, p < 0.001) and the inverse relationship with nuclear spacing (R2 > 0.17, p < 0.001).

Recent Advances in Attenuation Estimation

This chapter reviews some of the recent advances in the estimation of the local and the total attenuation, with an emphasis on reducing the bias and variance of the estimates. A special focus is put on describing the effect of power spectrum estimation on bias and variance, the introduction of regularization strategies, as well as on eliminating the need to use reference phantoms for compensating for system dependent effects.

Regularized Ultrasound Phantom-Free Local Attenuation Coefficient Slope (ACS) Imaging in Homogeneous and Heterogeneous Tissues

Attenuation maps or measurements based on the local attenuation coefficient slope (ACS) in quantitative ultrasound (QUS) have shown potential for the diagnosis of liver steatosis. In liver cancers, tissue abnormalities and tumors detected using ACS are also of interest to provide new image contrast to clinicians. Current phantom-based approaches have the limitation of assuming a comparable speed of sound between the reference phantom and insonified tissues. Moreover, these methods present the inconvenience for operators to acquire data on phantoms and patients. The main goal was to alleviate these drawbacks by proposing a methodology for constructing phantom-free regularized (PF-R) local ACS maps and investigate the performance in both homogeneous and heterogeneous media. The proposed method was tested on two tissue-mimicking media with different ACS constructed as homogeneous phantoms, side-by-side and top-to-bottom phantoms, and inclusion phantoms with different attenuations. Moreover, an in vivo proof-of-concept was performed on healthy, steatotic, and cancerous human liver datasets. Modifications brought to previous works include: 1) a linear interpolation of the power spectrum in the log scale; 2) the relaxation of the underlying hypothesis on the diffraction factor; 3) a generalization to nonhomogeneous local ACS; and 4) an adaptive restriction of frequencies to a more reliable range than the usable frequency range. Regularization was formulated as a generalized least absolute shrinkage and selection operator (LASSO), and a variant of the Bayesian information criterion (BIC) was applied to estimate the Lagrangian multiplier on the LASSO constraint. In addition, we evaluated the proposed algorithm when applying median filtering before and after regularization. Tests conducted showed that the PF-R yielded robust results in all tested conditions, suggesting potential for additional validation as a diagnosis method.

Spatially Variant Ultrasound Attenuation Mapping Using a Regularized Linear Least-Squares Approach

Quantitative ultrasound methods aim to estimate the acoustic properties of the underlying medium, such as the attenuation and backscatter coefficients, and have applications in various areas including tissue characterization. In practice, tissue heterogeneity makes the coefficient estimation challenging. In this work, we propose a computationally efficient algorithm to map spatial variations of the attenuation coefficient. Our proposed approach adopts a fast, linear least-squares strategy to fit the signal model to data from pulse-echo measurements. As opposed to existing approaches, we directly estimate the attenuation map, that is, the local attenuation coefficient at each axial location by solving a joint estimation problem. In particular, we impose a physical model that couples all these local estimates and combine it with a smooth regularization to obtain a smooth map. Compared to the conventional spectral log difference method and the more recent ALGEBRA approach, we demonstrate that the attenuation estimates obtained by our method are more accurate and better correlate with the ground-truth attenuation profiles over a wide range of spatial and contrast resolutions.

Robust Scatterer Number Density Segmentation of Ultrasound Images

Quantitative ultrasound (QUS) aims to reveal information about the tissue microstructure using backscattered echo signals from clinical scanners. Among different QUS parameters, scatterer number density is an important property that can affect the estimation of other QUS parameters. Scatterer number density can be classified into high or low scatterer densities. If there are more than ten scatterers inside the resolution cell, the envelope data are considered as fully developed speckle (FDS) and, otherwise, as underdeveloped speckle (UDS). In conventional methods, the envelope data are divided into small overlapping windows (a strategy here we refer to as patching), and statistical parameters, such as SNR and skewness, are employed to classify each patch of envelope data. However, these parameters are system-dependent, meaning that their distribution can change by the imaging settings and patch size. Therefore, reference phantoms that have known scatterer number density are imaged with the same imaging settings to mitigate system dependency. In this article, we aim to segment regions of ultrasound data without any patching. A large dataset is generated, which has different shapes of scatterer number density and mean scatterer amplitude using a fast simulation method. We employ a convolutional neural network (CNN) for the segmentation task and investigate the effect of domain shift when the network is tested on different datasets with different imaging settings. Nakagami parametric image is employed for multitask learning to improve performance. Furthermore, inspired by the reference phantom methods in QUS, a domain adaptation stage is proposed, which requires only two frames of data from FDS and UDS classes. We evaluate our method for different experimental phantoms and in vivo data.

Quantitative ultrasound imaging of soft biological tissues: a primer for radiologists and medical physicists

Quantitative ultrasound (QUS) aims at quantifying interactions between ultrasound and biological tissues. QUS techniques extract fundamental physical properties of tissues based on interactions between ultrasound waves and tissue microstructure. These techniques provide quantitative information on sub-resolution properties that are not visible on grayscale (B-mode) imaging. Quantitative data may be represented either as a global measurement or as parametric maps overlaid on B-mode images. Recently, major ultrasound manufacturers have released speed of sound, attenuation, and backscatter packages for tissue characterization and imaging. Established and emerging clinical applications are currently limited and include liver fibrosis staging, liver steatosis grading, and breast cancer characterization. On the other hand, most biological tissues have been studied using experimental QUS methods, and quantitative datasets are available in the literature. This educational review addresses the general topic of biological soft tissue characterization using QUS, with a focus on disseminating technical concepts for clinicians and specialized QUS materials for medical physicists. Advanced but simplified technical descriptions are also provided in separate subsections identified as such. To understand QUS methods, this article reviews types of ultrasound waves, basic concepts of ultrasound wave propagation, ultrasound image formation, point spread function, constructive and destructive wave interferences, radiofrequency data processing, and a summary of different imaging modes. For each major QUS technique, topics include: concept, illustrations, clinical examples, pitfalls, and future directions.