Anisotropy of the ultrasonic backscatter of myocardial tissue: II. Measurements in vivo

The purpose of this investigation was to determine the angular dependence of the backscatter from canine myocardial tissue in vivo and to compare it with the variation of backscatter over the cardiac cycle that has been recognized and reported previously. The backscatter was measured from regions of left ventricular wall in canine hearts in which the fibers of the muscle lay parallel to the surface of the heart and were oriented predominantly in a circumferential fashion. Because of technical considerations, the angle of insonification was varied systematically through two cycles in which the angle relative to the muscle fiber axes ranged from 60 degrees-120 degrees. Backscatter was maximum at angles of interrogation perpendicular to the myocardial fibers and minimum at those most acute (60 degrees) relative to the orientation of the fibers. The previously observed variation of integrated backscatter over the heart cycle was evident at each angle of interrogation. At end systole, the average maximum-to-minimum angular variation of integrated backscatter as 5.0 +/- 0.4 dB. At end diastole, the average maximum-to-minimum angular variation was 3.2 +/- 0.4 dB. Thus, even though angular dependence of the backscatter from tissues with directionally oriented structures is substantial, the anisotropy does not account for cardiac-cycle-dependent variation of backscatter. Accordingly, the angular dependence should be incorporated in approaches to quantitative tissue characterization with ultrasound.     

Measurements of the anisotropy of ultrasonic velocity in freshly excised and formalin-fixed myocardial tissue

The objective of this study was to quantify the anisotropy of ultrasonic velocity in freshly excised myocardial tissue and to examine the effects of formalin-fixation. Through-transmission radio-frequency-based measurements were performed on ovine and bovine myocardial specimens from 24 different hearts. A total of 81 specimens were obtained from specific locations within each heart to investigate the possibility of regional differences in anisotropy of velocity in the left ventricular wall and septum. No regional differences were observed for either lamb or cow myocardial specimens. In addition, no specific species-dependent differences were observed between ovine and bovine myocardium. Average values of velocity at room temperature for perpendicular and parallel insonification were 1556.9 +/- 0.6 and 1565.2 +/- 0.7 m/s (mean +/- standard error), respectively, for bovine myocardium (N=45) and 1556.3 +/- 0.6 and 1564.7 +/- 0.7 m/s for ovine myocardium (N=36). Immediately after measurements of freshly excised myocardium, ovine specimens were fixed in formalin for at least one month and then measurements were repeated. Formalin-fixation appears to increase the overall velocity at all angles of insonification and to increase the magnitude of anisotropy of velocity.     

Contraction-related variation in frequency dependence of acoustic properties of canine myocardium

Results of experiments performed in several laboratories indicate that contracting myocardium exhibits a cyclic variation of the magnitude of ultrasonic backscatter, with maxima occurring at end-diastole and minima at end-systole. The mechanisms responsible for this variation are not well understood. The purpose of the present study was to determine whether the frequency dependence of backscatter exhibits systematic variation throughout the cardiac cycle, analysis of which may facilitate improved understanding of biologic factors responsible for the cyclic variation of the magnitude of backscatter. In this study, the myocardial backscatter coefficient, as a function of frequency, was measured throughout the cardiac cycle in nine open-chest dogs. The frequency dependence of the backscatter coefficient was computed from a least-squares linear fit to log backscatter coefficient versus log frequency data. A cyclic variation of frequency dependence of backscatter was found with maximum near end-diastole (f2.6 +/- 0.1) and minimum near end-systole (f2.2 +/- 0.1), a significant variation (p less than 0.01). These results suggest that mechanisms responsible for the cyclic variation of backscatter may include changes in the effective size of the dominant scatterers throughout the cardiac cycle. An alternative explanation for the observed variation is an increase in the myocardial attenuation coefficient during systole followed by a decrease in diastole.     

Computer simulation model on ultrasonic myocardial backscatter

Myocardial ultrasonic tissue characterization with integrated backscatter (IB) has been studied in recent years. To evaluate the dynamic characters of the myocardium, dynamic tracing of each unit volume of myocardium is a key point. Thus, a 2D CVIB imaging algorithm based on an auto-tracing method was developed in our previous work. Where, an auto-tracking method was used in CVIB-weighted images so that the POI in each frame throughout the cardiac cycle is automatically traced to reflect the real myocardial tissue movement. It is necessary to validate the auto-tracing method. However, at present there are no appropriate models to support the process. In order to validate the auto-tracing method, two myocardium models have been established to simulate the short and long axis view of echocardiography. The motions of partially ischemic myocardium have been simulated. The models have been used to validate the 2D CVIB imaging methods in the detection of myocardial ischemia. The simulation results show that the auto-tracing algorithm is effective. The model developed in this work provides a tool for studying and developing technologies in myocardial behavior.     

Identification of human myocardial infarction in vitro based on the frequency dependence of ultrasonic backscatter.

It has been reported previously that acute and mature myocardial infarction in dogs can be differentiated in vitro and in vivo by ultrasonic tissue characterization based on measurement of the frequency dependence of ultrasonic backscatter. To characterize human infarction with an index of the frequency dependence of backscatter that could be obtained in patients, cylindrical biopsy specimens from 7 normal regions and 12 regions of infarction of 6 fixed, explanted human hearts in 2-deg steps around their entire circumference with a 5-MHz broadband transducer were insonified. One to six consecutive transmural levels were studied for each specimen. The dependence of apparent (uncompensated for attenuation or beam width) backscatter, /B(f)/2, on frequency (f) was computed from spectral analyses of radio-frequency data as /B(f)/2 = afn, where from theoretical considerations the magnitude of n decreases as scatterer size increases. Apparent integrated backscatter was computed as the average of /B(f)/2 from 3 to 7 MHz. The average value for n for normal tissue (0.9 +/- 0.1) exceeded that for tissue from regions of infarction (0.6 +/- 0.1; p less than 0.05). Infarct manifested a significant decrease of n from epicardial to endocardial levels (epi----mid----endo: 0.9----0.7----0.2; p less than 0.05) whereas normal tissue manifested similar values for n at each transmural level (0.8----1.1----0.9; p = NS). Average integrated backscatter across all transmural levels for infarct was significantly greater than for normal tissue (-48.3 +/- 0.5 vs -53.4 +/- 0.4 dB, infarct versus normal; p less than 0.05). The presence of fibrosis was associated with smaller values of n and greater integrated backscatter.(ABSTRACT TRUNCATED AT 250 WORDS)     

The role of cardiac tissue structure in defibrillation

The purpose of this paper is to investigate the relationship between cardiac tissue structure, applied electric field, and the transmembrane potential induced in the process of defibrillation. It outlines a general understanding of the structural mechanisms that contribute to the outcome of a defibrillation shock. Electric shocks defibrillate by changing the transmembrane potential throughout the myocardium. In this process first and foremost the shock current must access the bulk of myocardial mass. The exogenous current traverses the myocardium along convoluted intracellular and extracellular pathways channeled by the tissue structure. Since individual fibers follow curved pathways in the heart, and the fiber direction rotates across the ventricular wall, the applied current perpetually engages in redistribution between the intra- and extracellular domains. This redistribution results in changes in transmembrane potential (membrane polarization): regions of membrane hyper- and depolarization of extent larger than a single cell are induced in the myocardium by the defibrillation shock. Tissue inhomogeneities also contribute to local membrane polarization in the myocardium which is superimposed over the large-scale polarization associated with the fibrous organization of the myocardium. The paper presents simulation results that illustrate various mechanisms by which cardiac tissue structure assists the changes in transmembrane potential throughout the myocardium. (c) 1998 American Institute of Physics.     

Effect of the architecture of the left ventricle on the speed of the excitation wave in muscle fibers

It is known that preferential paths for the propagation of an electrical excitation wave in the human ventricular myocardium are associated with muscle fibers in tissue. The speed of the excitation wave along a fiber is several times higher than that across the direction of the fiber. To estimate the effect of the architecture and anisotropy of the myocardium of the left ventricle on the process of its electrical activation, we have studied the relation between the speed of the electrical excitation wave in a one-dimensional isolated myocardial fiber consisting of sequentially coupled cardiomyocytes and in an identical fiber located in the wall of a threedimensional anatomical model of the left ventricle. It has been shown that the speed of a wavefront along the fiber in the three-dimensional myocardial tissue is much higher than that in the one-dimensional fiber. The acceleration of the signal is due to the rotation of directions of fibers in the wall and to the position of the excitation wavefront with respect to the direction of this fiber. The observed phenomenon is caused by the approach of the excitable tissue with rotational anisotropy in its properties to a pseudoisotropic tissue.     

Anisotropy of ultrasonic velocity and elastic properties in normal human myocardium.

Measurements of ultrasonic quasilongitudinal velocity were made in the muscle fiber plane of excised human myocardium. Multiple adjacent planes across the left ventricular wall were interrogated to assess the transmural dependence of velocity. For each measurement plane, data were obtained in 2-deg increments through the full 360 deg relative to the myofibers. An approximate 1.3% magnitude of anisotropy was observed with maximum velocity along the muscle fibers and minimum velocity perpendicular to the muscle fibers. The known transmural shift in myofiber orientation was evidenced in the anisotropy of velocity as angular shifts between plots obtained from adjacent transmural planes within the same specimen. Measured values of velocity and density were used to estimate the effective C33 and C11 elastic constants of a thin layer of normal myocardium.     

Differentiation between acutely ischemic myocardium and zones of completed infarction in dogs on the basis of frequency-dependent backscatter

The goal of this work was to determine whether the frequency dependence of apparent backscatter coefficient (not corrected for attenuation within the myocardium) could differentiate completed, remote infarction from acute myocardial injury in vivo. Myocardial infarcts were produced in six dogs by coronary artery occlusion. One to 12 months later, acute ischemic injury was induced in each dog by ligation of a coronary artery that supplied a region of myocardium adjacent to the established infarct. Infarct, ischemic, and normal regions were interrogated with a 5-MHz, circular, 0.5-in. diam, broadband, focused, piezoelectric transducer mounted in a water-filled stand-off device placed against the exposed, beating heart. Apparent backscatter coefficients were measured over the range of frequencies from 3-7 MHz. The frequency dependence was obtained from the slope of log apparent backscatter coefficient versus log frequency. No significant difference in frequency dependence was found between normal and acutely ischemic myocardium for periods of up to 2 h of ischemia. In contrast, frequency dependence in regions of remote infarct (1.8 +/- 0.1, mean +/- standard error) was significantly lower than that in acutely ischemic or nonischemic regions (2.3 +/- 0.1) (p less than 0.01). These results suggest that remote myocardial infarction can be differentiated from acutely injured but still potentially salvageable myocardium in vivo on the basis of the frequency dependence of backscatter.     

Anisotropy of the ultrasonic backscatter of myocardial tissue: I. Theory and measurements in vitro

This research addresses the variations in the ultrasonic backscatter from specimens consisting of a suspension of approximately aligned cylindrical scatterers in a fluid medium as a function of the angle of propagation in the sample. Predictions of the angular dependence of backscatter based on the time-domain Born approximation described by Rose and Richardson [J. H. Rose and J. M. Richardson, J. Nondestr. Eval. 3, 45-53 (1982)] were compared with experimental measurements of the backscatter from both tissue-mimicking phantoms consisting of graphite fibers suspended in gelatin and from canine myocardial tissue. The angular dependence of the backscatter was predicted and measured to be maximum for propagation perpendicular to the cylinder axes and minimum for propagation parallel to the axes. Maximum to minimum (i.e., perpendicular to parallel) changes in the integrated backscatter were predicted to be between 5 and 10 dB in the phantom. The corresponding quantity measured in both the phantom and in canine myocardial tissue was approximately 6 dB.     

Spectral estimation for characterization of acoustic aberration

Spectral estimation based on acoustic backscatter from a motionless stochastic medium is described for characterization of aberration in ultrasonic imaging. The underlying assumptions for the estimation are: The correlation length of the medium is short compared to the length of the transmitted acoustic pulse, an isoplanatic region of sufficient size exists around the focal point, and the backscatter can be modeled as an ergodic stochastic process. The motivation for this work is ultrasonic imaging with aberration correction. Measurements were performed using a two-dimensional array system with 80 x 80 transducer elements and an element pitch of 0.6 mm. The f number for the measurements was 1.2 and the center frequency was 3.0 MHz with a 53% bandwidth. Relative phase of aberration was extracted from estimated cross spectra using a robust least-mean-square-error method based on an orthogonal expansion of the phase differences of neighboring wave forms as a function of frequency. Estimates of cross-spectrum phase from measurements of random scattering through a tissue-mimicking aberrator have confidence bands approximately +/- 5 degrees wide. Both phase and magnitude are in good agreement with a reference characterization obtained from a point scatterer.     

A two-dimensional CVIB imaging system with a speckle tracking algorithm

Quantitative ultrasound tissue characterization based on integrated backscatter (IB) has shown great potential in detecting myocardial ischemia. The magnitude of the cyclic variation in IB (CVIB) has been considered one promising parameter in assessing regional myocardial contractile performance. This lab has previously developed a novel ultrasonic fusion imaging method based on CVIB. However, the major problem for clinical applications of this technique is that the myocardial tissue could not be tracked effectively without cardiologist’s intervention. This paper introduced a speckle tracking method into the CVIB-weighted imaging system, called speckle tracking algorithm with adaptive window size (STAWAWS), to track myocardial tissue particle automatically. This method provides a way to obtain the particle’s positions frame by frame in a series of B-mode images. Then using the RF signals according to the particle’s positions the IB curve can be calculated to produce CVIB value. The method was applied on the experimental and clinical data cases’s analysis. The results of dog’s data processing showed that this method could eliminate the misunderstanding of myocardial ischemia especially near the endocardium. The results of clinical data suggested that this method had clinical significance in detecting ischemic myocardium. Though the CVIB-weighted images obtained by the use of this auto-tracking method can improve the accuracy of detecting myocardial ischemia, it is not real-time analysis and the clinical data cases are not sufficient. Further clinical validation is still needed in the future’ work.     

Extended three-dimensional impedance map methods for identifying ultrasonic scattering sites

The frequency-dependent ultrasound backscatter from tissues contains information about the microstructure that can be quantified. In many cases, the anatomic microstructure details responsible for ultrasonic scattering remain unidentified. However, their identification would lead to potentially improved methodologies for characterizing tissue and diagnosing disease from ultrasonic backscatter measurements. Recently, three-dimensional (3D) acoustic models of tissue microstructure, termed 3D impedance maps (3DZMs), were introduced to help to identify scattering sources [J. Mamou, M. L. Oelze, W. D. O'Brien, Jr., and J. F. Zachary, "Identifying ultrasonic scattering sites from 3D impedance maps," J. Acoust. Soc. Am. 117, 413-423 (2005)]. In the current study, new 3DZM methodologies are used to model and identify scattering structures. New processing procedures (e.g., registration, interpolations) are presented that allow more accurate 3DZMs to be constructed from histology. New strategies are proposed to construct scattering models [i.e., form factor (FF)] from 3DZMs. These new methods are tested on simulated 3DZMs, and then used to evaluate 3DZMs from three different rodent tumor models. Simulation results demonstrate the ability of the extended strategies to accurately predict FFs and estimate scatterer properties. Using the 3DZM methods, distinct FFs and scatterer properties were obtained for each tumor examined.     

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.     

Anisotropy of the ultrasonic attenuation in soft tissues: measurements in vitro

This study was designed to measure the ultrasonic attenuation within phantoms and tissue samples over a broad bandwidth and at many angles of incidence with respect to intrinsic orientations in order to elucidate both the frequency and angular dependence of the attenuation coefficient. Significant angular dependence, or anisotropy, of the attenuation was observed in canine myocardium (maximum to minimum ratio: 2.2 to 1) and a tissue mimicking phantom of oriented graphite fibers in gelatin (max to min: 2 to 1). In control studies, insignificant anisotropy was observed in the attenuation in canine liver samples and phantoms with graphite powder suspended in gelatin. Comparisons of the magnitude of variations of the oriented-fiber phantom to that predicted by a viscous relative motion model are presented.     

Analysis of microstructural alterations of normal and pathological breast tissue in vivo using the AR cepstrum

The mean scatterer spacing is considered to be an important parameter for describing ultrasonic scattering and characterization of biological tissue. Autoregressive models are widely used in parametric techniques for spectral estimation. In this paper, we describe the results of a careful examination of the mean scatterer spacing parameter in normal and pathological breast tissues in vivo using the autoregressive cepstrum. Our experimental results carried out at 4.5 MHz using weakly focused pulse-echo single element transducer show that the mean scatterer spacing in normal breast tissues in vivo is 1.25+/-0.21 mm whereas in several pathological breast tissues, it is between 0.82+/-0.10 and 1.09+/-0.07 mm. These results indicate good correlation with microstructure of breast tissue characterization, and hence the AR cepstrum holds promise that it could be used as an effective method for signal analysis of ultrasonic scattering and characterization of breast tissues scatterers.     

High frequency ultrasonic characterization of human vocal fold tissue

Recently, endolaryngeal sonography at frequencies ranging from 10 to 30 MHz has been found to be useful in diagnosing diseases of the vocal folds (VFs). However, image resolution can be further improved by ultrasound at higher frequencies, necessitating the measurement of high-frequency acoustic properties of VF tissue. The ultrasonic parameters of integrated backscatter, sound velocity, and frequency-dependent attenuation coefficient were measured in both the lamina propria (LP) and vocalis muscle (VM) of human VFs using a 47 MHz high-frequency ultrasonic transducer. The integrated backscatter was -173.44+/-6.14 (mean+/-s.d.) and -195.13+/-3.58 dB in the LP and VM, respectively, the sound velocity was 1667.68+/-44.9 and 1595.07+/-39.33 ms, and the attenuation coefficient at 47 MHz was 8.28+/-1.72 and 7.17+/-1.30 dBmm. The difference between these ultrasonic parameters may be attributed to variations in the structure and fiber concentrations in VF tissue. These results could serve as a useful clinical reference for the further development of high-frequency ultrasound devices for endolarynx sonography applications.     

Biochemical and anisotropical properties of tendons

Tendons are formed by dense connective tissue composed of an abundant extracellular matrix (ECM) that is constituted mainly of collagen molecules, which are organized into fibrils, fibers, fiber bundles and fascicles helicoidally arranged along the largest axis of the tendon. The biomechanical properties of tendons are directly related to the organization of the collagen molecules that aggregate to become a super-twisted cord. In addition to collagen, the ECM of tendons is composed of non-fibrillar components, such as proteoglycans and non-collagenous glycoproteins. The capacity of tendons to resist mechanical stress is directly related to the structural organization of the ECM. Collagen is a biopolymer and presents optical anisotropies, such as birefringence and linear dichroism, that are important optical properties in the characterization of the supramolecular organization of the fibers. The objective of this study was to present a review of the composition and organization of the ECM of tendons and to highlight the importance of the anisotropic optical properties in the study of alterations in the ECM.     

Anisotropy of ultrasonic propagation and scattering properties in fresh rat skeletal muscle in vitro

The anisotropy of frequency-dependent backscatter coefficient, attenuation, and speed of sound is assessed in fresh rat skeletal muscle within 5 h post-mortem. Excised rat semimembranosus and soleus muscles are measured in 37 degrees C Tyrode solution, with the muscle fibers at 90 degrees and 45 degrees orientations to the incident sound beam. Reflected and through transmission signals from either a 6- or 10-MHz focused transducer give frequency dependent information in the 4-14 MHz range. The attenuation coefficient in each muscle is consistently a factor of 2.0 +/- 0.4 lower for propagation perpendicular to the fibers than at 45 degrees, whereas speed of sound shows a much milder anisotropy, and is slightly faster for the 90 degrees orientation. The largest anisotropy is seen in the backscatter coefficient, most notably in the semimembranosus where the magnitude at 90 degrees is over an order of magnitude greater than at 45 degrees, with the frequency dependence in both cases giving a power law between 1.5 and 2.0.     

Characterization of dense bovine cancellous bone tissue microstructure by ultrasonic backscattering using weak scattering models

A weak scattering model was proposed for the ultrasonic frequency-dependent backscatter in dense bovine cancellous bone, using two autocorrelation functions to describe the medium: one with discrete homogeneities (spherical distribution of equal spheres) and another, which considers tissue as an inhomogeneous continuum (densely populated medium). The inverse problem to estimate trabecular thickness of bone tissue has been addressed. A combination of the two autocorrelation functions was required to closely approximate the backscatter from bovine bone with various microarchitecture, given that the shape of trabeculae ranges from a rodlike to a platelike shape. Because of the large variation in trabecular thickness, both at an intraspecimen and an interspecimen level, thickness distributions for individual trabeculae for each bone specimen were obtained, and dominant trabecular sizes were determined. Comparison of backscatter measurements to theoretical predictions indicated that there were more than one dominant trabecular sizes that scatter sound for most specimens. Linear regression, performed between dominant trabecular thickness and estimated correlation length, showed significant linear correlation (R(2)=0.81). Attenuation due to scattering by a continuous distribution of scatterers was predicted to be linear over a frequency range from 0.3 to 0.9 MHz, suggesting a possibility that scattering may be a significant source of attenuation.