High-Frame-Rate Volumetric Porcine Renal Vasculature Imaging

OBJECTIVE

The aim of this study was to assess the feasibility and imaging options of contrast-enhanced volumetric ultrasound kidney vasculature imaging in a porcine model using a prototype sparse spiral array.

METHODS

Transcutaneous freehand in vivo imaging of two healthy porcine kidneys was performed according to three protocols with different microbubble concentrations and transmission sequences. Combining high-frame-rate transmission sequences with our previously described spatial coherence beamformer, we determined the ability to produce detailed volumetric images of the vasculature. We also determined power, color and spectral Doppler, as well as super-resolved microvasculature in a volume. The results were compared against a clinical 2-D ultrasound machine.

RESULTS

Three-dimensional visualization of the kidney vasculature structure and blood flow was possible with our method. Good structural agreement was found between the visualized vasculature structure and the 2-D reference. Microvasculature patterns in the kidney cortex were visible with super-resolution processing. Blood flow velocity estimations were within a physiological range and pattern, also in agreement with the 2-D reference results.

CONCLUSION

Volumetric imaging of the kidney vasculature was possible using a prototype sparse spiral array. Reliable structural and temporal information could be extracted from these imaging results.

3-D Coherent Multi-Transducer Ultrasound Imaging with Sparse Spiral Arrays

Coherent multi-transducer ultrasound (CoMTUS) creates an extended effective aperture through the coherent combination of multiple arrays, which results in images with enhanced resolution, extended field-of-view, and higher sensitivity. The subwavelength localization accuracy of the multiple transducers required to coherently beamform the data is achieved by using the echoes backscattered from targeted points. In this study, CoMTUS is implemented and demonstrated for the first time in 3-D imaging using a pair of 256-element 2-D sparse spiral arrays, which keep the channel count low and limit the amount of data to be processed. The imaging performance of the method was investigated using both simulations and phantom tests. The feasibility of free-hand operation is also experimentally demonstrated. Results show that, in comparison to a single dense array system using the same total number of active elements, the proposed CoMTUS system improves spatial resolution (up to 10 times) in the direction where both arrays are aligned, contrast-to-noise-ratio (CNR, up to 46%), and generalized CNR (up to 15%). Overall, CoMTUS shows a narrower main lobe and higher contrast-to-noise ratio, which results in an increased dynamic range and better target detectability.

Surgeon eye lens dose monitoring in interventional neuroradiology, cardiovascular and radiology procedures

PURPOSE

This study investigated the radiation dose to surgeon eye lens for single procedure and normalised to exposure parameters for eight selected neuroradiology, cardiovascular and radiology interventional procedures.

METHODS

The procedures investigated were diagnostic study, Arteriovenous Malformations treatment (AVM) and aneurysm embolization for neuroradiology procedures, Coronary Angiography and Percutaneous Transluminal Coronary Angioplasty (CA-PTCA), Pacemaker and Implantable Cardioverter-Defibrillator implantation (PM-ICD), Endovascular Aortic Repair (EVAR) and Fenestrated Endovascular Aortic Repair (FEVAR) for cardiovascular and electrophysiology procedures. CT-guided lung biopsy was also monitored. All procedures were performed with table-mounted and ceiling-suspended shields (0.5 mm lead equivalent thickness), except for FEVAR and PM-ICD where only a table mounted shield was present, and CT-guided lung biopsy where no shield was used. Dose assessment was performed using a dosemeter positioned close to the most exposed eye of the surgeon, outside the protective eyewear.

RESULTS

The surgeon most exposed eye lens median H_p(3) equivalent dose for a single procedure, without protective eyewear contribution, was 18 μSv for neuroradiology diagnostic study, 62 μSv for AVM, 38 μSv for aneurysm embolization, 33 μSv for CA-PTCA, 39 μSv for PM-ICD, 49 μSv for EVAR, 2500 μSv for FEVAR, 153 μSv for CT-guided lung biopsy.

CONCLUSIONS

In interventional procedures, the 20 mSv/year dose limit for surgeon eye lens exposure might be exceeded if shields or protective eyewear are not used. Surgeon eye lens doses, normalised to single procedures and to exposure parameters, are a valuable tool for determining appropriate radiation protection measures and dedicated eye lens dosemeter assignment.

Design, Implementation, and Medical Applications of 2-D Ultrasound Sparse Arrays

An ultrasound sparse array consists of a sparse distribution of elements over a 2-D aperture. Such an array is typically characterized by a limited number of elements, which in most cases is compatible with the channel number of the available scanners. Sparse arrays represent an attractive alternative to full 2-D arrays that may require the control of thousands of elements through expensive application-specific integrated circuits (ASICs). However, their massive use is hindered by two main drawbacks: the possible beam profile deterioration, which may worsen the image contrast, and the limited signal-to-noise ratio (SNR), which may result too low for some applications. This article reviews the work done for three decades on 2-D ultrasound sparse arrays for medical applications. First, random, optimized, and deterministic design methods are reviewed together with their main influencing factors. Then, experimental 2-D sparse array implementations based on piezoelectric and capacitive micromachined ultrasonic transducer (CMUT) technologies are presented. Sample applications to 3-D (Doppler) imaging, super-resolution imaging, photo-acoustic imaging, and therapy are reported. The final sections discuss the main shortcomings associated with the use of sparse arrays, the related countermeasures, and the next steps envisaged in the development of innovative arrays.

Unfocused Field Analysis of a Density-Tapered Spiral Array for High-Volume-Rate 3-D Ultrasound Imaging

Spiral array transducers with a sparse 2-D aperture have demonstrated their potential in realizing 3-D ultrasound imaging with reduced data rates. Nevertheless, their feasibility in high-volume-rate imaging based on unfocused transmissions has yet to be established. From a metrology standpoint, it is essential to characterize the acoustic field of unfocused transmissions from spiral arrays not only to assess their safety but also to identify the root cause of imaging irregularities due to the array's sparse aperture. Here, we present a field profile analysis of unfocused transmissions from a density-tapered spiral array transducer (256 hexagonal elements, 220- [Formula: see text] element diameter, and 1-cm aperture diameter) through both simulations and hydrophone measurements. We investigated plane- and diverging-wave transmissions (five-cycle, 7.5-MHz pulses) from 0° to 10° steering for their beam intensity characteristics and wavefront arrival time profiles. Unfocused firings were also tested for B-mode imaging performance (ten compounded angles, -5° to 5° span). The array was found to produce unfocused transmissions with a peak negative pressure of 93.9 kPa at 2 cm depth. All transmissions steered up to 5° were free of secondary lobes within 12 dB of the main beam peak intensity. All wavefront arrival time profiles were found to closely match the expected profiles with maximum root-mean-squared errors of [Formula: see text] for plane wave (PW) and [Formula: see text] for diverging wave. The B-mode images showed good spatial resolution with a penetration depth of 22 mm in PW imaging. Overall, these results demonstrate that the density-tapered spiral array can facilitate unfocused transmissions below regulatory limits (mechanical index: 0.034; spatial-peak, pulse-average intensity: 0.298 W/cm2) and with suppressed secondary lobes while maintaining smooth wavefronts.

Sparse 2-D PZT-on-PCB Arrays With Density Tapering

Two-dimensional (2-D) arrays offer volumetric imaging capabilities without the need for probe translation or rotation. A sparse array with elements seeded in a tapering spiral pattern enables one-to-one connection to an ultrasound machine, thus allowing flexible transmission and reception strategies. To test the concept of sparse spiral array imaging, we have designed, realized, and characterized two prototype probes designed at 2.5-MHz low-frequency (LF) and 5-MHz high-frequency (HF) center frequencies. Both probes share the same electronic design, based on piezoelectric ceramics and rapid prototyping with printed circuit board substrates to wire the elements to external connectors. Different center frequencies were achieved by adjusting the piezoelectric layer thickness. The LF and HF prototype probes had 88% and 95% of working elements, producing peak pressures of 21 and 96 kPa/V when focused at 5 and 3 cm, respectively. The one-way -3-dB bandwidths were 26% and 32%. These results, together with experimental tests on tissue-mimicking phantoms, show that the probes are viable for volumetric imaging.

High Frame Rate Volumetric Imaging of Microbubbles Using a Sparse Array and Spatial Coherence Beamforming

Volumetric ultrasound imaging of blood flow with microbubbles enables a more complete visualization of the microvasculature. Sparse arrays are ideal candidates to perform volumetric imaging at reduced manufacturing complexity and cable count. However, due to the small number of transducer elements, sparse arrays often come with high clutter levels, especially when wide beams are transmitted to increase the frame rate. In this study, we demonstrate with a prototype sparse array probe and a diverging wave transmission strategy, that a uniform transmission field can be achieved. With the implementation of a spatial coherence beamformer, the background clutter signal can be effectively suppressed, leading to a signal to background ratio improvement of 25 dB. With this approach, we demonstrate the volumetric visualization of single microbubbles in a tissue-mimicking phantom as well as vasculature mapping in a live chicken embryo chorioallantoic membrane.

Experimental and simulation study of harmonic components generated by plane and focused waves

Although there is increasing interest in the use of plane waves (PW) in high-frame-rate imaging, not much experimental data is available about their behavior in terms of nonlinear propagation. This paper presents a detailed study of fundamental and harmonic components of the ultrasound beam associated to PW transmission from a linear array. Simulations and hydrophone measurements of PW propagation in water were performed and compared to the results obtained for focused waves (FWs) at various levels of peak negative pressure (PNP). Experimental results confirm that, at comparable PNP, the amplitudes of the harmonics reached by PWs are always higher, over extended regions, than those achieved with FW. For example, at MI = 0.2 the PW second harmonic turns out to be 9 dB higher at 25 mm depth (i.e. in the focal region), and 20 dB higher at 40 mm depth. Simulations additionally show that when ultrasound waves propagate through blood or muscle, the situation is in general reversed but, at low MI, the second harmonic amplitude can still be higher in PW than in FW. Furthermore, it is shown that increasing the array aperture size yields higher harmonic growth in PW compared to FW.

On the Depth-Dependent Accuracy of Plane-Wave-Based Vector Velocity Measurements With Linear Arrays

High-frame-rate vector Doppler methods are used to measure blood velocities over large 2-D regions, but their accuracy is often estimated over a short range of depths. This article thoroughly examines the dependence of velocity measurement accuracy on the target position. Simulations were carried out on flat and parabolic flow profiles, for different Doppler angles, and considering a 2-D vector flow imaging (2-D VFI) method based on plane wave transmission and speckle tracking. The results were also compared with those obtained by the reference spectral Doppler (SD) method. Although, as expected, the bias and standard deviation generally tend to worsen at increasing depths, the measurements also show the following. First, the errors are much lower for the flat profile (from ≈ -4 ± 3% at 20 mm to ≈ -17 ± 4% at 100 mm) than for the parabolic profile (from ≈ -4 ± 3% to ≈ -38 ±%). Second, only part of the relative estimation error is related to the inherent low resolution of the 2-D VFI method. For example, even for SD, the error bias increases (on average) from -0.7% (20 mm) to -17% (60 mm) up to -26% (100 mm). Third, conversely, the beam divergence associated with the linear array acoustic lens was found to have a great impact on the velocity measurements. By simply removing such lens, the average bias for 2-D VFI at 60 and 100 mm dropped down to -9.4% and -19.4%, respectively. In conclusion, the results indicate that the transmission beam broadening on the elevation plane, which is not limited by reception dynamic focusing, is the main cause of velocity underestimation in the presence of high spatial gradients.

Real-Time High Frame Rate Color Flow Mapping System

Plane wave (PW) transmission (TX) can be profitably used to improve the performance of color flow mapping (CFM) systems by increasing the autocorrelation ensemble length (EL) and/or the frame rate (FR). Although high-end scanners tend to include imaging schemes using PW TX and parallel receive beams, high frame rate (HFR) CFM has been so far experimentally implemented mostly through research platforms that transmit PWs and beamform/process the received channel data off-line. In this article, full real-time implementation of PW CFM with continuous-time clutter filtering and extended FR/EL is reported. The field-programmable gate arrays (FPGAs) and digital signal processors (DSPs) onboard the ULA-OP 256 research scanner were programmed to perform high-speed parallel beamforming and autocorrelation-based CFM processing, respectively. Different strategies were tested, in which the TX of PWs for CFM is either continuous or interleaved with the TX of packets of B-mode pulses. A fourth-order Chebyshev continuous-time high-pass filter with programmable cutoff frequency was implemented and its clutter rejection performance was positively compared with that obtained when operating on packet data. CFM FRs up to 575 were obtained. The possibility of programming the autocorrelation EL up to 64 permitted to detect flow with high sensitivity and accuracy (average relative errors down to 0.4% ± 8.4%). In vivo HFR movies are presented, showing the dynamics of flow in the common carotid artery, which highlight the presence of secondary flow components.

Real-Time 3-D Spectral Doppler Analysis With a Sparse Spiral Array

2-D sparse arrays may push the development of low-cost 3-D systems, not needing to control thousands of elements by expensive application-specific integrated circuits (ASICs). However, there is still some concern about their suitability in applications, such as Doppler investigation, which inherently involve poor signal-to-noise ratios (SNRs). In this article, a novel real-time 3-D pulsed-wave (PW) Doppler system, based on a 256-element 2-D spiral array, is presented. Coded transmission (TX) and matched filtering were implemented to improve the system SNR. Standard sonograms as well as multigate spectral Doppler (MSD) profiles, along lines that can be arbitrarily located in different planes, are presented. The performance of the system was assessed quantitatively on experimental data obtained from a straight tube flow phantom. An SNR increase of 11.4 dB was measured by transmitting linear chirps instead of standard sinusoidal bursts. For a qualitative assessment of the system performance in more realistic conditions, an anthropomorphic phantom of the carotid arteries was used. Finally, real-time B-mode and MSD images were obtained from healthy volunteers.

Virtual Real-Time for High PRF Multiline Vector Doppler on ULA-OP 256

The recent development of high-frame-rate (HFR) imaging/Doppler methods based on the transmission of plane or diverging waves has proposed new challenges to echographic data management and display. Due to the huge amount of data that need to be processed at very high speed, the pulse repetition frequency (PRF) is typically limited to hundreds hertz or few kilohertz. In Doppler applications, a PRF limitation may result unacceptable since it inherently translates to a corresponding limitation in the maximum detectable velocity. In this article, the ULA-OP 256 implementation of a novel ultrasound modality, called virtual real-time (VRT), is described. First, for a given HFR RT modality, the scanner displays the processed results while saving channel data into an internal buffer. Then, ULA-OP 256 switches to VRT mode, according to which the raw data stored in the buffer are immediately reprocessed by the same hardware used in RT. In the two phases, the ULA-OP 256 calculation power can be differently distributed to increase the acquisition frame rate or the quality of processing results. VRT was here used to extend the PRF limit in a multiline vector Doppler (MLVD) application. In RT, the PRF was maximized at the expense of the display quality; in VRT, data were reprocessed at a lower rate in a high-quality display format, which provides more detailed flow information. Experiments are reported in which the MLVD technique is shown capable of working at 16-kHz PRF, so that flow jet velocities higher up to 3 m/s can be detected.

Improving patient waiting time of centralized front office service in a regional hub hospital using the discrete event simulation model

BACKGROUND

Today, hospital rankings are based not only on basic clinical indicators, but even on quality service indicators such as patient waiting times. Improving these indicators is a very important issue for hospital management, so finding a solution to achieve it in a simple and effective way is one of the greatest goals.

OBJECTIVES

The aim of this article is to evaluate the use of a discrete event simulation model to improve healthcare processes and reduce waiting time of patients and hospital costs.

METHODS

The case study proposed in this paper is the reorganization of non-clinical front office operation for the patients (i.e. booking of exams, delivering medical reports, etc.) of the Careggi University Hospital of Florence, to optimize the utilization of the human resources and to improve performances of the process.

RESULTS

The development and validation of the model was made according to an analysis of real processes and data, pre and post implementation of model outcomes. The new organization shows a decrease of waiting times from an average value of 10 minutes and 37 seconds to 5 minutes and 57 seconds (-44%).

CONCLUSIONS

This paper shows that discrete event simulation could be a precise, cost-limited tool to optimize hospital processes and performance.

The effect of size range on ultrasound-induced translations in microbubble populations

Microbubble translations driven by ultrasound-induced radiation forces can be beneficial for applications in ultrasound molecular imaging and drug delivery. Here, the effect of size range in microbubble populations on their translations is investigated experimentally and theoretically. The displacements within five distinct size-isolated microbubble populations are driven by a standard ultrasound-imaging probe at frequencies ranging from 3 to 7 MHz, and measured using the multi-gate spectral Doppler approach. Peak microbubble displacements, reaching up to 10 μm per pulse, are found to describe transient phenomena from the resonant proportion of each bubble population. The overall trend of the statistical behavior of the bubble displacements, quantified by the total number of identified displacements, reveals significant differences between the bubble populations as a function of the transmission frequency. A good agreement is found between the experiments and theory that includes a model parameter fit, which is further supported by separate measurements of individual microbubbles to characterize the viscoelasticity of their stabilizing lipid shell. These findings may help to tune the microbubble size distribution and ultrasound transmission parameters to optimize the radiation-force translations. They also demonstrate a simple technique to characterize the microbubble shell viscosity, the fitted model parameter, from freely floating microbubble populations using a standard ultrasound-imaging probe.

A Comparison of Coherence-Based Beamforming Techniques in High-Frame-Rate Ultrasound Imaging With Multi-Line Transmission

One of the current challenges in ultrasound imaging is achieving higher frame rates, particularly in cardiac applications, where tracking the heart motion and other rapid events can provide potential valuable diagnostic information. The main drawback of ultrasound high-frame-rate strategies is that usually they partly sacrifice image quality in order to speed up the acquisition time. In particular, multi-line transmission (MLT), which consists in transmitting multiple ultrasound beams simultaneously in different directions, has been proven able to improve frame rates in echocardiography, unfortunately generating artifacts due to inter-beam crosstalk interferences. This work investigates the possibility to effectively suppress crosstalk artifacts in MLT while improving image quality by applying beamforming techniques based on backscattered signals spatial coherence. Several coherence-based algorithms (i.e., short-lag filtered-delay multiply and sum beamforming, coherence and generalized coherence factor, phase and sign coherence, and nonlinear beamforming with p th root compression) are implemented and compared, and their performance trends are evaluated when varying their design parameters. Indeed, experimental results of phantom and in vivo cardiac acquisitions demonstrate that this class of algorithms can provide significant benefits compared with delay and sum, well-suppressing artifacts (up to 48.5-dB lower crosstalk), and increasing image resolution (by up to 46.3%) and contrast (by up to 30 dB in terms of contrast ratio and 12.6% for generalized contrast-to-noise ratio) at the same time.

Doppler Ultrasound Technology for Fetal Heart Rate Monitoring: A Review

Fetal well-being is commonly assessed by monitoring the fetal heart rate (fHR). In clinical practice, the de facto standard technology for fHR monitoring is based on the Doppler ultrasound (US). Continuous monitoring of the fHR before and during labor is performed using a US transducer fixed on the maternal abdomen. The continuous fHR monitoring, together with simultaneous monitoring of the uterine activity, is referred to as cardiotocography (CTG). In contrast, for intermittent measurements of the fHR, a handheld Doppler US transducer is typically used. In this article, the technology of Doppler US for continuous fHR monitoring and intermittent fHR measurements is described, with emphasis on fHR monitoring for CTG. Special attention is dedicated to the measurement environment, which includes the clinical setting in which fHR monitoring is commonly performed. In addition, to understand the signal content of acquired Doppler US signals, the anatomy and physiology of the fetal heart and the surrounding maternal abdomen are described. The challenges encountered in these measurements have led to different technological strategies, which are presented and critically discussed, with a focus on the US transducer geometry, Doppler signal processing, and fHR extraction methods.

3-D Super-Resolution Ultrasound Imaging With a 2-D Sparse Array

High-frame-rate 3-D ultrasound imaging technology combined with super-resolution processing method can visualize 3-D microvascular structures by overcoming the diffraction-limited resolution in every spatial direction. However, 3-D super-resolution ultrasound imaging using a full 2-D array requires a system with a large number of independent channels, the design of which might be impractical due to the high cost, complexity, and volume of data produced. In this study, a 2-D sparse array was designed and fabricated with 512 elements chosen from a density-tapered 2-D spiral layout. High-frame-rate volumetric imaging was performed using two synchronized ULA-OP 256 research scanners. Volumetric images were constructed by coherently compounding nine-angle plane waves acquired at a pulse repetition frequency of 4500 Hz. Localization-based 3-D super-resolution images of two touching subwavelength tubes were generated from 6000 volumes acquired in 12 s. Finally, this work demonstrates the feasibility of 3-D super-resolution imaging and super-resolved velocity mapping using a customized 2-D sparse array transducer.

Spectral Doppler Measurements With 2-D Sparse Arrays

The 2-D sparse arrays, in which a few hundreds of elements are distributed on the probe surface according to an optimization procedure, represent an alternative to full 2-D arrays, including thousands of elements usually organized in a grid. Sparse arrays have already been used in B-mode imaging tests, but their application to Doppler investigations has not been reported yet. Since the sparsity of the elements influences the acoustic field, a corresponding influence on the mean frequency (Fm), bandwidth (BW), and signal-to-noise ratio (SNR) of the Doppler spectra is expected. This article aims to assess, by simulations and experiments, to what extent the use of a sparse rather than a full gridded 2-D array has an impact on spectral Doppler measurements. Parabolic flows were investigated by a 3 MHz, 1024-element gridded array and by a sparse array; the latter was obtained by properly selecting a subgroup of 256 elements from the full array. Simulations show that the mean Doppler frequency does not change between the sparse and the full array while there are significant differences on the BW (average reduction of 17.2% for the sparse array, due to different apertures of the two probes) and on the signal power (Ps) (22 dB, due to the different number of active elements). These results are confirmed by flow phantom experiments, which also highlight that the most critical difference between sparse and full gridded array in Doppler measurements is in terms of SNR (-16.8 dB).

High-Frame-Rate Tri-Plane Echocardiography With Spiral Arrays: From Simulation to Real-Time Implementation

Major cardiovascular diseases (CVDs) are associated with (regional) dysfunction of the left ventricle. Despite the 3-D nature of the heart and its dynamics, the assessment of myocardial function is still largely based on 2-D ultrasound imaging, thereby making diagnosis heavily susceptible to the operator's expertise. Unfortunately, to date, 3-D echocardiography cannot provide adequate spatiotemporal resolution in real-time. Hence, tri-plane imaging has been introduced as a compromise between 2-D and true volumetric ultrasound imaging. However, tri-plane imaging typically requires high-end ultrasound systems equipped with fully populated matrix array probes embedded with expensive and little flexible electronics for two-stage beamforming. This article presents an advanced ultrasound system for real-time, high frame rate (HFR), and tri-plane echocardiography based on low element count sparse arrays, i.e., the so-called spiral arrays. The system was simulated, experimentally validated, and implemented for real-time operation on the ULA-OP 256 system. Five different array configurations were tested together with four different scan sequences, including multi-line and planar diverging wave transmission. In particular, the former can be exploited to achieve, in tri-plane imaging, the same temporal resolution currently used in clinical 2-D echocardiography, at the expenses of contrast (-3.5 dB) and signal-to-noise ratio (SNR) (-8.7 dB). On the other hand, the transmission of planar diverging waves boosts the frame rate up to 250 Hz, but further compromises contrast (-10.5 dB), SNR (-9.7 dB), and lateral resolution (+46%). In conclusion, despite an unavoidable loss in image quality and sensitivity due to the limited number of elements, HFR tri-plane imaging with spiral arrays is shown to be feasible in real-time and may enable real-time functional analysis of all left ventricular segments of the heart.

Microbubble Radiation Force-Induced Translation in Plane-Wave Versus Focused Transmission Modes

Due to the primary radiation force, microbubble displacement has been observed previously in the focal region of single-element and array ultrasound probes. This effect has been harnessed to increase the contact between the microbubbles and targeted endothelium for drug delivery and ultrasound molecular imaging. In this study, microbubble displacements associated with plane-wave (PW) transmission are thoroughly investigated and compared to those obtained in focused-wave (FW) transmission over a range of pulse repetition frequencies, burst lengths (BLs), peak negative pressures, and transmission frequencies. In PW mode, the displacements, depending upon the experimental conditions, are in some cases consistently higher (e.g., by 28%, when the longest BL was used at PRF = 4 kHz), and the axial displacements are spatially more uniform compared to FW mode. Statistical analysis on the measured displacements reveals a slightly different frequency dependence of statistical quantities compared to transient peak microbubble displacements, which may suggest the need to consider the size range within the tested microbubble population.

Arterial wall shear rate response to reactive hyperaemia is markedly different between young and older humans

KEY POINTS

The vasodilatory response to reactive hyperaemia is impaired with advancing age, but it is unclear whether this is because of an altered wall shear rate (WSR) stimulus or an altered flow-mediated dilatation (FMD) response. Using new technology that allows detailed WSR measurement, we assessed the WSR-FMD response in healthy older people. Our data show that older people have a markedly altered and diminished WSR response to reactive hyperaemia compared to young people, but reduced WSR alone does not fully explain reduced FMD. In young people, WSR appears to be coupled to FMD but, by age ∼65 years, the arterial vasodilatory response has begun to uncouple from the WSR stimulus. These findings point to the importance and utility of comprehensively characterizing the WSR-FMD response when using reactive hyperaemia to assess vascular function, as well as giving new insight into the age-related alteration in vascular function.

ABSTRACT

The vasodilatory response to reactive hyperaemia is impaired with age, but it is unknown whether this is because of an altered wall shear rate (WSR) stimulus or an altered flow-mediated dilatation (FMD) response to the WSR stimulus. Inherent difficulties in measuring blood flow velocity close to the arterial wall have prevented detailed assessment of the WSR-FMD response. Using an enhanced multigate spectral Doppler ultrasound system (ultrasound advanced open platform), we aimed to produce new data on the WSR-FMD relationship in healthy older adults. Sixty healthy people, comprising 28 young (27.5 ± 5.5 years) and 32 older (64.9 ± 3.7 years) individuals, underwent FMD assessment. Raw data were post-processed using custom-designed software to obtain WSR and diameter parameters. The data revealed that older people have a much altered and diminished WSR response to reactive hyperaemia compared to younger people [e.g. WSR peak: 622 (571-673) vs. 443 (396-491) 1/s in young and older respectively; P < 0.05]. However, reduced WSR alone does not appear to fully explain the reduced FMD response in older people because associations between WSR and FMD were few and weak. This was in contrast to young adults, where associations were strong. We conclude that WSR during FMD is much altered and diminished in older people, and there appears to be an 'uncoupling' of WSR from FMD in older people that may reflect a loss of precision in the reactive hyperaemia stimulus-response relationship. These findings also point to the importance and utility of comprehensively characterizing the WSR-FMD response when using reactive hyperaemia to assess vascular function.

Continuous Simultaneous Recording of Brachial Artery Distension and Wall Shear Rate: A New Boost for Flow-Mediated Vasodilation

Vascular ultrasound has been extensively applied in the clinical setting to noninvasively assess the endothelial function by means of the so-called brachial artery flow mediated dilation (FMD). Despite the usefulness in large-scale epidemiological studies, this approach has revealed some pitfalls for assessing vascular physiology and health in individual subjects. Mainly, a reliable FMD examination should be based on the simultaneous and reliable measurement of both the stimulus, i.e., the wall shear rate (WSR), and the response, i.e., the diameter change. However, multiple technical, practical, and methodological challenges must be faced to meet this goal. In this work, we present the technical developments needed to implement a system to enable the extensive and reliable clinical ultrasound FMD examination. It integrates both a hardware part, i.e., an upgraded version of the ultrasound advanced open platform (ULA-OP), and a software part, i.e., a signal processing and data analysis platform. The system was applied for a two-center pilot clinical study on 35 young and healthy volunteers. Therefore, we present here the results of a statistical analysis on magnitude, time-course, and kinetic parameters of WSR and diameter trends that allowed us to accurately explore the vasodilatory response to the dynamic WSR changes. Our observations demonstrate that a direct and accurate estimation of WSR stimulus by multigate spectral Doppler allows understanding brachial artery vasodilatory response to reactive hyperemia. Drawing inferences on WSR stimulus from the diameter response along with an inaccurate estimation of WSR may cause further uncertainties for the accurate interpretation of the FMD response.

Plane-Wave Contrast Imaging: A Radiation Force Point of View

Radiation force is known to produce microbubble axial displacements by an amount that depends on the transmit burst frequency, amplitude, and length, as well as the pulse repetition frequency (PRF). In the standard focused-imaging mode, the actual PRF experienced by each microbubble is low, because it is of the order of the frame rate (i.e., usually tens of Hertz). In the plane-wave imaging mode, however, the actual PRF is considerably higher, as it is equivalent to the transmit PRF (kiloHertz range). Furthermore, the radiation pressure is expected to be almost uniform over the field of view, and typically lower than the peak pressure experienced in the focused transmit (TX) mode. We have experimentally investigated the possible effects of radiation force in the plane-wave mode. Here, we report on preliminary findings that show that the acoustic radiation force is negligible only at lower TX levels. At higher TX amplitudes, the bubble displacements due to radiation force are comparable to those obtained for focused waves at the same PRF. In addition, the radiation force is nearly uniform over the field of view and increases as the TX burst central frequency approaches the resonance frequency of size-isolated microbubbles.

Real-Time High-Frame-Rate Cardiac B-Mode and Tissue Doppler Imaging Based on Multiline Transmission and Multiline Acquisition

Cardiovascular diseases, the leading cause of death in the world, are often associated with the dysfunction of the left ventricle. Even if, in clinical practice, the myocardial function is often assessed through visual wall motion scoring on B-mode images, quantitative techniques have been introduced, e.g., ultrasound tissue Doppler imaging (TDI). However, this technique suffers from the limited frame rate of currently available imaging techniques that needs to be balanced with the field of view. High-frame-rate (HFR) cardiac imaging has been recently tested off-line by simultaneously transmitting multiple focused beams into different directions and acquiring raw channel data into a PC. Several image lines were then reconstructed from the echoes of each transmission (TX) event. The same approach has been used to increase the TDI frame rate without restricting the field of view. This paper demonstrates the real-time feasibility of multiline TX and acquisition methods for both HFR cardiac B-mode and TDI. These approaches have been implemented on the ULA-OP 256 research scanner, by taking care that the related resources were optimally exploited for these new applications. The obtainable performance in terms of image quality and frame rate has also been investigated. Experiments performed with a 128-element phased array probe show, for the first time, that real-time B-mode imaging is feasible at up to 1150 Hz without significant reduction in image quality or field of view. The implementation of a real-time TDI algorithm allowed obtaining TDI images with a frame rate of 288 Hz for a 90°-wide field of view. Finally, in vivo examples demonstrate the feasibility and the suitability of the method in clinical studies.

Ultrasound Open Platforms for Next-Generation Imaging Technique Development

Open platform (OP) ultrasound systems are aimed primarily at the research community. They have been at the forefront of the development of synthetic aperture, plane wave, shear wave elastography, and vector flow imaging. Such platforms are driven by a need for broad flexibility of parameters that are normally preset or fixed within clinical scanners. OP ultrasound scanners are defined to have three key features including customization of the transmit waveform, access to the prebeamformed receive data, and the ability to implement real-time imaging. In this paper, a formative discussion is given on the development of OPs from both the research community and the commercial sector. Both software- and hardware-based architectures are considered, and their specifications are compared in terms of resources and programmability. Software-based platforms capable of real-time beamforming generally make use of scalable graphics processing unit architectures, whereas a common feature of hardware-based platforms is the use of field-programmable gate array and digital signal processor devices to provide additional on-board processing capacity. OPs with extended number of channels (>256) are also discussed in relation to their role in supporting 3-D imaging technique development. With the increasing maturity of OP ultrasound scanners, the pace of advancement in ultrasound imaging algorithms is poised to be accelerated.

Experimental evaluation of ultrasound higher-order harmonic imaging with Filtered-Delay Multiply And Sum (F-DMAS) non-linear beamforming

Tissue Harmonic Imaging (THI) mode is currently one of the preferred choices by the clinicians for its ability to provide enhanced ultrasound images, thanks to the use of the second harmonic component of backscattered echoes. This paper aims at investigating whether the combination of THI with Filtered-Delay Multiply And Sum (F-DMAS) beamforming can provide further improvements in image quality. F-DMAS is a new non-linear beamformer, which, similarly to THI, is based on the use of the second harmonics of beamformed signals and is known to increase image contrast resolution and noise rejection. Thus, we have first compared the images obtained by using F-DMAS and the standard Delay And Sum (DAS) beamformers when only the second harmonics of the received signals was selected. Moreover, possible improvements brought about by other harmonic components generated by the combined use of the fundamental plus second harmonics and F-DMAS beamforming have been explored. Experimental results demonstrate that, as compared to standard harmonic imaging with DAS, THI and F-DMAS can be joined to improve the -20 dB lateral resolution up to 1 mm, the contrast ratio up to 12 dB on a cyst-phantom and up to 9 dB on in vivo images.

Real-Time Blood Velocity Vector Measurement Over a 2-D Region

Quantitative blood velocity measurements, as currently implemented in commercial ultrasound scanners, are based on pulsed-wave (PW) spectral Doppler and are limited to detect the axial component of the velocity in a single sample volume. On the other hand, vector Doppler methods produce angle-independent estimates by, e.g., combining the frequency shifts measured from different directions. Moreover, thanks to the transmission of plane waves, the investigation of a 2-D region is possible with high temporal resolution, but, unfortunately, the clinical use of these methods is hampered by the massive calculation power required for their real-time execution. In this paper, we present a novel approach based on the transmission of plane waves and the simultaneous reception of echoes from 16 distinct subapertures of a linear array probe, which produces eight lines distributed over a 2-D region. The method was implemented on the ULAO-OP 256 research scanner and tested both in phantom and in vivo. A continuous real-time refresh rate of 36 Hz was achieved in duplex combination with a standard B-mode at pulse repetition frequency of 8 kHz. Accuracies of -11% on velocity and of 2°on angle measurements have been obtained in phantom experiments. Accompanying movies show how the method improves the quantitative measurements of blood velocities and details the flow configurations in the carotid artery of a volunteer.

Multi Line Transmit Beamforming Combined With Adaptive Apodization

Increased frame rate is of high importance to cardiac diagnostic imaging as it enables examination of fast events during the cardiac cycle and improved quantitative analysis, such as speckle tracking. Multi-line transmission (MLT) is one of the methods proposed for this purpose. In contrast to the single-line transmission (SLT), where one focused beam is sent in each direction, MLT beams are simultaneously transmitted and focused in several (2,4,6..) directions improving the framerate accordingly. The simultaneous transmission is known to cause cross-talk artifacts due to the interference between the main-lobes and the side-lobes of the transmitted and received beams. Usually, the artifacts are attenuated using a Tukey window apodization, but the lateral resolution is degraded. Several other methods, such as minimum variance beamforming and filtered delay multiply and sum beamforming were proposed to deal with these artifacts.The assumption examined in this study is that a receive apodization can be chosen adaptively from a number of apodization windows in order to provide better artifact rejection and to increase the spatial resolution. The entire study was performed on experimental MLT dataset including wire and tissue mimicking phantoms, as well as in vivo cardiac data. The results demonstrate that application of a predefined apodization bank outperforms Tukey windowing alone, in terms of both resolution and receive crosstalk artifact rejection. Moreover, the achieved spatial resolution is superior to the non-apodized SLT, as measured from wire phantoms. The proposed method can also be combined with wider transmit beams, suitable for multi line acquisition.

Brachial artery vasodilatory response and wall shear rate determined by multigate Doppler in a healthy young cohort

Wall shear rate (WSR) is an important stimulus for the brachial artery flow-mediated dilation (FMD) response. However, WSR estimation near the arterial wall by conventional Doppler is inherently difficult. To overcome this limitation, we utilized multigate Doppler to accurately determine the WSR stimulus near the vessel wall simultaneously with the FMD response using an integrated FMD system [Ultrasound Advanced Open Platform (ULA-OP)]. Using the system, we aimed to perform a detailed analysis of WSR-FMD response and establish novel WSR parameters in a healthy young population. Data from 33 young healthy individuals (27.5 ± 4.9 yr, 19 females) were analyzed. FMD was assessed with reactive hyperemia using ULA-OP. All acquired raw data were postprocessed using custom-designed software to obtain WSR and diameter parameters. The acquired velocity data revealed that nonparabolic flow profiles within the cardiac cycle and under different flow states, with heterogeneity between participants. We also identified seven WSR magnitude and four WSR time-course parameters. Among them, WSR area under the curve until its return to baseline was the strongest predictor of the absolute ( R² = 0.25) and percent ( R² = 0.31) diameter changes in response to reactive hyperemia. For the first time, we identified mono- and biphasic WSR stimulus patterns within our cohort that produced different magnitudes of FMD response [absolute diameter change: 0.24 ± 0.10 mm (monophasic) vs. 0.17 ± 0.09 mm (biphasic), P < 0.05]. We concluded that accurate and detailed measurement of the WSR stimulus is important to comprehensively understand the FMD response and that this advance in current FMD technology could be important to better understand vascular physiology and pathology. NEW & NOTEWORTHY An estimation of wall shear rate (WSR) near the arterial wall by conventional Doppler ultrasound is inherently difficult. Using a recently developed integrated flow-mediated dilation ultrasound system, we were able to accurately estimate WSR near the wall and identified a number of novel WSR variables that may prove to be useful in the measurement of endothelial function, an important biomarker of vascular physiology and disease.

Architecture of an Ultrasound System for Continuous Real-Time High Frame Rate Imaging

High frame rate (HFR) imaging methods based on the transmission of defocused or plane waves rather than focused beams are increasingly popular. However, the production of HFR images poses severe requirements both in the transmission and the reception sections of ultrasound scanners. In particular, major technical difficulties arise if the images must be continuously produced in real-time, i.e., without any acquisition interruption nor loss of data. This paper presents the implementation of the real-time HFR-compounded imaging application in the ULA-OP 256 research platform. The beamformer sustains an average output sample rate of 470 MSPS. This allows continuously producing coherently compounded images, each of 64 lines by 1280 depths (here corresponding to 15.7 mm width and 45 mm depth, respectively), at frame rates up to 5.3 kHz. Imaging tests addressed to evaluate the achievable speed and quality performance were conducted on phantom. Results obtained by real-time compounding frames obtained with different numbers of steering angles between +7.5° and -7.5° are presented.

High Frame-Rate, High Resolution Ultrasound Imaging With Multi-Line Transmission and Filtered-Delay Multiply And Sum Beamforming

Multi-Line Transmission (MLT) was recently demonstrated as a valuable tool to increase the frame rate of ultrasound images. In this approach, the multiple beams that are simultaneously transmitted may determine cross-talk artifacts that are typically reduced, although not eliminated, by the use of Tukey apodization on both transmission and reception apertures, which unfortunately worsens the image lateral resolution. In this paper we investigate the combination, and related performance, of Filtered-Delay Multiply And Sum (F-DMAS) beamforming with MLT for high frame-rate ultrasound imaging. F-DMAS is a non-linear beamformer based on the computation of the receive aperture spatial autocorrelation, which was recently proposed for use in ultrasound B-mode imaging by some of the authors. The main advantages of such beamformer are the improved contrast resolution, obtained by lowering the beam side lobes and narrowing the main lobe, and the increased noise rejection. This study shows that in MLT images, compared to standard Delay And Sum (DAS) beamforming including Tukey apodization, F-DMAS beamforming yields better suppression of cross-talk and improved lateral resolution. The method's effectiveness is demonstrated by simulations and phantom experiments. Preliminary in vivo cardiac images also show that the frame rate can be improved up to 8-fold by combining F-DMAS and MLT without affecting the image quality.

Wideband 2-D Array Design Optimization With Fabrication Constraints for 3-D US Imaging

Ultrasound (US) 2-D arrays are of increasing interest due to their electronic steering capability to investigate 3-D regions without requiring any probe movement. These arrays are typically populated by thousands of elements that, ideally, should be individually driven by the companion scanner. Since this is not convenient, the so-called microbeamforming methods, yielding a prebeamforming stage performed in the probe handle by suitable custom integrated circuits, have so far been implemented in a few commercial high-end scanners. A possible approach to implement relatively cheap and efficient 3-D US imaging systems is using 2-D sparse arrays in which a limited number of elements can be coupled to an equal number of independent transmit/receive channels. In order to obtain US beams with adequate characteristics all over the investigated volume, the layout of such arrays must be carefully designed. This paper provides guidelines to design, by using simulated annealing optimization, 2-D sparse arrays capable of fitting specific applications or fabrication/implementation constraints. In particular, an original energy function based on multidepth 3-D analysis of the beam pattern is also exploited. A tutorial example is given, addressed to find the N _e elements that should be activated in a 2-D fully populated array to yield efficient acoustic radiating performance over the entire volume. The proposed method is applied to a 32 ×32 array centered at 3 MHz to select the 128, 192, and 256 elements that provide the best acoustic performance. It is shown that the 256-element optimized array yields sidelobe levels even lower (by 5.7 dB) than that of the reference 716-element circular and (by 10.3 dB) than that of the reference 1024-element array.

Wall Shear Rate Measurement: Validation of a New Method Through Multiphysics Simulations

Wall shear stress is known to affect the vessel endothelial function and to be related to important pathologies like the development of atherosclerosis. It is defined as the product of the blood viscosity by the blood velocity gradient at the wall position, i.e., the wall shear rate (WSR). The WSR measurement is particularly challenging in important cardiovascular sites, like the carotid bifurcation, because of the related complex flow configurations characterized by high spatial and temporal gradients, wall movement, and clutter noise. Moreover, accuracy of any method for WSR measurement can be effectively tested only if reliable gold standard WSR values, considering all the aforementioned disturbing effects, are available. Unfortunately, these requirements are difficult to achieve in a physical phantom, so that the accuracy test of the novel WSR measurement methods was so far limited to straight pipes and/or similar idealistic configurations. In this paper, we propose a new method for WSR measurement and its validation based on a mathematical model of the carotid bifurcation, which, exploiting fluid-structure simulations, is capable of reproducing realistic flow configuration, wall movement, and clutter noise. In particular, the profile near the wall, not directly measurable because affected by clutter, is estimated through a power-law fitting and compared with the gold standard provided by the model. In this condition, the WSR measurements featured an accuracy of ±20 %. A preliminary test on a volunteer confirmed the feasibility of the WSR method for in vivo application.

2-D Ultrasound Sparse Arrays Multidepth Radiation Optimization Using Simulated Annealing and Spiral-Array Inspired Energy Functions

Full matrix arrays are excellent tools for 3-D ultrasound imaging, but the required number of active elements is too high to be individually controlled by an equal number of scanner channels. The number of active elements is significantly reduced by the sparse array techniques, but the position of the remaining elements must be carefully optimized. This issue is faced here by introducing novel energy functions in the simulated annealing (SA) algorithm. At each iteration step of the optimization process, one element is freely translated and the associated radiated pattern is simulated. To control the pressure field behavior at multiple depths, three energy functions inspired by the pressure field radiated by a Blackman-tapered spiral array are introduced. Such energy functions aim at limiting the main lobe width while lowering the side lobe and grating lobe levels at multiple depths. Numerical optimization results illustrate the influence of the number of iterations, pressure measurement points, and depths, as well as the influence of the energy function definition on the optimized layout. It is also shown that performance close to or even better than the one provided by a spiral array, here assumed as reference, may be obtained. The finite-time convergence properties of SA allow the duration of the optimization process to be set in advance.

Coded Spectral Doppler Imaging: From Simulation to Real-Time Processing

Transmission of coded pulses and matched receive filtering can improve the ultrasound imaging penetration depth while preserving the axial resolution. This paper shows that the pulse compression technique may be integrated in a low-cost scanner to be profitably used also in spectral Doppler investigations. By operating on beamformed, demodulated, and down-sampled data in the frequency domain, a single digital signal processor is proved sufficient to perform both pulse compression and multigate spectral Doppler algorithms in real time. Simulations, phantom, and in vivo experiments demonstrate that the transmission of (2.5 or [Formula: see text] long) linear frequency-modulated chirps with bandwidths over the range 1.6-5.4 MHz, rather than of corresponding sine-burst pulses, provides signal-to-noise ratio (SNR) improvements very close to theory. Even in the presence of selective tissue attenuation, SNR gains up to 11 and 13.3 dB have been obtained for the short and the longer chirp, respectively. This may be important in clinical Doppler applications where the needed penetration depth is not achieved with sufficient SNR unless very long bursts are transmitted.

ULA-OP 256: A 256-Channel Open Scanner for Development and Real-Time Implementation of New Ultrasound Methods

Open scanners offer an increasing support to the ultrasound researchers who are involved in the experimental test of novel methods. Each system presents specific performance in terms of number of channels, flexibility, processing power, data storage capability, and overall dimensions. This paper reports the design criteria and hardware/software implementation details of a new 256-channel ultrasound advanced open platform. This system is organized in a modular architecture, including multiple front-end boards, interconnected by a high-speed (80 Gb/s) ring, capable of finely controlling all transmit (TX) and receive (RX) signals. High flexibility and processing power (equivalent to 2500 GFLOP) are guaranteed by the possibility of individually programming multiple digital signal processors and field programmable gate arrays. Eighty GB of on-board memory are available for the storage of prebeamforming, postbeamforming, and baseband data. The use of latest generation devices allowed to integrate all needed electronics in a small size ( 34 cm ×30 cm ×26 cm). The system implements a multiline beamformer that allows obtaining images of 96 lines by 2048 depths at a frame rate of 720 Hz (expandable to 3000 Hz). The multiline beamforming capability is also exploited to implement a real-time vector Doppler scheme in which a single TX and two independent RX apertures are simultaneously used to maintain the analysis over a full pulse repetition frequency range.

Wide-Angle Tissue Doppler Imaging at High Frame Rate Using Multi-Line Transmit Beamforming: An Experimental Validation In Vivo

Color tissue Doppler imaging (TDI) is a well-established methodology to assess local myocardial motion/deformation. Typically, a frame rate of ∼ 200 Hz can be achieved by imaging a narrow sector (∼ 30°, covering one cardiac wall) at moderate line density, using a dedicated pulse sequence and multi-line acquisition. However, a wide angle field-of-view is required in some clinical applications to image the whole left ventricle, which implies a drop in temporal resolution. Hereto, the aim of this study was to propose a novel imaging sequence using a multi-line transmit (MLT) beamforming approach to achieve high frame rate color TDI while preserving a wide field-of-view (i.e., 90° sector). To this end, a color MLT-TDI sequence achieving a frame rate of 208 Hz with a 90°-sector was implemented on an ultrasound experimental scanner interleaved with a conventional color TDI sequence achieving the same frame rate but only with a 22.5°-sector. Using this setup, the septal wall of 9 healthy volunteers was imaged and the corresponding velocity was extracted. The M-mode velocity images and the velocity profiles obtained from the MLT-TDI images presented physiologic patterns, very similar to those from conventional TDI. Moreover, for the peak systolic/diastolic velocities, good agreement and strong correlation between MLT-TDI and conventional TDI were found. The results thus demonstrate the feasibility of the novel MLT based TDI methodology to achieve high frame rate color TDI without compromising the field-of-view. This may open the opportunity to simultaneously assess regional myocardial function of the whole left ventricle at high temporal resolution.

Detection of Contrast Agents: Plane Wave Versus Focused Transmission

Ultrasound contrast agent (UCA) imaging provides a cost-effective diagnostic tool to assess tissue perfusion and vascular pathologies. However, excessive transmission (TX) levels may negatively impact both uniform diffusion and survival rates of contrast agents, limiting their density and thus their echogenicity. Contrast detection methods with both high sensitivity and low-contrast destruction rate are thus essential to maintain diagnostic capabilities. Plane-wave TX with a high number of compounding angles has been suggested to produce good quality images at pressure levels that do not destroy UCA. In this paper, we performed a quantitative evaluation of detection efficacy of flowing UCA with either traditional focused scanning or ultrafast plane-wave imaging. Amplitude modulation (AM) at nondestructive pressure levels was implemented in the ULA-OP ultrasound research platform. The influence of the number of compounding angles, peak-negative pressure, and flow speed on the final image quality was investigated. Results show that the images obtained by compounding multiple angled plane waves offer a greater contrast (up to a 12-dB increase) with respect to focused AM. This increase is attributed mainly to noise reduction caused by the coherent summation in the compounding step. Additionally, we show that highly sensitive detection is already achieved with a limited compounding number ( ), thus suggesting the feasibility of continuous contrast monitoring at a high frame rate. This capability is essential to properly detect contrast agents flowing at high speed, as an excessive angle compounding is shown to be destructive for the contrast signal, as the UCA motion quickly causes loss of correlation between consecutive echoes.

Plane-wave transverse oscillation for high-frame-rate 2-D vector flow imaging

Transverse oscillation (TO) methods introduce oscillations in the pulse-echo field (PEF) along the direction transverse to the ultrasound propagation direction. This may be exploited to extend flow investigations toward multidimensional estimates. In this paper, the TOs are coupled with the transmission of plane waves (PWs) to reconstruct high-framerate RF images with bidirectional oscillations in the pulse-echo field. Such RF images are then processed by a 2-D phase-based displacement estimator to produce 2-D vector flow maps at thousands of frames per second. First, the capability of generating TOs after PW transmissions was thoroughly investigated by varying the lateral wavelength, the burst length, and the transmission frequency. Over the entire region of interest, the generated lateral wavelengths, compared with the designed ones, presented bias and standard deviation of -3.3 ± 5.7% and 10.6 ± 7.4% in simulations and experiments, respectively. The performance of the ultrafast vector flow mapping method was also assessed by evaluating the differences between the estimated velocities and the expected ones. Both simulations and experiments show overall biases lower than 20% when varying the beam-to-flow angle, the peak velocity, and the depth of interest. In vivo applications of the method on the common carotid and the brachial arteries are also presented.

Density-tapered spiral arrays for ultrasound 3-D imaging

The current high interest in 3-D ultrasound imaging is pushing the development of 2-D probes with a challenging number of active elements. The most popular approach to limit this number is the sparse array technique, which designs the array layout by means of complex optimization algorithms. These algorithms are typically constrained by a few steering conditions, and, as such, cannot guarantee uniform side-lobe performance at all angles. The performance may be improved by the ungridded extensions of the sparse array technique, but this result is achieved at the expense of a further complication of the optimization process. In this paper, a method to design the layout of large circular arrays with a limited number of elements according to Fermat's spiral seeds and spatial density modulation is proposed and shown to be suitable for application to 3-D ultrasound imaging. This deterministic, aperiodic, and balanced positioning procedure attempts to guarantee uniform performance over a wide range of steering angles. The capabilities of the method are demonstrated by simulating and comparing the performance of spiral and dense arrays. A good trade-off for small vessel imaging is found, e.g., in the 60λ spiral array with 1.0λ elements and Blackman density tapering window. Here, the grating lobe level is -16 dB, the lateral resolution is lower than 6λ the depth of field is 120λ and, the average contrast is 10.3 dB, while the sensitivity remains in a 5 dB range for a wide selection of steering angles. The simulation results may represent a reference guide to the design of spiral sparse array probes for different application fields.

A real-time chirp-coded imaging system with tissue attenuation compensation

In ultrasound imaging, pulse compression methods based on the transmission (TX) of long coded pulses and matched receive filtering can be used to improve the penetration depth while preserving the axial resolution (coded-imaging). The performance of most of these methods is affected by the frequency dependent attenuation of tissue, which causes mismatch of the receiver filter. This, together with the involved additional computational load, has probably so far limited the implementation of pulse compression methods in real-time imaging systems. In this paper, a real-time low-computational-cost coded-imaging system operating on the beamformed and demodulated data received by a linear array probe is presented. The system has been implemented by extending the firmware and the software of the ULA-OP research platform. In particular, pulse compression is performed by exploiting the computational resources of a single digital signal processor. Each image line is produced in less than 20 μs, so that, e.g., 192-line frames can be generated at up to 200 fps. Although the system may work with a large class of codes, this paper has been focused on the test of linear frequency modulated chirps. The new system has been used to experimentally investigate the effects of tissue attenuation so that the design of the receive compression filter can be accordingly guided. Tests made with different chirp signals confirm that, although the attainable compression gain in attenuating media is lower than the theoretical value expected for a given TX Time-Bandwidth product (BT), good SNR gains can be obtained. For example, by using a chirp signal having BT=19, a 13 dB compression gain has been measured. By adapting the frequency band of the receiver to the band of the received echo, the signal-to-noise ratio and the penetration depth have been further increased, as shown by real-time tests conducted on phantoms and in vivo. In particular, a 2.7 dB SNR increase has been measured through a novel attenuation compensation scheme, which only requires to shift the demodulation frequency by 1 MHz. The proposed method characterizes for its simplicity and easy implementation.

Comparison of carotid artery blood velocity measurements by vector and standard Doppler approaches

Although severely affected by the angle dependency, carotid artery peak systolic velocity measurements are widely used for assessment of stenosis. In this study, blood peak systolic velocities in the common and internal carotid arteries of both healthy volunteers and patients with internal carotid artery stenosis were measured by two vector Doppler (VD) methods and compared with measurements obtained with the conventional spectral Doppler approach. Although the two VD techniques were completely different (using the transmission of focused beams and plane waves, respectively), the measurement results indicate that these techniques are nearly equivalent. The peak systolic velocities measured in 22 healthy common carotid arteries by the two VD techniques were very close (according to Bland-Altman analysis, the average difference was 3.2%, with limits of agreement of ± 8.6%). Application of Bland-Altman analysis to comparison of either VD technique with the spectral Doppler method provided a 21%-25% average difference with ± 13%-15% limits of agreement. Analysis of the results obtained from 15 internal carotid arteries led to similar conclusions, indicating significant overestimation of peak systolic velocity with the spectral Doppler method. Inter- and intra-operator repeatability measurements performed in a group of 8 healthy volunteers provided equivalent results for all of the methods (coefficients of variability in the range 2.7%-6.9%), even though the sonographers were not familiar with the VD methods. The results of this study suggest that the introduction of vector Doppler methods in commercial machines may finally be considered mature and capable of overcoming the angle-dependent overestimation typical of the standard spectral Doppler approach.

Profound increase in longitudinal displacements of the porcine carotid artery wall can take place independently of wall shear stress: a continuation report

The mechanisms underlying longitudinal displacements of the arterial wall, that is, displacements of the wall layers along the artery, and the resulting intramural shear strain remain largely unknown. We have already found that these displacements undergo profound changes in response to catecholamines. Wall shear stress, closely related to wall shear rate, represents the viscous drag exerted on the vessel wall by flowing blood. The aim of the work described here was to study possible relations between the wall shear rate and the longitudinal displacements. We investigated the carotid arteries of five anesthetized pigs in different hemodynamic situations using in-house developed non-invasive ultrasound techniques. The study protocol included administration of epinephrine, norepinephrine and β-blockade (metoprolol). No significant correlation between longitudinal displacement of the intima-media complex and wall shear rate was found. This result suggests that one or multiple pulsatile forces other than wall shear stress are also working along arteries, strongly influencing arterial wall behavior.

Accurate blood peak velocity estimation using spectral models and vector doppler

Ultrasound blood peak velocity estimates are routinely used for diagnostics, such as the grading of a stenosis. The peak velocity is typically assessed from the Doppler spectrum by locating the highest frequency detectable from noise. The selected frequency is then converted to velocity by the Doppler equation. This procedure contains several potential sources of error: the frequency selection is noise dependent and sensitive to the spectral broadening, which, in turn, is affected by the Doppler angle uncertainty. The result is, often, an inaccurate estimate. In this work we propose a new method that removes the aforementioned errors. The frequency is selected by exploiting a mathematical model of the Doppler spectrum that has recently been introduced. When a very large sample volume is used, which includes all the vessel section, the model is capable of predicting the exact threshold to be used without the need of broadening compensation. The angle ambiguity is solved by applying the threshold to the Doppler spectra measured from two different directions, according to the vector Doppler technique. The proposed approach has here been validated through Field II simulations, phantom experiments, and tests on volunteers by using defocused waves to insonify a large region from a linear array probe. A mean error lower than 1% and a mean coefficient of variability lower than 5% were measured in a variety of experimental conditions.