Low Variance Estimation of Backscatter Quantitative Ultrasound Parameters Using Dynamic Programming

Spatially Weighted Fidelity and Regularization Terms for Attenuation Imaging

Quantitative ultrasound (QUS) holds promise in enhancing diagnostic accuracy. For attenuation imaging, the regularized spectral log difference (RSLD) can generate accurate local attenuation maps. However, the performance of the method degrades when significant changes in backscatter amplitude occur. Variations in the technique were introduced involving a weighted approach to backscatter regularization, which, however, is not effective when changes in both attenuation and backscatter are present. This study introduces a novel approach that incorporates an L1-norm for backscatter regularization and spatially varying weights for both fidelity and regularization terms. The weights are calculated from an initial estimation of backscatter changes. Comparative analyses with simulated, phantom, and clinical data were performed. When changes in backscatter and attenuation occur, the proposed approach reduced the lowest root mean square error by up to 73%. It also improved the contrast-to-noise ratio (CNR) by a factor of 4.4 on average compared with previously available methods, considering the simulated and phantom data. In vivo results from healthy livers, thyroid nodules, and a breast tumor further confirm its effectiveness. In the liver, it is shown to be effective at reducing artifacts of attenuation images. In thyroid and breast tumors, the method demonstrated an enhanced CNR and better consistency of the attenuation measurements with the posterior acoustic enhancement. Overall, this approach offers promise for enhancing ultrasound attenuation imaging by helping differentiate tissue characteristics that may indicate pathology.

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.

Deep Estimation of Speckle Statistics Parametric Images

Quantitative Ultrasound (QUS) provides important information about the tissue properties. QUS parametric image can be formed by dividing the envelope data into small overlapping patches and computing different speckle statistics such as parameters of the Nakagami and Homodyned K-distributions (HK-distribution). The calculated QUS parametric images can be erroneous since only a few independent samples are available inside the patches. Another challenge is that the envelope samples inside the patch are assumed to come from the same distribution, an assumption that is often violated given that the tissue is usually not homogenous. In this paper, we propose a method based on Convolutional Neural Networks (CNN) to estimate QUS parametric images without patching. We construct a large dataset sampled from the HK-distribution, having regions with random shapes and QUS parameter values. We then use a well-known network to estimate QUS parameters in a multi-task learning fashion. Our results confirm that the proposed method is able to reduce errors and improve border definition in QUS parametric images.

Reverberation clutter signal suppression in ultrasound attenuation estimation using wavelet-based robust principal component analysis

Objective. Ultrasound attenuation coefficient estimation (ACE) has diagnostic potential for clinical applications such as quantifying fat content in the liver. Previously, we have proposed a system-independent ACE technique based on spectral normalization of different frequencies, called the reference frequency method (RFM). This technique does not require a well-calibrated reference phantom for normalization. However, this method may be vulnerable to severe reverberation clutter introduced by the body wall. The clutter superimposed on liver echoes may bias the estimation. Approach. We proposed to use robust principal component analysis, combined with wavelet-based sparsity promotion, to suppress the severe reverberation clutters. The capability to mitigate the reverberation clutters was validated through phantom and in vivo studies. Main Results. In the phantom studies with added reverberation clutters, higher normalized cross-correlation and smaller mean absolute errors were attained as compared to RFM results without the proposed method, demonstrating the capability to reconstruct tissue signals from reverberations. In a pilot patient study, the correlation between ACE and proton density fat fraction (PDFF), a measurement of liver fat by MRI as a reference standard, was investigated. The proposed method showed an improvement of the correlation (coefficient of determination, R = 0.82) as compared with the counterpart without the proposed method (R = 0.69). Significance: The proposed method showed the feasibility of suppressing the reverberation clutters, providing an important basis for the development of a robust ACE with large reverberation clutters.

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.

Identifying and overcoming limitations with in situ calibration beads for quantitative ultrasound

Ensuring the consistency of spectral-based quantitative ultrasound estimates in vivo necessitates accounting for diffraction, system effects, and propagation losses encountered in the tissue. Accounting for diffraction and system effects is typically achieved through planar reflector or reference phantom methods; however, neither of these is able to account for the tissue losses present in vivo between the ultrasound probe and the region of interest. In previous work, the feasibility of small titanium beads as in situ calibration targets (0.5-2 mm in diameter) was investigated. In this study, the importance of bead size for the calibration signal, the role of multiple echoes coming from the calibration bead, and sampling of the bead signal laterally through beam translation were examined. This work demonstrates that although the titanium beads naturally produce multiple reverberant echoes, time-windowing of the first echo provides the smoothest calibration spectrum for backscatter coefficient calculation. When translating the beam across the bead, the amplitude of the echo decreases rapidly as the beam moves across and past the bead. Therefore, to obtain consistent calibration signals from the bead, lateral interpolation is needed to approximate signals coming from the center of the bead with respect to the beam.

A multiparametric volumetric quantitative ultrasound imaging technique for soft tissue characterization

Quantitative ultrasound (QUS) offers a non-invasive and objective way to quantify tissue health. We recently presented a spatially adaptive regularization method for reconstruction of a single QUS parameter, limited to a two dimensional region. That proof-of-concept study showed that regularization using homogeneity prior improves the fundamental precision-resolution trade-off in QUS estimation. Based on the weighted regularization scheme, we now present a multiparametric 3D weighted QUS (3D QUS) method, involving the reconstruction of three QUS parameters: attenuation coefficient estimate (ACE), integrated backscatter coefficient (IBC) and effective scatterer diameter (ESD). With the phantom studies, we demonstrate that our proposed method accurately reconstructs QUS parameters, resulting in high reconstruction contrast and therefore improved diagnostic utility. Additionally, the proposed method offers the ability to analyze the spatial distribution of QUS parameters in 3D, which allows for superior tissue characterization. We apply a three-dimensional total variation regularization method for the volumetric QUS reconstruction. The 3D regularization involving N planes results in a high QUS estimation precision, with an improvement of standard deviation over the theoretical 1/N rate achievable by compounding N independent realizations. In the in vivo liver study, we demonstrate the advantage of adopting a multiparametric approach over the single parametric counterpart, where a simple quadratic discriminant classifier using feature combination of three QUS parameters was able to attain a perfect classification performance to distinguish between normal and fatty liver cases.

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.

Fast linear least-squares method for ultrasound attenuation and backscatter estimation

The ultrasonic attenuation and backscatter coefficients of tissues are relevant acoustic parameters due to their wide range of clinical applications. In this paper, a linear least-squares method for the estimation of these coefficients in a homogeneous region of interest based on pulse-echo measurements is proposed. The method efficiently fits an ultrasound backscattered signal model to the measurements in both the frequency and depth dimension simultaneously at a low computational cost. It is demonstrated that the inclusion of depth information has a positive effect particularly on the accuracy of the estimated attenuation. The sensitivity of the attenuation and backscatter coefficients' estimates to several predefined parameters such as the window length, window overlap and usable bandwidth of the spectrum is also studied. Comparison of the proposed method with a benchmark approach based on dynamic programming highlights better performance of our method in estimating these coefficients, both in terms of accuracy and computation time. Further analysis of the computation time as a function of the predefined parameters indicates our method's potential to be used in real-time clinical settings.

Ultrasound Scatterer Density Classification Using Convolutional Neural Networks and Patch Statistics

Quantitative ultrasound (QUS) can reveal crucial information on tissue properties, such as scatterer density. If the scatterer density per resolution cell is above or below 10, the tissue is considered as fully developed speckle (FDS) or underdeveloped speckle (UDS), respectively. Conventionally, the scatterer density has been classified using estimated statistical parameters of the amplitude of backscattered echoes. However, if the patch size is small, the estimation is not accurate. These parameters are also highly dependent on imaging settings. In this article, we adapt convolutional neural network (CNN) architectures for QUS and train them using simulation data. We further improve the network's performance by utilizing patch statistics as additional input channels. Inspired by deep supervision and multitask learning, we propose a second method to exploit patch statistics. We evaluate the networks using simulation data and experimental phantoms. We also compare our proposed methods with different classic and deep learning models and demonstrate their superior performance in the classification of tissues with different scatterer density values. The results also show that we are able to classify scatterer density in different imaging parameters with no need for a reference phantom. This work demonstrates the potential of CNNs in classifying scatterer density in ultrasound images.

Total attenuation compensation for backscatter coefficient estimation using full angular spatial compounding

The backscatter coefficient (BSC) quantifies the frequency-dependent reflectivity of tissues. Accurate estimation of the BSC is only possible with the knowledge of the attenuation coefficient slope (ACS) of the tissues under examination. In this study, the use of attenuation maps constructed using full angular spatial compounding (FASC) is proposed for attenuation compensation when imaging integrated BSCs. Experimental validation of the proposed approach was obtained using two cylindrical physical phantoms with off-centered inclusions having different ACS and BSC values than the background, and in a phantom containing an ex vivo chicken breast sample embedded in an agar matrix. With the phantom data, three different ACS maps were employed for attenuation compensation: (1) a ground truth ACS map constructed using insertion loss techniques, (2) the estimated ACS map using FASC attenuation imaging, and (3) a uniform ACS map with a value of 0.5 dBcm^{\protect \relax \special {t4ht=-}1}MHz^{\protect \relax \special {t4ht=-}1}, which is commonly used to represent attenuation in soft tissues. Comparable results were obtained when using the ground truth and FASC-estimated ACS maps in term of inclusion detectability and estimation accuracy, with averaged fractional error below 2.8 dB in both phantoms. Conversely, the use of the homogeneous ACS map resulted in higher levels of fractional error (>10 dB), which demonstrates the importance of an accurate attenuation compensation. The results with the ex vivo tissue sample were consistent with the observations using the physical phantoms, with the FASC-derived ACS map providing comparable BSC images to those formed using the ground truth ACS map and more accurate than those BSC images formed using a uniform ACS. These results suggest that BSCs can be reliably estimated using FASC when a self-consistent attenuation compensation stemming from prior estimation of an accurate ACS map is used.

Analytic Global Regularized Backscatter Quantitative Ultrasound

Although a variety of techniques have been developed to reduce the appearance of B-mode speckle, quantitative ultrasound (QUS) aims at extracting the hidden properties of the tissue. Herein, we propose two novel techniques to accurately and precisely estimate two important QUS parameters, namely, the average attenuation coefficient and the backscatter coefficient. Both the techniques optimize a cost function that incorporates data and continuity constraint terms, which we call AnaLytical Global rEgularized BackscatteR quAntitative ultrasound (ALGEBRA). We propose two versions of ALGEBRA, namely, 1-D- and 2-D-ALGEBRA. In 1-D-ALGEBRA, the regularized cost function is formulated in the axial direction, and the QUS parameters are calculated for one line of radio frequency (RF) echo data. In 2-D-ALGEBRA, the regularized cost function is formulated for the entire image, and the QUS parameters throughout the image are estimated simultaneously. This simultaneous optimization allows 2-D-ALGEBRA to "see" all the data before estimating the QUS parameters. In both the methods, we efficiently optimize the cost functions by casting it as a sparse linear system of equations. As a result of this efficient optimization, 1-D-ALGEBRA and 2-D-ALGEBRA are, respectively, 600 and 300 times faster than optimization using the dynamic programming (DP) method previously proposed by our group. In addition, the proposed technique has fewer input parameters that require manual tuning. Our results demonstrate that the proposed ALGEBRA methods substantially outperform least-square and DP methods in estimating the QUS parameters in phantom experiments.

Frequency-dependent attenuation reconstruction with an acoustic reflector

Attenuation of ultrasound waves varies with tissue composition, hence its estimation offers great potential for tissue characterization and diagnosis and staging of pathology. We recently proposed a method that allows to spatially reconstruct the distribution of the overall ultrasound attenuation in tissue based on computed tomography, using reflections from a passive acoustic reflector. This requires a standard ultrasound transducer operating in pulse-echo mode and a calibration protocol using water measurements, thus it can be implemented on conventional ultrasound systems with minor adaptations. Herein, we extend this method by additionally estimating and imaging the frequency-dependent nature of local ultrasound attenuation for the first time. Spatial distributions of attenuation coefficient and exponent are reconstructed, enabling an elaborate and expressive tissue-specific characterization. With simulations, we demonstrate that our proposed method yields a low reconstruction error of 0.04 dB/cm at 1 MHz for attenuation coefficient and 0.08 for the frequency exponent. With tissue-mimicking phantoms and ex-vivo bovine muscle samples, a high reconstruction contrast as well as reproducibility are demonstrated. Attenuation exponents of a gelatin-cellulose mixture and an ex-vivo bovine muscle sample were found to be, respectively, 1.4 and 0.5 on average, consistently from different images of their heterogeneous compositions. Such frequency-dependent parametrization could enable novel imaging and diagnostic techniques, as well as facilitate attenuation compensation of other ultrasound-based imaging techniques.

Segmentation of Ultrasound Images based on Scatterer Density using U-Net

Quantitative ultrasound can provide an objective estimation of different tissue properties, which may be used for tissue characterization and detection of abnormal tissue. The effective number of scatterers in different parts of a tissue is one of the important tissue properties that can be estimated by quantitative ultrasound techniques. The envelope echo is the signal which is usually used to estimate the scatterer density. In this study, we proposed using deep learning to estimate the effective number of scatterers. We generated 2000 simulated phantom data containing randomly distributed inclusions with three different values for number of scatterers per resolution cell. We used U-Net to segment the envelope data and to distinguish three different values of scatterer densities. We show that U-Net can discriminate different scattering regimes, particularly, when the difference between the number of scatterers is substantial. The overall accuracy of the network is 83.9%, and the average sensitivity and specificity among the three classes are 83.1% and 92.3% respectively. This study confirms the potential of deep learning framework in quantitative ultrasound and estimation of tissue properties using ultrasound images.

A Pilot Study on Scatterer Density Classification of Ultrasound Images Using Deep Neural Networks

Quantitative ultrasound estimates different intrinsic tissue properties, which can be used for tissue characterization. Among different tissue properties, the effective number of scatterers per resolution cell is an important parameter, which can be estimated by the echo envelope. Assuming the signal is stationary and coherent, if the number of scatterers per resolution cell is above approximately 10, envelope signal is considered to be fully developed speckle (FDS) and otherwise they are from low scatterer number density (LSND). Two statistical parameters named R and S are often calculated from envelope intensity to classify FDS from LSND. The main problem is that limited data from small patches often renders this classification inaccurate. Herein, we propose two techniques based on neural networks to estimate the effective number of scatterers. The first network is a multi-layer perceptron (MLP) that uses the hand-crafted features of R and S for classification. The second network is a convolutional neural network (CNN) that does not need hand-crafted features and instead utilizes spectrum and the envelope intensity directly. We show that the proposed MLP works very well for large patches wherein a reliable estimation of R and S can be made. However, its classification becomes inaccurate for small patches, where the proposed CNN provides accurate classifications.

Optimization of Ultrasound Backscatter Spectroscopy to Assess Neurotoxic Effects of Anesthesia in the Newborn Non-human Primate Brain

Studies in animal models have revealed that long exposures to anesthetics can induce apoptosis in the newborn and young developing brain. These effects have not been confirmed in humans because of the lack of a non-invasive, practical in vivo imaging tool with the ability to detect these changes. Following the successful use of ultrasound backscatter spectroscopy (UBS) to monitor in vivo cell death in breast tumors, we aimed to use UBS to assess the neurotoxicity of the anesthetic sevoflurane (SEVO) in a non-human primate (NHP) model. Sixteen 2- to 7-day-old rhesus macaques were exposed for 5 h to SEVO. Ultrasound scanning was done with a phased array transducer on a clinical ultrasound scanner operated at 10 MHz. Data consisting of 10-15 frames of radiofrequency (RF) echo signals from coronal views of the thalamus were obtained 0.5 and 6.0 h after initiating exposure. The UBS parameter "effective scatterer size" (ESS) was estimated by fitting a scattering form factor (FF) model to the FF measured from RF echo signals. The approach involved analyzing the frequency dependence of the measured FF to characterize scattering sources and selecting the FF model based on a χ² goodness-of-fit criterion. To assess data quality, a rigorous acceptance criterion based on the analysis of prevalence of diffuse scattering (an assumption in the estimation of ESS) was established. ESS changes after exposure to SEVO were compared with changes in a control group of five primates for which ultrasound data were acquired at 0 and 10 min (no apoptosis expected). Over the entire data set, the average measured FF at 0.5 and 6.0 h monotonically decreased with frequency, justifying fitting a single FF over the analysis bandwidth. χ² values of a (inhomogeneous continuum) Gaussian FF model were one-fifth those of the discrete fluid sphere model, suggesting that a continuum scatterer model better represents ultrasound scattering in the young rhesus brain. After application of the data quality criterion, only 5 of 16 subjects from the apoptotic group and 5 of 5 subjects from the control group fulfilled the acceptance criteria. All subjects in the apoptotic group that passed the acceptance criterion exhibited a significant ESS reduction at 6.0 h. These changes (-6.4%, 95% Interquartile Range: -14.3% to -3.3%) were larger than those in the control group (-0.8%, 95% Interquartile Range: -2.0% to 1.5%]). Data with a low prevalence of diffuse scattering corresponded to possibly biased results. Thus, ESS has the potential to detect changes in brain microstructure related to anesthesia-induced apoptosis.

Regularized Estimation of Effective Scatterer Size and Acoustic Concentration Quantitative Ultrasound Parameters Using Dynamic Programming

The objective of quantitative ultrasound (QUS) is to characterize tissue microstructure by parametrizing backscattered radiofrequency (RF) signals from clinical ultrasound scanners. Herein, we develop a novel technique based on dynamic programming (DP) to simultaneously estimate the acoustic attenuation, the effective scatterer size (ESS), and the acoustic concentration (AC) from ultrasound backscattered power spectra. This is achieved through two different approaches: (1) using a Gaussian form factor (GFF) and (2) using a general form factor (gFF) that is more flexible than the Gaussian form factor but involves estimating more parameters. Both DP methods are compared to an adaptation of a previously proposed least-squares (LSQ) method. Simulation results show that in the GFF approach, the variance of DP is on average 88%, 75% and 32% lower than that of LSQ for the three estimated QUS parameters. The gFF approach also yields similar improvements.

In vivo Ultrafast Quantitative Ultrasound and Shear Wave Elastography Imaging on Farm-Raised Duck Livers during Force Feeding

Shear wave elastography (speed and dispersion), local attenuation coefficient slope and homodyned-K parametric imaging were used for liver steatosis grading. These ultrasound biomarkers rely on physical interactions between shear and compression waves with tissues at both macroscopic and microscopic scales. These techniques were applied in a context not yet studied with ultrasound imaging, that is, monitoring steatosis of force-fed duck livers from pre-force-fed to foie gras stages. Each estimated feature presented a statistically significant trend along the feeding process (p values <10^-3). However, whereas a monotonic increase in the shear wave speed was observed along the process, most quantitative ultrasound features exhibited an absolute maximum value halfway through the process. As the liver fat fraction in foie gras is much higher than that seen clinically, we hypothesized that a change in the ultrasound scattering regime is encountered for high-fat fractions, and consequently, care has to be taken when applying ultrasound biomarkers to grading of severe states of steatosis.

Estimation of Backscatter Coefficients Using an In Situ Calibration Source

The objective of this article is to demonstrate the feasibility of estimating the backscatter coefficient (BSC) using an in situ calibration source. Traditional methods of estimating the BSC in vivo using a reference phantom technique do not account for transmission losses due to intervening layers between the ultrasonic source and the tissue region to be interrogated, leading to increases in bias and variance of BSC-based estimates. To account for transmission losses, an in situ calibration approach is proposed. The in situ calibration technique employs a titanium sphere that is well-characterized ultrasonically, biocompatible, and embedded inside the sample. A set of experiments was conducted to evaluate the embedded titanium spheres as in situ calibration targets for BSC estimation. The first experiment quantified the backscattered signal strength from titanium spheres of three sizes: 0.5, 1, and 2 mm in diameter. The second set of experiments assessed the repeatability of BSC estimates from the titanium spheres and compared these BSCs to theory. The third set of experiments quantified the ability of the titanium bead to provide an in situ reference spectrum in the presence of a lossy layer on top of the sample. The final set of experiments quantified the ability of the bead to provide a calibration spectrum over multiple depths in the sample. All experiments were conducted using an L9-4/38 linear array connected to a SonixOne system. The strongest signal was observed from the 2-mm titanium bead with the signal-to-noise ratio (SNR) of 11.6 dB with respect to the background speckle. Using an analysis bandwidth of 2.5-5.5 MHz, the mean differences between the experimentally derived BSCs and BSCs derived from the Faran theory were 0.54 and 0.76 dB using the array and a single-element transducer, respectively. The BSCs estimated using the in situ calibration approach without the layer and with the layer and using the reference phantom approach with the layer were compared to the reference phantom approach without the layer present. The mean differences in BSCs were 0.15, 0.73, and -9.69 dB, respectively. The mean differences of the BSCs calculated from data blocks located at depths that were either 30 pulse lengths above or below the actual bead depth compared to the BSC calculated at bead depth were -1.55 and -1.48 dB, respectively. The results indicate that an in situ calibration target can account for overlaying tissue losses, thereby improving the robustness of BSC-based estimates.

Validation of differences in backscatter coefficients among four ultrasound scanners with different beamforming methods

PURPOSE

The backscatter coefficient (BSC) indicates the absolute scatterer property of a material, independently of clinicians and system settings. Our study verified that the BSC differed among the scanners, transducers, and beamforming methods used for quantitative ultrasound analyses of biological tissues.

METHODS

Measurements were performed on four tissue-mimicking homogeneous phantoms containing spherical scatterers with mean diameters of 20 and 30 µm prepared at concentrations of 0.5 and 2.0 wt%, respectively. The BSCs in the different systems were compared using ultrasound scanners with two single-element transducers and five linear high- or low-frequency probes. The beamforming methods were line-by-line formation using focused imaging (FI) and parallel beam formation using plane wave imaging (PWI). The BSC of each system was calculated by the reference phantom method. The mean deviation from the theoretical BSC computed by the Faran model was analyzed as the benchmark validation of the calculated BSC.

RESULTS

The BSCs calculated in systems with different properties and beamforming methods well concurred with the theoretical BSC. The mean deviation was below ± 2.8 dB on average, and within the approximate standard deviation (± 2.2 dB at most) in all cases. These variations agreed with a previous study in which the largest error among four different scanners with FI beamforming was 3.5 dB.

CONCLUSION

The BSC in PWI was equivalent to those in the other systems and to those of FI beamforming. This result indicates the possibility of ultra-high frame-rate BSC analysis using PWI.

System-Independent Ultrasound Attenuation Coefficient Estimation Using Spectra Normalization

Ultrasound attenuation coefficient estimation (ACE) has diagnostic potential for clinical applications such as differentiating tumors and quantifying fat content in the liver. The two commonly used ACE methods in the ultrasound array imaging system are the spectral shift method and the reference-phantom-based methods. The spectra shift method estimates the central frequency downshift along depth, whereas the reference-phantom-based methods use a well-calibrated phantom to cancel system dependent effects in attenuation estimation. In this study, we propose a novel system-independent ACE technique based on spectra normalization of different frequencies. This technique does not require a reference phantom for normalization. The power of each frequency component is normalized by the power of an adjacent frequency component in the spectrum to cancel system-dependent effects, such as focusing and time gain compensation (TGC). This method is referred to as the reference frequency method (RFM), and its performance has been evaluated in phantoms and in vivo liver studies. The RFM technique can be applied to various transducer geometries (e.g., linear or curved arrays) with different beam patterns (e.g., focused or unfocused).