Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

  • Article
  • Published:

Ultrafast longitudinal imaging of haemodynamics via single-shot volumetric photoacoustic tomography with a single-element detector

Abstract

Techniques for imaging haemodynamics use ionizing radiation or contrast agents or are limited by imaging depth (within approximately 1 mm), complex and expensive data-acquisition systems, or low imaging speeds, system complexity or cost. Here we show that ultrafast volumetric photoacoustic imaging of haemodynamics in the human body at up to 1 kHz can be achieved using a single laser pulse and a single element functioning as 6,400 virtual detectors. The technique, which does not require recalibration for different objects or during long-term operation, enables the longitudinal volumetric imaging of haemodynamics in vasculature a few millimetres below the skin’s surface. We demonstrate this technique in vessels in the feet of healthy human volunteers by capturing haemodynamic changes in response to vascular occlusion. Single-shot volumetric photoacoustic imaging using a single-element detector may facilitate the early detection and monitoring of peripheral vascular diseases and may be advantageous for use in biometrics and point-of-care testing.

This is a preview of subscription content, access via your institution

Access options

Buy this article

Prices may be subject to local taxes which are calculated during checkout

Fig. 1: The PACTER system.
Fig. 2: Single-shot 3D reconstruction in PACTER.
Fig. 3: Spatiotemporal characterization of PACTER.
Fig. 4: PACTER analysis of mouse haemodynamics in vivo.
Fig. 5: PACTER analysis of human hand haemodynamics in vivo.
Fig. 6: PACTER analysis of haemodynamic changes in human foot vessels.

Similar content being viewed by others

Data availability

The data supporting the findings of this study are provided within the Article and its Supplementary Information. The raw and analysed datasets generated during the study are available for research purposes from the corresponding authors on reasonable request.

Code availability

The reconstruction code, the system-control software and the data-collection software are proprietary and used in licensed technologies, yet they are available from the corresponding author on reasonable request.

References

  1. Sanz, J. & Fayad, Z. A. Imaging of atherosclerotic cardiovascular disease. Nature 451, 953–957 (2008).

    Article  CAS  PubMed  Google Scholar 

  2. Forbes, J. M. & Cooper, M. E. Mechanisms of diabetic complications. Physiol. Rev. 93, 137–188 (2013).

    Article  CAS  PubMed  Google Scholar 

  3. Norgren, L. et al. Inter-society consensus for the management of peripheral arterial disease (TASC II). J. Vasc. Surg. 45, S5–S67 (2007).

    Article  PubMed  Google Scholar 

  4. Kadem, M. et al. Hemodynamic modeling, medical imaging, and machine learning and their applications to cardiovascular interventions. IEEE Rev. Biomed. Eng. 16, 403–423 (2023).

  5. Davies, P. F. Hemodynamic shear stress and the endothelium in cardiovascular pathophysiology. Nat. Clin. Pract. Cardiovasc. Med. 6, 16–26 (2009).

    Article  CAS  PubMed  Google Scholar 

  6. Fox, M. D. & Raichle, M. E. Spontaneous fluctuations in brain activity observed with functional magnetic resonance imaging. Nat. Rev. Neurosci. 8, 700–711 (2007).

    Article  CAS  PubMed  Google Scholar 

  7. Maurovich-Horvat, P. et al. Comprehensive plaque assessment by coronary CT angiography. Nat. Rev. Cardiol. 11, 390–402 (2014).

  8. Provost, J. et al. Simultaneous positron emission tomography and ultrafast ultrasound for hybrid molecular, anatomical and functional imaging. Nat. Biomed. Eng. 2, 85–94 (2018).

    Article  CAS  PubMed  Google Scholar 

  9. Fan, J. L. et al. High-speed volumetric two-photon fluorescence imaging of neurovascular dynamics. Nat. Commun. 11, 6020 (2020).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Shu, X., Beckmann, L. & Zhang, H. F. Visible-light optical coherence tomography: a review. J. Biomed. Opt. 22, 121707 (2017).

  11. Wang, C. et al. Continuous monitoring of deep-tissue haemodynamics with stretchable ultrasonic phased arrays. Nat. Biomed. Eng. 5, 749–758 (2021).

    Article  CAS  PubMed  Google Scholar 

  12. Hu, H. et al. A wearable cardiac ultrasound imager. Nature 613, 667–675 (2023).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Wang, L. V. & Hu, S. Photoacoustic tomography: in vivo imaging from organelles to organs. Science 335, 1458–1462 (2012).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Wang, L. V. & Yao, J. A practical guide to photoacoustic tomography in the life sciences. Nat. Methods 13, 627–638 (2016).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Weber, J., Beard, P. C. & Bohndiek, S. E. Contrast agents for molecular photoacoustic imaging. Nat. Methods 13, 639–650 (2016).

    Article  CAS  PubMed  Google Scholar 

  16. Zhang, H. F., Maslov, K. & Wang, L. V. In vivo imaging of subcutaneous structures using functional photoacoustic microscopy. Nat. Protoc. 2, 797–804 (2007).

    Article  CAS  PubMed  Google Scholar 

  17. Li, L. et al. Single-impulse panoramic photoacoustic computed tomography of small-animal whole-body dynamics at high spatiotemporal resolution. Nat. Biomed. Eng. 1, 0071 (2017).

    Article  PubMed  PubMed Central  Google Scholar 

  18. Xu, M. & Wang, L. V. Universal back-projection algorithm for photoacoustic computed tomography. Phys. Rev. E 71, 016706 (2005).

    Article  Google Scholar 

  19. Yao, J. et al. High-speed label-free functional photoacoustic microscopy of mouse brain in action. Nat. Methods 12, 407–410 (2015).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Cao, R. et al. Label-free intraoperative histology of bone tissue via deep-learning-assisted ultraviolet photoacoustic microscopy. Nat. Biomed. Eng. 7, 124–134 (2022).

  21. Shi, J. et al. High-resolution, high-contrast mid-infrared imaging of fresh biological samples with ultraviolet-localized photoacoustic microscopy. Nat. Photon. 13, 609–615 (2019).

    Article  CAS  Google Scholar 

  22. Lin, L. et al. Single-breath-hold photoacoustic computed tomography of the breast. Nat. Commun. 9, 2352 (2018).

    Article  PubMed  PubMed Central  Google Scholar 

  23. Na, S. et al. Massively parallel functional photoacoustic computed tomography of the human brain. Nat. Biomed. Eng. 6, 584–592 (2022).

    Article  CAS  PubMed  Google Scholar 

  24. Özbek, A., Deán-Ben, X. L. & Razansky, D. Optoacoustic imaging at kilohertz volumetric frame rates. Optica 5, 857–863 (2018).

    Article  PubMed  PubMed Central  Google Scholar 

  25. Wiskin, J. et al. Full wave 3D inverse scattering transmission ultrasound tomography in the presence of high contrast. Sci. Rep. 10, 20166 (2020).

  26. Duarte, M. F. et al. Single-pixel imaging via compressive sampling. IEEE Signal Process. Mag. 25, 83–91 (2008).

    Article  Google Scholar 

  27. Sun, B. et al. 3D computational imaging with single-pixel detectors. Science 340, 844–847 (2013).

    Article  CAS  PubMed  Google Scholar 

  28. Sun, M.-J. et al. Single-pixel three-dimensional imaging with time-based depth resolution. Nat. Commun. 7, 12010 (2016).

    Article  PubMed  PubMed Central  Google Scholar 

  29. Stellinga, D. et al. Time-of-flight 3D imaging through multimode optical fibers. Science 374, 1395–1399 (2021).

    Article  CAS  PubMed  Google Scholar 

  30. Kruizinga, P. et al. Compressive 3D ultrasound imaging using a single sensor. Sci. Adv. 3, e1701423 (2017).

    Article  PubMed  PubMed Central  Google Scholar 

  31. Luís Dean-Ben, X. & Razansky, D. Localization optoacoustic tomography. Light Sci. Appl. 7, 18004 (2018).

    Article  PubMed  PubMed Central  Google Scholar 

  32. Deán-Ben, X. L., Özbek, A., López-Schier, H. & Razansky, D. Acoustic scattering mediated single detector optoacoustic tomography. Phys. Rev. Lett. 123, 174301 (2019).

    Article  PubMed  Google Scholar 

  33. Hahamovich, E. et al. Single-detector 3D optoacoustic tomography via coded spatial acoustic modulation. Commun. Eng. 1, 25 (2022).

  34. Montaldo, G., Palacio, D., Tanter, M. & Fink, M. Time reversal kaleidoscope: a smart transducer for three-dimensional ultrasonic imaging. Appl. Phys. Lett. 84, 3879–3881 (2004).

    Article  CAS  Google Scholar 

  35. Montaldo, G., Palacio, D., Tanter, M. & Fink, M. Building three-dimensional images using a time-reversal chaotic cavity. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 52, 1489–1497 (2005).

    Article  PubMed  Google Scholar 

  36. Cox, B. T. & Beard, P. C. Photoacoustic tomography with a single detector in a reverberant cavity. J. Acoust. Soc. Am. 125, 1426–1436 (2009).

    Article  CAS  PubMed  Google Scholar 

  37. Brown, M. D. et al. Reverberant cavity photoacoustic imaging. Optica 6, 821–822 (2019).

  38. Li, Y. et al. Snapshot photoacoustic topography through an ergodic relay for high-throughput imaging of optical absorption. Nat. Photonics 14, 164–170 (2020).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Li, Y. et al. Multifocal photoacoustic microscopy using a single-element ultrasonic transducer through an ergodic relay. Light Sci. Appl. 9, 135 (2020).

  40. Li, L., Li, Y., Zhang, Y. & Wang, L. V. Snapshot photoacoustic topography through an ergodic relay of optical absorption in vivo. Nat. Protoc. 14, 164–170 (2021).

    Google Scholar 

  41. Zhao, Y. & Wang, L. V. Single-shot photoacoustic imaging with single-element transducer through a spatiotemporal encoder. J. Biomed. Opt. 28, 046004 (2023).

    Article  PubMed  PubMed Central  Google Scholar 

  42. Andersen, C. A. Noninvasive assessment of lower extremity hemodynamics in individuals with diabetes mellitus. J. Vasc. Surg. 52, 76S–80S (2010).

    Article  PubMed  Google Scholar 

  43. Sobieszczyk, P. & Beckman, J. Carotid artery disease. Circulation 114, e244–e247 (2006).

    Article  PubMed  Google Scholar 

  44. Vashist, S. K. et al. Emerging technologies for next-generation point-of-care testing. Trends Biotechnol. 33, 692–705 (2015).

  45. Marik, P. E. & Baram, M. Noninvasive hemodynamic monitoring in the intensive care unit. Crit. Care Clin. 23, 383–400 (2007).

    Article  CAS  PubMed  Google Scholar 

  46. Moore, C. L. & Copel, J. A. Point-of-care ultrasonography. N. Engl. J. Med. 364, 749–757 (2011).

  47. Trucco, A., Palmese, M. & Repetto, S. Devising an affordable sonar system for underwater 3-D vision. IEEE Trans. Instrum. Meas. 57, 2348–2354 (2008).

    Article  Google Scholar 

  48. Reigber, A. et al. Very-high-resolution airborne synthetic aperture radar imaging: signal processing and applications. Proc. IEEE 101, 759–783 (2013).

    Article  Google Scholar 

  49. Bernstein, E. F. Laser treatment of tattoos. Clin. Dermatol. 24, 43–55 (2006).

    Article  PubMed  Google Scholar 

  50. Zhou, Q., Lau, S., Wu, D. & Kirk Shung, K. Piezoelectric films for high frequency ultrasonic transducers in biomedical applications. Prog. Mater. Sci. 56, 139–174 (2011).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Zhou, Q. et al. Piezoelectric single crystal ultrasonic transducers for biomedical applications. Prog. Mater. Sci. 66, 87–111 (2014).

  52. Zhou, Q. et al. PMN-PT single crystal, high-frequency ultrasonic needle transducers for pulsed-wave Doppler application. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 54, 668–675 (2007).

    Article  PubMed  Google Scholar 

  53. Pugsley, M. K. & Tabrizchi, R. The vascular system: an overview of structure and function. J. Pharmacol. Toxicol. Methods 44, 333–340 (2000).

  54. Langham, M. C. et al. Evaluation of cuff-induced ischemia in the lower extremity by magnetic resonance oximetry. J. Am. Coll. Cardiol. 55, 598–606 (2010).

  55. Loenneke, J. P. et al. Effect of cuff type on arterial occlusion. Clin. Physiol. Funct. Imaging 33, 325–327 (2013).

    Article  PubMed  Google Scholar 

  56. Zhou, Y., Liang, J. & Wang, L. V. Cuffing-based photoacoustic flowmetry in humans in the optical diffusive regime. J. Biophoton. 9, 208–212 (2016).

    Article  Google Scholar 

  57. Sinex, J. E. Pulse oximetry: principles and limitations. Am. J. Emerg. Med. 17, 59–66 (1999).

    Article  CAS  PubMed  Google Scholar 

  58. Yang, J. et al. Photoacoustic assessment of hemodynamic changes in foot vessels. J. Biophoton. 12, e201900004 (2019).

    Article  Google Scholar 

  59. Choi, W. et al. Three-dimensional multistructural quantitative photoacoustic and US imaging of human feet in vivo. Radiology 303, 467–473 (2022).

    Article  PubMed  Google Scholar 

  60. Kragelj, R. et al. Parameters of postocclusive reactive hyperemia measured by near infrared spectroscopy in patients with peripheral vascular disease and in healthy volunteers. Ann. Biomed. Eng. 29, 311–320 (2001).

  61. de Mul, F. F. M., Morales, F., Smit, A. J. & Graaff, R. A model for post-occlusive reactive hyperemia as measured with laser-Doppler perfusion monitoring. IEEE Trans. Biomed. Eng. 52, 184–190 (2005).

    Article  PubMed  Google Scholar 

  62. Bioucas-Dias, J. M. & Figueiredo, M. A. T. A new TwIST: two-step iterative shrinkage/thresholding algorithms for image restoration. IEEE Trans. Image Process. 16, 2992–3004 (2007).

    Article  PubMed  Google Scholar 

  63. Beck, A. & Teboulle, M. A fast iterative shrinkage-thresholding algorithm for linear inverse problems. SIAM J. Imaging Sci. 2, 183–202 (2009).

    Article  Google Scholar 

  64. Zhu, Y. et al. Light emitting diodes based photoacoustic imaging and potential clinical applications. Sci. Rep. 8, 9885 (2018).

    Article  PubMed  PubMed Central  Google Scholar 

  65. Fatima, A. et al. Review of cost reduction methods in photoacoustic computed tomography. Photoacoustics 15, 100137 (2019).

    Article  PubMed  PubMed Central  Google Scholar 

  66. Ide, J. M. The velocity of sound in rocks and glasses as a function of temperature. J. Geol. 45, 689–716 (1937).

    Article  Google Scholar 

  67. Yao, J. & Wang, L. V. Photoacoustic microscopy: photoacoustic microscopy. Laser Photon. Rev. 7, 758–778 (2013).

    Article  Google Scholar 

  68. Wang, Y. et al. A robust and secure palm vessel biometric sensing system based on photoacoustics. IEEE Sens. J. 18, 5993–6000 (2018).

    Article  Google Scholar 

  69. Garrett, D. C., Xu, J., Ku, G. & Wang, L. V. Whole-body human ultrasound tomography. Preprint at https://doi.org/10.48550/arXiv.2307.00110 (2023).

  70. Fink, M. & de Rosny, J. Time-reversed acoustics in random media and in chaotic cavities. Nonlinearity 15, R1–R18 (2002).

    Article  Google Scholar 

  71. ANSI Z136.3-2018 American National Standard for Safe Use of Lasers in Health Care (American National Standards Institute, 2018).

  72. Jose, J. et al. Speed-of-sound compensated photoacoustic tomography for accurate imaging. Med. Phys. 39, 7262–7271 (2012).

    Article  PubMed  Google Scholar 

  73. Jerman, T., Pernuš, F., Likar, B. & Špiclin, Ž. Beyond Frangi: an improved multiscale vesselness filter. In Proc. SPIE Medical Imaging (eds Ourselin, S. & Styner, M. A.) 94132A (SPIE, 2015); https://doi.org/10.1117/12.2081147

Download references

Acknowledgements

We thank Y. Zhao for contributing to the universal calibration. This work was supported in part by National Institutes of Health grants R01 EB028277, U01 EB029823 and R35 CA220436 (Outstanding Investigator Award). The computations presented here were conducted at the Resnick High Performance Computing Center, a facility supported by the Resnick Sustainability Institute at the California Institute of Technology.

Author information

Authors and Affiliations

Authors

Contributions

Y. Zhang and L.V.W. conceived and designed the study. Y. Zhang, L.L., R.C. and K.M. built the imaging system. Y. Zhang developed the data acquisition program. P.H. developed the 3D reconstruction algorithm. Y. Zhang, L.L., R.C. and A.K. performed the experiments. Y. Zhang, P.H. and X.T. processed and analysed the data. Y. Zeng, L.J. and Q.Z. fabricated the ultrasonic transducer. L.V.W. supervised the study. All of the authors contributed to writing the paper.

Corresponding author

Correspondence to Lihong V. Wang.

Ethics declarations

Competing interests

L.V.W. has a financial interest in Microphotoacoustics, CalPACT and Union Photoacoustic Technologies. These companies did not provide support for this work. K.M. has a financial interest in Microphotoacoustics. All other authors declare no competing interests.

Peer review

Peer review information

Nature Biomedical Engineering thanks Xose Luis Dean-Ben, Wenfeng Xia and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

Additional information

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Extended data

Extended Data Fig. 1 PACTER of mouse haemodynamics in vivo.

a, 3D PACTER image of the abdominal vasculature of mouse 3. Norm., normalized. b, Cross-sectional 2D images corresponding to the yellow rectangle in a at four different time instances from the 4D PACTER datasets. t0 = 0.49 s. The white solid curve represents a two-term Gaussian fit of the vessels’ profile denoted by the yellow dashed line. Differences from the first image are highlighted in colour. c, PA amplitudes along the yellow dashed line (1D image) in b versus time, where the time instances in b are labelled with vertical grey lines. d, Blue solid and orange dash-dotted curves represent the centre positions and widths of the vessel on the left (based on the first term of the Gaussian fit in b) versus time. Shaded areas denote the standard deviations (n = 5). e, Fourier transforms of the centre positions and widths of the vessel in d, showing the respiratory frequency from the vessel widths only. Scale bars, 1 mm.

Extended Data Fig. 2 PACTER of human hand vasculature in vivo.

a,c–e, 3D PACTER images of the vasculature in a middle finger (a), an index finger (c), a proximal phalanx region (d), and a thenar region (e). Norm., normalized. b, Photograph of a human hand showing the imaged regions. Scale bars, 1 mm.

Extended Data Fig. 3 PACTER of human hand haemodynamics in vivo.

a, 3D PACTER image of the thenar vasculature of participant 1, in a region different from that in Fig. 5b. Norm., normalized. b, Maximum amplitude projections of the 3D volumes from the 4D PACTER datasets along the z axis in a at the time instances before, during, and after cuffing. t0 = 15.34 s. The solid lines flank the vessel under investigation. Differences from the t0 image are highlighted in colour. c, PA amplitudes along the vessel (1D image) in b versus time, where the time instances in b are labelled with vertical grey lines. The blue and orange arrows indicate peak responses in the occlusion and recovery phases, respectively. d, Positions of the blood front along the blood vessel during the occlusion (left) and recovery (right) phases in c. The left blue curve is an exponential fit with an occlusion rate of 0.6 ± 0.1 m/s, whereas the right orange curve is a linear fit showing the blood flow speed of 21.6 ± 7.9 m/s. e, Comparison between the durations of the occlusion and recovery phases in d. P < 0.001, calculated by the two-sample t-test. Error bars, means ± standard errors of the means (n = 9). Scale bars, 1 mm.

Extended Data Fig. 4 PACTER of haemodynamic changes in human foot vessels.

a, 3D PACTER images of the foot vessels of participant 3, in a region different from that in Fig. 6b, before (left) and after (right) vascular occlusion. Norm., normalized. b, Maximum amplitude projections of the 3D volumes from the 4D PACTER datasets along the z axis in a. t0 = 0 s. The blue and orange circles indicate regions of a vein and an artery, respectively. c, Difference between the two images in b. d, Relative PA signals of the venous and arterial regions indicated by the blue and orange circles in b and c. The shaded areas denote the standard deviations (n = 5). The arrows and vertical lines indicate the start times of occlusion and recovery. Scale bars, 1 mm.

Supplementary information

Supplementary Information

Supplementary notes, figures, tables and video captions.

Reporting Summary

Peer Review File

Supplementary Video 1

Principle and implementation of PACTER (with narration).

Supplementary Video 2

Simulation of PATER and PACTER signals.

Supplementary Video 3

4D PACTER image and its maximum amplitude projection of bovine blood flushing through an S-shaped tube.

Supplementary Video 4

4D PACTER images of bovine blood flowing through a tube at different speeds.

Supplementary Video 5

4D PACTER image and its maximum amplitude projection of bovine blood flowing through a tube with a speed of 272.5 mm s−1, captured at 1,000 volumes per second.

Supplementary Video 6

4D in vivo PACTER image of the abdominal vasculature of mouse 1.

Supplementary Video 7

4D in vivo PACTER image of the abdominal vasculature of mouse 2.

Supplementary Video 8

4D in vivo PACTER image of the abdominal vasculature of mouse 3.

Supplementary Video 9

4D in vivo PACTER image of the thenar vasculature of participant 1.

Supplementary Video 10

4D in vivo PACTER image of the thenar vasculature of participant 2.

Supplementary Video 11

4D in vivo PACTER image of the thenar vasculature of participant 1 in a different region.

Supplementary Video 12

4D in vivo PACTER image of the foot vessels of participant 3.

Supplementary Video 13

4D in vivo PACTER image of the foot vessels of participant 3 in a different region.

Supplementary Video 14

PACTER signals affected by temperature fluctuations.

Rights and permissions

Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, Y., Hu, P., Li, L. et al. Ultrafast longitudinal imaging of haemodynamics via single-shot volumetric photoacoustic tomography with a single-element detector. Nat. Biomed. Eng (2023). https://doi.org/10.1038/s41551-023-01149-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1038/s41551-023-01149-4

This article is cited by

Search

Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing