 Research
 Open Access
 Published:
Shortlag spatial coherence imaging using minimum variance beamforming on dual apertures
BioMedical Engineering OnLinevolume 18, Article number: 48 (2019)
Abstract
Background
Shortlag spatial coherence (SLSC) imaging, a newly proposed ultrasound imaging scheme, can offer a higher lesion detectability than conventional Bmode imaging. It requires a high focusing quality which can be satisfied by the synthetic aperture imaging mode. However, traditional nonadaptive synthesis for the SLSC still offers an unsatisfactory resolution. The spatial coherence estimation on the receive aperture cannot fully utilize the coherence information in twodimensional (2D) echo data.
Methods
To overcome these drawbacks, an improved SLSC scheme with adaptive synthesis on dual apertures is proposed in this paper. The minimum variance (MV) beamformer is applied in synthesizing both the receiving and transmitting apertures, while the SLSC function is estimated on both apertures as well. In this way, the resolution is enhanced by the MV implementation, while the coherence in dual apertures is fully utilized.
Results
Simulations, phantom experiments, and in vivo studies are conducted to evaluate the performance of the proposed method. Results demonstrate that the proposed method achieves the best performance in terms of the contrast ratio (CR), contrasttonoise ratio (CNR), and the speckle signaltonoise ratio (SNR). Specifically, compared with the delayandsum (DAS) method, the proposed method achieves 42.5% higher CR, 412.7% higher CNR, and 402.9% higher speckle SNR in simulations. The resolution is also better than the DAS and conventional SLSC beamformers.
Conclusions
The proposed method is a promising technique for improving the SLSC imaging quality and can provide better visualization for medical diagnosis.
Introduction
Shortlag spatial coherence (SLSC) imaging is a newly proposed beamforming technique for medical ultrasound imaging [1, 2]. Different from traditional beamforming methods, it compounds the spatial coherence rather than the magnitude of the array signal. Compared with conventional Bmode methods, the SLSC can better detect lesions inside tissues [1]. In addition, it can provide satisfactory contrast and lesion detectability even under a low signaltonoise (SNR) condition [3]. Previous research has proven that the SLSC imaging can be applied to clinical research and has the potential of obtaining better imaging quality [4,5,6].
The Van Cittert–Zernike theorem is the basis of the SLSC method [2, 7, 8]. It illustrates that the good performance of the SLSC can only be achieved when the transmitted wave is sufficiently focused. In consideration of this principle, a substantial number of conventional imaging modalities are not suitable for the SLSC application. For instance, the commonly used line scan mode can achieve sufficient focus around the focal depth. However, beams are highly diffracted in the nearfield and the deep regions, which means that the image quality inside these areas could suffer a severe degradation [9]. This limits the diagnostic value of the SLSC imaging. As another example, the acoustic wave in the plane wave imaging (PWI) is divergent. Consequently, the combination of the PWI and the SLSC results in poor imaging quality. To meet the requirements of the high focusing quality, it was suggested to apply the synthetic aperture (SA) focusing for the SLSC imaging [9]. Since multitransmission imaging modalities such as the SA imaging and the plane wave compounding (PWC) can naturally achieve the high focusing quality by the dynamic aperture focusing, the application of the SA imaging could improve the SLSC imaging quality mainly in terms of higher contrast [9].
Nevertheless, the inherent drawback of the SLSC imaging still persists. Despite the application of multitransmission modalities, the SLSC images still suffer a lower resolution than conventional Bmode ones [3]. To improve the resolution, one way is to modify the physical setup of the SLSC, such as changing the array aperture. Another alternative is to select an appropriate lag value, which is a vital parameter in the SLSC beamforming. Previous research confirmed that a large lag value could bring a comparably better resolution [2]. However, a small lag value is commonly adopted for a higher contrasttonoise ratio (CNR) to distinguish hyperechoic structures. In addition, the small lag value is also necessary for a high SNR [1, 3]. How to make a compromise between the resolution and the CNR is still a conundrum. To pursue a better SLSC imaging quality, several groups including ours proposed some improved algorithms from another perspective [10, 11]. Zhao et al. [11] proposed a novel method which adaptively synthesizes the transmit (Tx) aperture using a minimum variance (MV) beamformer and then calculates the spatial coherence function. According to previous simulations and experiments, the integration of MV and SLSC could bring better resolution and speckle performance, which give us inspiration for further study.
To improve the resolution of the SLSC, we propose an improved SLSC imaging method using adaptive MV beamforming on dual apertures (DAMV SLSC). In this method, the echo signals of the SA imaging are first recorded and written into a twodimensional (2D) echo data matrix. Then the MV method is applied in beamforming the Tx aperture, and the spatial coherence is estimated in the receive (Rx) aperture. Inspired by the jointtransmittingreceiving beamformer [12], we implement another MV SLSC process but this time, the Rx aperture is MV beamformed, while the SLSC calculation is through the Tx aperture. Finally, the shortlag coherences of both apertures are summed to generate the SLSC image. The innovation point of the method comes from two aspects. First, the MV weighting process is integrated with the SLSC imaging to improve the focusing quality. Although the combination of the MV and the coherence factor (CF) can be frequently found in previous research, the proposed method is conceptually different. Unlike the conventional Bmode MV beamformers with a CFbased postfilter, it is still an SLSC imaging scheme which adopts the MV algorithm for adaptive aperture synthesizing. Second, the spatial coherences are calculated in both the Tx and Rx apertures, which fully utilize the coherence information contained in the 2D data matrix. To evaluate the proposed method, simulations, experiments, and in vivo studies are conducted. Results demonstrate that the proposed method performs better in the resolution and contrast. Detailed analysis will follow in the sections on results and discussion.
The rest part of this paper is organized as follows: “Backgrounds” section presents the background of some previous methods. “Methods” section describes the algorithm of the proposed method and the experimental setup. In “Results” section, the simulated and experimental studies are illustrated in detail, including the comparison between different methods. Finally, in “Discussion” section, the improvements of the newly proposed method versus the performance of the conventional Bmode imaging method, as well as the advantages of the former, are discussed, and the conclusions drawn based on the findings of this study with the suggested scope for future study are given in “Conclusion” section.
Backgrounds
Synthetic aperture imaging model
The SA imaging originates from the radar system and is then applied into medical ultrasound imaging [13,14,15]. In this paper, the synthetic transmit aperture technique is used, which means that only one element fires for each transmitting event. After the propagation through the region of interest (ROI), the backscattered echoes are recorded by all elements. The firing events switch through all transmitting elements to synthesize the full aperture. Suppose an Melement array is adopted, we can acquire an \(M \times M\) size 2D echo data matrix after all data acquisitions:
where n is the time step and \(x_{i,j} \left( n \right)\) represents the echo data recorded by the jth element when using the ith transmitting element. An appropriate time delay is compensated into each channel. In the traditional delayandsum (DAS) algorithm, for each imaging point, all channel data are averaged through rows to generate lowresolution images at first. Then these images can be averaged for the final highresolution image:
where \(\omega_{\text{R}}\) is the receiving apodization window and \(\omega_{\text{T}}\) is the transmitting one. The \(n_{i,j}\) stands for the time index of \(x_{i,j} \left( n \right)\), which is calculated by the positions of the transmitting element, the receiving element, and the imaging point.
Shortlag spatial coherence algorithm
For echo signals of one single emission, the SLSC imaging algorithm differs a lot from the DAS beamformer. When a transducer array with M elements is adopted, the backscattered echo arriving at the array after a firing event could be expressed as
where \(x_{j} \left( n \right)\) represents the echo signal received by the jth element with an appropriated time delay for each channel, n is the time index, and \(\left( \cdot \right)^{\text{T}}\) stands for the matrix transpose. Then the mutual intensity or the socalled spatial coherence between two receiving elements could be estimated from their crosscorrelations:
where m is the distance or lag between two elements in terms of number of elements. The temporal sampling kernel \(n_{2}  n_{1}\) is usually shorter than one cycle [2, 3, 9]. The final SLSC value is calculated using the integral of the first q lags:
where q is the number of lags used for summation. An important parameter for the SLSC is Q = q/M × 100%, which is defined as the ratio between the lag number and the element number. Generally speaking, the lower the Q is, the less the highfrequency component is taken into account, which means the contrast would increase, while the resolution could be damaged, and vice versa [2]. Normally the parameter Q is set to range from 10% to 30%, in order to make a compromise between the resolution and contrast [2, 3].
The SLSC is based on the assumption that the focusing errors due to the acoustic clutter mainly influence the correlations of the “shortlag” elements which are closely separated [2]. This assumption illustrates why only signals of closed elements are used to compute the spatial coherence function. Another premise of the SLSC is that the hyperechoic target and the speckle region have a great difference in their spatial coherence function, which means that a high focusing quality is a prerequisite [2]. The conventional line scan ultrasound imaging mode could only ensure the satisfactory focusing quality near the focal zone, while the shortlag coherence function severely decreases in the region far away from the focus [9]. Hence, the SLSC imaging quality suffers a great degradation in these regions. As a solution, the SA imaging modality can help handle this drawback since the SA imaging dataset could be dynamically focused for each imaging point. In addition, on the basis of the acoustic reciprocity, the spatial coherence can be calculated in either the synthetic Tx aperture or the Rx aperture [16], which is this paper’s starting point.
Minimum variance beamformer
Adaptive beamformers replace the uniform weight used in the DAS beamformer with a datadependent weighting vector \(\varvec{w}\left( n \right)\) for the signal summation. The most representative one, MV beamformer, was originally proposed by Capon in 1969 [17]. It calculates the adaptive weights by minimizing the output energy subject to the constraint that the desired signal is distortionless. Under the incoherent premise between the signal and noise, it is theoretically equal to minimizing the output noise. The MV beamformed output is a product of the received signal vector and the weighting vector:
where \(\varvec{x}\left( n \right)\) is the echo signal in Eq. (3), \(\varvec{w}\left( n \right)\) is the MV weighting vector, and \(\left( \cdot \right)^{\text{H}}\) stands for the conjugate transpose. Then the MV beamforming process could be expressed by
where \(\varvec{a}\) stands for the steering vector, which can all be assumed as one vector after we apply the time delay to each channel. \(\varvec{R}\left( n \right)\) is the covariance matrix of the array data. Using the Lagrange method, the solution to (7) can be calculated as below:
In practice, the covariance matrix \(\varvec{R}\left( n \right)\) is usually unknown and should be estimated by the temporal array samples. A simple approximation could be \(\varvec{R}\left( n \right) = E\left[ {\varvec{x}\left( n \right)\varvec{x}^{\text{H}} \left( n \right)} \right]\) on the basis of the assumption that the signal and noise are uncorrelated. However, in practical ultrasound imaging, the signal and noise are usually highly correlated. To enhance the robustness, the spatial smoothing is proposed [18, 19]. It divides the array into several overlapped short subarrays and averages the covariance matrix through all subarrays. In this way, the onaxis signals can be decorrelated from offaxis ones. The spatial smoothing could be expressed as
where \(\varvec{x}_{l} = \left[ {x_{l} \left( n \right),x_{l + 1} \left( n \right), \ldots ,x_{l + L  1} \left( n \right)} \right]^{\text{T}}\) stands for the lth subarray. L is the subarray length and is usually a little smaller than half of the array size [19, 20]. Another frequently used method to improve the robustness is the diagonal loading. It adds a small constant to the original covariance matrix: \(\varvec{R}_{\text{DL}} \left( n \right) = \varvec{R}\left( n \right) + \varepsilon \cdot \varvec{I}\), which can make the covariance matrix nonsingular. Here \(\varepsilon\) is a small constant, and \(\varvec{I}\) is the identity matrix. The constant \(\varepsilon\) is often set to be \(\Delta \cdot {\text{trace}}\left( {\varvec{R}\left( n \right)} \right)\), while \(\Delta\) is usually less than 0.1 [19, 20].
Using the modified covariance matrix \(\varvec{R}_{\text{DL}} \left( n \right)\), the MV weighting vector \(\varvec{w}_{\text{MV}}\) can be computed using Eq. (8). Considering the spatial smoothing technique, Eq. (6) should be rewritten to form
Methods
Transmit aperture weighting for the SLSC imaging
The adaptive beamformer can be integrated with the SLSC process. Accordingly, a twostep signalprocessing procedure consisting of both the MV process and the SLSC imaging is logically proposed. The SA dataset in “Synthetic aperture imaging model” section, with the appropriate time delay compensation, can be obtained after data collection. The first step is to conduct the adaptive focusing process in the Tx aperture using the MV beamformer. After synthesizing the Tx aperture, the second step is to calculate the spatial coherence in the Rx aperture [11].
In the first step, the 2D echo data matrix is MV beamformed in the Tx dimension to synthesize the Rx aperture. To calculate the Tx–MV weights, the covariance matrix of the Tx aperture is averaged through all Rx array elements on the basis of the multiwave beamforming approach [21, 22]. In addition, since a temporal kernel is adopted in the SLSC algorithm, it is also necessary to adopt the temporal smoothing in the MV process, which averages the covariance matrix through the time index \(n_{1}\) to \(n_{2}\) (the temporal kernel). After all estimations, the covariance matrix should become as
where \(\varvec{x}_{{{\text{Tx}}_{j} ,l}} \left( n \right) = \left[ {x_{l,j} \left( n \right),x_{l + 1,j} \left( n \right), \ldots ,x_{l + L  1,j} \left( n \right)} \right]^{\text{T}}\) represents the lth Tx subarray of the jth Rx element data. Previous research has proven that the averaging process on both Tx aperture and Rx elements could help increase the accuracy of the estimation, which improves the imaging quality [12, 22]. Then the Tx–MV weight \(\varvec{w}_{{{\text{Tx}}  {\text{MV}}}} \left( n \right)\) can be calculated by Eq. (8) using the estimated \(\varvec{R}_{\text{Tx}} \left( n \right)\) with the appropriate diagonal loading. After the MV weighting process in the Tx aperture, the Txsynthesized Rx aperture can be written as \(\varvec{z}_{\text{Rx}} \left( n \right) = \left[ {z_{{{\text{Rx}}_{1} }} \left( n \right),z_{{{\text{Rx}}_{2} }} \left( n \right), \ldots ,z_{{{\text{Rx}}_{M} }} \left( n \right)} \right]\), where the Txweighted output of a single Rx channel \(z_{{{\text{Rx}}_{j} }} \left( n \right)\) can be calculated as
The second step is to estimate the spatial coherence through the Rx aperture, namely, the Rx spatial coherence function, \(\hat{R}_{\text{Rx}} \left( m \right)\):
Receive aperture weighting for the SLSC imaging
On the basis of the acoustic reciprocity, the adaptive weighting process can be implemented in both the Tx aperture and the Rx one [16]. In “Transmit aperture weighting for the SLSC imaging” section, the adaptive weighting process is done in the Tx aperture, and the spatial coherence function is calculated through the Rx aperture. Correspondingly, another MV–SLSC process can be implemented. However, this time the Rx aperture focusing process is first conducted using the MV beamformer, and then the SLSC function could be estimated through the Tx aperture. Similarly, the covariance matrix of the Rx aperture is estimated by averaging through not only all Rx subarrays, but also all Tx array elements:
where \(\varvec{x}_{{{\text{Rx}}_{i} ,l}} \left( n \right) = \left[ {x_{i,l} \left( n \right),x_{i,l + 1} \left( n \right), \ldots ,x_{i,l + L  1} \left( n \right)} \right]^{\text{T}}\) represents the lth Rx subarray of the ith Tx element data. The Rx–MV weight \(\varvec{w}_{{{\text{Rx}}  {\text{MV}}}} \left( n \right)\) can also be calculated by Eq. (8). After the MV weighting process in the Rx aperture, the Rxsynthesized Tx aperture can be written into \(\varvec{z}_{\text{Tx}} \left( n \right) = \left[ {z_{{{\text{Tx}}_{1} }} \left( n \right),z_{{{\text{Tx}}_{2} }} \left( n \right), \ldots ,z_{{{\text{Tx}}_{M} }} \left( n \right)} \right]\), where the Rxweighted output of a single Tx element \(z_{{{\text{Tx}}_{i} }} \left( n \right)\) can be calculated as
As mentioned earlier, the spatial coherence can also be estimated through the Tx aperture. The Tx spatial coherence function \(\hat{R}_{\text{Tx}} \left( m \right)\) can be calculated as follows:
SLSC image formation
The beamforming processes in “Transmit aperture weighting for the SLSC imaging” and “Receive aperture weighting for the SLSC imaging” sections are quite similar, or “symmetric” in a more precise way. If we transpose the SA data matrix in (1), and then repeat the MV SLSC process described in “Transmit aperture weighting for the SLSC imaging” section on the transposed matrix, we can find that both the procedure and the result are the same as described in the “Receive aperture weighting for the SLSC imaging” section. It does tally with the acoustic reciprocity mentioned before. There is no precedence in the calculation of the Rx spatial coherence or the Tx one. Both implementations are based on the original echo data (1). Reasonably both Rx and Tx spatial coherence functions can be used for the SLSC image formation.
Ordinarily, the SLSC image formation process can be realized through Eq. (5). In this paper, we adopt a modified version to accord with the aperture control technique [11]. The fnumber, defined as the ratio between the imaging depth and the effective aperture size, is commonly used to keep a consistent resolution through all depths in medical ultrasound imaging. In this paper’s SA implementation, since the MV weighting process is conducted in both apertures and so is the SLSC estimation process, the same fnumber value is adopted in both the Tx and Rx apertures in order to make two aperture sizes identical. Another reason is to meet the requirements of the reciprocity of the acoustics. An equal size of the transmitting and receiving apertures can maintain the same weight of two apertures in the SLSC calculation, which helps average and reduce the focusing errors in the two apertures. Since the aperture size varies through different imaging depths, the SLSC summation process should be weighted to maintain similar signal magnitudes in the whole imaging region. The modified version of Eq. (5) could be expressed as
Since q is correlated with the imaging depth, the averaging factor 1/q can help keep a consistent signal magnitude in different depths in the final SLSC image.
Implementation summary of the algorithm
To make the proposed scheme clear, we give a brief summary of the procedures here.

1.
Transform the received radio frequency (RF) data to the inphase and quadrature (IQ) domains with the appropriate time delay for each channel to obtain the echo data matrix (1).

2.
Adaptively synthesize the Tx aperture using the MV beamformer by Eqs. (11) and (12). Then estimate the SLSC function in the Rx aperture using Eq. (13).

3.
Accordingly, synthesize the Rx aperture by Eqs. (14) and (15) and estimate the Tx SLSC function using Eq. (16).

4.
Sum both the Rx and Tx spatial coherence function by Eq. (17) for the SLSC image formation.

5.
Repeat the procedure (1) to (4) for each imaging point.
After all beamforming process, the SLSC image could be directly displayed without a dynamic compression, which is different from the conventional Bmode imaging.
Experimental setup
Both the simulations and experiments were conducted to evaluate the performance of the proposed method. Simulated data were acquired using the Matlab simulation tool Field II [23, 24]. All phantom and in vivo data were generated from the Verasonics ultrasound platform (V1, Verasonics, Redmond, Washington).
In the Field II simulation, a 5MHz, 128element transducer array with a 0.3 mm pitch was adopted. The sampling rate was 40 MHz and a twocycle sinusoid was used for the excitation pulse. To study the point spread function (PSF), four point targets at z = 30 mm, 35 mm, 40 mm, and 45 mm were first simulated. A 10 dB zeromean white Gaussian noise (WGN) was added into each channel to simulate the background noise. In the second simulation, two 2.5mm radius circular anechoic cysts centered at z = 37 mm and z = 52 mm were simulated. The thickness of the cyst through the yaxis was set to be 8 mm. The speckle region was generated using randomly distributed and random amplitude scatterers with a total amplitude of 10 per wavelength cubic (\(\lambda^{3}\)).
In the phantom and in vivo experiments, a 5MHz, 128element linear transducer array with a 0.3 mm pitch (L114v, Verasonics, Redmond, Washington, USA) was used. The sampling frequency was originally set to 20 MHz, and the recorded RF data were then resampled at 40 MHz to maintain the same time accuracy with the simulations. All phantom data were obtained using a CIRS calibration phantom (Model 040GSE, Computerized Imaging Reference Systems Inc., Norfolk, Virginia, USA). In vivo human carotid artery data were acquired from a 28yearold male volunteer.
For the SA imaging modality, a spherical wave was transmitted by a single element and recorded by all elements for each emission. The fnumber was set to 1.4 using a rectangular window for both Tx and Rx apertures. A subarray length of L = 0.3M was used for the spatial smoothing process in the MV algorithm. The diagonal loading parameter \(\Delta\) was set to be 0.01. For the SLSC procedure, the ratio Q = q/M was set to be 15% and the temporal kernel was the same as one cycle [2, 3, 9]. These parameters are presented in Table 1.
Different results using the DAS, the MV, the SLSC and the proposed DAMV SLSC are shown together to make an exhaustive comparison. To quantify the performance of different methods, the fullwidth at halfmaximum (FWHM, defined as − 6 dB beam width for the mainlobe) was adopted for point targets, while the contrast ratio (CR), the CNR, and the speckle SNR were measured for cyst targets. The CR, CNR, and speckle SNR can be calculated by
where \(\mu_{\text{b}}\) and \(\mu_{\text{c}}\) are the mean signal magnitudes of the speckle and cyst regions, \(\sigma_{\text{b}}\) and \(\sigma_{\text{c}}\) are the standard deviations of the signal magnitude in the speckle and cyst regions, respectively [2].
Results
Simulated study
Figure 1 shows images of point targets using different methods. In Fig. 1a, sidelobes are visibly obvious with the DAS beamformer. Using the adaptive MV beamformer, sidelobes could be effectively suppressed in Fig. 1b. The original SLSC result in Fig. 1c seems to be worse than the DAS Bmode image since the points are wider. According to previous research, the wide mainlobe width could occur when the background noise level is very low [3]. As shown in Fig. 1d, the proposed DAMV SLSC can narrow the mainlobe when compared with conventional SLSC method. The sidelobes are also suppressed to a low degree.
For quantitative measurements, plots of lateral variations of different methods at the depth z = 30, 35, 40 and 45 mm are shown in Fig. 2 while statistical results of corresponding FWHMs are presented in Table 2. It is shown that our method achieves better performance than not only the SLSC but also the DAS, which confirms the resolution advantage of the proposed method.
Figure 3 shows images of the simulated anechoic cysts by different beamformers. It could be observed that noises are still at a high level inside the cysts of Fig. 3a, which mainly comes from the sidelobes of the background speckle region. The MV beamformer can slightly suppress the noises while the SLSC can offer an effective noise suppression. In the comparison between Fig. 3c, d, both show a good visualization of the anechoic cysts while the boundaries of the cysts in Fig. 3d are more distinct than in Fig. 3c. This indicates the advantage of integrating the MV beamformer with the SLSC algorithm. In Fig. 3d, noises are highly suppressed inside cysts, and the speckle regions are smoother with less dark spots. Generally speaking, the proposed method obtains a satisfactory performance in the cyst simulation.
We measure the CR, CNR, and speckle SNR for both cysts, and the numerical results are shown in Table 3. The corresponding speckle region marked by a black rectangle is selected at the same depth of the cyst marked by a white circle to avoid different attenuation through depths [11]. As shown in the results, Bmode MV images do not show much improvement in three metrics. The SLSC method could bring the significant increment mainly in the CNR. In comparison, the DAMV SLSC can increase three metrics and get a comparably better performance. Generally speaking, the noise inside cysts is remarkably suppressed by the proposed method, the CR has been notably increased. In addition, the speckle patterns are well preserved, which results in a higher speckle SNR.
Experimental phantoms
Figure 4 shows images of experimental phantoms consisting of two wire targets. The experimental results are similar to the simulated ones in “Experimental phantoms” section. The MV beamformer can obtain a narrower mainlobe width than the DAS beamformer. The speckle performance of the MV beamformer is also visually better than that of the DAS. For the SLSC beamformer, hyperechoic point targets are visually more detectable, which can be viewed from Fig. 4c, d. In comparison with the original SLSC method, the proposed method obtains better performance mainly in the mainlobe width.
Figure 5 plots the lateral variations across the depth z = 19 mm and 28 mm, while statistical results are given in Table 4. The MV gets the narrowest FWHM among all methods. In the SLSC series, the proposed method achieves a good performance just as it does in the point target simulation.
For the anechoic cyst phantom, results are shown in Fig. 6. The improvements of the proposed method mainly happen in less sidelobe noise inside cysts, clearer cyst margins and better speckle performance. Numerical results are presented in Table 5, which are calculated using the anechoic cyst marked by a white circle and its background speckles marked by a black rectangle. Similar to the anechoic cyst simulation, the MV shows no obvious differences in CR, CNR, and speckle SNR when compared with the DAS beamformer. With the SLSC method, the CNR is significantly increased, and the speckle SNR is also improved, which can be found in results of the SLSC and the DAMV SLSC. In the comparison between two SLSC methods, the DAMV SLSC obtains not only higher CR and CNR, but also better speckle SNR. This indicates that the proposed method is effective in suppressing background noises and enhancing the speckle performance.
In vivo study
In vivo images of a human carotid artery using different beamformers are presented in Fig. 7. The carotid artery has complicated anatomic structures, but the region inside the blood vessel has similar features with an anechoic cyst. Therefore, the same research process used in cyst simulations can be adopted to study the performance of these methods. As observed, noises inside the artery are obvious in Bmode images including the DAS one and the MV one. On the contrary, they are notably removed in the SLSC results. The artery wall can be better defined especially with the DAMV SLSC. The comparison between the SLSC and the DAMV SLSC demonstrates that the hyperechoic structures are more distinguishable with our method. The indistinct speckle regions in the original SLSC result are clearer with the proposed DAMV SLSC, as well, which improves the visualization of the anatomical structures.
Again for the quantitative assessment, the insideartery regions marked by a white circle and the background speckle marked by a black rectangle are adopted to calculate CRs, CNRs, and speckle SNRs. Corresponding results are given in Table 6. Compared with the original SLSC scheme, the proposed method brings an increased CR, along with a better CNR and speckle SNR.
Discussion
The proposed method is an optimization of the original SLSC method. The utilization of the MV beamforming in both the Tx and Rx aperture is effective in enhancing the imaging quality. From the simulated and experimental results, the proposed method obtains a comparably better performance than other beamformers. The resolution and contrast improvements result from two aspects. First, the implementation of the adaptive beamformer (MV) brings a better synthesis for both apertures than the conventional averaged compounding, which could further improve the performance of the following SLSC algorithm. The inherent adaptivity of the MV method mainly contributes to the resolution of the output images [20]. In addition, the MV helps remove the offaxis noise around hyperechoic targets, which results in a decline of the spatial coherence in these regions. Since the SLSC process sums all shortlag spatial coherence, the noise part in the final image is also suppressed. Second, the spatial coherence is measured in not only the Rx aperture but also the Tx one, which makes the proposed method different from others. The DAMV SLSC method is based on the acoustic reciprocity [16]. Thus, there is mutually dual between the beamforming processes in “Transmit aperture weighting for the SLSC imaging” and “Receive aperture weighting for the SLSC imaging” sections. The SLSC estimation in both apertures fully utilizes the coherence information in the synthetic aperture 2D echo data matrix. As a consequence, the contrast is further increased in comparison with the original SLSC scheme. In general, the high contrast along with a better resolution than conventional DAS and SLSC methods indicate that the proposed method could be a promising ultrasound imaging technique for better visualization of the anatomical structures and higher lesion detectability.
As for some parameters used in the proposed method, the subarray length L and the SLSC ratio Q are set to the same value for dual processes in order to accord with the acoustic reciprocity. Actually, using different values for the Tx and Rx aperture is also feasible since two processes in “Transmit aperture weighting for the SLSC imaging” and “Receive aperture weighting for the SLSC imaging” sections are mutually independent. Changing parameters for one aperture will not affect the beamforming process in another. It is a considerable issue that we can modify the final image by changing both apertures’ parameters separately. It might cause a degradation in the robustness but could bring advantages for specific situations. For instance, when the ROI contains several hyperechoic targets, enlarging L and decreasing Q for the Rx aperture can help achieve a higher contrast and better noise reduction around these targets, while the background speckle region could be damaged.
The computational complexity (CC) of the proposed method is a little higher than those of the MV beamformer and the original SLSC beamformer, since both the MV and the SLSC process are conducted in our method. The MV method has a complexity of \({\text{O}}\left( {L^{3} } \right)\) [25, 26], while the CC of the SLSC is comparably similar with that of the MV beamformer [2], also about \(O\left( {L^{3} } \right)\). The proposed scheme includes two MV beamforming processes and two SLSC estimations, for dual apertures, respectively. Thus, the CC is about four times the original MV methods, but still at the level of \(O\left( {L^{3} } \right)\). In consideration of the application of the graphics processing unit (GPU), it is possible to realize the realtime MV and SLSC algorithm [27,28,29]. We believe that our method also has the potential for realtime practical implementation.
There are still some unsolved issues in the proposed method. Above all, the integration of the MV and the SLSC could decrease the robustness in complicated situations. The hyperechoic cysts could be indistinguishable from the speckle region using the SLSC scheme. Another problem occurs in the clinical situation. The image enhancement of the proposed method could be influenced by tissue motions, channel noises, and phase aberration. Adopting a compensation technology for the tissue motion could be an effective alternative [30]. Future studies will focus on handling these drawbacks.
Conclusion
We develop a novel SLSC imaging method adopting the MV adaptive synthesis on dual apertures in this paper. The synthetic aperture echo dataset is adaptively beamformed in the Tx and Rx apertures, and the spatial coherence is also estimated through both apertures. Simulated, phantom and in vivo data are collected to demonstrate the advantages of the proposed method. Results show that our method can enhance the image resolution which is usually damaged by the original SLSC algorithm. More than this, the CR and speckle SNR are also increased due to the utilization of the adaptive focusing process. Therefore, we consider that the proposed method could be a promising SLSC imaging technique for the better imaging quality in medical ultrasound imaging.
Abbreviations
 CC:

computational complexity
 CF:

coherence factor
 CNR:

contrasttonoise ratio
 CR:

contrast ratio
 DAS:

delayandsum
 DAMV SLSC:

SLSC imaging using adaptive MV beamforming on dual apertures
 FWHM:

fullwidth at halfmaximum
 MV:

minimum variance
 PSF:

point spread function
 PWC:

plane wave compounding
 PWI:

plane wave imaging
 RF:

radio frequency
 ROI:

region of interest
 Rx:

receive
 SA:

synthetic aperture
 SLSC:

shortlag spatial coherence
 SNR:

signaltonoise ratio
 Tx:

transmit
References
 1.
Dahl JJ, Hyun D, Lediju M, Trahey GE. Lesion detectability in diagnostic ultrasound with shortlag spatial coherence imaging. Ultrason Imaging. 2011;33(2):119–33. https://doi.org/10.1177/016173461103300203.
 2.
Lediju MA, Trahey GE, Byram BC, Dahl JJ. ShortLag spatial coherence of backscattered echoes: imaging characteristics. IEEE Trans Ultrason Ferroelectr Freq Control. 2011;58(7):1377–88. https://doi.org/10.1109/tuffc.2011.1957.
 3.
Bell MAL, Dahl JJ, Trahey GE. Resolution and brightness characteristics of shortLag spatial coherence (SLSC) images. IEEE Trans Ultrason Ferroelectr Freq Control. 2015;62(7):1265–76. https://doi.org/10.1109/tuffc.2014.006909.
 4.
Bell MAL, Goswami R, Kisslo JA, Dahl JJ, Trahey GE. Shortlag spatial coherence imaging of cardiac ultrasound data: initial clinical results. Ultrasound Med Biol. 2013;39(10):1861–74. https://doi.org/10.1016/j.ultrasmedbio.2013.03.029.
 5.
Jakovljevic M, Trahey GE, Nelson RC, Dahl JJ. In vivo application of shortlag spatial coherence imaging in human liver. Ultrasound Med Biol. 2013;39(3):534–42. https://doi.org/10.1016/j.ultrasmedbio.2012.09.022.
 6.
Kakkad V, Dahl JJ, Ellestad S, Trahey G. In vivo application of shortlag spatial coherence and harmonic spatial coherence imaging in fetal ultrasound. Ultrason Imaging. 2014;37(2):101–16. https://doi.org/10.1177/0161734614547281.
 7.
Mallart R, Fink M. The Van CittertZernike theorem in pulse echo measurements. J Acoust Soc Am. 1991;90(5):2718–27. https://doi.org/10.1121/1.401867.
 8.
Walker WF, Trahey GE. Speckle coherence and implications for adaptive imaging. J Acoust Soc Am. 1997;101(4):1847–58. https://doi.org/10.1121/1.418235.
 9.
Bottenus N, Byram BC, Dahl JJ, Trahey GE. Synthetic aperture focusing for shortlag spatial coherence imaging. IEEE Trans Ultrason Ferroelectr Freq Control. 2013;60(9):1816–26. https://doi.org/10.1109/tuffc.2013.2768.
 10.
Chau G, Lavarello R, Dahl JJ. Shortlag spatial coherence weighted minimum variance beamformer for planewave images. In: IEEE International Ultrasonics Symposium, Tours, France, 18–21 September 2016:1–3. https://doi.org/10.1109/ultsym.2016.7728889.
 11.
Zhao J, Wang Y, Yu J, Guo W, Zhang S, Aliabadi S. Shortlag spatial coherence ultrasound imaging with adaptive synthetic transmit aperture focusing. Ultrason Imaging. 2017;39(4):224–39. https://doi.org/10.1177/0161734616688328.
 12.
Zhao J, Wang Y, Zeng X, Yu J, Yiu BYS, Yu ACH. Plane wave compounding based on a joint transmittingreceiving adaptive beamformer. IEEE Trans Ultrason Ferroelectr Freq Control. 2015;62(8):1440–52. https://doi.org/10.1109/tuffc.2014.006934.
 13.
Jensen JA, Nikolov SI, Gammelmark KL, Pedersen MH. Synthetic aperture ultrasound imaging. Ultrasonics. 2006;44:E5–15. https://doi.org/10.1016/j.ultras.2006.07.017.
 14.
Hemmsen MC, Rasmussen JH, Jensen JA. Tissue harmonic synthetic aperture ultrasound imaging. J Acoust Soc Am. 2014;136(4):2050–6. https://doi.org/10.1121/1.4893902.
 15.
Kortbek J, Jensen JA, Gammelmark KL. Sequential beamforming for synthetic aperture imaging. Ultrasonics. 2013;53(1):1–16. https://doi.org/10.1016/j.ultras.2012.06.006.
 16.
Bottenus N, Ustuner KF. Acoustic reciprocity of spatial coherence in ultrasound imaging. IEEE Trans Ultrason Ferroelectr Freq Control. 2015;62(5):852–61. https://doi.org/10.1109/tuffc.2014.006928.
 17.
Capon J. Highresolution frequencywavenumber spectrum analysis. Proc IEEE. 1969;57(8):1408–18. https://doi.org/10.1109/proc.1969.7278.
 18.
Shan TJ, Wax M, Kailath T. On spatial smoothing for directionofarrival estimation of coherent signals. IEEE Trans Acoust Speech Signal Process. 1985;33(4):806–11. https://doi.org/10.1109/tassp.1985.1164649.
 19.
Synnevag JF, Austeng A, Holm S. Adaptive beamforming applied to medical ultrasound imaging. IEEE Trans Ultrason Ferroelectr Freq Control. 2007;54(8):1606–13. https://doi.org/10.1109/tuffc.2007.431.
 20.
Synnevag JF, Austeng A, Holm S. Benefits of minimumvariance beamforming in medical ultrasound imaging. IEEE Trans Ultrason Ferroelectr Freq Control. 2009;56(9):1868–79. https://doi.org/10.1109/tuffc.2009.1263.
 21.
Jensen AC, Austeng A. An approach to multibeam covariance matrices for adaptive beamforming in ultrasonography. IEEE Trans Ultrason Ferroelectr Freq Control. 2012;59(6):1139–48. https://doi.org/10.1109/tuffc.2012.2304.
 22.
Rabinovich A, Friedman Z, Feuer A. Multiline acquisition with minimum variance beamforming in medical ultrasound imaging. IEEE Trans Ultrason Ferroelectr Freq Control. 2013;60(12):2521–31. https://doi.org/10.1109/tuffc.2013.2851.
 23.
Jensen JA, Svendsen NB. Calculation of pressure fields from arbitrarily shaped, apodized, and excited ultrasound transducers. IEEE Trans Ultrason Ferroelectr Freq Control. 1992;39(2):262–7. https://doi.org/10.1109/58.139123.
 24.
Jensen JA. Field: a program for simulating ultrasound systems. Med Biol Eng Comput. 1996;34(Suppl. 1, Pt. 1):351–3.
 25.
Zeng X, Chen C, Wang Y. Eigenspacebased minimum variance beamformer combined with Wiener postfilter for medical ultrasound imaging. Ultrasonics. 2012;52(8):996–1004. https://doi.org/10.1016/j.ultras.2012.07.012.
 26.
Golub GH, Van Loan CF. Matrix computations. Baltimore: Johns Hopkins University Press; 2012.
 27.
Asen JP, Buskenes JI, Colombo Nilsen IC, Austeng A, Holm S. Implementing Capon beamforming on a GPU for realtime cardiac ultrasound imaging. IEEE Trans Ultrason Ferroelectr Freq Control. 2014;61(1):76–85. https://doi.org/10.1109/tuffc.2014.6689777.
 28.
Yiu BYS, Yu ACH. GPUbased minimum variance beamformer for synthetic aperture imaging of the eye. Ultrasound Med Biol. 2015;41(3):871–83. https://doi.org/10.1016/j.ultrasmedbio.2014.11.005.
 29.
Hyun D, Trahey GE, Dahl JJ. In vivo demonstration of a realtime simultaneous Bmode/spatial coherence GPUbased beamformer. In: IEEE international ultrasonics symposium, Prague, Czech Republic, 21–25 July 2013:1280–3. https://doi.org/10.1109/ultsym.2013.0327.
 30.
Denarie B, Tangen TA, Ekroll IK, Rolim N, Torp H, Bjåstad T, Lovstakken L. Coherent plane wave compounding for very high frame rate ultrasonography of rapidly moving targets. IEEE Trans Med Imaging. 2013;32(7):1265–76. https://doi.org/10.1109/tmi.2013.2255310.
Authors’ contributions
Study concept and design (YQ); drafting of the manuscript (YQ); critical revision of the manuscript for important intellectual content (YQ, YW, JY, YG); obtained funding (YW, JY, YG); administrative, technical, and material support (YW); study supervision (YW). All authors read and approved the final manuscript.
Acknowledgements
This work was supported by the National Natural Science Foundation of China (61771143 and 81627804).
Competing interests
The authors declare that they have no competing interests.
Availability of data and materials
The simulated data and the experimental data used in this study are available from the corresponding author upon request after the publication of this article.
Consent for publication
Not applicable.
Ethics approval and consent to participate
Informed consent was obtained from the all individual participants included in the study.
Funding
This study was funded by the National Natural Science Foundation of China (61771143 and 81627804).
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Author information
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Received
Accepted
Published
DOI
Keywords
 Shortlag spatial coherence
 Synthetic aperture
 Minimum variance
 Adaptive beamforming
 Ultrasound imaging