Attenuation map or measurements based on local attenuation coefficient slope (ACS) in quantitative ultrasound (QUS) has shown potential for diagnosis of liver steatosis. In liver cancers, tissue ...abnormalities and tumors detected using ACS are also of interest to provide new image contrast to clinicians. Current phantom-based approaches have the limitation of assuming comparable speed of sound between the reference phantom and insonified tissues. Moreover, these methods present the inconvenience for operators to acquire data on phantoms as well as on patients. The main goal was to alleviate these drawbacks by proposing a methodology for constructing phantom-free regularized (PF-R) local ACS maps and investigate the performance in both homogeneous and heterogeneous media. The proposed method was tested on two tissue mimicking media with different ACS constructed as homogeneous phantoms, side-by-side and top-to-bottom phantoms, and inclusion phantoms with different attenuations. Moreover, an in-vivo proof-of-concept was performed on healthy, steatotic and cancerous human liver datasets. Modifications brought to previous works include: a) a linear interpolation of the power spectrum in log-scale; b) the relaxation of the underlying hypothesis on the diffraction factor; c) a generalization to nonhomogeneous local ACS; and d) an adaptive restriction of frequencies to a more reliable range than the usable frequency range. Regularization was formulated as a generalized LASSO, and a variant of the Bayesian Information Criterion (BIC) was applied to estimate the Lagrangian multiplier on the LASSO constraint. In addition, we evaluated the proposed algorithm when applying median filtering before and after regularization. Tests conducted showed that the PF-R yielded robust results in all tested conditions, suggesting potential for additional validation as a diagnosis method.
To develop a quantitative ultrasound (QUS)- and elastography-based model to improve classification of steatosis grade, inflammation grade, and fibrosis stage in patients with chronic liver disease in ...comparison with shear wave elastography alone, using histopathology as the reference standard.
This ancillary study to a prospective institutional review-board approved study included 82 patients with non-alcoholic fatty liver disease, chronic hepatitis B or C virus, or autoimmune hepatitis. Elastography measurements, homodyned K-distribution parametric maps, and total attenuation coefficient slope were recorded. Random forests classification and bootstrapping were used to identify combinations of parameters that provided the highest diagnostic accuracy. Receiver operating characteristic (ROC) curves were computed.
For classification of steatosis grade S0 vs. S1-3, S0-1 vs. S2-3, S0-2 vs. S3, area under the receiver operating characteristic curve (AUC) were respectively 0.60, 0.63, and 0.62 with elasticity alone, and 0.90, 0.81, and 0.78 with the best tested model combining QUS and elastography features. For classification of inflammation grade A0 vs. A1-3, A0-1 vs. A2-3, A0-2 vs. A3, AUCs were respectively 0.56, 0.62, and 0.64 with elasticity alone, and 0.75, 0.68, and 0.69 with the best model. For classification of liver fibrosis stage F0 vs. F1-4, F0-1 vs. F2-4, F0-2 vs. F3-4, F0-3 vs. F4, AUCs were respectively 0.66, 0.77, 0.72, and 0.74 with elasticity alone, and 0.72, 0.77, 0.77, and 0.75 with the best model.
Random forest models incorporating QUS and shear wave elastography increased the classification accuracy of liver steatosis, inflammation, and fibrosis when compared to shear wave elastography alone.
Quantitative ultrasound (QUS) techniques using radiofrequency (RF) backscattered signals have been used for tissue characterization of numerous organ systems. One approach is to use the magnitude and ...frequency dependence of backscatter echoes to quantify tissue structures. Another approach is to use first-order statistical properties of the echo envelope as a signature of the tissue microstructure. We propose a unification of these QUS concepts. For this purpose, a mixture of homodyned K-distributions is introduced to model the echo envelope, together with an estimation method and a physical interpretation of its parameters based on the echo signal spectrum. In particular, the total, coherent and diffuse signal powers related to the proposed mixture model are expressed explicitly in terms of the structure factor previously studied to describe the backscatter coefficient (BSC). Then, this approach is illustrated in the context of red blood cell (RBC) aggregation. It is experimentally shown that the total, coherent and diffuse signal powers are determined by a structural parameter of the spectral Structure Factor Size and Attenuation Estimator. A two-way repeated measures ANOVA test showed that attenuation ( p-value of 0.077) and attenuation compensation ( p-value of 0.527) had no significant effect on the diffuse to total power ratio. These results constitute a further step in understanding the physical meaning of first-order statistics of ultrasound images and their relations to QUS techniques. The proposed unifying concepts should be applicable to other biological tissues than blood considering that the structure factor can theoretically model any spatial distribution of scatterers.
The homodyned K-distribution appears naturally in the context of random walks and provides a useful model for the distribution of the received intensity in a wide range of non-Gaussian scattering ...configurations, including medical ultrasonics. An estimation method for the homodyned K-distribution based on the first moment of the intensity and two log-moments (XU method), namely the
and
-statistics previously studied in the special case of the K-distribution, is proposed as an alternative to a method based on the first three moments of the intensity (MI method) or the amplitude (MA method), and a method based on the signal-to-noise ratio (SNR), the skewness and the kurtosis of two fractional orders of the amplitude (labeled RSK method). Properties of the
and
statistics for the homodyned K-distribution are proved, except for one conjecture. Using those properties, an algorithm based on the bisection method for monotonous functions was developed. The algorithm has a geometric rate of convergence. Various tests were performed to study the behavior of the estimators. It was shown with simulated data samples that the estimations of the parameters 1/
and 1/(
+ 1) of the homodyned K-distribution are preferable to the direct estimations of the clustering parameter
and the structure parameter
(with respective relative root mean squared errors (RMSEs) of 0.63 and 0.13 as opposed to 1.04 and 4.37, when
= 1000). Tests on simulated ultrasound images with only diffuse scatterers (up to 10 per resolution cell) indicated that the XU estimator is overall more reliable than the other three estimators for the estimation of 1/
, with relative RMSEs of 0.79 (MI), 0.61 (MA), 0.53 (XU) and 0.67 (RSK). For the parameter 1/(
+ 1), the relative RMSEs were equal to 0.074 (MI), 0.075 (MA), 0.069 (XU) and 0.100 (RSK). In the case of a large number of scatterers (11 to 20 per resolution cell), the relative RMSEs of 1/
were equal to 1.43 (MI), 1.27 (MA), 1.25 (XU) and 1.33 (RSK), and the relative RMSEs of 1/(
+ 1) were equal to 0.14 (MI), 0.16 (MA), 0.17 (XU) and 0.20 (RSK). The four methods were also tested on simulated ultrasound images with a variable density of periodic scatterers to test images with a coherent component. The addition of noise on ultrasound images was also studied. Results showed that the XU estimator was overall better than the three other ones. Finally, on the simulated ultrasound images, the average computation times per image were equal to 6.0 ms (MI), 8.0 ms (MA), 6.8 ms (XU) and 500 ms (RSK). Thus, a fast, reliable, and novel algorithm for the estimation of the homodyned K-distribution was proposed.
Change in viscoelastic properties of biological tissues may often be symptomatic of dysfunction that can be correlated with tissue pathology. Shear wave (SW) elastography is an imaging method mainly ...used to assess stiffness but with the potential to measure viscoelasticity of biological tissues. This can enable tissue characterization and, thus, can be used as a marker to improve the diagnosis of pathological lesions. In this study, a frequency-shift method-based framework is presented for the reconstruction of viscosity by analyzing the spectral properties of acoustic radiation force-induced SWs. The aim of this study was to investigate the feasibility of viscosity reconstruction maps in homogeneous as well as heterogeneous samples. Experiments were performed in four in vitro phantoms, two ex vivo porcine liver samples, two ex vivo fatty duck liver samples, and one in vivo fatty goose liver. Successful viscosity maps were reconstructed in homogeneous and heterogeneous phantoms with embedded mechanical inclusions having different geometries. Quantitative values of viscosity obtained for two porcine liver tissues, two fatty duck liver samples, and one goose fatty liver were (mean ± SD) 0.61 ± 0.21 and 0.52 ± 0.35, 1.28 ± 0.54 and 1.36 ± 0.73, and 1.67 ± 0.70 Pa.s, respectively.
Objective myocardial contractility assessment during stress tests aims to improve the diagnosis of myocardial ischemia. Tissue Doppler imaging (TDI) or optical flow (OF) speckle tracking ...echocardiography has been used to quantify myocardial contractility at rest. However, this is more challenging during stress tests due to image decorrelation at high heart rates. Moreover, stress tests imply a high frame rate which leads to a limited lateral field-of-view. Therefore, a large lateral field-of-view robust ultrafast myocardial regularized lateral field-of-view robust ultrafast myocardial regularized OF-TDI principal strain estimator has been developed for high-frame-rate echocardiography of coherently compounded transmitted diverging waves. Feasibility and accuracy of the proposed estimator were validated with in vitro (using sonomicrometry as gold standard) and in vivo stress experiments. Compared with OF strain imaging, the proposed estimator improved the accuracy of principal major and minor strains during stress tests, with an average contrast-to-noise ratio improvement of 4.4 ± 2.7 dB ( p -value < 0.01). Moreover, there was a significant correlation and a very close agreement between the proposed estimator and sonomicrometry for tested heart rates between 60 and 180 beats per minute. The averages ± STD of R 2 and biases ± STD between them were 0.96 ± 0.04 ( p -value < 0.01) and 0.01 ± 0.03 % in the axial direction, respectively; and 0.94 ± 0.02 ( p -value < 0.01) and 0.04 ± 0.06 % in the lateral direction, respectively. These results suggest that the proposed estimator could be useful clinically to provide an accurate and quantitative 2D large lateral field-of-view myocardial strain assessment at high heart rates during stress echocardiography.
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.
An enhanced inflammatory response is a trigger to the production of blood macromolecules involved in abnormally high levels of erythrocyte aggregation.
This study aimed at demonstrating for the first ...time the clinical feasibility of a non-invasive ultrasound-based erythrocyte aggregation quantitative measurement method for potential application in critical care medicine.
Erythrocyte aggregation was evaluated using modeling of the backscatter coefficient with the Structure Factor Size and Attenuation Estimator (SFSAE). SFSAE spectral parameters W (packing factor) and D (mean aggregate diameter) were measured within the antebrachial vein of the forearm and tibial vein of the leg in 50 healthy participants at natural flow and reduced flow controlled by a pressurized bracelet. Blood samples were also collected to measure erythrocyte aggregation ex vivo with an erythroaggregometer (parameter S10).
W and Din vivo measurements were positively correlated with the ex vivoS10 index for both measurement sites and shear rates (correlations between 0.35-0.81, p < 0.05). Measurement at low shear rate was found to increase the sensitivity and reliability of this non-invasive measurement method.
We behold that the SFSAE method presents systemic measures of the erythrocyte aggregation level, since results on upper and lower limbs were highly correlated.
In ultrasound imaging, various statistical distributions have been proposed to model the first-order statistics of the amplitude of the echo envelope. We present an overview of these distributions ...based on their compound representation, which comprises three aspects: the modulated distribution (Rice or Nakagami); the modulating distribution (gamma, inverse Gaussian or even generalized inverse Gaussian); and the modulated parameters (the diffuse signal power with or without the coherent signal component or the coherent signal power). This unifying point of view makes the comparison of the various models conceptually easier. In particular, we discuss the implications of the modulated parameters on the mean intensity and the signal-to-noise ratio of the intensity in the case of a vanishing diffuse signal. We conclude that the homodyned K-distribution is the only model among the literature for which the parameters have a physical meaning that is consistent with the limiting case, although the other distributions may fit real data.
The Horn and Schunck (HS) method, which amounts to the Jacobi iterative scheme in the interior of the image, was one of the first optical flow algorithms. In this article, we prove the convergence of ...the HS method, whenever the problem is well-posed. Our result is shown in the framework of a generalization of the HS method in dimension
≥ 1, with a broad definition of the discrete Laplacian. In this context, the condition for the convergence is that the intensity gradients are not all contained in a same hyperplane. Two other articles (17 and 13) claimed to solve this problem in the case
= 2, but it appears that both of these proofs are erroneous. Moreover, we explain why some standard results on the convergence of the Jacobi method do not apply for the HS problem, unless
= 1. It is also shown that the convergence of the HS scheme implies the convergence of the Gauss-Seidel and SOR schemes for the HS problem.