Evaluation of the spectral fit algorithm as functions of frequency range and deltakaeff

Attenuation Compensation and Estimation

Estimating the loss of ultrasound signal with propagation depth as a function of frequency is essential for quantifying tissue properties. Specifically, ultrasound attenuation is used to correct for spectral distortion prior to estimating quantitative ultrasound parameters to assess the tissue. Ultrasound attenuation can also be used independently to characterize the tissue. In this chapter, we review the primary algorithms for estimating both the local attenuation within a region of interest as well as the total attenuation between a region of interest and an ultrasound source. The strengths and weaknesses of each algorithm are also discussed.

Added Value of Quantitative Ultrasound and Machine Learning in BI-RADS 4-5 Assessment of Solid Breast Lesions

The purpose of this study was to evaluate various combinations of 13 features based on shear wave elasticity (SWE), statistical and spectral backscatter properties of tissues, along with the Breast Imaging Reporting and Data System (BI-RADS), for classification of solid breast lesions at ultrasonography by means of random forests. One hundred and three women with 103 suspicious solid breast lesions (BI-RADS categories 4-5) were enrolled. Before biopsy, additional SWE images and a cine sequence of ultrasound images were obtained. The contours of lesions were delineated, and parametric maps of the homodyned-K distribution were computed on three regions: intra-tumoral, supra-tumoral and infra-tumoral zones. Maximum elasticity and total attenuation coefficient were also extracted. Random forests yielded receiver operating characteristic (ROC) curves for various combinations of features. Adding BI-RADS category improved the classification performance of other features. The best result was an area under the ROC curve of 0.97, with 75.9% specificity at 98% sensitivity.

Protocol for Robust In Vivo Measurements of Erythrocyte Aggregation Using Ultrasound Spectroscopy

Erythrocyte aggregation is a non-specific marker of acute and chronic inflammation. Although it is usual to evaluate this phenomenon from blood samples analyzed in laboratory instruments, in vivo real-time assessment of aggregation is possible with spectral ultrasound techniques. However, variable blood flow can affect the interpretation of acoustic measures. Therefore, flow standardization is required. Two techniques of flow standardization were evaluated with porcine and equine blood samples in Couette flow. These techniques consisted in either stopping the flow or reducing it. Then, the sensibility and repeatability of the retained method were evaluated in 11 human volunteers. We observed that stopping the flow compromised interpretation and repeatability. Conversely, maintaining a low flow provided repeatable measures and could distinguish between normal and high extents of erythrocyte aggregation. Agreement was observed between in vivo and ex vivo measures of the phenomenon (R² = 82.7%, p value < 0.0001). These results support the feasibility of assessing in vivo erythrocyte aggregation in humans by quantitative ultrasound means.

Scatterer size estimation using the center frequency assessed from ultrasound time domain data

Scatterer size estimation is useful when characterizing tissue using ultrasound. In all previous studies on scatterer size, the estimations are performed in the frequency domain and are thus subjected to a trade off in time-frequency resolution. This study focused on the feasibility of estimating scatterer size in the time domain using only the ultrasound center frequency, assuming a Gaussian-shaped pulse. A model for frequency normalization was derived and the frequency-dependent attenuation was compensated. Five phantoms with well-defined sizes of spherical glass beads were made and scanned with two different linear array transducers with variable center frequencies. A strong correlation (r = 0.99, p < 10^-19) between the backscattered center frequency and the product between the wave number and scatterer radius was demonstrated. On average the scatterer diameter was underestimated by 6% ± 24%. These results suggest that estimation of scatterer size is possible using only the center frequency assessed in the time domain.

Ultrasound imaging and characterization of biofilms based on wavelet de-noised radiofrequency data

The ability to non-invasively image and characterize bacterial biofilms in children during nasopharyngeal colonization with potential otopathogens and during acute otitis media would represent a significant advance. We sought to determine if quantitative high-frequency ultrasound techniques could be used to achieve that goal. Systematic time studies of bacterial biofilm formation were performed on three preparations of an isolated Haemophilus influenzae (NTHi) strain, a Streptococcus pneumoniae (Sp) strain and a combination of H. influenzae and S. pneumoniae (NTHi + Sp) in an in vitro environment. The process of characterization included conditioning of the acquired radiofrequency data obtained with a 15-MHz focused, piston transducer by using a seven-level wavelet decomposition scheme to de-noise the individual A-lines acquired. All subsequent spectral parameter estimations were done on the wavelet de-noised radiofrequency data. Various spectral parameters-peak frequency shift, bandwidth reduction and integrated backscatter coefficient-were recorded. These parameters were successfully used to map the progression of the biofilms in time and to differentiate between single- and multiple-species biofilms. Results were compared with those for confocal microscopy and theoretical evaluation of form factor. We conclude that high-frequency ultrasound may prove a useful modality to detect and characterize bacterial biofilms in humans as they form on tissues and plastic materials.

Improved algorithm for estimation of attenuation along propagation path using backscattered echoes from multiple sources

Accurately determining the attenuation along the propagation path leading to a region of interest could significantly improve diagnostic ultrasound tissue characterization since tissue characterization requires exact compensation for the frequency-dependent attenuation along the propagation path. In a previous study (JASA, 124:1367, 2008), it was shown that the total attenuation can be determined by using the backscattered echoes from multiple sources. The preliminary computer simulation results, had an average error between -0.3 and +0.2dB/MHz for the cases tested with a trend towards increasing error with increasing correlation length (i.e., characteristic size of the tissue microstructure of the scattering medium) and attenuation along the propagation path. Therefore, the goal of this study was to improve the accuracy of previously derived algorithm and reduce the dependence of the algorithm on correlation length and attenuation. In this study, the previous derivations were redone and the assumptions made by the algorithm regarding the scattering properties of the medium and the shape of the backscattered power spectrum were relaxed. The revised algorithm was then verified using computer simulations of five sources (6, 8, 10, 12, and 14MHz, 50% bandwidth) exposing a homogeneous tissue region. The simulated tissue had microstructure following a Gaussian spatial correlation function (i.e., exp(-0.827(ka(eff))(2)) where k is the wavenumber) with effective radii, a(eff), of 5-55microm (one size per simulated case) placed at a density of 250/mm(3) ( approximately 5 scatterers/resolution cell for 14MHz transducer). The attenuation of the tissue was also varied from 0.1 to 0.9dB/cm-MHz. The computer simulations demonstrated that the modifications significantly improved the accuracy of the algorithm resulting in average errors between -0.04 and 0.1dB/MHz which is three times better than the error performance of the original algorithm.

Estimating the total ultrasound attenuation along the propagation path by applying multiple filters to backscattered echoes from a single spherically focused source

Quantifying the correlation length of the tissue microstructure has shown potential for differentiating between benign and malignant tumors. To implement these advances in the clinic, the total frequency-dependent attenuation along the propagation path must be determined on a patient specific basis. Previously, an algorithm was developed to estimate this attenuation using echoes from multiple sources. In this study, the developed algorithm was extended to echoes from a single source by filtering the echoed signal into multiple frequency bands. This step was needed because it would be challenging to scan exactly the same tissue region using multiple sources in the clinic. Computer simulations and phantom experiments were conducted to verify the attenuation could be determined by filtering the echoes from a single source. The simulations utilized a spherically focused single-element source (5 cm focal length, f/4, 14 MHz center frequency, 50% bandwidth) exposing a homogeneous tissue region (Gaussian scattering structures with effective radii of 5 to 55 mum at a density of 250/mm(3), attenuation of 0.1 to 0.9 dB/cm.MHz). The phantom experiments utilized a spherically focused single-element source (5.08 cm focal length, f/4, 7.5 MHz center frequency) exposing a 0.5 dB/cm.MHz homogeneous glass bead phantom. The computer simulations and phantom experiment confirmed that the total attenuation along the propagation path can be determined by appropriately applying multiple filters to the backscattered echoes from a single source.

Ultrasound characterization of red blood cell aggregation with intervening attenuating tissue-mimicking phantoms

The analysis of the ultrasonic frequency-dependent backscatter coefficient of aggregating red blood cells reveals information about blood structural properties. The difficulty in applying this technique in vivo is due to the frequency-dependent attenuation caused by intervening tissue layers that distorts the spectral content of signals backscattered by blood. An optimization method is proposed to simultaneously estimate tissue attenuation and blood structure properties, and was termed the structure factor size and attenuation estimator (SFSAE). An ultrasound scanner equipped with a wide-band 25 MHz probe was used to insonify porcine blood sheared in both Couette and tubular flow devices. Since skin is one of the most attenuating tissue layers during in vivo scanning, four skin-mimicking phantoms with different attenuation coefficients were introduced between the transducer and the blood flow. The SFSAE gave estimates with relative errors below 25% for attenuations between 0.115 and 0.411 dBMHz and kR<2.08 (k being the wave number and R the aggregate radius). The SFSAE can be useful to examine in vivo and in situ abnormal blood conditions suspected to promote pathophysiological cardiovascular consequences.

Ultrasound attenuation estimation using backscattered echoes from multiple sources

The objective of this study was to devise an algorithm that can accurately estimate the attenuation along the propagation path (i.e., the total attenuation) from backscattered echoes. It was shown that the downshift in the center frequency of the backscattered ultrasound echoes compared to echoes obtained in a water bath was calculated to have the form Deltaf=mf(o)+b after normalizing with respect to the source bandwidth where m depends on the correlation length, b depends on the total attenuation, and f(o) is the center frequency of the source as measured from a reference echo. Therefore, the total attenuation can be determined independent of the scatterer correlation length by measuring the downshift in center frequency from multiple sources (i.e., different f(o)) and fitting a line to the measured shifts versus f(o). The intercept of the line gives the total attenuation along the propagation path. The calculations were verified using computer simulations of five spherically focused sources with 50% bandwidths and center frequencies of 6, 8, 10, 12, and 14 MHz. The simulated tissue had Gaussian scattering structures with effective radii of 25 mum placed at a density of 250 mm(3). The attenuation of the tissue was varied from 0.1 to 0.9 dB / cm-MHz. The error in the attenuation along the propagation path ranged from -3.5+/-14.7% for a tissue attenuation of 0.1 dB / cm-MHz to -7.0+/-3.1% for a tissue attenuation of 0.9 dB / cm-MHz demonstrating that the attenuation along the propagation path could be accurately determined using backscattered echoes from multiple sources using the derived algorithm.

Impact of local attenuation approximations when estimating correlation length from backscattered ultrasound echoes

Estimating the characteristic correlation length of tissue microstructure from the backscattered power spectrum could improve the diagnostic capability of medical ultrasound. Previously, size estimates were obtained after compensating for source focusing, the frequency-dependent attenuation along the propagation path (total attenuation), and the frequency-dependent attenuation in the scattering region (local attenuation). In this study, the impact of approximations of the local attenuation on the scatterer size estimate was determined using computer simulations and theoretical analysis. The simulations used Gaussian impedance distributions with an effective radius of 25 microm randomly positioned in a homogeneous half-space sonified by a spherically focused source (f/1 to f/4). The approximations of the local attenuation that were assessed neglected local attenuation (i.e., assume 0 dB/cm-MHz) neglected frequency dependence of the local attenuation, and assumed a finite frequency dependence (i.e., 0.5 dB/cm-MHz) independent of the true attenuation of the medium. Errors in the scatterer size estimate due to the local attenuation approximations increased with increasing window length, increasing true local attenuation and increasing f number. The most robust estimates were obtained when the local attenuation was approximated by a tissue-independent attenuation value that was greater than 70% of the largest attenuation expected in the tissue region of interest.