This protocol shows how to apply ultrafast ultrasound Doppler imaging to quantify blood flows. After a 1 s long acquisition, the experimenter has access to a movie of the full field of view with axial velocity values for each pixel every ≈0.3 ms (depending on the ultrasound time of flight).
The pulsed-Doppler effect is the main technique used in clinical echography to assess blood flow. Applied with conventional focused ultrasound Doppler modes, it has several limits. Firstly, a finely tuned signal filtering operation is needed to distinguish blood flows from surrounding moving tissues. Secondly, the operator must choose between localizing the blood flows or quantifying them. In the last two decades, ultrasound imaging has undergone a paradigm shift with the emergence of ultrafast ultrasound using unfocused waves. In addition to a hundredfold increase in framerate (up to 10000 Hz), this new technique also breaks the conventional quantification/localization trade-off, offering a complete blood flow mapping of the field of view and a simultaneous access to fine velocities measurements at the single-pixel level (down to 50 µm). This data continuity in both spatial and temporal dimensions strongly improves the tissue/blood filtering process, which results in an increase sensitivity to small blood flow velocities (down to 1 mm/s). In this method paper, we aim to introduce the concept of ultrafast Doppler as well as its main parameters. Firstly, we summarize the physical principles of unfocused wave imaging. Then, we present the Doppler signal processing main steps. Particularly, we explain the practical implementation of the critical tissue/blood flow separation algorithms and on the extraction of velocities from these filtered data. This theoretical description is supplemented by in vitro experiences. A tissue phantom embedding a canal with flowing blood-mimicking fluid is imaged with a research programmable ultrasound system. A blood flow image is obtained and the flow characteristics are displayed for several pixels in the canal. Finally, a review of in vivo applications is proposed, showing examples in several organs such as carotids, kidney, thyroid, brain and heart.
Ultrasound imaging is one of the most commonly used imaging techniques in clinical practice and research activities. The combination of ultrasound wave emission in the biological tissues followed by the recording of the backscattered echoes allows the reconstruction of anatomical images, the so-called “B-Mode”. This method is perfectly adapted for soft tissue imaging, such as biological tissues, which typically permit the penetration of ultrasound over several centimeters, with a propagation speed of ≈1540 m/s. Depending on the center frequency of the ultrasound probe, images with a resolution from 30 µm to 1 mm are obtained. Furthermore, it is well known that the motion of an acoustic source, affects the physical characteristics of the associated waves. Particularly, the link between the frequency shifts of a wave relative to the speed of its source is described as the Doppler effect1, whose simplest manifestation is the changing siren’s pitch of a moving ambulance. Ultrasound imaging has long used this physical effect to observe the moving red blood cells2, and it proposes a variety of imaging modes commonly labelled “Doppler imaging”. These modes enable the assessment of blood flows in very different applications and organs, such as brain, heart, kidney or peripheral arteries.
Remarkably, most of the currently available ultrasound systems rely on the same technology, referred to as conventional ultrasound. The underlying principles are the following: an acoustic beam insonifies the field of view and is swept along the ultrasound transducer aperture. For each position of the beam, the echoes are recorded and converted into a line of the final image. By progressively moving the beam along the transducer, the whole field of view can be imaged line-per-line (Figure 1, left panel). This strategy was well adapted to the electrical constraints and computing power prevailing until the beginning of the 21st century. Nonetheless, it has several drawbacks. Among these, the final framerate is limited to a few hundred images per second by the beam scanning process. In terms of blood flow, this relatively low framerate affects the maximum flow velocities that can be detected, which is dictated by the sampling criteria of Shannon-Nyquist3. Moreover, conventional Doppler must deal with a complex tradeoff. In order to assess the blood flow velocity in a particular region of interest (ROI), several echoes coming from that ROI have to be successively recorded. This implies that the ultrasound beam is temporarily maintained in a fixed position. The longer the echo ensemble, the better the velocity estimation will be for that ROI. However, to produce a complete image of the field of view, the beam must scan the medium. Therefore, one can sense the conflict between these two constraints: holding the beam to precisely assess the velocity along one line, or moving the beam to produce an image. The different conventional Doppler modes (i.e., Color Doppler or Pulse Wave Doppler) directly reflect this tradeoff. Typically, the Color Doppler produces a low-fidelity flow map used for localizing the vessels4, and the Pulse Wave Doppler is then used to accurately quantify the flow in a previously identified vessel5.
These two limitations (low framerate and localization/quantification tradeoff) are overcome with very high-framerate emerging techniques. Among these, the synthetic aperture approach6 or the multiline transmit technique can be cited7. In this study, we focus on the so-called Ultrafast ultrasound method. Introduced two decades ago8,9,10, this method also relies on the emission/reception of ultrasounds, but with a radically different pattern. Indeed, instead of using a scanning focused beam, ultrafast imaging uses plane wave or diverging waves, which are able to insonify the field of view with a single emission. Following that single emission, the associated electronics is also able to receive and process the huge number of echoes originating from the whole field of view. At the end, an image can be reconstructed from a single emission/reception pattern11 (Figure 1, right panel). These unfocused emissions can have a low signal to noise ratio (SNR) due to the spread of the acoustical energy. This can be tackled by emitting several titled plane-waves (or diverging waves with different sources) and by adding the resulting images. This method is named “coherent compounding”12. Two major consequences arise. Firstly, the framerate only depends on the ultrasound time of flight and can reach typical values from 1 to 10 kHz. Secondly, this ensures the data continuity in both spatial and temporal dimensions, also referred to as spatiotemporal coherence. The conventional localization/quantification tradeoff is thus broken. This combination of a high framerate and spatiotemporal coherence has a tremendous impact on the ability to detect blood flows with ultrasound. As compared to conventional ultrasound, ultrafast ultrasound provides full characterization of the blood flow3. Practically, the user has access to the velocity time course in every pixel of the image, for the whole duration of the acquisition (typically ≈1 s), with a timescale given by the framerate (typically, a framerate of 5 kHz for a temporal resolution of 200 µs). This high framerate makes the method suitable for a wide range of application such as fast flow in moving organs like heart chambers13 or myocardium with the coronary micro-perfusion14. Furthermore, it has been shown that its spatiotemporal coherence strongly improves its ability to separate slow blood flow from background moving tissues, therefore increasing the sensitivity to micro-vascular flow15. This capacity gives access to the micro vasculature of the brain in both animals16 and humans17.
Hence, ultrafast ultrasound is well suited to image blood flow in a variety of situations. It is restricted to soft biological tissues and will be strongly affected by the presence of hard interfaces such as bones, or gas cavity such as the lung. The tuning of the physical parameters of the ultrasound sequence allows the study of both slow (down to 1 mm/s11,16) and fast flows (up to several m/s). A tradeoff exists between the spatial resolution and the depth of penetration. Typically, a resolution of 50 µm can be achieved at the cost of a penetration around 5 mm. Conversely, the penetration can be extended to 15-20 cm at the cost of a resolution of 1 mm. It is worth noting that most ultrafast scanners such as the one used in this article only provide 2D images.
Here, we propose a simple protocol to introduce the concept of Ultrafast Doppler imaging, using a programmable research ultrasound scanner and Doppler phantom mimicking a vessel (artery or vein) embedded in biological tissue.
1. Doppler phantom preparation setup (Figure 2A)
2. Ultrafast ultrasound scanner setup (Figure 2A)
3. Ultrasound sequence programming
4. Probe positioning and data acquisition
5. Image reconstruction (Figure 2B)
6. Clutter filtering (Figure 2C)
NOTE: For steps 6-7, see the Matlab script provided in the Supplementary Material.
7. Flow visualization and velocity measurements (Figure 2C)
The quality of the acquisition and the post-processing is firstly assessed by visual inspection. The shape of the canal must be clearly visible in the power Doppler image, and the tissue area must appear dark. If the power Doppler signal is not restricted to the canal, it can mean that either the clutter filter step went wrong (SVD threshold is too low), or the probe experienced a strong movement during the acquisition.
After visual inspection, the study of the spectrogram inside the canal can provide good information on the success or failure of the experiment. The spectrogram should be one-sided (all the values above or below zeros). If the spectrogram is two-sided, aliasing is present. In that case, either the flow is too fast, or the PRF is too low.
If these quality criteria are met, the blood velocities can be extracted from any ROI in the image (Figure 2C). Tuning the size of the ROI allows more or less averaging of the signals. The velocity time course of a given ROI can then be used for multiple analysis such as the computation of resistivity indexes19, wall shear stress estimation20, reactive hyperemia quantification14 and much more21,22.
Figure 3A-D shows the transposition of this protocol to various in vivo applications. In particular, the neonate brain acquisition (Figure 3B) exhibits vessels with very different flow characteristics, from small cortical venules and arterioles to the major pericallosal artery. Figure 3D illustrates the ability of ultrafast Doppler to extract blood flow signal in a strongly moving organ such as the myocardium.
Figure 1: Conventional and Ultrafast Ultrasound Imaging. Legend: (Left) Conventional imaging with focused emission. (Right) Ultrafast imaging with plane wave emission. (adapted from Villemain et al.22). Please click here to view a larger version of this figure.
Figure 2: Ultrafast Doppler protocol workflow. (A) Experimental setup including the ultrafast-enabled scanner and the Doppler flow phantom. The dashed rectangle on the phantom indicates the footprint of the ultrasound transducer. (B) Automated data acquisition chain and post-processing triggered by a simple user button-press. (C) (Top) Extraction of blood flow signal and suppression of the tissue background noise (“clutter filter”) and display of the blood spectrogram in any ROI of the field of view. (Bottom) Spectrogram showing the blood velocity distribution in the ROI at different time points. The mean velocity in the ROI is traced in dashed green. Please click here to view a larger version of this figure.
Figure 3: Power Doppler images. Ultrafast Doppler acquisitions on several organs. (A) Adult transplanted kidney, (B) Sagittal view of a human neonate brain, (C) Adult thyroid, (D) Intramural coronary vasculature in open-chest swine experiments, (E) 3D directional power Doppler of the carotid artery and the jugular vein of a healthy volunteer (blue = descending flow, red = ascending flow). Several spectrograms are extracted for different ROI. (A-C are adapted from Baranger et al.18, D is adapted from Maresca et al.14, E is adapted from Provost et al.23). For each acquisition, the center frequency, number of angles, PRF and maximum depth was tuned according to the situation. The dynamic range for panel A, B and C are respectively -27, -35 and -30 dB. It was not provided for panel D and E. Please click here to view a larger version of this figure.
Supplementary Materials. Please click here to download these files.
Several variations are possible around the main frame of this protocol.
Hardware concerns
If the user supplies its custom host computer, the motherboard and the computer’s case must have an available PCI express slot. The CPU must also have enough PCIe lanes to handle all the devices.
Probe selection
The ultrasound probe (also named transducer) is chosen according to the spatial resolution needed and to the geometry of the field of view. The higher the center frequency of the probe, the better the spatial resolution but the shorter the imaging depth. Linear, curved or phased-array probes provide field of view of respectively rectangular, circular sector and flat-top sector shapes.
Angle dependency
Ultrafast Doppler shares the same constraint as conventional Doppler regarding the dependency to the blood flow angle. Indeed, the underlying Doppler effect only allows the detection of movement in the axial direction, meaning toward the probe surface, or away from the probe. Hence, only the axial components of the blood scatterers velocity vectors are actually measured. The true velocity can be recovered by manually indicating the local angle of the blood flow with the vertical axis, but this angle cannot always be correctly assessed. In an extreme case where the flow is perfectly orthogonal to the vertical depth axis, the Doppler effect cannot be used to reliably measure the blood velocity. More advanced techniques can measure multiple projections of the velocity vectors on several axes and can ultimately reconstruct the true velocity vector. These angle independent approaches are referred to as vector flow imaging8,9,10,24.
Aliasing
The protocol described in this manuscript has several limitations. Firstly, the Shannon-Nyquist sampling theorem states that the maximum measurable frequency in the sampled signal cannot exceed half of the framerate. With 3 angles and a PRF of 12 kHz, the framerate is 4 kHz. Consequently, we can derive from Eq. 3 that the maximum detectable axial velocity is 30 cm/s. Considering the angle of the canal in the phantom, this framerate allows the detection of flow with velocities up to 96 cm/s. Velocities above this threshold will appear aliased in the Doppler spectrogram. For the presented setup, the peak velocities were ranging between 95 cm and 8 cm/s.
Clutter filter optimization
Blood flow visualization strongly relies on the ability to separate the blood signals from the slowly moving tissue background. Depending on the respiration or the sonographer’s hand motion, the tissue can move with speeds similar to slow blood flow. Therefore, the so-called “clutter filter” stage aims to cancel tissue signals. The ability to detect slow blood flows only relies on the efficiency of this clutter filter stage. It has been shown that leveraging the spatiotemporal coherence of ultrafast ultrasound strongly increases the outcome of these filters. The singular value decomposition filter described by Demene et al.15 is widely used. Optimization of that method18 or more complex algorithms such as high-order SVD25, principal component pursuit26, independent component analysis27 or other low-rank decomposition28 can improve the quality of the filtered data. It is worth mentioning that in the presented in vitro setup, the only source of clutter is the hand motion of the operator. In vivo, many other factors such as respiration and arterial pulsatility are likely to induce clutter that is more predominant. In these cases, the advanced SVD filter described in this protocol becomes of prime importance.
Spectrogram interpretation
Spectrograms are the most common tool to study blood flow characteristics in both conventional and ultrafast Doppler ultrasound. For each time point, the spectrogram shows in grayscale the velocity distribution inside the considered ROI. Complex flows such as non-laminar flows will hence naturally broaden this spectrum. However, this link between spectral broadening and velocity distribution is only true to a certain extent. It has been shown by several studies that the spectrogram width, also call spectral broadening, is indeed linked to the velocity distribution in the ROI but also to several geometrical parameters of the imaging system (array width, angles, etc)29,30,31. Therefore, while the spectrogram of a steady laminar and homogeneous flow should be a thin, flat line, it shows in practice a certain width that does not reflect the velocity distribution, but rather the geometry of the imaging setup. This potential pitfall can lead to incorrect velocity measurements. It is recommended to consider the mean velocity inside the ROI (Figure 1C dashed green line) to avoid these effects32.
3D acquisitions
The present protocol was realized with a standard linear array transducer, resulting in 2D images. Nevertheless, 3D acquisitions can be performed, either by mechanically scanning the medium with a motorized linear probe16, or by using row-column arrays33 or matrix arrays34. The drawbacks to these new methods are the high computing costs and, for matrix probes, the need of specific scanners. An example of 3D acquisition is shown in Figure 3E.
Safety issues
Most of the research ultrafast ultrasound scanner are not approved for clinical use. It is the responsibility of the experimenter to comply with the standards prevailing in their country, both in term of electrical safety and acoustic output. For the latter, the FDA standards35 and the IEC 62127-1 International Standard36 must be considered.
Conclusion
In this paper, we have proposed a standard protocol to image blood flow with Ultrafast Doppler. By practicing on a calibrated flow phantom, the operator can eventually check the accuracy of their measurements. The protocol allows the user to design and execute an ultrafast ultrasound acquisition using plane-wave compounding. Ultimately, a post-processing framework is described and provides the first tool to display the blood flow profile in any region of interest of the image.
The authors have nothing to disclose.
We would like to thank Shreya Shah for her proofreading and advice.
Blood-mimicking fluid | CIRS Inc, Norfolk, Virginia, USA | 069DTF | |
Doppler flow phantom | CIRS Inc, Norfolk, Virginia, USA | ATS523A | |
Matlab | MathWorks, Natick, Massachusetts, United States | ||
Peristaltic pump / Doppler flow pump | CIRS Inc, Norfolk, Virginia, USA | 769 | Include tubings and pulse dampener |
Transducer adpter | Verasonics, Kirkland, Washington, USA | UTA 408-GE | |
Ultrafast ultrasound research scanner | Verasonics, Kirkland, Washington, USA | Vantage 256 | |
Ultrasound probe/transducer | GE Healthcare | GE 9L-D |