Estimating myocardial attenuation from M-mode ultrasonic backscatter

Evaluation method of the degree of red blood cell aggregation considering ultrasonic propagation attenuation by analyzing ultrasonic backscattering properties

PURPOSE

Red blood cell (RBC) aggregation is one of the main factors that determines blood viscosity and an important indicator for evaluating blood properties. As a noninvasive and quantitative method for diagnosing blood properties, our research group estimated the size of RBC aggregates by fitting the scattered power spectrum from the blood vessel lumen with the theoretical scattering characteristics to evaluate the degree of RBC aggregation. However, it was assumed that the propagation attenuation of ultrasound in the vascular lumen was the same regardless of whether RBCs were aggregated or not, which caused systematic errors in the estimated size.

METHODS

To improve the size estimation accuracy, we calculated and corrected the attenuation of the blood vessel lumen during RBC aggregation and non-aggregation. The attenuation in the blood vessel lumen was calculated with the spectra acquired from two different depths.

RESULTS

In the basic experiments using microparticles, the estimation accuracy decreased as the concentration increased in the case of the conventional method, but the estimated size tended to approach the true size irrespective of the concentration, removing the propagation attenuation component with the proposed method. In the in vivo experiment on the human hand dorsal vein, the size was estimated to be larger during RBC aggregation and smaller during non-aggregation using the proposed method.

CONCLUSION

These results suggest that the proposed method can provide precise size estimation by considering the propagation attenuation component regardless of differences in blood conditions such as RBC concentration and degree of aggregation.

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

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

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

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

Performance evaluation of the spectral centroid downshift method for attenuation estimation

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

Regional variation in the measured apparent ultrasonic backscatter of mid-gestational fetal pig hearts

The goal of this study was to characterize and compare regional backscatter properties of fetal hearts through measurements of the apparent integrated backscatter. Sixteen excised, formalin-fixed fetal pig hearts, representing an estimated 53 to 63 days of gestation, were investigated. Spatially localized measurements of integrated backscatter from these specimens were acquired using a 50 MHz single-element transducer. The apparent integrated backscatter measurements demonstrate different patterns of backscatter from the myocardium of the right ventricle compared with that of the left ventricle. These backscatter measurements appear to be consistent with the anisotropy of the fiber orientation observed in histologic assessment of the same specimens. For each of the 16 hearts, the apparent integrated backscatter from the right ventricular myocardium was larger than that from the left ventricular myocardium, exhibiting mean apparent backscatter values of -35.9 +/- 2.0 dB and -40.1 +/- 1.9 dB (mean +/- standard deviation; n = 16; p < 0.001), respectively. This study suggests that the intrinsic ultrasonic properties of the left and right ventricular myocardium are distinct in fetal pig hearts at mid-gestation.

Ultrasonic detection of the anisotropy of protein cross linking in myocardium at diagnostic frequencies

Increased myocardial stiffness in aging and diabetes that may result in pathologies such as diastolic dysfunction has been attributed, in part, to an increase in cross linking of extracellular matrix proteins such as collagen. With the development of new approaches to cardiovascular therapy, it becomes increasingly important to develop noninvasive approaches for monitoring changes in myocardial cross linking. The objective of this study was to use ultrasound at frequencies used in clinical echocardiography to measure changes in myocardial attenuation resulting from increased cross linking as a function of angle of insonification over a complete rotation. Through-transmission radiofrequency-based measurements were performed on 36 specimens from 12 freshly excised ovine hearts at room temperature, which were then fixed in formalin to induce protein cross linking prior to repeated measurements. For angles near perpendicular to the myofiber direction, the measured slope of attenuation increased from 0.52 +/- 0.07 dB/(cm MHz) (mean +/- one standard deviation) for freshly excised to 0.85 +/- 0.08 dB/(cm MHz) for formalin-fixed myocardium. In contrast, results for parallel insonification exhibit considerable overlap (1.88 +/- 0.17 for freshly excised and 1.75 +/- 0.19 dB/(cm MHz) for formalin-fixed myocardium). Results of this study suggest that the response of the extracellular collagenous matrix to changes in cross linking is directionally dependent. The anisotropy of ultrasonic attenuation thus may provide an approach for noninvasive monitoring of the extent and progression of myocardial disease associated with changes in protein cross linking. Accounting for effects due to anisotropy may be essential for the future detection of such changes using ultrasonic attenuation in vivo.

Anisotropy of the backscatter coefficient of formalin-fixed ovine myocardium

The objective of this study was to measure the backscatter coefficient of formalin-fixed myocardial tissue as a function of angle of insonification relative to the myocardial fiber direction. Backscatter measurements were performed on eight cylindrical formalin-fixed lamb myocardial specimens and compensated for attenuation and diffraction effects to determine the backscatter coefficient. The backscatter coefficient at 5 MHz was found to be maximum for insonification perpendicular to the predominant myofiber orientation and minimum for parallel insonification, with values of (17+/-14) and (1.2+/-0.7) x 10(-4) cm(-1) sr(-1) (mean+/-standard deviation), respectively.

Characterization of anisotropic myocardial backscatter using spectral slope, intercept and midband fit parameters

The specific myocardial structural components that contribute to the observed level of backscatter from the heart and its dependence on the angle of insonification have not been completely identified: The objectives of this study were to measure the anisotropy of backscatter from myocardium using the approach first introduced by Lizzi et al. [J Acoust Soc Am 73, 1366-1373 (1983)] and to use the extracted spectral parameters (spectral slope, intercept and midband fit) to characterize changes in the apparent scatterer size, spatial concentration and acoustic impedance properties as functions of the angle of insonification. Backscatter measurements were performed in vitro on eight cylindrical formalin-fixed lamb myocardial specimens using a 5 MHz focused transducer. The backscattered spectral data as a function of angle of insonification relative to the myocardial fiber direction were analyzed over the frequency range of 4 to 6 MHz. The spectral parameters describing features of backscatter were determined by applying a linear fit to attenuation-compensated normalized spectra. Results show that values of the spectral slope do not exhibit a significant dependence on the angle of insonification within uncertainties; however, the zero-frequency intercept showed clear anisotropy and was found to be a maximum for insonification perpendicular to the predominant myofiber orientation and a minimum for parallel insonification. A comparison of midband fit values at 5 MHz with attenuation-compensated integrated backscatter values showed excellent agreement for all angles of insonification. These data suggest that measurements of spectral slope, intercept, and midband fit can provide insights regarding aspects of tissue microstructure underlying the observed anisotropy of myocardial scattering properties. Measurements of the slope parameter suggest a very modest change in effective scatterer size with angle of insonification. However, the observed anisotropy in intercept and midband fit and apparent absence of anisotropy in the spectral slope suggests an angle of insonification dependence of acoustic concentration, the combination of effective spatial scatterer concentration and acoustic impedance properties, without a significant contribution from changes in effective scatterer size.

Ultrasonic radio-frequency spectrum analysis of normal brain tissue

Acoustic tissue properties can be estimated using texture and/or spectral parameter analysis. Spectral analysis is based on the rf-signals whose frequency-content is commonly neglected in conventional B-mode imaging. Attenuation and backscatter values of normal brain tissue were analyzed. Unprocessed rf-data of 20 patients were sampled intraoperatively after craniotomy using a modified conventional ultrasonic device (Hitachi CS 9600) and analyzed off-line by a custom-made software routine. Before parameter estimation, influences of the diffraction pattern were compensated by means of a correction function obtained using a tissue-mimicking phantom. Attenuation of white matter showed a linear frequency dependence with a slope of 0.94 +/- 0.13 dB cm(-1) MHz(-1). The spectral slope was determined using 10 distinct frequencies between 2.5 and 5.75 MHz. Backscattering properties were analyzed by determining the power spectral density (PSD) and a relative backscatter coefficient (rel BSC) against the values derived from the tissue-mimicking phantom. PSD and rel BSC values were frequency-dependent, with highest PSD values at the probe's center frequency (-75.69 +/- 8.26 dB V(2) Hz(-1)). The corresponding rel BSC value at 5 MHz was determined as 15.39 +/- 8.26 dB. Finally, backscatter coefficients (BSC) of brain tissue were computed using the known BSC of the phantom. The data provided in this study are meant to serve as a base for intended future characterization of brain tissue that potentially allows intraoperative differentiation between normal and pathologic areas and therefore provides the surgeon with additional information for defining the extent of resection in brain more precisely.

Attenuation estimation using spectral cross-correlation

Estimation of the local attenuation coefficient in soft tissue is important both for clinical diagnosis and for further analysis of ultrasound B-mode images. However, it is difficult to extract spectral properties in a small region of interest from noisy backscattered ultrasound radio frequency (RF) signals. Diffraction effects due to transducer beam focal properties also have to be corrected for accurate estimation of the attenuation coefficient. In this paper, we propose a new attenuation estimation method using spectral cross-correlation between consecutive power spectra obtained from the backscattered RF signals at different depths. Since the spectral cross-correlation method estimates the spectral shift by comparing the entire power spectra, it is more robust and stable to the spectral noise artifacts in the backscattered RF signals. A diffraction compensation technique using a reference phantom with a known attenuation coefficient value is also presented. Local attenuation coefficient estimates obtained using spectral cross-correlation are within 2.3% of the actual value with small estimation variances, as demonstrated in the simulation results.

Ultrasound imaging

Ultrasound imaging is now in very widespread clinical use. The most important underpinning technologies include transducers, beam forming, pulse compression, tissue harmonic imaging, contrast agents, techniques for measuring blood flow and tissue motion, and three-dimensional imaging. Specialized and emerging technologies include tissue characterization and image segmentation, microscanning and intravascular scanning, elasticity imaging, reflex transmission imaging, computed tomography, Doppler tomography, photoacoustics and thermoacoustics. Phantoms and quality assurance are necessary to maintain imaging performance. Contemporary ultrasonic imaging procedures seem to be safe but studies of bioeffects are continuing. It is concluded that advances in ultrasonic imaging have primarily been pushed by the application of physics and innovations in engineering, rather than being pulled by the identification of specific clinical objectives in need of scientific solutions. Moreover, the opportunities for innovation to continue into the future are both challenging and exciting.

Measurements of the anisotropy of ultrasonic attenuation in freshly excised myocardium

Echocardiography requires imaging of the heart with sound propagating at varying angles relative to the predominant direction of the myofibers. The degree of anisotropy of attenuation can significantly influence ultrasonic imaging and tissue characterization measurements in vivo. This study quantifies the anisotropy of attenuation of freshly excised myocardium at frequencies typical of echocardiographic imaging. Results show a significantly larger anisotropy than previously reported in specimens of locally unidirectional myofibers. Through-transmission radio frequency-based measurements were performed on specimens from 12 ovine and 12 bovine hearts. Although ovine hearts are closer in size to human, the larger size of bovine hearts offers the potential for specimens in which myofibers are more nearly unidirectionally aligned. The attenuation coefficient increased approximately linearly with frequency. The mean slope of attenuation with frequency was 3-4 times larger for propagation parallel than for perpendicular to the myofibers. At perpendicular insonification, slopes between ovine and bovine myocardium were approximately equal. However, attenuation in bovine specimens was larger for angles approaching parallel. The difference in results for parallel appears consistent with what might be expected from increased myofiber curvature associated with smaller lamb hearts. Quantitative knowledge of anisotropy of attenuation may be useful in understanding mechanisms underlying the interaction of ultrasound with myocardium.

Elastic stiffness coefficients (c11, C33, and C13) for freshly excised and formalin-fixed myocardium from ultrasonic velocity measurements

The goal of this study was to measure elastic stiffness coefficients of freshly excised and subsequently formalin-fixed myocardial tissue. Our approach was to measure the angle-dependent phase velocities associated with the propagation of a longitudinal ultrasonic wave (3-8 MHz) in ovine myocardium using phase spectroscopy techniques and to interpret the results in the context of orthotropic and transversely isotropic models describing the elastic properties of myocardium. The phase velocity results together with density measurements were used to obtain the elastic stiffness coefficients c11, c33, and c13 for both symmetries. The results for the elastic stiffness coefficients c11, c33, and c13 are the same for both symmetries. Measurements for freshly excised myocardium and the same tissue after a period of formalin fixation were compared to examine the impact of fixation on the elastic stiffness coefficients.