Doppler angle estimation using correlation

Phase shift variance imaging - a new technique for destructive microbubble imaging

The detection of microbubble contrast agents with ultrasound imaging techniques is the subject of ongoing research. Commonly, the nonlinear response of the agent is employed for detection. The performance of these techniques is, however, affected by nonlinear sound propagation. As an alternative, the change in echo response resulting from microbubble destruction can be employed to detect the agent. In this work, we propose a novel criterion for microbubble destruction detection that allows the rejection of tissue at a defined significance level even for highly echogenic structures in the presence of nonlinear propagation. Most clinical systems provide the hardware requirements for acquisitions consisting of multiple pulses transmitted at the same position, as used in Doppler imaging. Therefore, we develop a processing strategy that distinguishes contrast agent from other stationary or moving structures using these sequences. The proposed criterion is based on the variance of the phase shift of consecutive echoes in the sequence, which, in addition to tissue rejection, permits the distinction of motion from agent disruption. Phantom experiments are conducted to show the validity of the criterion and demonstrate the performance of the new method for contrast detection. Each detection series consists of 20 identical pulses at 9.5 MHz (4.7 MPa peak negative pressure) transmitted at a pulse repetition frequency of 5 kHz. The sequence is applied to phantoms under varied motion and flow conditions. As a first step toward molecular imaging, the technique is applied to microbubbles targeted to vascular endothelial growth factor receptor 2 (VEGFR2) in vitro. The results show a uniform rejection of the background signal while maintaining a contrast enhancement by more than 40 dB. The area under the receiver operating characteristics (ROC) curve is used as the performance metric for the separation of contrast agent and tissue signals, and values larger than 97% demonstrate that an excellent separation was achieved.

Volumetric Flow Measurement Using an Implantable CMUT Array

This paper describes volumetric-flow velocity measurement using an implantable capacitive micromachined ultrasonic transducer (CMUT) array. The array is comprised of multiple-concentric CMUT rings for ultrasound transmission and an outmost annular CMUT array for ultrasound reception. Microelectromechanical-system (MEMS) fabrication technology allows reception CMUT on this flowmeter to be implemented with a different membrane thickness and gap height than that of transmission CMUTs, optimizing the performance of these two different kinds of devices. The silicon substrate of this 2-mm-diameter CMUT ring array was bulk micromachined to approximately 80 to 100 μm thick, minimizing tissue disruption. The blood-flow velocity was detected using pulse ultrasound Doppler by comparing the demodulated echo ultrasound with the incident ultrasound. The demodulated ultrasound signal was sampled by a pulse delayed in time domain from the transmitted burst, which corresponds to detecting the signal at a specific distance. The flow tube/vessel diameter was detected through the time-flight delay difference from near and far wall reflections, which was measured from the ultrasound pulse echo. The angle between the ultrasound beam and the flow was found by using the cross-correlation from consecutive ultrasound echoes. Artificial blood flowing through three different polymer tubes was experimented with, while keeping the same volumetric flow rate. The discrepancy in flow measurement results between this CMUT meter and a calibrated laser Doppler flowmeter is less than 5%.

Aperture-domain processing and its applications in ultrasound imaging: a review

This paper provides a review of advanced algorithms for ultrasound image formation and signal processing that are based on aperture-domain data (i.e. the data recorded by individual channels prior to beam summation). First aperture-domain data are defined and their properties described, then two specific examples of phase-aberration correction and vector velocity estimation are presented. For phase-aberration correction, sidelobe-reduction techniques based on the coherence of the received aperture-domain data were tested with clinical breast data; the mean improvements in the contrast and contrast-to-noise ratios were 6.9 dB and 23.2 per cent, respectively. For flow estimation, a conventional scanner can only estimate the flow velocity parallel to the beam axis. The proposed flow estimation technique uses aperture-domain data for two-dimensional flow-velocity estimation. The experimental results demonstrate that the estimation errors for the proposed technique are 2.18 per cent and 18.11 per cent in the axial and lateral velocity components, respectively. Other applications in which aperture-domain data can be used are also discussed.

Photoacoustic flow measurements based on wash-in analysis of gold nanorods

In this study, photoacoustic flow measurement methods based on wash-in analysis are presented. These methods use the rod-to-sphere shape transformations of gold nanorods induced by pulsed-laser irradiation. Due to the shape dependence of the optical absorption of the gold nanorods, these shape transitions are associated with a change in the peak optical absorption wavelength. Pulsed-laser irradiation at the wavelength corresponding to the peak optical absorption of the original gold nanorods allows the particles that undergo shape changes to be viewed as "being destructed" by the laser irradiation at that wavelength, hence, flow information can be derived from the change in ultrasound intensity that is directly related to the wash-in rate of the gold nanorods and the laser intensity. Two flow estimation methods based on the wash-in analysis are described. The first method first applies high-energy laser pulses that induce shape changes in all the nanorods. A series of low-energy pulses then are applied to monitor the acoustic signal change as new nanorods flow into the region of interest. The second method uses single-energy laser pulses such that the "destruction" and "detection" are performed simultaneously. The simulation results show that it is valid to fit the time-intensity curves by exponential models. To demonstrate the validity of the proposed methods, an Nd:YAG pulsed laser operating at 1064 nm was used for optical irradiation, and a 1-MHz ultrasonic transducer was used for acoustic detection. Gold nanorods with a peak optical absorption at 1018 nm and a concentration of 0.26 nM were used to estimate flow velocities ranging from 0.35 to 2.83 mm/s. The linear regression results show that the correlation coefficients between the measured velocities and the true values are close to unity (> or = 0.94), thus demonstrating the feasibility of the proposed photoacoustic techniques for relative flow estimation.

Improved transverse flow estimation using differential maximum Doppler frequency

Conventional Doppler technique can only provide the axial component of the blood flow vector, which is actually a three dimensional (3-D) quantity. To acquire the complete flow vector, estimations of the other two velocity components are essential. For the two dimensional (2-D) Doppler-bandwidth-based transverse estimation, however, accuracy is generally limited because of the complex dependence of the Doppler spectral shape on the flow variation within the sample volume. Two factors that may lead to the Doppler spectral change were considered in this study. One is the position offset of the sample volume and the other is the length of the sample volume. Simulations were performed and experimental data were also collected. Results indicate that the position offset may result in severe underestimation of Doppler shift frequency. Consequently, Doppler bandwidth is overestimated when it is determined by the difference between Doppler shift frequency and maximum Doppler frequency. Compared with the position offset, influence of the length of sample volume on the Doppler bandwidth is minor. To overcome this problem, a novel method, which is based on the differential maximum Doppler frequency, is proposed. Specifically, two beams with different beam widths are simultaneously generated to observe the blood flow and the difference between the corresponding maximum Doppler frequencies is used to estimate the transverse velocity. It is demonstrated that the accuracy and stability of transverse estimation are significantly improved by the proposed method even when the position offset is present.

Estimating the blood velocity vector using aperture domain data

Most conventional blood flow estimation methods measure only the axial component of the blood velocity vector. In this study, we developed a new method for two-dimensional (2-D) velocity vector estimation in which time shifts resulting from blood motion are calculated for the individual channels using aperture domain data. This allows the construction of a time-shift profile along the array direction as a function of channel index, which is approximated by a first-order polynomial whose zeroth-order and first-order terms can be used to determine the axial and lateral velocity components, respectively. The efficacy of the proposed method was verified by simulations and experiments in which the transducer array had 64 elements, an aperture size of 1.96 cm, and a center frequency of 5 MHz. The flow velocity ranged from 5 to 35 cm/s and the Doppler angle ranged from 0 degrees to 90 degrees. The experimental results show that the accuracy of axial velocity estimation is higher for the new method than for the autocorrelation-based conventional method when the signal-to-noise ratio is larger than 0 dB. The mean estimation error for the axial velocity component is 2.18% for the new method, compared to 4.51% for the conventional method. The mean estimation error for the lateral velocity component is 15%, which is comparable to existing methods.

Vector-velocity estimation in swept-scan using a K-space approach

The swept-scan technique (i.e., continuously moving a single-crystal transducer during pulse-echo data acquisition) is used in high-frequency, ultrasonic flow imaging. Relative to the conventional step-scan technique, swept scanning improves the rate of data acquisition and enables near-real-time, high-frequency color flow mapping. However, the continuous transducer movement may have non-negligible effects on accuracy of velocity estimation. This paper introduces a spatial frequency domain (i.e., k-space) approach that quantifies the effects of both lateral and axial motions in a swept scan. It is shown that the k-space representation is equivalent to a Doppler-radio frequency (RF) frequency domain representation, and that transducer movement in the swept-scan technique results in a change in Doppler bandwidth. In addition, a vector velocity estimator is developed based on the proposed k-space approach. Both simulations and flow-phantom experiments were performed to evaluate the performance of the proposed vector velocity estimator. A 45-MHz transducer was scanned at 20 mm/s. The Doppler angle ranged from 29 degrees to 90 degrees, and the flow velocities ranged from 15 to 30 mm/s. The results show that the proposed k-space vector velocity estimator exhibited a mean error of 2.6 degrees for flow-direction estimation, with the standard deviation ranging from 2.2 degrees to 8.2 degrees. In comparison, for the conventional spectral-broadening-based vector velocity estimator ignoring the swept-scan effect, the mean error became 15 degrees and the standard deviations were from 2.7 degrees to 6.6 degrees.

Effects of swept scanning on velocity estimation

The swept-scan technique adopted in high-frequency ultrasound involves mechanically scanning a single-element transducer to acquire image data. Unlike conventional step scanning, where the image data are acquired at discrete positions, the swept-scan technique acquires the image data while the transducer is continuously moving. Such a scanning method is particularly advantageous for Doppler flow estimation because its frame rate is higher than that for the step-scan technique. However, the effects of the transducer motion on the accuracy of velocity estimation have not been studied comprehensively. This study employed a k-space approach to experimentally investigate the effects of swept scanning on both conventional Doppler axial velocity estimation and spectral-broadening-based lateral velocity estimation using a 45-MHz transducer. The results indicate that such effects must be corrected in order to obtain an accurate estimation of flow velocities.

Quantitative analysis of sinoatrial node using Doppler tissue images

A method has been developed to quantitatively analyze sinoatrial nodes (SAN) using Doppler tissue images (DTI). Doppler tissue images of SAN are acquired using an intracardiac catheter via the superior vena cava in an in vivo experiment. A sequence of DTI images of a SAN is obtained, and a complete cycle of the SAN excitation is observed. The tissue acceleration of the SAN is extracted and quantitatively analyzed. The estimated time-acceleration curve of the SAN exhibits remarkable similarity to the electrocardiogram curve. This is the first report on such finding. The experimental results show that the tissue movement of the SAN correlates with electrical cardiac activities and closely associates with the different phases of the cardiac cycle. This method has great potential in characterizing the local cardiac activities through the study of the conduct pathway.

A new method of ultrasound color flow mapping

Conventional ultrasound color flow mapping systems estimate and visualize only the axial velocity component. To obtain the transverse velocity component a modification of a multiple-beam method is proposed. The new two-dimensional color flow mapping system has a small size and consists of three transducers. The central transducer is an appodized and focused phased array. The other transducers are unfocused probes. Three transducers act as receivers and the central transducer operates as a transmitter. All receivers acquire rf scan lines that are then processed to estimate three axial velocity components using an autocorrelation method. These estimates are then combined to estimate the transverse velocity component, taking into account the geometric relationships among three transducers. Two algorithms for transverse velocity estimation are proposed. The first uses the Doppler angle estimate for calculation of the transverse velocity component. The other algorithm calculates the transverse velocity component directly from the axial components. The accuracy of the flow velocity estimators is estimated by simulations. Analysis of accuracy allows choosing the more effective algorithm for two-dimensional velocity estimation, which is insensitive to variations of the Doppler angle.

Doppler angle estimation using AR modeling

The transit time spectrum broadening effect has long been explored for Doppler angle estimation. Given acoustic beam geometry, the Doppler angle can be derived based on the mean Doppler frequency and the Doppler bandwidth. Spectral estimators based on the fast Fourier transform (FFT) are typically used. One problem with this approach is that a long data acquisition time is required to achieve adequate spectral resolution, with typically 32-128 flow samples being needed. This makes the method unsuitable for real-time two-dimensional Doppler imaging. This paper proposes using an autoregressive (AR) model to obtain the Doppler spectrum using a small number (e.g., eight) of flow samples. The flow samples are properly selected, then extrapolated to ensure adequate spectral resolution. Because only a small number of samples are used, the data acquisition time is significantly reduced and real-time, two-dimensional Doppler angle estimation becomes feasible. The approach was evaluated using both simulated and experimental data. Flows with various degrees of velocity gradient were simulated, with the Doppler angle ranging from 20 degrees to 75 degrees. The results indicate that the AR method generally provided accurate Doppler bandwidth estimates. In addition, the AR method outperformed the FFT method at smaller Doppler angles. The experimental data for Doppler angles, ranging from 33 degrees to 72 degrees, showed that the AR method using only eight flow samples had an average estimation error of 3.6 degrees, which compares favorably to the average error of 4.7 degrees for the FFT method using 64 flow samples. Because accurate estimates can be obtained using a small number of flow samples, it is concluded that real-time, two-dimensional estimation of the Doppler angle over a wide range of angles is possible using the AR method.

Doppler angle estimation of pulsatile flows using AR modeling

In quantitative ultrasonic flow measurements, the beam-to-flow angle (i.e., Doppler angle) is an important parameter. An autoregressive (AR) spectral analysis technique in combination with the Doppler spectrum broadening effect was previously proposed to estimate the Doppler angle. Since only a limited number of flow samples are used, real-time two-dimensional Doppler angle estimation is possible. The method was validated for laminar flows with constant velocities. In clinical applications, the flow pulsation needs to be considered. For pulsatile flows, the flow velocity is time-varying and the accuracy of Doppler angle estimation may be affected. In this paper, the AR method using only a limited number of flow samples was applied to Doppler angle estimation of pulsatile flows. The flow samples were properly selected to derive the AR coefficients and then more samples were extrapolated based on the AR model. The proposed method was verified by both simulations and in vitro experiments. A wide range of Doppler angles (from 3o degrees to 78 degrees) and different flow rates were considered. The experimental data for the Doppler angle showed that the AR method using eight flow samples had an average estimation error of 3.50 degrees compared to an average error of 7.08 degrees for the Fast Fourier Transform (FFT) method using 64 flow samples. Results indicated that the AR method not only provided accurate Doppler angle estimates, but also outperformed the conventional FFT method in pulsatile flows. This is because the short data acquisition time is less affected by the temporal velocity changes. It is concluded that real-time two-dimensional estimation of the Doppler angle is possible using the AR method in the presence of pulsatile flows. In addition, Doppler angle estimation with turbulent flows is also discussed. Results show that both the AR and FFT methods are not adequate due to the spectral broadening effects from the turbulence.

Time-intensity-based volumetric flow measurements: an in vitro study

Ultrasonic contrast agents have been used to assist blood flow measurements. Several contrast-specific flow measurement techniques have been proposed during the last few years. Among them, a method based on relative enhancement of the backscattered signal as a function of time is of particular interest. This method is also known as the time-intensity method. The method is based on the indicator-dilution theory, and the time-intensity curve is used to derive blood flow-related parameters such as the flow rate and the blood mixing volume. Previous in vitro studies done by other research groups were mainly based on a perfusion model or an artery model. Results showed that several parameters derived from the time-intensity curve had a good correlation with the flow rate under certain conditions. However, the studies did not focus on factors such as mixing volume, mixing chamber configuration and different types of mixing chamber. In this paper, dependence of the time-intensity curve is further studied. Specifically, two types of blood-mixing chambers were constructed. One was a spherical compartment phantom with two different sizes (260 and 580 mL) and different inflow/outflow configurations. The other was a perfusion phantom consisting of dialysis cartridges with the volume ranging from 114 to 351 mL. The time intensities were also measured at both the input and the output of the mixing chamber. A commercial agent (Levovist) and a self-made, albumin-based agent were used and the wash-out time constant and the mean transit time were derived for flow rates ranging from 500 to 1300 mL/min. For the perfusion phantom, results showed that the parameters had a good correlation with both the flow rate and the mixing volume. Results from the compartment phantom, on the other hand, indicated that the inflow/outflow configuration and the mixing size significantly affected the derived time constants. Potential applications of new volumetric flow estimation techniques based on both input and output intensities were also discussed.