1
|
Naji MA, Taghavi I, Schou M, Praesius SK, Hansen LN, Panduro NS, Andersen SB, Sogaard SB, Gundlach C, Kjer HM, Tomov BG, Thomsen EV, Nielsen MB, Larsen NB, Dahl AB, Sorensen CM, Jensen JA. Super-Resolution Ultrasound Imaging Using the Erythrocytes-Part II: Velocity Images. IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL 2024; 71:945-959. [PMID: 38857146 DOI: 10.1109/tuffc.2024.3411795] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2024]
Abstract
Super-resolution ultrasound imaging using the erythrocytes (SURE) has recently been introduced. The method uses erythrocytes as targets instead of fragile microbubbles (MBs). The abundance of erythrocyte scatterers makes it possible to acquire SURE data in just a few seconds compared with several minutes in ultrasound localization microscopy (ULM) using MBs. A high number of scatterers can reduce the acquisition time; however, the tracking of uncorrelated and high-density scatterers is quite challenging. This article hypothesizes that it is possible to detect and track erythrocytes as targets to obtain vascular flow images. A SURE tracking pipeline is used with modules for beamforming, recursive synthetic aperture (SA) imaging, motion estimation, echo canceling, peak detection, and recursive nearest-neighbor (NN) tracker. The SURE tracking pipeline is capable of distinguishing the flow direction and separating tubes of a simulated Field II phantom with 125-25- [Formula: see text] wall-to-wall tube distances, as well as a 3-D printed hydrogel micr-flow phantom with 100-60- [Formula: see text] wall-to-wall channel distances. The comparison of an in vivo SURE scan of a Sprague-Dawley rat kidney with ULM and micro-computed tomography (CT) scans with voxel sizes of 26.5 and [Formula: see text] demonstrated consistent findings. A microvascular structure composed of 16 vessels exhibited similarities across all imaging modalities. The flow direction and velocity profiles in the SURE scan were found to be concordant with those from ULM.
Collapse
|
2
|
Lin SC, Li PC. Fourier-Based Fast 3-D Ultrasound Imaging Using Row-Column-Addressed 2-D Arrays. IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL 2024; 71:85-101. [PMID: 38060356 DOI: 10.1109/tuffc.2023.3340507] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/10/2024]
Abstract
A Fourier-based fast 3-D ultrasound imaging method using row-column-addressed (RCA) 2-D arrays is presented. The row elements in an RCA array are activated sequentially, and all the column elements are used to receive. The obtained dataset is adapted to approximate to that obtained using a fully sampled array after a plane wave at a given incident angle is transmitted. In this way, the fast algorithm in plane-wave Fourier imaging (PWFI) can be applied to the adapted dataset. In addition, synthesizing multiple datasets based on multiple incident angles enables angular compounding, which improves the image quality. The proposed method was validated using computer simulations and physical-phantom experiments. The results show that the spatial resolution and contrast of the proposed method are comparable with those of its PWFI counterpart without requiring a fully sampled (FS) array. Compared with the delay-and-sum (DAS) method using the RCA array, the proposed method provides comparable spatial resolution but lower contrast; however, the computational complexity is significantly reduced from O(N4Nz) to O(WN2Nz log2(N2Nz)) , where N is the number of elements on each side of the RCA array, Nz is the number of voxels in the axial direction in the output image, and W is the number of compounding angles. For example, in the simulated results when the maximum compounding angle M is 5°, at a given point the lateral - 6-dB width provided by the proposed method is 0.241 mm (0.267 mm for DAS), the contrast ratio of a hyperechoic cyst is 8.87 dB (9.10 dB for DAS), the number of real number operations is reduced by a factor of 20.62, and the number of memory accesses is reduced by a factor of 47.21, both compared with DAS. This novel fast algorithm could facilitate the development of compact real-time 3-D imaging systems, especially when the channel count is high and a large field of view (FOV) is required.
Collapse
|
3
|
Jørgensen LT, Stuart MB, Jensen JA. Transverse oscillation tensor velocity imaging using a row-column addressed array: Experimental validation. ULTRASONICS 2023; 132:106962. [PMID: 36906961 DOI: 10.1016/j.ultras.2023.106962] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/29/2022] [Revised: 02/02/2023] [Accepted: 02/13/2023] [Indexed: 05/29/2023]
Abstract
Tensor velocity imaging (TVI) performance with a row-column probe was assessed for constant flow in a straight vessel phantom and pulsatile flow in a carotid artery phantom. TVI, i.e., estimating the 3-D velocity vector as a function of time and spatial position, was performed using the transverse oscillation cross-correlation estimator, and the flow was acquired with a Vermon 128+128 row-column array probe connected to a Verasonics 256 research scanner. The emission sequence used 16 emissions per image, and a TVI volume rate of 234 Hz was obtained for a pulse repetition frequency (fprf) of 15 kHz. The TVI was validated by comparing estimates of the flow rate through several cross-sections with the flow rate set by the pump. For the constant 8 mL/s flow in the straight vessel phantom with relative estimator bias (RB) and standards deviation (RSD) was found in the range of -2.18% to 0.55% and 4.58% to 2.48% in measurements performed with an fprf of 15, 10, 8, and 5 kHz. The pulsatile flow in the carotid artery phantom the was set to an average flow rate of 2.44 mL/s, and the flow was acquired with an fprf of 15, 10, and 8 kHz. The pulsatile flow was estimated from two measurement sites: one at a straight section of the artery and one at the bifurcation. In the straight section, the estimator predicted the average flow rate with an RB value ranging from -7.99% to 0.10% and an RSD value ranging from 10.76% to 6.97%. At the bifurcation, RB and RSD values were between -7.47% to 2.02% and 14.46% to 8.89%. This demonstrates that an RCA with 128 receive elements can accurately capture the flow rate through any cross-section at a high sampling rate.
Collapse
Affiliation(s)
- Lasse Thurmann Jørgensen
- Center for Fast Ultrasound Imaging, Department of Health Technology, Technical University of Denmark, DK-2800 Lyngby, Denmark.
| | - Matthias Bo Stuart
- Center for Fast Ultrasound Imaging, Department of Health Technology, Technical University of Denmark, DK-2800 Lyngby, Denmark
| | - Jørgen Arendt Jensen
- Center for Fast Ultrasound Imaging, Department of Health Technology, Technical University of Denmark, DK-2800 Lyngby, Denmark
| |
Collapse
|
4
|
Song H, Kang J, Boctor EM. Synthetic radial aperture focusing to regulate manual volumetric scanning for economic transrectal ultrasound imaging. ULTRASONICS 2023; 129:106908. [PMID: 36527822 PMCID: PMC10043828 DOI: 10.1016/j.ultras.2022.106908] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/30/2022] [Revised: 10/18/2022] [Accepted: 11/27/2022] [Indexed: 06/17/2023]
Abstract
In this paper, we present a volumetric transrectal ultrasound (TRUS) imaging under the presence of radial scanning angle disorientation (SAD) in a resource-limited diagnostic setting. Herein, we test our hypothesis that a synthetic radial aperture focusing (TRUS-rSAF) technique, in which a radial plane in target volume is reconstructed by coherent compounding of multiple transmittance/reception events, will reject a randomized SAD in a free-hand scanning setup based on external angular tracking. Based on an analytical model of the TRUS-rSAF technique, we first tested specific scenarios using a clinically available TRUS transducer under different SADs in a range of normal distributions (σ = 0.1°, 0.2°, 0.5°, 1°, 2°, and 5°). We found a benefit of the TRUS-rSAF technique for higher robustness when the SAD is contained within the radial synthetic aperture window, i.e., ±0.71° from a target scanning angle. However, no enhancement was found in spatial resolution because of the limited transmit beam field of the clinical TRUS transducer, limiting the synthetic aperture window. We further evaluated the TRUS-rSAF technique with a modified TRUS transducer for an extended synthetic aperture window to test whether higher spatial resolution and robustness to SAD can be obtained in the same evaluation setup. Widening of the synthetic aperture window (±3.54°, ± 5.91°, ± 8.27°, ± 10.63°, ± 12.99°, ± 15.35°) resulted in proportional enhancements of spatial resolution, but it also progressively built up sidelobe artifacts due to randomized synthesis with limited phase cancellations. The results suggest the need for careful calibration of the TRUS-rSAF technique to enable TRUS imaging with free-hand radial scanning and external angle tracking in resource-limited settings.
Collapse
Affiliation(s)
- Hyunwoo Song
- Department of Computer Science, Whiting School of Engineering, the Johns Hopkins University, Baltimore, MD 21218, USA
| | - Jeeun Kang
- Laboratory for Computational Sensing and Robotics, Whiting School of Engineering, the Johns Hopkins University, Baltimore, MD 21218, USA.
| | - Emad M Boctor
- Department of Computer Science, Whiting School of Engineering, the Johns Hopkins University, Baltimore, MD 21218, USA; Laboratory for Computational Sensing and Robotics, Whiting School of Engineering, the Johns Hopkins University, Baltimore, MD 21218, USA.
| |
Collapse
|
5
|
Jensen JA, Schou M, Jorgensen LT, Tomov BG, Stuart MB, Traberg MS, Taghavi I, Oygaard SH, Ommen ML, Steenberg K, Thomsen EV, Panduro NS, Nielsen MB, Sorensen CM. Anatomic and Functional Imaging Using Row-Column Arrays. IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL 2022; 69:2722-2738. [PMID: 35839193 DOI: 10.1109/tuffc.2022.3191391] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Row-column (RC) arrays have the potential to yield full 3-D ultrasound imaging with a greatly reduced number of elements compared to fully populated arrays. They, however, have several challenges due to their special geometry. This review article summarizes the current literature for RC imaging and demonstrates that full anatomic and functional imaging can attain a high quality using synthetic aperture (SA) sequences and modified delay-and-sum beamforming. Resolution can approach the diffraction limit with an isotropic resolution of half a wavelength with low sidelobe levels, and the field of view can be expanded by using convex or lensed RC probes. GPU beamforming allows for three orthogonal planes to be beamformed at 30 Hz, providing near real-time imaging ideal for positioning the probe and improving the operator's workflow. Functional imaging is also attainable using transverse oscillation and dedicated SA sequence for tensor velocity imaging for revealing the full 3-D velocity vector as a function of spatial position and time for both blood velocity and tissue motion estimation. Using RC arrays with commercial contrast agents can reveal super-resolution imaging (SRI) with isotropic resolution below [Formula: see text]. RC arrays can, thus, yield full 3-D imaging at high resolution, contrast, and volumetric rates for both anatomic and functional imaging with the same number of receive channels as current commercial 1-D arrays.
Collapse
|
6
|
Song H, Kang J, Boctor EM. A novel design framework of synthetic radial aperture focusing for volumetric transrectal ultrasound imaging. JOURNAL OF COMPUTATIONAL DESIGN AND ENGINEERING 2022; 9:1852-1865. [PMID: 36268473 PMCID: PMC9563629 DOI: 10.1093/jcde/qwac083] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/06/2022] [Revised: 04/19/2022] [Accepted: 05/14/2022] [Indexed: 06/16/2023]
Abstract
In this paper, we present a novel design framework of synthetic radial aperture focusing for three-dimensional (3D) transrectal ultrasound imaging (TRUS-rSAF), in which multiple transmittance/reception events at different scanning angles are synthesized to reconstruct a radial plane in the target volume, securing high spatial resolution and texture uniformity. A theory-based design approach has not been available to push the envelope of the 3D rSAF technique. Herein, a closed-form analytical description of the TRUS-rSAF method is presented for the first time, effectively delineating spatial resolution and grating lobe positions in the radial dimension of a TRUS transducer. We demonstrate a solid optimization workflow based on the theoretical foundation to improve its spatiotemporal resolution, grating lobe artifacts, and signal-to-noise ratio. A specific design criterion was considered to outperform a clinical 3D TRUS imaging as a reference (TRUS-REF), where each radial plane is reconstructed with a single transmittance/reception event using a motorized actuator. The optimized TRUS-rSAF method significantly enhanced spatial resolution up to 50% over the TRUS-REF method while providing clinically effective temporal resolution (2-8 volume/sec) with negligible grating lobe artifacts. The results indicate that the proposed design approach would enable a novel TRUS imaging solution in clinics.
Collapse
Affiliation(s)
- Hyunwoo Song
- Department of Computer Science, Whiting School of Engineering, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | | | | |
Collapse
|
7
|
Mozaffarzadeh M, Verschuur DJE, Verweij MD, de Jong N, Renaud G. Accelerated 2-D Real-Time Refraction-Corrected Transcranial Ultrasound Imaging. IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL 2022; 69:2599-2610. [PMID: 35797321 DOI: 10.1109/tuffc.2022.3189600] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
In a recent study, we proposed a technique to correct aberration caused by the skull and reconstruct a transcranial B-mode image with a refraction-corrected synthetic aperture imaging (SAI) scheme. Given a sound speed map, the arrival times were calculated using a fast marching technique (FMT), which solves the Eikonal equation and, therefore, is computationally expensive for real-time imaging. In this article, we introduce a two-point ray tracing method, based on Fermat's principle, for fast calculation of the travel times in the presence of a layered aberrator in front of the ultrasound probe. The ray tracing method along with the reconstruction technique is implemented on a graphical processing unite (GPU). The point spread function (PSF) in a wire phantom image reconstructed with the FMT and the GPU implementation was studied with numerical synthetic data and experiments with a bone-mimicking plate and a sagittally cut human skull. The numerical analysis showed that the error on travel times is less than 10% of the ultrasound temporal period at 2.5 MHz. As a result, the lateral resolution was not significantly degraded compared with images reconstructed with FMT-calculated travel times. The results using the synthetic, bone-mimicking plate, and skull dataset showed that the GPU implementation causes a lateral/axial localization error of 0.10/0.20, 0.15/0.13, and 0.26/0.32 mm compared with a reference measurement (no aberrator in front of the ultrasound probe), respectively. For an imaging depth of 70 mm, the proposed GPU implementation allows reconstructing 19 frames/s with full synthetic aperture (96 transmission events) and 32 frames/s with multiangle plane wave imaging schemes (with 11 steering angles) for a pixel size of [Formula: see text]. Finally, refraction-corrected power Doppler imaging is demonstrated with a string phantom and a bone-mimicking plate placed between the probe and the moving string. The proposed approach achieves a suitable frame rate for clinical scanning while maintaining the image quality.
Collapse
|
8
|
Jorgensen LT, Traberg MS, Stuart MB, Jensen JA. Performance Assessment of Row-Column Transverse Oscillation Tensor Velocity Imaging Using Computational Fluid Dynamics Simulation of Carotid Bifurcation Flow. IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL 2022; 69:1230-1242. [PMID: 35133963 DOI: 10.1109/tuffc.2022.3150106] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
In this work, the accuracy of row-column tensor velocity imaging (TVI), i.e., 3-D vector flow imaging (VFI) in 3-D space over time, is quantified on a complex, clinically relevant flow. The quantification is achieved by transferring the flow simulated using computational fluid dynamics (CFD) to a Field II simulation environment, and this allows for a direct comparison between the actual and estimated velocities. The carotid bifurcation flow simulations were performed with a peak inlet velocity of 80 cm/s, nonrigid vessel walls, and a flow cycle duration of 1.2 s. The flow was simulated from two observation angles, and it was acquired using a 3-MHz 62+62 row-column addressed array (RCA) at a pulse repetition frequency ( fprf ) of 10 and 20 kHz. The tensor velocities were obtained at a frame rate of 208.3 Hz, at fprf = 10 kHz , and the results from two velocity estimators were compared. The two estimators were the directional transverse oscillation (TO) cross correlation estimator and the proposed autocorrelation estimator. Linear regression between the actual and estimated velocity components yielded, for the cross correlation estimator, an R 2 value in the range of 0.89-0.91, 0.46-0.77, and 0.91-0.97 for the x -, y -, and z -components, and 0.87-0.89, 0.40-0.83, and 0.91-0.96 when using the autocorrelation estimator. The results demonstrate that an RCA can, with just 62 receive channels, measure complex 3-D flow fields at a high volume rate.
Collapse
|