Login| Sign Up| Help| Contact|

Patent Searching and Data


Title:
COHERENCE GATED PHOTOACOUSTIC REMOTE SENSING (CG-PARS)
Document Type and Number:
WIPO Patent Application WO/2019/145764
Kind Code:
A1
Abstract:
A coherence gated photoacoustic remote sensing system for imaging a subsurface structure in a sample with optical resolution may include an excitation beam source configured to generate an excitation beam that induces ultrasonic signals in the sample at an excitation location; an interrogation team source configured to generate an interrogation team incident on the sample at an interrogation location, a portion of the interrogation beam returning from the sample that is indicative of the generated ultrasonic signals, the interrogation beam being a low-coherent beam; an optical system that focuses the excitation beam onto the sample at an excitation location, and focuses the interrogation beam onto the sample at an interrogation location, at least the interrogation location being below the surface of and within the sample; and a low coherence interferometer that isolates a returning portion of the interrogation beam that corresponds to an interrogation event of the sample.

Inventors:
ZEMP ROGER (CA)
HAJI REZA PARSIN (CA)
BELL KEVAN (CA)
Application Number:
PCT/IB2018/057585
Publication Date:
August 01, 2019
Filing Date:
September 28, 2018
Export Citation:
Click for automatic bibliography generation   Help
Assignee:
ILLUMISONICS INC (CA)
International Classes:
G01N21/17; A61B5/00; G01B9/02; G01N21/47
Domestic Patent References:
WO2009055705A22009-04-30
WO2014036405A22014-03-06
Foreign References:
DE102010012809A12011-09-29
CN103048271A2013-04-17
US20140185055A12014-07-03
US20140185055A12014-07-03
US20120200845A12012-08-09
US20160113507A12016-04-28
US20170215738A12017-08-03
US20140275942A12014-09-18
US9335253B22016-05-10
Other References:
CEDRIC BLATTER ET AL: "Intrasweep phase-sensitive optical coherence tomography for noncontact optical photoacoustic imaging", OPTICS LETTERS, OPTICAL SOCIETY OF AMERICA, US, vol. 37, no. 21, 1 November 2012 (2012-11-01), pages 4368 - 4370, XP001579227, ISSN: 0146-9592, [retrieved on 20121017], DOI: 10.1364/OL.37.004368
ZHIHUA DING ET AL: "Real-time phase-resolved optical coherence tomography and optical Doppler tomography", OPTICS EXPRESS, vol. 10, no. 5, 11 March 2002 (2002-03-11), pages 236 - 1181, XP055532518, DOI: 10.1364/OE.10.000236
KEVAN L. BELL ET AL: "Coherence-gated photoacoustic remote sensing microscopy", OPTICS EXPRESS, vol. 26, no. 18, 29 August 2018 (2018-08-29), pages 23689, XP055531971, DOI: 10.1364/OE.26.023689
HAJI REZA ET AL., LIGHT: SCIENCE & APPLICATIONS, vol. 6, 2017, pages 16278
Download PDF:
Claims:
CLAIMS

1. A coherence gated photoacoustic remote sensing system for imaging a subsurface structure in a sample with optical resolution, comprising:

an excitation beam source configured to generate an excitation beam that induces ultrasonic signals in tire sample at an excitation location;

an interrogation beam source configured to generate an interrogation beam incident on the sample at an interrogation location, a portion of the interrogation beam returning from the sample that is indicative of the generated ultrasonic signals, the interrogation beam being a low-coherent beam;

an optical system that focuses the excitation beam onto the sample at an excitation location, and focuses the interrogation beam onto the sample at an interrogation location, at least the interrogation location being below the surface of and within the sample; and

a low coherence interferometer that isolates a returning portion of the interrogation team that corresponds to an interrogation event of the sample.

2. The system of claim 1 , further comprising a reference beam source configured to generate a reference beam that travels along a reference path, and wherein the lo " coherence interferometer isolates the returning portion using the reference beam.

3. Tire system of claim 2, wherein the reference beam source is configured to generate one or more additional reference beams that ar e phase shifted relative to the reference team, and wherein the low coherence interferometer isolates the returning portion using the reference beam and the one or more additional reference beams.

4. The system of claim 3, wherein the one or more additional reference beams are phased shifted by at least one of a different path length, one or more wave plates, and one or more circulators.

5. The system of claim 3, wherein the one or more additional reference beams are defected either in parallel or serially with the reference beam.

6. The system of claim 1, wherein the excitation beam and the interrogation team are pulsed or intensity-modulated.

7. The system of claim 1 , wherein the excitation location and the interrogation location are each below the surface of and within the sample.

8. The system of claim 1, wherein at least one of the excitation location and the interrogation location are within 1 mm of the surface of the sample.

9. The system of claim 1, wherein at least one of the excitation location and the interrogation location are gr eater than 1 pm below the surface of the sample.

10. The system of claim 1 , wherein tire excitation location and the interrogation location are focal points that are at least partially overlapping.

11. The system of claim 1. further comprising a processes: that calculates an image of the sample based on the returning portion of the interrogation beam.

12. The system of claim 1 , wherein the interrogation beam has pulses that are sufficiently short that acoustic propagation is negligible.

13. The system of claim 1, wherein, for each detection location, the system applies an excitation beam with more than one frequency, bandwidth, phase shift or combination thereof

14. Tire system of claim 1, wherein the optical system interrogates each interrogation location the sample in a non-excited state and after an excitation beam excites the sample.

15. The system of claim 1, wherein the excitation beam source is configured to generate one or more excitation beams that excites the sample with a plurality of frequencies, a plurality of bandwidths or combinations thereof.

16. A method of using the system of claim 1 during:

functional imaging during brain surgery1;

assessing internal bleeding and cauterization verification;

imaging perfusion sufficiency of organs and organ transplants;

imaging angiogenesis around islet transplants;

imaging of skin-grafts; imaging of tissue scaffolds and biomaterials to evaluate vascularization and/or immune rejection;

imaging to aid microsurgery; or

procedures for guidance to avoid cuting critical blood vessels and nerves.

17. A method of using the system of claim 1 in combination with fluorescence microscopy two-photon and confoeal fluorescence microscopy. Coherent- Raman-

Stokes microscopy, Raman microscopy, or Optical coherence tomography.

18. The method of claim 17, further including performing microcirculation imaging or performing blood oxygenation parameter imaging with the system.

19. An endoscope comprising the system of claim 1.

20. A surgical microscope comprising the system of claim 1.

21. A method of remote sensing a sample, comprising the steps of:

providing a coherence gated photoacoustic remote sensing system comprising an excitation beam and an interrogation beam, the interrogation beam being a low-coherent beam:,

causing the excitation beam to induce ultrasonic signals in the sample at an excitation location;

causing the interrogation beam to interrogate the sample at an interrogation location, wherein a portion of the interrogation beam returns from the sample that is indicative of the generated ultrasonic signals, the interrogation location being below the surface of and within the sample;

using a low coherence interferometer to isolate a. returning portion of the interrogation beam to achieve an interrogation event of the sample.

22. Tire method of claim 21. further comprising providing a reference beam that travels along a reference path and wherein the low coherence interferometer isolates the returning portion using the reference beam.

23. The method of claim 22, further comprising the step of prowling one or more additional reference beams that are phase shifted relative to the reference beam, and wherein the low coherence interferometer isolates the returning portion using the reference beam and the one or more additional reference beams.

24. The method of claim 23, wherein the one or more additional reference beams are phased shifted by at least one of a different path length, one or more wave plates, and one or more circulators.

25. The method of claim 23, wherein tire one or more additional reference beams are detected either hi parallel or serially with the reference beam.

26. The method of claim 21 , wherein the excitation beam and the interrogation beam ar e pulsed or intensity-modulated

27. The method of claim 21, wherein the excitation location and the interrogation location are each below the surface of and within the sample.

28. The method of claim 21, wherein at least one of the excitation location and the interrogation location are within 1 mm of the surface of the sample.

29. The method of claim 21, wherein at least one of the excitation location and the interrogation location are greater than 1 p below the surface of the sample.

30. Tiie method of claim 21 , wherein the excitation location and the interrogation location are focal points that are at least partially overlapping.

31. The method of claim 21 , further comprising the step of calculating an image of the sample based on the returning portion of the interrogation beam.

32. The method of claim 21, wherein the interrogation beam has pulses that are sufficiently short that acoustic propagation is negligible.

33. The method of claim 21, wherein, for each detection Location, the excitation beam is operated to provide with more than one frequency, bandwidth phase shift, or combination thereof

34. The method of claim 21, further comprising the step of interrogating each interrogation location in a non-excited state and after the excitation beam excites the sample.

35. The method of claim 21, wherein the excitation beam comprises one or more excitation beams that excites the sample with a plurality of frequencies, a plurality of bandwidths or combinations thereof.

Description:
COHERENCE GATED PHOTOACOUSTIC REMOTE SENSING (CG-PARS)

CROSS-REFERENCE TO RELATED APPLiCATION(S)

[0001] This patent application claims benefit of priority under 35 U.S.C. § 119 to U.S. Provisional Patent Application No. 62/622,816, filed January 26, 2018, the entirety of which is incorporated herein by reference.

FIELD

[0002] This relates to the field of optical imaging and, in particular, to a laser-based method and system for non-contact imaging of biological tissue in vivo, ex vivo, or in vitro.

BACKGROUND

[0003] The entireties of the U.S. Patents and Patent Publications set forth herein are expressly incorporated by reference.

[0004] Photoacoustic imaging is an emerging hybrid imaging technology providing optical contrast with high spatial resolution. Nanosecond or picosecond laser pulses fired into tissue launch thenno-elastic-induced acoustic waves, winch are detected and reconstructed to form high-resolution images.

[0005] Photoacoustic imaging has been developer! into multiple embodiments, primarily including photoacoustic tomography (PAT), photoacoustic microscopy (PAM) which is sometimes referred to as acousiio-resolurion photoacoustic microscopy (AR-PAM), and optical-resolution photoacoustic microscopy (OR-PAM). In PAT, signals are collected from multiple transducer locations and reconstructed to form a tomographic image in a way similar to ultrasound (US) or X-ray computed tomography (CT). One of the differences between PAT and the other two modalities is that an assumption must be made about the sample in order to facilitate reconstruction; typically, tins involves assuming the acoustic propagation velocity within the sample in PAM, typically, a single element focused high-frequency ultrasound transducer is used to collect photoacoustic signals providing acoustic focusing. This transducer, along with the excitation beam may be scanned laterally about the sample to perform volumetric imaging. Both PAT and PAM are typically implemented using an unfocused excitation beam. Both modalities provide acoustic-limited resolution and have penetration depth limited by surface optical exposure limits and acoustic atenuation. OR- PAM, typically, utilizes both optical and acoustic focusing providing further improved resolution (~3um) at further reduced penetration depths (~lmm) now limited by fundamental light transport, that is, the distance which optical focus can be reasonably maintained hi all 7 three embodiments. the acoustic signal is typically collected through an acoustically coupled transducer or other acoustic- or acousto-optic resonator. In all cases the photoacoustic signal can be recorded for various positions to form a 2D or 3D photoacoustic image representing the optical absorption in the sample at the excitation wavelength. The amplitude of the various recorded peaks implies the local optical absorption, and the relative time delay infers the depth from the time required for acoustic propagation.

[0006] Photoacoustic microscopy has shown significant potential for imaging vascular· structures from macro-vessels to micro-vessels. It has also show'll great promise for functional and molecular imaging, including imaging of nanoparticle contrast agents and imaging of gene expression. Multi-wavelength photoacoustic imaging has been used for spectral unmixiiig, such as mapping of blood oxygen saturation, by using known oxy- and deoxy-liemogiobin molar extinction spectra. Since conventional photoacoustic imaging requires acoustic coupling to the sample tire technique is inappropriate for many clinical applications such as wound healing, bum diagnostics, surgery, and many endoscopic procedures. Here, physical contact, coupling, or immersion is undesirable or impractical. Some non-contact photoacoustic detection strategies have been reported.

[0007] However until recently no technique has demonstrated practical non-contact in vivo microscopy in reflection mode with eonfbcal resolution and optical absorption as the contrast mechanism. Most previous approaches detected surface oscillations with interferometric methods wliich have suffered from poor sensitivity' and have been ineffective for high quality' in vivo imaging. One example of a low-coherence interferometry method for sensing photoacoustic signals was proposed in (Gurton et al., US Patent Publication No. 2014/0185055) to be combined with an optical coherence tomography' (OCT) system, resulting in 30pm lateral resolution. Another system is described in (Rousseau et al, US Patent Publication No. 2012/0200845) entitled“Biological Tissue Inspection Method and System”, which describes a noncontact photoacoustic imaging system for in vivo or ex vivo, non-contact imaging of biological tissue without the need for a coupling agent. Other systems use a fiber based interferometer with optical amplification to detect photoacoustic signals and form photoacoustic images of phantoms with acoustic (not optical) resolution. However these systems sutler from a poor signal-fo-noise ratio. Furthermore, in vivo imaging was not demonstrated, and optical-resolution excitation was not demonstrated.

[0008] A recently reported photoacoustic technology known as photoacoustic remote sensing (PARS) microscopy' (Haji Reza et al., US Patent Publication No. 2016.0113507, and Haji Reza et al., US Patent Publication No. 2017/0215738) has been able to solve many of a these sensitivity issues through its detection mechanism. PARS utilizes the elasto-optic effect in which the large photoacoustic initial pressures generate nontrivial modulations in the local refractive ides of a material. By co-focusing a continuous wave interrogation beam with the exci ation spot, the back-reflected time- varying intensity of die interrogation beam encodes information regarding this elasto-optic modulation which in turn implies the magnitude of the generated photoacoustic initial pressure, which is directly related to the local optical absorption in the sample at the excitation spot. PARS has thus far demonstrated improved sensitivity and resolution characteristics over conventional contact-based OR-PAM, with lateral resolutions car-par with confocal microscopy (~60Qnm). However, in some examples, depth sensitivity can be improved. Since PARS may be solely sensitive to the large initial photoacoustic pressures near the excitation spot, time-domain information is not indicative of depth. This may require three dimensional optical scanning when recording 3D volumes. Since PARS has been implemented, in some examples, using a low-coherence superhmhnescent diode (SLD) as the detection source, some advantages may be gained by implementing a low-coherence interferometer.

[0009] Optical coherence tomography (OCT) provides a means of capturing depth- resolved optical scattering information from a sample. This is generally accomplished by the use of low-coherence interferometry. TWO common embodiments of the technique involve a time-domain approach, known as time-domain optical coherence tomography (TD-OCT), and a frequency-domain approach, known as frequency-domain optical coherence tomography (FD-OCT) or spectral-domain optical coherence tomography (SD-OCT). TD-OCT generally is implemented with a single broadband continuous-wave interrogation source which is split into a sample- and reference-path, where the total path length of tire reference-path is scanned such that low-coherence interferometry is performed at various depths along the sample-path. This modality still may necessitate a 3D voxel-based scan for capturing for volumes. SD- OCT generally is implemented with either a broadband source, or a modulated frequency source, where imaging is commonly performed with a fixed reference-path length and depth information is acquired through Fourier transform of the collected spectral data. Here, volumetric scanning only may necessitate lateral scanning as full depth-resolved information is collected with a single acquisition event. There has been a great body of work within the OCT field to provide quantitative optical absorption measurement. This is tire particular interest within the opthalrnie imaging community, which requires oxygen saturation measurement about the fimdus of the eye. While there have been several notable works on tins topic, the current approach is still incapable of direct optical absorption measurement (unlike photoacousfic modalities). Rather optical absorption must be inferred through the use of a visible probe source which can greatly limit the penetration depth into the sample. The resulting OCT image is fit to optical extinction curves providing optical absoiption. It would be beneficial to the biomedical imaging community to offer an improved optical absorption modality.

[0010] There have been several notable attempts to provide a multi-modality implementation of non-PARS-based non-contact photoacoustics and OCT. These include but are not limited to (Wang, US Patent Publication No. 2014/0185055, Johnson et ah, US Patent Publication No. 2014/0275942, and Ode, US Patent No. 9335253). However, all of these works do not provide the same method of operation presented here in that they simply provide separately a non-contact PAT and OCT system. The proposed approach is not to be confused with previous OCT-based photoacoustic detection methods which aimed to detect propagated acoustic waves manifesting themselves as subtle oscillations at the sample outer surface. Instead tire proposed approach locally detects optical-absorptiaii-mdueed initial pressures directly at their sub-surface origins. Additionally, the photoacoustic component of each is specifically analogous to a PAT system in that lateral tomographic reconstruction is required and acoustic-resolution is provided

[0011] Given these complementary properties between PARS and OCT, there would be a dear benefit towar ds augmenting PARS with var ious coherence-gated detection mechanisms. However, for reasons which will be discussed in further sections, a great deal of technical challenges arise with these implementation which are addressed in this disclosure.

SUMMARY

[0012] According to an aspect, there is provided a coherence-gated photoacoustic remote sensing system (CG-PARS) for imaging a subsurface structure in the sample known as coherence-enhanced photoacoustic remote sensing (CEPARS) microscopy, which provides significant axial-resolution characteristics over conventional PARS. Tins may be accomplished through the addition of a tow-coherence interferometer between the sample- path and a newly included reference-path, wherein by virtue of low-coherence interferometry, signals which are associated with path lengths significantly longer or shorter than the reference-path length (when compared with the coherence-length of the broadband interrogation source) are rejected. This may comprise an excitation beam configured to generate ultrasonic signals in the sample-path at an excitation location; an interrogation beam incident on the sample at the excitation location, a portion of the interrogation beam returning from the sample that is indicative of the generated ultrasonic signals; a single reference-path, or multiple reference-paths which may provide various phase offsets, or an optical quadrature detector: an optical combiner to compare the back-reflected sample beam with the single reference, or multiple combiners to compare the back-reflected sample beam with the multiple reference-paths; single, or multiple detectors for collecting the combined beams; and a processing unit for interpreting collected results.

[0013] According to another aspect, there is provided an endoscopic CEPARS which may provide significant axial-resolution characteristics over conventional endoscopic PARS. This may comprise of a fiber optic cable having an input end and a detection end; an excitation beam coupled to the input into the input end of the optical fiber configured to generate ultrasonic signals in the sample-path at an excitation location; an interrogation beam coupled into the input end of the optical fiber incident on the sample at the excitation location, a portion of the interrogation beam returning from the sample back along the optical fiber that is indicative of the generated ultrasonic signals; a single reference-path, or multiple reference-paths, which may provide various phase offsets; an optical combiner to compare the back-reflected sample beam with tire single reference, or multiple combiners to compare the back-reflected sample beam with the multiple reference-paths; single, or multiple detectors for collecting the combined beams; and a processing unit for interpreting collected results.

[0014] According to another aspect, there is provided a CG-PARS system for imaging a subsurface structure in the sample known as spectral-domain coherence-gate photoacoustic remote sensing (SDCG-PARS) microscopy which provides the ability to image Ml depth- resolved optical-absorption within a sample within a single rapid pulse-train drastically improving imaging speeds over conventional PARS, and aforementioned CEPARS. This is accomplished through the addition of a low-coherence interferometer between the sample- path and a reference-path, a detector capable of detecting the spectral content of the combined reference- and sample-paths, and the addition of a rapid (< 100ns) interrogation mechanism such as a pulsed interrogation source, a rapidly modulated continuous-wave (CW) source, photodiode array, rapid shudder, etc. This allows for acquisition of the depth- resolved scattering profile both before, and directly after the sample has undergone photoaeoustie excitation. The difference between these two scattering profiles being indicative of the optical absorption. This may comprise an excitation beam configured to generate ultrasonic signals in the sample-path at an excitation location; an interrogation beam incident on the sample at the excitation location, a portion of the interrogation beam returning from the sample, where the spectrum is indicative of the generated ultrasonic signals; a reference-path which may provide various phase offsets; and optical combiner to compare the back-reflected sample beam with the reference beam; a spectrum detector, which by its own virtue, or virtue of other components is capable of short interrogation times (<lO0ns); and a processing unit for interpreting collected results.

[0015] According to another aspect, there is provided an endoscopic SDCG-PARS which provides fell depth-resolved acquisitions. Ibis comprises a fiber optic cable having an input end and a detection end; an excitation beam coupled to the input into die input end of the optical fiber configured to generate ultrasonic signals in the sample-path at an excitation location: an interrogation beam coupled into the input end of the optical fiber incident on the sample at the excitation location, a portion of the interrogation beam returning from the sample back along the optical fiber where the spectrum is indicative of the generate ultrasonic signals; a reference-path which may provide various phase offsets; and optical combiner to compare the back-reflected sample beam with the reference beam; a spectrum detector, which by its own virtue, or virtue of other components is capable of short interrogation times (<lG0ns); and a processing unit for interpreting collected results.

[0016] For other embodiments of CEPARS and SDCG-PARS, the excitation source may comprise of a single or multiple sources which are pulsed, or CW and modulated. Excitation sources may be narrow-band and may cover a wide range of wavelengths or broadban individually providing wider spectra. This variety of excitation spectral content provides a means of implementing absorption-contrast spectral unmixing of the various target species in a sample. The interrogation source may likewise comprise of a single or multiple sources which are pulsed, or CW and modulated. Interrogation sources may be narrow-band and may cover a wide range of wavelengths or broadband individually providing wider spectra. This variety of interrogation spectral content provides a means of controlling the extinction (thereby the penetration) of the interrogation beam and a means of controlling tire effective coherence-length which dictates the axial-resolution of the device. The optical beam combiner may comprise of an optical coupler such as a beam-splitting cube for bulk optical implementation or a fiber coupler for fiber-based implementation, or some variety of interferometer such as a bulk- or fiber-based Miehekon interferometer common path interferometer (using specially designed interferometer objective lenses), Fizeau interferometer, Ramsey interferometer, Fabry-Perot interferometer or Mach-Zehnder interferometer. Scanning of the interrogation location may be performed through optical scanning, such a galvo-mirror, MEMS mirror, resonant scanner, polygon scanner etc., or through mechanical scanning of either the optics or the sample using single- or multiple-axes linear, or rotational stages. Extraction of relevant signal data may be performed in a solely programmatic implementation, to a relevant circuit-based processor, or through some combination of the two.

[0017] The CEPARS may be implemented using a single reference-path where phase variation is contained within a polarization state (such as circular· polarization) or may require that multiple acquisitions be performed, or may’ be implemented using multiple reference-paths which inherently provide phase variation through the use of different path lengths. Detection of the various combined beams may be performed by some manner of optical intensity detector such as a photodiode, balanced photodiode, avalanche photodiode, etc., CCD, EMCCD, iCCD, CMOS, etc., or an array of aforementioned detectors.

[0018] The SDCG-PARS interrogation may be implemented using either a pulsed source or a CW source which is modulated when using some form of sampie-and-hokl detector array such as a CCD, EMCCD, iCCD etc., or may be implemented using a CW source when using some form of rapid optical switching such as a shutter or optical switch, or when using a higher bandwidth detector array such as a photodiode, balanced photodiode, avalanche photodiode, etc..

[0019] The CEPARS is distinct from time-domain optical coherence tomography (TD- OCT) in that it; (1) may include the use of a pulsed excitation laser, and (2) may be sensitive to optical absorption contrast. CEPARS may necessitate die use of at least two optical beams such that one acts to excite the sample and the oiher acts to detect perturbations in the sample. Tire CEPARS system may be distinct from PARS in that it may include: (1 } one or more reference paths, (2) a means of separating the in-phase (sample with no delay reference) and quadrature (sample with delayed reference) beams, and (3) a means of detecting at least two of these beams.

[0020] The SDCG-PARS may be distinct from spectral-domain optical coherence tomography (SD-GCT) and PARS in that may include; (1) the use of a pulsed excitation laser (2) the use of a pulsed interrogation laser, or a rapidly modulated continuous-wave laser, or a continuous-wave laser along with the use of a gated camera exposure to detect signals on a sufficiently short timescale such that acoustic propagation is negligible, (3) a system to subtract the depth-resolved scatterer distributions before and immediately after the excitation pulse, and that it may require (4) at least two distinct interrogation events per acquisition location such that the difference between acquisitions infers depth-resolved optical absorption distribution. SDCG-PARS may necessitate the use of at least two optical beams such that one S acts to excite the sample and the other acts to detect perturbations in the sample.

[0021] Other aspects will be apparent from the description and claims below. In other aspects, the aspects described herein may be combined together in any reasonable combination as will be recognized by those skilled in the art.

[0022] A coherence gated photoacoustic remote sensing system for imaging a subsurface structure in a sample with optical resolution may include an excitation beam source configured to gener ate an excitation beam that induces ultrasonic signals i the sample at an excitation location; an interrogation beam source configured to generate an interrogation beam incident on the sample at an interrogation location, a portion of the interrogation beam returning from die sample that is indicative of the generated ultrasonic signals, the interrogation beam being a low-coherent beam; an optical system that focuses the excitation beam onto the sample at an excitation location, and focuses the interrogation beam onto the sample at an interrogation location, at least the interrogation location being below the surface of and within the sample; and a low coherence interferometer that isolates a returning portion of the interrogation beam that corresponds to an interrogation event of the sample.

[0023] The system may include a reference beam source configured to generate a reference beam that travels along a reference path, and wherein the low coherence interferometer isolates the returning portion using the reference beam. The reference beam source is configured to generate one or more additional reference beams that are phase shifted relative to the reference beam, and wherein the low coherence interferometer isolates the returning portion using the reference beam and the one or more additional reference beams. One or more additional reference beams are phased shifted by at least one of a different path length, one or more wave plates, and one or more circulators. The one or more additional reference beams are detected either in parallel or serially with the reference beam. The excitation beam and the interrogation beam are pulsed or intensity-modulated. The excitation location and the interrogation location are each below the surface of and within the sample. At least one of the excitation location find the interrogation location are within 1 ram of the surface of the sample. At least one of the excitation location and the interrogation location are greater than 1 pm below the surface of the sample. The excitation location and the interrogation location are focal points that are at least partially overlapping. The system includes a processor that calculates an image of the sample based on the returning portion of the interrogation beam. Tlie interrogation team has pulses that are sufficiently short that acoustic propagation is negligible. For each detection location, the system applies an excitation beam with more than one frequency, bandwidth, phase shift, or combination thereof. The optical system interrogates each interrogation location the sample in a non-excited slate and after an excitation beam excites die sample. The excitation beam source is configured to generate one or more excitation beams that excites the sample with a plurality of frequencies, a plurality of bandwidths or combinations thereof.

[0024] A method of using the system may include functional imaging during brain surgery; assessing internal bleeding and cauterization verification; imaging perfusion sufficiency of organs and organ transplants; imaging angiogenesis around islet transplants; imaging of skin- grafts; imaging of tissue scaffolds and biomaterials to evaluate vascularization and ' or immune rejection: imaging to aid microsurgery; or procedures for guidance to avoid cutting critical Wood vessels and nerves. A method of using the system of claim may be combined with fiuorescenc microscopy, two-photon and confocal fluorescence microscopy, Cohereni- Anti-Raman- Stokes microscopy, Raman microscopy, or Optical coherence tomography.

[0025] Tiie method may include performing microcirculation imaging or performing bloo oxygenation parameter imaging with the system

[0026] An endoscope may include the system.

[0027] A surgical microscope may include the system.

[0028] A method of remote sensing a sample may comprise the steps of: providing a coherence gated photoacoustic remote sensing system comprising an excitation beam and an interrogation beam, the interrogation beam being a low-coherent beam; causing the excitation beam to induce ultrasonic signals in the sample at an excitation location; causing the interrogation team to interrogate the sample at. an interrogation location, wherein a portion of the interrogation beam returns from the sample that is indicative of the generated ultrasonic signals the interrogation location being below the surface of an within the sample; using a low coherence interferometer to isolate a returning portion of the interrogation beam to achieve an interrogation even t of the sample.

[0029] Tiie method further comprises providing a reference beam that travels along a reference path, and wherein the low coherence interferometer isolates the returning portion using the reference beam. The method further comprises the ste of providing one or more additional reference beams that are phase shifted relative to the reference beam, and wherein the low coherence interferometer isolates the returning portion using the reference beam and the one or more additional reference beams. Tire one or more additional reference teams are phased shifted by at least one of a different path length, one or more wave plates, and one or more circulators. The one or more additional reference beams are detected either in parallel or serially with the reference beam The excitatio beam and the interrogation team are pulsed or intensity-modulated. The excitation location and the interrogation location are each below the surface of and within the sample. At least one of the excitation location and the interrogation location are within 1 mm of the surface of the sample. At least one of the excitation location and the interrogation location are greater than 1 p below the surface of the sample. The excitation location and the interrogation location are focal points that are at least partially overlapping. The method further comprises the step of calculating an image of the sample based on the returning portion of the interrogation beam. The interrogation beam has pulses that are sufficiently short that acoustic propagation is negligible. For each detection location, the excitation beam is operated to provide with more than one frequency, bandwidth, phase shift, or combination thereof. The method further comrpises the step of interrogating each interrogation location in a non-excited state and after the excitation beam excites the sample. The excitation beam comprises one or more excitation beams that excites the sample with a plurality of frequencies, a plurality of bandwidths or combinations thereof.

BRIEF DESCRIPTION OF THE DRAWINGS

[0030] These and other features will become more apparent from the following description in which reference is made to the appended drawings, the drawings are for tire purpose of illustration only and are not intended to be in any way limiting, wherein:

[0031] In this patent document the word "comprising" is used in its non-limiting sense to mean that items following the word are included, but items not specifically mentioned are not excluded A reference to an element by the indefinite article "a" does not requires that there be one and only one of the elements.

[0032] The scope of the following claims should not he limited by the preferred embodiments set forth in the examples above and in the drawings, but should be given the broadest interpretation consistent with the description as a whole.

[0033] Fig. 1 depicts a schematic overview of the excitation pathway.

[0034] Fig. 2 depicts a schematic overvie ' of the interrogation pathway.

[0035] Fig. 3 depicts a schematic view implementation of optical sources.

[0036] Fig. 4 depicts a schematic view of yet another implementation of optical sources.

[0037] Fig. 5 depicts a schematic view of an implementation of a beam combiner.

[0038] Fig. 6 depicts a schematic view of yet another implementation of a beam combiner.

[0039] Fig. 7 is a graphical illustration of the PARS mechanism.

[0040] Fig. 8 is a graphical illustration of CEPARS signals.

[0041] Fig. 9 depicts an imaging process flow diagram for CEPARS.

[0042] Fig. 10 is a schematic view of an example system layout for a CEPARS (Parallel).

[0043] Fig. 1 J is a schematic view of another example system layout for a CEPARS

(Parallel).

[0044] Fig. 12 is a schematic view of yet another example system layout for a CEPARS (Serial).

[0045] Fig. 13 is a schematic view of an endoscopic example system layout for a CEPARS.

[0046] Fig. 14 is a graphical illustration of an outline of the SDCG-PARS detection mechanism, primarily highlighting the relative time with which key process is ar e carried out.

[0047] Eig. 15 is a graphical illustration and enlargement of an example of a SDCG- PARS spectrum both before and after photoacoustic excitation.

[0048] Fig. 16 is an imaging process flow diagram for SDCG-PARS

[0049] Fig. 17 is a schematic view of an example system layout for a SDCG-PARS.

[0050] Erg. 18 is a schematic view of an example endoscopic system layout for a SDCG-

PARS

[0051] Fig. 19 is a schematic view of a system Layout for a CEPARS with a quadrature interferometer.

DETAILED DESCRIPTION

[0052] Fig. 1 shows a high-level overview of the excitation path. This primarily consists of an optical excitation source (1), an optical scanning system (2), and focusing optics (3) such as an objective lens which focuses the light onto the sample (4). The purpose of the excitation path is to direct the excitation source onto the sample to produce photoacoustic excitation within the sample.

[0053] Fig. 2 shows a high-level overview of the interrogation path. In general this consists of an optical interrogation source (5). an optical combiner (6), an optical reference path (7), an optical detector (8), and is directed onto the same optical scanning system (2), focusing optics (3), and sample (4) as in Fig. 1. The primary purpose of the interrogation path is to direct a portion of the interrogation somce onto the sample, another portion being directed to the reference path to provide a desired reference path length, then to combine the beam from the sample path and the reference path as to perform low-coherence interferometry at the beam combiner. These combined optical signals are then processed appropriately at the detector to extract desired information.

[0054] Fig. 3 shows one possible implementation of the (1) excitation source, or the (5) interrogation source which consists of one or more pulsed or modulated optical radiation sources (101) of one or more optical wavelength (1,2,...,N) which are fiber coupled (102) together at their respective outputs. The optical fibers may be of any type such as multimode, single mode, polarization maintaining, nonlinear etc.

[0055] Fig. 4 shows another possible implementation of the ( 1) excitation source, or the (5) interrogation source which consists of one or more pulsed or modulated optical radiation sources (101) of one or more optical wavelengths (1 ,2, ... ) which are coupled together through tree-space optics (103) such as beam combiners or diehroic minors.

[0056] Fig. 5 shows one possible implementation of the (6) beam combiner which consists of a (11) fiber-based device such as a fiber-based interferometer or fiber-based coupler.

[0057] Fig. 6 shows another possible implementation of tire (6) beam combiner which consists of a (10) free-space optical beam combiner in a Michelson interferometer layout. Note that other free-space interferometer layouts may be used such as common path interferometer (using specially designed interferometer objective lenses), Fizeau interferometer, Ramsey interferometer, Fabry-Perot interferometer and Mach-Zehnder interferometer.

[0058] Fig. 7 highlights aspects of the PARS mechanism. Upon absorption of a sufficiently short excitation pulse (such that thermal and stress confinement conditions are met, typically shorter than 100ns) rapid heating will occur proportional to the local optical absorption at the excitation wavelength. This heating will in turn produce significant pressures known as photoacoustic initial pressures through thermo-elastic expansion following p 0 a; h A Gfm b where h L is a conversion efficiency factor, G is a material property known as the Griineissen parameter, f is the focal fiiience of the excitation beam and p a is the optical absorption of the medium at the given excitation wavelength. These pressures can be substantial easily surpassing !OOMPa for excitation pubes within the ANSI optical exposure limits. These nontrivial pressures can produce modulations 6n in the local refractive-index n 0 through the elasto-optic effect following n * = ¾ 4- dh = i½ + (fciioP 0 ) /(2pv ). where n * is the new refractive-index profile, t is the elasto-optic coefficient, p is the mass density and v s Is the acoustic propagation velocity" in certain PARS imaging embodiments, these refractive-index modulations can be measured using a continuous-wave interrogation beam which is co-focused to the excitation location. This is detected as a total intensity measured on a photodiode such that all phase information from the interrogation spectrum is rejected. Tins can be simply represented as a change in reflectivity AR from the interrogation locatio which will be the difference between the perturbed reflectivity' and the unperturbed reflectivity" R =

between the two media and i½ such that for small perturbations dh (which are themselves proportional to the optical absorption p a ) we get the approximate relationship AR cc Sntyty— n 2 (Haji Reza et al., Light: Science & Applications volume 6, page 16278 (2017), the entirety" of which is incorporated by reference herein). One interpretation of this result is that the intensity" reflectivity from an excited interface relates directly to the inherent scattering contrast (¾— n 2 ) and to the optical absorption.

[0059] In CEPARS, it may be desirable to exclude signals which have originated far from the focus. Previously, with conventional PARS embodiments, the axial characteristics were solely provided by the optical section defined by the focusing optics. However, it was found experimentally that axial performance can easily be far worse than this value. To improve tins, CEPARS may add low-coherence interferometry such that signals which have originated from a path length significantly longer or shorter (defined by the coherence length of the interrogation source) than the reference path length will be excluded. In other words, signals which have originated from a path length that is more than a threshold amount different than the reference path length may be excluded. However, tins will lead to ambiguity" within the received signal as the two paths may still provide a signal which has undergone some amount of deconstruetive interference. To combat this, CEPARS captures several (at least two) low- coherence interferometry signals which involve different reference path lengths. One example would be to compare half of the sample signal with one reference path, and the other half of the sample signal with a reference path where the phase has been offset by p /2 . For complete characterization of the received signal at least four components with appropriate phase offset such as of 0, tt/2, p, and 3p/2 are required following from quadrature interferomehy. This would allow for extraction of both an in-phase and quadrature signal simultaneously by rejecting uadesired self-interference effects and reference-path signals such that phase-derived ambiguity can be eliminated.

[0060] Fig. 8 (CEPARS Signals) shows an example of the above described signals. If we assume a single optical seatterer is placed at some location in the sample path (E s (t,v)), and the mean reference path length is scanned about that same distance then the following two signals are acquired: for interference between reference path I (E r (t, v)) (no added delay) the corresponding measured intensity signal would be processed to provide

E r dv— I c likewise, for interference between reference path 2 v )) (p /2 added delay) the corresponding measured intensity signal would be processed to provide I q (t) ¾ j jE s + E rZ l z dv— I c where I c is a calibration intensity, v is the optical frequency and E s , E r , E r2 are considered to have wide spectral content. Note that this is approximate as it assumes small self-interference effects within the sample. These signals then undergo high- pass filtering to remove the remaining mean signal offset, they are rectified, then finally their squares are summed producing the final time-domain signal Sig(t). These steps are highlighted in Fig. 9. To capture Sig(t) without approximation or calibration, a full- quadrature detection may be implemented providing a measure for example of w hich corresponds to sample-reference path delays of 0, p/Z, p. and 3p/2 respectively. From here the complete low-coherence optical quadrature can be determined as Sig{t) = ( — .) 2 + This process allows for acquisition of the low- coherence information within a rapid time scale. This is in contrast to other low-coherence methods such as time-domain optical-coherence tomography (TD-OCT) winch commonly may necessitate axial scanning across the seatterer to properly characterize the sample. Such a TD-OCT approach would be ineffective for capturing the PARS mechanism lar gely due to issues involving acquisition time and repeatability concerns.

[0061] Fig. 10 highlights one possible implementation of CEPARS. A polarized interrogation source (1001) is fed into a beam-splitter (1008) which directs a portion of the beam towards the sample path and mother portion towards the reference mirror (1005). The sample path of the interrogation is then combined with the excitation path using an appropriate dicbroic minor (1009). The two beams are then directed onto the sample (1022) using a set of scanning minors (1019) and an objective lens (1020). Here, scanning can also be performed «sing a mechanical scanning stage (1021) to overcome field of view limitations of the objective. The reference path passes through an eighth wave plate (1006) twice providing a circular polarized state where the total path length is controlled by position of the reference minor. This circular polarized state will inherently provide the two desired reference phases. The linear polarized sample path returning from the sample is then combined with the circular reference path at the beam-splitter. Excess excitation light which lias been transmitted through tire dicbroic mirror is further rejected by the use of a narrow filter (1010). Finally the two polarization states are split using a polarized beam-splitter (1013) and individual detection is then performed. Since this device will inherently be sensitive to polarization-dependent scattering in the sample, it may also be necessary to first characterize the given interrogation location with the reference path blocked such that the relative received values can be appropriately adjusted.

[0062] Fig. 1 1 highlights another possible implementation of CEPARS. This implementation features primarily fiber-based optics and takes advantage of a randomly polarized interrogation source to avoid polarization-dependent sensitivity at the sample. Tire interrogations source (1101) is spit (1110) between tire reference and sample pathways as before. Here, the reference path is further split (1 114) to provide the desired added phase offsets. Polarization-independent circulators (1 113, 1115, 1116) then redirect the reference paths (Rl, R2) towards respective beam couplers (1106, 1107) where they are combined with the sample path components (S 1 , S2).

[0063] Fig. 12 highlights another possible implementation of CEPARS. This implementation features serial acquisition as opposed to those represente and Figs. 10 an 11, which utilize a parallel capture. Serial CEPARS only may necessitate a single low- coherence interferometer, but may require multiple acquisitions. Moreover, subsequent acquisitions must be performed with a varied reference path-length. For example, a dual acquisition might take one acquisition with a n/2 phase offset relative to a first acquisition provided by 7 a piezo-mounted mirror (1205). ha this manner, the in-phase and quadrature data can still be captured.

[0064] Fig.13 highlights yet another possible implementation of CEPARS. This implementation features a serial acquisition as that presented in Fig. 12. However, rather than focusing directly onto a sample with ffee-space optics, here of the excitation and sample-path of the interrogation beams are coupled into a fiber wliich is fed through an endoscopic probe. At the distal end. optical focusing is provided by a GRIN lens (1327), and optical scanning is provided by a set of MEMS mirrors (1319). This represents a compact implementation capable of accessing challenging locales.

[0065] Fig. 19 highlights yet another possible implementation of CEPARS. This implementation utilizes a full optical-quadrature detection path. Unlike other and simpler described architectures, this embodiment may not require additional calibration, may not require assumption of small self-interference terms, and may not require multiple acquisition events providing more complete characterization of the tissue. The detection pathways include an interrogation source (1901), which is polarized ( 1905) and split (1903). The sample path transmits through a polarization-sensitive splitter (1 23), is circularly polarized (by a quarter waveplate 1925), combined with the excitation pathway (at dichroic mirror 1926) and directed to the sample. The back-reflected portion is converted back to a linear polarization state (at quarter waveplate 1925), has remaining excitation removed by a filter (1924) and the light is again passed through a linear polarizer (1922) to ensure a clean polarization state. The reference path consists of a similar non-reciprocal pathways using a quarter-wave plate (1910) and PBS (1911). A dispersion cell (1909) can be added to compensate for sample-path dispersion. The length of this path can be controlled by changing the reference mirror (1908) position for .appropriate depth selection within the sample. Tins light is circularly polarized (by quarter waveplate 1912) contributing a tt/2 phase shift along one axis and recombined with the sample pathway in a non-polarizing splitter (1917). These two paths winch are composed to multiple polarization states are further separated in two PBSs (1916, 1921) yielding the desired combinations of sample-path and reference-path phase delays for full-quadrature detection across four sensors (1913, 1915, 1918, 1920). The sensors 1913, 1915, 1918, and 1920 may be optica! sensors, such as, e.g., a single photodiode, array of photodiodes, CCD etc. Then the collected data will be processes to extract the PARS modulation quadrature information.

[0066] In some embodiments of SDCG-PARS, one goal is to provide a full depth- resolved optical absorption profile of a sample without necessitating axial optical scanning. Conceptually, this is similar to how SD-OCT is operated. However the techniques are highly distinct from each other. First it. is assumed that the optical section can be considered a collection of ideal reflectors at some spatial distribution (along the z direction) such that is can Ire represented as r E (z). By interrogating the sample with a range of optical frequencies. commonly implemented as either a swept source or a stationary broadband spectrum source, a respective reflection spectrum can be collected. This involves combining the back-reflected light from the sample with a reference such that the interference fringes now encode the locations of the optical scatterers within the optical section. Recovery of the spatial reflection distribution then simply involves performing a frequency transform on the collected spectrum. Since tire PARS mechanism involves modulation of the optical scattering properties within a sample where these modulations correspond to locations of optical absorption, by comparing the distribution of scatterers both before and directly after photoaeoustie excitation by an execution pulse, the difference at a given location will now correspond to the PARS-modulated regions winch are optically absorbing. However, although high-bandwidth detectors are ideal for such a task they may prove highly impractical for implementation, and as such there is a requirement for a means of providing these two distinct interrogations. One proposed method is the use of a short (<100ns) pulsed, or modulated interrogation laser which can effectively force a short interrogation time on a lower bandwidth detector such as a CCD array by reducing the amount of time back-reflected light from the sample will be incident on the array. This method allows for proper control over the relative timing of the excitation and interrogation pulses and the duration of the interrogation time.

[0067] Fig. 14 shows an example of the relative timing between the reflectivity ' properties of a given wavelength in the sample, and the excitation and interrogation pulses. The second interrogation pulse which corresponds to be excited sample, must be timed such that it takes frill advantage of the perturbed sample. This exact timing will vary significantly given all the available parameters such as the sample under consideration, the time evolution char acteristics of die excitation, and the time revolution characteristics of the interrogation. In general, the rising edge of the interrogation will be less than Ips from the rising edge of the excitation. Tire duration of the intenogation pulse will also be less than Ips.

[0068] Fig. 15 shows an example of two collected spectra. One of the spectra is associated with an unperturbed interrogation event, the other is associated with an excited interrogation event. The small differences Dh between the spectra are associated with the PARS-modulated regions.

[0069] Fig. 16 shows a flow chart of the collection and processing involved with SDCG- PARS. The two collected spectra are first deconvolved with the original spectral content S(v). Here, other processing steps may be taken to reduce the effects of noise and other non- desirable effects. The spectra are then transformed back into a physical distribution representing the relative strength of optical scattering at a given depth r s (z). Tire envelope of each scattering distribution is taken, then the two envelopes are subtracted from each other to form a SDCG-PARS A-line. One of fee two original envelopes can also be used to produce a conventional SD-OCT A-line.

[0070] Fig. 17 shows an example system of a fiber-based SDCG-PARS. A pulsed interrogation source (1701) is split (by a splitter 1703) such that a portion is collected at a detector (1704) to characterize pulse-to-puise consistency. The oilier portion is split into a sample path and a reference path. The reference path is directe on a reference minor (1711) such that the total pat length will be appropriately similar to the total sample path length facilitating low-coherence interferometry. The sample path is combined with the pulsed excitation source through a multiplexer (1713). The two beams are then scanned along the surface of the sample with a set of galvanometer mirrors (1725) and an appropriate objective lens (1726). The hack-reflected light from the reference path and fire sample path are then combined in a fiber coupler ( 1706) such that they interfere with each other. This resulting light is then fed into a CCD-based spectr ometer (1705) for detecting of the spectra.

[0071] Fig. 18 shows another example of a SDCG-PARS, here with an endoscopic implementation lire primary difference between this and the previous embodiment is that after the multiplexer (1813), the combined beams are fed into an endoscopic casing (1812). Positioning of the final focus is now controlled by the use of a GRIN lens (1815) at the distal end of the fiber which is focusing through a MEMS mirror (1816) which provides lateral scanning of the interrogation spot on fee sample (1817).

[0072] It will be apparent that other examples may be designed wife different components to achieve similar results. Other alternatives may include various coherence length sources, use of balanced photodetectors, interrogation-beam modulation, incorporation of optical amplifiers in the return signal path, etc.

[0073] During in vivo imaging experiments, no agent or ultrasound coupling medium are required. However fee target can be prepared wife water or any liquid such as oil before non- contact imaging session. No special holder or immobilization is required to hold the target d ing imaging sessions.

[0074] Other advantages that are inherent to the structure will be apparent to those skilled in fee art. The embodiments described herein are illustrative and not intended to limit the scope of the claims, which are to be interpreted in light of the specification as a whole.

[0075] The excitation beam may be any pulsed or modulated source of electromagnetic radiation including lasers or other optical sources. In one example, a nanosecond-pulsed laser was used. The excitation beam may be set to any wavelength suitable for taking advantage of optical (or other electromagnetic) absorption of the sample. The source may· be monochromatic or polychromatic.

[0076] The interrogation beam may be any pulsed, or modulated source of electromagnetic radiation including lasers or other optical sources. Any wavelength can be used for interrogation purpose depending on the application.

[0077] CG-PARS may use any interferometry designs such as common path interferometer (using specially designed interferometer objective lenses), Miche!son interferometer, Fizeau interferometer, Ramsey interferometer, Fabry-Perot interferometer, Mach-Zehnder interferometer, and optical-quadrature detection. The basic principle is that phase (and maybe amplitude) oscillations in the probing receiver beam can be detected using interferometry and detected at AC, RF or ultrasonic frequencies using various detectors.

[0078] In one example, both excitation and receiver beams may be combined and scanned. In this way, photoacoustic excitations may Ire sensed in the same area as they·' are generated and where they are the largest. Other arrangements may also be used, including keeping the receiver beam fixed while scanning the excitation beam or vice versa. Galvanometers, MEMS mirrors, polygon scanners, and stepper/DC motors may be used as a means of scanning the excitation beam, probe/reeeiver beam or both.

[0079] The excitation beam and sensing/ receiver beam can be combined using dichroic minors, prisms, beam splitters, polarizing beam sphiters etc. They can also be focused using different optical paths

[0080] The reflected light may be collected by photodiodes, avalanche photodiodes, phototubes, photomultipliers, CMOS cameras, CCD cameras (including EM-CCD, intensified-CCDs, back-thinned and cooled CCDs), etc. The detected signal may be amplifie by an RF amplifier, lock-in amplifier, trans-impedance amplifier, or other amplifier configuration. Also different methods may be used in order to filter the excitation beam from the receiver beam before detection. CG-PARS may use optical amplifiers to amplify detected light.

[0081] A table top, handheld, endoscopic, surgical microscope, or ophthalmic CG-PARS system may be constructed based on principles known in the ait. CG-PARS may be used for A-. B- or C- scan images for in vivo, ex vivo or phantom studies.

[0082] CG-PARS may be optimized in order to take advantage of a multi-focus design for improving the depth-of-fbcus of 2D and 3D GR-CG-PARS imaging. The chromatic aberration in the collimating and objective lens pair may be harnessed to refocus light from a fiber into the object so that each wavelength is focused at a slightly different depth location. Using these wavelengths simultaneously may be used to improve the depth of field and signal to noise ratio (SNR) of CG-PARS images. During CG-PARS imaging, depth scanning by wavelength tuning may be performed.

[0083] The CG-PARS system may be combined with other imaging modalities such as fluorescence microscopy, two-photon and confocal fluorescence microscopy. Coherent- Anti- Raman-Stokes microscopy, Raman microscopy, Optical coherence tomography, other photoacoustic and ultrasound systems, etc. This systems could be combined by designing an optical combiner to integrate the optical paths of each systems. Also a processor to synchronise the results if necessary and analyse the results either separately or in combination. These integrated modalities can bring complementary imaging contrast. This could permit imaging of the microcirculation, blood oxygenation parameter imaging, and imaging of other molecularly-specifrc targets simultaneously, a potentially important task that is difficult to implement with only fluorescence based microscopy methods. A multi- wavelength visible laser source may also be implemented to generate photoaconstic signals for functional or structural imaging.

[0084] Polarization analyzers may be used to decompose detected light into respective polarization states. Tire light detected in each polarization state may proride information about ultrasound-tissue interaction.

[0085] APPLICATIONS

[0086] It will be understood that the system described herein may be used in various ways, such as those purposes described above, and also ma be used in other ways to take advantage of the aspects described above. A non-exhaustive list of applications is discussed below.

[0087] The system may be used for imaging angiogenesis for different pre-clinical tumor models.

[0088] Tire system may also be used for clinical imaging of micro- and macro-circulation and pigmented cells, which may find use for applications such as in (1 ) the eye, potentially augmenting or replacing fluorescein angiography: (2) imaging dermatological lesions including melanoma, basal cell carcinoma, hemangioma, psoriasis eczema, dermatitis, imaging Mohs surgery, imaging to verify tumor margin resections; (3) peripheral vascular disease; (4) diabetic and pressure ulcers; (5) bum imaging; (6) plastic surgery an microsurgery; (7) imaging of circulating tumor cells, especially melanoma cells; (8) imaging lymph node angiogenesis; (9) imaging response to photodynamic therapies including those with vascular ablative mechanisms; (10) imaging response to chemotherapeutics including anti-angiogenic drags; (11) imaging response to radiotherapy.

[0089] The system may be useful in estimating oxygen saturation using multi-wavelength photoacoustic excitation and CG-PARS detection and applications including: ( 1 ) estimating venous oxygen saturation where pulse oximetry cannot be used including estimating cerebrovenous oxygen saturation and central venous oxygen saturation. This could potentially replace catheterization procedures which can be risky, especially 7 iu small children and infants.

[0090] Oxygen flux and oxygen consumption may also be estimated by using CG-PARS imaging to estimate oxygen saturation, and an auxiliary 7 method to estimate blood flow in vessels flowing into and out of a region of tissue.

[0091] The system may also have some gastroenterological applications, such as imaging vascular beds an depth of invasion hi Barrett’s esophagus and colorectal cancels. Depth of invasion is key to prognosis and metabolic potential. Gastroenterological applications may be combined or piggy-backed off of a clinical endoscope and the miniaturized CG-PARS system may be designed either as a standalone endoscope or fit within the accessory 7 channel of a clinical endoscope.

[0092] The system may 7 have some surgical applications, such as functional imaging during brain surgery, use for assessment of internal bleeding and cauterization verification, imaging perfusion sufficiency of organs and organ transplants, imaging angiogenesis around islet transplants, imaging of skin-grafts, imaging of tissue scaffolds and biomaterials to evaluate vascularization and immune rejection imaging to aid microsurgery, guidance to avoid cutting critical blood vessels and nerves.

[0093] Other examples of applications ma include CG-PARS imaging of contrast agents in clinical or pre-clinical applications: identification of sentinel lymph nodes: non- or minimally-invasive identification of tumors in lymph nodes; imaging of genetically-encoded reporters such as tyrosinase, chromoproteins, fluorescent proteins for pre-clinical or clinical molecular imaging applications; imaging actively or passively targeted optically absorbing nanoparticles for molecular imaging; and imaging of blood clots and potentially staging the age of the clots.

[0094] In some embodiments, any suitable technology, such as, e.g., OCT, can be used for surface topology (for constant- or variable-depth focusing for photoacoustic remote sensing technologies) before imaging with CG-PARS.

[0095] In at least some embodiments, systems of the present disclosure may include variable-focal-length lenses (including voice-eoil-driven, MEMS-based, piezoelectric-based, and tunable acoustic gradient lenses). Furthermore, systems of the present disclosure may include double-clad fiber couplers for both OCT and PARS microscopy (including CG- PARS) to deliver excitation light (and/or interrogation light) from a single-mode fiber to the sample, hut collect interrogation light using the multi-mode cladding of the double-clad fiber. Systems of the present disclosure also may Ire used with angiography or Doppler.

[0096] Embodiments of the present disclosure may include one or more of the following advantages:

1. The proposed CG-PARS provides depth-dependent contrast which is directly proportional to optical absorption of the excitation laser. For example, CW CE-PARS extracts modulated components of signals using high-pass or bandpass filters. Pulsed detection systems associated with Pulsed CE-PARS or SD-CG-PARS uses differences in detected signals with and without excitation pulses.

2. The coherence length of the source is preferably shorter than the depih-of- focus of the interrogation beam into the sample, and more preferably significantly shorter. In tills way, improved depth resolution can be achieved by use of coherence-gating.

3. The proposed SD-CG-PARS system incorporates a spectrometer and is able to detect enveloped A-scans with and without excitation pulses (or with different pulse energies). The system uses a processor for extracting differences in the enveloped A- scans with and without excitation pulse (or with different pulse energies).

4. in the proposed CE-PARS system, there may be two or more interferometers, or a method for sequentially interrogating with two or more successive reference path phase shifts and a processor for combining serial or parallel acquisitions to extract temporal modulations of tlie envelope signal.

5. The proposed CG-PARS methods uses OCT signals to detect refractive index changes associated with initial pressures and uses at least two acquisitions (either in serial or parallel with multiple detectors). In SD-CG-PARS, an A-scan OCT envelope acquisition is obtained with and without and excitation pulse, where each A-scan is acquired with a spectrometer. In CE-PARS we in-phase an quadrature components of the signal are acquired.

6. As noted, the SD-CG PARS method uses a spectrometer. Additionally, SD-CG-PARS may be used to detect enveloped OCT A-seans with and without excitation pulses (or with different pulse energies). The phase in the detected signal may be removed to form an envelope. For SD-CG-PARS, a processor may be used to extract differences in the enveloped A-seans with and without excitation pulses (or with different pulse energies).

[0097] Certain examples of remote sensing systems may be described as follows:

1. A Spectral-Domain Coherence-Gated PARS Tomography (SD-CG-PARS Tomogr aphy) system having:

a. A pulsed excitation electromagnetic radiation source

b. A low-coherence interrogation light source, the coherence length being the principal determinant of the depth resolution. Typically, the interrogation wavelengths and the excitation wavelengths are spectrally distinct, but in an optional embodiment, the excitation and interrogation sources could be one and the same. c. A combiner to combine the pulsed excitation beam and the interrogation beam to enable eo- scanning of both beams

d. Focusing iens(es) for focusing the respective or combined beams and for collecting interrogation light horn the sample.

e. An interferometer, having a splitter to split the interrogation beam into a reference path and a signal path, the reference path having an adjustable path-length, and flie signal path returning collected signals back to interfere with reference path light.

f. A light analysis module consisting of a spectrometer (with various types of dispersive elements: gratings, prisms, etc) and detector arrays (CCD, CMOS, photodiode array).

g. A temporal gating system to ensure that optical signals recorded after the excitatio pulses are read out within a short (< teas of nanoseconds or <lus) time scale before acoustic waves propagate far from then origin. Specifically, the acoustic distance-of-propagation over the interrogation readout time should not be significantly greater than the desired axial or lateral spatial resolution. This temporal gating could be accomplished using (1) a (fs-us-scale pulsed interrogation source and pulse-sequencer and acquisition electronics carefully timed to read out signals within nanoseconds after the excitation source. (2) an optical or electronic shutter with nanosecond-scale response times to enable the capture of ONLY the interrogation light within the desired temporal window (3) fast photodiode array to electronically capture time domain signals from each element and capturing only the first T time samples.

h. A pulse sequencer and acquisition system for forming at least two OCT A-sean hues per scan position: one with an excitation pulse and one without an excitation puke OR one with a different excitation pulse energy than another.

i. Optional reference photodiode measurement subsystem to account for pulse-fo-pulse variations of the excitation source and for variations in the interrogation source.

j. Optional programmable controller and actuator to adjust the reference patMength between scans or to adjust the desired depth-sectioning.

fc Optional filter to reject excitation laser wavelengths from being detected by tire spectrometer detectors.

1. A processor for processing the OCT RF A-sean lines (with- and without excitation pukes or with differing pulse energies and optionally with an without reference pathiengtir shifts) to form CG-PARS A-Scans with contrast proportional to optical absorption at each depth position. One such processor embodiment includes forming the envelope of each OCT A-sean and subtracting the envelopes of A-seans with and without excitation laser pulses. This strategy has the advantage of eliminating unwanted phase-noise sensitivity· but will still capture refractive index changes associated with photoacoustic initial pressures. m. A processing system to lender and display OCT and CG-PARS images

2. A coherence-enhanced PARS (CE-PARS) microscopy system having:

a. A pulsed excitation light source

b. A low-coherence inteirogatioii laser, the coherence length being the principal determinant of the depth resolution. Typically, the interrogation wavelengths and the excitation wavelengths are spectrally distinct but in an optional embodiment, the excitation and interrogation sources could be one and tire same.

c. A combiner to combine the pulsed excitation beam and the interrogation beam to enable co-scanning of both beams

d. Focusing iens(es) for focusing the respective or combined beams and for collecting interrogation light from the sample.

e. An interferometer, having a splitter to split the interrogation beam into a reference path and a. signal path, file reference path having an adjustable path-length, and the signal path returning collected signals back to interfere with reference path light

f. Light detection module{s) including associated optional amplifiers and filters, for example, consisting of photodiode(s) or balanced photodiode(s). Filters may be included to reject DC scattered light and collect only the modulated component in tire ca.se of CW interrogation beams. See description of module for pulsed interrogation light below.

g. A method for acquiring effective inphase- and quadrature complex envelope signals from the interfering ligh using one of two methods (1) serially, by performing a point-scan, lateral-scan, depth-scan or C-scan then adjusting the reference pathiengfh by TT/2 phase then scanning again. (2) in parallel by using an additional interferometer with a reference path differing by p/2 from the reference path of the other interferometer. Tins parallel interferometer may be implemente with separate optical paths or as a common-path configuration. Note that this quadrature-sampling scheme offers the flexibility of C-scanning or en-face scanning at a particular depth gating (or depth range) without requiring acquisition of complete depth-scans (A-scans) to create an effective PARS image within a precise depth- section. If an A-scan is acquired, there must be an excitation pulse for every depth sample in file A-scan line which could lead to unwanted persistent laser exposure compared to the scanned beam case.

li A processor for estimating the envelope or specifically, the magnitude of the complex envelope signal for cases with and without an excitation pulse (or with excitation pulses of different strengths).

i. A processor for processing the OCT RF envelope signals (with- and without excitation pulses or with differing pulse energies) to form CE-PARS signals with contrast proportional to optical absorption at each scan position. One such processor embodiment includes forming the envelope of each OCT signal and subtracting the envelopes with and without excitation laser pulses. This strategy lias the advantage of eliminating unwanted phase-noise sensitivity but will still capture refractive index changes associated with photoacoustic initial pressures.

j . A temporal gating system to ensure that optical signals recorded after the excitation pulses are read out within a short (< tens of nanoseconds) time-scale before acoustic waves propagate far from their origin. Specifically, the acoustic distance-of-propagaiion over the interrogation readout time should not be significantly greater than the desired axial or lateral spatial resolution. This temporal gating could be accomplished using (1) a nanosecond-scale pulsed interrogation source and pulse-sequencer and acquisition electronics eareMly timed to read out signals within nanoseconds after the excitation source. (2) an optical or electronic shutter with nanosecond-scale response times to enable the capture of ONLY the interrogation light within the desired temporal window (3) acquiring the photodiode signal as a function of time then sampling only the first few tens to hundreds of ns OR (4) using an analog or digital peak detector to extract the peak (envelope) signal or peak-fo-peak (envelope) signal for each pulse.

k. A processing system to render and display OCT and CG-PARS images.

3. A pulsed interrogation detection subsystem which involves capturing an interrogation pulsed signal from the sample (with or without reference beam interference) both with or without an excitation pulse (or with differing pulse energies) and subtraction of the respective signals or estimating then- relative difference normalized to t e OCT signal without an excitation pulse present. This can be done by recording amplified photodiode signals with an analog-to-digiial converter and doing the subtraction (and optionally division) operations digitally. It can also be done with analog electronics

4. A functional imaging system involving sequential pulses using (1) different excitation wavelengths or (2) different pulse widths (e.g. ps pulses and ns pulses). In the ease of both (1) and (2) the PARS initial pressure signals are proportional to optical absorption and detected optically using interrogation beams using the CG-PARS systems described above, or using previously described interferometric or non- interferometric PARS systems.

[0098] Examples of methods of for remote sensing may be described as follows:

a. (SDCG-PARS) A method for interrogating the optical properties of a sample which comprises :

A method of generating photoacoustic signals within a sample;

A low-coherence interferometer used to detect photoacoustic signals;

A method for directing light towards a sample at a given location;

A method of collecting light from a sample at a given location:

An optical spectrum detector;

A processor for collecting multiple optical spectra; and

A processor for extracting differences between multiple optical spectra.

I. The method of statement a., wherein the method of generating photoacoustic signals within a sample comprises a narrowband or broadband electromagnetic source which is one of a pulsed source, or a continuous-wave source which is intensify modulated.

i. The method of I., wherein the portion of tire excitation source is detected by photodiode to account for pulse-to-pulse variations.

II. The method of statement a., wherein the low-coherence interferometer comprises a broadband electromagnetic source which is one of a pulsed source, or a continuous-wave source which is intensify modulated, a method of splitting this beam into a reference path and a sample path, and a method of combining the beams returning from the reference path and the sample path.

i. The method of statement II, wherein the optical spectrum detector comprises one or more dispersive elements (gratings, prisms, etc) and one or more detector arrays (CCD, CMOS, photodiode, etc.).

ii. The method of statement II., wherein the portion of the interrogation source is detected by photodiode to account for pulse-to-pulse variations.

III. The method of statement a, wherein the low-coherence interferometer comprises a broadband continuous-wave source electromagnetic source, a method of splitting this beam into a reference path and a sample path, and a method of combining the beams returning from the reference path and die sample path.

i. The method of statement HI, wherein the optica! spectrum detector comprises one or more dispersive elements (gratings, prisms, etc) and one or more high-bandwidth detector arrays (photodiode, avalanche photodiode, etc.).

ii. The method of statement HI, wherein the portion of the interrogation source is detected by photodiode to account for power and variations.

IV. The method of statement a. , wherein the method for directing towards and from a sample comprises of an optical scanner (one or more of Galvanometer mirrors, resonant minors, MEMS mirrors, polygon scanner, etc.), focusing optic subsystem (objective lens, reflective objective, parabolic mirror, GRIN lens, and a system of optical filters to reject excitation wavelengths along the detection path.

V The method of statement a., wherein the method for directing towards and from a sample comprises of an light guide (optical fiber, double clad fiber, optical fiber bundle, etc.), an optical scanner (one or more of Galvanometer mirrors, resonant minors, MEMS mirrors, polygon scanner, etc.), and focusing optic subsystem (objective lens, reflective objective, parabolic minor, GRIN lens) and a system of optical filters to reject excitation wavelengths along the detection path.

VI. The method of statement a, wherein the processors for collecting multiple optical spectra and for extracting differences between multiple optical spectra are implemented as electronic devices.

b. (Parallel CEPARS) A method for interrogating the optical properties of a sample which comprises:

A method of generating photoacoustic signals within a sample:

Two or more optical low-coherence interferometers used to detect photoacoustic signals;

A method for directing light towards a sample at a given location;

A method of collecting light from a sample at a given location;

A processor to combine data channels from the interferometers; and

A processor to extract temporal modulations of the envelope signal.

I. The method of statement h, wherein the method of generating photoacoustic signals within a sample comprises a narrowband or broadband electromagnetic source which is one of a pulsed source, or a continuous-wave source which is intensity modulated.

i. The method of statement I, wherein the portion of the excitation source is detected by photodiode to account for pulse-to-pulse variations.

P. The method of statement b, wherein the method for directing towards and from a sample comprises of an optical scanner (one or more of Galvanometer mirrors, resonant minors, MEMS mirrors, polygon scanner, etc.), focusing optic subsystem (objective lens, reflective objective, parabolic mirror, GRIN lens, and a system of optical filters to reject excitation wavelengths along the detection path.

I. The method of statement b, wherein the method for directing towards and from a sample compulses of an light guide (optical fiber, double clad fiber, optical fiber bundle, etc.), an optical scanner (one or more of Galvanometer mirrors, resonant minors, MEMS mirrors, polygon scanner, etc.), and focusing optic subsystem (objective lens, reflective objective, parabolic minor, GRIN lens), and a system of optical filters to reject excitation wavelengths along the detection path. c. (serial QSCG-PARS) A method tor interrogating the optical properties of a sample which comprises:

A method of generating photoacoustic signals within a sample;

A low-coherence interferometer used to detect photoacoustic signals where reference phase must he adjusted between sequential acquisitions;

A method for directing light towards a sample at a given location;

A method of collecting light h orn a sample at a given location;

A method of acquisition which necessitates multiple acquisitions;

A processor to combine serial data channels from the interferometer; and

A processor to extract temporal modulations of the envelope signal.

I. The method of statement c, wherein the method of generating photoacoustic signals within a sample comprises a narrowband or broadband electromagnetic source which is one of a pulsed source, or a continuous-wave source which is intensity modulated ii. The method of statement I, wherein the portion of the excitation source is detected by photodiode to account for puLse-io-pulse variations .

II. Tire method of statement c, wherein the method for directing towards and from a sample comprises of an optical scanner (one or more of Galvanometer mirrors, resonant mirrors, MEMS minors, polygon scanner, etc.), focusing optic subsystem (objective lens, reflective objective, parabolic mirror, GRIN lens and a system of optical filters to reject excitation wavelengths along the detection path.

HI. The method of statement c, wherein the method for directing towards and from a sample comprises of an light guide (optical fiber, double clad fiber, optical fiber bundle, etc,), an optical scanner (one or more of Galvanometer minors, resonant mirrors, MEMS mirrors, polygon seamier, etc.), and focusing optic subsystem (objective lens, reflective objective, parabolic mirror, GRIN lens), and a system of optical filters to reject excitation wavelengths along the detection path.

d. (Quadrature CEPARS) A method for interrogating the optical properties of a sample which comprises:

A method of generating photoacoustic signals within a sample; An optical quadrature detector;

A method for directing light towards a sample at a given location;

A method of collecting light from a sample at a given location;

A processor to combine dat channels from the quadrature detector; and

A processor to extract temporal modulations of the envelope signal.

I. The method of statement cL wherein the method of generating photoacoustic signals within a sample comprises a narrowband or broadband electromagnetic source which is one of a pulsed source or a continuous-wave source which is intensity modulated i. The method of statement I, wherein the portion of the excitation source is detected by photodiode to account for pulse-to-pnlse variations.

P. The method of statement d., wherein the method for directing towards and fr om a sample comprises of an optical scanner (one or more of Galvanometer mirrors, resonant minors, MEMS mirrors, polygon scanner, etc,), focusing optic subsystem (objective lens, reflective objective, parabolic mirror, GRIN lens, and a system of optical filters to reject excitation wavelengths along the detection path.

PI. The method of statement d, wherein the method for directing towards and from a sample comprises of an light guide (optical fiber, double clad fiber, optical fiber bundle, etc.), an optical scanner (one or more of Galvanometer mirrors, resonant mirrors, MEMS minors, polygon scanner, etc.), and focusing optic subsystem (objective lens, reflective objective, parabolic mirror, GRIN lens), and a system of optical filters to reject excitation wavelengths along the detection path.