3-D Super-Resolution Ultrasound (SR-US) Imaging using a 2-D Sparse Array with High Volumetric Imaging Rate

3-D Super-Resolution Ultrasound (SR-US) Imaging
using a 2-D Sparse Array
with High Volumetric Imaging Rate

Sevan Harput, Kirsten Christensen-Jeffries, Jemma Brown, Jiaqi Zhu, Ge Zhang, Chee Hau Leow,
Matthieu Toulemonde, Alessandro Ramalli, Enrico Boni, Piero Tortoli,
Robert J. Eckersley, Chris Dunsby, and Meng-Xing Tang

ULIS Group, Department of Bioengineering, Imperial College London, London, SW7 2BP, UK
Biomedical Engineering Department, Division of Imaging Sciences, King’s College London, SE1 7EH, London, UK
Department of Information Engineering, University of Florence, 50139 Florence, IT
Lab. on Cardiovascular Imaging & Dynamics, Dept. of Cardiovascular Sciences, KU Leuven, Leuven, Belgium
Department of Physics and the Centre for Pathology, Imperial College London, London, SW7 2AZ, UK
These authors contributed equally to this work
E-mail: S.Harput@imperial.ac.uk, Mengxing.Tang@imperial.ac.uk

Super-resolution ultrasound imaging has been so far achieved in 3-D by mechanically scanning a volume with a linear probe, by co-aligning multiple linear probes, by using multiplexed 3-D clinical ultrasound systems, or by using 3-D ultrasound research systems. In this study, a 2-D sparse array was designed with 512 elements according to a density-tapered 2-D spiral layout and optimized to reduce the sidelobes of the transmitted beam profile. High frame rate volumetric imaging with compounded plane waves was performed using two synchronized ULA-OP256 systems. Localization-based 3-D super-resolution images of two touching sub-wavelength tubes were generated from a 120 second acquisition.

I Introduction

By localizing spatially isolated microbubbles through multiple frames, super-resolution ultrasound (SR-US) images can be generated. Even at ultrasonic frequencies in the low MHz range, a localization precision of a few micrometers can be achieved [1, 2]. In the absence of sample motion, it is this localization precision that determines the maximum achievable resolution in super-resolution images. If motion is present and subsequently corrected post-acquisition, then the motion correction accuracy can also limit the achievable spatial resolution [3, 4]. Researchers demonstrated the use of 2-D SR-US imaging in many studies using microbubbles [5, 6, 7, 8, 9, 10, 11, 12, 13] and nanodroplets [14, 15, 16]. All of these studies were based on imaging with 1-D probes that can only display 2-D slices of a 3-D structure, thus making the volumetric observations more challenging. In addition to this, 2-D ultrasound imaging introduces more limitations in localization based super-resolution imaging. First, super-resolution cannot be achieved in the elevational direction. Second, out-of-plane motion can only be compensated for movements smaller than the elevational beamwidth of the transducer. However, with the implementation of 3-D super-resolution ultrasound imaging, diffraction limited resolution can be overcome in every direction and 3-D motion can be tracked over larger scales.

Although 3-D super-resolution imaging has not been achieved with a 2-D imaging probe with a high volumetric rate, several studies have extended super-resolution into the third dimension. O’Reilly and Hynynen used a subset of a hemispherical transcranial therapy array to generate 3-D super-resolution images of a spiral tube phantom through an ex vivo human skullcap [17]. Desailly et al. used 2 parallel series of transducers to image microfluidic channels and obtained 3-D super-localization by fitting parallel parabolas in the elevation direction [18]. Errico et al. performed a coronal scan of an entire rat brain and generated 2-D super-resolution slices by using a micro-step motor and a linear array [19]. Lin et al. mechanically scanned a rat FSA tumor using a linear array mounted on a motorized precision motion stage to generate 3-D super-resolution images by calculating the maximum intensity projection from all 2-D slices [20]. Christensen-Jeffries et al. generated volumetric 3-D super-resolution images of cellulose tubes at the overlapping imaging region of two orthogonal linear arrays at the focus [21].

Fig. 1: Simulated ultrasound field propagation from the sparse array computed at 6, 12 and 18 mm depths.

Using a full 2-D array for 3-D SR-US imaging requires a system with very large number of independent channels, the design of which might be impractical due to the high cost, complexity, and large data size. Several technologies have been developed to use a large number of active elements with fewer scanner channels in order to simplify ultrasound systems and probes. In this study, a density-tapered sparse array method was chosen instead of a full 2-D array to reduce the number of channels and hence the amount of data while maintaining the frame rate. This approach is similar to previous studies on minimally redundant 2-D arrays [22] and sparse 2-D arrays [23, 24, 25, 26, 27], but uses more number of elements to improve transmit power and receive sensitivity. Our method significantly differs from row-column addressing and multiplexing approaches since it maintains simultaneous access to all probe elements through independent channels. The sparse array was designed specifically for high volumetric rate 3-D super-resolution ultrasound imaging based on a density-tapered spiral layout [28].

Ii Materials and Methods

Ii-a 2-D Sparse Array

Fig. 2: Layout of the 2-D sparse array. Empty rows are introduced due to manufacturing limitations and not related to the density-tapered 2-D spiral method.

A sparse array was designed by selecting 512 elements from a element 2-D matrix array as shown in Fig. 2. Row number 9, 18, and 27 are intentionally left blank for wiring. The method to select the location of sparse array elements is based on the density-tapered 2-D spiral layout. This method arranges the elements according to Fermat’s spiral seeds with spatial density modulation and reduces the sidelobes of the transmitted beam profile. This deterministic, aperiodic, and balanced positioning procedure guarantees uniform performance over a wide range of imaging angles [28].

To verify the feasibility of the sparse array for 3-D imaging with plane waves, Field II simulations were performed [29, 30]. [29]. Ultrasound propagation from the sparse array is simulated at different depths as shown in Fig. 1. Closer to the transducer the radiated plane wave has a tail, which makes it not suitable for plane wave imaging. However, far from the transducer, the ultrasound field becomes uniform and the residual field behind the wavefront is diminished, which is suitable for plane wave imaging.

After the in silico feasibility test, the 2-D sparse array (Vermon S.A., Tours, France) was specifically fabricated for 3-D SR-US with an element size of m, center frequency of 3.7 MHz and a bandwidth of 60%.

Ii-B Experimental Setup

Two ULA-OP256 [31, 32] systems were synchronized to transmit 9 plane waves steered within a range of degrees in the lateral and elevational directions from the 512 selected elements. Each compounded volume consisted of 9 volumetric datasets acquired in 3.6 ms with a pulse repetition frequency of 2500 Hz. This compounded plane wave imaging scheme was used to measure the super-localization precision of the system before the experiments, which was found to be 18 m.

Two 200 m tubes arranged in a double helix shape were imaged during the experiments as shown in Fig. 3. A 1:1000 diluted Sonovue (Bracco S.p.A, Milan, Italy) solution was flowed through both tubes with a constant flow rate that produced an average microbubble velocity of 10 mm/s. A total of 3000 volumetric ultrasound frames were acquired in 120 seconds. Beamforming and volumetric reconstruction were performed offline. SVD was used to separate the microbubble signals from the echoes originating from the tube and the assembly holding the tube. Localization of isolated microbubbles was performed on every acquired volume to generate the 3-D super-resolved volumes [33].

Fast data acquisition was used with a PRF of 2500 Hz to ensure minimum artefacts on compounded volumes due to moving microbubbles. Before the next fast acquisition, there was a 40 ms interval to allow the microbubbles to travel approximately half the size of the B-mode PSF. This short pause between acquisitions was introduced to reduce the data redundancy and maximize the total acquisition duration for the allocated memory size.

Fig. 3: Optical Images of two 200 m cellulose tubes arranged in a double helix pattern.

Iii Results & Discussion

Fig. 4 shows the 3-D visualization of microvessels overlaid on 2-D maximum-intensity-projection (MIP) slices, where it is not possible to visualize two 200 m tubes. The full-width-half-maximum (FWHM) of the 3-D B-mode point-spread-function (PSF) was measured by using linear interpolation as 793, 772, and 499 m in X, Y, Z directions respectively [34].

Fig. 5 shows the 3-D SR-US images of the sub-wavelength structures, where the imaging wavelength is 404 m in water at 25C. A total of 2319 microbubbles were localized within the 3000 volumes after compounding. Due to the large number of localizations, the 3-D structure of the tubes cannot be clearly visualized in a single 2-D image. To improve the visualization, 3-D SR-US images are plotted with depth information color-coded in the image from different viewing angles. The FWHM of a single tube, which was measured over a 2 mm length, appeared as 193 m at the widest point and 81% of the super-localizations were within a diameter of 200 m.

Volumetric imaging with 2-D sparse arrays can be a reasonable choice instead of using a full matrix array or a multiplexed array. In comparison to full matrix arrays, the proposed solution is less expensive; it can achieve the same volumetric acquisition speed since all elements of 2-D spiral array are continuously connected to every channel in the system; and it has a comparable spatial resolution since the bandwidths and aperture sizes are almost the same. In comparison to 2-D multiplexed arrays, it has the flexibility to modify the number of compoundings depending on the application. When the acquired ultrasound volume data is a combination of multiple transmissions and acquisitions, intra-volume motion artefacts will be generated. This can be a problem for multiplexed 3-D ultrasound systems for some applications.

Iv Conclusion

The main limitation of localization-based SR-US imaging performed in 2-D is the lack of super-resolution in the elevation direction. In this study, this issue was addressed by using a bespoke 2-D sparse array, which achieved super-resolution in the elevational direction and sufficient SNR for 3-D SR-US imaging. Two 200 m, smaller than half wavelength, tubes arranged in a double helix shape were resolved in 3-D using the 2-D spiral array probe.

Fig. 4: Isosurface of the 3-D ultrasound image is plotted in white at -15, -20, and -25 dB levels with degrading transparency. 2-D maximum intensity projections with a 30 dB dynamic range is overlaid on the volumetric image.
Fig. 5: Two images zoom show the structure from different viewing angles. Depth-encoded colorscale is added to improve the visualization of two tubes arranged in a double helix shape. The corresponding optical images are presented as an inset.


This work was supported mainly by the EPSRC under Grant EP/N015487/1 and EP/N014855/1, in part by the King’s College London (KCL) and Imperial College London EPSRC Centre for Doctoral Training in Medical Imaging (EP/L015226/1), in part by the Wellcome EPSRC Centre for Medical Engineering at KCL (WT 203148/Z/16/Z), in part by the Department of Health through the National Institute for Health Research comprehensive Biomedical Research Center Award to Guy’s and St Thomas’ NHS Foundation Trust in partnership with KCL and King’s College Hospital NHS Foundation Trust, in part by the Graham-Dixon Foundation and in part by NVIDIA GPU grant.


  • [1] O. M. Viessmann, R. J. Eckersley, K. Christensen-Jeffries, M.-X. Tang, and C. Dunsby, “Acoustic super-resolution with ultrasound and microbubbles,” Phys. Med. Biol, vol. 58, pp. 6447–6458, 2013.
  • [2] Y. Desailly, J. Pierre, O. Couture, and M. Tanter, “Resolution limits of ultrafast ultrasound localization microscopy,” Phys. Med. Biol, vol. 60, pp. 8723–8740, 2015.
  • [3] S. Harput, K. Christensen-Jeffries, Y. Li, J. Brown, R. J. Eckersley, C. Dunsby, and M.-X. Tang, “Two stage sub-wavelength motion correction in human microvasculature for ceus imaging,” in IEEE International Ultrasonics Symposium (IUS), 2017, pp. 1–4.
  • [4] S. Harput, K. Christensen-Jeffries, J. Brown, Y. Li, K. J. Williams, A. H. Davies, R. J. Eckersley, C. Dunsby, and M. Tang, “Two-stage motion correction for super-resolution ultrasound imaging in human lower limb,” IEEE Trans. Ultrason., Ferroelectr., Freq. Control, vol. 65, no. 5, pp. 803–814, 2018.
  • [5] K. Christensen-Jeffries, R. J. Browning, M.-X. Tang, C. Dunsby, and R. J. Eckersley, “In vivo acoustic super-resolution and super-resolved velocity mapping using microbubbles,” IEEE Trans. Med. Imag., vol. 34, no. 2, pp. 433–440, 2015.
  • [6] D. Ackermann and G. Schmitz, “Detection and tracking of multiple microbubbles in ultrasound b-mode images,” IEEE Trans. Ultrason., Ferroelectr., Freq. Control, vol. 63, no. 1, pp. 72–82, 2016.
  • [7] A. Bar-Zion, C. Tremblay-Darveau, O. Solomon, D. Adam, and Y. C. Eldar, “Fast vascular ultrasound imaging with enhanced spatial resolution and background rejection,” IEEE Trans. Med. Imag., vol. 36, pp. 169–180, 2017.
  • [8] J. Foiret, H. Zhang, T. Ilovitsh, L. Mahakian, S. Tam, and K. W. Ferrara, “Ultrasound localization microscopy to image and assess microvasculature in a rat kidney,” Scientific Reports, vol. 7, no. 13662, pp. 1–12, 2017.
  • [9] S. Harput, K. Christensen-Jeffries, J. Brown, R. J. Eckersley, C. Dunsby, and M.-X. Tang, “Localisation of multiple non-isolated microbubbles with frequency decomposition in super-resolution imaging,” in IEEE International Ultrasonics Symposium (IUS), 2017, pp. 1–4.
  • [10] O. Couture, V. Hingot, B. Heiles, P. Muleki-Seya, and M. Tanter, “Ultrasound localization microscopy and super-resolution: A state of the art,” IEEE Trans. Ultrason., Ferroelectr., Freq. Control, vol. 65, no. 8, pp. 1304–1320, 2018.
  • [11] P. Song, J. D. Trzasko, A. Manduca, R. Huang, R. Kadirvel, D. F. Kallmes, and S. Chen, “Improved super-resolution ultrasound microvessel imaging with spatiotemporal nonlocal means filtering and bipartite graph-based microbubble tracking,” IEEE Trans. Ultrason., Ferroelectr., Freq. Control, vol. 65, no. 2, pp. 149–167, 2018.
  • [12] T. Opacic, S. Dencks, B. Theek, M. Piepenbrock, D. Ackermann, A. Rix, T. Lammers, E. Stickeler, S. Delorme, G. Schmitz, and F. Kiessling, “Motion model ultrasound localization microscopy for preclinical and clinical multiparametric tumor characterization,” Nature Communications, vol. 9, no. 1527, pp. 1–13, 2018.
  • [13] T. Ilovitsh, A. Ilovitsh, J. Foiret, B. Z. Fite, and K. W. Ferrara, “Acoustical structured illumination for super-resolution ultrasound imaging,” Communications Biology, vol. 1, no. 3, 2018.
  • [14] G. P. Luke, A. S. Hannah, and S. Y. Emelianov, “Super-resolution ultrasound imaging in vivo with transient laser-activated nanodroplets,” Nano Letters, vol. 16, pp. 2556–2559, 2016.
  • [15] H. Yoon, K. A. Hallam, C. Yoon, and S. Y. Emelianov, “Super-resolution imaging with ultrafast ultrasound imaging of optically triggered perfluorohexane nanodroplets,” IEEE Trans. Ultrason., Ferroelectr., Freq. Control, vol. 65, no. pre-print, pp. 1–1, 2018.
  • [16] G. Zhang, S. Harput, S. Lin, K. Christensen-Jeffries, C. H. Leow, J. Brown, C. Dunsby, R. J. Eckersley, and M.-X. Tang, “Acoustic wave sparsely activated localization microscopy (awsalm): Super-resolution ultrasound imaging using acoustic activation and deactivation of nanodroplets,” Applied Physics Letters, vol. 113, no. 1, p. 014101, 2018.
  • [17] M. A. O’Reilly and K. Hynynen, “A super-resolution ultrasound method for brain vascular mapping,” Medical Physics, vol. 40, no. 110701, 2013.
  • [18] Y. Desailly, O. Couture, M. Fink, and M. Tanter, “Sono-activated ultrasound localization microscopy,” Applied Physics Letters, vol. 103, no. 174107, 2013.
  • [19] C. Errico, J. Pierre, S. Pezet, Y. Desailly, Z. Lenkei, O. Couture, and M. Tanter, “Ultrafast ultrasound localization microscopy for deep super-resolution vascular imaging,” Nature, vol. 527, pp. 499–507, 2015.
  • [20] F. Lin, S. E. Shelton, D. Espindola, J. D. Rojas, G. Pinton, and P. A. Dayton, “3-d ultrasound localization microscopy for identifying microvascular morphology features of tumor angiogenesis at a resolution beyond the diffraction limit of conventional ultrasound,” Theranostics, vol. 7, no. 1, pp. 196–204, 2017.
  • [21] K. Christensen-Jeffries, J. Brown, P. Aljabar, M.-X. Tang, C. Dunsby, and R. J. Eckersley, “3-d in vitro acoustic super-resolution and super-resolved velocity mapping using microbubbles,” IEEE Trans. Ultrason., Ferroelectr., Freq. Control, vol. 64, no. 10, pp. 1478–1486, 2017.
  • [22] M. Karaman, I. O. Wygant, O. Oralkan, and B. T. Khuri-Yakub, “Minimally redundant 2-d array designs for 3-d medical ultrasound imaging,” IEEE Trans. Med. Imag., vol. 28, no. 7, pp. 1051–1061, 2009.
  • [23] A. Austeng and S. Holm, “Sparse 2-d arrays for 3-d phased array imaging - design methods,” IEEE Trans. Ultrason., Ferroelectr., Freq. Control, vol. 49, no. 8, pp. 1073–1086, 2002.
  • [24] B. Diarra, M. Robini, P. Tortoli, C. Cachard, and H. Liebgott, “Design of optimal 2-d nongrid sparse arrays for medical ultrasound,” IEEE Trans. Biomed. Eng, vol. 60, no. 11, pp. 3093–3102, 2013.
  • [25] E. Roux, A. Ramalli, P. Tortoli, C. Cachard, M. Robini, and H. Liebgott, “2-d ultrasound sparse arrays multidepth radiation optimization using simulated annealing and spiral-array inspired energy functions,” IEEE Trans. Ultrason., Ferroelectr., Freq. Control, vol. 63, no. 12, pp. 2138–2149, 2016.
  • [26] E. Roux, A. Ramalli, H. Liebgott, C. Cachard, M. Robini, and P. Tortoli, “Wideband 2-d array design optimization with fabrication constraints for 3-d us imaging,” IEEE Trans. Ultrason., Ferroelectr., Freq. Control, vol. 64, no. 1, pp. 108–125, 2017.
  • [27] E. Roux, F. Varray, L. Petrusca, C. Cachard, P. Tortoli, and H. Liebgott, “Experimental 3-d ultrasound imaging with 2-d sparse arrays using focused and diverging waves,” Scientific Reports, vol. 8, no. 9108, pp. 1–12, 2018.
  • [28] A. Ramalli, E. Boni, A. S. Savoia, and P. Tortoli, “Density-tapered spiral arrays for ultrasound 3-d imaging,” IEEE Trans. Ultrason., Ferroelectr., Freq. Control, vol. 62, no. 8, pp. 1580–1588, 2015.
  • [29] J. Jensen and N. B. Svendsen, “Calculation of pressure fields from arbitrarily shaped, apodized, and excited ultrasound transducers,” IEEE Trans. Ultrason., Ferroelectr., Freq. Control, vol. 39, pp. 262–267, 1992.
  • [30] J. Jensen, “Field: A program for simulating ultrasound systems,” in Medical & Biological Engineering & Computing, vol. 34, no. 1, 1996, pp. 351–353.
  • [31] E. Boni, L. Bassi, A. Dallai, F. Guidi, V. Meacci, A. Ramalli, S. Ricci, and P. Tortoli, “Ula-op 256: A 256-channel open scanner for development and real-time implementation of new ultrasound methods,” IEEE Trans. Ultrason., Ferroelectr., Freq. Control, vol. 63, no. 10, pp. 1488–1495, 2016.
  • [32] E. Boni, L. Bassi, A. Dallai, V. Meacci, A. Ramalli, M. Scaringella, F. Guidi, S. Ricci, and P. Tortoli, “Architecture of an ultrasound system for continuous real-time high frame rate imaging,” IEEE Trans. Ultrason., Ferroelectr., Freq. Control, vol. 64, no. 9, pp. 1276–1284, 2017.
  • [33] K. Christensen-Jeffries, S. Harput, J. Brown, P. N. T. Wells, P. Aljabar, C. Dunsby, M.-X. Tang, and R. J. Eckersley, “Microbubble axial localization errors in ultrasonic super-resolution imaging,” IEEE Trans. Ultrason., Ferroelectr., Freq. Control, vol. 64, no. 11, pp. 1644–1654, 2017.
  • [34] S. Harput, J. McLaughlan, D. M. Cowell, and S. Freear, “New performance metrics for ultrasound pulse compression systems,” in IEEE International Ultrasonics Symposium (IUS), 2014, pp. 440–443.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description