Multi-spectral photoacoustic elasticity tomography

Yubin Liu; Zhen Yuan

doi:10.1364/BOE.7.003323

1. Introduction

Photoacoustic tomography (PAT) is a robust biomedical imaging method that can offer the structural and functional information of biological tissues with excellent resolution and high contrast [1–7]. PAT imaging techniques have been successfully applied to the early detection of cancer, probing the brain function and examining vascular and skin diseases [8–11]. More importantly, recent work also shows that PAT can reconstruct the tissue mechanical properties including the acoustic velocity, the elastic modulus and the temperature [12–15], and optical and physiological properties such as the optical absorption and scattering coefficients, the deoxyhemoglobin (HbR) and oxyhemoglobin (HbO₂) concentrations by using spectrally resolved photoacoustic (PA) measurements [6,16–19]. In particular, the uniqueness regarding simultaneous reconstruction of different chromophore concentrations and acoustic velocity has been resolved for multi-spectral PAT [18], which can reveal spatially resolved quantitative physiological and molecular information by exploiting the known spectral characteristics of specific chromophores. The generated physiological properties of biological tissues including the concentrations of HbR and HbO₂ and water (H2O) content are essential for an accurate diagnostic decision-making in lesion detection and image-guided cancer treatment.

Interestingly, the bulk elastic modulus of normal and diseased tissues have been explored by ultrasonography, where the ranges of the elastic moduli from soft tissues span over as much as four orders of the magnitudes [20,21]. Recently, PAT has been implemented to characterize the elastic moduli of phantoms or biological tissues using single-wavelength measurements. However, most of work conducted can only reconstruct the elastic moduli-related parameters, which cannot achieve the truly quantitative and spectrally resolved PA elasticity tomography (PAET) [14]. In this study, multi-spectral PAET is proposed to directly reconstruct the chromophore concentrations and quantitative elastic modulus of biological tissues simultaneously. The capability of reconstructing the physiology and elastic properties by using multi-spectral PAET, paves a new avenue for better differentiating benign from malignant lesions since the elastic contrast between diseased and healthy tissues is very high [21].

2. Methods and materials

2.1 The multi-spectral PAET reconstruction method

To formulate the reconstruction model for multi-spectral PAET, the PA wave generation and propagation in an acoustic coupling media are described by the basic Newton’s law of motion Eq. (1), equation of continuity (2) and the thermal elastic Eq. (3) in terms of thermal confinement [1],

ρ (r) \frac{\partial}{\partial t} V (r, t) = - \nabla p (r, t) .

\nabla \cdot V (r, t) = - \frac{1}{ρ (r) v_{s}^{2} (r)} \frac{\partial}{\partial t} p (r, t) + β \frac{\partial}{\partial t} T (r, t) .

ρ (r) C_{p} \frac{\partial}{\partial t} T (r, t) = H (r, t) .

where V is particle velocity, v_s is the acoustic velocity, t is time, r is the spatial location,

ρ

is the density of the media, T and p are the temperature and pressure, H is the source term which can be written as

H = Ψ I (t)

, C_p is the specific heat,

Ψ

is the light absorbed energy density, and I(t) is the temporal illumination function. Considering Eqs. (1)-(3) and eliminating V, we get

ρ (r) \nabla \cdot (\frac{1}{ρ (r)} \nabla p (r, t)) - \frac{1}{v_{s}^{2} (r)} \frac{\partial^{2}}{\partial t^{2}} p (r, t) = - \frac{β}{C_{p}} \frac{\partial}{\partial t} H (r, t) .

If the bulk elastic modulus K is denoted as

K = ρ v_{s}^{2}

, Eq. (4) is written,

\nabla (\frac{1}{ρ (r)} \nabla p (r, t)) - \frac{1}{K} \frac{\partial^{2}}{\partial t^{2}} p (r, t) = - \frac{β}{ρ C_{p}} \frac{\partial}{\partial t} H (r, t) .

If a homogeneous elastic reference medium is assumed with density

ρ = ρ_{0}

, Eq. (5) is written,

\nabla^{2} p (r, t) - \frac{ρ_{0}}{K} \frac{\partial^{2}}{\partial t^{2}} p (r, t) = - \frac{β}{C_{p}} \frac{\partial}{\partial t} H (r, t) .

Denoting the following Fourier transform form for acoustic pressure,

P (r, ω) = \int_{- \infty}^{+ \infty} p (r, t) \exp (- i ω t) d t .

And taking the Fourier transform on variable t of Eq. (7), we obtain,

\nabla^{2} p (r, ω) + \frac{K_{0}}{K} k_{0}^{2} p (r, ω) = i \frac{k_{0} β}{C_{p}} v_{0} Ψ (r) .

where

k_{0} = ω / v_{0}

is the wave number described by the angular frequency

ω

and the speed of acoustic wave v₀ in a reference medium, K₀ is the bulk elastic modulus of the reference medium and the light absorbed energy density

Ψ

is equal to the product of the optical absorption coefficient

μ_{a}

and photon density

Φ

. Equation (8) can be further written as

\nabla^{2} p (r, ω) + O k_{0}^{2} p (r, ω) = i \frac{k_{0} β}{C_{p}} v_{0} Ψ (r) .

where O is a coefficient that depends on the bulk modulus K and

O = K_{0} / K

.

In multi-spectral PAT, the frequency-domain Helmholtz wave equation is expresses as

\nabla^{2} p (r, ω, λ) + O k_{0}^{2} p (r, ω, λ) = i \frac{k_{0} β}{C_{p}} v_{0} μ (r, ω) Φ (r, λ)

where λ is the wavelength of the incident light. According to the Beer’s law, the wavelength-dependent tissue absorption coefficient is written [16,18],

μ_{a} (λ) = \sum_{i = 1}^{} ε_{i} (λ) c_{i} .

in which

c_{i}

is the concentration and

ε_{i} (λ)

is the extinction coefficient of the ith chromophore at wavelength λ . Consequently, Eq. (10) is rewritten,

\nabla^{2} p (r, ω, λ) + O k_{0}^{2} p (r, ω, λ) = i \frac{k_{0} β}{C_{p}} v_{0} \sum_{i = 1}^{} ε_{i} (λ) c_{i} Φ (r, λ) .

Equation (12) is the forward equation for multi-spectral PAET and the finite element discretization form of Eq. (12) is [16,18],

A p = B

The inverse solution can be obtained by solving the following inverse equation:

(J^{T} J + ξ I) Δ χ = J^{T} (p^{o} - p^{c}) .

in which

Δ χ = {(\begin{matrix} Δ c_{l, 1} & Δ c_{l, 2} & ... & Δ c_{l, n} & Δ O_{1} & Δ O_{2} & ... & Δ O_{n} \end{matrix})}^{T}

is the update vector for l chromophores, elasticity modulus-related parameter O and N is the node number of the finite element mesh;

ξ

is the regularization parameter determined by combined Marquardt and Tikhonov regularization schemes [16];

p_{i}^{o}

and

p_{i}^{c}

are measured and calculated data for i = 1,2,…M boundary locations and are written as for each acoustic frequency ω within each incident optical wavelength λ,

p^{o} = {(p_{1}^{o}, p_{2}^{o}, \dots, p_{M}^{o})}^{T} | (ω, λ), p^{c} = {(p_{1}^{c}, p_{2}^{c}, \dots, p_{M}^{c})}^{T} | (ω, λ) .

The Jacobian matrix, J is denoted as

J = [{\tilde{J}}_{O, λ, ω}, {\tilde{J}}_{c_{i}, λ, ω}]

, where

{\tilde{J}}_{O, λ, ω}

and

{\tilde{J}}_{c_{i}, λ, ω}

represents the Jacobian submatrix for the elastic modulus and different chromophores, respectively [15]. In consideration of

Ψ = μ_{a} Φ

, the elements in Jacobian matrix

{\tilde{J}}_{c_{i}, λ, ω}

are determined by

\frac{\partial p_{i} (ω, λ)}{\partial {(c_{l})}_{, j}} = \sum_{k = 1}^{N} (\frac{\partial p_{i} (ω, λ)}{\partial Ψ_{k} (λ)} \frac{\partial Ψ_{k} (λ)}{\partial μ_{a, j} (λ)} ε_{l} (λ) \begin{matrix} ) & (i = 1, 2... M; j = 1, 2... N; l = 1, 2... C h r o m) \end{matrix} .

in which Chrom is the number of chromophores and the derivatives

\partial Ψ / \partial μ_{a}

in Eq. (3) are further denoted as

\frac{\partial Ψ_{k} (λ)}{\partial μ_{a}_{j} (λ)} = (\begin{matrix} Φ_{j} + μ_{a, j} (\partial Φ_{j} / \partial μ_{a, j}) & i f \begin{matrix} \begin{matrix} (k = j) \end{matrix} \end{matrix} \\ μ_{a, k} (\partial Φ_{k} / \partial μ_{a, j}) & i f \begin{matrix} (k \neq j) \end{matrix} \end{matrix} | (λ) .

The sensitivity of $\partial p (ω, λ) / \partial Ψ (λ)$ and $\partial Φ / \partial μ_{a}$ in each iteration can be calculated from Eq. (10) and the following photon diffusion equation for each wavelength, respectively,

\nabla \cdot D (r, λ) \nabla Φ (r, λ) - μ_{a} (r, λ) Φ (r, λ) = - S (r, λ) .

in which S(r) is the light source term and D(r) is the optical diffusion coefficient and assumed as a constant in this study.

Likewise, the elements in Jacobian submatrix ${\tilde{J}}_{O, λ, ω}$ are determined by [22],

[A] {\frac{\partial p (ω, λ)}{\partial O_{j}}} = - [\frac{\partial A}{\partial O_{j}}] {p (ω, λ)} .

The Jacobian matrix can be calculated through the following steps using the adjoint method:

First, we define a MxN matrix $Ψ$ , and let $Ψ$ satisfy the following relationship:

[Ψ] [A] = [Δ_{d}] .

where the vector

Δ_{d}

has the unit value at the measurement sites/nodes and zero at other nodes. Then we left multiply Eq. (19) by

Ψ

, which yields

[Ψ] [A] {\frac{\partial p (λ)}{\partial O_{j}}} = - [Ψ] [\frac{\partial A}{\partial O_{j}}] {p (λ)} .

Inspecting Eq. (20) into Eq. (21), we can immediately find that the left hand side of Eq. (21) actually gives the corresponding elements in the Jacobian submatrix

{\tilde{J}}_{O, λ, ω}

.

As such, in multi-spectral PAET, the image formation task is to update chromophores concentration and bulk elastic modulus distributions via the iterative solution of Eqs. (13), (14) and (18) so that a weighted sum of the squared difference between computed and measured acoustic pressures can be minimized.

2.2 The multi-spectral PAET imaging systems

For our home-made PAT imaging system at the Faculty of Health Sciences of the University of Macau, a pulsed light from an Nd:YAG laser with OPO based on the multi-wavelength excitation (wavelength range from 680 to 1064 nm; pulse duration: 5-10 ns; frequency rate: 20 Hz; Surelite I-10, Continuum) was used to illuminate the phantom/biological tissues via an optical subsystem and generate acoustic signals. A transducer (1MHz central frequency; OLYMPUS NDT) and the phantom were immersed into a water tank. A rotary stage rotated the transducer relative to the center of the tank as shown in Fig. 1. The incident optical fluence was controlled at 10 mJ/cm², and the diameter of the incident laser beam was 2.0 cm. The complex wavefield signal was amplified by a Pulser/Receiver (5073R, OLYMPUS).

Fig. 1 Schematic of the PAT system. BS: beam splitter; PC: personal computer.

Download Full Size | PDF

2.3 Simulation and in vitro experimental tests

For simulation test 1, a circular background with the diameter of 3cm contained three circular targets (2mm in radius each as shown in Fig. 2(a)), where each target had different contrasts in physiological properties (Hb, HbO₂, and H₂O) and bulk elastic modulus (K). The chromophore concentrations and K used were provided in Table 1 and six optical wavelengths (633, 670, 723, 805, 854and 896nm) were utilized to generate the measurements for the multi-spectral PAET. A total of 120 detectors were equally distributed along the boundary of the circular background. The extinction coefficient for each chromophore was adopted from the website at http://omlc.ogi.edu/spectra/index.html [23].

Fig. 2 (a) Test geometry. The dimension units are in millimeters. (b) The exact (left column) and recovered (right column) images using 6 optical wavelengths. The first to the third rows of Fig. 2(b) denote the Hb, HBO₂ and H₂O, respectively. The axes (left and bottom) illustrate the spatial scale, in mm, whereas the color scale (right) denote HbO₂ and HbR in μM, and H₂O in percentage.

Download Full Size | PDF

Table 1. Chromophore Concentrations and K for the targets in simulation test 1

View Table | View all tables in this article

The test geometry for the second simulation was the same as simulation test 1, in which we assumed the chromophore of the targets only consisted of HbR. Two inclusions (top right and middle bottom targets) had different HbR contrast levels while two inclusions (top left and middle bottom targets) had the same contrast level for K. The HbR concentrations and K used for the second simulation test were given in Table 2. The measurements from two wavelengths (670nm and 896nm) were used for the multi-wavelength PAET reconstruction and the extinction coefficients were also taken from the website at http://omlc.ogi.edu/spectra/index.html [23]. In this study, the constant Gruneisen function $Γ = c_{0}^{2} β / C_{p}$ ( $c_{0}$ = 1495 m/s), ${μ^{'}}_{s}$ = 1.0 mm^-1, and K₀ = 2.23GPa were used for all the tests.

Table 2. HbR concentration and K of the targets in simulation test 2

View Table | View all tables in this article

In addition, for the first in vitro experiment, we embedded a square chicken breast (length: 3mm; width:2mm) into a 3.5-cm-diameter solid cylindrical phantom. We then immersed the chicken breast-bearing solid phantom into the water tank. Two optical wavelengths (532 and 680nm) were used for the measurements. The phantom materials used consisted of Intralipid as scatterer and India ink as absorber with Agar powder (1-2%) for solidifying the Intralipid and India ink solution.

For the second in vitro experiment, we placed a phantom that contained the porcine liver tissue into the water. The solid cylindrical phantom with the diameter of 3.5cm consisted of Intralipid as the scatterer and India ink as the absorber with Agar powder (1-2%) for solidifying the Intralipid and India ink solution. In this test, we used the dual-wavelength of 700nm and 750nm for the measurements and the molar extinction coefficients of porcine liver were adopted from reference [24].

3. Results and discussion

In the following sections we show the reconstruction results that demonstrate the feasibility of the multi-spectral PAET reconstruction method. The multi-spectral PAET approach is examined using several simulations and in vitro experiments based on the reconstruction model and imaging systems mentioned above.

Figures 2 and 3 showed the recovered results for simulation case 1, in which the figures on the left column of Fig. 2 displayed the exact locations of the targets including HbO₂, HbR, H₂O and K whereas the recovered images were presented on the right Figs. 2(b) and 3 corresponding to HbO₂, HbR, H₂O and K, respectively. It was discovered from the recovered results in Figs. 2 and 3 that the crosstalk errors between the chromophore concentrations and bulk elastic modulus were effectively resolved by using our multi-spectral PAET reconstruction method. Meanwhile, the sizes and the quantitative physiological and elastic properties of the targets were quantified by estimating the full width half maximum of the recovered quantitative properties, which were provided in Table 3. Interestingly, for the reconstructed sizes, we discovered that the errors ranged from 5% to 7% whereas for the recovered quantitative elastic or physiology properties, the errors were less than 12.5%.

Fig. 3 The exact (left side) and recovered (right side) K images using 6 optical wavelengths for the simulation test 1 based on the test geometry in Fig. 2(a). The axes (left and bottom) illustrate the spatial scale, in mm, whereas the color scale (right) denotes K in GPa.

Download Full Size | PDF

Table 3. The reconstructed chromophore concentrations and K for the targets in simulation test 1

View Table | View all tables in this article

For the second simulation test, the aim was to examine if the targets that had different HbR and elastic modulus contrast levels can be effectively recovered in terms of the position, size and quantitative physiology and elastic properties. The reconstruction results in Fig. 4 indicated that the multi-spectral PAET algorithm was very robust for simultaneous recoveries of tissue physiological and elastic properties. In particular, multi-spectral PAET can effectively minimize the crosstalk errors between HbR and K and improve the image quality for multi-parameter reconstructions. In addition, by estimating the full width half maximum of the recovered quantitative properties in Fig. 4, we discovered that the reconstructed errors for the target sizes as well as the quantitative HbR concentration and K were less than 10%.

Fig. 4 Quantitative HbR (left) and K (right) images recovered by our multi-spectral PAET. The axes (left and bottom) illustrate the spatial scale, in mm, whereas the color scale (right) denote HbR in μM and K in GPa.

Download Full Size | PDF

The first in vitro test was performed to explore the merits of the developed multi-spectral PAET reconstruction method. For this experiment, we used the chicken breast as the target, as plotted in Fig. 5(a). In addition, the optical properties and molar extinction coefficients of the chicken breast were provided in Table 4. In particular, we used balance scales to measure the mass of a piece of chicken breast, and then placed the chicken breast inside the water to measure its volume. The measured mass density of the chicken breast was 1045kg/m³. The imaging results were provided in Fig. 5, in which Figs. 5(b) and 5(c) showed the recovered HbR concentration and bulk elastic modulus K, respectively. Likewise, by estimating the full width half maximum of the quantitative properties profiles in Fig. 5, we found the reconstructed HbR concentration of the chicken breast was nearly 10 μM while the corresponding peak value of K was 2.5GPa. In particular, the reconstructed size of the chicken breast in length was 3.1mm, which was in good agreement with its true size of 3.0mm.

Fig. 5 The Photograph of the chicken breast used for the first in vitro experimental test (a), the recovered HbR concentration image (b), the recovered K image (c), the HbR profile plotted along the center of the target (d), and the K profile plotted along the center of the target (e). The axes (left and bottom) illustrate the spatial scale, in mm, whereas the color scale (right) denote HbR in μM and K in GPa.

Download Full Size | PDF

Table 4. Optical properties and extinction coefficients of the chicken breast

View Table | View all tables in this article

For the second in vitro experimental test, the liver samples were imaged by using our PAT imaging system and the recovered K and HbR concentration were presented in Fig. 6. By estimating the full width half maximum of the quantitative properties profiles in Fig. 6, the reconstructed HbR concentration of subovate porcine liver was around 241.7 μM while the recovered bulk elasticity modulus K was about 2.6 GPa.

Fig. 6 The Photograph of the porcine liver tissue used for the second in vitro experimental test (a), the recovered HbR concentrations image (b), the recovered K (c), the HbR profile plotted along the center of the target (d), and the K profile plotted along the center of the target (e). The axes (left and bottom) illustrate the spatial scale, in mm, whereas the color scale (right) denote HbR in μM and K in GPa.

Download Full Size | PDF

In addition, it should be pointed out that to access the recovered quantitative properties of the targets by using the full width half maximum method, six profiles plotted along the center of the target were extracted and used for further analysis. For each profile, the quantitative property at the center of the target was considered as the reconstructed result. Then we calculated the mean value of the reconstructed results from the six profiles, which was considered as the recovered value for the target. It should also be pointed out that for the in vivo case, it may be very challenging to recover the quantitative elastic and physiological properties due to the highly heterogeneous distributions of optical properties from different biological tissues. In particular, for the in vivo test, the absorption spectra are very hard to be accurately determined. However, it is possible to recover the quantitative parameters if the reconstruction regions are focused on the tumor or lesion areas.

Further, to experimentally examine the reconstruction accuracy of the developed reconstruction method, it is essential to compare the recovered elastic moduli with the actual ones from different tissues (Table 5). It is noted that the actual bulk elastic modulus K was determined as $K = ρ v_{s}^{2}$ . Because the density $ρ$ was measured in this study for different experimental tests, the actual bulk elastic modulus can easily be calculated by using the following longitudinal velocities for the different tissues [21]. We found the reconstruction error between the recovered bulk elastic modulus and actual one was less than 15%.

Table 5. Average longitudinal velocities for the different phantoms [21]

View Table | View all tables in this article

In summary, our reconstruction results indicated that the quantification of chromophore concentrations and elastic bulk modulus is feasible by using the developed PAET imaging method when spectral priors are incorporated into the inverse reconstruction. In particular, images that are generated from simulated and in vitro data are presented, in which both the physiological and K properties images can be reconstructed based on our multi-spectral PAET reconstruction algorithm. Consequently, the spectrally resolved PAET approach provides an efficient means of concurrent reconstruction of multiple parameters including different chromophore concentrations and K. More importantly, it is clear from the reconstructed results that the crosstalk issue for recovery of multiple parameters are successfully resolved by using multi-spectral PA measurements. Specifically, multi-spectral PAET can effectively minimize the crosstalk error and improve the image quality for all the test cases.

According to the previous work, we know that the ranges of the elastic moduli of soft tissues span over as much as four orders of magnitudes [21]. Consequently, the recovered elastic property by using PAET can be applied to the diagnosis and therapeutic treatment of diseases that are characterized by changes in the elastic property of biological tissues. In addition, we discovered that the PAET image reconstruction methods are also able to quantitatively reconstruct absorbing objects of different sizes and different elastic modulus contrast levels. Nonetheless, the ability of reconstructing the elastic property image by using PA measurements may provide us a novel tool to quantify physiological function, disease progression, or response to intervention.

Acknowledgments

This study was supported by grants MYRG2014-00093-FHS, MYRG 2015-00036-FHS and MYRG2016-00110-FHS from the University of Macau in Macau and grants FDCT 026/2014/A1 and FDCT 025/2015/A1 from the Macao government.

References and links

1. L. V. Wang and H. Wu, Biomedical Optics: Principles and Imaging (Wiley, 2007).

2. G. Paltauf and H. Schmidt-Kloiber, “Pulsed optoacoustic characterization of layered media,” J. Appl. Phys. 88(3), 1624–1631 (2000). [CrossRef]

3. A. A. Oraevsky, S. L. Jacques, and F. K. Tittel, “Measurement of tissue optical properties by time-resolved detection of laser-induced transient stress,” Appl. Opt. 36(1), 402–415 (1997). [CrossRef] [PubMed]

4. R. J. Zemp, “Quantitative photoacoustic tomography with multiple optical sources,” Appl. Opt. 49(18), 3566–3572 (2010). [CrossRef] [PubMed]

5. J. Ripoll and V. Ntziachristos, “Quantitative point source photoacoustic inversion formulas for scattering and absorbing media,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 71(3), 031912 (2005). [CrossRef] [PubMed]

6. B. Cox, S. Arridge, K. Kostli, and P. Beard, “2D quantitative photoacoustic image reconstruction of absorption distributions in scattering medium using a simple iterative method,” Appl. Opt. 45, 1866–1875 (2006). [CrossRef] [PubMed]

7. C. Xu, P. Kumavor, A. Aguirre, and Q. Zhu, “Quantitative recovery of absorption coefficient using DOT-assisted photoacoustic tomography for breast imaging,” BsuD in OSA Biomedical Optics Topical Meeting, Miami, FL, 2010. [CrossRef]

8. S. Yang, D. Xing, Y. Lao, D. Yang, L. Zeng, L. Xiang, and W. R. Chen, “Noninvasive monitoring of traumatic brain injury and post-traumatic rehabilitation with laser-induced photoacoustic imaging,” Appl. Phys. Lett. 90(24), 243902 (2007). [CrossRef]

9. J. T. Oh, M. L. Li, H. F. Zhang, K. Maslov, G. Stoica, and L. V. Wang, “Three-dimensional imaging of skin melanoma in vivo by dual-wavelength photoacoustic microscopy,” J. Biomed. Opt. 11(3), 034032 (2006). [CrossRef] [PubMed]

10. A. Agarwal, S. Huang, M. O’donnell, K. C. Day, M. Day, N. Kotov, and S. Ashkenazi, “Targeted gold nanorod contrast agent for prostate cancer detection by photoacoustic imaging,” J. Appl. Phys. 102(6), 064701 (2007). [CrossRef]

11. X. Wang, W. W. Roberts, P. L. Carson, D. P. Wood, and J. B. Fowlkes, “Photoacoustic tomography: a potential new tool for prostate cancer,” Biomed. Opt. Express 1(4), 1117–1126 (2010). [CrossRef] [PubMed]

12. Z. Yuan, Q. Zhang, and H. Jiang, “Simultaneous reconstruction optical and acoustic properties by quantitative photoacoustic tomography,” Opt. Express 14, 6749–6754 (2006). [CrossRef] [PubMed]

13. J. Shah, S. Park, S. Aglyamov, T. Larson, L. Ma, K. Sokolov, K. Johnston, T. Milner, and S. Y. Emelianov, “Photoacoustic imaging and temperature measurement for photothermal cancer therapy,” J. Biomed. Opt. 13(3), 034024 (2008). [CrossRef] [PubMed]

14. M. S. Singh and H. Jiang, “Elastic property attributes to photoacoustic signals: an experimental phantom study,” Opt. Lett. 39(13), 3970–3973 (2014). [CrossRef] [PubMed]

15. Z. Yuan, X. Li, and L. Xi, “Listening to light scattering in turbid media: quantitative optical scattering imaging using photoacoustic measurements with one-wavelength illumination,” J. Opt. 16, 065301 (2014).

16. Z. Yuan and H Jiang, “Quantitative photoacoustic tomography: recovery of optical absorption coefficient maps of heterogeneous medium,” Appl. Phys. Lett. 88, 231101 (2006).

17. Y. Liu, H. Jiang, and Z. Yuan, “Two-scheme for quantitative photoacoustic tomography based on Monte Carlo simulation,” Med. Phys. 43(7), 1–11 (2016). [CrossRef] [PubMed]

18. Z. Yuan and H. Jiang, “Simultaneous recovery of tissue physiological and acoustic properties and the criteria for wavelength selection in multispectral photoacoustic tomography,” Opt. Lett. 34(11), 1714–1716 (2009). [CrossRef] [PubMed]

19. S. Li, B. Montcel, Z. Yuan, W. Liu, and D. Vray, “Multigrid-based reconstruction algorithm for quantitative photoacoustic tomography,” Biomed. Opt. Express 6(7), 2424–2434 (2015). [CrossRef] [PubMed]

20. F. A. Duck, Physical properties of tissues: a comprehensive reference book (Academic press, 2013).

21. T. Glozman and H. Azhari, “A method for characterization of tissue elastic properties combining ultrasonic computed tomography with elastography,” J. Ultrasound Med. 29(3), 387–398 (2010). [PubMed]

22. H. Jiang, Z. Yuan, and X. Gu, “Spatially varying optical and acoustic property reconstruction using finite-element-based photoacoustic tomography,” J. Opt. Soc. Am. A 23(4), 878–888 (2006). [CrossRef] [PubMed]

23. http//omlc.ogi.edu/spectra/index.html.

24. J. P. Ritz, A. Roggan, C. Isbert, G. Müller, H. J. Buhr, and C. T. Germer, “Optical properties of native and coagulated porcine liver tissue between 400 and 2400 nm,” Lasers Surg. Med. 29(3), 205–212 (2001). [CrossRef] [PubMed]

Objects position	HbO₂ (μM)	HbR (μM)	H₂O (%)	K (GPa)
Background	6	2	20	2.2
Top right	6	10	20	2.2
Top left	6	2	70	2.8
Bottom middle	20	2	20	2.8

Objects position	HbR (μM)	K(GPa)
Background	2	2.2
Top right	20	2.2
Top left	2	4.0
Bottom middle	10	4.0

Targets	HbO₂ (μM)	HbR (μM)	H₂O (%)	K (GPa)	Sizes errors (%)
Top right	5	11	22	2.4	5
Top left	5	2	80	2.9	7
Bottom middle	21	2	22	2.9	5

Wavelengths (nm)	Target positions	Molar extinction coefficient (cm⁻¹/M)	Absorptioncoefficient (cm⁻¹)
532	Background Chicken breast	40580	0.01 0.8
680	Background Chicken breast	2407	0.01 0.046

Objects position	HbO₂ (μM)	HbR (μM)	H₂O (%)	K (GPa)
Background	6	2	20	2.2
Top right	6	10	20	2.2
Top left	6	2	70	2.8
Bottom middle	20	2	20	2.8

Multi-spectral photoacoustic elasticity tomography

Abstract

1. Introduction

2. Methods and materials

2.1 The multi-spectral PAET reconstruction method

2.2 The multi-spectral PAET imaging systems

2.3 Simulation and in vitro experimental tests

3. Results and discussion

Acknowledgments

References and links

Cited By

Figures (6)

Tables (5)

Equations (21)

Biomedical Optics Express