Coherent plane wave compounding for very high frame rate ultrasonography of rapidly moving targets

Optimum Parameters in Ultrasound Coherent Plane Wave Compounding for High LPWV Estimation: Validation on Phantom and Feasibility in 10 Subjects

OBJECTIVES

The aim of this study is to determine the optimum and fine values of the number and transmission angles of tilted plane waves for coherent plane-wave compounding (CPWC)-based high local pulse wave velocity (LPWV) estimation.

METHODS

A Verasonics system incorporating a linear array probe L14-5/38 with 128 elements and a pulsatile pump, CompuFlow1000, were used to acquire radio frequency data of 3, 5, 7, and 9 tilted plane wave sequences with angle intervals from 0° to 12° with a coarse interval increment step of 1°, and the angle intervals from 0° to 2° with a fine interval increment step of 0.25° from a carotid vessel phantom with the LPWV of 13.42 ± 0.90 m/s. The mean value, standard deviation, and coefficients of variation (CV) of the estimated LPWVs were calculated to quantitatively assess the performance of different configurations for CPWC-based LPWV estimation. Ten healthy human subjects of two age groups were recruited to assess the in vivo feasibility of the optimum parameter values.

RESULTS

The CPWC technique with three plane waves (PRF of 12 kHz corresponding to a frame rate of 4000 Hz) with an interval of 0.75° had LPWVs of 13.52 ± 0.08 m/s with the lowest CV of 1.84% on the phantom, and 5.49 ± 1.46 m/s with the lowest CV of 12.35% on 10 subjects.

CONCLUSIONS

The optimum parameters determined in this study show the best repeatability of the LPWV measurements with a vessel phantom and 10 healthy subjects, which support further studies on larger datasets for potential applications.

Ultrafast 3-D Transcutaneous Super Resolution Ultrasound Using Row-Column Array Specific Coherence-Based Beamforming and Rolling Acoustic Sub-aperture Processing: In Vitro, in Rabbit and in Human Study

OBJECTIVE

This study aimed to realise 3-D super-resolution ultrasound imaging transcutaneously with a row-column array which has far fewer independent electronic channels and a wider field of view than typical fully addressed 2-D matrix arrays. The in vivo image quality of the row-column array is generally poor, particularly when imaging non-invasively. This study aimed to develop a suite of image formation and post-processing methods to improve image quality and demonstrate the feasibility of ultrasound localisation microscopy using a row-column array, transcutaneously on a rabbit model and in a human.

METHODS

To achieve this, a processing pipeline was developed which included a new type of rolling window image reconstruction, which integrated a row-column array specific coherence-based beamforming technique with acoustic sub-aperture processing. This and other processing steps reduced the 'secondary' lobe artefacts, and noise and increased the effective frame rate, thereby enabling ultrasound localisation images to be produced.

RESULTS

Using an in vitro cross tube, it was found that the procedure reduced the percentage of 'false' locations from ∼26% to ∼15% compared to orthogonal plane wave compounding. Additionally, it was found that the noise could be reduced by ∼7 dB and the effective frame rate was increased to over 4000 fps. In vivo, ultrasound localisation microscopy was used to produce images non-invasively of a rabbit kidney and a human thyroid.

CONCLUSION

It has been demonstrated that the proposed methods using a row-column array can produce large field of view super-resolution microvascular images in vivo and in a human non-invasively.

Adaptive higher-order singular value decomposition clutter filter for ultrafast Doppler imaging of coronary flow under non-negligible tissue motion

BACKGROUND AND OBJECTIVE

With the development of advanced clutter-filtering techniques by singular value decomposition (SVD) and leveraging favorable acquisition settings such as open-chest imaging by a linear high-frequency probe and plane waves, several studies have shown the feasibility of cardiac flow measurements during the entire cardiac cycle, ranging from coronary flow to myocardial perfusion. When applying these techniques in a routine clinical setting, using transthoracic ultrasound imaging, new challenges emerge. Firstly, a smaller aperture is needed that can fit between ribs. Consequently, diverging waves are employed instead of plane waves to achieve an adequate field of view. Secondly, to ensure imaging at a larger depth, the maximum pulse repetition frequency has to be reduced. Lastly, in comparison to the open-chest scenario, tissue motion induced by the heartbeat is significantly stronger. The latter complicates substantially the distinction between clutter and blood signals.

METHODS

This study investigates a strategy to overcome these challenges by diverging wave imaging with an optimal number of tilt angles, in combination with dedicated clutter-filtering techniques. In particular, a novel, adaptive, higher-order SVD (HOSVD) clutter filter, which utilizes spatial, temporal, and angular information of the received ultrasound signals, is proposed to enhance clutter and blood separation.

RESULTS

When non-negligible tissue motion is present, using fewer tilt angles not only reduces the decorrelation between the received waveforms but also allows for collecting more temporal samples at a given ensemble duration, contributing to improved Doppler performance. The addition of a third angular dimension enables the application of HOSVD, providing greater flexibility in selecting blood separation thresholds from a 3-D tensor. This differs from the conventional threshold selection method in a 2-D spatiotemporal space using SVD. Exhaustive threshold search has shown a significant improvement in Contrast and Contrast-to-Noise ratio for Power Doppler images filtered with HOSVD compared to the SVD-based clutter filter.

CONCLUSION

With the improved settings, the obtained Power Doppler images show the feasibility of measuring coronary flow under the influence of non-negligible tissue motion in both in vitro and ex vivo.

Improved Spatiotemporal Resolution in Echocardiography Using Mixed Geometry Imaging Sequences

Cardiac ultrasound seeks to image the most dynamic environment in the body-the moving heart. Many modern ultrasound imaging techniques address the tradeoff between spatial and temporal resolution using either narrow focused beams or with broad beam, synthetic aperture (SA) sequences that have been shown to suffer from motion artifacts. Retrospective encoding for conventional ultrasound sequences (REFoCUS) unifies the processing of these various geometric sequences, but the motion sensitivity of this approach has yet to be investigated. We hypothesize that a "mixed sequence" enabled by the REFoCUS method incorporating several beam geometries may better resolve cardiac motion over a wide field of view (FOV) and at a high frame rate. First, the motion sensitivity of REFoCUS was evaluated in simulation for several focused and broad transmit profiles. Focused transmissions resolve both lateral and axial motion much more effectively than broad transmissions, with performance similar to conventional beamforming techniques. Second, a mixed sequence was designed that insonifies the full field-of-view with plane wave (PW) transmissions and key moving targets with focused transmissions. This mixed sequence was tested in simulation and in vivo and was used to image the heart as well as the liver, a low-motion control. By combining a sparse PW sequence ( n=60 ) with a small group of targeted focused transmissions ( n = 10), the anterior mitral valve leaflet (AML) at its peak observed velocity was better resolved. We believe that mixed sequences have strong potential to resolve cardiac motion at clinically relevant frame rates.

Windowed Radon Transform for Robust Speed-of-Sound Imaging With Pulse-Echo Ultrasound

In recent years, methods estimating the spatial distribution of tissue speed of sound with pulse-echo ultrasound are gaining considerable traction. They can address limitations of B-mode imaging, for instance in diagnosing fatty liver diseases. Current state-of-the-art methods relate the tissue speed of sound to local echo shifts computed between images that are beamformed using restricted transmit and receive apertures. However, the aperture limitation affects the robustness of phase-shift estimations and, consequently, the accuracy of reconstructed speed-of-sound maps. Here, we propose a method based on the Radon transform of image patches able to estimate local phase shifts from full-aperture images. We validate our technique on simulated, phantom and in-vivo data acquired on a liver and compare it with a state-of-the-art method. We show that the proposed method enhances the stability to changes of beamforming speed of sound and to a reduction of the number of insonifications. In particular, the deployment of pulse-echo speed-of-sound estimation methods onto portable ultrasound devices can be eased by the reduction of the number of insonifications allowed by the proposed method.

On the Role of Coherent Plane Wave Compounding in Shear Wave Elasticity Imaging: The Convolution Effect and Its Implications

OBJECTIVE

The clinical applicability of shear wave elasticity imaging (SWEI) has been confounded by its appreciable inter-system variability and unsatisfactory sensitivity. While SWEI relies on plane wave imaging (PWI) to achieve real-time rendering, it has been rarely noticed that PWI can affect SWEI's performance. This work is aimed at demonstrating that the use of coherent plane wave compounding (CPWC) can be a factor causing SWEI's underperformance.

METHODS

We presented a model to formally describe the slow-time behavior of CPWC in motion tracking. This model reveals that CPWC introduces temporal convolution on the observed motion, making the motion sampling process a low-pass filter (LPF). For validation, shear waves were produced in a phantom in the same way but sampled via PWI using different compounding numbers (CN) and pulse repetition frequencies (PRF), with the obtained signals compared with the inferences drawn from our model. Similar experiments were performed to reconstruct two small targets in the phantom in order to appraise the impact of CPWC on SWEI's sensitivity.

DISCUSSION

The validation experiment shows that the measurements match well with the model inferences, which verifies the LPF nature of CPWC. The phantom study also shows that either increasing CN or decreasing PRF can cause the loss of high-frequency motion information, leading to blurred target delineation by SWEI.

CONCLUSION

The convolution effect can help understand the variability of SWEI. Researchers should beware this effect when working on SWEI standardization. Clinicians using SWEI should also be cautious because this effect makes it harder to identify small lesions.

Simultaneous High-Frame-Rate Acoustic Plane-Wave and Optical Imaging of Intracranial Cavitation in Polyacrylamide Brain Phantoms during Blunt Force Impact

Blunt and blast impacts occur in civilian and military personnel, resulting in traumatic brain injuries necessitating a complete understanding of damage mechanisms and protective equipment design. However, the inability to monitor in vivo brain deformation and potential harmful cavitation events during collisions limits the investigation of injury mechanisms. To study the cavitation potential, we developed a full-scale human head phantom with features that allow a direct optical and acoustic observation at high frame rates during blunt impacts. The phantom consists of a transparent polyacrylamide material sealed with fluid in a 3D-printed skull where windows are integrated for data acquisition. The model has similar mechanical properties to brain tissue and includes simplified yet key anatomical features. Optical imaging indicated reproducible cavitation events above a threshold impact energy and localized cavitation to the fluid of the central sulcus, which appeared as high-intensity regions in acoustic images. An acoustic spectral analysis detected cavitation as harmonic and broadband signals that were mapped onto a reconstructed acoustic frame. Small bubbles trapped during phantom fabrication resulted in cavitation artifacts, which remain the largest challenge of the study. Ultimately, acoustic imaging demonstrated the potential to be a stand-alone tool, allowing observations at depth, where optical techniques are limited.

Optimization of the Tracking Beam Sequence in Harmonic Motion Imaging

Harmonic motion imaging (HMI) is an ultrasound elastography technique that estimates the viscoelastic properties of tissues by inducing localized oscillatory motion using focused ultrasound (FUS). The resulting displacement, assumed to be inversely proportional to the tissue local stiffness, is estimated using an imaging array based on RF speckle tracking. In conventional HMI, this is accomplished with plane-wave (PW) imaging, which inherently suffers from low lateral resolution. Coherent PW compounding (PWC) leverages spatial and temporal resolution using synthetic focusing in transmit. In this study, we introduced focused imaging with parallel tracking in HMI and compared parallel tracking of various transmit F-numbers (F/2.6, 3, 4, and 5) qualitatively and quantitatively with PW and PWC imaging at various compounded angle ranges (6°, 12°, and 18°). An in silico model of a 56-kPa spherical inclusion (diameter: 3.6 mm) embedded in a 5.3-kPa background and a 5.3-kPa elastic phantom with cylindrical inclusions (Young's moduli: 22-56 kPa, diameters: 2.0-8.6 mm) were imaged to assess different tracking beam sequences. Speckle biasing in displacement estimation associated with parallel tracking was also investigated and concluded to be negligible in HMI. Parallel tracking in receive (Rx) resulted in 2%-7% and 8%-12% increase compared to PW imaging ( ) in HMI contrast and contrast-to-noise ratio in silico and phantoms. Focused imaging with parallel tracking in Rx was concluded to be most robust among PW and PWC imaging for displacement estimation, and its preclinical feasibility was demonstrated in postsurgical human cancerous breast tissue specimens and in vivo murine models of breast cancer.

The Effect of Spatial Velocity Gradients on Block-Matching Accuracy for Ultrasound Velocimetry

OBJECTIVE

Block matching serves as the foundation for ultrasound velocimetry techniques such as blood speckle tracking and echo-particle image velocimetry. Any spatial velocity gradients (SVGs) inside a block-matching pair will result in tracking error, due to both the finite block size and the ultrasound point-spread-function. We assess, using an in silico sinusoidal flow phantom, the effect of SVG magnitude and beam-to-flow angle on block-matching bias and precision. Secondarily we assess the effect that SVGs have on velocimetry bias when using angled plane-wave compounding.

METHODS

The magnitude and angle of SVGs were varied by adjusting the wavelength and direction of a sinusoidal flow profile. Scatterers displaced by this flow profile were used for simulating ultrasound radio frequency data at discrete time points. After beamforming, the 2-D flow field was estimated using block matching. Two imaging sequences were tested, a single plane-wave and a three-angled plane-wave.

RESULTS

Smaller sinusoidal flow wavelengths resulted in increased bias and reduced precision, revealing an inverse relationship between sinusoidal flow wavelength and tracking error, with median errors ranging from 69%-90% for the smallest flow wavelengths (highest SVGs) down to 3%-5% for the largest (lowest SVGs). The SVG angle was also important, in which lateral SVGs (with axially oriented flows) resulted in significant speckle decorrelation and high tracking errors in regions with high SVGs. Conversely, axial SVGs (laterally oriented flow) experienced higher bias in the peak velocity regions of the flow profile. Coherent compounding resulted in higher velocity errors than using a single transmission for lateral SVGs but not for axial SVGs.

CONCLUSION

The highest SVGs that could be measured with ≤10% error was when the sinusoidal flow wavelength was less than 20 times the ultrasound pulse wavelength. The clinical significance is that the high SVGs present in high kinetic energy flows, such as severe carotid stenosis and aortic regurgitation, will limit the ability to accurately quantify the velocities in these flow structures.

A KL Divergence-Based Loss for In Vivo Ultrafast Ultrasound Image Enhancement with Deep Learning

Ultrafast ultrasound imaging, characterized by high frame rates, generates low-quality images. Convolutional neural networks (CNNs) have demonstrated great potential to enhance image quality without compromising the frame rate. However, CNNs have been mostly trained on simulated or phantom images, leading to suboptimal performance on in vivo images. In this study, we present a method to enhance the quality of single plane wave (PW) acquisitions using a CNN trained on in vivo images. Our contribution is twofold. Firstly, we introduce a training loss function that accounts for the high dynamic range of the radio frequency data and uses the Kullback-Leibler divergence to preserve the probability distributions of the echogenicity values. Secondly, we conduct an extensive performance analysis on a large new in vivo dataset of 20,000 images, comparing the predicted images to the target images resulting from the coherent compounding of 87 PWs. Applying a volunteer-based dataset split, the peak signal-to-noise ratio and structural similarity index measure increase, respectively, from 16.466 ± 0.801 dB and 0.105 ± 0.060, calculated between the single PW and target images, to 20.292 ± 0.307 dB and 0.272 ± 0.040, between predicted and target images. Our results demonstrate significant improvements in image quality, effectively reducing artifacts.

Automatic Probe Localization in Multiaperture Ultrasound by Radon Domain Signal Matching

In multiaperture ultrasound, several ultrasound probes with different insonification angles are combined to increase the field of view and angular coverage of image structures. A full reconstruction incorporating all possible combinations of transmitting and receiving probes has been shown to improve resolution, contrast, and angular coverage beyond what can be achieved by the registration of single images from different probes. A major challenge in multiaperture imaging is the correct determination of relative probe locations. A registration based on the content of images from different probes is challenging due to the decorrelation of image structures and speckle with increasing angle between the probes. We propose a probe localization method for plane-wave ultrasound that uses solely the receive dataset of a nontransmitting probe. The localization is performed by signal tracking in the Radon domain. To demonstrate that the method does not rely on common structures in the individual images, we show that a satisfying localization can be performed in pure speckle for angles, where the speckle patterns have completely decorrelated. The method shows potential for real-time probe localization in free-hand multiprobe ultrasound imaging or for flexible and wearable multiarray combination of multiple capacitive micromachined (CMUT)-based systems in the future.

Optimized Virtual Sources Distributions for 3-D Ultrafast Diverging Wave Compounding Imaging: A Simulation Study

Ultrafast ultrasound imaging allows observing rapid phenomena; combined with 3-D imaging it has the potential to provide a more accurate analysis of organs which leads, in the end, to better diagnosis. Coherent compounding using diverging waves is commonly used to reconstruct high-quality images on large volumes while keeping the frame rate high enough to allow dynamic analysis. In practice, the virtual sources (VSs) that drive the diverging waves are often distributed in a deterministic way: following a regular grid, concentric rings, and spirals. Even though those deterministic distributions can offer various tradeoffs in terms of imaging performance, other distributions can be considered to improve imaging performance. It is herein suggested to look at alternative VSs distributions for optimizing the lateral resolution and the secondary lobes level (SLL) on several point spread functions (PSFs) by means of a multiobjective genetic algorithm. The optimization framework has led to seven pseudo-irregular distributions of VSs distributions that have not yet been found in the literature. An analysis of the imaging performance with a simulated phantom shows that these new distributions offer different tradeoffs between lateral resolution and contrast, respectively, measured on point-like reflectors and anechoic cysts. As an example, one of these optimized distributions improves the lateral resolution by 16% and gives equivalent contrast values on cysts and PSF isotropy properties, when compared to a concentric-rings-based distribution.

Frequency-Dependent F-Number Suppresses Grating Lobes and Improves the Lateral Resolution in Coherent Plane-Wave Compounding

Ultrafast imaging modes, such as coherent plane-wave compounding (CPWC), increase image uniformity and reduce grating lobe artifacts by dynamic receive apertures. The focal length and the desired aperture width maintain a given ratio, which is called the F -number. Fixed F -numbers, however, exclude useful low-frequency components from the focusing and reduce the lateral resolution. Herein, this reduction is avoided by a frequency-dependent F -number. This F -number derives from the far-field directivity pattern of a focused aperture and can be expressed in closed form. The F -number, at low frequencies, widens the aperture to improve the lateral resolution. The F -number, at high frequencies, narrows the aperture to avoid lobe overlaps and suppress grating lobes. Phantom and in vivo experiments with a Fourier-domain beamforming algorithm validated the proposed F -number in CPWC. The lateral resolution, which was measured by the median lateral full-widths at half-maximum of wires, improved by up to 46.8% and 14.9% in a wire and a tissue phantom, respectively, in comparison to fixed F -numbers. Grating lobe artifacts, which were measured by the median peak signal-to-noise ratios of wires, reduced by up to 9.9 dB in comparison to the full aperture. The proposed F -number thus outperformed F -numbers that were recently derived from the directivity of the array elements.

ARU-GAN: U-shaped GAN based on Attention and Residual connection for super-resolution reconstruction

Plane-wave ultrasound imaging technology offers high-speed imaging but lacks image quality. To improve the image spatial resolution, beam synthesis methods are used, which often compromise the temporal resolution. Herein, we propose ARU-GAN, a super-resolution reconstruction model based on residual connectivity and attention mechanisms, to address this issue. ARU-GAN comprises a Full-scale Skip-connection U-shaped Generator (FSUG) with an attention mechanism and a Residual Attention Patch Discriminator (RAPD). The former captures global and local features of the image by using full-scale skip-connections and attention mechanisms. The latter focuses on changes in the image at different scales to enhance its discriminative ability at the patch level. ARU-GAN was trained using a combined loss function on the Plane-Wave Imaging Challenge in Medical Ultrasound (PICMUS) 2016 dataset, which includes three types of targets: point targets, cyst targets, and in-vivo targets. Compared to Coherent Plane-Wave Compounding (CPWC), ARU-GAN achieved a reduction in Full Width at Half Maximum (FWHM) by 5.78%-20.30% on point targets, improved Contrast (CR) by 7.59-11.29 percentage points, and Contrast to Noise Ratio (CNR) by 30.58%-45.22% on cyst targets. On in-vivo target, ARU-GAN improved the Peak Signal-to-Noise Ratio (PSNR) by 11.94%, the Complex-Wavelet Structural Similarity Index Measurement (CW-SSIM) by 17.11%, and the Normalized Cross Correlation (NCC) by at least 2.17% compared to existing deep learning methods. In conclusion, ARU-GAN is a promising model for the super-resolution reconstruction of plane-wave medical ultrasound images. It provides a novel solution for improving image quality, which is essential for clinical practice.

An improved Fourier Ptychography algorithm for ultrasonic array imaging

Inspired by the optical imaging algorithm, the Fourier Ptychography (FP) algorithm is adopted to improve the resolution of ultrasonic array imaging. In the FP algorithm, the steady-state spectrum is utilized to recover the high-resolution ultrasonic images. Meanwhile, the parameters of FP algorithm are empirical, which can affect the imaging quality of ultrasonic array. Then the particle swarm optimization (PSO) algorithm is used to optimize the parameters of FP algorithm to further improve the imaging quality of ultrasonic array. The tungsten imaging experiments and pig eye imaging experiments are conducted to demonstrate the feasibility and effectiveness of the developed algorithm. In addition, the proposed algorithm and the coherent wave superposition (CWS) algorithm are both based on single plane wave (SPW) algorithms and they are then compared. The results show that the CWS algorithm and FP algorithm have good longitudinal and lateral resolutions, respectively. The particle swarm optimization-based FP (PSOFP) imaging algorithm has both excellent lateral and longitudinal resolutions. The average lateral resolution of PSOFP imaging algorithm is improved by 34.47% compared with CWS imaging algorithm in the tungsten wires experiments, and the lateral boundary structure width of the lens is improved by 49.48% in the pig eye experiments. The proposed algorithm can effectively improve the ultrasonic imaging quality for medical application.

Doppler and Pair-Wise Optical Flow Constrained 3D Motion Compensation for 3D Ultrasound Imaging

Volumetric (3D) ultrasound imaging using a 2D matrix array probe is increasingly developed for various clinical procedures. However, 3D ultrasound imaging suffers from motion artifacts due to tissue motions and a relatively low frame rate. Current Doppler-based motion compensation (MoCo) methods only allow 1D compensation in the in-range dimension. In this work, we propose a new 3D-MoCo framework that combines 3D velocity field estimation and a two-step compensation strategy for 3D diverging wave compounding imaging. Specifically, our framework explores two constraints of a round-trip scan sequence of 3D diverging waves, i.e., Doppler and pair-wise optical flow, to formulate the estimation of the 3D velocity fields as a global optimization problem, which is further regularized by the divergence-free and first-order smoothness. The two-step compensation strategy is to first compensate for the 1D displacements in the in-range dimension and then the 2D displacements in the two mutually orthogonal cross-range dimensions. Systematical in-silico experiments were conducted to validate the effectiveness of our proposed 3D-MoCo method. The results demonstrate that our 3D-MoCo method achieves higher image contrast, higher structural similarity, and better speckle patterns than the corresponding 1D-MoCo method. Particularly, the 2D cross-range compensation is effective for fully recovering image quality.

Smart Sensors and Microtechnologies in the Precision Medicine Approach against Lung Cancer

BACKGROUND AND RATIONALE

The therapeutic interventions against lung cancer are currently based on a fully personalized approach to the disease with considerable improvement of patients' outcome. Alongside continuous scientific progresses and research investments, massive technologic efforts, innovative challenges, and consolidated achievements together with research investments are at the bases of the engineering and manufacturing revolution that allows a significant gain in clinical setting.

AIM AND METHODS

The scope of this review is thus to focus, rather than on the biologic traits, on the analysis of the precision sensors and novel generation materials, as semiconductors, which are below the clinical development of personalized diagnosis and treatment. In this perspective, a careful revision and analysis of the state of the art of the literature and experimental knowledge is presented.

RESULTS

Novel materials are being used in the development of personalized diagnosis and treatment for lung cancer. Among them, semiconductors are used to analyze volatile cancer compounds and allow early disease diagnosis. Moreover, they can be used to generate MEMS which have found an application in advanced imaging techniques as well as in drug delivery devices.

CONCLUSIONS

Overall, these issues represent critical issues only partially known and generally underestimated by the clinical community. These novel micro-technology-based biosensing devices, based on the use of molecules at atomic concentrations, are crucial for clinical innovation since they have allowed the recent significant advances in cancer biology deciphering as well as in disease detection and therapy. There is an urgent need to create a stronger dialogue between technologists, basic researchers, and clinicians to address all scientific and manufacturing efforts towards a real improvement in patients' outcome. Here, great attention is focused on their application against lung cancer, from their exploitations in translational research to their application in diagnosis and treatment development, to ensure early diagnosis and better clinical outcomes.

Improving lateral resolution and contrast by combining coherent plane-wave compounding with adaptive weighting for medical ultrasound imaging

Due to the severe lateral lobe artifact by coherent plane-wave compounding (CPWC) and the low signal-to-noise ratio of radiofrequency (RF) data collected from the plane wave, the adaptive beamforming methods based on focused wave imaging (FWI) are improper to be directly applied to CPWC. To obtain a high-quality image with high resolution and contrast, this study combined the threshold phase coherence factor (THR-PCF) with the reconstructed covariance matrix minimum variance (RCM-MV) and then proposed a novel CPWC-based adaptive beamforming algorithm, THR-PCF + RCM-MV. The simulation, phantom, and in-vivo experiments were performed to investigate the performance of the proposed methods in comparison with the CPWC and the classical adaptive methods including the minimum variance (MV), generalized coherence factor (GCF) and their combination GCF + MV. The simulation results demonstrated that the THR-PCF + RCM-MV beamformer improved contrast ratio (CR) by 28.14%, contrast noise ratio (CNR) by 22.01%, speckle signal-to-noise ratio (s_SNR) by 23.58%, generalized contrast-to-noise ratio (GCNR) by 0.3%, and the full width at half maximum (FWHM) by 43.38% on average, compared with the GCF + MV method. The phantom experimental results showed a better performance of the THR-PCF + RCM-MV beamformer with an average improvement by 21.95% in CR, 2.62% in s_SNR, and 48.64% in FWHM compared with the GCF + MV. Meanwhile, the results showed that the image quality of the near and far fields was enhanced by the THR-PCF + RCM-MV. The in-vivo imaging results showed that our new method had potential for clinical application. In conclusion, the lateral resolution and contrast of medical ultrasound imaging could be improved greatly with our proposed method.

A study of wide unfocused wavefront for convex-array ultrasound imaging

Ultrafast ultrasound imaging modalities have attracted a lot of attention in the ultrasound community. It breaks the compromise between the frame rate and the region of interest by insonifying the whole medium with wide unfocused waves. Coherent compounding can be performed to enhance the image quality at a cost of frame rate. Ultrafast imaging has wide clinical applications, such as vector Doppler imaging and shear elastography. On the other hand, the use of unfocused waves is still marginal with convex-array transducers. For convex array, plane wave imaging is limited by the complicated transmission delay calculation, limited field-of-view, and inefficient coherent compounding. In this article, we study three wide unfocused wavefronts, namely, lateral virtual-source defined diverging wave imaging (latDWI), tilt virtual-source defined diverging wave imaging (tiltDWI), and Archimedean-spiral-based imaging (AMI) for convex-array imaging using the full-aperture transmission. The analytical monochromatic wave solutions to this three imaging are given. The mainlobe width and grating lobe position are given explicitly. Theoretical -6 dB beamwidth and synthetic transmit field response are studied. Simulation studies are carried on with the point targets and hypoechoic cysts. Time-of-flight formulas are given explicitly for beamforming. The conclusions are in good agreement with the theory: latDWI provides the finest lateral resolution but generates the severest axial lobe level for scatterers with large obliquities (i.e., for scatterers located at the image border) which degrades the image contrast. This effect gets worsen as the compound number increases. The tiltDWI and AMI give a very close performance on resolution and image contrast. AMI displays better contrast with a small compound number.

Point detection in textured ultrasound images

Detection of point scatterers in textured ultrasound images can be challenging. This paper investigates how four multilook methods can improve the detection. We analyze many images with known point scatterer locations and randomly textured backgrounds. The normalized matched filter (NMF) and multilook coherence factor (MLCF) methods are normalized methods that do not require any texture correction prior to detection analysis. They are especially propitious when optimal texture correction of the ultrasound images is difficult to obtain. The results show significant improvement in detection performance when the MLCF method is weighted with the prewhitened and texture corrected image. The method can be applied even when we do not have prior knowledge about the optimal prewhitening limits. The multilook methods NMF and NMF weighted (NMFW) are very favorable methods to apply on images where acoustic noise dominates the speckle background.

High Frame Rate Ultrasound Imaging by Means of Tensor Completion: Application to Echocardiography

High frame rate ultrasound (US) imaging enables the monitoring of fast-moving organs. In echocardiography, this is especially needed due to the existence of rapidly moving structures, such as the heart valves. In the last two decades, various methods have been proposed to improve the frame rate. Here, we propose a novel method, based on binary coding patterns (BCPs) and tensor completion (TC), to increase the temporal resolution (i.e., frame rate) in the preprocessing stage of conventional focused ultrasound imaging (CFUI). The rationale behind our proposal is to perform, at first, the beamforming of a fraction of the scan lines, randomly selected in each frame based on BCP. Then, we reconstruct the missing scan lines through TC. The latter is an effective technique for recovering missing information from a low-rank tensor, based on a small number of observations using rank minimization. Following our approach, reducing the transmissions events needed to generate an image, the frame rate is increased by the same proportion. We have applied the proposed technique to a pre-beamformed radio frequency (RF) echocardiographic dataset. Our results show that we can improve the frame rate by a factor from 3 to 4, while keeping the structural similarity (SSIM) of the reconstructed tensor and the original one at values higher than 0.98.

High frequency ultrasound vibrational shear wave elastography for preclinical research

Preclinical evaluation of novel therapies using models of cancer is an important tool in cancer research, where imaging can provide non-invasive tools to characterise the internal structure and function of tumours. The short propagation paths when imaging tumours and organs in small animals allow the use of high frequencies for both ultrasound and shear waves, providing the opportunity for high-resolution shear wave elastography and hence its use for studying the heterogeneity of tissue elasticity, where heterogeneity may be a predictor of tissue response. Here we demonstrate vibrational shear wave elastography (VSWE) using a mechanical actuator to produce high frequency (up to 1000 Hz) shear waves in preclinical tumours, an alternative to the majority of preclinical ultrasound SWE studies where an acoustic radiation force impulse is required to create a relatively low-frequency broad-band shear-wave pulse. We implement VSWE with a high frequency (17.8 MHz) probe running a focused line-by-line ultrasound imaging sequence which as expected was found to offer improved detection of 1000 Hz shear waves over an ultrafast planar wave imaging sequence in a homogenous tissue-mimicking phantom. We test the VSWE in anex vivotumour xenograft, demonstrating the ability to detect shear waves up to 10 mm from the contactor position at 1000 Hz. By reducing the kernel size used for shear wave speed estimation to 1 mm we are able to produce shear wave speed images with spatial resolution of this order. Finally, we present VSWE data from xenograft tumoursin vivo, demonstrating the feasibility of the technique in mice under isoflurane sedation. Mean shear wave speeds in the tumours are in good agreements with those reported by previous authors. Characterising the frequency dependence of shear wave speed demonstrates the potential to quantify the viscoelastic properties of tumoursin vivo.

Multi Spherical Wave Imaging Method Based on Ultrasonic Array

The imaging range of traditional plane wave imaging is usually limited by the directivity of the plane wave. In this paper, a multi spherical wave imaging method based on an ultrasonic array is proposed, which radiates both compression and shear waves in a solid medium to form the multi spherical wave. Firstly, excitation characteristics of the multi spherical wave are analyzed theoretically and the calculation method of echo delay of multi spherical wave imaging is derived. Multi spherical wave imaging is compared with conventional ultrasonic plane wave imaging by designing experiments. Compared with ultrasonic plane wave imaging, multi spherical wave imaging is not limited to the size of the transducer and can greatly improve the detection range. In addition, compared with the multi plane wave imaging method, the multi spherical wave imaging algorithm is relatively simple, fast, and has high application value.

Improving plane wave ultrasound imaging through real-time beamformation across multiple arrays

Ultrasound imaging is a widely used diagnostic tool but has limitations in the imaging of deep lesions or obese patients where the large depth to aperture size ratio (f-number) reduces image quality. Reducing the f-number can improve image quality, and in this work, we combined three commercial arrays to create a large imaging aperture of 100 mm and 384 elements. To maintain the frame rate given the large number of elements, plane wave imaging was implemented with all three arrays transmitting a coherent wavefront. On wire targets at a depth of 100 mm, the lateral resolution is significantly improved; the lateral resolution was 1.27 mm with one array (1/3 of the aperture) and 0.37 mm with the full aperture. After creating virtual receiving elements to fill the inter-array gaps, an autoregressive filter reduced the grating lobes originating from the inter-array gaps by − 5.2 dB. On a calibrated commercial phantom, the extended field-of-view and improved spatial resolution were verified. The large aperture facilitates aberration correction using a singular value decomposition-based beamformer. Finally, after approval of the Stanford Institutional Review Board, the three-array configuration was applied in imaging the liver of a volunteer, validating the potential for enhanced resolution.

Independent Component Analysis Filter for Small Vessel Contrast Imaging During Fast Tissue Motion

Suppressing tissue clutter is an essential step in blood flow estimation and visualization, even when using ultrasound contrast agents. Blind source separation (BSS)-based clutter filter for high-framerate ultrasound imaging has been reported to perform better in tissue clutter suppression than the conventional frequency-based wall filter and nonlinear contrast pulsing schemes. The most notable BSS technique, singular value decomposition (SVD) has shown compelling results in cases of slow tissue motion. However, its performance degrades when the tissue motion is faster than the blood flow speed, conditions that are likely to occur when imaging the small vessels, such as in the myocardium. Independent component analysis (ICA) is another BSS technique that has been implemented as a clutter filter in the spatiotemporal domain. Instead, we propose to implement ICA in the spatial domain where motion should have less impact. In this work, we propose a clutter filter with the combination of SVD and ICA to improve the contrast-to-background ratio (CBR) in cases where tissue velocity is significantly faster than the flow speed. In an in vitro study, the range of fast tissue motion velocity was 5-25 mm/s and the range of flow speed was 1-12 mm/s. Our results show that the combination of ICA and SVD yields 7-10 dB higher CBR than SVD alone, especially in the tissue high-velocity range. The improvement is crucial for cardiac imaging where relatively fast myocardial motions are expected.

Point Detection in Ultrasound Using Prewhitening and Multilook Optimization

We investigate methods to improve the detection of point scatterers in ultrasound imaging using the standard delay-and-sum (DAS) image as our starting point. An optimized whitening transform can increase the spatial resolution of the image. By splitting an image's frequency spectrum into many subsets using the multilook technique, we can exploit the coherent properties of a point scatterer. We present three new multilook methods and evaluate their effect on point detection. The performances are compared to DAS using synthetic aperture Field II simulations of a point scatterer in uniform speckle background. The results show that optimized prewhitening of the images can significantly improve the point detection. The multilook methods have the potential to improve the detection performance further when a sufficient number of looks are used. If prior knowledge about the optimal spectrum limits is unavailable and a nonoptimal prewhitening is applied, applying that the new multilook methods can considerably improve the point detection.

CNN-Based Image Reconstruction Method for Ultrafast Ultrasound Imaging

Ultrafast ultrasound (US) revolutionized biomedical imaging with its capability of acquiring full-view frames at over 1 kHz, unlocking breakthrough modalities such as shear-wave elastography and functional US neuroimaging. Yet, it suffers from strong diffraction artifacts, mainly caused by grating lobes, sidelobes, or edge waves. Multiple acquisitions are typically required to obtain a sufficient image quality, at the cost of a reduced frame rate. To answer the increasing demand for high-quality imaging from single unfocused acquisitions, we propose a two-step convolutional neural network (CNN)-based image reconstruction method, compatible with real-time imaging. A low-quality estimate is obtained by means of a backprojection-based operation, akin to conventional delay-and-sum beamforming, from which a high-quality image is restored using a residual CNN with multiscale and multichannel filtering properties, trained specifically to remove the diffraction artifacts inherent to ultrafast US imaging. To account for both the high dynamic range and the oscillating properties of radio frequency US images, we introduce the mean signed logarithmic absolute error (MSLAE) as a training loss function. Experiments were conducted with a linear transducer array, in single plane-wave (PW) imaging. Trainings were performed on a simulated dataset, crafted to contain a wide diversity of structures and echogenicities. Extensive numerical evaluations demonstrate that the proposed approach can reconstruct images from single PWs with a quality similar to that of gold-standard synthetic aperture imaging, on a dynamic range in excess of 60 dB. In vitro and in vivo experiments show that trainings carried out on simulated data perform well in experimental settings.

Experimental Study of Aperiodic Plane Wave Imaging for Ultrafast 3-D Ultrasound Imaging

OBJECTIVE

Although plane wave imaging (PWI) with multiple plane waves (PWs) steered at different angles enables ultrafast three-dimensional (3-D) ultrasonic imaging, there is still a challenging tradeoff between image quality and frame rate. To address this challenge, we recently proposed the aperiodic PWI (APWI) with mathematical analysis and simulation study. In this paper, we demonstrate the feasibility of APWI and evaluate the performance with phantom and in vivo experiments.

METHODS

APWI with a concentric ring angle pattern (APWI-C) and APWI with a sunflower pattern (APWI-S) are evaluated. For experimental verification of the methods, the experimental results are compared with simulation results in terms of the mainlobe-to-sidelobe ratio. In addition, the performance of APWI is compared with that of conventional PWI by using a commercial phantom. To examine the potential for clinical use of APWI, a gallstone-mimicking phantom study and an in vivo carotid artery experiment are also conducted.

RESULTS

In the phantom study, the APWI methods provide a contrast ratio approximately 23 dB higher than that of PWI. In a gallstone mimicking experiment, the proposed methods yield 3-D rendered stone images more similar to the real stones than PWI. In the in vivo carotid artery images, APWI reduces the clutter artifacts inside the artery.

CONCLUSION

Phantom and in vivo studies show that the APWI enhances the contrast without compromising the spatial resolution and frame rate.

SIGNIFICANCE

This study experimentally demonstrates the feasibility and advantage of APWI for ultrafast 3-D ultrasonic imaging.

XDoppler: Cross-Correlation of Orthogonal Apertures for 3D Blood Flow Imaging

Row column addressing (RCA) transducers have the potential to provide volumetric imaging at ultrafast frame rate using a low channel count over a large field of view. In previous works we have shown that vascular imaging of large arteries as well as functional neuroimaging of the rat brain were feasible using RCA orthogonal plane wave imaging (OPW), but these applications required to transmit many plane waves, significantly reducing the frame rate. In this study, we introduce XDoppler imaging, a novel method to increase the performances of RCA flow imaging by taking advantage of the blood spatial decorrelation statistics combined with the limited spatial overlap of the point spread functions (PSF) of the two orthogonal apertures of the RCA transducer. We provide at first a theoretical basis to understand how the correlation operation reduces the sidelobe level. Then, we demonstrate both in vitro and in vivo in the human carotid artery and in the rat brain that XDoppler provides a significant gain in contrast-to-noise ratio (CNR) (between 3 and 6 dB depending on the application) compared to OPW. This improvement also leads to a sensitivity increase in the rat brain as more blood vessels are detected by XDoppler imaging.

Use of Intertransmission Coherence for Haze Artifact Suppression in Cardiovascular Synthetic Aperture Ultrasound Imaging

Synthetic aperture (SA) beamforming is a principal technology of modern medical ultrasound imaging. In that the use of focused transmission provides superior signal-to-noise ratio (SNR) and is promising for cardiovascular diagnosis at the maximum imaging depth of about 160 mm. But there is a pitfall in increasing the frame rate to more than 80 frames per second (frames/s) without image degradation by the haze artifact produced when the transmit foci (SA virtual sources) placed within the imaging field. We hypothesize that the source of this artifact is a grating lobe caused by coarse (decimated) multiple transmission and manifesting in the low brightness region in the accelerated-frame-rate images. We propose an intertransmission coherence factor (ITCF) method suppressing haze artifacts caused by coarse-pitch multiple transmission. The method is expected to suppress the image blurring because the SA grating lobe signal is less coherent than the main lobe signals. We evaluated an ITCF algorithm for suppressing the grating artifact when the transmission pitch is up to four times larger than the normal pitch (equivalent to 160 frames/s). In in-vitro and in-vivo experiments, we demonstrated the strong relation of haze artifact with the grating lobe due to the coarse-pitch transmission. Then, we confirmed that the ITCF method suppresses the haze artifact of a human heart by 15 dB while preserving the spatial resolution. The ITCF method combined with focused transmission SA beamforming is a valid method for getting cardiovascular ultrasound B-mode images without making a compromise in the trade-off relationship between the frame rate and the SNR.

Simultaneous Noise Suppression and Incoherent Artifact Reduction in Ultrafast Ultrasound Vascular Imaging

Ultrasound vascular imaging based on ultrafast plane wave imaging and singular value decomposition (SVD) clutter filtering has demonstrated superior sensitivity in blood flow detection. However, ultrafast ultrasound vascular imaging is susceptible to electronic noise due to the weak penetration of unfocused waves, leading to a lower signal-to-noise ratio (SNR) at larger depths. In addition, incoherent clutter artifacts originating from strong and moving tissue scatterers that cannot be completely removed create a strong mask on top of the blood signal that obscures the vessels. Herein, a method that can simultaneously suppress the background noise and incoherent artifacts is proposed. The method divides the tilted plane or diverging waves into two subgroups. Coherent spatial compounding is performed within each subgroup, resulting in two compounded data sets. An SVD-based clutter filter is applied to each data set, followed by a correlation between the two data sets to produce a vascular image. Uncorrelated noise and incoherent artifacts can be effectively suppressed with the correlation process, while the coherent blood signal can be preserved. The method was evaluated in wire-target simulations and phantom, in which around 7-10-dB SNR improvement was shown. Consistent results were found in a flow channel phantom with improved SNR by the proposed method (39.9 ± 0.2 dB) against conventional power Doppler (29.1 ± 0.6 dB). Last, we demonstrated the effectiveness of the method combined with block-wise SVD clutter filtering in a human liver, breast tumor, and inflammatory bowel disease data sets. The improved blood flow visualization may facilitate more reliable small vessel imaging for a wide range of clinical applications, such as cancer and inflammatory diseases.

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.

Quantitative assessment of ensemble coherency in contrast-free ultrasound microvasculature imaging

Purpose

Contrast‐free visualization of microvascular blood flow (MBF) using ultrasound can play a valuable role in diagnosis and detection of diseases. In this study, we demonstrate the importance of quantifying ensemble coherence for robust MBF imaging. We propose a novel approach to quantify ensemble coherence by estimating the local spatiotemporal correlation (LSTC) image, and evaluate its efficacy through simulation and in vivo studies.

Methods

The in vivo patient studies included three volunteers with a suspicious breast tumor, 15 volunteers with a suspicious thyroid tumor, and two healthy volunteers for renal MBF imaging. The breast data displayed negligible prior motion and were used for simulation analysis involving synthetically induced motion, to assess its impact on ensemble coherency and motion artifacts in MBF images. The in vivo thyroid data involved complex physiological motion due to its proximity to the pulsating carotid artery, which was used to assess the in vivo efficacy of the proposed technique. Further, in vivo renal MBF images demonstrated the feasibility of using the proposed ensemble coherence metric for curved array‐based MBF imaging involving phase conversion. All ultrasound data were acquired at high imaging frame rates and the tissue signal was suppressed using spatiotemporal clutter filtering. Thyroid tissue motion was estimated using two‐dimensional normalized cross correlation‐based speckle tracking, which was subsequently used for ensemble motion correction. The coherence of the MBF image was quantified based on Casorati correlation of the Doppler ensemble.

Results

The simulation results demonstrated that an increase in ensemble motion corresponded with a decrease in ensemble coherency, which reciprocally degraded the MBF images. Further the data acquired from breast tumors demonstrated higher ensemble coherency than that from thyroid tumors. Motion correction improved the coherence of the thyroid MBF images, which substantially improved its visualization. The proposed coherence metrics were also useful in assessing the ensemble coherence for renal MBF imaging. The results also demonstrated that the proposed coherence metric can be reliably estimated from downsampled ensembles (by up to 90%), thus allowing improved computational efficiency for potential applications in real‐time MBF imaging.

Conclusions

This pilot study demonstrates the importance of assessing ensemble coherency in contrast‐free MBF imaging. The proposed LSTC image quantified coherence of the Doppler ensemble for robust MBF imaging. The results obtained from this pilot study are promising, and warrant further development and in vivo validation.

In Vivo Comparison of Multiline Transmission and Diverging Wave Imaging for High-Frame-Rate Speckle-Tracking Echocardiography

High-frame-rate (HFR) speckle-tracking echocardiography (STE) assesses myocardial function by quantifying motion and deformation at high temporal resolution. Among the proposed HFR techniques, multiline transmission (MLT) and diverging wave (DW) imaging have been used in this context both being characterized by specific advantages and disadvantages. Therefore, in this article, we directly contrast both approaches in an in vivo setting while operating at the same frame rate (FR). First, images were recorded at baseline (resting condition) from healthy volunteers and patients. Next, additional acquisitions during stress echocardiography were performed on volunteers. Each scan was contoured and processed by a previously proposed 2-D HFR STE algorithm based on cross correlation. Then, strain curves and their end-systolic (ES) values were extracted for all myocardial segments for further statistical analysis. The baseline acquisitions did not reveal differences in estimated strain between the acquisition modes ( ); myocardial segments ( ); or an interaction between imaging mode and depth ( ). Similarly, during stress testing, no difference ( p = 0.7 ) was observed for the two scan sequences, stress levels or an interaction sequence-stress level ( p = 0.94 ). Overall, our findings show that MLT and DW compoundings give comparable HFR STE strain values and that the choice for using one method or the other may thus rather be based on other factors, for example, system requirements or computational cost.

Tensor Velocity Imaging With Motion Correction

This article presents a motion compensation procedure that significantly improves the accuracy of synthetic aperture tensor velocity estimates for row-column arrays. The proposed motion compensation scheme reduces motion effects by moving the image coordinates with the velocity field during summation of low-resolution volumes. The velocity field is estimated using a transverse oscillation cross-correlation estimator, and each image coordinate's local tensor velocity is determined by upsampling the field using spline interpolation. The motion compensation procedure is validated using Field II simulations and flow measurements acquired using a 3-MHz row-column addressed probe and the research scanner SARUS. For a peak velocity of 25 cm/s, a pulse repetition frequency of 2 kHz, and a beam-to-flow angle of 60°, the proposed motion compensation procedure was able to reduce the relative bias from -27.0% to -9.4% and the standard deviation from 8.6% to 8.1%. In simulations performed with a pulse repetition frequency of 10 kHz, the proposed method reduces the bias in all cases with beam-to-flow angles of 60° and 75° and peak velocities between 10 and 150 cm/s.

Resolution and Speckle Reduction in Cardiac Imaging

Cardiac imaging depends on clear visualization of several different structural and functional components to determine left ventricular and overall cardiac health. Ultrasound imaging is confounded by the characteristic speckle texture resulting from subwavelength scatterers in tissues, which is similar to a multiplicative noise on underlying tissue structure. Reduction of this texture can be achieved through physical means, such as spatial or frequency compounding, or through adaptive image processing. Techniques in both categories require a tradeoff of resolution for speckle texture reduction, which together contribute to overall image quality and diagnostic value. We evaluate this tradeoff for cardiac imaging tasks using spatial compounding as an exemplary speckle reduction method. Spatial compounding averages the decorrelated speckle patterns formed by views of a target from multiple subaperture positions to reduce the texture at the expense of active aperture size (and, in turn, lateral resolution). We demonstrate the use of a novel synthetic aperture focusing technique to decompose harmonic backscattered data from focused beams to their aperture-domain spatial frequency components to enable combined transmit and receive compounding. This tool allows the evaluation of matched data sets from a single acquisition over a wide range of spatial compounding conditions. We quantified the tradeoff between resolution and texture reduction in an imaging phantom and demonstrated improved lesion detectability with increasing levels of spatial compounding. We performed a cardiac ultrasound on 25 subjects to evaluate the degree of compounding useful for diagnostic imaging. Of these, 18 subjects were included in both qualitative and quantitative analysis. We found that compounding improved detectability of the endocardial border according to the generalized contrast-to-noise ratio in all cases, and more aggressive compounding made further improvements in ten out of 18 cases. Three expert reviewers evaluated the images for their usefulness in several diagnostic tasks and ranked four compounding conditions ("none," "low," "medium," and "high"). Contrary to the quantitative metrics that suggested the use of high levels of compounding, the reviewers determined that "low" was usually preferred (77.9%), while "none" or "medium" was selected in 21.2% of cases. We conclude with a brief discussion of the generalization of these results to other speckle reduction methods using the imaging phantom data.

Retrospective Transmit Beamforming and Coherent Plane-Wave Compounding for Microvascular Doppler Imaging: A Comparison Study

Imaging blood flow in small vessels is of great clinical value for evaluating increased vascularization, potentially related to angiogenesis in cancer or inflammation processes in musculoskeletal disease. Using a traditional duplex imaging approach, a major challenge in color Doppler imaging is the limited amount of samples available for clutter filtering. Coherent plane-wave compounding (CPWC) enables a continuous high frame rate acquisition and improved image quality due to dynamic transmit focusing. However, the presence of moving scatterers in the image can lead to a loss in signal-to-noise ratio (SNR) and contrast. In this study, typical CPWC sequences for low-flow imaging were compared with retrospective transmit beamforming (RTB) sequences with similar frame rates and transmit power. The comparison was based on resolution, contrast, and SNR, using a stationary phantom, a flow phantom, a thread phantom, and in vivo recordings of blood vessels in the thyroid and kidney. A model was developed to estimate the difference in SNR between RTB and CPWC in the presence of static and moving scatterers while varying the transmit sequence parameters. The model predicted that RTB may yield an increased SNR compared with CPWC, especially for flow imaging, where the SNR difference reached 6 dB for a maximum velocity of 15 cm/s. The measured SNR values were in agreement with the predicted values, both in the case of stationary scatterers and for the flow phantom. We further demonstrated that reducing beam density to increase frame rate is associated with spatial undersampling (stripe) artifacts for RTB and grating lobes for CPWC. Both phantom and in vivo results indicate that transmit focusing may be beneficial in a low-flow imaging setup that, combined with adaptive clutter filtering, can yield superior microvascular imaging.

CNN-Based Ultrasound Image Reconstruction for Ultrafast Displacement Tracking

Thanks to its capability of acquiring full-view frames at multiple kilohertz, ultrafast ultrasound imaging unlocked the analysis of rapidly changing physical phenomena in the human body, with pioneering applications such as ultrasensitive flow imaging in the cardiovascular system or shear-wave elastography. The accuracy achievable with these motion estimation techniques is strongly contingent upon two contradictory requirements: a high quality of consecutive frames and a high frame rate. Indeed, the image quality can usually be improved by increasing the number of steered ultrafast acquisitions, but at the expense of a reduced frame rate and possible motion artifacts. To achieve accurate motion estimation at uncompromised frame rates and immune to motion artifacts, the proposed approach relies on single ultrafast acquisitions to reconstruct high-quality frames and on only two consecutive frames to obtain 2-D displacement estimates. To this end, we deployed a convolutional neural network-based image reconstruction method combined with a speckle tracking algorithm based on cross-correlation. Numerical and in vivo experiments, conducted in the context of plane-wave imaging, demonstrate that the proposed approach is capable of estimating displacements in regions where the presence of side lobe and grating lobe artifacts prevents any displacement estimation with a state-of-the-art technique that relies on conventional delay-and-sum beamforming. The proposed approach may therefore unlock the full potential of ultrafast ultrasound, in applications such as ultrasensitive cardiovascular motion and flow analysis or shear-wave elastography.

Optimum Speckle Tracking Based on Ultrafast Ultrasound for Improving Blood Flow Velocimetry

Speckle tracking using optimum comparison frames (STO) is proposed to improve the blood flow velocity profile (BFVP) estimation based on ultrafast ultrasound with coherent plane-wave compounding. The optimum comparison frames are as far as possible from the reference frame image while possessing a speckle correlation above a given threshold. The correlation thresholds for different kernel sizes are determined via an experiment based on a vascular-mimicking phantom. In in vitro experiments with different peak velocities of the flow ranging from 0.38 to 1.18 m/s, the proposed STO method with three kernel sizes ( 0.46 × 0.46 , 0.31 × 0.69 , and 0.92 × 0.92 mm²) is used for the BFVP estimations. The normalized root mean square errors (NRMSEs) between the estimated and theoretical BFVPs are calculated and compared with the results based on the speckle tracking using adjacent-frame images. For the three kernel sizes, the mean relative decrements in the STO-based NRMSEs are 46.6%, 44.7%, and 52.9%, and the standard deviations are 36.8%, 37.6%, and 35.9%, respectively. The STO method is also validated by in vivo experiments using rabbit iliac arteries with contrast agents. With parabolic curves fitting to the mean velocity estimates, the average relative increments for the STO-based R² (coefficients of determination) are 7.22% and 6.25% for kernel sizes of 0.46 × 0.46 and 0.31 × 0.69 mm², respectively. In conclusion, the STO method improves the BFVP measurement accuracy, whereby accurate diagnosis information can be acquired for clinical applications.

Space Optimized Plane Wave Imaging for Fast Ultrasonic Inspection with Small Active Aperture: Simulation and Experiment

Plane wave imaging (PWI) is attracting more attention in industrial nondestructive testing and evaluation (NDT&E). To further improve imaging quality and reduce reconstruction time in ultrasonic imaging with a limited active aperture, an optimized PWI algorithm was proposed for rapid ultrasonic inspection, with the comparison of the total focusing method (TFM). The effective area of plane waves and the space weighting factor were defined in order to balance the amplitude of the imaging area. Experiments were carried out to contrast the image quality, with great agreement to the simulation results. Compared with TFM imaging, the space-optimized PWI algorithm demonstrated a wider dynamic detection range and a higher defects amplitude, where the maximum defect amplitude attenuation declined by 6.7 dB and average attenuation on 12 defects decreased by 3.1 dB. In addition, the effects of plane wave numbers on attenuation and reconstruction time were focused on, achieving more than 10 times reduction of reconstruction times over TFM.

High contrast power Doppler imaging in side-viewing intravascular ultrasound imaging via angular compounding

The ability to assess likelihood of plaque rupture can determine the course of treatment in coronary artery disease. One indicator of plaque vulnerability is the development of blood vessels within the plaque, or intraplaque neovascularization. In order to visualize these vessels with increased sensitivity in the cardiac catheterization lab, a new approach for imaging blood flow in small vessels using side-viewing intravascular ultrasound (IVUS) is proposed. This approach based on compounding adjacent angular acquisitions was evaluated in tissue mimicking phantoms and ex vivo vessels. In phantom studies, the Doppler CNR increased from 3.3 ± 1.0 to 13 ± 2.6 (conventional clutter filtering) and from 1.9 ± 0.15 to 7.5 ± 1.1 (SVD filtering) as a result of applying angular compounding. When imaging flow at a rate of 5.6 mm/s in 200 µm tubes adjacent to the lumen of ex vivo porcine arteries, the Doppler CNR increased from 5.3 ± 0.95 to 7.2 ± 1.3 (conventional filtering) and from 23 ± 3.3 to 32 ± 6.7 (SVD filtering). Applying these strategies could allow increased sensitivity to slow flow in side-viewing intravascular ultrasound imaging.

Impact of Aperture, Depth, and Acoustic Clutter on the Performance of Coherent Multi-Transducer Ultrasound Imaging

Transducers with a larger aperture size are desirable in ultrasound imaging to improve resolution and image quality. A coherent multi-transducer ultrasound imaging system (CoMTUS) enables an extended effective aperture through the coherent combination of multiple transducers. In this study, the discontinuous extended aperture created by CoMTUS and its performance for deep imaging and through layered media are investigated by both simulations and experiments. Typical image quality metrics—resolution, contrast and contrast-to-noise ratio—are evaluated and compared with a standard single probe imaging system. Results suggest that the image performance of CoMTUS depends on the relative spatial location of the arrays. The resulting effective aperture significantly improves resolution, while the separation between the arrays may degrade contrast. For a limited gap in the effective aperture (less than a few centimetres), CoMTUS provides benefits to image quality compared to the standard single probe imaging system. Overall, CoMTUS shows higher sensitivity and reduced loss of resolution with imaging depth. In general, CoMTUS imaging performance was unaffected when imaging through a layered medium with different speed of sound values and resolution improved up to 80% at large imaging depths.

A Novel 2-D Speckle Tracking Method for High-Frame-Rate Echocardiography

Speckle tracking echocardiography (STE) is a clinical tool to noninvasively assess regional myocardial function through the quantification of regional motion and deformation. Even if the time resolution of STE can be improved by high-frame-rate (HFR) imaging, dedicated HFR STE algorithms have to be developed to detect very small interframe motions. Therefore, in this article, we propose a novel 2-D STE method, purposely developed for HFR echocardiography. The 2-D motion estimator consists of a two-step algorithm based on the 1-D cross correlations to separately estimate the axial and lateral displacements. The method was first optimized and validated on simulated data giving an accuracy of ~3.3% and ~10.5% for the axial and lateral estimates, respectively. Then, it was preliminarily tested in vivo on ten healthy volunteers showing its clinical applicability and feasibility. Moreover, the extracted clinical markers were in the same range as those reported in the literature. Also, the estimated peak global longitudinal strain was compared with that measured with a clinical scanner showing good correlation and negligible differences (-20.94% versus -20.31%, p -value = 0.44). In conclusion, a novel algorithm for STE was developed: the radio frequency (RF) signals were preferred for the axial motion estimation, while envelope data were preferred for the lateral motion. Furthermore, using 2-D kernels, even for 1-D cross correlation, makes the method less sensitive to noise.

Ultrafast Endoscopic Ultrasonography With Circular Array

Rapid development of ultrafast ultrasound imaging has led to novel medical ultrasound applications, including shear wave elastography and super-resolution vascular imaging. However, these have yet to incorporate endoscopic ultrasonography (EUS) with a circular array, which provides a wider view in the alimentary canal than traditional linear and convex arrays. A coherent diverging wave compounding (CDWC) imaging method was proposed for ultrafast EUS imaging and implemented on a custom circular array. In CDWC, virtual acoustic point sources are allocated and virtually insonified diverging waves from each source are achieved by adjusting all circular array elements' emission time delays. Diverging waves emitted from different virtual sources are coherently compounded, generating synthetic transmit focusing at every location in the image plane. As the field of view of the circular array is centrally symmetric, all virtual sources are equidistantly distributed on a concentric circle of radius r . To achieve the highest frame rate possible with image quality comparable to that obtained with the traditional multi-focus imaging method, the effects of various radii r and virtual source quantities on the compounded image quality were theoretically analyzed and experimentally verified. Simulation, phantom, and ex-vivo experiments were conducted with an 8 MHz, 124-element circular array, with a 5.35 mm radius. When 16 virtual sources were used with r=1.605 mm, image quality comparable to that obtained with the multi-focus approach was achieved at a frame rate of 1000 frames/s. This demonstrates the feasibility of the proposed ultrafast EUS imaging method and promotes further development of multi-functional EUS devices.

High-Frame-Rate Color Doppler Echocardiography: A Quantitative Comparison of Different Approaches

Ultrasound color Doppler imaging (CDI) provides a map of the axial blood flow velocities in a 2-D/3-D region of interest. While CDI is clinically effective for a qualitative analysis of abnormal blood flows, e.g., for valvular disease in cardiology, it is in limited use for quantitative measures, mainly hampered by low frame rate and measurement bias. These limitations can be reduced by different approaches toward high-frame-rate (HFR) imaging at the expense of reduced image quality and penetration depth. The aim of this study was to compare the impact of different HFR sequences on CDI quantitatively. Different cardiac scan sequences, including diverging waves and multiline transmission, were designed, implemented on a research system, and compared in terms of patient safety parameters, image quality, and penetration depth. Furthermore, in vivo images were acquired and compared for healthy volunteers. Results showed that the HFR techniques spread artifacts on larger areas than the standard single-line scans (> +50%). In addition, due to patient safety limitations, they reduce the penetration depth up to -5 cm. On the other hand, the HFR techniques provide comparable velocity estimates (relative difference <6%) and enhance the time resolution of the color Doppler images, achieving frame rates up to 625 Hz in continuous acquisition.

A unified deep network for beamforming and speckle reduction in plane wave imaging: A simulation study

Plane Wave Imaging is a fast imaging method used in ultrasound, which allows a high frame rate, but with compromised image quality when a single wave is used. In this work a learning-based approach was used to obtain improved image quality. The entire process of beamforming and speckle reduction was embedded in a single deep convolutional network, and trained with two types of simulated data. The network architecture was designed based on traditional physical considerations of the ultrasonic image formation pipe. As such, it includes beamforming with spatial matched filters, envelope detection, and a speckle reduction stage done in log-signal representation, with all stages containing trainable parameters. The approach was tested on the publicly available PICMUS datasets, achieving axial and lateral full-width-half-maximum (FWHM) resolution values of 0.22 mm and 0.35 mm respectively, and a Contrast to Noise Ratio (CNR) metric of 16.75 on the experimental datasets.

A United Sign Coherence Factor Beamformer for Coherent Plane-Wave Compounding with Improved Contrast

In this study, we present a united sign coherence factor beamformer for coherent plane-wave compounding (CPWC). CPWC is capable of reaching an image quality comparable to the conventional B-mode with a much higher frame rate. Conventional coherence factor (CF) based beamformers for CPWC are based on one-dimensional (1D) frameworks, either in the spatial coherence dimension or angular coherence dimension. Both 1D frameworks do not take into account the coherence information of the dimensions of each other. In order to take full advantage of the radio-frequency (RF) data, this paper proposes a united framework containing both spatial and angular information for CPWC. A united sign coherence factor beamformer (uSCF), which combines the conventional sign coherence factor (SCF) and the united framework, is introduced in the paper as well. The proposed beamformer is compared with the conventional 1D SCF beamformers (spatial and angular dimension beamformers) using simulation, phantom and in vivo studies. In the in vivo images, the proposed method improves the contrast ratio (CR) and generalized contrast-to-noise ratio (gCNR) by 197% and 20% over CPWC. Compared with other 1D methods, uSCF also shows an improved contrast and lateral resolution on all datasets.

Ultrafast Radial Modulation Imaging

Radial modulation imaging improves the detection of microbubbles at high frequency using a dual ultrasonic excitation. However, the synchronization between the imaging pulses is nontrivial because microbubbles need to be interrogated in the compression and the rarefaction phase, and the time-delay difference from dispersion has to be corrected. To address these issues, we propose the use of ultrafast radial modulation imaging (uRMI). In this technique, a beat frequency between the modulation pulse (around 1 MHz) and the ultrafast pulse-repetition frequency was exploited to separate microbubbles from tissue phantom in vitro. This led to a modulated images' set in the spectral domain of the slow time that may then be demodulated through a digital lock-in amplifier to retrieve the contrast image. Ultrafast RMI, applied on a flow phantom with microbubbles, provided a contrast-to-tissue ratio from 7.2 to 14.8 dB at 15 MHz. For flow speed lower than 0.05 mL/min, uRMI (16 dB) provided a better contrast-to-tissue ratio than other techniques: singular value decomposition spatiotemporal filter (11 dB), amplitude modulation (9 dB), or microbubbles disruption (6 dB). This technique may then be suitable to improve the detection of targeted microbubbles, in ultrasound molecular imaging applications, and the detection of extremely slow microbubbles moving in the finest vessels in ultrasound localization microscopy.