# Patent

## Title
Quantification of optical absorption coefficients using acoustic spectra in photoacoustic tomography

## Abstract

Accurately quantifying optical absorption coefficient using acoustic spectra of photoacoustic signals. Optical absorption is closely associated with many physiological parameters, such as the concentration and oxygen saturation of hemoglobin, and it can be used to quantify the concentrations of non-fluorescent molecules. A sample is illuminated by, for example, a pulsed laser and following the absorption of optical energy, a photoacoustic pressure is generated via thermo-elastic expansion. The acoustic waves then propagate and are detected by a transducer. The optical absorption coefficient of the sample is quantified from spectra of the measured photoacoustic signals. Factors, such as system bandwidth and acoustic attenuation, may affect the quantification but are canceled by dividing the acoustic spectra measured at multiple optical wavelengths.

## Background
CROSS REFERENCE TO RELATED APPLICATIONS
This application is a continuation of U.S. application Ser. No. 13/637,897, filed Oct. 24, 2012 which was the National Stage of International Application No. PCT/US2011/031823, filed Apr. 8, 2011, which claims the benefit of U.S. Provisional Patent Application No. 61/322,605 filed Apr. 9, 2010, all of which are hereby incorporated by reference in their entirety.


STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH & DEVELOPMENT
This invention was made with government support under grant R01 EB008085, awarded by the U.S. National Institutes of Health. The government has certain rights in the invention.


BACKGROUND
The embodiments described herein relate generally to photoacoustic scanning methods and apparatus and, more particularly, to a method and a system for determining an optical absorption coefficient by an object using acoustic spectra.
Total, oxygenated, and deoxygenated hemoglobin concentrations ([HbT], [HbO2], and [HbR], respectively) are fundamental pathophysiological parameters in biomedicine. For example, abnormally low [HbT] may be caused by loss of blood, nutritional deficiency, chemotherapy, inflammation, kidney failure or bone marrow problems, while abnormally high [HbT] may be related to exposure to high altitude, smoking, dehydration and tumors. Blood oxygen saturation (sO2), which is defined as [HbO2] divided by [HbT], is vital in understanding brain hemodynamics in response to sensory stimulations, monitoring healing of burns and wounds, and evaluating the effectiveness of chemotherapy and radiotherapy on tumors. Several techniques have been developed to quantify hemoglobin concentration and sO2 in vivo, including near-infrared spectroscopy (NIRS), blood oxygen level dependent (BOLD) contrast magnetic resonance imaging (MRI), electron paramagnetic resonance imaging (EPRI), positron emission tomography (PET), and single photon emission computed tomography (SPECT). However, all of these modalities have disadvantages. For example, at least some of these modalities have poor spatial resolution, relative quantification, and undesirable contrast agent injection. Photoacoustic (PA) tomography (PAT) has already demonstrated its ability to monitor biological hemodynamic functions without using exogenous contrast agents. Quantitative PAT is challenging, because compensating for the fluence in quantitative in vivo photoacoustic tomography is difficult, and factors, such as the tissue acoustic attenuation and the imaging system bandwidth, also affect the quantification accuracy.
BRIEF DESCRIPTION
In one aspect, a photoacoustic imaging method includes illuminating an object using a light beam emitted by a light source, and detecting a pressure wave emitted by the object using an acoustic transducer, wherein the pressure wave is induced by the object in response to the light beam (e.g., light pulse). The object is illuminated with at least two light pulses having different wavelengths, and the absorption coefficient of the object is determined from the detected pressure waves. The method also includes calculating an amount of optical energy that is absorbed by the object from the light beam, based at least on the pressure.



BRIEF DESCRIPTION OF THE DRAWINGS
The embodiments described herein may be better understood by referring to the following description in conjunction with the accompanying drawings.
FIG. 1 is a schematic diagram of an exemplary optical resolution photoacoustic microscopy (OR-PAM) system.
FIG. 2 is a diagram of an exemplary experimental setup of a sample for OR-PAM.
FIG. 3 is a set of graphs that illustrate an example of an object spectrum O(ωλ), a system dependent response H(ω), and a tissue-related acoustic attenuation effect a(ω) acquired via the OR-PAM system of FIGS. 1 and 2.
FIG. 4 is a graph illustrating acoustic spectra of a number of photoacoustic signals.
FIG. 5 is a graph illustrating a comparison of photoacoustic spectra of two objects.
FIG. 6 is a graph illustrating quantified absorption coefficients of a plurality of samples.
FIG. 7 is a graph illustrating acoustic spectra of a photoacoustic signal with and without an optical phantom layer.
FIG. 8 is a photoacoustic maximum amplitude projection image of a blood vessel.
FIG. 9 is a graph that illustrates a ratio of acoustic spectra at two different optical wavelengths.
FIG. 10 is a diagram that illustrates raster scanning an acoustic transducer to obtain three-dimensional images via photoacoustic microscopy.
FIG. 11 is a diagram that illustrates light diffusion at or near a top surface of a blood vessel.
FIG. 12 is a diagram that illustrates a characterization of light on the top surface of the blood vessel of FIG. 11 as an isotropic point source.
FIG. 13 is a diagram that illustrates an alternative embodiment of an OR-PAM that includes a plurality of acoustic transducers arranged into an array.
FIG. 14 is a flow chart of a method of determining an absorption coefficient of a volume of an object.



DETAILED DESCRIPTION
While the making and using of various embodiments of the present disclosure are discussed in detail below, it should be appreciated that the present disclosure provides many applicable inventive concepts that can be embodied in a wide variety of specific contexts. The specific embodiments discussed herein are merely illustrative of specific ways to make and use the disclosure and do not delimit the scope of the disclosure.
To facilitate the understanding of the embodiments described herein, a number of terms are defined below. The terms defined herein have meanings as commonly understood by a person of ordinary skill in the areas relevant to the present disclosure. Terms such as “a,” “an,” and “the” are not intended to refer to only a singular entity, but rather include the general class of which a specific example may be used for illustration. The terminology herein is used to describe specific embodiments of the disclosure, but their usage does not delimit the disclosure, except as outlined in the claims.
To be consistent with the commonly used terminology, whenever possible, the terms used herein will follow the definitions recommended by the Optical Society of America (OCIS codes).
In some embodiments, the term “photoacoustic microscopy” refers generally to a photoacoustic imaging technology that detects pressure waves generated by light absorption in the volume of a material (such as biological tissue) and propagated to the surface of the material. In other words, photoacoustic microscopy is a method for obtaining three-dimensional images of the optical contrast of a material by detecting acoustic or pressure waves traveling from the object. Imaging resolution is on the micrometer scale.
In some embodiments, the term “photoacoustic tomography” also refers to a photoacoustic imaging technology that detects acoustic or pressure waves generated by light absorption in the volume of a material (such as biological tissue) and propagated to the surface of the material. The emphasis is sometimes on photoacoustic computed tomography, i.e., cross-sectional or three-dimensional photoacoustic imaging based on computer reconstruction, although the most general definition of photoacoustic tomography encompasses photoacoustic microscopy.
In some embodiments, the terms “reflection mode” and “transmission mode” refer generally to a laser photoacoustic microscopy system that employs the detection of acoustic or pressure waves transmitted from the volume of their generation to the optically irradiated surface and a surface that is opposite to, or substantially different from, the irradiated surface, respectively.
In some embodiments, the term “time-resolved detection” refers generally to the recording of the time history of a pressure wave with a temporal resolution sufficient to reconstruct the pressure wave profile.
In some embodiments, the term “photoacoustic waves” refers generally to pressure waves produced by light absorption.
In photoacoustic imaging, a sample is illuminated, usually by a pulsed laser, and following the absorption of optical energy, an initial pressure is generated via thermo-elastic expansion. The photoacoustic waves then propagate and are detected by acoustic or ultrasonic sensors. The strength of the initial pressure is directly proportional to the absorbed optical energy in the tissue and, therefore, multi-wavelength photoacoustic measurements may yield optical absorption spectral information. Since ultrasound scattering is much weaker than optical scattering in biological tissues, photoacoustic tomography is capable of high resolution imaging at depths beyond the optical transport mean free path. Moreover, because photoacoustic imaging inherently exploits optical absorption contrast, it usually has a higher sensitivity to optical absorption than other optical imaging technologies. By measuring the optical absorption spectrum with photoacoustic imaging, sO2 may be quantified in the same way near infrared spectroscopy (NIRS) quantifies sO2, except with high spatial resolution and sensitivity.
Photoacoustic images are spatial mappings of the absorbed optical energy A(r) (Jm−3), which is the product of the optical absorption coefficient μa(r, λ) (m −1) and the fluence F(r) (Jm−2). To obtain the intrinsic tissue property μa(r, λ), there is a need to compensate for the extrinsic quantity F(r). Since the optical parameters for tissues are usually heterogeneous and unknown, the fluence varies from case to case and is difficult to model. As a consequence, compensating for the fluence in quantitative in vivo photoacoustic tomography is a challenge.
Currently, fluence compensation can be done invasively or non-invasively. The invasive method includes positioning an optical absorber with a known spectrum close to the region of interest, and then normalizing the measured photoacoustic signals from the objects with the amplitudes of the photoacoustic signals from the optical absorber having a known spectrum. Alternatively, the incident fluence may be estimated by measuring light transmission through a sample of excised tissue of the same type as the tissue overlaying a region of interest. The non-invasive method involves solving the radiative transfer equation (RTE) and photoacoustic wave equations with iterative algorithms. Both of these methods are based on the linear relationship between the local fluence incident on the blood vessel, acquired either by experimental measurements or numerical simulation, and the peak amplitude of the photoacoustic signal.
A temporal profile of the photoacoustic signal is used to quantify the optical absorption coefficients. This temporal method is self-calibrating since it depends on the relative temporal profile rather than the absolute amplitude of the photoacoustic signal. Therefore, the temporal profile is less dependent on changes in the optical properties of overlying tissues. However, since the temporal profile is distorted by various factors, such as the limited bandwidth of acoustic detectors and frequency-dependent acoustic attenuation of the sample or region of interest, linear translation of the temporal profile to the optical absorption may be inaccurate.
As will be described in greater detail below, embodiments of the present disclosure provide a method for quantifying the optical absorption coefficient using the acoustic spectra of the photoacoustic signals. By carefully investigating the factors involved in generating the acoustic spectrum, the effects of detector bandwidth and acoustic attenuation are eliminated, as shown later. At least some embodiments of the method are self-calibrating since it deals only with the relative change in various acoustic frequencies. The acoustic spectrum S(ω) of the received photoacoustic signal depends on three factors: 1) the ‘real’ object spectrum O(ω, λ) measured with unit fluence, which is related to the target object's shape, size, optical properties, and fluence incident directions; 2) the system dependent response H(ω), which is the Fourier transform of the photoacoustic signal from an ideal point absorber measured with this system without acoustic attenuation in the tissue; and 3) the tissue related acoustic attenuation effect a(ω), which is related to the acoustic properties of the tissue that lies between the target objects (e.g., region of interest) and the detector. Based on a system linearity assumption, the above factors lead to the expression S(ω, λ)=F(λ)O(ω, λ)H(ω)a(ω). The last two terms remain unchanged when samples are measured with the same system under the same condition, and therefore are usually cancellable. An example is where light at various optical wavelengths is used to excite one blood vessel. Therefore, by dividing the photoacoustic acoustic spectrum measured at one optical wavelength by the spectrum measured at another wavelength, the system dependent effects and the acoustic attenuation effect may be eliminated. As such, the absolute value of μa may be quantified with this method even though F varies with the optical wavelength. By contrast, previous methods can quantify only the relative value of μa.
As a feasibility study, this idea has been implemented and validated using one form of photoacoustic imaging, optical resolution photoacoustic microscopy (OR-PAM), where the object spectrum O(ω, λ) may be relatively easily modeled. Referring to FIG. 1, in OR-PAM, photoacoustic A-scan signals are acquired through time-resolved acoustic detection, and three-dimensional images are formed by raster scanning the acoustic transducer 102 along the transverse plane, as shown in FIG. 1. The axial resolution of the system depended on the transducer bandwidth (centered at 50 MHz with 80% bandwidth), while the lateral resolution relied on optical focusing, which can reach the theoretical optical diffraction limit. For the system of FIG. 1, the axial and lateral resolutions were quantified to be approximately 15 micrometers (μm) and approximately 5 μm, respectively. Therefore, the surface of blood vessels with a diameter of greater than approximately 30 μm may be roughly treated as a flat surface. In this case, the acoustic spectrum of the generated photoacoustic signal is only related to the optical penetration depth. If F0 is used to denote the incident fluence on the surface of the blood vessel, the fluence inside the blood vessel obeys Beer's law and may be written as F(z)=F0 exp(−μaz), as shown in FIG. 2. Here, the reduced scattering coefficient is much less than the absorption coefficient, because the anisotropy factor is so close to 1 in blood in the optical spectral region that was used (around 585 nm); therefore, it is neglected. The photoacoustic signal generated by the object is expressed using Equation (1):

O(t, λ)=μa(λ)exp [−μa(λ)ct]  Eq. (1)

where c is the speed of sound in the biological tissue. Performing a Fourier transformation of Eq. (1) leads to Equation (2):





  
    
      
        
          
          
            O
            ⁡
            
              (
              
                ω
                ,
                λ
              
              )
            
          
          
        
        =
        
          1
          
            
              
                
                  [
                  
                    ω
                    /
                    
                      
                        μ
                        a
                      
                      ⁡
                      
                        (
                        λ
                        )
                      
                    
                  
                  ]
                
                2
              
              +
              
                c
                2
              
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        
          (
          2
          )
        
      
    
  




Moreover, if the photoacoustic signals of the blood vessel are measured at two optical wavelengths, the ratio of the spectra of the photoacoustic signals may be written as shown in Equation (3):





  
    
      
        
          
            S
            ⁡
            
              (
              
                ω
                ,
                
                  λ
                  1
                
              
              )
            
          
          
            S
            ⁡
            
              (
              
                ω
                ,
                
                  λ
                  2
                
              
              )
            
          
        
        =
        
          
            
              
                F
                ⁡
                
                  (
                  
                    λ
                    1
                  
                  )
                
              
              ⁢
              
                O
                ⁡
                
                  (
                  
                    ω
                    ,
                    
                      λ
                      1
                    
                  
                  )
                
              
              ⁢
              
                H
                ⁡
                
                  (
                  ω
                  )
                
              
              ⁢
              
                a
                ⁡
                
                  (
                  ω
                  )
                
              
            
            
              
                F
                ⁡
                
                  (
                  
                    λ
                    2
                  
                  )
                
              
              ⁢
              
                O
                ⁡
                
                  (
                  
                    ω
                    ,
                    
                      λ
                      2
                    
                  
                  )
                
              
              ⁢
              
                H
                ⁡
                
                  (
                  ω
                  )
                
              
              ⁢
              
                a
                ⁡
                
                  (
                  ω
                  )
                
              
            
          
          =
          
            
              
                F
                ⁡
                
                  (
                  
                    λ
                    1
                  
                  )
                
              
              ⁢
              
                
                  
                    
                      [
                      
                        ω
                        /
                        
                          
                            μ
                            a
                          
                          ⁡
                          
                            (
                            
                              λ
                              
                                2
                                ⁢
                                
                                  
                                
                              
                            
                            )
                          
                        
                      
                      ]
                    
                    2
                  
                  +
                  
                    c
                    2
                  
                
              
            
            
              
                F
                ⁡
                
                  (
                  
                    λ
                    2
                  
                  )
                
              
              ⁢
              
                
                  
                    
                      [
                      
                        ω
                        /
                        
                          
                            μ
                            a
                          
                          ⁡
                          
                            (
                            
                              λ
                              1
                            
                            )
                          
                        
                      
                      ]
                    
                    2
                  
                  +
                  
                    c
                    2
                  
                
              
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        
          (
          3
          )
        
      
    
  




Furthermore, the absolute values of μa(λ1), μa(λ2) and F(λ1)/F(λ2) may be derived by fitting the ratio of Eq. (3). The assumption of laser beam collimation within the blood vessel for Eqs. (1)-(3) is valid within and far from the optical focal zone. Otherwise, the 3D Gaussian beam profile should be considered. Here, the Rayleigh range of the Gaussian beam is approximately 30 μm, which is comparable to the optical penetration depth in blood around 585 nm.

FIG. 1 is a schematic diagram of an exemplary OR-PAM system, FIG. 2 is a diagram of an exemplary experimental setup, and FIG. 3 is a set of graphs that illustrate an example of O(ω, λ), H(ω), and a(ω). The system includes a dye laser 104 pumped by a Nd: YLF laser 106 that is used as the irradiation source. The laser beam from the dye laser 104 passes through a condenser lens 140, is spatially filtered by a pinhole 108, is redirected by a mirror 142, and is then focused by an objective lens 110. Sonic and ultrasonic focusing is achieved through a plano-concave lens 112. The optical objective lens 110 and 50 MHz acoustic transducer 102 are confocally configured via a correction lens 122, a right angle prism 124, and a silicone oil layer 126. A water tank 128 and membrane 130 separate the acoustic lens 112 from the sample 132. Volumetric images are generated by data acquisition device 114 through a combination of time-resolved detection of the photoacoustic waves with a two-dimensional raster scanning in the transverse plane. A scanner controller 116 and scanner 118 (i.e., two dimensional motivator platform) provide rasterizing under the control of data acquisition device 114, and the data acquisition device 114 receives data from acoustic transducer 102 via amplifier 120. The data acquisition device 114 provides a trigger signal to the pump laser 106 and the scanner controller 116 and also provides a clock signal to the scanner controller 116. The trigger and clock signals initiate and control the timing of light pulses from the pump laser 106 and dye laser 104 in coordination with movement of the scanner 118.
As shown in FIG. 2, the optical focus 202 of the system of FIG. 1 is much smaller than a targeted blood vessel whose top surface 204 within the optical focal diameter can therefore be approximated as a plane. An overlying tissue 206 covers the blood vessel top surface 204, and blood 208 lies beneath the blood vessel top surface 204 in the setup of FIGS. 1 and 2. The optical fluence within the blood vessel decays exponentially with depth at a rate of the optical absorption coefficient. FIG. 3 is an example of the object spectrum O(ω, λ), the system dependent response H(ω), and the tissue related acoustic attenuation effect a(ω).
In a phantom study, original black ink was diluted with water in six ratios ranging from 1:1 to 1:6. The original and diluted ink samples were sequentially placed in a container, sealed with plastic membrane, and then the container was placed in a water tank. Photoacoustic A-line signals were acquired from these samples. The acoustic spectra of the photoacoustic signals are shown in FIG. 4. Compared with the spectrum of the photoacoustic signal from the original ink sample, the spectra of the photoacoustic signals from the diluted ink samples are “shifted” to lower frequencies. Light penetrated deeper in lower concentration ink samples, and the corresponding photoacoustic signals decay more slowly in the time domain. Therefore, the spectra contain more low-frequency components. By dividing the measured spectra of any two ink samples frequency-by-frequency, as shown in FIG. 5, the absorption coefficients of both samples may be determined by fitting the resultant ratio curve with Eq. (3). Because parts of the spectra (e.g., 0-5 MHz and greater than 70 MHz in the exemplary system) are unreliable due to the limited-band detection, they are not used for the fitting. The quantified absorption coefficients of all seven samples and their theoretical values are plotted in FIG. 6. To demonstrate whether the recovered absorption coefficients are independent of acoustic attenuation and optical fluence, three of the ink samples were covered with an identical layer of optical phantom ( ̃1.5 mm 2% Agar, 0.1% intralipid, 1% black ink). The spectra of the photoacoustic signals from one ink sample with and without this layer are shown in FIG. 7. The spectral profiles differ because of the acoustic attenuation, while the spectral magnitudes differ owing to the optical fluence attenuation. Since the acoustic properties of the layer added between the samples and the detector are the same for the three ink samples, the acoustic attenuation may be cancelled by taking the ratio of the acoustic spectra of photoacoustic signals from any two covered ink samples. The quantified absorption coefficients of these samples agree with the expected values as shown in FIG. 6.
In an in vivo experiment, a region measuring approximately 1 square millimeter (mm2) in a nude mouse ear was imaged with two optical wavelengths of 561 nanometers (nm) and 570 nm. FIG. 8 shows the photoacoustic maximum amplitude projection (MAP) image acquired with an optical wavelength of 570 nm, which is an oxygen insensitive absorption wavelength of hemoglobin. Each point in the MAP image records the maximum value of a Hilbert transformed photoacoustic A-scan. Two vessels marked with V1 and V2 in FIG. 8 were selected for a quantitative study. The A-scans acquired within these two vessels were properly aligned and then averaged. For each vessel, the acoustic spectrum measured at 561 nm was divided point by point by the acoustic spectrum measured at 570 nm, and the absorption coefficients were acquired by fitting the ratio with Eq. (3), as shown in FIG. 9. The [HbT], [HbO2], and [HbR], together with the sO2 values were calculated based on the quantified optical absorption coefficients at the two optical wavelengths, as shown below in Eqs.(4) and (5) and Table 1.




  
    
      
        
          [
          
            
              
                
                  [
                  
                    HbO
                    2
                  
                  ]
                
              
            
            
              
                
                  [
                  HbR
                  ]
                
              
            
          
          ]
        
        =
        
          
            
              
                64500
                2.303
              
              ⁡
              
                [
                
                  
                    
                      
                        
                          ɛ
                          ox
                        
                        ⁡
                        
                          (
                          
                            λ
                            1
                          
                          )
                        
                      
                    
                    
                      
                        
                          ɛ
                          de
                        
                        ⁡
                        
                          (
                          
                            λ
                            1
                          
                          )
                        
                      
                    
                  
                  
                    
                      
                        
                          ɛ
                          ox
                        
                        ⁡
                        
                          (
                          
                            λ
                            2
                          
                          )
                        
                      
                    
                    
                      
                        
                          ɛ
                          de
                        
                        ⁡
                        
                          (
                          
                            λ
                            2
                          
                          )
                        
                      
                    
                  
                
                ]
              
            
            
              -
              1
            
          
          ⁡
          
            [
            
              
                
                  
                    
                      μ
                      a
                    
                    ⁡
                    
                      (
                      
                        λ
                        1
                      
                      )
                    
                  
                
              
              
                
                  
                    
                      μ
                      
                        a
                        ⁢
                        
                          
                        
                      
                    
                    ⁡
                    
                      (
                      
                        λ
                        2
                      
                      )
                    
                  
                
              
            
            ]
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        
          (
          4
          )
        
      
    
  
  
    
      
        
          sO
          2
        
        =
        
          
            [
            
              HbO
              2
            
            ]
          
          
            
              [
              
                HbO
                2
              
              ]
            
            +
            
              [
              HbR
              ]
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        
          (
          5
          )
        
      
    
  




In Eq.(4), λ1 and λ2 are the two wavelengths, and εox and εde are the known molar extinction coefficients of oxy- and deoxyhemoglobin, respectively. According to the sO2 values, V1 and V2 were identified to be an arteriole-venule pair. The incident fluence ratio at the two optical wavelengths F(λ1)/F(λ2) was also quantified for both vessels. In this special case, the two vessels are embedded at a similar depth, and the optical and the acoustic properties of the overlying tissue are comparable. Here, the quantified fluence ratios turned out to be the same for V1 and V2. If the wavelength-dependent fluence variations are ignored by simply assuming F(λ1)/F(λ2)=1, as shown in the last column of Table 1, the quantified sO2 values become inaccurate by approximately 8% and 11% for the artery and the vein, respectively.
















TABLE 1







μa(λ1)
μa(λ2)

[HbT]
[HbO2]
[HbR]

sO2



(cm−1)
(cm−1)
F(λ2)/F(λ1)
(g/L)
(g/L)
(g/L)
sO2
F2/F1 = 1
























V1
143 ± 3
188 ± 4
0.96 ± 0.01
110.6 ± 8.1
106.2 ± 4.3
 4.4 ± 3.8
0.96 ± 0.04
0.88


(Artery)


V2
159 ± 4
186 ± 5
0.96 ± 0.01
110.2 ± 9.2
 77.1 ± 4.9
33.1 ± 4.3
0.70 ± 0.07
0.62


(Vein)









Embodiments of the method may be applied in the optical diffusive regime with photoacoustic microscopy, whose lateral resolution relies on the acoustic focus (e.g.,  ̃50 μm at 50 MHz acoustic frequency). As shown in FIG. 10, photoacoustic A-scan signals are acquired through time-resolved acoustic detection, and three-dimensional images are formed by raster scanning an acoustic transducer 402 and optical components (e.g., optical illumination source 404 and mirror 406) along a transverse plane (e.g., transverse to the direction of the optical illumination striking a surface 408 of a sample 410). In the embodiment of FIG. 10, the surface 408 of the sample 410 is illuminated with an annular ring having a dark center 412 as shown in the call out from the main figure. The surface of blood vessels with sufficiently large diameter (e.g., greater than 300 μm for the 50-μm lateral resolution) may be approximately treated as a flat surface. As shown in FIG. 11, it may generally be assumed that the light is completely diffused (e.g., see diffused light 504) when it reaches a top surface 502 of a blood vessel within tissue 410. Blood 506 is located below the blood vessel top surface 502, overlying tissue 508 is located above the blood vessel top surface 502, and the acoustic focus 510 narrows at the surface of the blood vessel 502 as shown in FIG. 11. Moreover, as shown in FIG. 12, light at each point (e.g., point 602) on the top surface of the blood vessel 502 can be seen as an isotropic point source for each point (e.g., point 604) in a thin layer of blood 606, and the fluence F(z) in the blood vessel can be expressed as shown in Equation (6):




  
    
      
        
          F
          ⁡
          
            (
            z
            )
          
        
        =
        
          
            π
            ⁢
            
              
            
            ⁢
            
              R
              2
            
            ⁢
            
              ∫
              
                
                  ∫
                  
                    S
                    ′
                  
                
                ⁢
                
                  
                    exp
                    ⁡
                    
                      (
                      
                        
                          -
                          
                            μ
                            a
                          
                        
                        ⁢
                        
                          
                            
                              x
                              ′2
                            
                            +
                            
                              y
                              ′2
                            
                            +
                            
                              z
                              ′2
                            
                          
                        
                      
                      )
                    
                  
                  ⁢
                  
                    d
                    
                      x
                      ′
                    
                  
                  ⁢
                  
                    d
                    
                      y
                      ′
                    
                  
                
              
            
          
          =
          
            
              F
              0
            
            ⁢
            
              π
              2
            
            ⁢
            
              R
              2
            
            ⁢
            
              
                ∫
                0
                ∞
              
              ⁢
              
                
                  exp
                  ⁡
                  
                    (
                    
                      
                        -
                        
                          μ
                          a
                        
                      
                      ⁢
                      
                        
                          u
                          +
                          
                            z
                            2
                          
                        
                      
                    
                    )
                  
                
                ⁢
                
                  
                
                ⁢
                
                  d
                  u
                
              
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        
          (
          6
          )
        
      
    
  




where R is the radius of the acoustic focus 510. The object spectrum O(ω, λ) may then be calculated using Fourier transformation, where z is first converted to time t through z=ct. Similar to the case of OR-PAM, the system-dependent response H(ω) and the tissue related acoustic attenuation effect a(ω) are canceled by dividing, at each acoustic frequency, the acoustic spectra measured at two optical wavelengths. FIG. 13 shows another embodiment of this method with an acoustic transducer array 702. With a synthetic aperture focusing technique (SAFT), a virtual acoustic focus704 may be formed at any position in the field of view by applying proper time delays to each element of the array system. FIG. 13 includes diffused light 706 on the surface of blood vessel 710 within tissue 708. Using the same process as described above, the optical absorption coefficients may be extracted from the acoustic spectra measured at multiple optical wavelengths. Although a linear array is illustrated in the figure, arrays of other shapes—such as a semicircle, a hemisphere, and a 2D plane—can be used as well.

To increase the accuracy of this method, an acoustic transducer with an appropriate bandwidth is selected. The acoustic spectrum of the photoacoustic signal is related to the light penetration depth. Therefore, the central frequency of the transducer should vary with the penetration depth to achieve the best signal to noise ratio (SNR). Therefore, the central frequency of the transducer should match the penetration depth to maximize signal to noise ratio (SNR). Moreover, SNR is usually low at high acoustic frequencies due to acoustic attenuation. O(ω, λ), H(ω), and a(ω) are all band-limited, and in an exemplary embodiment, H(ω) is chosen to match O(ω, λ) and a(ω).
Referring to FIG. 14, a method of determining an absorption coefficient of a volume of an object begins with a focusing element (e.g., an objective or lens) focuses a first light pulse from a light source (e.g., a laser) into a volume of an object or sample at 802. At 804, a transducer detects a first acoustic spectrum emitted by the object in response to the volume of the object receiving the first light pulse. At 806, the focusing element focuses a second light pulse from the light source into the volume of the object, and at 808, the transducer detects a second acoustic spectrum emitted by the object in response to the volume of the object receiving the second light pulse. In operation, the transducer receives an acoustic wave from the object and performs a Fourier transform to generate the acoustic spectrum. It is contemplated that the Fourier transform may be performed by a component separate from the transducer component receiving the acoustic waveform. At 810, a controller determines an absorption coefficient from the detected first and second acoustic spectra by dividing the first acoustic spectrum by the second acoustic spectrum frequency by frequency yielding a normalized acoustic spectrum. The normalized acoustic spectrum is fitted to an ideal curve of an absorption coefficient to determine the absorption coefficient of the volume of the object. The first and second light pulses have different wavelengths, and the method may be repeated for multiple volumes of the object and to obtain an image of the object.
It will be understood that the particular embodiments described herein are shown by way of illustration and not as limitations of the disclosure. The principal features of this disclosure may be employed in various embodiments without departing from the scope of the disclosure. Those of ordinary skill in the art will recognize numerous equivalents to the specific procedures described herein. Such equivalents are considered to be within the scope of this disclosure and are covered by the claims.
All of the compositions and/or methods disclosed and claimed herein may be made and/or executed without undue experimentation in light of the present disclosure. While the compositions and methods of this disclosure have been described in terms of the embodiments included herein, it will be apparent to those of ordinary skill in the art that variations may be applied to the compositions and/or methods and in the steps or in the sequence of steps of the method described herein without departing from the concept, spirit, and scope of the disclosure. All such similar substitutes and modifications apparent to those skilled in the art are deemed to be within the spirit, scope, and concept of the disclosure as defined by the appended claims.
It will be understood by those of skill in the art that information and signals may be represented using any of a variety of different technologies and techniques (e.g., data, instructions, commands, information, signals, bits, symbols, and chips may be represented by voltages, currents, electromagnetic waves, magnetic fields or particles, optical fields or particles, or any combination thereof). Likewise, the various illustrative logical blocks, modules, circuits, and algorithm steps described herein may be implemented as electronic hardware, computer software, or combinations of both, depending on the application and functionality. Moreover, the various logical blocks, modules, and circuits described herein may be implemented or performed with a general purpose processor (e.g., microprocessor, conventional processor, controller, microcontroller, state machine or combination of computing devices), a digital signal processor (“DSP”), an application specific integrated circuit (“ASIC”), a field programmable gate array (“FPGA”) or other programmable logic device, discrete gate or transistor logic, discrete hardware components, or any combination thereof designed to perform the functions described herein. Similarly, steps of a method or process described herein may be embodied directly in hardware, in a software module executed by a processor, or in a combination of the two. A software module may reside in RAM memory, flash memory, ROM memory, EPROM memory, EEPROM memory, registers, hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art. Although preferred embodiments of the present disclosure have been described in detail, it will be understood by those skilled in the art that various modifications can be made therein without departing from the spirit and scope of the disclosure as set forth in the appended claims.
A controller, computing device, or computer, such as described herein, includes at least one or more processors or processing units and a system memory. The controller typically also includes at least some form of computer readable media. By way of example and not limitation, computer readable media may include computer storage media and communication media. Computer storage media may include volatile and nonvolatile, removable and non-removable media implemented in any method or technology that enables storage of information, such as computer readable instructions, data structures, program modules, or other data. Communication media typically embody computer readable instructions, data structures, program modules, or other data in a modulated data signal such as a carrier wave or other transport mechanism and include any information delivery media. Those skilled in the art should be familiar with the modulated data signal, which has one or more of its characteristics set or changed in such a manner as to encode information in the signal. Combinations of any of the above are also included within the scope of computer readable media.
This written description uses examples to disclose the disclosure, including the best mode, and also to enable any person skilled in the art to practice the disclosure, including making and using any devices or systems and performing any incorporated methods. The patentable scope of the disclosure is defined by the claims, and may include other examples that occur to those skilled in the art. Such other examples are intended to be within the scope of the claims if they have structural elements that do not differ from the literal language of the claims, or if they include equivalent structural elements with insubstantial differences from the literal language of the claims.

## Claims

1. A method of quantifying an optical absorption coefficient of an object, said method comprising:
directing a first light pulse emitted by a light source into a volume of the object, wherein the first light pulse has a first wavelength;
receiving at an acoustic transducer a first acoustic wave emitted by the object in response to receiving the first light pulse;
performing on a controller a Fourier transformation on the received first acoustic wave to generate a first acoustic spectrum;
directing a second light pulse emitted by the light source into the volume of the object, wherein the second light pulse has a second wavelength;
receiving at the acoustic transducer a second acoustic wave emitted by the object in response to receiving the second light pulse;
performing on the controller a Fourier transformation on the received second acoustic wave to generate a second acoustic spectrum; and
quantifying the absorption coefficient of the volume of the object, wherein the absorption coefficient is a function of the first acoustic spectrum and the second acoustic spectrum.

2. The method of claim 1 wherein quantifying the absorption coefficient of the volume of the object comprises dividing the first acoustic spectrum by the second acoustic spectrum frequency by frequency to generate a normalized acoustic spectrum, and fitting the normalized acoustic spectrum to an absorption coefficient curve.

3. The method of claim 1 wherein the light source comprises a first light source for emitting the first light pulse and a second light source for emitting the second light pulse.

4. The method of claim 1 wherein the first wavelength is different from the second wavelength.

5. The method of claim 1 further comprising disregarding a portion of the first acoustic spectrum and the second acoustic spectrum outside of a system bandwidth, wherein the system bandwidth is at least partially determined as a function of a bandwidth of the acoustic transducer.

6. The method of claim 1 further comprising determining at least one of HbT, HbO2, HbR, and sO2 values as a function of the quantified absorption coefficient.

7. A system for quantifying an optical absorption coefficient of a volume of an object, said system comprising:
a light source for emitting a first light pulse having a first wavelength and a second light pulse having a second wavelength;
an optical element for directing the first light pulse and the second light pulse into a volume of the object;
an acoustic transducer for receiving a first acoustic wave emitted by the object in response to the volume of the object absorbing the first light pulse and a second acoustic wave emitted by the object in response to the volume of the object absorbing the second light pulse; and
a controller for performing a Fourier transformation on the received first acoustic wave and second acoustic wave to generate a first acoustic spectrum and a second acoustic spectrum, respectively, and quantifying the absorption coefficient of the volume of the object as a function of the first acoustic spectrum and the second acoustic spectrum.

8. The system of claim 7 wherein the controller is further configured to quantify the absorption coefficient of the volume of the object by dividing the first acoustic spectrum by the second acoustic spectrum, frequency by frequency, to generate a normalized acoustic spectrum and to fit the normalized acoustic spectrum to an absorption coefficient curve.

9. The system of claim 7 wherein the light source comprises a first light source for emitting the first light pulse and a second light source for emitting the second light pulse, and wherein the first wavelength is different from the second wavelength.

10. The system of claim 7 wherein the controller is further configured to disregard a portion of the first acoustic spectrum and the second acoustic spectrum outside of a system bandwidth, wherein said system bandwidth is at least partially determined as a function of a bandwidth of the acoustic transducer.

11. The system of claim 7 wherein the controller is configured to quantify at least one of HbT, HbO2, HbR, and sO2 values as a function of the quantified absorption coefficient.

12. A system for generating an image of an object and quantifying an optical absorption coefficient of a volume of an object, said system comprising:
a light source for emitting a first light pulse having a first wavelength and a second light pulse having a second wavelength;
an optical element for directing the first light pulse and the second light pulse into a volume of a plurality of volumes of the object;
an acoustic transducer for receiving a first acoustic wave produced by the object in response to the volume of the object absorbing the first light pulse and a second acoustic wave produced by the object in response to the volume of the object absorbing the second light pulse;
a controller for performing a Fourier transformation on the received first acoustic wave and second acoustic wave to generate a first acoustic spectrum and a second acoustic spectrum, respectively, quantifying the absorption coefficient of the volume of the object as a function of the first acoustic spectrum and the second acoustic spectrum and for reconstructing an image from acoustic spectra detected from the plurality of volumes of the object; and
a scanner for altering a spatial relationship between the object and the optical element after the first acoustic wave and the second acoustic wave are detected by the acoustic transducer for the volume such that the first light pulse and the second light pulse are directed into another volume of the plurality of volumes of the object.

13. The system of claim 12 wherein the controller is configured to quantify the absorption coefficient of the volume of the object by dividing the first acoustic spectrum by the second acoustic spectrum, frequency by frequency, to generate a normalized acoustic spectrum and fits the normalized acoustic spectrum to an absorption coefficient curve.

14. The system of claim 12 wherein the light source comprises a first light source for emitting the first light pulse and a second light source for emitting the second light pulse, and wherein the first wavelength is different from the second wavelength.

15. The system of claim 12 wherein the controller is configured to disregard a portion of the first acoustic spectrum and the second acoustic spectrum outside of a system bandwidth, wherein said system bandwidth is at least partially determined as a function of a bandwidth of the transducer.

16. The system of claim 12 wherein the controller is configured to quantify at least one of HbT, HbO2, HbR, and sO2 values as a function of the quantified absorption coefficient.