1.

Introduction

Over the past few years, photoacoustic microscopy (PAM) has been proven to be capable of structural, functional, molecular, and metabolic imaging.^1–4 In PAM, the object is illuminated by a short-pulsed laser beam. Following the absorption of light, the increased temperature generates an initial pressure rise, which propagates as photoacoustic (PA) waves and is detected by an ultrasonic transducer.⁵ Because the initial pressure is directly proportional to the absorbed optical energy density $A(\overrightarrow{r})$ ($\mathrm{J}·{\mathrm{m}}^{-3}$) in the tissue, which is the product of the optical fluence $F(\overrightarrow{r})$ ($\mathrm{J}·{\mathrm{m}}^{-2}$) and the optical absorption coefficient ${\mu }_{\mathrm{a}}(\overrightarrow{r})$ (${\mathrm{m}}^{-1}$), multiwavelength PA measurements can provide spectral information of optical absorption. Moreover, because it is essentially background free for blood-vessel imaging, PAM has better imaging contrast than other back-scattering-based optical imaging modalities. So far, PAM has shown its viability in detecting many intrinsic contrasts in biological tissues, such as hemoglobin,⁶ bilirubin,⁷ DNA and RNA in nuclei,⁸ lipid,⁹ cytochromes,¹⁰ and melanin.¹¹ To obtain the intrinsic contrast property ${\mu }_{\mathrm{a}}(\overrightarrow{r})$, we need to compensate for the extrinsic quantity $F(\overrightarrow{r})$. However, the optical fluence is usually unknown because of the light attenuation in tissue.¹² Therefore, it is generally challenging to provide quantitative PA studies of, for example, oxygen metabolism.

Currently, the optical fluence $F(\overrightarrow{r})$ can be compensated for either invasively¹³ or noninvasively.^14–17 In the invasive method, an optical absorber with a priori known absorption spectrum is placed close to the target of interest and used as a reference. PA amplitude from the target of interest is normalized by the measured PA amplitude of this reference absorber. A noninvasive method involves solving both optical diffusion and PA wave equations. However, the inverse problems are computationally intensive and ill-posed due to optical diffusion. In another noninvasive method, the temporal profiles of PA signals are used to determine the absolute absorption coefficient.^18–21 However, the acoustic attenuation and the limited bandwidth of the transducer may distort the PA signals. Thus, the direct fitting method may be inaccurate. To avoid this issue, our group has proposed to use acoustic spectra to determine the absolute absorption coefficient.^22,23

In this paper, we present another quantitative approach to measure absolute absorber concentrations by statistical analyses of PA signals. To obtain a fluence-independent detection, the particle count in the detection volume needs to be small enough so that PA signal fluctuation due to the particle Brownian motion is dominant compared to other fluctuating sources, such as laser intensity fluctuation, electronic thermal noise, and photon shot noise. This method is analogous to fluorescence correlation spectroscopy and dynamic light scattering (also known as photon correlation spectroscopy). The former has been widely used to measure the fluorescence intensity fluctuations, while the latter has been used to measure the scattered light intensity fluctuations.

The paper is organized as follows. We first developed a model to describe the relationship between the measured PA signal and the particle count. The principle was then demonstrated experimentally. After samples of varied microsphere concentrations were studied, red blood cells (RBCs) were measured to show that our method can be applied to biological samples. Finally, we performed in vivo experiments to demonstrate the potential of our method in biological applications.

2.

Methods and Materials

Particle fluctuation in a detection volume caused by Brownian motion is governed by Poisson distribution.²⁴ Thus, the variance of the particle count ($N_{\mathrm{p}}$) equals its mean, i.e.,

Because the thermal noise is zero-mean, the mean and variance of $A_{\mathrm{PD}}$ are derived from Eq. (2) as

Based on Eqs. (4) and (5), we have

Equation (6) shows a parabolic relationship between $\mathrm{Var}(A_{\mathrm{PD}})$ and $\mathrm{E}(A_{\mathrm{PD}})$.

Similar to the analysis of PD, the PA amplitude ($A_{\mathrm{PA}}$) can be expressed as

Based on Eq. (7), the mean and the variance of $A_{\mathrm{PA}}$ are

Substituting Eq. (1) into Eq. (9), we have

The ratio between the mean squared and the variance of $A_{\mathrm{PA}}$ is the power-based signal-to-noise ratio (SNR):

Defining ${\mathrm{SNR}}_{\mathrm{th}}=\frac{{\mathrm{E}}^2(A_{\mathrm{PA}})}{\mathrm{Var}(n_{\mathrm{PA}})}$, and substituting it into Eq. (12), we have

Equation (14) can be discussed in two different cases. In case 1, when particles are dense enough so that $\mathrm{E}(N_{\mathrm{p}})\gg {\alpha }^{-2}$. we have

${\mathrm{SNR}}_{\mathrm{PA}}$ approaches ${\alpha }^{-2}$ when ${\mathrm{SNR}}_{\mathrm{th}}$ increases. In case 2, when particles are sparse enough so that $\mathrm{E}(N_{\mathrm{p}})\gg {\alpha }^{-2}$, we have

${\mathrm{SNR}}_{\mathrm{PA}}$ approaches $\mathrm{E}(N_{\mathrm{p}})$ when ${\mathrm{SNR}}_{\mathrm{th}}$ increases. Using Eq. (16) or more generally Eq. (14), one can measure the expected value of the absolute particle count. In addition, if $\mathrm{E}(N_{\mathrm{p}})\ll {\mathrm{SNR}}_{\mathrm{th}}$, Eq. (16) reduces to

Thus, the measurement outcome is independent of the laser fluence. Although our model is derived for static particle suspension, it is also valid for flowing particles because the temporal statistics of particle count fluctuations in a static medium are equivalent to the spatial statistics of the particle count fluctuations in a flowing medium.²⁴

We validated this idea using our optical-resolution PAM (OR-PAM) system shown in Fig. 1.²⁵ A tunable dye laser (CBR-D, Sirah, GmbH, Grevenbroich and Germany), pumped by a Nd:YLF laser (INNOSAB, Edgewave, GmbH, Würselen and Germany) at a 1 kHz pulse repetition rate, provides nanosecond laser pulses. The laser beam is filtered by a 50-μm diameter pinhole and focused into the sample by a microscope objective lens (numerical aperture is 0.1). An acoustic-optical beam combiner, composed by two prisms sandwiching a layer of silicone oil, provides the acoustic-optical coaxial alignment. Acoustic detection is achieved with an ultrasonic transducer (V214-BC, Panametrics-NDT, Inc., Waltham) and a spherically focused acoustic lens, which is placed confocally with the optical objective lens. The detected PA signal is then amplified, digitized, and stored. Unlike confocal and two-photon microscopy, depth scanning is not needed in OR-PAM for volumetric imaging. A three-dimensional image is obtained by raster scanning the sample and piecing together the depth-resolved one-dimensional images (A-lines). So far, this OR-PAM system has achieved a lateral resolution of $\sim 5\mu \text{m}$ and an axial resolution of $\sim 15\mu \text{m}$ with an imaging depth of $\sim 1\mathrm{mm}$.²⁵ Thus, the detection volume is about $500{\mu \text{m}}^3$ based on $1/e$ excitation beam width. In all experiments, plastic tubing with inner diameter of 300 μm was used to flow the liquid sample. A fast flow speed of $9\mathrm{mm}/\mathrm{s}$ was applied to the fluid in the tubing to avoid potential heat aggregation or photobleaching, and to make sure the multiple PA signals were completely independent.

Fig. 1

Optical-resolution photoacoustic microscopy schematic. BS, beam splitter; ConL, condenser lens; CorL, correction lens; EP, eyepiece; ND, neutral density filter; OL, objective lens; PD, photodiode; PH, pinhole; RAP, right-angle prism; RhP, rhomboid prism; SOL, silicone oil layer; UT, ultrasonic transducer; and WT, water tank.

To calculate $\mathrm{Var}(A_{\mathrm{PA}})$ and $\mathrm{E}(A_{\mathrm{PA}})$, for each particle count, we measured 1500 A-lines with the same incident laser energy. Representative $A$-lines and corresponding PA amplitudes are shown in Figs. 2(a) and 2(b), respectively. Thus, $\mathrm{Var}(A_{\mathrm{PA}})$ and $\mathrm{E}(A_{\mathrm{PA}})$ can be obtained by calculating the variance and mean of the PA amplitudes. The background of each A-line inside the dashed box in Fig. 2(a) was used to calculate the electronic thermal noise $n_{\mathrm{PA}}$. Then we tuned the neutral density filter, shown in Fig. 1, to change the incident laser energy. Thus, different ${\mathrm{SNR}}_{\mathrm{PA}}$ and ${\mathrm{SNR}}_{\mathrm{th}}$ could be obtained. Meanwhile, for each PA A-line, we monitored the incident energy by using a PD. Similarly, we could also calculate $\mathrm{Var}(A_{\mathrm{PD}})$ and $\mathrm{E}(A_{\mathrm{PD}})$ after 1500 measurements. By selecting PA A-lines with corresponding $\mathrm{E}(A_{\mathrm{PD}})$ only in a predetermined range, the photon count fluctuation $\alpha$ can be controlled. If $\mathrm{E}(N_{\mathrm{p}})\ll {\alpha }^{-2}$, the absolute particle count can be quantified based on Eqs. (16) or (17).

Fig. 2

Representative PA signals (a) and PA amplitudes (b) of multiple A-lines with the same incident laser energy for each particle count.

Lysed bovine blood was used for the experiment with large particle counts. It is estimated that there are about one billion hemoglobin molecules in the detection volume,²⁶ which satisfies the condition that $\mathrm{E}(N_{\mathrm{p}})$ is much larger than ${\alpha }^{-2}$ for Eq. (15). For small particle counts, defibrinated whole bovine blood was used instead for the experiment. In whole blood, RBCs, instead of hemoglobin molecules, constitute the particles moving in and out of the detection volume. Since there are only a few RBCs in the detection volume, $\mathrm{E}(N_{\mathrm{p}})$ will be much smaller than ${\alpha }^{-2}$, satisfying the condition for Eq. (16).

To demonstrate quantitative particle count measurement, we prepared samples with different particle concentrations. Since we require that particle count fluctuation be dominant in the PA signal fluctuation, two different sizes of red dyed microspheres were used to make the samples: a 0.5-μm diameter particle with a stock concentration of $\sim 3.64\times {10}^{11}$ particles per milliliter (2.5% w/w), and a 1-μm diameter particle with a stock concentration of $\sim 4.55\times {10}^{10}$ particles per milliliter (2.5% w/w). The smaller particle stock solution was diluted by two and four times, while the larger particle stock solution was diluted to eight and 16 times the stock concentration of the smaller particle solution.

We used defibrinated bovine blood in two experiments: an ex vivo experiment and an in vivo mimicking experiment. For the ex vivo experiment, defibrinated bovine blood with four different concentrations was prepared. The original undiluted defibrinated bovine blood ($150\mathrm{g}/\mathrm{L}$) was serially diluted to 3/4, 1/2, and 1/4 of the original concentration. Since there were too few RBCs in the original resolution voxel of OR-PAM, we increased the detection volume to $1.2\times {10}^{-5}{\mathrm{mm}}^3$ by defocusing the light beam. For the simulated in vivo experiment, the original undiluted defibrinated bovine blood ($150\mathrm{g}/\mathrm{L}$) was serially diluted to 1/2 and 1/4 of the original concentration.

3.

Results

3.1.

Model Demonstration

Based on Eqs. (6) and (11), $\mathrm{Var}(A_{\mathrm{PD}})$ and $\mathrm{Var}(A_{\mathrm{PA}})$ parabolically depend on $\mathrm{E}(A_{\mathrm{PD}})$ and $\mathrm{E}(A_{\mathrm{PA}})$, respectively. In addition, $\mathrm{E}(A_{\mathrm{PA}})$ should be proportional to $\mathrm{E}(A_{\mathrm{PD}})$. In the experiment using lysed bovine blood, the result shows a clear parabolic relationship between $\mathrm{Var}(A_{\mathrm{PD}})$ and $\mathrm{E}(A_{\mathrm{PD}})$ [Fig. 3(a)]. The photon count fluctuation coefficient $\alpha$ was fitted to be 5.67%, close to the laser specification. The similar result between $\mathrm{Var}(A_{\mathrm{PA}})$ and $\mathrm{E}(A_{\mathrm{PA}})$ is shown in Fig. 3(b). The coefficient $\alpha$ was fitted to be 5.76%, very close to the one for PD. In addition, Fig. 3(c) shows that there is a good linearity between $\mathrm{E}(A_{\mathrm{PA}})$ and $\mathrm{E}(A_{\mathrm{PD}})$.

Fig. 3

Relationships between (a) $\mathrm{Var}(A_{\mathrm{PD}})$ and $\mathrm{E}(A_{\mathrm{PD}})$, (b) $\mathrm{Var}(A_{\mathrm{PA}})$ and $\mathrm{E}(A_{\mathrm{PA}})$, (c) $\mathrm{E}(A_{\mathrm{PA}})$ and $\mathrm{E}(A_{\mathrm{PD}})$ as estimated from the experimental data. Fitting equation for (a): $\mathrm{Var}(A_{\mathrm{PD}})={\alpha }^2\times {\mathrm{E}}^2(A_{\mathrm{PD}})+2631$, where $\alpha =0.0567$, coefficient of determination $R^2=0.9994$. Fitting equation for (b): $\mathrm{Var}(A_{\mathrm{PA}})={\alpha }^2\times {\mathrm{E}}^2(A_{\mathrm{PA}})+7508$, where $\alpha =0.0576$ and $R^2=0.9974$. Fitting equation for (c): $\mathrm{E}(A_{\mathrm{PA}})=0.6656\times \mathrm{E}(A_{\mathrm{PD}})-190.7$, $R^2=0.9964$. (a) and (b) show that when the particle count is sufficiently large, both $\mathrm{Var}(A_{\mathrm{PA}})$ and $\mathrm{Var}(A_{\mathrm{PD}})$ follow a parabolic relationship with $\mathrm{E}(A_{\mathrm{PA}})$ and $\mathrm{E}(A_{\mathrm{PD}})$, respectively. The coefficients $\alpha$ for PD and PA are almost identical both due to the photon count fluctuation. (c) shows a good linearity between $\mathrm{E}(A_{\mathrm{PA}})$ and $\mathrm{E}(A_{\mathrm{PD}})$.

As mentioned above, we can regroup the laser pulses based on the PD readings. By doing this, the laser intensity fluctuation coefficient $\alpha$ can be artificially controlled to a predetermined value. Figure 4(a) shows that with large particle counts, ${\mathrm{SNR}}_{\mathrm{PA}}$ is dominated by the light fluctuations when ${\mathrm{SNR}}_{\mathrm{th}}$ increases, consistent with Eq. (15). Here, the neutral density filter was adjusted to change the laser fluence and thus to change the ${\mathrm{SNR}}_{\mathrm{th}}$. On the other hand, when the particle counts are small in the detection volume, the measurement results with high ${\mathrm{SNR}}_{\mathrm{th}}$ only rely on the particle counts [Figs. 4(b)–4(d)], consistent with Eq. (16). Although the photon count fluctuation coefficients are varied, the measurement results stay the same. Thus, the fluence-independent measurement is achieved with the conditions of high ${\mathrm{SNR}}_{\mathrm{th}}$ and small particle counts.

Fig. 4

Measurements for different photon count fluctuation coefficients. (a) For larger particle counts, the ${\mathrm{SNR}}_{\mathrm{PA}}$ at sufficiently high ${\mathrm{SNR}}_{\mathrm{th}}$ depends only on photon count fluctuation. (b)–(d) If the particle count is small, the ${\mathrm{SNR}}_{\mathrm{PA}}$ at sufficiently high ${\mathrm{SNR}}_{\mathrm{th}}$ depends primarily on the particle count. (b) $\alpha =0.0066$. (c) $\alpha =0.0284$. (d) $\alpha =0.0587$. The dots are the measurement data and the solid curves are the fitted results based on Eq. (14).

3.2.

Measurement of Particle Counts

Next, we applied our method to quantitatively measure five different particle concentration samples. Figure 5(a) shows measured results with five different preset particle counts of 440, 220, 110, 55, and 27 in the detection volume. The experimental results agree well with the preset values. In addition, we observed that the measured particle count does not change with the laser fluence any longer when the ${\mathrm{SNR}}_{\mathrm{th}}$ is high enough. The fitted coefficient $\mathrm{E}(N_{\mathrm{p}})$ in the model was used to quantify the absolute particle count for each sample. As shown in Fig. 5(b), the experimental results agree well with the preset particle counts in the detection volume.

Fig. 5

Quantitatively measuring particle concentration from different samples with the same light fluctuation. (a) The preset particle counts are 440 (purple), 220 (blue), 110 (red), 55 (green), and 27 (black). Solid lines: theoretical fitting based on Eq. (16). (b) Measured particle count versus the preset value.

3.3.

Measurement of RBC Counts Ex Vivo

We also measured RBCs to validate the feasibility of our method in biological applications. Figure 6(a) illustrates the measured results with the fitting curve for four different preset RBC counts of 400, 300, 200, and 100 in the detection volume. We used the fitted coefficient $\mathrm{E}(N_{\mathrm{p}})$ to represent the counts of RBCs and compared with the corresponding preset RBC counts [Fig. 6(b)]. The results demonstrate that our method is fluence-independent and can be used to quantitatively measure biological samples.

Fig. 6

Measured RBC counts under the same light fluctuation. (a) The preset particle counts are 400 (blue), 300 (red), 200 (green), and 100 (black). Solid lines: theoretical fitting based on Eq. (16). (b) Measured particle count versus the preset value.

3.4.

Measurement of RBC Counts in the Simulated In Vivo Case

For in vivo application, additional signal fluctuation mainly comes from the perturbation of laser fluence by micro-flows in the light path. To mimic the in vivo situation while maintaining the flexibility of varying the RBC concentration, we covered targets made of black tape and bovine-blood-filled tubing with a living mouse ear, as shown in Figs. 7(a)–7(c).

Fig. 7

Measurement of RBC counts in tubing covered by a living mouse ear. (a) Photo of the experimental setup. (b) Close-up of the region within the red dashed box in (a), showing the blood-filled tubing and black-tape targets covered by the living mouse ear. (c) A typical cross-sectional PA image showing that the tubing and part of the black tape were covered by the mouse ear. (d) Measurement of PA amplitude fluctuation from black tape at 40 random spots covered by the mouse ear.

We first measured the signal fluctuations from black tape, which did not show Brownian motion. As shown in Fig. 7(d), measurements were made at 40 random spots on the black tape. For all the spots, the PA signal fluctuations were around 2%. As described before, the power-based signal-to-noise ratio (${\mathrm{SNR}}_{\mathrm{PA}}$) was used to denote the measured particle counts. Since a 2% fluctuation corresponds to an ${\mathrm{SNR}}_{\mathrm{PA}}$ of $4\times {10}^{-4}$, it is negligible. The experimental results show that the in vivo PA signal fluctuation due to other sources rather than Brownian motion is insignificant.

Furthermore, we measured RBC counts in the same setup. Figure 8(a) illustrates the measured results, with fitting curves for three different preset RBC counts of 1080, 540, and 270 in the detection volume. We used the fitted coefficient $\mathrm{E}(N_{\mathrm{p}})$ to represent the RBC counts. As shown in Fig. 8(b), the experimental results agree with the preset RBC counts. The result shows that our method can be used to measure RBC counts in an in vivo mimicking situation using OR-PAM.

Fig. 8

RBC count measurement with the blood-filled tubing covered by a mouse ear. (a) The preset RBC counts in the voxel were 1080 (red), 540 (green), and 270 (black). (b) Measured versus preset RBC counts.

3.5.

Measurement of RBC Counts In Vivo

To show the potential for actual in vivo applications of our method, we performed two more demonstrations, which correspond to anemia caused by blood loss and by RBC aggregation, respectively. In each in vivo experiment on RBC count, we maintained the laser fluence constant (see Sec. 4). In the first demonstration, saline was injected into mouse tail vein to mimic blood loss (blood dilution), as shown in Fig. 9(a). The same vein was monitored to quantify the RBC count changes. Figure 9(b) shows a PA image of the vein, where the black cross denotes the monitoring spot. We slowly injected $\sim 0.03\mathrm{mL}$ saline into the tail vein during a 3 min window. As shown in Figs. 9(c) and 9(d), the measured RBC count was reduced by injected saline. As soon as we stopped injection, the RBC count quickly returned to its normal value.

Fig. 9

In vivo RBC count changes in a mouse tail vein when injecting saline. (a) Animal setup in the experiment. The imaged area was downstream of the saline injection spot. (b) PA MAP image of part of the tail in vivo. The large vessel in the middle is the mouse tail vein where we injected the saline. The black cross denotes the monitoring spot. (c) Time course of PA amplitude. The red dashed box highlights a clear decrease in the PA signal mean value. (d) RBC counts measured by our method before, during, and after saline were injected into the tail vein.

In the other demonstration, to induce RBC aggregation, we used $\sim 0.3\mathrm{mL}$ 3% dextran-70 solution to mix with 3 mL original RBC solution. The diluted original and aggregated RBC solutions were imaged ex vivo by a commercial microscope, as shown in Fig. 10. The results show that we can induce RBC aggregation successfully ex vivo and probably in vivo.

Fig. 10

RBC aggregation induced by 3% dextran-70 solution. The concentration of RBCs was the same in the two images. (a) Original and (b) aggregated RBCs imaged by a commercial microscope.

Then we performed an in vivo experiment to observe RBC aggregation. The same dextran-70 solution was injected into a mouse through the tail vein, and the mouse ear was used to monitor the RBC aggregation, as shown in Fig. 11(a). Figure 11(b) shows the PA image of the mouse ear, where the red cross denotes the RBC monitoring spot. As shown in Figs. 11(c) and 11(d), an increase of PA signal fluctuation and corresponding RBC count decrease were observed about 50 min after the dextran-70 solution injection. Based on the results of our ex vivo experiment, we think that this decrease in RBC count was from RBC aggregation, induced by the mixing of dextran-70 solution with blood.

Fig. 11

RBC aggregation measured in vivo. (a) Animal mounting in the experiment. Dextran-70 solution was injected into the mouse through the tail vein. The mouse ear was imaged and monitored over the time. (b) PA MAP image of the mouse ear. The red cross denotes where we monitored the RBC count. (c) Time course of PA amplitude. The red dashed box highlights a clear increase in the PA signal fluctuation. (d) Time course of RBC count measured by our method. The selected red dashed box window highlights a decrease in RBC count.

4.

Discussion

Both the theoretical model and experimental study showed that our method is applicable when the particle count in the detection volume is sufficiently small such that the particle Brownian motion dominates the PA signal fluctuation. If we can control the photon count fluctuation coefficient within 1%, the maximum count of particles that we can measure is on the order of 1000 ($\ll 1/{\alpha }^2$). By contrast, the minimum measurable count of particles is limited by ${\mathrm{SNR}}_{\mathrm{th}}$. If the count of particles in the detection volume is too small, noise reduces the measuring accuracy significantly. Either particles with larger absorption cross sections or greater detection volumes can reduce the minimum. Based on our experimental results [Figs. 4(b)–4(d)], we conclude that our method can measure as few as three RBCs in the detection volume.

As shown in Sec. 3 (Fig. 5), our measured particle counts are slightly larger than the preset absolute counts estimated based on the detection volume of $500{\mu \text{m}}^3$. This is most likely due to the inaccurate estimation of the detection volume. We calculated the detection volume based on the $1/e$ detection sensitivity field. However, particles not in this volume can still generate PA signals, which will contribute to the total PA signals. If we take those particles into account, the actual detection volume becomes greater.

In addition, as shown in Figs. 4(a), 5(a), and 6(a), when ${\mathrm{SNR}}_{\mathrm{th}}$ becomes high enough, the measured particle counts are independent of the ${\mathrm{SNR}}_{\mathrm{th}}$ and laser fluence. Although we examined the complete trend in this study with varying ${\mathrm{SNR}}_{\mathrm{th}}$, a measurement with a single, high enough ${\mathrm{SNR}}_{\mathrm{th}}$ is sufficient to measure the fluence-independent particle counts in the detection volume, as we have done in the in vivo experiment.

In our experiments, we chose a flowing medium for the following two reasons. First, we focused on measuring counts of RBCs, which usually flow constantly in vivo. Second, using a flowing medium can avoid potential heat accumulation and photobleaching, which, however, can be mitigated by decreasing the pulse repetition rate. For example, the heat diffusion time is about 0.2 ms in our experiment (characteristic dimension: 5 μm; thermal diffusivity: $1.3\times {10}^{-3}{\mathrm{cm}}^2/\mathrm{s}$).⁵ Thus, a laser repetition rate of 1 kHz causes minimal heat accumulation, which is further reduced by blood perfusion.

Finally, we may be able to apply our method in deep tissue imaging modalities such as acoustic resolution PAM and PA computed tomography where fluence compensation is even more challenging. However, compared with OR-PAM, their greater detection volume may require even lower particle concentrations. Although it may be difficult to detect absolute RBC counts, we can use our method to monitor the absolute counts of other cells in blood, such as intact or labeled white blood cells,²⁷ which are much sparser in the blood stream.

5.

Conclusion

To the best of our knowledge, this is the first study using PAM for fluence-independent particle count measurement by statistically analyzing PA signals. This method has been experimentally verified by using different kinds of samples, including RBCs. We also demonstrated the feasibility of our method in the in vivo experiment. Therefore, our method shows promise for quantitatively measuring biological samples in vivo.