Motion compensation for intravascular ultrasound palpography

Rupture of vulnerable plaques in coronary arteries is the major cause of acute coronary syndromes. Most vulnerable plaques consist of a thin fibrous cap covering an atheromous core. These plaques can be identified using intravascular ultrasound (IVUS) palpography, which measures radial strain by cross-correlating RF signals at different intraluminal pressures. Multiple strain images (i.e., partial palpograms) are averaged per heart cycle to produce a more robust compounded palpogram. However, catheter motion due to cardiac activity causes misalignment of the RF signals and thus of the partial palpograms, resulting in less valid strain estimates. To compensate for in-plane catheter rotation and translation, we devised four methods based on block matching. The global rotation block matching (GRBM) and contour mapping (CMAP) methods measure catheter rotation, and local block matching (LBM) and catheter rotation and translation (CRT) estimate displacements of local tissue regions. These methods were applied to nine in vivo pullback acquisitions, made with a 20 MHz phased-array transducer. We found that all these methods significantly increase the number of valid strain estimates in the partial and compounded palpograms (P < 0.008). The best method, LBM, attained an average increase of 17% and 15%, respectively. Implementation of this method should improve the information coming from IVUS palpography, leading to better vulnerable plaque detection.     

Decorrelation characteristics of transverse blood flow along an intravascular array catheter

In recent years, a new method to measure transverse blood flow based on the decorrelation of radio frequency (RF) signals has been introduced. In this paper, we investigated the decorrelation characteristics of transverse blood flow measurement using an intravascular ultrasound (IVUS) array catheter by means of computer modeling. Blood was simulated as a collection of randomly located point scatterers. Moving this scattering medium transversally across the acoustical beam represented flow. First-order statistics were evaluated, and the signal-to-noise ratio from the signals was measured. The correlation coefficient method was used to present the results. The decorrelation patterns for RF and for RF-envelope signals were studied. The decorrelation patterns from the RF signals were in good agreement with those obtained from theoretical beam profiles. This agreement suggests that the decorrelation properties of an IVUS array catheter for measuring quantitative transverse blood flow can be assessed by measuring the ultrasound beam. A line of point scatterers, moved transversally across the acoustical beam (line spread function), can determine this decorrelation behaviour.     

A deconvolution filter for improvement of time-delay estimation in elastography

In elastography, tissue under investigation is compressed, and the resulting strain is estimated from the gradient of displacement estimates. Therefore, it is important to accurately estimate the displacements (time-delay) for good quality elastograms. A principal source of error in time-delay estimation in elastography is the decorrelation of the echo signal due to tissue compression (decorrelation noise). Temporal stretching of the postcompression signals has been shown to reduce the decorrelation noise at small strains. In this article, we present a deconvolution filter that reduces the decorrelation even further when applied in conjunction with signal stretching. The performance of the proposed filter is evaluated using simulated data.     

Decorrelation-based blood flow velocity estimation: effect of spread of flow velocity,linear flow velocity gradients,and parabolic flow

In recent years, a new method to measure transverse blood flow, based on the decorrelation of the radio frequency (RF) signals has been developed. In this paper, we investigated the influence of nonuniform flow on the velocity estimation. The decorrelation characteristics of transverse blood flow using an intravascular ultrasound (IVUS) array catheter are studied by means of computer modeling. Blood was simulated as a collection of randomly located point scatterers; moving this scattering medium transversally across the acoustical beam represented flow. First-order statistics were evaluated, and the signal-to-noise ratio from the signals were measured. The correlation coefficient method was used to present the results. Three velocity profiles were simulated: random spread of blood-flow velocity, linear blood-flow velocity gradient, and parabolic blood-flow. Radio frequency and envelope signals were used to calculate the decorrelation pattern. The results were compared to the mean decorrelation pattern for plug blood-flow. The RF signals decorrelation patterns were in good agreement with those obtained for plug blood flow. Envelope decorrelation patterns show a close agreement with the one for plug blood flow. For axial blood flow, there is a discrepancy between decorrelation patterns. The results presented here suggest that the decorrelation properties of an IVUS array catheter for measuring quantitative transverse blood flow probably will not be affected by different transverse blood-flow conditions     

Enhancement of echo-signal correlation in elastography using temporal stretching

Echo-signal decorrelation due to tissue compression is a significant source of error in tissue displacement estimates obtained using crosscorrelation. Tissue displacement estimates are used to compute strain values for imaging the elasticity of biological soft tissues. The correlation coefficient between the pre- and post-compression echo rf signals reduces rapidly with signal decorrelation due to increased compression. Miniscule reductions in the value of the correlation coefficient can have a significant impact on the performance of the strain estimator as illustrated by the strain filter. Reducing the rate of signal decorrelation using temporal stretching (which improves the value of the correlation coefficient), significantly improves the performance of the strain filter. The reduction in the rate of signal decorrelation with the subsequent increase in the correlation coefficient using temporal stretching is discussed in this paper. Theoretical, simulation and experimental results quantify the enhancement in the value of the correlation coefficient attained with temporal stretching.     

Quantitative blood speed imaging with intravascular ultrasound

Previously, we presented a method of real-time arterial color flow imaging using an intravascular ultrasound (IVUS) imaging system, where real-time RF A-scans were processed with an FIR (finite-impulse response) filter bank to estimate relative blood speed. Although qualitative flow measurements are clinically valuable, realizing the full potential of blood flow imaging requires quantitative flow speed and volume measurements in real time. Unfortunately, the rate of RF echo-to-echo decorrelation is not directly related to scatterer speed in a side-looking IVUS system because the elevational extent of the imaging slice varies with range. Consequently, flow imaging methods using any type of decorrelation processing to estimate blood speed without accounting for spatial variation of the radiation pattern will have estimation errors that prohibit accurate comparison of speed estimates from different depths. The FIR filter bank approach measures the rate of change of the ultrasound signal by estimating the slow-time spectrum of RF echoes. A filter bank of M bandpass filters is applied in parallel to estimate M components of the slow-time DFT (discrete Fourier transform). The relationship between the slow-time spectrum, aperture diffraction pattern, and scatterer speed is derived for a simplified target. Because the ultimate goal of this work is to make quantitative speed measurements, we present a method to map slow time spectral characteristics to a quantitative estimate. Results of the speed estimator are shown for a simulated circumferential catheter array insonifying blood moving uniformly past the array (i.e., plug flow) and blood moving with a parabolic profile (i.e., laminar flow)     

An adaptive strain estimator for elastography

Elastography is based on the estimation of strain due to applied tissue compression. In conventional elastography, strain is computed from the gradient of the displacement estimates between gated pre- and postcompression echo signals. Gradient-based estimation methods are known to be susceptible to noise. In elastography, in addition to the electronic noise, a principal source of estimation error is the decorrelation of the echo signal as a result of tissue compression (decorrelation noise). Temporal stretching of postcompression signals previously was shown to reduce the decorrelation noise. In this paper, we introduce a novel estimator that uses the stretch factor itself as an estimator of the strain. It uses an iterative algorithm that adaptively maximises the correlation between the pre- and postcompression echo signals by appropriately stretching the latter. We investigate the performance of this adaptive strain estimator using simulated and experimental data. The estimator has exhibited a vastly superior performance compared with the conventional gradient-based estimator.     

Motion artifact reduction for IVUS-based thermal strain imaging

Thermal strain imaging (TSI) using intravascular ultrasound (IVUS) has the potential to identify lipid pools within rupture-prone arterial plaques and serve as a valuable supplement to current IVUS systems in diagnosing acute coronary syndromes. The major challenge for in vivo application of TSI will be cardiac motion, including bulk motion and tissue deformation. Simulations based on an artery model, including a lipid-filled plaque, demonstrate that effective bulk motion compensation can be achieved within a certain motion range using spatial interpolation. We also propose a practical imaging scheme to minimize mechanical strains caused by tissue deformation based on a linear least squares fitting strategy. This scheme was tested on clinical data by artificially superimposing thermal displacements corresponding to different temperature rises. Results suggest a 1-2 degrees C temperature rise is required to detect lipids in an atherosclerotic plaque in vivo.     

Estimates of echo correlation and measurement bias in acoustic radiation force impulse imaging

Acoustic radiation force impulse (ARFI) imaging is a novel imaging modality in which pulses from a diagnostic ultrasound scanner are used to displace tissue and track its motion. The region displaced has lateral and elevational dimensions of similar scale to the ultrasound beams used to track the motion. Therefore, there is a range of tissue displacements present within the tracking beam, leading to decorrelation of the echo signal. Expressions are derived for the expected value of the displacement estimate and the cross-correlation at the expected displacement. Numerical simulations confirm the analytical model.     

Echo decorrelation from displacement gradients in elasticity and velocity estimation

Several ultrasonic techniques for the estimation of blood velocity, tissue motion and elasticity are based on the estimation of displacement through echo time-delay analysis. A common assumption is that tissue displacement is constant within a short observation time used for time delay estimation (TDE). The precision of TDE is mainly limited by noise sources corrupting the echo signals. In addition to electronic and quantization noise, a substantial source of TDE error is the decorrelation of echo signals because of displacement gradients within the observation time. The authors present a theoretical model that describes the mean changes of the crosscorrelation function as a function of observation time and displacement gradient. The gradient is assumed to be small and uniform within the observation time; the decorrelation introduced by the lateral and elevational displacement components is assumed to be small compared with the decorrelation caused by the axial component. The decorrelation model predicts that the expected value of the crosscorrelation function is a low-pass filtered version of the autocorrelation function (i.e., the crosscorrelation obtained without gradients). The filter is a function of the axial gradient and the observation time. This theoretical finding is corroborated experimentally. Limitations imposed by decorrelation in displacement estimation and potential uses of decorrelation in medical ultrasound are discussed.     

A fundamental limit on delay estimation using partially correlatedspeckle signals

Delay estimation is used in ultrasonic imaging to estimate blood or soft tissue motion, to measure echo arrival time differences for phase aberration correction, and to estimate displacement for tissue elasticity measurements. In each of these applications delay estimation is performed using speckle signals which are at least partially decorrelated relative to one another. Delay estimates which utilize such data are subject to large errors known as false peaks and smaller magnitude errors known as jitter. While false peaks can sometimes be removed through nonlinear processing, jitter errors place a fundamental limit on the performance of delay estimation techniques. The authors apply the Cramer-Rao Lower Bound to derive an analytical expression which predicts the magnitude of jitter errors incurred when estimating delays using radio frequency (RF) data from speckle targets. The analytical expression presented includes the effects of signal decorrelation due to physical processes, corruption by electronic noise, and a number of other factors. Simulation results are presented which show that the performance of the normalized cross correlation algorithm closely matches theoretical predictions. These results indicate that for poor signal to noise ratios (0 dB) a small improvement in signal to noise ratio can dramatically reduce jitter magnitude. At high signal to noise ratios (30 dB) small amounts of signal decorrelation can significantly increase the magnitude of jitter errors     

Visualizing the radial and circumferential strain distribution within vessel phantoms using synthetic-aperture ultrasound elastography

Noninvasive elastography (NIVE) produces elastograms that are difficult to interpret because NIVE visualizes strain in the transducer coordinate system. In this paper, we hypothesized that transforming normal and shear strain elastograms to the vessel coordinate system will produce better strain elastograms. To corroborate this hypothesis, we acquired synthetic-aperture (SA) ultrasound data from simulated and physical vessel phantoms. In both studies, SA echo frames were reconstructed from data acquired with a sparse transducer array. The simulation study was performed with homogeneous and heterogenous phantoms, but in the experimental study we used a modified ultrasound scanner to acquire SA data from homogeneous (n = 1) and heterogeneous (n = 3) vessel phantoms. Axial and lateral displacements were estimated by performing two-dimensional cross-correlation analysis on the beamformed RF echo frames. We generated radial and circumferential strain elastograms by transforming normal and shear strain elastograms to the vessel coordinate system. The results revealed: 1) radial and circumferential strain elastograms acquired from simulated data had a relative root mean squared error on the order of 0.1%; 2) experimentally acquired radial and circumferential strain elastograms had elastographic contrast-to-noise ratio (CNRe) between 10 and 40 dB, and elastographic signal-to-noise ratio (SNRe) between 10 and 35 dB, depending on the number of active transmission elements employed during imaging; 3) radial and circumferential strain elastograms produced with fewer than 8 active transmission elements were inferior to those computed with a greater number of active elements; and 4) plaques were evident in the strain elastograms, except in those obtained with the sparsest transducer array. This study demonstrated that a syntheticaperture ultrasound system could visualize radial and circumferential strain noninvasively.     

Frequency-domain-based strain estimation and high-frame-rate imaging for quasi-static elastography

In freehand elastography, quasi-static tissue compression is applied through the ultrasound probe, and the corresponding axial strain is estimated by calculating the time shift between consecutive echo signals. This calculation typically suffers from a poor signal-to-noise ratio or from the decorrelation between consecutive echoes resulting from an erroneous axial motion impressed by the operator. This paper shows that the quality of elastograms can be improved through the integration of two distinct techniques in the strain estimation procedure. The first technique evaluates the displacement of the tissue by analyzing the phases of the echo signal spectra acquired during compression. The second technique increases the displacement estimation robustness by averaging multiple displacement estimations in a high-frame-rate imaging system, while maintaining the typical elastogram frame-rate. The experimental results, obtained with the Ultrasound Advanced Open Platform (ULA-OP) and a cyst phantom, demonstrate that each of the proposed methods can independently improve the quality of elastograms, and that further improvements are possible through their combination.     

Estimation of displacement location for enhanced strain imaging

Ultrasonic strain imaging usually begins with displacement estimates computed using finite-length sections of RF ultrasound signals. Amplitude variations in the ultrasound are known to perturb the location at which the displacement estimate is valid. If this goes uncorrected, it is a significant source of estimation noise, which is amplified when displacement fields are converted into strain images. We present a study of this effect based on theoretical analysis and practical experiments. A correction method based on the analysis is tested on phase zero and correlation coefficient strain imaging, and compared to the amplitude compression techniques of earlier studies. We also test adaptive strain estimation to provide a benchmark, but the performance of our new method matches or surpasses this benchmark under normal scanning conditions. Furthermore, the new correction is suitable for real time applications owing to its extreme computational simplicity.     

Strain compounding: a new approach for speckle reduction

A new compounding technique for reducing speckle brightness variations is proposed. This method exploits the decorrelation between signals under different strain states. The different strain states can be created using externally applied forces such as the ones used in sonoelastography. Such forces produce three-dimensional tissue motion. By correcting only the in-plane (i.e., axial and lateral) motion, the images under different strain states have similar characteristics except for speckle appearance caused by the uncorrected out-of-plane (i.e., elevational) motion. Additional speckle decorrelation is also introduced through tissue motion correction caused by the change of effective in-plane sample volume geometry. Therefore, these images can be combined for speckle reduction with less degradation in in-plane spatial resolution than conventional approaches. In this paper, three-dimensional tissue motion under various strain conditions were simulated. It was found that significant speckle decorrelation existed at strains achievable in some clinical situations. Experiments were also conducted to test efficacy of this approach. Pulse-echo data from a gelatin-based phantom were acquired using a 5-MHz, single crystal transducer, and both conventional and compound B-mode images were formed. Results indicated that speckle brightness variations were reduced, and detectability of low contrast objects was enhanced. Performance limitations and fundamental differences between the proposed technique and existing techniques are discussed     

Estimation and reduction of decorrelation effect due to tissue lateral displacement in elastography

In cross-correlation based elastography, the quality of the strain image is degraded by the distortion of echo waveforms due to tissue axial and lateral displacement. To study the effects of tissue lateral displacement on echo decorrelation, a tissue axial stretching model is developed and a concept called correlation signal-to-noise ratio (CSNR) is introduced to quantify the decorrelation effect due to tissue lateral displacement. A computer simulation based on the tissue stretching model is carried out to study the influence of several important elastographic parameters on echo decorrelation due to tissue lateral displacement. Finally, guided by the CSNR concept, a 2-D spatial comprehensive cross-correlation method is proposed to reduce the decorrelation noise. Results indicate that CSNR can be used as a quality indicator of elastography and the 2-D spatial comprehensive cross-correlation method can effectively reduce the decorrelation noise while slightly decreasing the lateral resolution of the strain image     

The impact of various imaging parameters on ultrasonic displacementand velocity estimates

The accuracy of displacement and velocity data in ultrasonic motion detection systems depends on a combination of ultrasonic imaging parameters. These include magnitude and direction of target motion, target region dimensions, scattering media, ultrasonic frequency of interrogation, digital sampling frequency, and signal type (envelope detected or RF). Because the impact of scattering media in particular has heretofore received little or no attention, we provide experimental results and computer analysis to evaluate the influence of different scattering media on the accuracy of ultrasonic displacement and velocity estimates using porcine liver, porcine muscle, and woolen sea sponge samples. Our experimental results show that for identical target dimensions and displacements, the accuracy of ultrasonic displacement and velocity estimates in porcine muscle samples are substantially higher than in porcine liver samples. Analysis of experimentally derived autocovariance curves for each tissue type indicates that the improvement in accuracy for muscle samples is not, in fact, due to differences in the speckle characteristics for each tissue type. The improvement is attributed to the presence of well-defined and resolvable image structures from muscle and the lack of such resolvable structures in porcine liver tissue. We provide a summary of the factors impacting ultrasonic displacement and velocity measurements, and discuss why and how a combination of one or more imaging parameters affects these measurements     

A time-efficient and accurate strain estimation concept for ultrasonic elastography using iterative phase zero estimation

In ultrasonic elastography, the exact estimation of temporal displacements between two signals is the key to estimating strain. An algorithm was previously proposed that estimates these displacements using phase differences of the corresponding base-band signals. A major advantage of these algorithms compared with correlation techniques is the computational efficiency. In this paper, an extension of the algorithm is presented that iteratively takes into account the time shifts of the signals to overcome the problems of aliasing and accuracy in the estimation of the phase shift. Thus, it can be proven that the algorithm is equivalent to the search of the maximum of the correlation function. Furthermore, a robust logarithmic compression is proposed that only compresses the envelope of the signal. This compression does not introduce systematic errors and significantly reduces decorrelation noise. The resulting algorithm is a computationally simple and very fast alternative to conventional correlation techniques, and the accuracy of strain images is improved.     

Speckle tracking methods for ultrasonic elasticity imaging using short-time correlation

In ultrasound elasticity imaging, strain decorrelation is a major source of error in displacements estimated using correlation techniques. This error can be significantly decreased by reducing the correlation kernel. Additional gains in signal-to-noise ratio (SNR) are possible by filtering the correlation functions prior to displacement estimation. Tradeoffs between spatial resolution and estimate variance are discussed, and estimation in elasticity imaging is compared to traditional time-delay estimation. Simulations and experiments on gel-based phantoms are presented. The results demonstrate that high resolution, high SNR strain estimates can be computed using small correlation kernels (on the order of the autocorrelation width of the ultrasound signal) and correlation filtering.     

Theoretical quality assessment of myocardial elastography with in vivo validation

Myocardial elastography (ME), a radio frequency (RF)-based speckle tracking technique with one-dimensional (1-D) cross correlation and novel recorrelation methods in a 2-D search was proposed to estimate and fully image 2-D transmural deformation field and to detect abnormal cardiac function. A theoretical framework was first developed in order to evaluate the performance of 2-D myocardial elastography based on a previously developed 3-D finite-element model of the canine left ventricle. A normal (control) and an ischemic (left-circumflex, LCx) model, which more completely represented myocardial deformation than a kinematic model, were considered. A 2-D convolutional image formation model was first used to generate RF signals for quality assessment of ME in the normal and ischemic cases. A 3-D image formation model was further developed to investigate the effect of the out-of-plane motion on the 2-D, in-plane motion estimation. Both orthogonal, in-plane displacement components (i.e., lateral and axial) between consecutive RF frames were iteratively estimated. All the estimated incremental 2-D displacements from end-diastole (ED) to end-systole (ES) were then accumulated to acquire the cumulative 2-D displacements, which were further used to calculate the cumulative 2-D systolic finite strains. Furthermore, the cumulative systolic radial and circumferential strains, which were angle- and frame-rate independent, were obtained from the 2-D finite-strain components and imaged in full view to detect the ischemic region. We also explored the theoretical understanding of the limitations of our technique for the accurate depiction of disease and validated it in vivo against tagged magnetic resonance imaging (tMRI) in the case of a normal human myocardium in a 2-D short-axis (SA) echocardiographic view. The theoretical framework succeeded in demonstrating that the 2-D myocardial elastography technique was a reliable tool for the complete estimation and depiction of the in-plane myocardial deformation field as well as for accurate identification of pathological mechanical function using established finite-element, left-ventricular canine models. In a preliminary study, the 2-D myocardial elastography was shown capable of imaging myocardial deformation comparable to equivalent tMRI estimates in a clinical setting.