Photoacoustic imaging explained

Photoacoustic imaging

Photoacoustic imaging or optoacoustic imaging is a biomedical imaging modality based on the photoacoustic effect. Non-ionizing laser pulses are delivered into biological tissues and part of the energy will be absorbed and converted into heat, leading to transient thermoelastic expansion and thus wideband (i.e. MHz) ultrasonic emission. The generated ultrasonic waves are detected by ultrasonic transducers and then analyzed to produce images. It is known that optical absorption is closely associated with physiological properties, such as hemoglobin concentration and oxygen saturation.[1] As a result, the magnitude of the ultrasonic emission (i.e. photoacoustic signal), which is proportional to the local energy deposition, reveals physiologically specific optical absorption contrast. 2D or 3D images of the targeted areas can then be formed.[2]

Biomedical imaging

See main article: Biomedical imaging. The optical absorption in biological tissues can be due to endogenous molecules such as hemoglobin or melanin, or exogenously delivered contrast agents. As an example, Fig. 2 shows the optical absorption spectra of oxygenated hemoglobin (HbO2) and deoxygenated hemoglobin (Hb) in the visible and near infrared region.[3] Since blood usually has orders of magnitude higher absorption than surrounding tissues, there is sufficient endogenous contrast for photoacoustic imaging to visualize blood vessels. Recent studies have shown that photoacoustic imaging can be used in vivo for tumor angiogenesis monitoring, blood oxygenation mapping, functional brain imaging, skin melanoma detection, methemoglobin measuring, etc.[2]

Δf Primary contrast Δz δz δx Speed
Hz mm μm μm Mvx/s
Photoacoustic microscopy50 MOptical absorption315450.5
Photoacoustic tomography5 MOptical absorption507007000.5
Confocal microscopyFluorescence, scattering0.23-200.3-310-100
Two-photon microscopyFluorescence0.5-1.01-100.3-310-100
Optical coherence tomography300 TOptical scattering1-20.5-101-1020-4.000
Scanning laser acoustic microscopy300 MUltrasonic scattering1-2202010
Acoustic microscopy50 MUltrasonic scattering2020-10080-1600.1
Ultrasonography5 MUltrasonic scattering603003001
Table 1. Comparison of contrast mechanisms, penetration depth (Δz), axial resolution (δz), lateral resolution (δx=δy) and imaging speed of confocal microscopy, two-photon microscopy, optical coherence tomography (300 THz), ultrasound microscopy (50 MHz), ultrasound imaging (5 MHz), photoacoustic microscopy (50 MHz), and photoacoustic tomography (3.5 MHz). Speeds in megavoxel per second of non-parallel techniques.

Two types of photoacoustic imaging systems, photoacoustic/thermoacoustic computed tomography (also known as photoacoustic/thermoacoustic tomography, i.e., PAT/TAT) and photoacoustic microscopy (PAM), have been developed. A typical PAT system uses an unfocused ultrasound detector to acquire the photoacoustic signals, and the image is reconstructed by inversely solving the photoacoustic equations. A PAM system, on the other hand, uses a spherically focused ultrasound detector with 2D point-by-point scanning, and requires no reconstruction algorithm.

Photoacoustic computed tomography

General equation

Given the heating function

H(\vec{r},t)

, the generation and propagation of photoacoustic wave pressure

p(\vec{r},t)

in an acoustically homogeneous inviscid medium is governed by
2p(\vec{r},t)-1
2
v
s
\nabla
\partial2
\partial{t2
}p(\vec,t)=-\frac\fracH(\vec,t) \qquad \qquad \quad \quad (1),

where

vs

is the speed of sound in medium,

\beta

is the thermal expansion coefficient, and

Cp

is the specific heat capacity at constant pressure. Eq. (1) holds under thermal confinement to ensure that heat conduction is negligible during the laser pulse excitation. The thermal confinement occurs when the laser pulsewidth is much shorter than the thermal relaxation time.[4]

The forward solution of Eq. (1) is given by

\left.p(\vec{r},t)=\beta
4\piCp

\int

d\vec{r'
}
\frac \right|_ \qquad \quad \,\,\,\,(2).

In stress confinement, which occurs when the laser pulsewidth is much shorter than the stress relaxation time,[4] Eq. (2) can be further derived as

p(\vec{r},t)=1
4\pi
2
v
s
\partial
\partialt

\left[

1
vst

\intd\vec{r'}p0(\vec{r'})\delta\left(t-

|\vec{r
-\vec{r'}|}{v

s}\right)\right]    (3),

where

p0

is the initial photoacoustic pressure.

Universal reconstruction algorithm

In a PAT system, the acoustic pressure is detected by scanning an ultrasonic transducer over a surface that encloses the photoacoustic source. To reconstruct the internal source distribution, we need to solve the inverse problem of equation (3) (i.e. to obtain

p0

). A representative method applied for PAT reconstruction is known as the universal backprojection algorithm.[5] This method is suitable for three imaging geometries: planar, spherical, and cylindrical surfaces.

The universal back projection formula is

where

\Omega0

is the solid angle subtended by the entire surface

S0

with respect to the reconstruction point

\vec{r}

inside

S0

, and

Simple system

A simple PAT/TAT/OAT system is shown in the left part of Fig. 3. The laser beam is expanded and diffused to cover the whole region of interest. Photoacoustic waves are generated proportional to the distribution of optical absorption in the target, and are detected by a single scanned ultrasonic transducer. A TAT/OAT system is the same as PAT except that it uses a microwave excitation source instead of a laser. Although single-element transducers have been employed in these two systems, the detection scheme can be extended to use ultrasound arrays as well.

Biomedical applications

Intrinsic optical or microwave absorption contrast and diffraction-limited high spatial resolution of ultrasound make PAT and TAT promising imaging modalities for wide biomedical applications:

Brain lesion detection

Soft tissues with different optical absorption properties in the brain can be clearly identified by PAT.[6]

Hemodynamics monitoring

Since HbO2 and Hb are the dominant absorbing compounds in biological tissues in the visible spectral range, multiple wavelength photoacoustic measurements can be used to reveal the relative concentration of these two chromophores.[6] [7] Thus, the relative total concentration of hemoglobin (HbT) and the hemoglobin oxygen saturation (sO2) can be derived. Therefore, cerebral hemodynamic changes associated with brain function can be successfully detected with PAT.

Breast cancer diagnosis

By utilizing low scattered microwave for excitation, TAT is capable of penetrating thick (several cm) biological tissues with less than mm spatial resolution.[8] Since cancerous tissue and normal tissue have about the same responses to radio frequency radiation, TAT has limited potential in early breast cancer diagnosis.

Photoacoustic microscopy

See main article: Photoacoustic microscopy. The imaging depth of photoacoustic microscopy is mainly limited by the ultrasonic attenuation. The spatial (i.e. axial and lateral) resolutions depend on the ultrasonic transducer used. An ultrasonic transducer with high central frequency and broader bandwidth are chosen to obtain high axial resolution. The lateral resolution is determined by the focal diameter of the transducer. For instance, a 50 MHz ultrasonic transducer provides 15 micrometre axial and 45 micrometre lateral resolution with ~3 mm imaging depth.

Photoacoustic microscopy has multiple important applications in functional imaging: it can detect changes in oxygenated/deoxygenated hemoglobin in small vessels.[9] [10]

Other applications

Photoacoustic imaging was introduced recently in the context of artwork diagnostics with emphasis on the uncovering of hidden features such as underdrawings or original sketch lines in paintings. Photoacoustic images, collected from miniature oil paintings on canvas, illuminated with a pulsed laser on their reverse side, revealed clearly the presence of pencil sketch lines coated over by several paint layers.[11]

Advances in photoacoustic imaging

Photoacoustic imaging has seen recent advances through the integration of deep learning principles and compressed sensing. For more information on the deep learning applications in photoacoustic imaging, see Deep learning in photoacoustic imaging.

See also

External links

Notes and References

  1. A. Grinvald. Functional architecture of cortex revealed by optical imaging of intrinsic signals . Nature . 324. 361–364 . 1986 . 10.1038/324361a0 . 3785405 . 6095. etal. 1986Natur.324..361G. 4328958 .
  2. M. Xu . L.H. Wang . Photoacoustic imaging in biomedicine . Review of Scientific Instruments . 77 . 4. 041101–041101–22 . 2006 . 10.1063/1.2195024. 2006RScI...77d1101X.
  3. http://omlc.ogi.edu/spectra/ Optical Properties Spectra
  4. Book: L.H. Wang . H.I. Wu . Biomedical Optics. Wiley . 2007. 978-0-471-74304-0.
  5. M. Xu. Universal back-projection algorithm for photoacoustic-computed tomography. Physical Review E. 71. 016706. 2005. 10.1103/PhysRevE.71.016706 . 15697763. 1. etal. 2005PhRvE..71a6706X. 1969.1/180492. free.
  6. X. Wang. Non-invasive laser-induced photoacoustic tomography for structural and functional imaging of the brain in vivo. Nature Biotechnology. 21. 803–806. 2003. 10.1038/nbt839 . 12808463 . 7. 2961096. etal.
  7. X. Wang. Non-invasive imaging of hemoglobin concentration and oxygenation in the rat brain using high-resolution photoacoustic tomography. Journal of Biomedical Optics. 11. 024015. 2006. 10.1117/1.2192804 . 16674205 . 2. etal. 2006JBO....11b4015W. 9488754 .
  8. G. Ku. Thermoacoustic and photoacoustic tomography of thick biological tissues toward breast imaging. Technology in Cancer Research and Treatment. 4. 559–566. 2005 . 5 . 16173826. etal. 10.1177/153303460500400509. 1969.1/181686. 15782118. free.
  9. Yao. Junjie. Wang. Lihong V.. 2013-01-31. Photoacoustic microscopy. Laser & Photonics Reviews. 7. 5. 758–778. 10.1002/lpor.201200060. 24416085. 3887369. 2013LPRv....7..758Y. 1863-8880.
  10. Zhang. Hao F. Maslov. Konstantin. Stoica. George. Wang. Lihong V. 2006-06-25. Functional photoacoustic microscopy for high-resolution and noninvasive in vivo imaging. Nature Biotechnology. 24. 7. 848–851. 10.1038/nbt1220. 16823374. 912509. 1087-0156.
  11. Tserevelakis. George J.. Vrouvaki. Ilianna. Siozos. Panagiotis. Melessanaki. Krystallia. Hatzigiannakis. Kostas. Fotakis. Costas. Zacharakis. Giannis. 2017-04-07. Photoacoustic imaging reveals hidden underdrawings in paintings. Scientific Reports. En. 7. 1. 747. 10.1038/s41598-017-00873-7. 28389668. 5429688. 2045-2322. 2017NatSR...7..747T.