Comparison of ultrasound attenuation and backscatter estimates in layered tissue-mimicking phantoms among three clinical scanners

Ultrasound backscatter coefficient for fat quantification is affected by the measurement depth

PURPOSE

It has been reported that the estimate of ultrasound attenuation coefficient (AC) is affected by depth of measurement, with linear decrease of values with depth. It is unknown whether backscatter coefficient (BSC) has similar behavior.

METHODS

This retrospective study was performed with Sequoia US system equipped with ultrasound derived fat fraction (UDFF) algorithm (Siemens Healthineers, Issaquah, WA, USA) that combines BSC with AC. UDFF was obtained positioning upper edge of the region of interest at 1.5,2,3,4,5 cm below liver capsule. BSC data were extracted from UDFF offline. A fractional polynomial regression, which selects the best model considering the polynomial development of the variables of interest, was used. Covariates included were age, sex, skin-to-liver-capsule distance, stiffness. Distance was included as linear factor or with a power ranging from - 2 to 3, and the best fitting model was chosen according to partial F test. Body mass index (BMI) was not included because of collinearity with skin-to-liver capsule distance.

RESULTS

104 individuals (56 females; age: 57.9 ± 13.0 years; BMI: 29.0 ± 6.5 kg/m2; skin-to-liver-capsule distance: 2.3 ± 0.7 cm; liver stiffness: 7.5 ± 5.5 kiloPascal) were studied. Best fitting model for BSC included a combination of depth as linear factor and square root. BSC showed a decrease of - 13.98 dB/cm-steradian for each logarithmic increase of 1 cm depth (coefficient: - 13.98; 95% CI: - 21.016; - 5.379; p = .001). Skin-to-liver-capsule distance and stiffness also were independent predictors of BSC.

CONCLUSIONS

The estimation of the BSC in the liver exhibits a depth dependence that significantly affects results. A standardized acquisition protocol is needed to compare results and reliably assess changes over time.

Confounders of Ultrasound Attenuation Imaging in a Linear Probe Using the Canon Aplio i800 System: A Phantom Study

There have been studies showing attenuation imaging (ATI) with ultrasound as an approach to diagnose liver diseases such as steatosis or cirrhosis. So far, this technique has only been used on a convex probe. The goal of the study was to investigate the feasibility of ATI measurements using the linear array on a canon Aplio i800 scanner on certified phantoms. Three certified liver tissue attenuation phantoms were measured in five different positions using a linear probe. The effects of positioning and depth were explored and compared. The values were compared to the certified expected value for each phantom as well as the different measurement values for each measurement position. The ATI measurements on phantoms showed significant effect for the different probe positions and region of interest (ROI) depths. Values taken in the center with the probe perpendicular to the phantom were closest to certified values. Median values at 2.5-4.5 cm depth for phantoms 1 and 2 and 0.5-2.5 cm for phantom 3 were comparable with certified values. Measurements taken at a depth greater than 6 cm in any position were the least representative of the certified values (p-value < 0.01) and had the widest range throughout the different sessions. ATI measurements can be performed with the linear probe in phantoms; however, careful consideration should be given to depth dependency, as it can significantly affect measurement values. Remaining measurements at various depths within the 0.5-6.0 cm range showed deviation from the certified values of approximately 25%.

Quantitative evaluation of rat sciatic nerve degeneration using high-frequency ultrasound

In this study, we evaluated the utility of using high-frequency ultrasound to non-invasively track the degenerative process in a rat model of peripheral nerve injury. Primary analyses explored spatial and temporal changes in quantitative backscatter coefficient (BSC) spectrum-based outcomes and B-mode textural outcomes, using gray level co-occurrence matrices (GLCMs), during the progressive transition from acute to chronic injury. As secondary analyses, correlations among GLCM and BSC spectrum-based parameters were evaluated, and immunohistochemistry were used to suggest a structural basis for ultrasound outcomes. Both mean BSC spectrum-based and mean GLCM-based measures exhibited significant spatial differences across presurgical and 1-month/2-month time points, distal stumps enclosed proximity to the injury site being particularly affected. The two sets of parameters sensitively detected peripheral nerve degeneration at 1-month and 2-month post-injury, with area under the receiver operating charactersitic curve > 0.8 for most parameters. The results also indicated that the many BSC spectrum-based and GLCM-based parameters significantly correlate with each other, and suggested a common structural basis for a diverse set of quantitative ultrasound parameters. The findings of this study suggest that BSC spectrum-based and GLCM-based analysis are promising non-invasive techniques for diagnosing peripheral nerve degeneration.

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.

Liver Fat Quantification With Ultrasound: Depth Dependence of Attenuation Coefficient

OBJECTIVES

The primary aim was to estimate the influence of various depths on ultrasound attenuation coefficient (AC) of multiple vendors in the liver. The secondary aim was to evaluate the impact of region of interest (ROI) size on AC measurements in a subset of participants.

METHODS

This Institutional Review Board (IRB)-approved Health Insurance Portability and Accountability Act (HIPAA)-compliant retrospective study was carried out in two centers using AC-Canon and AC-Philips algorithms and extracting AC-Siemens values from ultrasound-derived fat fraction algorithm. Measurements were performed positioning ROI upper edge (3 cm size) at 2, 3, 4, 5 cm from the liver capsule with AC-Canon and AC-Philips and at 1.5, 2, 3 cm with Siemens algorithm. In a subset of participants, measurements were obtained with 1 and 3 cm ROI size. Univariate and multivariate linear regression models and Lin's concordance correlation coefficient (CCC) were used for statistical analysis as appropriate.

RESULTS

Three different cohorts were studied. Sixty-three participants (34 females; mean age: 51 ± 14 years) were studied with AC-Canon, 60 (46 females; mean age: 57 ± 11 years) with AC-Philips, and 50 (25 females; 61 ± 13 years) with AC-Siemens. There was a decrease in AC values per 1 cm increase in depth in all. In multivariable analysis, the coefficient was -0.049 (-0.060; -0.038 P < .001) with AC-Canon, -0.058 (-0.066; -0.049 P < .001) with AC-Philips and -0.081 (-0.112; -0.050 P < .001) with AC-Siemens. AC values with 1 cm ROI were significantly higher than those obtained with 3 cm ROI at all depths (P < .001) but the agreement between AC values obtained with different ROI size was excellent (CCC 0.82 [0.77-0.88]).

CONCLUSIONS

There is depth dependence in AC measurement that affects results. A standardized protocol with fixed ROI's depth and size is needed.

US Quantification of Liver Fat: Past, Present, and Future

Fatty liver disease has a high and increasing prevalence worldwide, is associated with adverse cardiovascular events and higher long-term medical costs, and may lead to liver-related morbidity and mortality. There is an urgent need for accurate, reproducible, accessible, and noninvasive techniques appropriate for detecting and quantifying liver fat in the general population and for monitoring treatment response in at-risk patients. CT may play a potential role in opportunistic screening, and MRI proton-density fat fraction provides high accuracy for liver fat quantification; however, these imaging modalities may not be suited for widespread screening and surveillance, given the high global prevalence. US, a safe and widely available modality, is well positioned as a screening and surveillance tool. Although well-established qualitative signs of liver fat perform well in moderate and severe steatosis, these signs are less reliable for grading mild steatosis and are likely unreliable for detecting subtle changes over time. New and emerging quantitative biomarkers of liver fat, such as those based on standardized measurements of attenuation, backscatter, and speed of sound, hold promise. Evolving techniques such as multiparametric modeling, radiofrequency envelope analysis, and artificial intelligence-based tools are also on the horizon. The authors discuss the societal impact of fatty liver disease, summarize the current state of liver fat quantification with CT and MRI, and describe past, currently available, and potential future US-based techniques for evaluating liver fat. For each US-based technique, they describe the concept, measurement method, advantages, and limitations. © RSNA, 2023 Online supplemental material is available for this article. Quiz questions for this article are available through the Online Learning Center.

Liver Fibrosis, Fat, and Iron Evaluation with MRI and Fibrosis and Fat Evaluation with US: A Practical Guide for Radiologists

Quantitative imaging biomarkers of liver disease measured by using MRI and US are emerging as important clinical tools in the management of patients with chronic liver disease (CLD). Because of their high accuracy and noninvasive nature, in many cases, these techniques have replaced liver biopsy for the diagnosis, quantitative staging, and treatment monitoring of patients with CLD. The most commonly evaluated imaging biomarkers are surrogates for liver fibrosis, fat, and iron. MR elastography is now routinely performed to evaluate for liver fibrosis and typically combined with MRI-based liver fat and iron quantification to exclude or grade hepatic steatosis and iron overload, respectively. US elastography is also widely performed to evaluate for liver fibrosis and has the advantage of lower equipment cost and greater availability compared with those of MRI. Emerging US fat quantification methods can be performed along with US elastography. The author group, consisting of members of the Society of Abdominal Radiology (SAR) Liver Fibrosis Disease-Focused Panel (DFP), the SAR Hepatic Iron Overload DFP, and the European Society of Radiology, review the basics of liver fibrosis, fat, and iron quantification with MRI and liver fibrosis and fat quantification with US. The authors cover technical requirements, typical case display, quality control and proper measurement technique and case interpretation guidelines, pitfalls, and confounding factors. The authors aim to provide a practical guide for radiologists interpreting these examinations. ^© RSNA, 2023 See the invited commentary by Ronot in this issue. Quiz questions for this article are available in the supplemental material.

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.

Lithium Niobate Piezoelectric Micromachined Ultrasonic Transducers for high data-rate intrabody communication

In recent years, there has been an increased interest in continuous monitoring of patients and their Implanted Medical Devices (IMDs) with different wireless technologies such as ultrasounds. This paper demonstrates a high data-rate intrabody communication link based on Lithium Niobate (LN) Piezoelectric Micromachined Ultrasonic Transducers (pMUTs). The properties of the LN allow to activate multiple flexural mode of vibration with only top electrodes. When operating in materials like the human tissue, these modes are merging and forming a large communication bandwidth. Such large bandwidth, up to 400 kHz, allows for a high-data rate communication link for IMDs. Here we demonstrate a full communication link in a tissue phantom with a fabricated LN pMUT array of 225 elements with an area of just 3 by 3 mm square, showing data-rates up to 800 kbits/s, starting from 3.5 cm and going up to 13.5 cm, which covers the vast majority of IMDs.

This paper demonstrates Lithium Niobate based Piezoelectric Micromachined Ultrasonic Transducers. The devices allow activation of multiple vibration modes, achieve a high communication link, and monitor an implanted medical device up to 13.5 cm inside the body.

US Attenuation for Liver Fat Quantification: An AIUM-RSNA QIBA Pulse-Echo Quantitative Ultrasound Initiative

Nonalcoholic fatty liver disease (NAFLD) is the most common chronic liver disease worldwide, with an estimated prevalence of up to 30% in the general population and higher in people with type 2 diabetes. The assessment of liver fat content is essential to help identify patients with or who are at risk for NAFLD and to follow their disease over time. The American Institute of Ultrasound in Medicine-RSNA Quantitative Imaging Biomarkers Alliance Pulse-Echo Quantitative Ultrasound Initiative was formed to help develop and standardize acquisition protocols and to better understand confounding factors of US-based fat quantification. The three quantitative US parameters explored by the initiative are attenuation, backscatter coefficient, and speed of sound. The purpose of this review is to present the current state of attenuation imaging for fat quantification and to provide expert opinion on examination performance and interpretation. US attenuation methods that need further study are outlined.

Bottom layer absorption coefficients extraction from two-layer phantoms based on crossover point in diffuse reflectance

SIGNIFICANCE

Numerous optical imaging and spectroscopy techniques are used to study the tissue-optical properties; the majority of them are limited in information regarding the penetration depth. A simple, safe, easily applicable diagnostic technique is required to get deeper tissue information in a multilayer structure.

AIM

A fiber-based diffuse reflectance (DR) technique is used to extract and quantify the bottom layer absorption coefficients in two-layer (2L) tissue-mimicking solid phantoms. We determine the Indian black ink concentrations in a deep-hidden layer that is sandwiched between agar and silicone-based phantom layers.

APPROACH

A fiber-based DR experiment was performed to study the optical properties of the tissue at higher penetration depth, with different fiber core diameters and a constant numerical aperture (0.5 NA). The optimal core diameter of the fiber was chosen by measuring solid phantoms. In 2L phantoms, the thickness of the top layer was kept 5.5 mm with a constant absorption and reduced scattering coefficients (μa  =  0.045  mm  -  1 and μs  '    =  2.622  mm  -  1), whereas the absorption coefficients of the bottom layers were varied from 0.014 to 0.037  mm  -  1 keeping the μs  '   the same as the top layer. A unique crossover point (Cp) was found in the DR intensity profile against distance. We examined the slope before and after the Cp. These two slopes indicate the difference between the optical properties of the top and bottom layers. Our technique got further verification, as we successfully determined the Cp with different Indian black ink concentrations, placed at the junction between the agar and silicone-based phantom layers.

RESULTS

The DR measurements were applied to 2L phantoms. Two different slopes were found in 2L phantoms compared to the one-layer (optical properties equal to the top layer of 2L). We extracted the slopes before and after the Cp in the 2L phantoms. The calculated absorption coefficients before the Cp were 0.014  ±  0.0004, 0.022  ±  0.0003, 0.028  ±  0.0003, and 0.036  ±  0.0014  mm  -  1, and the absorption coefficients after the Cp were 0.019  ±  0.0013, 0.013  ±  0.0004, 0.014  ±  0.0006, and 0.031  ±  0.0001  mm  -  1, respectively. The calculated absorption coefficients before the Cp were in good agreement with the optical properties of the bottom layer. The calculated absorption coefficients after the Cp were not the same as the top layer. Our DR system successfully determines the crossover points 12.14  ±  0.11 and 11.73  ±  0.15  mm for 70% and 100% ink concentrations placed at the junction of the agar and silicone layers.

CONCLUSIONS

In a 2L tissue structure, the Cp depends on the absorption coefficients of top and bottom layers and the thickness of the top layer. With the help of the Cp and the absorption coefficients, one can determine the thickness of the top layer or vice versa. The slope value before the Cp in the DR profile allowed us to determine the absorption properties of the bottom layer instead of having the average behavior of the 2L phantom in the far detection range (11.0 to 17.0 mm).

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.

Reference Phantom Method for Ultrasonic Imaging of Thin Dynamic Constructs

Quantitative ultrasound has a great potential for the non-destructive evaluation of tissue engineered constructs, where the local attenuation and the integrated backscatter coefficient (IBC) can be used for monitoring the development of biological processes. The local determination of both parameters can be achieved using the reference phantom method (RPM). However, its accuracy can be affected when evaluating constructs of evolving sound speed, attenuation and thickness, for example, when evaluating biodegradable hydrogels developing neocartilage. To assess the feasibility of using the RPM under such dynamic conditions while employing a 50-MHz transducer, we conducted a series of experiments on 3-mm-thick acellular hydrogels laden with microspheres. The ultrasonic evaluation procedure used was validated by detecting and compensating for large attenuation variations occurring in the construct, up to 20-fold with respect to the reference phantom, with estimations errors below 1%. We found that sound speed mismatch does not affect the local attenuation estimation, but causes a strong diffraction effect by reducing the backscatter intensity. Such intensity reduction was compensated by determining the IBC percentage change (IBC_Δ) as function of sound speed mismatch with respect to the reference phantom (ΔSS), with the equation IBC_Δ = (0.63 ± 0.07) ΔSS + (8.54 ± 0.76) 10^-3 ΔSS². The investigated ΔSS interval was up to 120 m/s and using two different concentrations of microspheres, with estimation errors below 7% relative to the construct's actual IBC. Finally, we found that the spectral difference method is sufficient to measure within a few millimetres in depth mismatch, and when combining with sound speed mismatch, we found negligible additional effects. These results pave the way for the use of a generic reference phantom for the evaluation of thin dynamic constructs, thus simplifying the need for using different phantoms depending on the construct's properties.

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.

Attenuation Coefficient Parameter Computations for Tissue Composition Assessment of Carotid Atherosclerotic Plaque in Vivo

Quantitative ultrasound has been used to assess carotid plaque tissue composition. Here, we compute the attenuation coefficient (AC) in vivo with the optimum power spectral shift estimator (OPSSE) and reference phantom method (RPM), extract AC parameters and form parametric maps. Differences between OPSSE and RPM AC parameters are computed. Relationships between AC parameters, surgical scores and histopathology assessments are examined. Kendall's τ correlations between OPSSE AC and surgical scores are significant, including those between cholesterol and Standard Deviation (adjusted p = 0.038); thrombus and Minimum (adjusted p = 0.002), Maximum (adjusted p = 0.021) and Standard Deviation (adjusted p = 0.001); ulceration and Average (adjusted p = 0.033), Median (unadjusted p = 0.013), Maximum (unadjusted p = 0.039) and Mode (adjusted p = 0.009). The strongest correlations with histopathology are percentage cholesterol and Median OPSSE (unadjusted p = 0.007); percentage hemorrhage and Minimum OPSSE (adjusted p < 0.001); hemosiderin score and Median OPSSE (adjusted p = 0.010); and percentage calcium and Percentage Non-physical RPM Pixels (unadjusted p = 0.014). Kruskal-Wallis H and Dunn's post hoc tests have the ability to distinguish between groups (p < 0.05). Results suggest AC parameters may assist in vivo evaluation of carotid plaque vulnerability.

Quantitative ultrasound radiofrequency data analysis for the assessment of hepatic steatosis using the controlled attenuation parameter as a reference standard

Purpose

This study was aimed to investigate the value of quantitative ultrasound (US) parameters from radiofrequency (RF) data analysis for assessing hepatic steatosis, using controlled attenuation parameter (CAP)-based steatosis grades as the reference standard.

Methods

We analyzed 243 participants with both B-mode liver US with RF data acquisition and CAP measurements. On B-mode US images, hepatic steatosis was visually scored (0/1/2/3, none/mild/moderate/severe), and the hepatorenal index (HRI) was calculated. From the RF data analysis, the tissue scatter-distribution imaging parameter (TSI-p) and tissue attenuation imaging parameter (TAI-p) of the liver parenchyma were measured. US parameters were correlated with CAP-based steatosis grades (S0/1/2/3, none/mild/moderate/severe) and their diagnostic performance was evaluated using receiver operating characteristic (ROC) curve analysis. Multivariate linear regression analysis was performed to identify determinants of TSI-p and TAI-p.

Results

Participants were classified as having S0 (n=152), S1 (n=54), S2 (n=14), and S3 (n=23) on CAP measurements. TSI-p and TAI-p were significantly correlated with steatosis grades (ρ =0.593 and ρ=-0.617, P<0.001 for both). For predicting ≥S1, ≥S2, and S3, the areas under the ROC curves (AUCs) of TSI-p were 0.827/0.914/0.917; TAI-p, 0.844/0.914/0.909; visual scores, 0.659/0.778/0.794; and HRI, 0.629/0.751/0.759, respectively. TSI-p and TAI-p had significantly higher AUCs than did visual scores or HRI for ≥S1 or ≥S2 (P≤0.003). In the multivariate analysis, the transient elastography-based fibrosis grade (P=0.034) and steatosis grade (P<0.001) were independent determinants of TSI-p, while steatosis grade (P<0.001) was an independent determinant of TAI-p.

Conclusion

TSI-p and TAI-p derived from US RF data may be useful for detecting hepatic steatosis and assessing its severity.

Quantitative Ultrasound Monitoring of Breast Tumour Response to Neoadjuvant Chemotherapy: Comparison of Results Among Clinical Scanners

Quantitative ultrasound (QUS) techniques have been demonstrated to detect cell death in vitro and in vivo. Recently, multi-feature classification models have been incorporated into QUS texture-feature analysis methods to increase further the sensitivity and specificity of detecting treatment response in locally advanced breast cancer patients. To effectively incorporate these analytic methods into clinical applications, QUS and texture-feature estimations should be independent of data acquisition systems. The study here investigated the consistencies of QUS and texture-feature estimation techniques relative to several factors. These included the ultrasound system properties, the effects of tissue heterogeneity and the effects of these factors on the monitoring of response to neoadjuvant chemotherapy. Specifically, tumour-response-detection performance based on QUS and texture parameters using two clinical ultrasound systems was compared. Observed variations in data between the systems were small and the results exhibited good agreement in tumour response predictions obtained from both ultrasound systems. The results obtained in this study suggest that tissue heterogeneity was a dominant feature in the parameters measured with the two different ultrasound systems; whereas differences in ultrasound system beam properties only exhibited a minor impact on texture features. The McNemar statistical test performed on tumour response prediction results from the two systems did not reveal significant differences. Overall, the results in this study demonstrate the potential to achieve reliable and consistent QUS and texture-based analyses across different ultrasound imaging platforms.

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.

A Quantitative Ultrasound-Based Multi-Parameter Classifier for Breast Masses

This manuscript reports preliminary results obtained by combining estimates of two or three (among seven) quantitative ultrasound (QUS) parameters in a model-free, multi-parameter classifier to differentiate breast carcinomas from fibroadenomas (the most common benign solid tumor). Forty-three patients scheduled for core biopsy of a suspicious breast mass were recruited. Radiofrequency echo signal data were acquired using clinical breast ultrasound systems equipped with linear array transducers. The reference phantom method was used to obtain system-independent estimates of the specific attenuation (ATT), the average backscatter coefficients, the effective scatterer diameter (ESD) and an effective scatterer diameter heterogeneity index (ESDHI) over regions of interest within each mass. In addition, the envelope amplitude signal-to-noise ratio (SNR), the Nakagami shape parameter, m, and the maximum collapsed average (maxCA) of the generalized spectrum were also computed. Classification was performed using the minimum Mahalanobis distance to the centroids of the training classes and tested against biopsy results. Classification performance was evaluated with the area under the receiver operating characteristic (ROC) curve. The best performance with a two-parameter classifier used the ESD and ESDHI and resulted in an area under the ROC curve of 0.98 (95% confidence interval [CI]: 0.95-1.00). Classification performance improved with three parameters (ATT, ESD and ESDHI) yielding an area under the ROC curve of 0.999 (0.995-1.000). These results suggest that system-independent QUS parameters, when combined in a model-free classifier, are a promising tool to characterize breast tumors. A larger study is needed to further test this idea.

Quantitative ultrasound imaging of cell-laden hydrogels and printed constructs

In the present work we have revisited the application of quantitative ultrasound imaging (QUI) to cellular hydrogels, by using the reference phantom method (RPM) in combination with a local attenuation compensation algorithm. The investigated biological samples consisted of cell-laden collagen hydrogels with PC12 neural cells. These cell-laden hydrogels were used to calibrate the integrated backscattering coefficient (IBC) as a function of cell density, which was then used to generate parametric images of local cell density. The image resolution used for QUI and its impact on the relative IBC error was also investigated. Another important contribution of our work was the monitoring of PC12 cell proliferation. The cell number estimates obtained via the calibrated IBC compared well with data obtained using a conventional quantitative method, the MTS assay. Evaluation of spectral changes as a function of culture time also provided additional information on the cell cluster size, which was found to be in close agreement with that observed by microscopy. Last but not least, we also applied QUI on a 3D printed cellular construct in order to illustrate its capabilities for the evaluation of bioprinted structures. STATEMENT OF SIGNIFICANCE: While there is intensive research in the areas of polymer science, biology, and 3D bio-printing, there exists a gap in available characterisation tools for the non-destructive inspection of biological constructs in the three-dimensional domain, on the macroscopic scale, and with fast data acquisition times. Quantitative ultrasound imaging is a suitable characterization technique for providing essential information on the development of tissue engineered constructs. These results provide a detailed and comprehensive guide on the capabilities and limitations of the technique.

Quantitative Ultrasound Biomarkers Based on Backscattered Acoustic Power: Potential for Quantifying Remodeling of the Human Cervix during Pregnancy

As pregnancy progresses, the cervix remodels from a rigid structure to one pliable enough to allow delivery of a fetus, a process that involves progressive disorganization of cervical microstructure. Quantitative ultrasound biomarkers that may detect this process include those derived from the backscattered echo signal, namely, acoustic attenuation and backscattered power loss. We recently reported that attenuation and backscattered power loss are affected by tissue anisotropy and heterogeneity in the ex vivo cervix. In this study, we compared attenuation and backscattered power difference in a group of women in early pregnancy (first trimester) with those in a group in late pregnancy (third trimester). We observed a significant decrease in the backscattered power difference in late as compared with early pregnancy, suggesting decreased microstructural organization in late pregnancy, a finding that is consistent with animal models of cervical remodeling. In contrast, we found no difference in attenuation between the time points. These results suggest that the backscattered power difference, but perhaps not attenuation, may be a useful clinical biomarker of cervical remodeling.

Quantitative Hepatic Fat Quantification in Non-alcoholic Fatty Liver Disease Using Ultrasound-Based Techniques: A Review of Literature and Their Diagnostic Performance

Non-alcoholic fatty liver disease is a condition that is characterized by the presence of >5% fat in the liver and affects more than one billion people worldwide. If adequate and early precautions are not taken, non-alcoholic fatty liver disease can progress to cirrhosis and death. The current reference standard for detecting hepatic steatosis is a liver biopsy. However, because of the potential morbidity associated with liver biopsies, non-invasive imaging biomarkers have been extensively investigated. Magnetic resonance imaging-based methods have proven accuracy in quantifying liver steatosis; however, these techniques are costly and have limited availability. Ultrasound-based quantitative imaging techniques are increasingly utilized because of their widespread availability, ease of use and relative cost-effectiveness. Several ultrasound-based liver fat quantification techniques have been investigated, including techniques that measure changes in the acoustic properties of the liver caused by the presence of fat. In this review, we focus on quantitative ultrasound approaches and their diagnostic performance in the realm of non-alcoholic fatty liver disease.

Power Spectrum Consistency among Systems and Transducers

Use of the reference phantom method for computing acoustic attenuation and backscatter is widespread. However, clinical application of these methods has been limited by the need to acquire reference phantom data. We determined that the data acquired from 11 transducers of the same model and five clinical ultrasound systems of the same model produce equivalent estimates of reference phantom power spectra. We describe that the contribution to power spectral density variance among systems and transducers equals that from speckle variance with 59 uncorrelated echo signals. Thus, when the number of uncorrelated lines of data is small, speckle variance will dominate the power spectral density estimate variance introduced by different systems and transducers. These results suggest that, at least for this particular transducer and imaging system combination, one set of reference phantom calibration data is highly representative of the average among equivalent transducers and systems that are in good working order.

Low Variance Estimation of Backscatter Quantitative Ultrasound Parameters Using Dynamic Programming

One of the main limitations of ultrasound imaging is that image quality and interpretation depend on the skill of the user and the experience of the clinician. Quantitative ultrasound (QUS) methods provide objective, system-independent estimates of tissue properties, such as acoustic attenuation and backscattering properties of tissue, which are valuable as objective tools for both diagnosis and intervention. Accurate and precise estimation of these properties requires correct compensation for intervening tissue attenuation. Prior attempts to estimate intervening-tissue attenuation based on minimizing cost functions that compared backscattered echo data to models have resulted in limited precision and accuracy. To overcome these limitations, in this paper, we incorporate the prior information of piecewise continuity of QUS parameters as a regularization term into our cost function. We further propose to calculate this cost function using dynamic programming (DP), a computationally efficient optimization algorithm that finds the global optimum. Our results on tissue-mimicking phantoms show that DP substantially outperforms a published least squares method in terms of both estimation bias and variance.

Anisotropy and Spatial Heterogeneity in Quantitative Ultrasound Parameters: Relevance to the Study of the Human Cervix

Imaging biomarkers based on quantitative ultrasound can offer valuable information about properties that inform tissue function and behavior such as microstructural organization (e.g., collagen alignment) and viscoelasticity (i.e., compliance). For example, the cervix feels softer as its microstructure remodels during pregnancy, an increase in compliance that can be objectively quantified with shear wave speed and therefore shear wave speed estimation is a potential biomarker of cervical remodeling. Other proposed biomarkers include parameters derived from the backscattered echo signal, such as attenuation and backscattered power loss, because such parameters can provide insight into tissue microstructural alignment and organization. Of these, attenuation values for the pregnant cervix have been reported, but large estimate variance reduces their clinical value. That said, parameter estimates based on the backscattered echo signal may be incorrect if assumptions they rely on, such as tissue isotropy and homogeneity, are violated. For that reason, we explored backscatter and attenuation parameters as potential biomarkers of cervical remodeling via careful investigation of the assumptions of isotropy and homogeneity in cervical tissue. Specifically, we estimated the angle- and spatial-dependence of parameters of backscattered power and acoustic attenuation in the ex vivo human cervix, using the reference phantom method and electronic steering of the ultrasound beam. We found that estimates are anisotropic and spatially heterogeneous, presumably because the tissue itself is anisotropic and heterogeneous. We conclude that appropriate interpretation of imaging biomarkers of cervical remodeling must account for tissue anisotropy and heterogeneity.

High-Frequency Quantitative Ultrasound for Imaging Prostate Cancer Using a Novel Micro-Ultrasound Scanner

Currently, biopsies guided by transrectal ultrasound (TRUS) are the only method for definitive diagnosis of prostate cancer. Studies by our group suggest that quantitative ultrasound (QUS) could provide a more sensitive means of targeting biopsies and directing focal treatments to cancer-suspicious regions in the prostate. Previous studies have utilized ultrasound signals at typical clinical frequencies, i.e., in the 6-MHz range. In the present study, a 29-MHz, TRUS, micro-ultrasound system and transducer (ExactVu micro-ultrasound, Exact Imaging, Markham, Canada) was used to acquire radio frequency data from 163 patients immediately before 12-core biopsy procedures, comprising 1956 cores. These retrospective data are a subset of data acquired in an ongoing, multisite, 2000-patient, randomized, clinical trial (clinicaltrials.gov NCT02079025). Spectrum-based QUS estimates of effective scatter diameter (ESD), effective acoustic concentration (EAC), midband (M), intercept (I) and slope (S) as well as envelope statistics employing a Nakagami distribution were used to train linear discriminant classifiers (LDCs) and support vector machines (SVMs). Classifier performance was assessed using area-under-the-curve (AUC) values obtained from receiver operating characteristic (ROC) analyses with 10-fold cross validation. A combination of ESD and EAC parameters resulted in an AUC value of 0.77 using a LDC. When Nakagami-µ or prostate-specific antigen (PSA) values were added as features, the AUC value increased to 0.79. SVM produced an AUC value of 0.77, using a combination of envelope and spectral QUS estimates. The best classification produced an AUC value of 0.81 using an LDC when combining envelope statistics, PSA, ESD and EAC. In a previous study, B-mode-based scoring and evaluation using the PRI-MUS protocol produced a maximal AUC value of 0.74 for higher Gleason-score values (GS >7) when read by an expert. Our initial results with AUC values of 0.81 are very encouraging for developing a new, predominantly user-independent, prostate-cancer, risk-assessing tool.

Quantifying Backscatter Anisotropy Using the Reference Phantom Method

Acoustic properties can be exploited to infer and evaluate tissue microstructure. However, common assumptions are that the medium of interest is homogeneous and isotropic, and that its underlying physical properties cause diffuse scattering. In this paper, we describe how we developed and tested novel parameters designed to address isotropy/anisotropy in backscattered echo signal power in complex biological tissues. Specifically, we explored isotropy/anisotropy in backscattered power in isotropic phantoms (spherical glass beads), an anisotropic phantom (dialysis phantom with rodlike fibers), and an in vivo human tissue with well-described anisotropy (bicep muscle). Our approach uses the reference phantom method to compensate for system transfer and diffraction losses when electronically beamsteering a linear array transducer. We define three parameters to quantify the presence and orientation of anisotropic scatterers, as well as address magnitude of anisotropy. We found that these parameters can detect and sense the degree of anisotropy in backscatter in both phantoms and bicep muscle. Bias of the summary anisotropy parameters, induced through a speed of sound mismatch of sample media and reference phantom, was less than 0.2 dB if the speed of sound was within ±20 m/s of the sample media. In summary, these new parameters may be useful for testing the assumption of isotropy as well as providing more detailed information about the underlying microstructural sources of backscatter in complex biological tissues.

Monitoring Microwave Ablation of Ex Vivo Bovine Liver Using Ultrasonic Attenuation Imaging

Thermal ablation of soft tissue changes the tissue microstructure and, consequently, induces changes in its acoustic properties. Although B-mode ultrasound provides high-resolution and high-frame-rate images of ablative therapeutic procedures, it is not particularly effective at delineating boundaries of ablated regions because of poor contrast in echogenicity between ablated and surrounding normal tissue. Quantitative ultrasound techniques can provide quantitative estimates of acoustic properties, such as backscatter and attenuation coefficients, and differentiate ablated and unablated regions more effectively, with the potential for monitoring minimally invasive thermal therapies. In this study, a previously introduced attenuation estimation method was used to create quantitative attenuation coefficient maps for 11 microwave ablation procedures performed on refrigerated ex vivo bovine liver. The attenuation images correlate well with the pathologic images of the ablated region. The mean attenuation coefficient for regions of interest drawn inside and outside the ablated zones were 0.9 (±0.2) and 0.45 (±0.15) dB/cm/MHz, respectively. These estimates agree with reported values in the literature and establish the usefulness of non-invasive attenuation imaging for monitoring therapeutic procedures in the liver.

Lower Bound on Estimation Variance of the Ultrasonic Attenuation Coefficient Using the Spectral-Difference Reference-phantom Method

Ultrasonic attenuation is one of the primary parameters of interest in Quantitative Ultrasound (QUS). Non-invasive monitoring of tissue attenuation can provide valuable diagnostic and prognostic information to the physician. The Reference Phantom Method (RPM) was introduced as a way of mitigating some of the system-related effects and biases to facilitate clinical QUS applications. In this paper, under the assumption of diffuse scattering, a probabilistic model of the backscattered signal spectrum is used to derive a theoretical lower bound on the estimation variance of the attenuation coefficient using the Spectral-Difference RPM. The theoretical lower bound is compared to simulated and experimental attenuation estimation statistics in tissue-mimicking (TM) phantoms. Estimation standard deviation (STD) of the sample attenuation in a region of interest (ROI) of the TM phantom is measured for various combinations of processing parameters, including Radio-Frequency (RF) data block length (i.e., window length) from 3 to 17 mm, RF data block width from 10 to 100 A-lines, and number of RF data blocks per attenuation estimation ROI from 3 to 10. In addition to the Spectral-Difference RPM, local attenuation estimation for simulated and experimental data sets was also performed using a modified implementation of the Spectral Fit Method (SFM). Estimation statistics of the SFM are compared to theoretical variance predictions from the literature.1 Measured STD curves are observed to lie above the theoretical lower bound curves, thus experimentally verifying the validity of the derived bounds. This theoretical framework benefits tissue characterization efforts by isolating processing parameter ranges that could provide required precision levels in estimation of the ultrasonic attenuation coefficient using Spectral Difference methods.

In Vivo Estimation of Attenuation and Backscatter Coefficients From Human Thyroids

Fine-needle aspiration (FNA) remains the gold standard for the diagnosis of thyroid cancer. However, currently, a large number of FNA biopsies result in negative or undetermined diagnosis, which suggests that better noninvasive tools are needed for the clinical management of thyroid cancer. Spectral-based quantitative ultrasound (QUS) characterizations may offer a better diagnostic management as previously demonstrated in mouse cancer models ex vivo. As a first step toward understanding the potential of QUS markers for thyroid disease management, this paper deals with the spectral-based QUS estimation of healthy human thyroids in vivo. Twenty volunteers were inspected by a trained radiologist using two ultrasonic imaging systems, which allowed them to acquire radio-frequency data spanning the 3-16-MHz frequency range. Estimates of attenuation coefficient slope (ACS) using the spectral logarithmic difference method had an average value of [Formula: see text]) with a standard deviation of [Formula: see text]. Estimates of backscatter coefficient (BSC) using the reference-phantom method had an average value of [Formula: see text] over the useful frequency range. The intersubject variability when estimating BSCs was less than 1.5 dB over the analysis frequency range. Further, the effectiveness of three scattering models (i.e., fluid sphere, Gaussian, and exponential form factors) when fitting the experimentally estimated BSCs was assessed. The exponential form factor was found to provide the best overall goodness of fit ( R² = 0.917), followed by the Gaussian ( R² = 0.807) and the fluid-sphere models ( R² = 0.752). For all scattering models used in this study, average estimates of the effective scatterer diameter were between 44 and 56 μm. Overall, an excellent agreement in the estimated attenuation and BSCs with both scanners was exhibited.

Optimum Diffraction-Corrected Frequency-Shift Estimator of the Ultrasonic Attenuation Coefficient

The ultrasonic attenuation coefficient is an important parameter that has been studied extensively in Quantitative Ultrasound and Tissue Characterization. There are various methods described in the literature that estimate this parameter by measuring either spectral difference (i.e., decay) or spectral shift of the backscattered echo signal. Under ideal conditions, i.e., in the absence of abrupt changes in tissue backscattering, Spectral Difference methods can produce estimates with high accuracy and precision. On the other hand, diffraction-corrected Spectral Shift methods (e.g., the Hybrid method) are better suited for application in practical settings using clinical ultrasound scanners. However, current Spectral Shift methods use inefficient frequency shift estimators that ultimately degrade the quality of attenuation coefficient estimates. In this paper, a probabilistic model of the backscattered radiofrequency (RF) echo is used to derive the Cramér-Rao lower bound (CRLB) on estimation variance of the spectral centroid. Next, an efficient correlation-based shift estimator is presented that achieves the CRLB. Used in conjunction with a well-characterized reference phantom to correct for diffraction and other system-related effects, this estimator greatly improves the accuracy and precision of Spectral- Shift attenuation estimation. A theoretical analysis of this method is provided, and its performance is quantitatively compared with that of the Hybrid method using simulated and experimental phantom studies. A minimum of 3-fold reduction in the standard deviation of attenuation coefficient estimates is observed using the new method.

Review of Quantitative Ultrasound: Envelope Statistics and Backscatter Coefficient Imaging and Contributions to Diagnostic Ultrasound

Conventional medical imaging technologies, including ultrasound, have continued to improve over the years. For example, in oncology, medical imaging is characterized by high sensitivity, i.e., the ability to detect anomalous tissue features, but the ability to classify these tissue features from images often lacks specificity. As a result, a large number of biopsies of tissues with suspicious image findings are performed each year with a vast majority of these biopsies resulting in a negative finding. To improve specificity of cancer imaging, quantitative imaging techniques can play an important role. Conventional ultrasound B-mode imaging is mainly qualitative in nature. However, quantitative ultrasound (QUS) imaging can provide specific numbers related to tissue features that can increase the specificity of image findings leading to improvements in diagnostic ultrasound. QUS imaging can encompass a wide variety of techniques including spectral-based parameterization, elastography, shear wave imaging, flow estimation, and envelope statistics. Currently, spectral-based parameterization and envelope statistics are not available on most conventional clinical ultrasound machines. However, in recent years, QUS techniques involving spectral-based parameterization and envelope statistics have demonstrated success in many applications, providing additional diagnostic capabilities. Spectral-based techniques include the estimation of the backscatter coefficient (BSC), estimation of attenuation, and estimation of scatterer properties such as the correlation length associated with an effective scatterer diameter (ESD) and the effective acoustic concentration (EAC) of scatterers. Envelope statistics include the estimation of the number density of scatterers and quantification of coherent to incoherent signals produced from the tissue. Challenges for clinical application include correctly accounting for attenuation effects and transmission losses and implementation of QUS on clinical devices. Successful clinical and preclinical applications demonstrating the ability of QUS to improve medical diagnostics include characterization of the myocardium during the cardiac cycle, cancer detection, classification of solid tumors and lymph nodes, detection and quantification of fatty liver disease, and monitoring and assessment of therapy.

Beyond Cervical Length: A Pilot Study of Ultrasonic Attenuation for Early Detection of Preterm Birth Risk

The purpose of this study was to determine whether cervical ultrasonic attenuation could identify women at risk of spontaneous preterm birth. During pregnancy, women (n = 67) underwent from one to five transvaginal ultrasonic examinations to estimate cervical ultrasonic attenuation and cervical length. Ultrasonic data were obtained with a Zonare ultrasound system with a 5- to 9-MHz endovaginal transducer and processed offline. Cervical ultrasonic attenuation was lower at 17-21 wk of gestation in the SPTB group (1.02 dB/cm-MHz) than in the full-term birth groups (1.34 dB/cm-MHz) (p = 0.04). Cervical length was shorter (3.16 cm) at 22-26 wk in the SPTB group than in the women delivering full term (3.68 cm) (p = 0.004); cervical attenuation was not significantly different at this time point. These findings suggest that low attenuation may be an additional early cervical marker to identify women at risk for SPTB.

Noninvasive Diagnosis of Nonalcoholic Fatty Liver Disease and Quantification of Liver Fat Using a New Quantitative Ultrasound Technique

BACKGROUND & AIMS

Liver biopsy analysis is the standard method used to diagnose nonalcoholic fatty liver disease (NAFLD). Advanced magnetic resonance imaging is a noninvasive procedure that can accurately diagnose and quantify steatosis, but is expensive. Conventional ultrasound is more accessible but identifies steatosis with low levels of sensitivity, specificity, and quantitative accuracy, and results vary among technicians. A new quantitative ultrasound (QUS) technique can identify steatosis in animal models. We assessed the accuracy of QUS in the diagnosis and quantification of hepatic steatosis, comparing findings with those from magnetic resonance imaging proton density fat fraction (MRI-PDFF) analysis as a reference.

METHODS

We performed a prospective, cross-sectional analysis of a cohort of adults (N = 204) with NAFLD (MRI-PDFF, ≥5%) and without NAFLD (controls). Subjects underwent MRI-PDFF and QUS analyses of the liver on the same day at the University of California, San Diego, from February 2012 through March 2014. QUS parameters and backscatter coefficient (BSC) values were calculated. Patients were assigned randomly to training (n = 102; mean age, 51 ± 17 y; mean body mass index, 31 ± 7 kg/m(2)) and validation (n = 102; mean age, 49 ± 17 y; body mass index, 30 ± 6 kg/m(2)) groups; 69% of patients in each group had NAFLD.

RESULTS

BSC (range, 0.00005-0.25 1/cm-sr) correlated with MRI-PDFF (Spearman ρ = 0.80; P < .0001). In the training group, the BSC analysis identified patients with NAFLD with an area under the curve value of 0.98 (95% confidence interval, 0.95-1.00; P < .0001). The optimal BSC cut-off value identified patients with NAFLD in the training and validation groups with 93% and 87% sensitivity, 97% and 91% specificity, 86% and 76% negative predictive values, and 99% and 95% positive predictive values, respectively.

CONCLUSIONS

QUS measurements of BSC can accurately diagnose and quantify hepatic steatosis, based on a cross-sectional analysis that used MRI-PDFF as the reference. With further validation, QUS could be an inexpensive, widely available method to screen the general or at-risk population for NAFLD.

The mechanical role of the cervix in pregnancy

Appropriate mechanical function of the uterine cervix is critical for maintaining a pregnancy to term so that the fetus can develop fully. At the end of pregnancy, however, the cervix must allow delivery, which requires it to markedly soften, shorten and dilate. There are multiple pathways to spontaneous preterm birth, the leading global cause of death in children less than 5 years old, but all culminate in premature cervical change, because that is the last step in the final common pathway to delivery. The mechanisms underlying premature cervical change in pregnancy are poorly understood, and therefore current clinical protocols to assess preterm birth risk are limited to surrogate markers of mechanical function, such as sonographically measured cervical length. This is what motivates us to study the cervix, for which we propose investigating clinical cervical function in parallel with a quantitative engineering evaluation of its structural function. We aspire to develop a common translational language, as well as generate a rigorous integrated clinical-engineering framework for assessing cervical mechanical function at the cellular to organ level. In this review, we embark on that challenge by describing the current landscape of clinical, biochemical, and engineering concepts associated with the mechanical function of the cervix during pregnancy. Our goal is to use this common platform to inspire novel approaches to delineate normal and abnormal cervical function in pregnancy.

Performance evaluation of the spectral centroid downshift method for attenuation estimation

Estimation of frequency-dependent ultrasonic attenuation is an important aspect of tissue characterization. Along with other acoustic parameters studied in quantitative ultrasound, the attenuation coefficient can be used to differentiate normal and pathological tissue. The spectral centroid downshift (CDS) method is one the most common frequencydomain approaches applied to this problem. In this study, a statistical analysis of this method's performance was carried out based on a parametric model of the signal power spectrum in the presence of electronic noise. The parametric model used for the power spectrum of received RF data assumes a Gaussian spectral profile for the transmit pulse, and incorporates effects of attenuation, windowing, and electronic noise. Spectral moments were calculated and used to estimate second-order centroid statistics. A theoretical expression for the variance of a maximum likelihood estimator of attenuation coefficient was derived in terms of the centroid statistics and other model parameters, such as transmit pulse center frequency and bandwidth, RF data window length, SNR, and number of regression points. Theoretically predicted estimation variances were compared with experimentally estimated variances on RF data sets from both computer-simulated and physical tissue-mimicking phantoms. Scan parameter ranges for this study were electronic SNR from 10 to 70 dB, transmit pulse standard deviation from 0.5 to 4.1 MHz, transmit pulse center frequency from 2 to 8 MHz, and data window length from 3 to 17 mm. Acceptable agreement was observed between theoretical predictions and experimentally estimated values with differences smaller than 0.05 dB/cm/MHz across the parameter ranges investigated. This model helps predict the best attenuation estimation variance achievable with the CDS method, in terms of said scan parameters.

Effective scatterer diameter estimates for broad scatterer size distributions

Acoustic form factors have been used to model the frequency dependence of acoustic scattering in phantoms and tissues. This work demonstrates that a broad range of scatterer sizes, individually well represented by Faran theory or a Gaussian form factor, is not accurately described by a single effective scatterer from either of these models. Contributions from a distribution of discrete scatterer sizes for two different form factor functions (Gaussian form factors and scattering functions from Faran's theory) were calculated and linearly combined. Composite form factors created from Gaussian distributions of scatterer sizes centered at 50 µm with standard deviations of up to σ = 40 µm were fit to each scattering model between 2 and 12 MHz. Scatterer distributions were generated using one of two assumptions: the number density of the scatterer diameter distribution was Gaussian distributed, or the volume fraction of each scatterer diameter in the distribution was Gaussian distributed. Each simulated form factor was fit to a single-diameter form factor model for Gaussian and exponential form factors. The mean-squared error (MSE) between the composite simulated data and the best-fit single-diameter model was smaller with an exponential form factor model, compared with a Gaussian model, for distributions with standard deviations larger than 30% of the centroid value. In addition, exponential models were shown to have better ability to distinguish between Faran scattering model-based distributions with varying center diameters than the Gaussian form factor model. The evidence suggests that when little is known about the scattering medium, an exponential scattering model provides a better first approximation to the scattering correlation function for a broad distribution of spherically symmetric scatterers than when a Gaussian form factor model is assumed.

Experimental evaluation of spectral-based quantitative ultrasound imaging using plane wave compounding

Quantitative ultrasound (QUS) based on backscatter coefficient (BSC) estimation has shown potential for tissue characterization. Beamforming using plane wave compounding has advantages for echographic, Doppler, and elastographic imaging; however, to date, plane wave compounding has not been experimentally evaluated for the purpose of BSC estimation. In this study, two BSC-derived parameters (i.e., the BSC midband fit and intercept) were estimated from experimental data obtained using compound plane wave beamforming. For comparison, QUS parameters were also estimated from data obtained using both fixed focus and dynamic receive beamforming. An ultrasound imaging system equipped with a 9-MHz center frequency, 64-element array was used to collect data up to a depth of 45 mm. Two gelatin phantoms with randomly distributed 20-μm inclusions with a homogeneous scatterer concentration and a two-region scatterer concentration were used for assessing the precision and lateral resolution of QUS imaging, respectively. The use of plane wave compounding resulted in accurate QUS estimation (i.e., bias in the BSC parameters of less than 2 dB) and relatively constant lateral resolution (i.e., BSC midband fit 10% to 90% rise distance ranging between 1.0 and 1.5 mm) throughout a 45 mm field of view. Although both fixed focus and dynamic receive beamforming provided the same performance around the focal depth, the reduction in SNR away from the focus resulted in a reduced field of view in the homogeneous phantom (i.e., only 28 mm). The lateral resolution also degraded away from the focus, with up to a 2-fold and 10-fold increase in the rise distance at 20 mm beyond the focal depth for dynamic receive and fixed focus beamforming, respectively. These results suggest that plane wave compounding has the potential to improve the performance of spectral-based quantitative ultrasound over other conventional beamforming strategies.

Signal to noise ratio comparisons for ultrasound attenuation slope estimation algorithms

PURPOSE

Attenuation imaging has a promising role in the detection of tissue abnormalities. The authors have previously compared three different frequency domain ultrasound attenuation estimation methods, for accuracy and bias. The mean estimated attenuation value in a region of interest has been the determining factor of how well a method performs; however, the noise level has not been quantified for attenuation estimated using different methods.

METHODS

The authors compare three different frequency domain ultrasound attenuation estimation methods [the reference phantom method (RPM), the centroid downshift method (CEN), and the hybrid method (HYB)] using the signal to noise ratio (SNR) metric. Both simulated and experimental tissue-mimicking phantoms are used in the performance comparison study, evaluating the impact of the variation in acoustical properties.

RESULTS

For attenuation estimation in a tissue-mimicking phantom with a known attenuation coefficient of 0.5 dB/cm/MHz, all the three methods estimated the attenuation coefficient to be ≈ 0.49 dB/cm/MHz for a transmit center frequency of 6 MHz, however, the signal to noise ratio obtained was found to be 8.5, 5.7, and 2.2 for the HYB, RPM, and CEN methods, respectively. These results demonstrate the need for the SNR metric in the comparison of different algorithms and to evaluate the impact of varying different ultrasound system and tissue parameters.

CONCLUSIONS

In this paper, the authors demonstrate that although the estimated mean attenuation value with a region of interest may be closely estimated using different methods, the signal to noise ratio obtained of the estimates can vary significantly. The centroid downshift method presented with the lowest signal-to-noise ratio of the methods compared. The hybrid method was the least susceptible to changes in the acoustical properties and provided unbiased attenuation coefficient estimates with the highest signal-to-noise ratios.

Real-time implementation of a dual-mode ultrasound array system: in vivo results

A real-time dual-mode ultrasound array (DMUA) system for imaging and therapy is described. The system utilizes a concave (40-mm radius of curvature) 3.5 MHz, 32 element array, and modular multichannel transmitter/receiver. The system is capable of operating in a variety of imaging and therapy modes (on transmit) and continuous receive on all array elements even during high-power operation. A signal chain consisting of field-programmable gate arrays and graphical processing units is used to enable real time, software-defined beamforming and image formation. Imaging data, from quality assurance phantoms as well as in vivo small- and large-animal models, are presented and discussed. Corresponding images obtained using a temporally-synchronized and spatially-aligned diagnostic probe confirm the DMUA's ability to form anatomically-correct images with sufficient contrast in an extended field of view around its geometric center. In addition, high-frame rate DMUA data also demonstrate the feasibility of detection and localization of echo changes indicative of cavitation and/or tissue boiling during high-intensity focused ultrasound exposures with 45-50 dB dynamic range. The results also show that the axial and lateral resolution of the DMUA are consistent with its f(number) and bandwidth with well-behaved speckle cell characteristics. These results point the way to a theranostic DMUA system capable of quantitative imaging of tissue property changes with high specificity to lesion formation using focused ultrasound.

Method for estimating total attenuation from a spatial map of attenuation slope for quantitative ultrasound imaging

Estimating total ultrasound attenuation from backscatter data is essential in the field of quantitative ultrasound (QUS) because of the need to compensate for attenuation when estimating the backscatter coefficient and QUS parameters. This work uses a reference phantom method of attenuation estimation to create a spatial map of attenuation slope (AS) from backscatter radio-frequency (RF) data of three phantoms and a rat mammary adenocarcinoma tumor (MAT). The attenuation maps show changes in attenuation between different regions of the phantoms and the MAT tumor. Analyses of the attenuation maps of the phantoms suggest that the AS estimates are in good quantitative agreement with the known values for the phantoms. Furthermore, estimates of total attenuation from the attenuation maps are likewise in good quantitative agreement with known values.