Research Article | Open Access
Zhijie Dong, Jihun Kim, Chengwu Huang, Matthew R. Lowerison, U-Wai Lok, Shigao Chen, Pengfei Song, "Three-Dimensional Shear Wave Elastography Using a 2D Row Column Addressing (RCA) Array", BME Frontiers, vol. 2022, Article ID 9879632, 11 pages, 2022. https://doi.org/10.34133/2022/9879632
Three-Dimensional Shear Wave Elastography Using a 2D Row Column Addressing (RCA) Array
Objective. To develop a 3D shear wave elastography (SWE) technique using a 2D row column addressing (RCA) array, with either external vibration or acoustic radiation force (ARF) as the shear wave source. Impact Statement. The proposed method paves the way for clinical translation of 3D SWE based on the 2D RCA, providing a low-cost and high volume rate solution that is compatible with existing clinical systems. Introduction. SWE is an established ultrasound imaging modality that provides a direct and quantitative assessment of tissue stiffness, which is significant for a wide range of clinical applications including cancer and liver fibrosis. SWE requires high frame rate imaging for robust shear wave tracking. Due to the technical challenges associated with high volume rate imaging in 3D, current SWE techniques are typically confined to 2D. Advancing SWE from 2D to 3D is significant because of the heterogeneous nature of tissue, which demands 3D imaging for accurate and comprehensive evaluation. Methods. A 3D SWE method using a RCA array was developed with a volume rate up to 2000 Hz. The performance of the proposed method was systematically evaluated on tissue-mimicking elasticity phantoms and in an in vivo case study. Results. 3D shear wave motion induced by either external vibration or ARF was successfully detected with the proposed method. Robust 3D shear wave speed maps were reconstructed for phantoms and in vivo. Conclusion. The high volume rate 3D imaging provided by the 2D RCA array provides a robust and practical solution for 3D SWE with a clear pathway for future clinical translation.
Shear wave elastography (SWE) is an established ultrasound imaging modality that found significant impact in a wide range of clinical applications including cancer and liver fibrosis [1–4]. SWE estimates the speed of fast propagating shear waves in tissue, which is directly related to tissue stiffness (i.e., shear modulus) and provides an important biomarker for tissue health [5–11]. Shear waves are typically generated by either using the acoustic radiation force (ARF) from the push beam of an ultrasound transducer  or by using an external mechanical vibration coupled to the tissue . Due to the fast propagating speed of shear waves (e.g., 1–10 m/s in most soft tissues), a high imaging frame rate (e.g., 1000-2000 Hz) is typically required for robust SWE . While such high imaging speed is now widely available on mainstream ultrasound scanners powered by software beamformers for 2D imaging, a viable 3D imaging solution that supports such high imaging speed (i.e., >1000 Hz volume rate) remains elusive. As such, current SWE implementations are largely limited to 2D, although 3D SWE has many known advantages over 2D because of the anisotropic tissue mechanical properties and the resulting complex shear wave propagations in 3D . Furthermore, 2D imaging is particularly ill-posed for external vibration-induced shear waves because of the complexities associated with the shear wave source and various boundary conditions .
To date, several solutions for 3D SWE have been proposed. One category uses mechanical scanning of 1D transducers (e.g., wobblers) to acquire 2D SWE images (usually ARF-based) at different elevational slices followed by 3D reconstruction. These methods are not true 3D SWE given that shear wave detection and tracking are still performed in 2D (i.e., within each 2D imaging slice), which does not carry information of tissue viscoelasticity orthogonal to the 2D imaging plane. However, wobbler-based 3D SWE is commercially available and readily accessible to clinicians worldwide .
The other main category of 3D SWE methods is based on 2D matrix arrays, with either external vibration or ARF as the shear wave source. Although a previous study reported the use of ARF for 3D SWE on a element 2D array , in general, it is challenging to implement such an approach on 2D arrays with multiplexing , as multiplexed 2D arrays cannot sustain the long-duration, high-voltage push pulse for ARF. In addition, ARF-induced shear waves are transient and wideband (e.g., 75–500 Hz) , which requires a high volume rate imaging for tracking (e.g., 1000–2000 Hz). Due to the high channel count and expensive computational cost associated with 3D beamforming, it is challenging to reach such high volume rate with 3D SWE.
On the other hand, external vibration-based methods are more feasible for 3D SWE based on 2D arrays because they do not use ARF and the shear wave signal is typically lower in frequency and narrower in bandwidth [21–23]. In particular, for methods using continuous vibrations, thanks to the cyclic nature of the resulting shear wave signal, one can use subapertures of the transducer to acquire subvolumes of the shear wave signal separately and then stitch the subvolumes of shear wave signals together to synthesize the full 3D volume [24, 25]. The cyclic nature of the continuous vibration-induced shear waves can be further exploited to enable 3D SWE on a conventional ultrasound scanner with a very low imaging volume rate. For example, Huang et al.  proposed a sub-Nyquist sampling technique that allows tracking of 3D shear waves with an 88.9 Hz volume rate. However, one common issue with these subvolume and sub-Nyquist sampling-based methods is that multiple data acquisitions are necessary to synthesize the data for the full 3D volume, which makes these methods susceptible to artifacts from tissue motion.
To address these limitations, here we propose a new 3D SWE method based on the 2D row column addressing (RCA) array. Different from the design of a conventional 2D matrix array (e.g., elements), RCA arrays use orthogonally arranged, bar-shaped elements (e.g., elements) to reduce the fabrication and computational complexities of 3D imaging. Thanks to the low channel count, RCA arrays are compatible with most commercial ultrasound scanners and support ultrafast, 3D volumetric imaging with thousands of Hertz of volume rate [27, 28]. Meanwhile, thanks to the absence of multiplexing, RCA arrays can sustain the push pulses for ARF-based 3D SWE. As such, 2D RCA arrays present an enticing solution for robust 3D SWE. Although several studies have presented methods for 3D ultrafast blood flow and super-resolution imaging based on RCA [29, 30], to the best of our knowledge, there have been no studies reporting the use of RCA for 3D SWE. In this paper we will first present a study that investigates different plane wave compounding schemes based on the RCA for 3D shear wave tracking, followed by phantom and in vivo studies that demonstrate 3D SWE using the RCA for both external vibration- and ARF-induced shear waves.
2.1. Imaging Sequence for 3D Shear Wave Detection
Three compounding plane wave imaging methods were studied based on the RCA array. The first scheme is compounding RC, as shown in Figure 1(a), which uses row elements (R) to transmit plane waves with different steering angles and column elements (C) to receive (i.e., the RC scheme). Similarly, the CR scheme uses columns (C) to transmit and rows (R) to receive. The compounding RC+CR combines RC and CR and achieves symmetrical PSF (see Imaging Sequence for 3D Shear Wave Detection in Materials and Methods). The specifications of the sequences are listed in Table 1.
Figure 1(b) shows the simulation results of the three transmission schemes (compounding RC, compounding CR, and compounding RC+CR) at the volume rate of 400 Hz as listed in Table 1. Due to the symmetrical layout of the row elements and column elements, the PSF of the RC scheme is identical to that of the CR scheme except the image is rotated by 90°. Because the compounding RC+CR scheme is a combination of the RC scheme and CR scheme, the PSF is symmetrical in the plane. Since the RC+CR scheme used a smaller angular aperture (for or plane) as compared to RC or CR to maintain the same volume rate, it had weaker transmit focusing and a wider main lobe (i.e., worse spatial resolution). The average full width at half maximum (FWHM) along -direction was calculated using 30 point targets located at depths between 20 and 50 mm. Figure 1(c) shows the similar simulation results of the three transmission schemes at the volume rate of 2000 Hz. Due to the limited number of compounding angles (e.g., 7 angles), the main lobe along the transmit focusing direction was wide. The compounding RC+CR scheme had a better B-mode spatial resolution (measured at -6 dB) in both - and -direction than either RC or CR alone.
The same 3D SWE detection sequences were then used to image a multi-purpose multi-tissue phantom, as shown in Figures 2(a) and 2(b). The FWHM was measured for the wire target located at 41 mm depth and averaged along the -direction. Similar to the results obtained from the simulation, the RC scheme or the CR scheme had better B-mode resolution than the compounding RC+CR scheme at 400 Hz, while the compounding RC+CR scheme had better resolution at 2000 Hz.
Another important metric for 3D SWE is the signal-to-noise ratio (SNR) for shear wave detection, which affects the quality of the detected shear wave signal, especially when the motion is weak such as in the case of ARF-induced shear waves. Figures 2(c) and 2(d) show the measured signals, noises, and SNRs of the three sequences as a function of imaging depth. Due to the symmetry of RC and CR schemes, the SNR performances were identical for these two schemes. Because RC+CR scheme used a smaller angular aperture that distributes more acoustic energy within the field of view (FOV), SNR was higher for RC+CR than either RC or CR when the volume rate was 400 Hz (see Figure 2(c)).
2.2. External Vibration-Based 3D SWE
After studying the RCA-based detection sequence for 3D SWE, we then performed the 3D shear wave study using a homogenous phantom and an elasticity phantom (see External Vibration-Based 3D SWE in Materials and Methods). The imaging sequence details are reported in Table 2 for each application. Figure 3(a) shows the detected 3D shear wave motion at 100 Hz in the homogenous phantom (i.e., the multi-purpose multi-tissue phantom) with a volumetric imaging rate of 400 Hz. Figure 3(b) shows the reconstructed shear wave speed (SWS) volume using the 3D local frequency estimation (LFE) method, and Figure 3(c) shows the corresponding histogram of SWS from the full volume. The average SWS estimated using the proposed 3D SWE method was m/s, while the SWS measured by the LOGIQ E10 ultrasound system (General Electric Healthcare, Wauwatosa, WI, USA) was 2.98 m/s. The nominal SWS of the phantom according to the literature  is m/s. The measurement using the GE system was conducted with an L2-9-D linear array (9 MHz; General Electric Healthcare, Wauwatosa, WI, USA) and the ARF-based 2D elastography mode . The underestimation of SWS using the proposed method may be explained by the lower shear wave frequency (i.e., 100 Hz) from using the external vibrator.
For the elasticity phantom study, a stiff spherical lesion with Young’s modulus of 64.9 kPa (4.65 m/s) was imaged, and Young’s modulus of the background is 18 kPa (2.45 m/s). Figure 4(a) shows representative volumetric frames of detected shear wave motion (200 Hz) as a function of time. Figure 4(b) shows the reconstructed SWS volume of the elasticity phantom with clear visualization of the spherical lesion; the artifact located at the left top corner is probably due to a boundary effect. The measured average SWS of the background was 2.67 m/s with the GE system (Figure. S1) and m/s with the proposed method; the SWS of the stiff target was 4.94 m/s with the GE system (Figure. S1) and m/s with the proposed 3D SWE method.
To evaluate the in vivo performance of the proposed method, an in vivo study was performed in the biceps muscle of a healthy volunteer using another RCA transducer (Daxsonics Ultrasound Inc., Halifax, Canada; central frequency at 7 MHz). The imaging sequence details are reported in Table 2, and the experimental setup is shown in Figure 5(a). Figure 5(b) shows the detected shear wave motion. The estimated SWS using the local cross-correlation method is 3.29 m/s, and the SWS measured with the GE system is 3.84 m/s.
2.3. ARF-Based 3D SWE
Finally, we applied the proposed method to detect the ARF-based shear wave signal. Because of the relatively low element sensitivity of the RCA array used in this study (see Figure. S2), it was challenging to use the same array to generate shear waves. This is only a limitation with the particular RCA array used in this study and does not necessarily apply to other RCA probes. Due to this limitation, we used a 1D linear array transducer (L7-4, Philips Healthcare, Andover, MA, USA) to produce the push beam for shear wave generation for the RCA to detect. Figure 6(a) shows the detected shear wave motion in the slice at 1.5 ms using compounding RC and compounding RC+CR scheme, as listed in Table 2. The detected shear wave using compounding RC scheme has better signal quality compared to the RC+CR scheme, although the RC+CR scheme achieves better B-mode spatial resolution at 2000 Hz (see Figure 2(b)). Figure 6(c) shows the detected shear wave motion at six consecutive time points using the compounding RC scheme at a volume rate of 2000 Hz. The ARF-induced cylindrical-shaped shear wave by the L7-4 linear array can be clearly visualized. The estimated SWS of the homogenous phantom using the random sample consensus (RANSAC) method was 3.05 m/s, as shown in Figure 6(b). The SWS measured using the GE system was 2.75 m/s.
Currently, 3D SWE is challenged by various technical limitations such as the limited volume rate of 3D imaging systems and the high fabrication and computation costs associated with conventional 2D matrix arrays. These limitations remain the major challenge for 3D elastography in clinical practice. In this article, we introduced a 3D SWE method based on the RCA array with shear waves induced by either external vibration or ARF. The phantom studies demonstrate that the proposed method is capable of robust 3D shear wave tracking at high volume rates (e.g., 2000 Hz) without the need for repeated shear wave generation and detection. The in vivo study shows that our method is capable of in vivo shear wave imaging by showing the feasibility of external vibration-based 3D SWE in the arm muscles of a healthy volunteer.
The main advantage of RCA-based 3D SWE is that the RCA array is compatible with most commercial ultrasound systems because of its lower channel count, and the RCA array has lower fabrication costs compared to conventional 2D matrix arrays . For example, our RCA array includes elements with a mm2 footprint, while for a 2D matrix array, elements (i.e., 16,384 channels) would be needed for the same footprint. Meanwhile, a distinct advantage of the RCA array is that it is feasible to achieve a high imaging volume rate (e.g., 1000–2000 Hz) that is challenging for conventional 2D matrix arrays to achieve. Moreover, because 2D RCA arrays do not require multiplexing, they can better sustain the long push pulses for ARF-based SWE. In theory, there should be no technical barriers for RCA to transmit long-duration push pulses.
In this study, three different transmission schemes (i.e., compounding RC, compounding CR, and compounding RC+CR) were evaluated in silico and in vitro. We found that with a fixed angular pitch (i.e., 1°) and a relatively low volume rate (e.g., 400 Hz), the compounding RC scheme or compounding CR scheme achieved better B-mode spatial resolution as compared to the compounding RC+CR scheme, which uses a smaller angular aperture due to the limited number of transmissions (Figure 2(a)). Meanwhile, higher energy gain was achieved by the RC+CR scheme (Figure 2(c)) because of smaller angular aperture, which showed better shear wave detection SNR. In this study, the compounding RC scheme was used for shear wave detection with external vibration due to the better B-mode resolution, as listed in Table 2, which is preferred for detecting high-quality shear waves to reconstruct small objects like the spherical lesion. However, the resolution difference between RC and RC+CR was not significant (Figure 1(b)). For the ARF-based SWE, both compounding RC and compounding RC+CR scheme were investigated, and the detection using compounding RC achieved better shear wave signal quality (Figure 6(a)).
The tradeoff between spatial resolution and SNR needs to be balanced for different applications associated with different volume rates. For example, to detect a small lesion, shear waves with small wavelengths are necessary, and a sequence with high B-mode resolution is preferred because of the more robust spatial sampling. However, if the SNR is more important than spatial resolution (e.g., weak shear wave signal from deep tissues), then a sequence with higher SNR is preferred. For example, one can use the RC+CR scheme or the RC scheme with a reduced angular aperture to distribute maximum acoustic energy within the FOV. Alternatively, the number of transmissions (i.e., number of compounding angles) can be increased, and multiple acquisitions with repeated pushes or vibrations can be used to further boost shear wave SNR. Other advanced transmission schemes can also be applied to improve SNR such as the Hadamard encoding method .
Since RCA is only capable of one-way focusing, the imaging quality is always compromised when the number of compounding angles is reduced to allow for higher imaging volume rate. In our study, a 2000 Hz volume rate (with 6 or 7 compounding angles) was used to detect the ARF-induced shear waves, which lead to noisy shear wave signals (Figure 6) that were challenging for SWS estimation methods such as LFE and cross-correlation. To address this issue, multiple acquisitions with time averaging can be used to improve the shear SNR, although in practice for human imaging a breath hold may be necessary to avoid motion artifacts. Recently several RCA-based adaptive beamforming methods were proposed to improve the RCA imaging quality by taking advantage of the RCA array layout [33, 34]. However, these methods generally involve multiplication between the RC and CR data, which distorts the phase information and is not suitable for detecting shear wave motion.
Apart from compounding plane wave imaging, other transmission schemes can also be used on the 2D RCA array for 3D SWE. For example, synthetic aperture imaging  can offer optimal spatial resolution (see Figure. S3 (d)-(f)) and line-by-line focused beam imaging (see Figure. S3 (g) and (h)) can offer optimal SNR performance. Wide beam imaging is another option that offers both high shear wave detection SNR and spatial resolution (see Figure. S3 (i) and (j)). However, one major limitation of these three methods is that a large number of transmissions is needed, which results in a volume rate that may be too low to be used for 3D SWE. Nevertheless, these sequences may be useful for external vibration-based SWE where shear wave frequency is low. Alternatively one can use repeated cycles of shear wave generation and detection to synthesize shear wave data with these low volume rate sequences .
Both external vibration-based SWE and ARF-based SWE were studied. There are many advantages of external vibration-induced shear waves such as deeper penetration (when using low-frequency vibrations) and lower requirement on imaging frame rate thanks to the cyclic nature of the shear wave. However, ARF-based shear waves have broader bandwidths, which correspond to a higher spatial resolution in SWS images. A wide-bandwidth shear wave is also beneficial for applications involved with shear wave dispersion analysis . In addition, ARF-based SWE is more convenient to implement in clinical practice. In general, if the detection volume rate is sufficient (e.g., 2000 Hz) and ARF is adequate to generate robust shear wave signal, then ARF-based 3D SWE is a preferred choice in practice.
The estimated shear wave speeds of the elasticity phantom are lower than the velocities derived from the GE system, and this underestimation can be explained by the lower shear wave frequency associated with the external vibration compared to the ARF-induced shear waves by the GE system , and it may be corrected for by increasing the vibration frequency.
2D SWE is limited to capturing the complex shear wave propagation in 3D. 2D imaging is particularly ill-posed for the external vibration-based SWE when the propagating direction is not confined to 2D in-plane, under which case SWS estimation will be biased. Another benefit of the 3D SWE using RCA arrays is that high temporal resolution could be achieved (e.g., 100 ms used for the acquisition of the external vibration-based 3D SWE); thus, the proposed method is robust to the artifacts from tissue motion.
In this study, for the first time to the best of our knowledge, RCA-based 3D SWE with external vibration or ARF was developed and verified on phantoms and in vivo with the clinic-ready ultrasound imaging system, and the quantitative measurements were compared with a state-of-the-art clinical 2D ultrasound SWE system. The proposed method can provide 3D elastography at a high volume rate (e.g., 2000 Hz) with short acquisition time (e.g., tens of milliseconds) and fast processing speed by GPU-based computation acceleration. Ongoing studies are being conducted to use a new probe for ARF-based 3D SWE and evaluate its performance for in vivo imaging.
4. Materials and Methods
4.1. Imaging Sequence for 3D Shear Wave Detection
To detect the shear waves generated by either ARF or external vibration, a high volume rate imaging with adequate imaging quality is necessary. The RCA array includes two groups of orthogonally arranged 1D linear arrays (i.e., rows and columns), which transmit and receive in an alternate fashion for 3D imaging (e.g., rows transmit, and columns receive, and vice versa). One limitation of such an arrangement is that RCA arrays only support one-way focusing. Figure 1(a) shows an example of transmission of steered plane waves in the plane (using the row elements) and receive of the plane waves in the plane (with the column elements). In this case, only transmit focusing can be achieved in the plane (i.e., there is no receive focusing), and only receive focusing can be achieved in the plane (i.e., there is no x-z transmit focusing).
Among various transmission schemes, such as compounding plane wave imaging [14, 37], synthetic aperture imaging , and line-by-line focused beam imaging, compounding plane wave imaging provides the best tradeoff for high volume rate scanning because it presents a balanced performance between imaging quality and imaging frame rate. However, due to the limitation of one-way focusing, the point spread function (PSF) of either scheme is asymmetric (e.g., elongated along the -direction when using the RC scheme). A common approach to improve the spatial resolution and obtain a symmetric PSF is to use the RC+CR scheme. However, this method reduces the imaging volume rate by half because it requires twice as many pulse-echo data acquisition cycles. To maintain the same volume rate as the RC or CR scheme, the RC+CR scheme uses half the number of compounding angles in either or plane, and the angular aperture will be reduced by half with a fixed angular pitch.
To study the impact of compounding RC, compounding CR, and compounding RC+CR schemes for 3D SWE (e.g., spatial resolution and SNR), three corresponding sequences with the same volume rate of 400 Hz or 2000 Hz and a fixed angular pitch of 1° were designed and tested to identify the suitable transmission method for robust 3D shear wave detection and tracking. The 400 Hz volume rate is more appropriate for shear wave detection with external vibration, and the 2000 Hz volume is appropriate for the ARF-based shear wave detection. Table 1 summarizes the configurations of the studied sequences. Both simulation and phantom studies were carried out to evaluate the imaging performance of each sequence. Simulations were conducted using the Verasonics simulator with multiple point targets located at different depths (e.g., 20–50 mm). The simulation settings were specified to match the transducer parameters which were used for the phantom study. For the phantom study, a Verasonics Vantage 256 system (Verasonics Inc., Kirkland, WA, USA) and a element 2D RCA array transducer (Vermon SA, Tours, France; central frequency at 5.2 MHz (measured at -20 dB)) were used for the acquisitions. The transmit frequency was 5.2 MHz with a sampling frequency of 20.8 MHz. A multi-purpose multi-tissue ultrasound phantom (Model 040GSE, CIRS Inc., Norfolk, VA, USA) was used for the resolution and SNR evaluation.
The FWHM of the point target was used to evaluate the spatial resolution. To quantify the experimental SNR performance (i.e., for shear wave detection), 100 volumes were acquired from a homogeneous phantom and a region of the phantom with uniform speckle patterns which were used for SNR estimation. The image noise is calculated as the time standard deviation of the volumetric image . where is the mean signal over time and is the number of volumes. The signal, noise, and SNR were averaged along the - and -direction to obtain an accurate profile of the SNR as a function of the depth , and the SNR is given by
Based on the quantitative evaluation of the resolution and SNR in silico and in vitro, proper transmission schemes were selected and used for the rest of the studies.
4.2. External Vibration-Based 3D SWE
An electromechanical vibrator (Bruel & Kjaer VTS Ltd., Royston, UK) was used to generate harmonic shear waves with different frequencies (e.g., 100 Hz–200 Hz) from a function generator (33210A, Keysight Technologies Inc., Santa Rosa, CA, USA). The phantom experimental setup is illustrated in Figure 7(a), and the angle between the vibrator and the phantom was approximately 60 degrees. The sinusoidal signal from the function generator was amplified by an amplifier (XLS 2502, Crown Audio, Los Angeles, CA) before being sent to the vibrator. A metal rod with a bar-shaped tip was used to conduct the vibration from the vibrator to the phantom. The vibrator was synchronized with the ultrasound system. Both the multi-purpose phantom (040GSE) and an elasticity phantom (Model 049, CIRS Inc., Norfolk, VA, USA) were utilized to test the 3D SWE imaging performance.
To avoid aliasing, the volume rate of the detection sequence was set to four times the shear wave frequency, as listed in Table 2. For the homogenous phantom study (i.e., the multi-purpose multi-tissue phantom), a 100 Hz shear wave was generated and detected with a scanning rate of 400 Hz. For the heterogeneous phantom study (i.e., the elasticity phantom), a stiff spherical lesion (10 mm diameter) was imaged. Because the resolution of the shear wave imaging is fundamentally limited by the wavelength of the shear wave, to resolve the spherical lesion from the background with better resolution, a 200 Hz shear wave was used and detected by an 800 Hz RCA detection sequence.
In this study, we developed a GPU-based beamforming method for fast ultrasound data reconstruction. In-phase/quadrature (IQ) data were beamformed and then used to calculate the particle velocity induced by shear waves using the autocorrelation method . A five-pixel spatial window (2.5 wavelengths) along the -direction was used for averaging in the autocorrelation calculation, and a spatial median filter (2.5 wavelengths in three directions) was used to denoise of the calculated shear wave motion. To alleviate shear wave interference inside the phantom from external vibration, a 3D directional filter was applied to the shear wave data to extract shear wave signals that propagate away from the source . 3D SWS map of the homogenous phantom was reconstructed using the 3D local frequency estimation (LFE) method [41, 42], and the 3D SWS maps of the elasticity phantom and in vivo study were reconstructed using the 3D local cross-correlation methods  with a step size of 20 pixels (10 wavelengths) and an averaging window of 40 pixels (20 wavelengths).
4.3. ARF-Based 3D SWE
The experimental setup is depicted in Figure 7(b). The two transducers were carefully aligned and operated by two Verasonics Vantage systems, which were synchronized. A homemade homogenous tissue-mimicking phantom was sandwiched between the two transducers. The L7-4 transmitted the push beam (push beam central frequency, 4 MHz; F-number, 1; focal depth, 3 cm; and push length, 600 us) immediately after the RCA started data acquisition. Similarly, 1D autocorrelation was used for shear wave motion calculation. Due to the low sensitivity of the RCA array (see Figure. S2) and the high volume rate used for detection (i.e., 2000 Hz), the shear wave signal was noisy due to low SNR from the limited number of transmissions (e.g., 7 transmissions), especially along with the transmit focusing direction (e.g., -direction) when only the RC scheme was applied. Furthermore, the ARF-based shear wave motion is much weaker than the external vibration. To best estimate the SWS of the phantom from noisy shear wave signal, the RANSAC  method was used to fit the shear wave motion and calculate SWS.
The data used to support the findings of this study are available from the corresponding author upon request.
Opinions, interpretations, conclusions, and recommendations are those of the author and are not necessarily endorsed by the Department of Defense.
Conflicts of Interest
Pengfei Song, Shigao Chen, and Mayo Clinic have patents and commercial interests for the shear wave elastography technologies presented in this paper, which have been licensed. The remaining authors have no conflicts of interest to declare.
Z.D., S.C., and P.S. generated the idea and designed the experiments. Z.D., J.K., C.H., U.L., S.C., and P.S. developed and applied the shear wave motion detection and shear wave speed estimation algorithms. Z.D., J.K., C.H., and M. R. L. performed the experiments. All authors were involved in the data analysis, critical revision of the manuscript, and discussion.
This work was supported by the Department of Defense (DoD) through the Breast Cancer Research Program (BCRP) under Award Nos. W81XWH-21-1-0062 and W81XWH-21-1-0063.
Figure. S1: shear wave speed maps of the (a) stiffer spherical lesion and (b) background of the elasticity phantom using the GE LOGIC E10 system and a 1D linear array. The measurement was conducted in ARF-based 2D SWE mode. Figure. S2: mechanical index of the RCA array as a function of the input voltage. The measurements were performed using a needle hydrophone in the water tank, and the RCA transmitted a focused beam using 64 elements () or 128 elements (), and the focal depth is 25 mm. Figure. S3: reconstructed volumetric images (slice view) of the phantom using four different transmission schemes (compounding plane wave imaging, synthetic aperture imaging (SA), focused beam imaging (FB), and wide beam imaging (WB)). (a) Compounding plane wave imaging with RC scheme with an angular pitch of 1°, angular aperture of 15°, 31 compounding angles, and 400 Hz volume rate. (b) Compounding plane wave imaging with CR scheme with an angular pitch of 1°, angular aperture of 15°, 31 compounding angles, and 400 Hz volume rate. (c) Compounding plane wave imaging with RC+CR scheme with an angular pitch of 1°, angular aperture of 15°, 62 compounding angles, and 200 Hz volume rate. (d) SA imaging with 128 transmissions and 100 Hz volume rate. one element was used for each transmission. (e) SA imaging with 128 transmissions and 100 Hz volume rate, virtual source (a subaperture of 12 elements) was used for each transmission. (f) SA imaging with 64 transmissions and 200 Hz volume rate, virtual source (a subaperture of 64 elements) was used for each transmission. (g) FB imaging with 64 transmissions and 200 Hz volume rate, 32 elements were used for each focused beam with a F-number of 4. (h) FB imaging with 128 transmissions and 100 Hz volume rate, 32 elements were used for each focused beam with a F-number of 4. (i) WB imaging with 30 transmissions and 400 Hz volume rate, 50 elements were used for each subregion compounding. (j) WB imaging with 60 transmissions and 200 Hz volume rate, 50 elements were used for each subregion compounding. (Supplementary Materials)
- T. Shiina, K. R. Nightingale, M. L. Palmeri et al., “WFUMB guidelines and recommendations for clinical use of ultrasound elastography: part 1: basic principles and terminology,” Ultrasound in Medicine & Biology, vol. 41, no. 5, pp. 1126–1147, 2015.
- R. G. Barr, K. Nakashima, D. Amy et al., “WFUMB guidelines and recommendations for clinical use of ultrasound elastography: part 2: breast,” Ultrasound in Medicine & Biology, vol. 41, no. 5, pp. 1148–1160, 2015.
- R. G. Barr, S. R. Wilson, D. Rubens, G. Garcia-Tsao, and G. Ferraioli, “Update to the society of radiologists in ultrasound liver elastography consensus statement,” Radiology, vol. 296, no. 2, pp. 263–274, 2020.
- G. Ferraioli, P. Parekh, A. B. Levitov, and C. Filice, “Shear wave elastography for evaluation of liver fibrosis,” Journal of Ultrasound in Medicine, vol. 33, no. 2, pp. 197–203, 2014.
- A. P. Sarvazyan, O. V. Rudenko, S. D. Swanson, J. B. Fowlkes, and S. Y. Emelianov, “Shear wave elasticity imaging: a new ultrasonic technology of medical diagnostics,” Ultrasound in Medicine & Biology, vol. 24, no. 9, pp. 1419–1435, 1998.
- L. Sandrin, M. Tanter, S. Catheline, and M. Fink, “Shear modulus imaging with 2-D transient elastography,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 49, no. 4, pp. 426–435, 2002.
- J. Bercoff, M. Tanter, and M. Fink, “Supersonic shear imaging: a new technique for soft tissue elasticity mapping,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 51, no. 4, pp. 396–409, 2004.
- J.-L. Gennisson, T. Deffieux, M. Fink, and M. Tanter, “Ultrasound elastography: principles and techniques,” Diagnostic and Interventional Imaging, vol. 94, no. 5, pp. 487–495, 2013.
- P. Song, M. W. Urban, A. Manduca, H. Zhao, J. F. Greenleaf, and S. Chen, “Comb-push ultrasound shear elastography (CUSE) with various ultrasound push beams,” IEEE Transactions on Medical Imaging, vol. 32, no. 8, pp. 1435–1447, 2013.
- P. Song, H. Zhao, A. Manduca, M. W. Urban, J. F. Greenleaf, and S. Chen, “Comb-push ultrasound shear elastography (CUSE): a novel method for two-dimensional shear elasticity imaging of soft tissues,” IEEE Transactions on Medical Imaging, vol. 31, no. 9, pp. 1821–1832, 2012.
- R. M. Sigrist, J. Liau, A. El Kaffas, M. C. Chammas, and J. K. Willmann, “Ultrasound elastography: review of techniques and clinical applications,” Theranostics, vol. 7, no. 5, pp. 1303–1329, 2017.
- K. Nightingale, S. McAleavey, and G. Trahey, “Shear-wave generation using acoustic radiation force: in vivo and ex vivo results,” Ultrasound in Medicine & Biology, vol. 29, no. 12, pp. 1715–1723, 2003.
- L. Sandrin, B. Fourquet, J. M. Hasquenoph et al., “Transient elastography: a new noninvasive method for assessment of hepatic fibrosis,” Ultrasound in Medicine & Biology, vol. 29, no. 12, pp. 1705–1713, 2003.
- G. Montaldo, M. Tanter, J. Bercoff, N. Benech, and M. Fink, “Coherent plane-wave compounding for very high frame rate ultrasonography and transient elastography,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 56, no. 3, pp. 489–506, 2009.
- S. Rosenzweig, M. Palmeri, N. Rouze et al., “Comparison of concurrently acquired in vivo 3D ARFI and SWEI images of the prostate,” in 2012 IEEE International Ultrasonics Symposium, pp. 97–100, Dresden, Germany, 2012.
- H. Zhao, P. Song, D. D. Meixner et al., “External vibration multi-directional ultrasound shearwave elastography (EVMUSE): application in liver fibrosis staging,” IEEE Transactions on Medical Imaging, vol. 33, no. 11, pp. 2140–2148, 2014.
- S. H. Lee, J. M. Chang, W. H. Kim et al., “Differentiation of benign from malignant solid breast masses: comparison of two-dimensional and three-dimensional shear-wave elastography,” European Radiology, vol. 23, no. 4, pp. 1015–1026, 2013.
- J.-l. Gennisson, J. Provost, T. Deffieux et al., “4-D ultrafast shear-wave imaging,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 62, no. 6, pp. 1059–1065, 2015.
- Q. Huang and Z. Zeng, “A review on real-time 3D ultrasound imaging technology,” BioMed research international, vol. 2017, 20 pages, 2017.
- M. L. Palmeri and K. R. Nightingale, “Acoustic radiation force-based elasticity imaging methods,” Interface focus, vol. 1, no. 4, pp. 553–564, 2011.
- P.-Y. Chen, T.-H. Yang, L.-C. Kuo, H.-Y. Hsu, F.-C. Su, and C.-C. Huang, “Evaluation of hand tendon elastic properties during rehabilitation through high-frequency ultrasound shear elastography,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 68, no. 8, pp. 2716–2726, 2021.
- C.-C. Weng, P.-Y. Chen, D. Chou, C.-C. Shih, and C.-C. Huang, “High frequency ultrasound elastography for estimating the viscoelastic properties of the cornea using lamb wave model,” IEEE Transactions on Biomedical Engineering, vol. 68, no. 9, pp. 2637–2644, 2020.
- P.-Y. Chen, T.-H. Yang, L.-C. Kuo, C.-C. Shih, and C.-C. Huang, “Characterization of hand tendons through high-frequency ultrasound elastography,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 67, no. 1, pp. 37–48, 2020.
- M. Wang, B. Byram, M. Palmeri, N. Rouze, and K. Nightingale, “Imaging transverse isotropic properties of muscle by monitoring acoustic radiation force induced shear waves using a 2-D matrix ultrasound array,” IEEE Transactions on Medical Imaging, vol. 32, no. 9, pp. 1671–1684, 2013.
- M. Wang, B. Byram, M. Palmeri, N. Rouze, and K. Nightingale, “On the precision of time-of-flight shear wave speed estimation in homogeneous soft solids: initial results using a matrix array transducer,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 60, no. 4, pp. 758–770, 2013.
- C. Huang, P. Song, D. C. Mellema et al., “Three-dimensional shear wave elastography on conventional ultrasound scanners with external vibration,” Physics in Medicine & Biology, vol. 65, no. 21, article 215009, 2020.
- C. E. Morton and G. R. Lockwood, “Theoretical assessment of a crossed electrode 2-D array for 3-D imaging,” in IEEE Symposium on Ultrasonics, 2003, vol. 1, pp. 968–971, Honolulu, HI, USA, 2003.
- H. Bouzari, M. Engholm, S. I. Nikolov, M. B. Stuart, E. V. Thomsen, and J. A. Jensen, “Imaging performance for two row–column arrays,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 66, no. 7, pp. 1209–1221, 2019.
- J. Sauvage, J. Porée, C. Rabut et al., “4D functional imaging of the rat brain using a large aperture row-column array,” IEEE Transactions on Medical Imaging, vol. 39, no. 6, pp. 1884–1893, 2020.
- J. A. Jensen, M. L. Ommen, S. H. Oygard et al., “Three-dimensional super-resolution imaging using a row–column array,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 67, no. 3, pp. 538–546, 2020.
- P. Song, M. Macdonald, R. Behler et al., “Two-dimensional shear-wave elastography on conventional ultrasound scanners with time-aligned sequential tracking (TAST) and comb-push ultrasound shear elastography (CUSE),” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 62, no. 2, pp. 290–302, 2015.
- P. Gong, P. Song, and S. Chen, “Hadamard-encoded multipulses for contrast-enhanced ultrasound imaging,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 64, no. 11, pp. 1674–1683, 2017.
- J. Hansen-Shearer, M. Lerendegui, M. Toulemonde, and M.-X. Tang, “High contrast ultrafast 3D ultrasound imaging using row column specific frame multiply and sum,” 2021, https://arxiv.org/abs/2102.07746.
- A. Bertolo, J. Sauvage, M. Tanter, M. Pernot, and T. Deffieux, “XDoppler: cross-correlation of orthogonal apertures for 3D blood flow imaging,” IEEE Transactions on Medical Imaging, vol. 40, no. 12, pp. 3358–3368, 2021.
- C. Amador, M. W. Urban, S. Chen, and J. F. Greenleaf, “Loss tangent and complex modulus estimated by acoustic radiation force creep and shear wave dispersion,” Physics in Medicine & Biology, vol. 57, no. 5, pp. 1263–1282, 2012.
- M. W. Urban, B. Qiang, P. Song, I. Z. Nenadic, S. Chen, and J. F. Greenleaf, “Investigation of the effects of myocardial anisotropy for shear wave elastography using impulsive force and harmonic vibration,” Physics in Medicine & Biology, vol. 61, no. 1, pp. 365–382, 2016.
- M. Flesch, M. Pernot, J. Provost et al., “4D in vivo ultrafast ultrasound imaging using a row-column addressed matrix and coherently-compounded orthogonal plane waves,” Physics in Medicine & Biology, vol. 62, no. 11, pp. 4571–4588, 2017.
- E. Mace, G. Montaldo, B.-F. Osmanski, I. Cohen, M. Fink, and M. Tanter, “Functional ultrasound imaging of the brain: theory and basic principles,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 60, no. 3, pp. 492–506, 2013.
- C. Kasai, K. Namekawa, A. Koyano, and R. Omoto, “Real-time two-dimensional blood flow imaging using an autocorrelation technique,” IEEE Transactions on sonics and ultrasonics, vol. 32, no. 3, pp. 458–464, 1985.
- T. Deffieux, J.-L. Gennisson, J. Bercoff, and M. Tanter, “On the effects of reflected waves in transient shear wave elastography,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 58, no. 10, pp. 2032–2035, 2011.
- H. Knutsson, C.-F. Westin, and G. Granlund, “Local multiscale frequency and bandwidth estimation,” in Proceedings of the 1st IEEE International Conference on Image Processing, vol. 1, pp. 36–40, Austin, Texas, 1994.
- A. Manduca, R. Muthupillai, P. J. Rossman, J. F. Greenleaf, and R. L. Ehman, “Image processing for magnetic-resonance elastography,” in Medical Imaging 1996: Image Processing, 1996, vol. 2710, pp. 616–623, International Society for Optics and Photonics.
- M. A. Fischler and R. C. Bolles, “Random sample consensus,” Communications of the ACM, vol. 24, no. 6, pp. 381–395, 1981.
Copyright © 2022 Zhijie Dong et al. Exclusive Licensee Suzhou Institute of Biomedical Engineering and Technology, CAS. Distributed under a Creative Commons Attribution License (CC BY 4.0).