A point process approach to assess the frequency dependence of ultrasound backscattering by aggregating red blood cells

Contrast analysis in ultrafast ultrasound blood flow imaging of jugular vein

PURPOSE

The contrasts of flowing blood in in vitro experiments using porcine blood and in vivo measurements of human jugular veins were analyzed to demonstrate that the hemorheological property was dependent on the shear rate.

METHODS

Blood samples (45% hematocrit) suspended in saline or plasma were compared with examine the difference in viscoelasticity. Ultrafast plane-wave imaging at an ultrasonic center frequency of 7.5 MHz was performed on different steady flows in a graphite-agar phantom. Also, in vivo measurement was performed in young, healthy subjects and patients with diabetes. A spatiotemporal matrix of beamformed radio-frequency data was used for the singular value decomposition (SVD) clutter filter. The clutter-filtered B-mode image was calculated as the amplitude envelope normalized at the first frame in the diastolic phase to evaluate contrast. The shear rate was estimated as the velocity gradient perpendicular to the lateral axis.

RESULTS

Although nonaggregated erythrocytes at a high shear rate exhibited a low echogenicity, the echogenicity in the plasma sample overall increased due to erythrocyte aggregation at a low shear rate. In addition, the frequency of detection of specular components, defined as components beyond twice the standard deviation of a contrast map obtained from a clutter-filtered B-mode image, increased in the porcine blood at a high shear rate and the venous blood in healthy subjects versus patients with diabetes.

CONCLUSION

The possibility of characterizing hemorheological properties dependent on the shear rate and diabetes condition was indicated using ultrafast plane-wave imaging with an SVD-based clutter filter.

Effect of Clutter Filter in High-Frame-Rate Ultrasonic Backscatter Coefficient Analysis

High-frame-rate imaging with a clutter filter can clearly visualize blood flow signals and provide more efficient discrimination with tissue signals. In vitro studies using clutter-less phantom and high-frequency ultrasound suggested a possibility of evaluating the red blood cell (RBC) aggregation by analyzing the frequency dependence of the backscatter coefficient (BSC). However, in in vivo applications, clutter filtering is required to visualize echoes from the RBC. This study initially evaluated the effect of the clutter filter for ultrasonic BSC analysis for in vitro and preliminary in vivo data to characterize hemorheology. Coherently compounded plane wave imaging at a frame rate of 2 kHz was carried out in high-frame-rate imaging. Two samples of RBCs suspended by saline and autologous plasma for in vitro data were circulated in two types of flow phantoms without or with clutter signals. The singular value decomposition was applied to suppress the clutter signal in the flow phantom. The BSC was calculated using the reference phantom method, and it was parametrized by spectral slope and mid-band fit (MBF) between 4-12 MHz. The velocity distribution was estimated by the block matching method, and the shear rate was estimated by the least squares approximation of the slope near the wall. Consequently, the spectral slope of the saline sample was always around four (Rayleigh scattering), independently of the shear rate, because the RBCs did not aggregate in the solution. Conversely, the spectral slope of the plasma sample was lower than four at low shear rates but approached four by increasing the shear rate, because the aggregations were presumably dissolved by the high shear rate. Moreover, the MBF of the plasma sample decreased from -36 to -49 dB in both flow phantoms with increasing shear rates, from approximately 10 to 100 s^-1. The variation in the spectral slope and MBF in the saline sample was comparable to the results of in vivo cases in healthy human jugular veins when the tissue and blood flow signals could be separated.

Pulse-Echo Technique to Compensate for Laminate Membrane Transmission Loss in Phantom-Based Ultrasonic Attenuation Coefficient Measurements

OBJECTIVES

Accurately measuring the attenuation coefficient (AC) of reference phantoms is critical in clinical applications of quantitative ultrasound. Phantom AC measurement requires proper compensation of membrane transmission loss. Conventional methods require separate membrane samples to obtain membrane transmission loss. Unfortunately, separate membrane samples are often unavailable. A pulse-echo approach is proposed herein to compensate for membrane transmission loss without requiring separate membrane samples.

METHODS

The proposed method consists of the following steps. First, the insertion loss, caused by phantom attenuation and membrane transmission loss, is measured. Second, the membrane reflection coefficient is measured. Third, the unknown acoustic parameters of the membrane and phantom material are estimated by fitting theoretical reflection coefficient to the measured one. Finally, the fitted parameters are used to estimate membrane transmission loss and phantom AC. The proposed method was validated through k-Wave simulations and phantom experiments. Experimental AC measurements were repeated on 5 distinct phantoms by 2 operators to assess the repeatability and reproducibility of the proposed method. Five transducers were used to cover a broad bandwidth (0.7-16 MHz).

RESULTS

The acquired AC in the simulations had a maximum error of 0.06 dB/cm-MHz for simulated phantom AC values ranging from 0.5 to 1 dB/cm-MHz. The acquired AC in the experiments had a maximum error of 0.045 dB/cm-MHz for phantom AC values ranging from 0.28 to 1.48 dB/cm-MHz. Good repeatability and cross-operator reproducibility were observed with a mean coefficient of variation below 0.054.

CONCLUSION

The proposed method simplifies phantom AC measurement while providing satisfactory accuracy and precision.

Extracting Quantitative Ultrasonic Parameters from the Backscatter Coefficient

The ultrasonic backscatter coefficient (BSC) is a fundamental quantitative ultrasound (QUS) parameter that contains rich information about the underlying tissue. Deriving parameters from the BSC is essential for fully utilizing the information contained in BSC for tissue characterization. In this chapter, we review two primary approaches for extracting parameters from the BSC versus frequency curve: the model-based approach and the model-free approach, focusing on the model-based approach, where a scattering model is fit to the observed BSC to yield model parameters. For this approach, we will attempt to unite commonly used models under a coherent theoretical framework. We will focus on the underlying assumptions and conditions for various BSC models. Computer code is provided to facilitate the use of some of the models. The strengths and weaknesses of various models are also discussed.

Power laws prevail in medical ultrasound

Major topics in medical ultrasound rest on the physics of wave propagation through tissue. These include fundamental treatments of backscatter, speed of sound, attenuation, and speckle formation. Each topic has developed its own rich history, lexicography, and particular treatments. However, there is ample evidence to suggest that power law relations are operating at a fundamental level in all the basic phenomena related to medical ultrasound. This review paper develops, from literature over the past 60 years, the accumulating theoretical basis and experimental evidence that point to power law behaviors underlying the most important tissue-wave interactions in ultrasound and in shear waves which are now employed in elastography. The common framework of power laws can be useful as a coherent overview of topics, and as a means for improved tissue characterization.

Ultrasound Scattering From Cell-Pellet Biophantoms and Ex Vivo Tumors Provides Insight Into the Cellular Structure Involved in Scattering

The histologically identifiable cellular structure(s) involved in ultrasonic scattering is(are) yet to be uniquely identified. The study quantifies six possible cellular scattering parameters, namely, cell and nucleus radii and their respective cell and nucleus volume fractions as well as a combination of cell and nucleus radii and their volume fraction. The six cellular parameters are each derived from four cell lines (4T1, JC, LMTK, and MAT) and two tissue types (cell-pellet biophantom and ex vivo tumor). Optical histology and quantitative ultrasound (QUS), both independent approaches, are used to yield these cellular parameters. QUS scatterer parameters are experimentally determined using two ultrasonic scattering models: the spherical Gaussian model (GM) and the structure factor model (SFM) to yield insight about scattering from nuclei only and cells only. GM is a classical ultrasonic scattering model to evaluate QUS parameters and is well adapted for diluted media. SFM is adapted for dense media to estimate reasonably well scatterer parameters of cellular structures from ex vivo tissue. Nucleus and cell radii and volume fractions are measured optically from histology. They were used as inputs to calculate BSC for scattering from cells, nuclei, and both cells and nuclei. The QUS-derived scatterers (radii and volume fractions) distributions were then compared to the optical histology scatterer parameters derived from these calculated BSCs. The results suggest scattering from cells only (LMTK and MAT) or both cells and nuclei (4T1 and JC) for cell-pellet biophantoms and scattering from nuclei only for tumors.

Ultrasound simulation of blood with different red blood cell aggregations and concentrations

BACKGROUND

Considerable progress of ultrasound simulation on blood has enhanced the characterizing of red blood cell (RBC) aggregation.

OBJECTIVE

A novel simulation method aims at modeling the blood with different RBC aggregations and concentrations is proposed.

METHODS

The modeling process is as follows: (i) A three-dimensional scatterer model is first built by a mapping with a Hilbert space-filling curve from the one-dimensional scatterer distribution. (ii) To illustrate the relationship between the model parameters and the RBC aggregation level, a variety of blood samples are prepared and scanned to acquire their radiofrequency signals in-vitro. (iii) The model parameters are determined by matching the Nakagami-distribution characteristics of envelope signals simulated from the model with those measured from the blood samples.

RESULTS

Nakagami metrics m estimated from 15 kinds of blood samples (hematocrits of 20%, 40%, 60% and plasma concentrations of 15%, 30%, 45%, 60%, 75%) are compared with metrics estimated by their corresponding models (each with different eligible parameters). Results show that for the three hematocrit levels, the mean and standard deviation of the root-mean-squared deviations of m are 0.27 ± 0.0026, 0.16 ± 0.0021, 0.12 ± 0.0018 respectively.

CONCLUSION

The proposed simulation model provides a viable data source to evaluate the performance of the ultrasound-based methods for quantifying RBC aggregation.

Statistical modeling of ultrasound signals related to the packing factor of wave scattering phenomena for structural characterization

The two-dimensional homodyned K-distribution has been widely used to model the echo envelope of ultrasound radio frequency (RF) signals in the field of medical ultrasonics. The main contribution of this work is to present a theoretical framework for supporting this model of the echo envelope and statistical models of the RF signals and their Hilbert transform in the case in which the scatterers' positions may be dependent. In doing so, the law of large numbers, Lyapounov's central limit theorem, and the Berry-Esseen theorem are being used. In particular, the proposed theoretical framework supports a previous conjecture relating the scatterer clustering parameter of the homodyned K-distribution to the packing factor W, which is related to the spatial organization of the scatterers, appearing in statistical physics or backscatter coefficient modeling. Simulations showed that the proposed modeling is valid for a number of scatterers and packing factors varying by steps of 2 from 1 to 21 and 1 to 11, respectively. The proposed framework allows, in principle, the detection of the structural information taking place at a scale smaller than the wavelength based solely on the statistical analysis of the RF signals or their echo envelope, although this goal was previously achieved based on the spectral analysis of ultrasound signals.

Classification of red blood cell aggregation using empirical wavelet transform analysis of ultrasonic radiofrequency echo signals

Grading red blood cell (RBC) aggregation is important for the early diagnosis and prevention of related diseases such as ischemic cardio-cerebrovascular disease, type II diabetes, deep vein thrombosis, and sickle cell disease. In this study, a machine learning technique based on an adaptive analysis of ultrasonic radiofrequency (RF) echo signals in blood is proposed, and its feasibility for classifying RBC aggregation is explored. Using an adaptive empirical wavelet transform (EWT) analysis, the ultrasonic RF signals are decomposed into a series of empirical mode functions (EMFs); then, dominant empirical mode functions (DEMFs) are selected from the series. Six statistical characteristics, including the mean, variance, median, kurtosis, root mean square (RMS), and skewness are calculated for the locally normalized DEMFs, aiming to form primary feature vectors. Random forest (RDF) and support vector machine (SVM) classifiers are trained with the given feature vectors to obtain prediction models for RBC classification. Ultrasonic RF echo signals are acquired from five groups of six types of porcine blood samples with average numbers of aggregated RBCs of 1.04, 1.20, 1.83, 2.31, 2.72, and 4.28, respectively, to test the classification performance of the proposed method. The best subset with regard to the variance, kurtosis, and RMS is determined according to the maximum accuracy based on the RDF and SVM classifiers. The classification accuracies are 84.03 ± 3.13% for the RDF classifier, and 85.88 ± 2.99% for the SVM classifier. The mean classification accuracy of the SVM classifier is 1.85% better than that of the RDF classifier. In conclusion, the machine learning method is useful for the discrimination of varying degrees of RBC aggregation, and has potential for use in characterizing and monitoring the RBC aggregation in vessels.

A study on understanding the physical mechanism of change in ultrasonic envelope statistical property during temperature elevation

PURPOSE

Our previous studies demonstrate that the variation in ultrasonic envelope statistics is correlated with the temperature change inside scattering media. This variation is identified as the change in the scatterer structure during thermal expansion or contraction. However, no specific evidence has been verified to date. This study numerically reproduces the change in the scatterer distribution during thermal expansion or contraction using finite element simulations and also investigates how the situation is altered by different material properties.

METHODS

The material properties of a linear elastic solid depend on the thermal expansion coefficient, thermal conductivity, specific heat, and initial scatterer number density. Three-dimensional displacements, calculated in the simulation, were sequentially used to update the positions of the randomly distributed scatterers. Ultrasound signals from the scatterer distribution were generated by simulating a 7.5-MHz linear array transducer whose specifications were the same as those in the experimental measurements of several phantoms and excised porcine livers. To represent the change in the envelope statistical feature, the absolute value of the ratio change in the logarithmic Nakagami (NA) parameter, Δ m , at each time was calculated as a value normalized with the initial NA parameter.

RESULTS

The change in the scatterer number density relates to the volume change during temperature elevation. The magnitude of the Δ m shift against the temperature change increases depending on the higher thermal expansion coefficient. In contrast, the relationship between Δ m and the scatterer number density is similar with any material property. Additionally, the changes in Δ m obtained by several experimental phantoms with low to high scatterer number densities are comparable with the numerical simulation results.

CONCLUSIONS

The change in Δ m is indirectly related to the change in the scatterer number density owing to the volume change during thermal expansion or contraction.

Quantifying scattering from dense media using two-dimensional impedance maps

A better understanding of ultrasound scattering in a three-dimensional (3D) medium can provide more accurate methods for ultrasound tissue characterization. The possibility of using two-dimensional impedance maps (2DZMs) based on correlation coefficients has shown promise in the case of isotropic and sparse medium [Luchies and Oelze, J. Acoust. Soc. Am. 139, 1557-1564 (2016)]. The present study investigates the use of 2DZMs in order to quantify 3D scatterer properties of dense media from two-dimensional (2D) histological slices. Two 2DZM approaches were studied: one based on the correlation coefficient and the other based on the 2D Fourier transform of 2DZMs. Both 2DZM approaches consist in estimating the backscatter coefficient (BSC) from several 2DZMs, and then the resulting BSC was fit to the theoretical polydisperse structure factor model to yield 3D scatterer properties. Simulation studies were performed to evaluate the ability of both 2DZM approaches to quantify scattering of a 3D medium containing randomly distributed polydisperse spheres or monodisperse ellipsoids. Experimental studies were also performed using the histology photomicrographs obtained from HT29 cell pellet phantoms. Results demonstrate that the 2DZM Fourier transform-based approach was more suitable than the correlation coefficient-based approach for estimating scatterer properties when using a small number of 2DZMs.

Speckle from branching vasculature: dependence on number density

Purpose: Recent theories examine the role of the fractal branching vasculature as a primary site of Born scattering from soft normal tissues. These derivations postulate that the first-order statistics of speckle from soft tissue, such as the liver, thyroid, and prostate, will follow a Burr distribution with a power law parameter that can be related back to the underlying power law, which governs the branching network. However, the issue of scatterer spacing, or the number of cylindrical vessels per sample volume of the interrogating pulse, has not been directly addressed. Approach: Speckle statistics are examined with a 3D simulation that varies the number density and the governing power law parameter of an ensemble of different sized cylinders. Several in vivo liver scans are also analyzed for confirmation across different conditions. Results: The Burr distribution is found to be an appropriate model for the histogram of amplitudes from speckle regions, where the parameters track the underlying power law and scatterer density conditions. These results are also tested in a more general model of rat liver scans in normal versus abnormal conditions, and the resulting Burr parameters are also found to be appropriate and sensitive to underlying scatterer distributions. Conclusions: These preliminary results suggest that the classical Burr distribution may be useful in the quantification of scattering of ultrasound from soft vascularized tissues and as a tool in tissue characterization.

Ultrasonic backscattering and microstructure in sheared concentrated suspensions

Quantitative ultrasound techniques based on the parametrization of the backscatter coefficient (BSC) are used to characterize concentrated particle suspensions. Specifically, a scattering model is fit to the measured BSC and the fit parameters can provide local suspension properties. The scattering models generally assume an isotropic microstructure (i.e., spatial organization) of the scatterers, whereas the sheared concentrated suspensions can develop an anisotropic microstructure. This paper studied the influence of the shear-induced anisotropic microstructure of concentrated suspensions on the ultrasonic backscattering. Experiments were conducted on suspensions of polymethylmetacrylate spheres (5.8 μm in radius) sheared in a Couette flow device to obtain anisotropic microstructure and then mixed by hand to obtain isotropic microstructure. Experimental structure factors that are related to the spatial distribution of sphere positions were obtained by comparing the BSCs of one concentrated and one diluted suspension. Finally, Stokesian dynamics numerical simulations of sheared concentrated suspensions are used to determine the pair correlation function, which is linked to the Fourier transform of the structure factor. The experimental structure factors are found to be in good agreement with numerical simulations. The numerical simulation demonstrates that the angular-dependent BSCs and structure factors are caused by the shear-induced anisotropic microstructure within the suspension.

The first order statistics of backscatter from the fractal branching vasculature

The issue of speckle statistics from ultrasound images of soft tissues such as the liver has a long and rich history. A number of theoretical distributions, some related to random scatterers or fades in optics and radar, have been formulated for pulse-echo interference patterns. This work proposes an alternative framework in which the dominant echoes are presumed to result from Born scattering from fluid-filled vessels that permeate the tissue parenchyma. These are modeled as a branching, fractal, self-similar, multiscale collection of cylindrical scatterers governed by a power law distribution relating to the number of branches at each radius. A deterministic accounting of the echo envelopes across the scales from small to large is undertaken, leading to a closed form theoretical formula for the histogram of the envelope of the echoes. The normalized histogram is found to be related to the classical Burr distribution, with the key power law parameter directly related to that of the number density of vessels vs diameter, frequently reported in the range of 2 to 4. Examples are given from liver scans to demonstrate the applicability of the theory.

Quantitative Measurement and Evaluation of Red Blood Cell Aggregation in Normal Blood Based on a Modified Hanai Equation

The aggregation of red blood cells (RBCs) in normal blood (non-coagulation) has been quantitatively measured by blood pulsatile flow based on multiple-frequency electrical impedance spectroscopy. The relaxation frequencies f_c under static and flowing conditions of blood pulsatile flow are utilized to evaluate the RBC aggregation quantitatively with the consideration of blood flow factors (RBC orientation, deformation, thickness of electrical double layer (EDL)). Both porcine blood and bovine blood are investigated in experiments, for the reason that porcine blood easily forms RBC aggregates, while bovine blood does not. The results show that the relaxation frequencies f_c of porcine blood and bovine blood present opposite performance, which indicates that the proposed relaxation frequency f_c is efficient to measure RBCs aggregation. Furthermore, the modified Hanai equation is proposed to quantitatively calculate the influence of RBCs aggregation on relaxation frequency f_c. The study confirms the feasibility of a high speed, on-line RBC aggregation sensing method in extracorporeal circulation systems.

A Method for Stereological Determination of the Structure Function From Histological Sections of Isotropic Scattering Media

The frequency-dependent ultrasonic backscatter coefficient (BSC) from tissues, a fundamental parameter estimated by quantitative ultrasound (QUS) techniques, contains microstructure information useful for tissue characterization. To extract the microstructure information from the BSC, the tissue under investigation is often modeled as a collection of discrete scatterers embedded in a homogeneous background. From a discrete scatterer point of view, the BSC is dependent on not only the properties of individual scatterers relative to the background but also the scatterer spatial arrangement [described by the structure function (SF)]. Recently, the 2-D SF was computed from histological tissue sections, and was shown to be related to the volumetric SF extracted from QUS measurements. In this paper, a stereological method is proposed to extract the volumetric (3-D) SF from 2-D histological tissue sections. Simulations and experimental cell pellet biophantom studies were conducted to evaluate the proposed method. Simulation results verified the proposed method. Experimental results showed that the volumetric SF extracted using the proposed method had a significantly better agreement with the QUS-extracted SF than did the 2-D SF extracted in the previous study. The proposed stereological approach provides a useful tool for predicting the SF from histology.

High frequency ultrasound imaging and simulations of sea urchin oocytes

High frequency ultrasound backscatter signals from sea urchin oocytes were measured using a 40 MHz transducer and compared to numerical simulations. The Faran scattering model was used to calculate the ultrasound scattered from single oocytes in suspension. The urchin oocytes are non-nucleated with uniform size and biomechanical properties; the backscatter from each cell is similar and easy to simulate, unlike typical nucleated mammalian cells. The time domain signal measured from single oocytes in suspension showed two distinct peaks, and the power spectrum was periodic with minima spaced approximately 10 MHz apart. Good agreement to the Faran scattering model was observed. Measurements from tightly packed oocyte cell pellets showed similar periodic features in the power spectra, which was a result of the uniform size and consistent biomechanical properties of the cells. Numerical simulations that calculated the ultrasound scattered from individual oocytes within a three dimensional volume showed good agreement to the measured signals and B-scan images. A cepstral analysis of the signal was used to calculate the size of the cells, which was 78.7 μm (measured) and 81.4 μm (simulated). This work supports the single scattering approximation, where ultrasound is discretely scattered from single cells within a bulk homogeneous sample, and that multiple scattering has a negligible effect. This technique can be applied towards understanding the complex scattering behaviour from heterogeneous tissues.

Structure Function Estimated From Histological Tissue Sections

Ultrasonic scattering is determined by not only the properties of individual scatterers but also the correlation among scatterer positions. The role of scatterer spatial correlation is significant for dense medium, but has not been fully understood. The effect of scatterer spatial correlation may be modeled by the structure function as a frequency-dependent factor in the backscatter coefficient (BSC) expression. The structure function has been previously estimated from the BSC data. The aim of this study is to estimate the structure function from histology to test if the acoustically estimated structure function is indeed caused by the scatterer spatial distribution. Hematoxylin and eosin stained histological sections from dense cell pellet biophantoms were digitized. The scatterer positions were determined manually from the histological images. The structure function was calculated from the extracted scatterer positions. The structure function obtained from histology showed reasonable agreement in the shape but not in the amplitude, compared with the structure function previously estimated from the backscattered data. Fitting a polydisperse structure function model to the histologically estimated structure function yielded relatively accurate cell radius estimates ([Formula: see text]). Furthermore, two types of mouse tumors that have similar cell size and shape but distinct cell spatial distributions were studied, where the backscattered data were shown to be related to the cell spatial distribution through the structure function estimated from histology. In conclusion, the agreement between acoustically estimated and histologically estimated structure functions suggests that the acoustically estimated structure function is related to the scatterer spatial distribution.

Quantitative Characterization of Tissue Microstructure in Concentrated Cell Pellet Biophantoms Based on the Structure Factor Model

Quantitative ultrasound (QUS) methods based on the backscatter coefficient (BSC) are typically model-based. The BSC is estimated from experiments and is fit to a model. The fit parameters are often termed QUS estimates and are used to characterize the scattering properties of the tissue under investigation. Nevertheless, for physical interpretation of QUS estimates to be accurate, the scattering model chosen must also be accurate. The goal of this work was to investigate the use of the structure factor model (SFM) to take into account coherent scattering from high volume fractions of scatterers. The study focuses on comparing the performance of two sparse models (fluid-filled sphere and Gaussian) and one concentrated model (SFM) to estimate QUS parameters from simulations and cell pellet biophantoms with a range of scatterer volume fractions. Results demonstrated the superiority of the SFM for all investigated volume fractions (i.e., from 0.006 to 0.30). In particular, the sparse models underestimated scatterer size and overestimated acoustic concentration when the volume fraction was greater than 0.12. In addition, the SFM has the ability to provide the volume fraction and the relative impedance contrast (instead of only the acoustic concentration provided by the sparse models), which could have a great benefit for tissue characterization. This study demonstrates that the SFM could prove to be an invaluable tool for QUS and could help to more accurately characterize tissue from ultrasound measurements.

High-Frequency Quantitative Ultrasound Spectroscopy of Excised Canine Livers and Mouse Tumors Using the Structure Factor Model

Three scattering models were examined for characterizing ex vivo canine livers and HT29 mouse tumors in the 10-38- and the 15-42-MHz frequency bandwidth, respectively. The spherical Gaussian model (SGM) and the fluid sphere model (FSM) that were examined are suitable for dealing with sparse media, whereas the structure factor model (SFM) is adapted for characterizing concentrated media. For the canine livers, the scatterer radius and the acoustic concentration estimated with the three models were similar and matched well the nuclear structures obtained from histological analysis (with relative errors less than 7%). These results show that the livers could be considered as a diluted medium and that the nuclei in liver could be a dominant source of scattering. For the homogeneous mouse tumors, containing mostly viable HT29 cells, scatterer radius and volume fraction estimated with the SFM showed good agreement with the whole cell structures obtained from histological analysis (with relative errors less than 15%), whereas the sparse models (the SGM and the FSM) gave no consistent quantitative ultrasound parameters. This suggests that the viable HT29 cell areas have densely packed cellular content and that the whole HT29 cell could be responsible for scattering. For the heterogeneous tumors, the hyperechogenic zones observed in the B-mode images were linked to the presence of small necrotic areas surrounded by viable HT29 cells. Comparison between sparse and concentrated models shows that these hyperechogenic zones could be considered as a concentrated medium.

Coherent and incoherent ultrasound backscatter from cell aggregates

The effective medium theory (EMT) was recently developed to model the ultrasound backscatter from aggregating red blood cells [Franceschini, Metzger, and Cloutier, IEEE Trans. Ultrason. Ferroelectr. Freq. Control 58, 2668-2679 (2011)]. The EMT assumes that aggregates can be treated as homogeneous effective scatterers, which have effective properties determined by the aggregate compactness and the acoustical characteristics of the cells and the surrounding medium. In this study, the EMT is further developed to decompose the differential backscattering cross section of a single cell aggregate into coherent and incoherent components. The coherent component corresponds to the squared norm of the average scattering amplitude from the effective scatterer, and the incoherent component considers the variance of the scattering amplitude (i.e., the mean squared norm of the fluctuation of the scattering amplitude around its mean) within the effective scatterer. A theoretical expression for the incoherent component based on the structure factor is proposed and compared with another formulation based on the Gaussian direct correlation function. This theoretical improvement is assessed using computer simulations of ultrasound backscatter from aggregating cells. The consideration of the incoherent component based on the structure factor allows us to approximate the simulations satisfactorily for a product of the wavenumber times the aggregate radius kr_ag around 2.

Unifying Concepts of Statistical and Spectral Quantitative Ultrasound Techniques

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.

Experimental application of ultrafast imaging to spectral tissue characterization

Ultrasound ultrafast imaging (UI) allows acquisition of thousands of frames per second with a sustained image quality at any depth in the field of view. Therefore, it would be ideally suited to obtain good statistical sampling of fast-moving tissues using spectral-based techniques to derive the backscatter coefficient (BSC) and associated quantitative parameters. In UI, an image is formed by insonifying the medium with plane waves steered at different angles, beamforming them and compounding the resulting radiofrequency images. We aimed at validating, experimentally, the effect of these beamforming protocols on the BSC, under both isotropic and anisotropic conditions. Using UI techniques with a linear array transducer (5-14 MHz), we estimated the BSCs of tissue-mimicking phantoms and flowing porcine blood at depths up to 35 mm with a frame rate reaching 514 Hz. UI-based data were compared with those obtained using single-element transducers and conventional focusing imaging. Results revealed that UI compounded images can produce valid estimates of BSCs and effective scatterer size (errors less than 2.2 ± 0.8 and 0.26 ± 0.2 dB for blood and phantom experiments, respectively). This work also describes the use of pre-compounded UI images (i.e., steered images) to assess the angular dependency of circulating red blood cells. We have concluded that UI data sets can be used for BSC spectral tissue analysis and anisotropy characterization.

Structure function for high-concentration biophantoms of polydisperse scatterer sizes

Ultrasonic backscattering coefficient (BSC) has been used extensively to characterize tissue. In most cases, sparse scatterer concentrations are assumed. However, many types of tissues have dense scattering media. This study addresses the problem of dense media scattering by taking into account the correlation among scatterers using the structure functions. The effect of scatterer polydispersity on the structure functions is investigated. Structure function models based on polydisperse scatterers are theoretically developed and experimentally evaluated against the structure functions obtained from cell pellet biophantoms. The biophantoms were constructed by placing live cells of known concentration in coagulation media to form a clot. The BSCs of the biophantoms were estimated using single-element transducers over the frequency range from 11 to 105 MHz. Experimental structure functions were obtained by comparing the BSCs of two cell concentrations. The structure functions predicted by the models agreed with the experimental structure functions. Fitting the models yielded cell radius estimates that were consistent with direct light microscope measures. The results demonstrate the role of scatterer position correlation on dense media scattering, and the significance of scatterer polydispersity on structure functions. This work may lead to more accurate modeling of ultrasonic scattering in dense medium for improved tissue characterization.

In vivo observation of the hypo-echoic "black hole" phenomenon in rat arterial bloodstream: a preliminary Study

The "black hole," a hypo-echoic hole at the center of the bloodstream surrounded by a hyper-echoic zone in cross-sectional views, has been observed in ultrasound backscattering measurements of blood with red blood cell aggregation in in vitro studies. We investigated whether the phenomenon occurs in the in vivo arterial bloodstream of rats using a high-frequency ultrasound imaging system. Longitudinal and cross-sectional ultrasound images of the rat common carotid artery (CCA) and abdominal aorta were obtained using a 40-MHz ultrasound system. A high-frame-rate retrospective imaging mode was employed to precisely examine the dynamic changes in blood echogenicity in the arteries. When the imaging was performed with non-invasive scanning, blood echogenicity was very low in the CCA as compared with the surrounding tissues, exhibiting no hypo-echoic zone at the center of the vessel. Invasive imaging of the CCA by incising the skin and subcutaneous tissues at the imaging area provided clearer and brighter blood echo images, showing the "black hole" phenomenon near the center of the vessel in longitudinal view. The "black hole" was also observed in the abdominal aorta under direct imaging after laparotomy. The aortic "black hole" was clearly observed in both longitudinal and cross-sectional views. Although the "black hole" was always observed near the center of the arteries during the diastolic phase, it dissipated or was off-center along with the asymmetric arterial wall dilation at systole. In conclusion, we report the first in vivo observation of the hypo-echoic "black hole" caused by the radial variation of red blood cell aggregation in arterial bloodstream.

Structure factor model for understanding the measured backscatter coefficients from concentrated cell pellet biophantoms

Ultrasonic backscatter coefficient (BSC) measurements were performed on K562 cell pellet biophantoms with cell concentrations ranging from 0.006 to 0.30 in the 10-42 MHz frequency bandwidth. Three scattering models, namely, the fluid-filled sphere model (FFSM), the particle model (PM), and the structure factor model (SFM), were compared for modeling the scattering from an ensemble of concentrated cells. A parameter estimation procedure was developed in order to estimate the scatterer size and relative impedance contrast that could explain the measured BSCs from all the studied cell concentrations. This procedure was applied to the BSC data from K562 cell pellet biophantoms in the 10-42 MHz frequency bandwidth and to the BSC data from Chinese hamster ovary cell pellet biophantoms in the 26-105 MHz frequency bandwidth given in Han, Abuhabsah, Blue, Sarwate, and O'Brien [J. Acoust. Soc. Am. 130, 4139-4147 (2011)]. The data fitting quality and the scatterer size estimates show that the SFM was more suitable than the PM and the FFSM for modeling the responses from concentrated cell pellet biophantoms.

Comparison of three scattering models for ultrasound blood characterization

Ultrasonic backscattered signals from blood contain frequency-dependent information that can be used to obtain quantitative parameters reflecting the aggregation level of red blood cells (RBCs). The approach is based on estimating structural aggregate parameters by fitting the spectrum of the backscattered radio-frequency echoes from blood to an estimated spectrum considering a theoretical scattering model. In this study, three scattering models were examined: a new implementation of the Gaussian model (GM), the structure factor size estimator (SFSE), and the new effective medium theory combined with the structure factor model (EMTSFM). The accuracy of the three scattering models in determining mean aggregate size and compactness was compared by 2-D and 3-D computer simulations in which RBC structural parameters were controlled. Two clustering conditions were studied: 1) the aggregate size varied and the aggregate compactness was fixed in both 2-D and 3-D cases, and 2) the aggregate size was fixed and the aggregate compactness varied in the 2-D case. For both clustering conditions, the EMTSFM was found to be more suitable than GM and SFSE for characterizing RBC aggregation.

Experimental assessment of four ultrasound scattering models for characterizing concentrated tissue-mimicking phantoms

Tissue-mimicking phantoms with high scatterer concentrations were examined using quantitative ultrasound techniques based on four scattering models: The Gaussian model (GM), the Faran model (FM), the structure factor model (SFM), and the particle model (PM). Experiments were conducted using 10- and 17.5-MHz focused transducers on tissue-mimicking phantoms with scatterer concentrations ranging from 1% to 25%. Theoretical backscatter coefficients (BSCs) were first compared with the experimentally measured BSCs in the forward problem framework. The measured BSC versus scatterer concentration relationship was predicted satisfactorily by the SFM and the PM. The FM and the PM overestimated the BSC magnitude at actual concentrations greater than 2.5% and 10%, respectively. The SFM was the model that better matched the BSC magnitude at all the scatterer concentrations tested. Second, the four scattering models were compared in the inverse problem framework to estimate the scatterer size and concentration from the experimentally measured BSCs. The FM did not predict the concentration accurately at actual concentrations greater than 12.5%. The SFM and PM need to be associated with another quantitative parameter to differentiate between low and high concentrations. In that case, the SFM predicted the concentration satisfactorily with relative errors below 38% at actual concentrations ranging from 10% to 25%.

Forward problem study of an effective medium model for ultrasound blood characterization

The structure factor model (SFM) is a scattering model developed to simulate the backscattering coefficient (BSC) of aggregated red blood cells (RBCs). However, the SFM can hardly be implemented to estimate the structural aggregate parameters in the framework of an inverse problem formulation. A scattering model called the effective medium theory combined with the SFM (EMTSFM) is thus proposed to approximate the SFM. The EMTSFM assumes that aggregates of RBCs can be treated as individual homogeneous scatterers, which have effective properties determined by the acoustical characteristics and concentration of RBCs within aggregates. The EMTSFM parameterizes the BSC by three indices: the aggregate radius, the concentration of RBCs with- in aggregates (the aggregate compactness), and the systemic hematocrit. The goodness of fit of the EMTSFM approximation in comparison with the SFM was then examined. Based on a 2-D study, the EMTSFM was found to approximate the SFM with relative errors less than 30% for a product of the wavenumber times the mean aggregate radius krΛκ <; 1.32. The main contribution of this work is the parameterization of the BSC with the RBC aggregate compactness, which is of relevance in clinical hemorheology because it reflects the binding energy between RBCs.

Computational models of distributed aberration in ultrasound breast imaging

Two methods for simulation of ultrasound wavefront distortion are introduced and compared with aberration produced in simulations using digitized breast tissue specimens and a conventional multiple time-shift screen model. In the first method, aberrators are generated using a computational model of breast anatomy. In the second method, 10 to 12 irregularly shaped, strongly scattering inclusions are superimposed on the multiple-screen model to create a screen-inclusion model. Linear 2-D propagation of a 7.5-MHz planar, pulsed wavefront through each aberrator is computed using a first-order k-space method. The anatomical and screen-inclusion models reproduce two characteristics of arrival-time fluctuations observed in simulations using the digitized specimens that are not represented in simulations using the multiple-screen model: non-Gaussian first-order statistics and sharp changes in the rms arrival-time fluctuation as a function of propagation distance. The anatomical and screen-inclusion models both produce energy- level fluctuations similar to the digitized specimens, but the anatomical model more closely matches the pulse-shape distortion produced by the specimens. Both aberration models can readily be extended to 3-D, and the screen-inclusion model has the advantage of simplicity of implementation. Both models should enable more rigorous evaluation of adaptive focusing algorithms than is possible using conventional time-shift screen models.

Generalized poisson 3-D scatterer distributions

This paper describes a simple, yet powerful ultrasound scatterer distribution model. The model extends a 1-D generalized Poisson process to multiple dimensions using a Hilbert curve. The model is intuitively tuned by spatial density and regularity parameters which reliably predict the first and second-order statistics of varied synthetic imagery.

High-frequency ultrasound backscattering by blood: analytical and semianalytical models of the erythrocyte cross section

This paper proposes analytical and semianalytical models of the ultrasonic backscattering cross section (BCS) of various geometrical shapes mimicking a red blood cell (RBC) for frequencies varying from 0 to 90 MHz. By assuming the first-order Born approximation and by modeling the shape of a RBC by a realistic biconcave volume, different scattering behaviors were identified for increasing frequencies. For frequencies below 18 MHz, a RBC can be considered a Rayleigh scatterer. For frequencies less than 39 MHz, the general concept of acoustic inertia tensor is introduced to describe the variation of the BCS with the frequency and the incidence direction. For frequencies below 90 MHz, ultrasound backscattering by a RBC is equivalent to backscattering by a cylinder of height 2 microm and diameter 7.8 microm. These results lay the basis of ultrasonic characterization of RBC aggregation by proposing a method that distinguishes the contribution of the individual RBC acoustical characteristics from collective effects, on the global blood backscattering coefficient. A new method of data reduction that models the frequency dependence of the ultrasonic BCS of micron-sized weak scatterers is also proposed. Applications of this method are in tissue characterization as well as in hematology.

Anatomy and flow in normal and ischemic microvasculature based on a novel temporal fractal dimension analysis algorithm using contrast enhanced ultrasound

Strategies for improvement of blood flow by promoting new vessel growth in ischemic tissue are being developed. Recently, contrast-enhanced ultrasound (CEU) imaging has been used to assess tissue perfusion in models of ischemia-related angiogenesis, growth-factor mediated angiogenesis, and tumor angiogenesis. In these studies, microvascular flow is measured in order to assess the total impact of adaptations at different vascular levels. High-resolution methods for imaging larger vessels have been developed in order to derive "angiograms" of arteries, veins, and medium to large microvessels. We describe a novel method of vascular bed (microvessel and arterial) characterization of vessel anatomy and flow simultaneously, using serial measurement of the fractal dimension (FD) of a temporal sequence of CEU images. This method is proposed as an experimental methodology to distinguish ischemic from nonischemic tissue. Moreover, an improved approach for extracting the FD unique to this application is introduced.

Effect of red cell clustering and anisotropy on ultrasound blood backscatter: a Monte Carlo study

When flowing at a low shear rate, blood appears hyperechogenic on ultrasound B-scans. The formation of red blood cell (RBC) aggregates that also alters blood viscosity is the microscopic mechanism explaining this acoustical phenomenon. In this study, Monte Carlo simulations were performed to predict how RBC clustering increases ultrasound scattering by blood. A bidimensional Gibbs-Markov random point process parameterized by the adhesion energy epsilon and an anisotropy index nu was used to describe RBC positions for a hematocrit H = 40%. The frequency dependence of the backscattering coefficient chi(f) was computed using Born approximation. The backscattering coefficient chi0 at 5 MHz and the spectral slopes n(x) and n(y) (chi alpha f(nx) or f(ny)) measured, respectively, when the insonification is parallel and perpendicular with the RBC cluster axis were then extracted. Under isotropic conditions, chi0 increased up to 7 dB with epsilon and n(x) = n(y) decreased from 4.2 to 3.4. Under anisotropic conditions, the backscattering was stronger perpendicularly to aggregate axis, resulting in n(x) < n(y). The anisotropy in scattering appeared more pronounced when epsilon or nu increased. These two dimensional results generally predict that low-frequency blood backscatter is related to cluster dimension, and higher-frequency properties are affected by finer morphological features as anisotropy. This numerically establishes that ultrasound backscatter spectroscopy on a large frequency range is pertinent to characterize in situ hemorheology.

Non-Gaussian statistics and temporal variations of the ultrasound signal backscattered by blood at frequencies between 10 and 58 MHz

Very little is known about the blood backscattering behavior and signal statistics following flow stoppage at frequencies higher than 10 MHz. Measurements of the radio frequency (rf) signals backscattered by normal human blood (hematocrit = 40%, temperature = 37 degrees C) were performed in a tube flow model at mean frequencies varying between 10 and 58 MHz. The range of increase of the backscattered power during red blood cell (RBC) rouleau formation was close to 15 dB at 10 and 36 MHz, and dropped, for the same blood samples, below 8 dB at 58 MHz. Increasing the frequency from 10 to 58 MHz raised the slope of the power changes at the beginning of the kinetics of aggregation, and could emphasize the non-Gaussian behavior of the rf signals interpreted in terms of the K and Nakagami statistical models. At 36 and 58 MHz, significant increases of the kurtosis coefficient, and significant reductions of the Nakagami parameter were noted during the first 30 s of flow stoppage. In conclusion, increasing the transducer frequency reduced the magnitude of the backscattered power changes attributed to the phenomenon of RBC aggregation, but improved the detection of rapid growth in aggregate sizes and non-Gaussian statistical behavior.

Modeling the frequency dependence (5-120 MHz) of ultrasound backscattering by red cell aggregates in shear flow at a normal hematocrit

The frequency dependence of the ultrasound signal backscattered by blood in shear flow was studied using a simulation model. The ultrasound backscattered signal was computed with a linear model that considers the characteristics of the ultrasound system and tissue acoustic properties. The tissue scattering properties were related to the position and shape of the red blood cells (RBCs). A 2D microrheological model simulated the RBC dynamics in a Couette shear flow system. This iterative model, described earlier [Biophys. J. 82, 1696-1710 (2002)], integrates the hydrodynamic effect of the flow, as well as adhesive and repulsive forces between RBCs. RBC aggregation was simulated at 40% hematocrit and shear rates of 0.05-2 s(-1). The RBC aggregate sizes ranged, on average, from 3.3 RBCs at 2 s(-1) to 33.5 cells at 0.05 s(-1). The ultrasound backscattered power was studied at frequencies between 5-120 MHz and insonification angles between 0-180 degrees. At frequencies below approximately 30 MHz, the ultrasound backscattered power increased as the shear rate was decreased and the size of the aggregates was raised. A totally different scattering behavior was noted above 30 MHz. Typical spectral slopes of the backscattered power (log-log scale) between 5-25 MHz equaled 3.8, whereas slopes down to 0.6 were measured at 0.05 s(-1), between 40-60 MHz. The ultrasound backscattered power was shown to be angle dependent at low frequencies (5-25 MHz). The anisotropy persisted at high frequencies (>25 MHz) for small aggregates (at 2 s(-1)). In conclusion, this study sheds some light on the blood backscattering behavior with an emphasis on the non-Rayleigh regime. Additional experimental studies may be necessary to validate the simulation results, and to fully understand the relation between the ultrasound backscattered power, level of RBC aggregation, shear rate, frequency, and insonification angle.

A new high resolution color flow system using an eigendecomposition-based adaptive filter for clutter rejection

We present a new signal processing strategy for high frequency color flow mapping in moving tissue environments. A new application of an eigendecomposition-based clutter rejection filter is presented with modifications to deal with high blood-to-clutter ratios (BCR). Additionally, a new method for correcting blood velocity estimates with an estimated tissue motion profile is detailed. The performance of the clutter filter and velocity estimation strategies is quantified using a new swept-scan signal model. In vivo color flow images are presented to illustrate the potential of the system for mapping blood flow in the microcirculation with external tissue motion.

A new high resolution color flow system using an eigendecomposition-based adaptive filter for clutter rejection

We present a new signal processing strategy for high frequency color flow mapping in moving tissue environments. A new application of an eigendecomposition-based clutter rejection filter is presented with modifications to deal with high blood-to-clutter ratios (BCR). Additionally, a new method for correcting blood velocity estimates with an estimated tissue motion profile is detailed. The performance of the clutter filter and velocity estimation strategies is quantified using a new swept-scan signal model. In vivo color flow images are presented to illustrate the potential of the system for mapping blood flow in the microcirculation with external tissue motion.