Subsample volume processing of doppler ultrasound signals

Processing of Doppler signals produced by pulsed Doppler systems is based on the assumption that the phase of the received high frequency ultrasound signals changes linearly with depth. However, the random spatial distribution of scatterers is not in accordance with this basic assumption. Consequently, averaging of the demodulated signal over an observation window, covering a few periods of the received signal, does not improve the estimate for the instantaneous quadrature components of the Doppler signal originating from a given depth. Hence, the accuracy of the Doppler velocity estimate is independent of the length of the observation window employed. However, splitting the observation window in subsample volumes, each with a length of one period at the emission frequency, and combining the Doppler signals of the subsample volumes at the last stage of signal processing, i.e., mean Doppler frequency estimation using the autocorrelation technique, results in a considerable reduction of the variance of the velocity estimate. Using a computer simulation of the signal processing involved, it is demonstrated that with subsample volume processing the variance of the velocity estimate attains the same variance as is expected for the RF cross correlation technique.     

Estimation of Doppler shift frequency using selected phase information for high frame rate color flow mapping

Purpose We describe a new approach to processing signals used to estimate the Doppler shift frequency in high frame-rate color flow mapping with fewer pulse transmissions. When an ultrasound pulse is transmitted to a large number of scatterers, the echoes from the scatterers overlap and interfere with one another. This interference causes the phase of the received echo signal to fluctuate, thus disturbing the estimated shift in Doppler frequency. The technique proposed here eliminates this disturbed phase information, leaving the remaining information for use in estimating the shift in Doppler frequency. The instantaneous frequency of the echo signal can serve as an index of the influence of interference.Methods To test this technique in vivo we used radio-frequency echo signals from the carotid artery for simulation and evaluated the error of the estimated Doppler shift frequency in several cases.Conclusion Performance was enhanced when the number of pulses transmitted was limited and this technique was used.This article is a translation of the original that was published in Jpn J Med Ultrasonics 2001;28:J15–23     

A simulation of transit time effects in Doppler ultrasound signals

A signal model is proposed which can be used to study frequency extraction techniques for Doppler ultrasound. The signal is based on the physics of the Doppler process and depends on a sliding window used to average a set of independent Gaussian random numbers. This window is related to the shape of the sample volume for the Doppler pulse and depends on the Doppler angle. Simulation results compare favorably with results from flow experiments in terms of the variance of the estimated Doppler shift, the shape of the power spectra and the behavior of the signals with respect to Burg autoregressive power spectra. A potential use of the signal in the study of spectral analysis techniques is presented.     

Comparison of the performance of the RF cross correlation and doppler autocorrelation technique to estimate the mean velocity of simulated ultrasound signals

In pulsed Doppler systems the received RF (radio frequency) signal is multiplied by a quadrature reference signal and subsequently averaged over a short depth range to obtain a sample of the complex Doppler signal. The mean frequency of the sampled Doppler signal, obtained with the autocorrelation function, reflects the mean velocity of the scatterers moving through the sample volume. An alternative is to evaluate the two-dimensional cross correlation function of a short segment of the RF signals over subsequent lines, giving the mean velocity of the scatterers. Both methods of velocity estimation were applied to computer-generated RF signals with varying RF bandwidth, signal-to-noise ratio, and mean and width of the imposed velocity distribution. The length of the RF signal segment and the number of lines for velocity estimation (package length) affects the accuracy of the velocity estimate. It can be concluded that the cross correlation technique behaves superiorly especially for a low velocity dispersion. Furthermore, the standard deviation of the velocity estimate decreases for an increasing sample volume length and package length, while the performance of the conventional Doppler technique is rather independent of the length of the sample volume. The difference between both techniques decreases for a greater package length or for signals simulating a wide velocity distribution.     

Effects of transducer beam geometry and flow velocity profile on the Doppler power spectrum: a theoretical study

A theoretical model is used to show how the Doppler spectrum for various axisymmetric velocity profiles is affected by beam misalignment and incomplete insonation. Results are presented for both circular and square beam geometries. Moreover, a closed-form expression is derived for the power spectral density received by an on-axis transducer with a Gaussian beam profile. It is shown that the error incurred in measuring the mean Doppler frequency with such a profile will generally be bounded by the results for the circular and square beam geometries. The effects of an ideal high-pass filter on the mean Doppler frequency and the backscattered Doppler power are examined. It is shown that such a filter can introduce large differences in the measured systolic to diastolic power ratios. Finally, theoretical expressions and results are presented for the spectral broadening index (SBI), normalized spectral variance (NSV), coefficient of kurtosis (CK), the coefficient of skewness (CS) as functions of the axisymmetric velocity profile shape assuming complete uniform insonation.     

Influence of spectral broadening on continuous wave Doppler ultrasound spectra: a geometric approach

A model is presented that enables the detailed effects of spectral broadening to be calculated for a continuous wave (CW) Doppler system by using geometric boundary arguments. The model assumes a uniform distribution of isotropic scatterers and treats the transmitter and receiver crystals as incremental sources and receivers. Detailed results for rectangular and circular geometries are presented in order to provide a physical understanding of the manner in which spectral broadening arises. Results are given for the circular geometry, to illustrate the manner in which the received spectrum is affected by the transducer size and distance from the vessel.     

Subspace-Based Blood Power Spectral Capon Combined with Wiener Postfilter to Provide a High-Quality Velocity Waveform with Low Mathematical Complexity

In Doppler analysis, the power spectral density (PSD), which accounts for the axial velocity distribution of the blood scatterers, is estimated. The conventional spectral estimator is Welch's method, which suffers from frequency leakage at small observation window length. The performance of adaptive techniques such as blood power Capon (BPC) has been promising at the cost of higher computation complexity. Reducing the computational complexity while retaining the benefits of BPC would be necessary for real-time implementation. The purpose of the work described here was to investigate whether it is possible to decrease the computation load in BPC and still obtain acceptable results. The computation complexity in BPC is owing primarily to the matrix inversion required for computing the PSD estimate. We here propose the subspace blood power Capon technique, which employs a data covariance matrix with reduced number of rows in estimation of the weight vector. In maximum velocity estimation in the spectra, the signal noise slope intersection envelop estimator that makes use of the integrated power spectrum is employed. The evaluations are made based on both simulated and in vivo data. The results indicate that it is possible to reduce the order of complexity to almost 12.25% at the cost of 2.31% and 2.24% increases in the relative standard deviation and relative bias of the estimates. Moreover, the Wiener post-filter as a post-weighting factor, which will be multiplied by the final weight vector of the spectral estimator, estimates the power of the desired signal and the power of the interference plus noise to improve the contrast. The proposed estimator has exhibited a promising performance at beam-to-flow angles of 45°, 60° and 75°. Furthermore, the robust performance of the proposed estimator against variation in the flow rate is also documented.     

Selective transmission of a focused Doppler ultrasound beam through a plastic layer

Laboratory test objects are widely used in Doppler ultrasound (US). Although the acoustic properties of in vitro materials are usually known, they are unlikely to match each other, or their in vivo counterparts, exactly. We conducted theoretical and experimental studies of a focused ultrasound beam as it passes from one fluid, through an intervening plastic layer at an oblique angle, and then into a different fluid. Dual mode propagation may occur (i.e., both longitudinal and shear waves can propagate in the plastic layer). Our calculations show that the power transmitted by either mode drops very rapidly to zero at certain critical angles. A range of angles of incidence exists within a focused beam and this, combined with the highly angle-dependent power transmission behaviour, can produce major distortions of Doppler data. These may persist even when the beam axis is not oriented exactly at the critical angle. The total power transmitted depends on all the wave speeds, may involve mode conversion, and is a very complicated function of the angle of incidence. This study reports a practical method for the calculation of power transmission though a plastic layer, and shows how the resulting power vs. angle graph can be used to avoid artefacts in in vitro Doppler studies.     

Differences in definity and optison microbubble destruction rates at a similar mechanical index with different real-time perfusion systems.

The purpose of this study was to determine microbubble responses to different pulse sequence schemes that exist on low mechanical index (MI) real-time perfusion imaging systems using either intravenous albumin-coated (Optison) or lipid-encapsulated (Definity) microbubbles. A tissue-mimicking phantom was created that permitted insonation of microbubbles at 3 cm (near field) and 9 cm (far field) from the diagnostic transducer face. Differences in effluent microbubble concentration were measured after they passed through vessels being insonified with pulse sequence schemes that transmitted alternating polarity (pulse inversion Doppler), alternating amplitude (power modulation), or both (contrast pulse sequencing) at a similar MI, frame rate, and transmit frequency. Normalized contrast signal intensity within a recirculating chamber was also measured in the near and far field. Pulse inversion Doppler produced less initial normalized contrast signal intensity and greater destruction rates than amplitude varying pulse sequence schemes like power modulation or contrast pulse sequencing at both the 0.1- and 0.2-MI settings. These differences indicate that the same MI setting on different real-time perfusion imaging techniques will produce different microbubble responses.     

Coded excitation in TCD ultrasound systems to improve axial resolution

The radiofrequency (RF) signal from a transcranial Doppler (TCD) ultrasound system allows us to track the motion of an embolus in a 2-D space of time and depth. The technique is limited by the narrow bandwidth (long duration) of the transmitted pulse because this provides poor axial resolution. This study aimed to assess whether implementing coded excitation and pulse compression in a TCD system would be a feasible means of increasing the bandwidth and hence improving the axial resolution, without the need for reducing the pulse length and increasing peak power levels. Embolic signals were collected in vitro from a flow phantom using a TCD system, which alternately transmitted coded and noncoded pulses. This allowed the same event to be investigated using the two different types of processing. Quantitative and qualitative measures were used to compare the two processing methods. The in vitro results were promising. They showed that the axial resolution could be improved, on average, by a factor of 6.6, using a pulse length of 13 micros and a chirp bandwidth of 0.8 MHz. This was reduced to a factor of 6 when a temporal bone sample was placed between the transducer and the flow phantom. Qualitatively the journey of an embolus was easier to track using the pulse compressed signal. Sonograms could be generated from smaller receive gates using the pulse compressed signal while still achieving an adequate signal-to-noise ratio. We found that it is both feasible and beneficial to implement coded excitation and pulse compression in a TCD system.     

Implementation of spectral width Doppler in pulsatile flow measurements

In this paper, we present an automatic beam-vector (Doppler) angle and flow velocity measurement method and implement it in pulsatile flow measurements using a clinical Doppler ultrasound system. In current clinical Doppler ultrasound flow velocity measurements, the axis of the blood vessel needs to be set manually on the B-scan image to enable the estimation of the beam-vector angle and the beam-vector angle corrected flow velocity (the actual flow velocity). In this study, an annular array transducer was used to generate a conical-shaped and symmetrically focused ultrasound beam to measure the flow velocity vectors parallel and perpendicular to the ultrasound beam axis. The beam-vector angle and flow velocity is calculated from the mode frequency (f(d)) and the maximum Doppler frequency (f(max)) of the Doppler spectrum. We develop a spectrum normalization algorithm to enable the Doppler spectrum averaging using the spectra obtained within a single cardiac cycle. The Doppler spectrum averaging process reduces the noise level in the Doppler spectrum and also enables the calculation of the beam-vector angle and flow velocity for pulsatile flows to be measured. We have verified the measurement method in vivo over a wide range of angles, from 52 degrees to 80 degrees, and the standard deviations of the measured beam-vector angles and flow velocities in the carotid artery are lower than 2.2 degrees and 12 cm/s (about 13.3%), respectively.     

Measurement of blood perfusion in tissue using Doppler ultrasound

A diagnostic tool for noninvasive evaluation of microcirculatory blood flow using continuous-wave CW Doppler ultrasound is presented. In this study, the properties of this method are investigated both theoretically and experimentally. The method utilizes a nondirectional CW Doppler flowmeter. Blood perfusion in tissue is shown to be proportional to the integral integral of fS (f)df where S(f) is the Doppler power spectrum and f is the Doppler frequency. The instrumentation needed to implement the method is described. Using an experimental flow model it is demonstrated that the above integral is proportional to the product between the number of scatterers in the sample volume of the Doppler probe and the mean speed of these scatterers. This is true even for low flow velocities (down to 1 mm/s). The results from in-vivo measurements on tissues in the finger, and the calf demonstrate that the method can monitor changes in the blood perfusion. It also shows the present limitations of the method due to movement artefacts.     

Doppler ultrasound tracking instrument for monitoring blood flow velocity

Doppler ultrasound (US) is potentially a valuable method for monitoring changes of blood flow velocity over a period of many minutes or even hours, but is seldom used in this way. One difficulty that may have contributed to this is the problem of maintaining a fixed geometry between the US beam and the blood vessel. A method of improving the success of monitoring might be to actively steer the US beam so as to maintain an adequate signal even when small displacements of the transducer occur. We have designed and built a prototype system for this purpose. The system comprises a continuous-wave phased-array transducer controlled by a purpose-built Doppler unit. The system constantly evaluates the quality of the returning Doppler signal in terms of total power and signal-to-noise ratio (SNR) (evaluated by assessing the quality of derived envelope signals), and steers the ultrasonic beam in a manner so as to improve the signal, should this be necessary. The system was tested in vitro, where the automatic tracking of the Doppler signal doubled the effective beam width of the transducer. Further developments that increase sensitivity and steering range should result in US Doppler systems that are better suited to long-term monitoring.     

Auditory phenomena during transcranial Doppler insonation of the basilar artery

During transcranial Doppler insonation of the basilar artery, 15 normal subjects heard high-frequency tones similar to tinnitus. The pitch of the tones equated with the PRF, and intensity equated with loudness. Alteration of pulse length did not change the tones heard, and the effect was highly dependent on probe position and direction. This is an example of ultrasonography acting as a sensory stimulus. The exact mechanism and implications of this effect will need to be clarified by further investigation. Although this phenomenon is not proved to be harmful, it would seem prudent to use the minimum transmitted power necessary to obtain a satisfactory signal from the basilar artery.     

An improved ultrasound simulation model: use in evaluating log versus linear processing for lesion detection

This paper reports the development of an improved three-dimensional computer simulation model for evaluation of ultrasonic imaging systems. This model was used to successfully evaluate a signal processing method for improving lesion detection in ultrasound imaging. Linear processing of the rf signal amplitudes from a limited region of tissue was compared with the logarithmic compression employed by most commercial scanners. Two lesions were simulated by spherical distributions of scatterers having backscatter coefficients greater than the scatterers in the surrounding medium. Linear processing improved the differential contrast by a factor of about two. The simulation is based on the three-dimensional distribution of acoustic frequency spectra in a transducer beam and integration of scattered pulses from a corresponding three-dimensional array of scatterers. The simulation reported in previous papers depended upon physical measurement of the impulse response of a transducer. An original contribution described briefly herein, and in more detail in a companion article, is the addition of a model of the transducer's pulse waveform generation. Another new addition is the definition of a specific lesion detection task for objective assessment of a change in image quality following perturbation of some system parameter.     

彩色Doppler能量图血流定量的实验研究

本研究的目的在于建立应用Ｄｏｐｐｌｅｒ能量图进行血流量测定的方法。应用标准预设的超声Ｄｏｐｐｌｅｒ能量图测定体外模拟循环系统管道内流体的Ｄｏｐｐｌｅｒ散射能量，观察一定液体浓度及Ｄｏｐｐｌｅｒ增益水平时管道横截面上Ｄｏｐｐｌｅｒ能量信号变化。利用计算机辅助图像处理求得各水平的管道截面上平均能量强度（彩色亮度），并与实际流量对比。结果显示在一定的增益浓度水平，管道截面上平均能量强度与实际流量之间存在直线相关关系（Ｐ值＜０．０１）。本研究证明在适当的图像设置条件下Ｄｏｐｐｌｅｒ能量图能较准确地反映血流量变化的信息。此方法在组织灌流、肿瘤血管定量分析等领域有潜在的应用价值。     

Effects of atherosclerotic lesions on ultrasonic beam and cw doppler signals

Doppler ultrasound has been shown to be a useful clinical tool in the noninvasive detection of vascular stenoses, occlusion, or plaques. The presence of atherosclerotic lesions in blood vessel walls can distort the Doppler ultrasonic signal, producing inverted or biphasic tracings. This observation, in fact, has been used frequently as a criterion for the diagnosis of plaque formation and calcification of these lesions. To define the source of this distortion, we have examined the effects of various types of atherosclerotic lesions on the transmitted ultrasonic beam generated by a continuous wave (CW) Doppler ultrasonic probe. These results are compared to the Doppler flow tracings obtained by a CW ultrasonic Doppler flowmeter from blood vessels with overlaying lesions in a mock flow system. It is concluded that severe distortion in the form of beam refraction or diffraction can be produced by calcified lesions, resulting in biphasic or inverted flow signals.     

Dependence of power Doppler image on a high pass filter instrumented in ultrasound machine

The aim of this study was to examine the influence of flow velocity, beam incident angle and a wall motion filter on the intensity of power Doppler image of steady flow in vitro. Power Doppler images of flow were recorded using a 7.5-MHz linear transducer with the fixed repetition frequency, gain and frame rate settings at the velocity levels of 15, 20, 35 and 65 cm/s with the beam incident angles of 30° and 60°. The power (P) of received Doppler signal was digitized to 25 gray values (G) for imaging using a wall motion filter with the transfer function: G = a × log P + b (a, b = constant values). Image intensity was measured at the center of the flow image off-line and compared with the flow velocity measured by conventional pulsed wave Doppler. Below the Nyquist limit, G increased as velocity increased, and G with the incident angle 30° was lower than those with 60°. Overall relationship between the Doppler shift frequency and the image intensity exhibited an excellent correlation when fit to the theoretical curve derived from the filter property (r = 0.97). The signal intensity of power Doppler image of flow depends upon the mean Doppler frequency shift or the flow velocity and was clearly affected by the wall motion filter.     

超声多普勒功率及其在检测颅内血流分布上的应用

本文结合背向散射超声多普勒信号模型对超声多普勒功率的概念及其特性进行了阐述，并提出了一种用超声多普勒功率来检测颅内血管血流分布的方法。实验结果证明，该方法能在噪声背景下灵敏准确地反应沿声束指向一维空间各个位置是否存在血流以及血流强度的大小。这一方法已被用于经颅多普勒（ＴＣＤ）脑血流分析仪中，解决了现有经颅多普勒检查时只能凭经验盲目搜寻颅内血管的问题，并提供了颅内血管分布的相对位置信息以帮助临床医生更方便准确地判定所检测的是哪一根血管。     

Deviations from Rayleigh statistics in ultrasonic speckle

The statistics of speckle patterns in ultrasound images have potential for tissue characterization. In "fully developed speckle" from many random scatterers, the amplitude is widely recognized as possessing a Rayleigh distribution. This study examines how scattering populations and signal processing can produce non-Rayleigh distributions. The first order speckle statistics are shown to depend on random scatterer density and the amplitude and spacing of added periodic scatterers. Envelope detection, amplifier compression, and signal bandwidth are also shown to cause distinct changes in the signal distribution.