Estimation of mean frequency and variance of ultrasonic Doppler signal by using second-order autoregressive model

In order to estimate the mean frequency and variance of the diagnostic ultrasound Doppler signal in the presence of clutter noise, a new estimator using a second-order autoregressive (AR) model, called the AR estimator, is proposed. The sampled signal that contains information of both the Doppler signal and clutter is described by the second-order AR model with two poles. The mean frequency and variance of a unidirectional Doppler signal can be estimated, respectively, from the phase and the magnitude of the pole, with larger phase between the two poles. If the clutter is not completely rejected, all conventional estimators, including the autocorrelation (AC) estimator, result in erroneous estimations for the mean frequency and variance of the Doppler signal, whereas the AR estimator gives an accurate estimation. In the absence of clutter, however, the performance of both the AC and AR estimators are similar. If the blood flows in both directions in a sample volume and the clutter is rejected to the extent that it no longer obscures the Doppler signal, the proposed method can estimate simultaneously the mean frequencies and variances of both the forward and reverse blood flows. The performance of the proposed AR estimator was compared with that of the AC estimator by computer simulations and experiments, and it was found that when the number of available sampled data is small, the AR estimator does not require the use of a clutter filter, which simplifies Doppler signal detection.     

Real-time adaptive clutter rejection filtering in color flow imaging using power method iterations

We propose a new algorithm for real-time, adaptive-clutter-rejection filtering in ultrasound color flow imaging (CFI) and related techniques. The algorithm is based on regression filtering using eigenvectors of the signal correlation matrix as a basis for representing clutter, a method that previously has been considered too computationally demanding for real-time processing in general CFI applications. The data acquisition and processing scheme introduced allows for a more localized sampling of the clutter statistics and, therefore, an improved clutter attenuation for lower filter orders. By using the iterative power method technique, the dominant eigenvalues and corresponding eigenvectors of the correlation matrix can be estimated efficiently, rendering real-time operation feasible on desktop computers. A new adaptive filter order algorithm is proposed that successfully estimates the proper dimension of the clutter basis, previously one of the major drawbacks of this clutter-rejection technique. The filter algorithm performance and computational demands has been compared to that of conventional clutter filters. Examples have been included which confirms that, by adapting the clutter-rejection filter to estimates of the clutter-signal statistics, improved attenuation of the clutter signal can be achieved in normal as well as more excessive cases of tissue movement and acceleration.     

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.     

Maximum likelihood blood velocity estimator incorporating properties of flow physics

The aspect of correlation among the blood velocities in time and space has not received much attention in previous blood velocity estimators. The theory of fluid mechanics predicts this property of the blood flow. Additionally, most estimators based on a cross-correlation analysis are limited on the maximum velocity detectable. This is due to the occurrence of multiple peaks in the cross-correlation function. In this study a new estimator (CMLE), which is based on correlation (C) properties inherited from fluid flow and maximum likelihood estimation (MLE), is derived and evaluated on a set of simulated and in vivo data from the carotid artery. The estimator is meant for two-dimensional (2-D) color flow imaging. The resulting mathematical relation for the estimator consists of two terms. The first term performs a cross-correlation analysis on the signal segment in the radio frequency (RF)-data under investigation. The flow physic properties are exploited in the second term, as the range of velocity values investigated in the cross-correlation analysis are compared to the velocity estimates in the temporal and spatial neighborhood of the signal segment under investigation. The new estimator has been compared to the cross-correlation (CC) estimator and the previously developed maximum likelihood estimator (MLE). The results show that the CMLE can handle a larger velocity search range and is capable of estimating even low velocity levels from tissue motion. The CC and the MLE produce incorrect velocity estimates due to the multiple peaks, when the velocity search range is increased above the maximum detectable velocity. The root-mean square error (RMS) on the velocity estimates for the simulated data is on the order of 7 cm/s (14%) for the CMLE, and it is comparable to the RMS for the CC and the MLE. When the velocity search range is set to twice the limit of the CC and the MLE, the number of incorrect velocity estimates are 0, 19.1, and 7.2% for the CMLE, CC, and MLE, respectively. The ability to handle a larger search range and estimating low velocity levels was confirmed on in vivo data.     

Clutter filter design for ultrasound color flow imaging

For ultrasound color flow images with high quality, it is important to suppress the clutter signals originating from stationary and slowly moving tissue sufficiently. Without sufficient clutter rejection, low velocity blood flow cannot be measured, and estimates of higher velocities will have a large bias. The small number of samples available (8 to 16) makes clutter filtering in color flow imaging a challenging problem. In this paper, we review and analyze three classes of filters: finite impulse response (FIR), infinite impulse response (IIR), and regression filters. The quality of the filters was assessed based on the frequency response, as well as on the bias and variance of a mean blood velocity estimator using an autocorrelation technique. For FIR filters, the frequency response was improved by allowing a non-linear phase response. By estimating the mean blood flow velocity from two vectors filtered in the forward and backward direction, respectively, the standard deviation was significantly lower with a minimum phase filter than with a linear phase filter. For IIR filters applied to short signals, the transient part of the output signal is important. We analyzed zero, step, and projection initialization, and found that projection initialization gave the best filters. For regression filters, polynomial basis functions provide effective clutter suppression. The best filters from each of the three classes gave comparable bias and variance of the mean blood velocity estimates. However, polynomial regression filters and projection-initialized IIR filters had a slightly better frequency response than could be obtained with FIR filters     

Clutter metric based on the Cramer-Rao lower bound on automatic target recognition

This is a performance evaluation on the implementation of the Cramer-Rao lower bound (CRLB) for background clutter measurement on automatic target recognition (ATR). In essence the background clutter evaluation problem for ATR is consistent with the deterministic parameter estimation problem. Thus, useful concepts and theories of deterministic parameter estimation can be introduced into the investigation of background clutter. In this paper, the CRLB is employed as a metric for clutter measurement. Requirements needed for this application are analyzed, and the approach for obtaining the CRLB of a scene image is produced. The flexibility of the CRLB metric is analyzed. Discussion and comparison are made on the relationship between the CRLB metric and the Sims signal-to-clutter metric. Finally, we illustrate how this metric defines the potential for false alarms by determining the correspondence level between a target and background through the application of the CRLB.     

Clutter rejection filters in color flow imaging: a theoretical approach

A general class of linear clutter rejection filters is described, covering the commonly used filter types including FIR/IIR filters with linear initialization, as well as regression filters, where the clutter component is estimated by least square curve fitting. The filter can be described by a complex valued matrix, and a frequency response is defined. However, in contrast to a time invariant filter, the general linear filter may create frequency components which are not present in the input signal. This produces bias in the velocity and velocity spread estimates. It is shown that the clutter filter effect on the autocorrelation estimates can be described by a frequency domain transfer function, but unlike time invariant filters, the transfer function is different for each temporal lag of the autocorrelation function. Using a two dimensional (axial and temporal dimension) model of the received signal, the bias in velocity and velocity spread is quantified, both for the autocorrelation algorithm and the time shift cross-correlation estimator. Theoretical expressions, as well as numerical examples are given.     

Estimation of stochastic environment force for master–slave robotic system

The aim of this work is to obtain the maximum likelihood estimate (MLE) of controller-gain parameters \(\widehat{K} \) of the slave robot to determine the stochastic environment force. This is accomplished by measuring the joint positions of master and slave for a known master torque using stochastic difference equation. Here, the environmental force is modelled as a zero-mean white Gaussian random process. Therefore, the joint probability distribution function (pdf) of the slave angle over a given time duration can be computed as a function of the parameters ‘K’. This pdf is maximized with respect to ‘K’ to obtain the MLE of controller-gain parameters. Subsequently, convergence analysis of error in the estimates is performed. Also, an expression of the Cramer–Rao lower bound (CRLB) is derived to measure accuracy of the estimation. Comparison of CRLB with variance of MLE supports that our estimates are asymptotically efficient. The estimation performance is validated analytically and through simulations carried out on a two-link master–slave robotic system.     

A new wideband spread target maximum likelihood estimator for blood velocity estimation. I. Theory

The derivation and theoretical evaluation of new wideband maximum-likelihood strategies for the estimation of blood velocity using acoustic signals are presented. A model for the received signal from blood scatterers, using a train of short wideband pulses, is described. Evaluation of the autocorrelation of the signal based on this model shows that the magnitude, periodicity, and phase of the autocorrelation are affected by the mean scatterer velocity and the presence of a velocity spread target. New velocity estimators are then derived that exploit the effect of the scatterer velocity on both the signal delay and the shift in frequency. The wideband range spread estimator is derived using a statistical model of the target. Based on the point target assumption, a simpler wideband maximum-likelihood estimator is also obtained. These new estimation strategies are analyzed for their local and global performance. Evaluation of the Cramer-Rao bound shows that the bound on the estimator variance is reduced using these estimators, in comparison with narrowband strategies. In order to study global accuracy, the expected estimator output is evaluated, and it is determined that the width of the mainlobe is reduced. In addition, it is shown that the height of subsidiary velocity peaks is reduced through the use of these new estimators.     

A new time-domain narrowband velocity estimation technique for Doppler ultrasound flow imaging. II. Comparative performance assessment

For pt.I see ibid., vol.45, no.4, pp.939-54 (1998). The statistical performance of the new 2-D narrowband time-domain root-MUSIC blood velocity estimator described previously is evaluated using both simulated and flow phantom wideband (50% fractional bandwidth) ultrasonic data. Comparisons are made with the standard 1-D Kasai estimator and two other wideband strategies: the time domain correlator and the wideband point maximum likelihood estimator. A special case of the root-MUSIC, the "spatial" Kasai, is also considered. Simulation and flow phantom results indicate that the root-MUSIC blood velocity estimator displays a superior ability to reconstruct spatial blood velocity information under a wide range of operating conditions. The root-MUSIC mode velocity estimator can be extended to effectively remove the clutter component from the sample volume data. A bimodal velocity estimator is formed by processing the signal subspace spanned by the eigenvectors corresponding to the two largest eigenvalues of the Doppler correlation matrix. To test this scheme, in vivo common carotid flow complex Doppler data was obtained from a commercially available color flow imaging system. Velocity estimates were made using a reduced form of this data corresponding to higher frame rates. The extended root-MUSIC approach was found to produce superior results when compared to both 1- and 2-D Kasai-type estimators that used initialized clutter filters. The results obtained using simulated, flow phantom, and in vivo data suggest that increased sensitivity as well as effective clutter suppression can be achieved using the root-MUSIC technique, and this may be particularly important for wideband high frame rate imaging applications.     

An extended autocorrelation method for estimation of blood velocity

The conventional autocorrelation method (AM) to estimate the blood velocity for color flow imaging (CFI) is based on the phase estimation of the autocorrelation function. In this paper, a new extended autocorrelation method (EAM) that uses both phase and magnitude of the two dimensional (depth and temporal direction) autocorrelation function for estimating the blood velocity is presented. It is shown that the EAM has similar performance as the cross-correlation method (CCM). Both of them have smaller estimation variance than the AM and have the ability to estimate velocities beyond the Nyquist velocity, but the EAM is more computationally efficient than the CCM. A 2-D blood flow signal with rectilinear velocity including the transit time effect has also been simulated and the results are presented in this paper. For comparison, the EAM and the CCM have been applied to the simulated signals in which the flow velocities are up to four times the Nyquist velocity. The EAM has been further verified by experimental RF data from the subclavian artery     

Estimation of a dense velocity field based on the statistics of dynamic speckle

This article presents a new technique for flow velocity estimation from ultrasound image sequences. The method is based on the analysis of the temporal statistics of the speckle pattern in motion. We demonstrate that the biased local temporal variance (LTV) of a single pixel within an image of dynamic speckle is related to velocity. This allows us to estimate the total velocity magnitude without the requirement of neither block matching nor autocovariance estimation. A new estimator, asymptotically without bias, called LTV is presented. Results conducted on experimental B mode sequences (40 MHz) of blood mimicking fluid with calibrated velocities are presented. Performances of the estimator are studied and results show good agreement with the statistical model. Magnitude of the 3D velocity vector in the range of 0.1–2 mm/s have been estimated with a standard deviation error of less than 12%. The validity of the method and its limitations are then discussed. © 2010 Wiley Periodicals, Inc. Int J Imaging Syst Technol, 20, 268‐276, 2010     

Single-ensemble-based eigen-processing methods for color flow imaging--Part I. The Hankel-SVD filter

Because of their adaptability to the slow-time signal contents, eigen-based filters have shown potential in improving the flow detection performance of color flow images. This paper proposes a new eigen-based filter called the Hankel-SVD filter that is intended to process each slowtime ensemble individually. The new filter is derived using the notion of principal Hankel component analysis, and it achieves clutter suppression by retaining only the principal components whose order is greater than the clutter eigen-space dimension estimated from a frequency based analysis algorithm. To assess its efficacy, the Hankel-SVD filter was first applied to synthetic slow-time data (ensemble size: 10) simulated from two different sets of flow parameters that model: 1) arterial imaging (blood velocity: 0 to 38.5 cm/s, tissue motion: up to 2 mm/s, transmit frequency: 5 MHz, pulse repetition period: 0.4 ms) and 2) deep vessel imaging (blood velocity: 0 to 19.2 cm/s, tissue motion: up to 2 cm/s, transmit frequency: 2 MHz, pulse repetition period: 2.0 ms). In the simulation analysis, the post-filter clutter-to- blood signal ratio (CBR) was computed as a function of blood velocity. Results show that for the same effective stopband size (50 Hz), the Hankel-SVD filter has a narrower transition region in the post-filter CBR curve than that of another type of adaptive filter called the clutter-downmixing filter. The practical efficacy of the proposed filter was tested by application to in vivo color flow data obtained from the human carotid arteries (transmit frequency: 4 MHz, pulse repetition period: 0.333 ms, ensemble size: 10). The resulting power images show that the Hankel-SVD filter can better distinguish between blood and moving-tissue regions (about 9 dB separation in power) than the clutter-downmixing filter and a fixed-rank multi ensemble-based eigen-filter (which showed a 2 to 3 dB separation).     

Performance analysis of the maximum a posteriori probability direction finding

Abstract

In this paper, the mean square error performance of the maximum a posteriori (MAP) probability direction finding by sensor array in terms of its Cramer‐Rao lower bound (CRLB) is analyzed. Based on the principle of Bayesian estimator, a log posteriori probability function is formed when the a priori knowledge of location is given. The Fisher information matrix (FIM) is found accordingly. It shows that the CRLB of the MAP estimator is much lower than that of maximum likelihood technique, especially when the element SNR is low and/or the number of snapshots is small. In addition, the CRLB remains at a relatively low level in terms of the variance of DOA of sources. It also shows that the location variance dominates the behavior of the MAP direction finder when locations of sources are Gaussian distributed.     

Approximate Cramer-Rao bound on Doppler error in correlation-processing relatively narrowband noise radar

Limitations on accuracy of Doppler estimation in continuous-wave noise radar with correlation processing are studied. Second-order properties of output of the correlation receiver are evaluated and an approximate Cramer-Rao bound on errors of Doppler measurement is derived. The accuracy of Doppler measurements is found to be affected by the following factors: power spectral density of noise signal, frequency response of the low-pass filter in correlator, observation time, velocity of the target and signal to noise ratio. It is shown that the random nature of the transmitted signal induces additional fluctuations at the output of correlator which limit the accuracy even in the infinite signal to noise case. Qualitative extension of the results to a case covering multiple targets and clutter is made. It is argued that the performance will decrease and that increasing transmitted power may not provide significant improvement when clutter is present.     

A comparison of the performance of time-delay estimators in medical ultrasound

Time-delay estimation (TDE) is a common operation in ultrasound signal processing. In applications such as blood flow estimation, elastography, phase aberration correction, and many more, the quality of final results is heavily dependent upon the performance of the time-delay estimator implemented. In the past years, several algorithms have been developed and applied in medical ultrasound, sonar, radar, and other fields. In this paper we analyze the performances of the widely used normalized and non-normalized correlations, along with normalized covariance, sum absolute differences (SAD), sum squared differences (SSD), hybrid-sign correlation, polarity-coincidence correlation, and the Meyr-Spies method. These techniques have been applied to simulated ultrasound radio frequency (RF) data under a variety of conditions. We show how parameters, which include center frequency, fractional bandwidth, kernel window size, signal decorrelation, and signal-to-noise ratio (SNR) affect the quality of the delay estimate. Simulation results also are compared with a theoretical performance limit set by the Cramer-Rao lower bound (CRLB). Results show that, for high SNR, high signal correlation, and large kernel size, all of the algorithms closely match the theoretical bound, with relative performances that vary by as much as 20%. As conditions degrade, the performances of various algorithms differ more significantly. For signals with a correlation level of 0.98, SNR of 30 dB, center frequency of 5 MHz with a fractional bandwidth of 0.5, and kernel size of 2 /spl mu/s, the standard deviation of the jitter error is on the order of few nanoseconds. Normalized correlation, normalized covariance, and SSD have an approximately equal jitter error of 2.23 ns (the value predicted by the CRLB is 2.073 ns), whereas the polarity-coincidence correlation performs less well with a jitter error of 2.74 ns.     

Single-ensemble-based eigen-processing methods for color flow imaging--Part II. The matrix pencil estimator

Parametric spectral estimators can potentially be used to obtain flow estimates directly from raw slow-time ensembles whose clutter has not been suppressed. We present a new eigen-based parametric flow estimation method called the matrix pencil, whose principles are based on a matrix form under the same name. The presented method models the slow-time signal as a sum of dominant complex sinusoids in the slow-time ensemble, and it computes the principal Doppler frequencies by using a generalized eigen-value problem-formulation and matrix rank reduction principles. Both fixed rank (rank-one, rank-two) and adaptive-rank matrix pencil flow estimators are proposed, and their potential applicability to color flow signal processing is discussed. For the adaptive-rank estimator, the nominal rank was defined as the minimum eigen-structure rank that yields principal frequency estimates with a spread greater than a prescribed bandwidth. In our initial performance evaluation, the fixed-rank matrix pencil estimators were applied to raw color flow data (transmit frequency: 5 MHz; pulse repetition period: 0.175 ms; ensemble size: 14) acquired from a steady flow phantom (70 cm/s at centerline) that was surrounded by rigid-tissue-mimicking material. These fixed-rank estimators produced velocity maps that are well correlated with the theoretical flow profile (correlation coefficient: 0.964 to 0.975). To facilitate further evaluation, the matrix pencil estimators were applied to synthetic slow-time data (transmit frequency: 5 MHz; pulse repetition period: 1.0 ms; ensemble size: 10) modeling flow scenarios without and with tissue motion (up to 1 cm/s). The bias and root-mean-squared error of the estimators were computed as a function of blood-signal-to-noise ratio and blood velocity. The matrix pencil flow estimators showed that they are comparatively less biased than most of the existing frequency-based flow estimators like the lagone autocorrelator.     

Adaptive clutter filtering based on sparse component analysis in ultrasound color flow imaging

An adaptive method based on the sparse component analysis is proposed for stronger clutter filtering in ultrasound color flow imaging (CFI). In the present method, the focal underdetermined system solver (FOCUSS) algorithm is employed, and the iteration of the algorithm is based on weighted norm minimization of the dependent variable with the weights being a function of the preceding iterative solutions. By finding the localized energy solution vector representing strong clutter components, the FOCUSS algorithm first extracts the clutter from the original signal. However, the different initialization of the basis function matrix has an impact on the filtering performance of FOCUSS algorithms. Thus, 2 FOCUSS clutter- filtering methods, the original and the modified, are obtained by initializing the basis function matrix using a predetermined set of monotone sinusoids and using the discrete Karhunen-Loeve transform (DKLT) and spatial averaging, respectively. Validation of 2 FOCUSS filtering methods has been performed through experimental tests, in which they were compared with several conventional clutter filters using simplistic simulated and gathered clinical data. The results demonstrate that 2 FOCUSS filtering methods can follow signal varying adaptively and perform clutter filtering effectively. Moreover, the modified method may obtain the further improved filtering performance and retain more blood flow information in regions close to vessel walls.     

Clutter filters adapted to tissue motion in ultrasound color flow imaging

The quality of ultrasound color flow images is highly dependent on sufficient attenuation of the clutter signals originating from stationary and slowly moving tissue. Without sufficient clutter rejection, the detection of low velocity blood flow will be poor, and the velocity estimates will have a large bias. In some situations, e.g., when imaging the coronary arteries or when the operator moves the probe in search for small vessels, there is considerable movement of tissue. It has been suggested that clutter rejection can be improved by mixing down the signal with an estimate of the mean frequency prior to high pass filtering. In this paper, we compare this algorithm with several other adaptive clutter filtering algorithms using both experimental data and simulations. We found that realistic accelerations of the tissue have a large effect on the clutter rejection. The best results were obtained by mixing down the signal with non-constant phase increments estimated from the signal. This adapted the filter to a possibly accelerated tissue motion and produced a significant improvement in clutter rejection