Analysis of Two Quantitative Ultrasound Approaches

Ultrasound normalized cumulative residual entropy imaging: Theory, methodology, and application

BACKGROUND AND OBJECTIVE

Ultrasound information entropy imaging is an emerging quantitative ultrasound technique for characterizing local tissue scatterer concentrations and arrangements. However, the commonly used ultrasound Shannon entropy imaging based on histogram-derived discrete probability estimation suffers from the drawbacks of histogram settings dependence and unknown estimator performance. In this paper, we introduced the information-theoretic cumulative residual entropy (CRE) defined in a continuous distribution of cumulative distribution functions as a new entropy measure of ultrasound backscatter envelope uncertainty or complexity, and proposed ultrasound CRE imaging for tissue characterization.

METHODS

We theoretically analyzed the CRE for Rayleigh and Nakagami distributions and proposed a normalized CRE for characterizing scatterer distribution patterns. We proposed a method based on an empirical cumulative distribution function estimator and a trapezoidal numerical integration for estimating the normalized CRE from ultrasound backscatter envelope signals. We presented an ultrasound normalized CRE imaging scheme based on the normalized CRE estimator and the parallel computation technique. We also conducted theoretical analysis of the differential entropy which is an extension of the Shannon entropy to a continuous distribution, and introduced a method for ultrasound differential entropy estimation and imaging. Monte-Carlo simulation experiments were performed to evaluate the estimation accuracy of the normalized CRE and differential entropy estimators. Phantom simulation and clinical experiments were conducted to evaluate the performance of the proposed normalized CRE imaging in characterizing scatterer concentrations and hepatic steatosis (n = 204), respectively.

RESULTS

The theoretical normalized CRE for the Rayleigh distribution was π/4, corresponding to the case where there were ≥10 randomly distributed scatterers within the resolution cell of an ultrasound transducer. The theoretical normalized CRE for the Nakagami distribution decreased as the Nakagami parameter m increased, corresponding to that the ultrasound backscattered statistics varied from pre-Rayleigh to Rayleigh and to post-Rayleigh distributions. Monte-Carlo simulation experiments showed that the proposed normalized CRE and differential entropy estimators can produce a satisfying estimation accuracy even when the size of the test samples is small. Phantom simulation experiments showed that the proposed normalized CRE and differential entropy imaging can characterize scatterer concentrations. Clinical experiments showed that the proposed ultrasound normalized CRE imaging is capable to quantitatively characterize hepatic steatosis, outperforming ultrasound differential entropy imaging and being comparable to ultrasound Shannon entropy and Nakagami imaging.

CONCLUSION

This study sheds light on the theory and methodology of ultrasound normalized CRE. The proposed ultrasound normalized CRE can serve as a new, flexible quantitative ultrasound envelope statistics parameter. The proposed ultrasound normalized CRE imaging may find applications in quantified characterization of biological tissues. Our code will be made available publicly at https://github.com/zhouzhuhuang.

Multimodality quantitative ultrasound envelope statistics imaging based support vector machines for characterizing tissue scatterer distribution patterns: Methods and application in detecting microwave-induced thermal lesions

Ultrasound envelope statistics imaging, including ultrasound Nakagami imaging, homodyned-K imaging, and information entropy imaging, is an important group of quantitative ultrasound techniques for characterizing tissue scatterer distribution patterns, such as scatterer concentrations and arrangements. In this study, we proposed a machine learning approach to integrate the strength of multimodality quantitative ultrasound envelope statistics imaging techniques and applied it to detecting microwave ablation induced thermal lesions in porcine liver ex vivo. The quantitative ultrasound parameters included were homodyned-K α which is a scatterer clustering parameter related to the effective scatterer number per resolution cell, Nakagami m which is a shape parameter of the envelope probability density function, and Shannon entropy which is a measure of signal uncertainty or complexity. Specifically, the homodyned-K log10(α), Nakagami-m, and horizontally normalized Shannon entropy parameters were combined as input features to train a support vector machine (SVM) model to classify thermal lesions with higher scatterer concentrations from normal tissues with lower scatterer concentrations. Through heterogeneous phantom simulations based on Field II, the proposed SVM model showed a classification accuracy above 0.90; the area accuracy and Dice score of higher-scatterer-concentration zone identification exceeded 83% and 0.86, respectively, with the Hausdorff distance <26. Microwave ablation experiments of porcine liver ex vivo at 60-80 W, 1-3 min showed that the SVM model achieved a classification accuracy of 0.85; compared with single log10(α),m, or hNSE parametric imaging, the SVM model achieved the highest area accuracy (89.1%) and Dice score (0.77) as well as the smallest Hausdorff distance (46.38) of coagulation zone identification. We concluded that the proposed multimodality quantitative ultrasound envelope statistics imaging based SVM approach can enhance the capability to characterize tissue scatterer distribution patterns and has the potential to detect the thermal lesions induced by microwave ablation.

A magnetic phantom technique for investigating structural effects on quantitative ultrasound parameters

Media that contain ultrasound scatterers arranged in a regular spatial distribution can be considered as structured. Structural effects affect quantitative ultrasound parameters that reflect the microstructure properties. Prior studies examined structural effects using simulations or phantoms with fixed microarchitecture, focusing on a limited set of ultrasound parameters, with limited attention given to their underlying physical significance. This study aims to investigate the concordance of the physical interpretations of multiple quantitative ultrasound parameters experimentally by introducing a phantom type with an adjustable microarchitecture. The phantom consists of an aqueous solution containing superparamagnetic microspheres, acting as scatterers. The spatial arrangement of the magnetic particles is modified by applying an external magnetic field, therefore changing the degree of structure of the phantom. Quantitative ultrasound parameters are estimated in three different configurations: the magnetic field intensity is varied over time, strength, and orientation. In each experiment, the backscatter coefficient and the envelope quantitative ultrasound parameters are successfully extracted (R2 ≈ 0.94). Their physical interpretations are supported by microphotographs and geometrical considerations through concordant hypotheses. This study paves the way for the use of magnetic phantoms. This methodology could be followed to validate theoretical scattering models and the physical meanings of quantitative ultrasound parameters.

Enhanced identification of women at risk for preterm birth via quantitative ultrasound: a prospective cohort study

BACKGROUND

Historically, clinicians have relied on medical risk factors and clinical symptoms for preterm birth risk assessment. In nulliparous women, clinicians may rely solely on reported symptoms to assess for the risk of preterm birth. The routine use of ultrasound during pregnancy offers the opportunity to incorporate quantitative ultrasound scanning of the cervix to potentially improve assessment of preterm birth risk.

OBJECTIVE

This study aimed to investigate the efficiency of quantitative ultrasound measurements at relatively early stages of pregnancy to enhance identification of women who might be at risk for spontaneous preterm birth.

STUDY DESIGN

A prospective cohort study of pregnant women was conducted with volunteer participants receiving care from the University of Illinois Hospital in Chicago, Illinois. Participants received a standard clinical screening followed by 2 research screenings conducted at 20±2 and 24±2 weeks. Quantitative ultrasound scans were performed during research screenings by registered diagnostic medical sonographers using a standard cervical length approach. Quantitative ultrasound features were computed from calibrated raw radiofrequency backscattered signals. Full-term birth outcomes and spontaneous preterm birth outcomes were included in the analysis. Medically indicated preterm births were excluded from the analysis. Using data from each visit, logistic regression with Akaike information criterion feature selection was conducted to derive predictive models for each time frame based on historical clinical and quantitative ultrasound features. Model evaluations included a likelihood ratio test of quantitative ultrasound features, cross-validated receiver operating characteristic curve analysis, sensitivity, and specificity.

RESULTS

On the basis of historical clinical features alone, the best predictive model had an estimated receiver operating characteristic area under the curve of 0.56±0.03. By the time frame of Visit 1, a predictive model using both historical clinical and quantitative ultrasound features provided a modest improvement in the area under the curve (0.63±0.03) relative to that of the predictive model using only historical clinical features. By the time frame of Visit 2, the predictive model using historical clinical and quantitative ultrasound features provided significant improvement (likelihood ratio test, P<.01), with an area under the curve of 0.69±0.03.

CONCLUSION

Accurate identification of women at risk for spontaneous preterm birth solely through historical clinical features has been proven to be difficult. In this study, a history of preterm birth was the most significant historical clinical predictor of preterm birth risk, but the historical clinical predictive model performance was not statistically significantly better than the no-skill level. According to our study results, including quantitative ultrasound yields a statistically significant improvement in risk prediction as the pregnancy progresses.

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.

A Highly Flexible Piezoelectric Ultrasonic Sensor for Wearable Bone Density Testing

Driven by the loss of bone calcium, the elderly are prone to osteoporosis, and regular routine checks on bone status are necessary, which mainly rely on bone testing equipment. Therefore, wearable real-time healthcare devices have become a research hotspot. Herein, we designed a high-performance flexible ultrasonic bone testing system using axial transmission technology based on quantitative ultrasound theory. First, a new rare-earth-element-doped PMN-PZT piezoelectric ceramic was synthesized using a solid-state reaction, and characterized by X-ray diffraction and SEM. Both a high piezoelectric coefficient d33 = 525 pC/N and electromechanical coupling factors of k33 = 0.77, kt = 0.58 and kp = 0.63 were achieved in 1%La/Sm-doped 0.17 PMN-0.47 PZ-0.36 PT ceramics. Combining a flexible PDMS substrate with an ultrasonic array, a flexible hardware circuit was designed which includes a pulse excitation module, ultrasound array module, amplification module, filter module, digital-to-analog conversion module and wireless transmission module, showing high power transfer efficiency and power intensity with values of 35% and 55.4 mW/cm2, respectively. Finally, the humerus, femur and fibula were examined by the flexible device attached to the skin, and the bone condition was displayed in real time on the mobile client, which indicates the potential clinical application of this device in the field of wearable healthcare.

Predicting Spontaneous Pre-term Birth Risk Is Improved When Quantitative Ultrasound Data Are Included With Historical Clinical Data

OBJECTIVE

Predicting women at risk for spontaneous pre-term birth (sPTB) has been medically challenging because of the lack of signs and symptoms of pre-term birth until interventions are too late. We hypothesized that prediction of the sPTB risk level is enhanced when using both historical clinical (HC) data and quantitative ultrasound (QUS) data compared with using only HC data. HC data defined herein included birth history prior to that of the current pregnancy as well as, from the current pregnancy, a clinical cervical length assessment and physical examination data.

METHODS

The study population included 248 full-term births (FTBs) and 26 sPTBs. QUS scans (Siemens S2000 and MC9-4) were performed by registered diagnostic medical sonographers using a standard cervical length approach. Two cervical QUS scans were conducted at 20 ± 2 and 24 ± 2 wk of gestation. Multiple QUS features were evaluated from calibrated raw radiofrequency backscattered ultrasonic signals. Two statistical models designed to determine sPTB risk were compared: (i) HC data alone and (ii) combined HC and QUS data. Model comparisons included a likelihood ratio test, cross-validated receiver operating characteristic area under the curve, sensitivity and specificity. The study's birth outcomes were only FTBs and sPTBs; medically induced pre-term births were not included.

DISCUSSION

Combined HC and QUS data identified women at risk of sPTB with better AUC (0.68, 95% confidence interval [CI]: 0.57-0.78) compared with HC data alone (0.53, 95% CI: 0.40-0.66) and HC data + cervical length at 18-20 wk of gestation (average AUC = 0.51, 95% CI: 0.38-0.64). A likelihood ratio test for significance of QUS features in the classification model was highly statistically significant (p < 0.01).

CONCLUSION

Even with only 26 sPTBs among 274 births, value was added in predicting sPTB when QUS data were included with HC data.

Multiplexed Ultrasound Imaging Using Spectral Analysis on Gas Vesicles

Current advances in ultrasound imaging techniques combined with the next generation contrast agents such as gas vesicles (GV) revolutionize the visualization of biological tissues with spatiotemporal precision. In optics, fluorescent proteins enable understanding of molecular and cellular functions in biological systems due to their multiplexed imaging capability. Here, a panel of GVs is investigated using mid-band fit (MBF) spectral imaging to realize multiplexed ultrasound imaging to uniquely visualize locations of different types of stationary GVs. The MBF spectral imaging technique demonstrates that stationary clustered GVs are efficiently localized and distinguished from unclustered GVs in agarose gel phantom and 3D vessel structures are visualized in ex vivo mouse liver specimens. Mouse macrophages serve as carriers of clustered and unclustered GVs and multiplexing beacons to report cells' spatial locations by emitting distinct spectral signals. 2D MBF spectral images are reconstructed, and pixels in these images are classified depending on MBF values by comparing predetermined filters that predict the existence of cells with clustered and unclustered GVs. This pseudo-coloring scheme clearly distinguishes the locations of two classes of cells like pseudo-color images in fluorescence microscopy.

Ultrasound Backscatter Envelope Statistics Parametric Imaging for Liver Fibrosis Characterization: A Review

Early detection and diagnosis of liver fibrosis is of critical importance. Currently the gold standard for diagnosing liver fibrosis is biopsy. However, liver biopsy is invasive and associated with sampling errors and can lead to complications such as bleeding. Therefore, developing noninvasive imaging techniques for assessing liver fibrosis is of clinical value. Ultrasound has become the first-line tool for the management of chronic liver diseases. However, the commonly used B-mode ultrasound is qualitative and can cause interobserver or intraobserver difference. Ultrasound backscatter envelope statistics parametric imaging is an important group of quantitative ultrasound techniques that have been applied to characterizing different kinds of tissue. However, a state-of-the-art review of ultrasound backscatter envelope statistics parametric imaging for liver fibrosis characterization has not been conducted. In this paper, we focused on the development of ultrasound backscatter envelope statistics parametric imaging techniques for assessing liver fibrosis from 1998 to September 2019. We classified these techniques into six categories: constant false alarm rate, fiber structure extraction technique, acoustic structure quantification, quantile-quantile probability plot, the multi-Rayleigh model, and the Nakagami model. We presented the theoretical background and algorithms for liver fibrosis assessment by ultrasound backscatter envelope statistics parametric imaging. Then, the specific applications of ultrasound backscatter envelope statistics parametric imaging techniques to liver fibrosis evaluation were reviewed and analyzed. Finally, the pros and cons of each technique were discussed, and the future development was suggested.

Assessment of Hepatic Steatosis in Nonalcoholic Fatty Liver Disease by Using Quantitative US

Background Advanced confounder-corrected chemical shift-encoded MRI-derived proton density fat fraction (PDFF) is a leading parameter for fat fraction quantification in nonalcoholic fatty liver disease (NAFLD). Because of the limited availability of this MRI technique, there is a need to develop and validate alternative parameters to assess liver fat. Purpose To assess relationship of quantitative US parameters to MRI PDFF and to develop multivariable quantitative US models to detect hepatic steatosis and quantify hepatic fat. Materials and Methods Adults with known NAFLD or who were suspected of having NAFLD were prospectively recruited between August 2015 and February 2019. Participants underwent quantitative US and chemical shift-encoded MRI liver examinations. Liver biopsies were performed if clinically indicated. The correlation between seven quantitative US parameters and MRI PDFF was evaluated. By using leave-one-out cross validation, two quantitative US multivariable models were evaluated: a classifier to differentiate participants with NAFLD versus participants without NAFLD and a fat fraction estimator. Classifier performance was summarized by area under the receiver operating characteristic curve and area under the precision-recall curve. Fat fraction estimator performance was evaluated by correlation, linearity, and bias. Results Included were 102 participants (mean age, 52 years ± 13 [standard deviation]; 53 women), 78 with NAFLD (MRI PDFF ≥ 5%). A two-variable classifier yielded a cross-validated area under the receiver operating characteristic curve of 0.89 (95% confidence interval: 0.82, 0.96) and an area under the precision-recall curve of 0.96 (95% confidence interval: 0.93, 0.99). The cross-validated fat fraction predicted by a two-variable fat fraction estimator was correlated with MRI PDFF (Spearman ρ = 0.82 [P < .001]; Pearson r = 0.76 [P < .001]). The mean bias was 0.02% (P = .97), and 95% limits of agreement were ±12.0%. The predicted fat fraction was linear with MRI PDFF (R ² = 0.63; slope, 0.69; intercept, 4.3%) for MRI PDFF of 34% or less. Conclusion A multivariable quantitative US approach yielded excellent correlation with MRI proton density fat fraction for hepatic steatosis assessment in nonalcoholic fatty liver disease. © RSNA, 2020 Online supplemental material is available for this article.