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
Access Nature and 54 other Nature Portfolio journals
Get Nature+, our best-value online-access subscription
$29.99 / 30 days
cancel any time
Subscribe to this journal
Receive 12 digital issues and online access to articles
$99.00 per year
only $8.25 per issue
Buy this article
- Purchase on Springer Link
- Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout
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
Sanz, J. & Fayad, Z. A. Imaging of atherosclerotic cardiovascular disease. Nature 451, 953–957 (2008).
Forbes, J. M. & Cooper, M. E. Mechanisms of diabetic complications. Physiol. Rev. 93, 137–188 (2013).
Norgren, L. et al. Inter-society consensus for the management of peripheral arterial disease (TASC II). J. Vasc. Surg. 45, S5–S67 (2007).
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).
Davies, P. F. Hemodynamic shear stress and the endothelium in cardiovascular pathophysiology. Nat. Clin. Pract. Cardiovasc. Med. 6, 16–26 (2009).
Fox, M. D. & Raichle, M. E. Spontaneous fluctuations in brain activity observed with functional magnetic resonance imaging. Nat. Rev. Neurosci. 8, 700–711 (2007).
Maurovich-Horvat, P. et al. Comprehensive plaque assessment by coronary CT angiography. Nat. Rev. Cardiol. 11, 390–402 (2014).
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).
Fan, J. L. et al. High-speed volumetric two-photon fluorescence imaging of neurovascular dynamics. Nat. Commun. 11, 6020 (2020).
Shu, X., Beckmann, L. & Zhang, H. F. Visible-light optical coherence tomography: a review. J. Biomed. Opt. 22, 121707 (2017).
Wang, C. et al. Continuous monitoring of deep-tissue haemodynamics with stretchable ultrasonic phased arrays. Nat. Biomed. Eng. 5, 749–758 (2021).
Hu, H. et al. A wearable cardiac ultrasound imager. Nature 613, 667–675 (2023).
Wang, L. V. & Hu, S. Photoacoustic tomography: in vivo imaging from organelles to organs. Science 335, 1458–1462 (2012).
Wang, L. V. & Yao, J. A practical guide to photoacoustic tomography in the life sciences. Nat. Methods 13, 627–638 (2016).
Weber, J., Beard, P. C. & Bohndiek, S. E. Contrast agents for molecular photoacoustic imaging. Nat. Methods 13, 639–650 (2016).
Zhang, H. F., Maslov, K. & Wang, L. V. In vivo imaging of subcutaneous structures using functional photoacoustic microscopy. Nat. Protoc. 2, 797–804 (2007).
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).
Xu, M. & Wang, L. V. Universal back-projection algorithm for photoacoustic computed tomography. Phys. Rev. E 71, 016706 (2005).
Yao, J. et al. High-speed label-free functional photoacoustic microscopy of mouse brain in action. Nat. Methods 12, 407–410 (2015).
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).
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).
Lin, L. et al. Single-breath-hold photoacoustic computed tomography of the breast. Nat. Commun. 9, 2352 (2018).
Na, S. et al. Massively parallel functional photoacoustic computed tomography of the human brain. Nat. Biomed. Eng. 6, 584–592 (2022).
Özbek, A., Deán-Ben, X. L. & Razansky, D. Optoacoustic imaging at kilohertz volumetric frame rates. Optica 5, 857–863 (2018).
Wiskin, J. et al. Full wave 3D inverse scattering transmission ultrasound tomography in the presence of high contrast. Sci. Rep. 10, 20166 (2020).
Duarte, M. F. et al. Single-pixel imaging via compressive sampling. IEEE Signal Process. Mag. 25, 83–91 (2008).
Sun, B. et al. 3D computational imaging with single-pixel detectors. Science 340, 844–847 (2013).
Sun, M.-J. et al. Single-pixel three-dimensional imaging with time-based depth resolution. Nat. Commun. 7, 12010 (2016).
Stellinga, D. et al. Time-of-flight 3D imaging through multimode optical fibers. Science 374, 1395–1399 (2021).
Kruizinga, P. et al. Compressive 3D ultrasound imaging using a single sensor. Sci. Adv. 3, e1701423 (2017).
Luís Dean-Ben, X. & Razansky, D. Localization optoacoustic tomography. Light Sci. Appl. 7, 18004 (2018).
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).
Hahamovich, E. et al. Single-detector 3D optoacoustic tomography via coded spatial acoustic modulation. Commun. Eng. 1, 25 (2022).
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).
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).
Cox, B. T. & Beard, P. C. Photoacoustic tomography with a single detector in a reverberant cavity. J. Acoust. Soc. Am. 125, 1426–1436 (2009).
Brown, M. D. et al. Reverberant cavity photoacoustic imaging. Optica 6, 821–822 (2019).
Li, Y. et al. Snapshot photoacoustic topography through an ergodic relay for high-throughput imaging of optical absorption. Nat. Photonics 14, 164–170 (2020).
Li, Y. et al. Multifocal photoacoustic microscopy using a single-element ultrasonic transducer through an ergodic relay. Light Sci. Appl. 9, 135 (2020).
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).
Zhao, Y. & Wang, L. V. Single-shot photoacoustic imaging with single-element transducer through a spatiotemporal encoder. J. Biomed. Opt. 28, 046004 (2023).
Andersen, C. A. Noninvasive assessment of lower extremity hemodynamics in individuals with diabetes mellitus. J. Vasc. Surg. 52, 76S–80S (2010).
Sobieszczyk, P. & Beckman, J. Carotid artery disease. Circulation 114, e244–e247 (2006).
Vashist, S. K. et al. Emerging technologies for next-generation point-of-care testing. Trends Biotechnol. 33, 692–705 (2015).
Marik, P. E. & Baram, M. Noninvasive hemodynamic monitoring in the intensive care unit. Crit. Care Clin. 23, 383–400 (2007).
Moore, C. L. & Copel, J. A. Point-of-care ultrasonography. N. Engl. J. Med. 364, 749–757 (2011).
Trucco, A., Palmese, M. & Repetto, S. Devising an affordable sonar system for underwater 3-D vision. IEEE Trans. Instrum. Meas. 57, 2348–2354 (2008).
Reigber, A. et al. Very-high-resolution airborne synthetic aperture radar imaging: signal processing and applications. Proc. IEEE 101, 759–783 (2013).
Bernstein, E. F. Laser treatment of tattoos. Clin. Dermatol. 24, 43–55 (2006).
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).
Zhou, Q. et al. Piezoelectric single crystal ultrasonic transducers for biomedical applications. Prog. Mater. Sci. 66, 87–111 (2014).
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).
Pugsley, M. K. & Tabrizchi, R. The vascular system: an overview of structure and function. J. Pharmacol. Toxicol. Methods 44, 333–340 (2000).
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).
Loenneke, J. P. et al. Effect of cuff type on arterial occlusion. Clin. Physiol. Funct. Imaging 33, 325–327 (2013).
Zhou, Y., Liang, J. & Wang, L. V. Cuffing-based photoacoustic flowmetry in humans in the optical diffusive regime. J. Biophoton. 9, 208–212 (2016).
Sinex, J. E. Pulse oximetry: principles and limitations. Am. J. Emerg. Med. 17, 59–66 (1999).
Yang, J. et al. Photoacoustic assessment of hemodynamic changes in foot vessels. J. Biophoton. 12, e201900004 (2019).
Choi, W. et al. Three-dimensional multistructural quantitative photoacoustic and US imaging of human feet in vivo. Radiology 303, 467–473 (2022).
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).
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).
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).
Beck, A. & Teboulle, M. A fast iterative shrinkage-thresholding algorithm for linear inverse problems. SIAM J. Imaging Sci. 2, 183–202 (2009).
Zhu, Y. et al. Light emitting diodes based photoacoustic imaging and potential clinical applications. Sci. Rep. 8, 9885 (2018).
Fatima, A. et al. Review of cost reduction methods in photoacoustic computed tomography. Photoacoustics 15, 100137 (2019).
Ide, J. M. The velocity of sound in rocks and glasses as a function of temperature. J. Geol. 45, 689–716 (1937).
Yao, J. & Wang, L. V. Photoacoustic microscopy: photoacoustic microscopy. Laser Photon. Rev. 7, 758–778 (2013).
Wang, Y. et al. A robust and secure palm vessel biometric sensing system based on photoacoustics. IEEE Sens. J. 18, 5993–6000 (2018).
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).
Fink, M. & de Rosny, J. Time-reversed acoustics in random media and in chaotic cavities. Nonlinearity 15, R1–R18 (2002).
ANSI Z136.3-2018 American National Standard for Safe Use of Lasers in Health Care (American National Standards Institute, 2018).
Jose, J. et al. Speed-of-sound compensated photoacoustic tomography for accurate imaging. Med. Phys. 39, 7262–7271 (2012).
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
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
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
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.
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.
About this article
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
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41551-023-01149-4
This article is cited by
-
Fast capturing of deep blood flow
Nature Biomedical Engineering (2023)