PORTABLE DEVICE FOR QUANTITATIVE MEASUREMENT OF TISSUE AUTOREGULATION AND NEUROVASCULAR COUPLING USING EEG, METABOLISM, AND BLOOD FLOW DIAGNOSTICS

Abstract

The present invention is directed to a portable device for quantitative measurement of tissue autoregulation and neurovascular coupling via portable measurement of blood flow, oxygenation, metabolism, and/or EEG signals and methods for using said device. The device may comprise a body and a plurality of legs pivotably attached to the body. The plurality of legs may comprise at least one reference electrode leg and at least one measurement electrode leg for electrical measurement, and an optical detection fiber leg and at least one optical source fiber leg for optical blood flow, oxygenation, and metabolism measurement. The present invention is additionally directed to a portable device for blood flow measurement and therapeutic photobiomodulation. The device may comprise a body and a plurality of legs. The plurality of legs may comprise at least one optical detection fiber leg and at least one optical source fiber leg, and at least one leg for therapeutic photobiomodulation.

Claims

1. A system for quantitative intracranial measurement of cerebral blood flow, oxygenation, metabolism, and autoregulation, the system comprising: a. a device body; b. one or more light sources, extending from the device body and configured to be positioned in proximity to a head of a subject, wherein one or more of the light sources are configured to emit a coherent light signal; c. one or more detectors, extending from the device body and configured to be positioned in proximity to the head and to detect one or more backscattered signals; d. a microprocessor, operatively connected to the one or more light sources and to the one or more detectors; and e. a memory component, operatively connected to the microprocessor, wherein the microprocessor is capable of executing instructions held in the memory component, the memory component comprising instructions for decoupling components of the one or more backscattered signals by: i. differentiating between components of the one or more backscattered signals that are due to different layers of the head; ii. determining a dynamic perfusion metric using the one or more backscattered signals; iii. determining a tissue absorption coefficient using the one or more backscattered signals; iv. determining a tissue reduced scattering coefficient using the one or more backscattered signals; v. calculating a value of an absolute perfusion metric, using the dynamic perfusion metric, the tissue absorption coefficient, and the tissue reduced scattering coefficient; vi. calculating a value of an absolute metabolic metric, using the absolute perfusion metric, the tissue absorption coefficient, and the tissue reduced scattering coefficient; and vii. calculating a quantitative value of cerebral autoregulation, using the absolute values of the perfusion metric and the metabolic metric; thereby providing for quantitative intracranial measurement of cerebral blood flow, oxygenation, metabolism, and autoregulation.

2. The system of claim 1, additionally comprising an electroencephalography (EEG) electrode extending from the device body, the electrode configured to allow for co-localized EEG monitoring, wherein the system is configured to detect an EEG signal and thereby allows for the evaluation of neurovascular coupling.

3. The system of claim 1, wherein the system continuously calculates the quantitative value of the autoregulation metric in real time.

4. The system of claim 1, wherein the system is non-invasive, detects an intrinsic optical signal, and does not require any exogenous analyte or contrast agent.

5. A device for quantitative intracranial measurement of a brain metric, the device comprising: a. a device body; b. one or more light sources, extending from the device body and configured to be positioned in proximity to a head of a subject, wherein the one or more light sources are configured to emit one or more light signals for measurement of a dynamic perfusion metric, a tissue absorption metric and a tissue scattering metric; c. one or more detectors configured to detect one or more backscattered light signals, the detectors extending from the device body and configured to be positioned in proximity to the head; wherein the backscattered light signals allow for determination of an absolute value of the brain metric using the dynamic perfusion metric, the tissue absorption metric and the tissue scattering metric.

6. The device of claim 5, wherein the device comprises two or more detectors with different source-detector separations, and wherein the different source-detector separations enable decoupling of signals from the skull and signals from the brain.

7. The device of claim 5, wherein the brain metric is indicative of brain perfusion, oxygenation, metabolism, or of cerebral edema.

8. The device of claim 5, wherein the brain metric comprises cerebral metabolic rate of oxygen (CMRO.sub.2), cerebral blood flow (CBF), tissue concentration of deoxy-hemoglobin (ctHb), tissue concentration of oxygenated hemoglobin (ctHbO.sub.2), tissue oxygenation (StO.sub.2), tissue water content, tissue lipid content, tissue reduced scattering coefficient, tissue scattering amplitude, tissue scattering slope, tissue reflectance, or any combination thereof.

9. The device of claim 5, wherein the dynamic perfusion metric, the tissue absorption metric, or the tissue scattering metric provides information on neuronal injury, edema, sickle cell disease, depolarization, ischemia, hypoxia, metabolic injury, impaired autoregulation, or a combination thereof.

10. The device of claim 5, wherein the light signals comprise separate coherent and modulated light signals, or individual coherent modulated light signals, and wherein the coherent light signals allow for measurement of the dynamic perfusion metric and the modulated light signals allow for measurement of the tissue absorption metric and the tissue scattering metric.

11. The device of claim 5, additionally comprising an electroencephalography (EEG) electrode configured to allow for co-localized EEG monitoring.

12. The device in claim 11, wherein concurrent EEG monitoring and optical monitoring allow for measurement of a time offset between corresponding neurologic and hemodynamic activity.

13. The device of claim 5, additionally comprising a therapeutic probe configured to administer a colocalized therapy, wherein the therapy comprises colocalized transcranial light therapy, photobiomodulation, magnetic stimulation, electrical stimulation, or another therapy.

14. The device of claim 13, wherein determination of the absolute value of the brain metric, or a change in the value of the brain metric, is configured to guide the therapy in real time.

15. The device of claim 5, wherein the device is configured to use diffuse correlation spectroscopy (DCS), frequency-domain diffuse optical spectroscopy (FD-DOS), or a combination thereof, at one or more wavelengths in a range comprising the visible, near-infrared (NIR), and short-wave infrared (SWIR) regimes.

16. The device of claim 5, wherein the sources and detectors are attached to the device body via a plurality of extendable support fibers which may be retracted into the device body.

17. The device of claim 5, wherein two or more of the detectors have different source-detector separations, and wherein the source-detector separations are measured by a distance sensor, a camera, a ruler, calipers, or by a patch or apparatus which guides placement of the sources and detectors.

18. A device for quantitative subdermal measurement of a tissue metric, the device comprising: a. a device body; b. one or more light sources, extending from the device body and configured to be positioned in proximity to a body surface of a subject, wherein the one or more light sources are configured to emit a coherent light signal and a modulated light signal; and c. two or more detectors, extending from the device body and configured to be positioned in proximity to the body surface, the detectors configured to detect one or more backscattered light signals; wherein the modulated light signal allows for decoupling of a plurality of components of the backscattered light signals, the components comprising: i. a tissue absorption component; ii. a tissue scattering component; and iii. a dynamic perfusion component; wherein the decoupled components of the backscattered light signals allow for determination of a value of the tissue metric.

19. The device of claim 18, wherein the value of the tissue metric comprises an absolute value.

20. The device of claim 18, wherein the value of the tissue metric provides information on the perfusion or metabolism of an organ, is indicative of tissue autoregulation, or allows for comparative analysis of the autoregulation of two or more body parts.

Description

BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWINGS

[0022] The patent application or application file contains at least one drawing executed in color. Copies of this patent or patent application publication with color drawings will be provided by the Office upon request and payment of the necessary fee.

[0023] The features and advantages of the present invention will become apparent from a consideration of the following detailed description presented in connection with the accompanying drawings in which:

[0024] FIG. 1A shows a diagram of an embodiment of a portable EEG and blood flow measurement device.

[0025] FIG. 1B shows a diagram of an embodiment of a portable EEG and blood flow measurement device comprising an optical filter and custom correlator.

[0026] FIG. 1C shows an additional embodiment where the correlator could be either a hardware or a software correlator and the data processing can either be on-board or off-board, and the analysis of the data may be done in real-time.

[0027] FIG. 1D shows an additional embodiment where multiple detector fibers can be used to detect optical signals from the tissue. The use of multiple detector fibers allows for multiple source-detector separations (e.g., r1 and r2), which allows the device to probe different depths (z) beneath the surface of the tissue to interrogate different regions of the brain or other bodily tissue and/or better separate the signal from the brain or other bodily tissue from that of the skin, subcutaneous tissue, and/or bone. This embodiment can also be applied to all other embodiments described in this invention.

[0028] FIG. 2 shows a diagram of a leg of a portable EEG and blood flow measurement device.

[0029] FIG. 3 shows a flow chart of EEG measurement instructions of a signal processing component.

[0030] FIG. 4 shows an additional embodiment where the data analysis algorithm can be on-board or off-board and can include a seizure detection algorithm or other algorithm assisting with diagnosis or ongoing treatment (e.g. management of medically-induced coma, ischemia detection, etc.).

[0031] FIG. 5A shows a flow chart of optical blood flow measurement instructions of a signal processing component.

[0032] FIG. 5B shows an additional embodiment where the data analysis algorithm can be on-board or off-board.

[0033] FIG. 6A shows a diagram of a portable EEG, blood flow, and oximetry measurement device comprising additional legs for oximetry measurements.

[0034] FIG. 6B shows a diagram of a portable EEG, blood flow, and oximetry measurement device comprising a leg capable of serving as the source fiber for measuring both oximetry and blood flow.

[0035] FIG. 6C shows a diagram of a portable EEG, blood flow, and oximetry measurement device comprising an electrical-optical transducer and an optical-electrical transducer.

[0036] FIG. 6D shows a diagram of a portable EEG, blood flow, and oximetry measurement device comprising a plurality of movable clamps for retracting and extracting fibers.

[0037] FIG. 7A shows a flow chart of oximetry measurement instructions of a signal processing component.

[0038] FIG. 7B shows an additional embodiment where the data analysis algorithm can be on-board or off-board.

[0039] FIG. 8 shows a diagram of a portable blood flow measurement, oximetry measurement, and therapeutic biomodulation device.

[0040] FIG. 9 shows a flow chart of a therapeutic photobiomodulation method of the present invention.

[0041] FIG. 10 shows a diagram of a portable blood flow measurement, oximetry measurement, and therapeutic biomodulation device comprising a leg for transmitting therapeutic electromagnetic stimulation.

[0042] FIG. 11 shows a flow chart of a therapeutic electromagnetic stimulation method of the present invention.

[0043] FIG. 12 shows an illustration of the concomitant measurement of electrical and optical signals, where retractable fibers, multiple source-detector separations, and separation of tissue absorption and scattering at each time point enable correction of the optical data for the effects of the scalp and skull to obtain quantitatively accurate brain blood flow, hemoglobin concentration, oxygenation, metabolism, and edema index, using a small form factor pen instrument.

[0044] FIG. 13 shows an illustration of the measurement of absolute values of hemodynamic/metabolic parameters concomitant with EEG parameters. Note: In contrast to the present invention, previous products do not feature the combined flow+oxygenation+EEG with correction for varying tissue scattering to enable calculation of absolute values of CBF, ctHbO.sub.2 ctHb, StO.sub.2, and CMRO.sub.2. Additionally, this enables the measurement of flow-metabolism autoregulation metrics.

[0045] FIG. 14 shows an illustration of the measurement of tissue scattering changes (edema, cortical spreading depolarization, etc.) concomitant with hemodynamic and EEG changes.

[0046] FIG. 15 shows one embodiment that is a multimodal platform (not to scale) for Laser Speckle Imaging (LSI) and multispectral Spatial Frequency Domain Imaging (SFDI) of the rat brain. A craniectomy (6 mm4 mm area) is performed to provide direct access to the brain for optical imaging. For SFDI, light-emitting diodes (LEDs) of 655 nm, 730 nm, and 850 nm are sequentially sent into a spatial light modulator that acts as a projector to send spatially-modulated patterns of light onto the brain. A scientific CMOS camera detects the backscattered light. For LSI, an 809 nm laser illuminates the brain with coherent light, and the remitted speckle pattern is captured at 60 fps with a CCD camera.

[0047] FIGS. 16A-C show an embodiment that illustrates how venous and parenchymal regions of interest (ROIs) demonstrate different hemodynamics in response to cardiac arrest. 16A shows images of blood flow (speckle flow index; SFI) and deoxyhemoglobin concentration (ctHb), measured from the rat brain using LSI and SFDI, respectively. ROIs over the parenchyma (P) and a large vein (V) are labeled in the ctHb image. 16B shows that ctHb (black dots) increases in these ROIs during the initial minute of asphyxia. The sigmoidal fit (red line) and linear fit (blue circles) to the measured data are used to calculate the parameter in Eq. 6. 16C shows that the rate of change of tissue ctHb (dctHb/dt) during the initial minute of asphyxia is higher in the venous ROI than parenchyma ROI.

[0048] FIG. 17 shows an embodiment showing absolute CMRO.sub.2 (M O.sub.2/min) maps of a 6 mm4 mm region of the rat brain at different time points during a CA/CPR experiment. Metabolic activity increases as anesthesia is being washed out (between Baseline and Start Ischemia), followed by a sharp decrease during ischemia. Following CPR, CMRO.sub.2 recovers to anesthesia-free baseline level (3 min post-CPR), subsequently increases to values higher than baseline (5-8 min post-CPR), and then declines to values approaching anesthetized baseline level once cerebral electrical activity resumes (12 min post-CPR). Large vessels (dark blue) have been removed from the CMRO.sub.2 images to signify that the oxygen metabolism measured is occurring in the parenchyma.

[0049] FIG. 18 shows an embodiment showing how optical properties affect calculation of CMRO.sub.2. Comparison of SFI (light blue), CMRO.sub.2 calculated using SFI (dark blue), and CMRO.sub.2 calculated using v.sub.c (red), which corrects for the effects of tissue optical properties on SFI. For ease of comparison, the three curves are normalized to their value at a point near the end of the washout period (t4 min). This correction reveals differences in the observed rates of change in CMRO.sub.2 during reperfusion and resumption of ECoG bursting, suggesting the need to take optical properties into account even for rCMRO.sub.2 measurements.

[0050] FIGS. 19A-C show an embodiment demonstrating how cerebral electric activity after resuscitation from cardiac arrest is associated with baseline (pre-cardiac arrest) CMRO.sub.2. 19A shows time courses of absolute CMRO.sub.2 for two representative subjects. The subject with the higher baseline absolute CMRO.sub.2 (blue) had higher ECoG Information Quantity (IQ) 90 min post-CPR than the one with the lower baseline absolute CMRO.sub.2 (red). 19B shows a scatter plot illustrating correlation (Spearman; r=0.63; p=0.076) between baseline absolute CMRO.sub.2 and cerebral electrical activity (ECoG IQ) 90 min post-CPR. Blue dots are subjects with mild (5 min) ischemia duration, and red dots are subjects with prolonged (7 min) ischemia duration. These results suggest that the pre-injury, baseline value of CMRO.sub.2 may be predictive of cerebral electrical recovery following CA and CPR. 19C shows that subjects in the High IQ group (ECoG IQ>0.75) had significantly higher baseline absolute CMRO.sub.2 than subjects in the Low IQ group (ECoG IQ<0.75) 90 min post-CPR (p=0.016).

[0051] FIG. 20 shows an embodiment in which absolute CMRO.sub.2 values (M O.sub.2/min) measured with the combined SFDI+LSI optical system are in with the range of absolute CMRO.sub.2 values reported in previous preclinical studies using different imaging modalities (MRI/MRS, PET). Each LSI+SFDI data point represents an individual subject (n=10) imaged in this report. Each data point for MRI/MRS and PET represents the reported average CMRO.sub.2 of the subjects measured in individual studies.

[0052] FIG. 21 shows an embodiment where blood flow and oxygenation data are combined into a derivation of a different equation for absolute tissue metabolic rate of oxygen without the need to induce a zero-flow ischemic perturbation in the subject. By removing the need for the zero-flow condition, this embodiment shows significant additional clinical translational potential by enabling non-invasive measurement of tissue metabolism in physiological units without perturbing the subject. This embodiment involves multiplying the tissue deoxy-hemoglobin concentration (ctHb) by the volume of tissue being oxygenated per unit time to calculate tissue metabolic rate of oxygen. The volume of tissue oxygenated per unit time is determined by multiplying the Brownian diffusion coefficient Db (the area of tissue perfused per unit time) by the mean penetration depth of light in the tissue (determined from the measured absorption and scattering coefficients using a light transport model). The light transport model used in the penetration depth calculation was a diffusion model in this embodiment, but it could also be a Monte Carlo model, finite element model, or other type of analytical, semi-empirical, numerical, computational, or other type of model.

[0053] FIG. 22A shows an embodiment where the Brownian diffusion coefficient Db is calculated to quantify cerebral blood flow in physiological units (mm.sup.2/sec) at multiple time points throughout a cardiac arrest and resuscitation experiment in a rat. The calculation of Db is performed using a method similar to that in FIG. 13, where the measured blood flow index is corrected for tissue absorption and scattering coefficients.

[0054] FIG. 22B shows an embodiment where the mean penetration depth of unmodulated (f.sub.x=0) near-infrared light in tissue is estimated using a diffusion model with the measured tissue absorption and scattering coefficients, for use in the metabolic rate of oxygen calculation described in FIG. 21.

[0055] FIG. 22C shows an embodiment where the mean penetration depth of modulated (f.sub.x0.3/mm) near-infrared light in tissue is estimated using a diffusion model with the measured tissue absorption and scattering coefficients, for use in the metabolic rate of oxygen calculation described in FIG. 21.

[0056] FIG. 23 shows an embodiment where the cerebral metabolic rate of oxygen (CMRO2) in physiological units (micro-moles of oxygen per minute) is calculated (using the method described in FIG. 21) at multiple time points throughout a cardiac arrest and resuscitation experiment in a rat.

[0057] FIG. 24 shows an embodiment illustrating a significant correlation between the slope of the deoxy-hemoglobin concentration during entry into cardiac arrest in a rat model and the baseline CMRO2 calculated with the method described in FIG. 21.

[0058] FIG. 25 shows an embodiment illustrating that the mean baseline CMRO2 (over the last minute prior to start of asphyxia in our rat model of cardiac arrest), calculated with the method described in FIG. 21, correlates with recovery of cerebral electrical activity (ECoG Information Quantity; IQ) 90 minutes after resuscitation from cardiac arrest.

[0059] FIG. 26 shows an embodiment illustrating that the mean baseline Brownian diffusion coefficient Db of cerebral blood flow (over the last minute prior to start of asphyxia in our rat model of cardiac arrest), calculated with the method described in FIG. 22A, correlates with recovery of cerebral electrical activity (ECoG Information Quantity; IQ) 90 minutes after resuscitation from cardiac arrest.

[0060] FIG. 27 shows an embodiment where the measured blood flow index can be corrected for tissue absorption and scattering (using a method similar to that in FIG. 13) to obtain a spatially-resolved map of either the Brownian diffusion coefficient Db or the directed-flow speed vc. It is important to note that this method can also be used to map either or both of these coefficients as a function of time, and that these coefficients can be extracted simultaneously, along with weighting coefficients characterizing the relative contributions of each type of flow to the measured blood flow index. Types of flow in this model could include other flow equations besides those described here for diffuse and directed flow. FIG. 27 also shows an embodiment where the directed flow speed vc is used in a modified version of the absolute CMRO2 equation described in FIG. 21. In this modified equation, Db is replaced with vc and the mean penetration depth of the light in tissue is squared to maintain units of micro-moles of oxygen per unit time for the calculated CMRO2. As in FIG. 21, this CMRO2 equation does not require a zero-flow ischemic perturbation to the tissue. This equation can also be modified further to include a combination of diffuse flow, directed flow, and other types of flow within the same model.

DETAILED DESCRIPTION OF THE INVENTION

[0061] As used herein, the term intracranial measurement refers to the measurement of tissue (such as brain tissue) inside the skull, by a probe which may be positioned outside the skull (i.e. the measurement is taken through the skull in a non-invasive manner).

[0062] In one embodiment, the present invention features a system for quantitative intracranial measurement of cerebral blood flow, oxygenation, metabolism, autoregulation, or a combination thereof. As a non-limiting example, the system may comprise: a device body; one or more light sources; one or more detectors, a microprocessor, and a memory component. The light sources and the detectors may extend from the device body and be configured to be positioned in proximity to a head of a subject. As a non-limiting example, the light sources and the detectors may extend from the device body so as to pass through the hair of the subject and contact the skin surface at a plurality of points in a measurement area. In another embodiment, the light sources and detectors may not extend from the device body, but instead be integrated within an end of the device body. In some embodiments, the microprocessor may be operatively connected to the one or more light sources, the one or more detectors, or a combination thereof. In further embodiments, the memory component may be operatively connected to the microprocessor, and the microprocessor may be capable of executing instructions held or stored in the memory component. According to preferred embodiments, one or more of the light sources may be configured to emit a coherent light signal. In selected embodiments, the system may be configured to detect and decouple one or more backscattered signals via the detectors. In one embodiment, the memory component may comprise instructions for decoupling components of the one or more backscattered signals. As a non-limiting example, the system may be configured to: differentiate between components of the one or more backscattered signals which are due to different layers of the head; measure or determine a dynamic perfusion metric; measure or determine a tissue absorption coefficient; measure or determine a tissue reduced scattering coefficient; calculate a value of an absolute perfusion metric, using the dynamic perfusion metric, the tissue absorption coefficient, and the tissue reduced scattering coefficient; calculate a value of an absolute metabolic metric, using the absolute perfusion metric, the tissue absorption coefficient, and the tissue reduced scattering coefficient; and calculate a quantitative value of cerebral autoregulation, using the absolute values of the perfusion metric and the metabolic metric. As another non-limiting example, the instructions for decoupling components of the one or more backscattered signals may comprise differentiating between components of the one or more backscattered signals that are due to different layers of the head; determining a dynamic perfusion metric using the one or more backscattered signals; determining a tissue absorption coefficient using the one or more backscattered signals; determining a tissue reduced scattering coefficient using the one or more backscattered signals; calculating a value of an absolute perfusion metric, using the dynamic perfusion metric, the tissue absorption coefficient, and the tissue reduced scattering coefficient; calculating a value of an absolute metabolic metric, using the absolute perfusion metric, the tissue absorption coefficient, and the tissue reduced scattering coefficient; calculating a quantitative value of cerebral autoregulation, using the absolute values of the perfusion metric and the metabolic metric; or a combination thereof, thereby providing for quantitative intracranial measurement of cerebral blood flow, oxygenation, metabolism, and autoregulation As used herein, the terms tissue absorption coefficient and tissue reduced scattering coefficient may have both hemodynamic and non-hemodynamic components, and they may have tissue-based (e.g., brain, skin, muscle, heart, lung, etc.) and non-tissue-based (e.g., implantable or injectable tissue sensors such as implantable probes, nanoparticles, quantum dots, etc.) components.

[0063] In some embodiments, the instructions for decoupling the components of the one or more backscattered signals may be executed within the device. In other embodiments, the instructions for decoupling the components of the one or more backscattered signals may be executed external to the device. As a non-limiting example, any of the devices of the present invention may comprise a wired or wireless connection to an external processor, such that the data processing may be accomplished external to the device itself. As one non-limiting example, the system may include wireless transmitters and receivers to enable communication between the device and an external data processor. In some embodiments, processing of the data obtained by the device may be enhanced via a machine learning algorithm.

[0064] In one embodiment, the quantitative value of cerebral autoregulation may be calculated by dividing the quantitative value of the absolute perfusion metric by the quantitative value of the absolute metabolic metric. Alternatively, the quantitative value of cerebral autoregulation may be calculated by dividing the rate of change of an absolute or relative perfusion metric by the rate of change of an absolute or relative metabolic metric. Alternatively, the quantitative value of cerebral autoregulation may be calculated by creating a coordinate in a multi-dimensional space in which the x-axis is the absolute perfusion metric, the y-axis is the absolute metabolic metric, and other axes (dimensions) may include the rate of change in perfusion, the rate of change in metabolism, and the absolute values and/or rates of change in oxygenated hemoglobin concentration, deoxygenated hemoglobin concentration, total hemoglobin concentration, tissue oxygenation, tissue scattering coefficient, and tissue water content. Additionally, these autoregulation metrics are not limited to the brain and can be applied to other organs so long as a perfusion metric and metabolic metric is obtainable. Moreover, in all organs, including but not limited to the brain, a perturbation (e.g. pharmacologic intervention or modification of the oxygen or carbon dioxide level) can be introduced to determine the change in autoregulation metrics at baseline and during a perturbation to improve assessment of autoregulation or lack thereof.

[0065] In some embodiments, the system may additionally include one or more electroencephalography (EEG) electrodes, which may extend from the device body and may allow for co-localized EEG monitoring. As a non-limiting example, the system may be configured to detect an EEG signal and evaluate neurovascular coupling. Alternatively, detection of the EEG signal may allow for quantification of cerebral electrical activity via parameters including, but not limited to, root-mean-squared (RMS) intensity, information quantity (IQ) and other entropy-based measures, burst frequency, coherence, and phase-amplitude coupling, either in a specific sub-band (or a ratio of two sub-bands) or over a wider range of frequencies, using either values at individual time points or changes in these values over a period of time. A highly robust EEG signal, especially with certain quantification metrics, can signify a highly activate brain, which in turn would be expected to have a high blood flow, suggesting adequate neurovascular coupling. On the other hand, if neural activity and blood flow are not matched accordingly, this may suggest neurovascular decoupling and the possibility of ischemia or overperfusion, both of which can be pathologic.

[0066] In one embodiment the system may continuously calculate the quantitative value of the autoregulation metric in real time. As used herein, the term real time refers to a very short time delay (e.g., from less than a second to several minutes, depending on the embodiment) between data acquisition and display of measured parameters. This may advantageously allow for rapid monitoring, because measuring these types of parameters in absolute quantitative units enables characterization of the tissue in a very short period of time, without need for a perturbation challenge that would require more time to administer and more stress to the patient. Immediate assessment of a patient by a first-responder being called for a medical emergency (e.g. stroke, trauma, hemorrhage, cardiac arrest, etc.) would highly benefit from real-time monitoring to enable healthcare providers to make an accurate diagnosis and thereby mobilize an accurate and rapid treatment. This may also advantageously allow for real-time feedback to inform clinical treatment by monitoring the effects of a clinical intervention (e.g., raising or lowering blood pressure, administering oxygen) during and immediately after it is performed. In preferred embodiments, the system is non-invasive, detects an intrinsic optical signal, and does not require any exogenous analyte or contrast agent. Here, an intrinsic signal is defined as a signal that is inherent to the tissue.

[0067] According to one embodiment, the present invention features a device for quantitative intracranial measurement of a brain metric. As a non-limiting example, the device may comprise: a device body; one or more light sources; and one or more detectors. The device may additionally comprise an operatively connected microprocessor and an operatively connected memory component capable of executing instructions held or stored in the memory component. In some embodiments the light sources and the detectors may extend from the device body and be positioned in proximity to a head of a subject. In some embodiments, the light sources may be configured to emit a coherent light signal (either spatially or temporally coherent) and a modulated light signal (including, but not limited to, intensity/amplitude modulation, phase modulation, frequency modulation, temporal modulation, and/or spatial modulation). The coherent light signal may allow for measurement of a dynamic perfusion metric and the modulated light signal may allow for measurement of a tissue absorption metric and a tissue scattering metric. The detectors may be configured to detect one or more backscattered light signals. In a preferred embodiment, the backscattered light signals may allow for determination of an absolute value of the brain metric using the dynamic perfusion metric, the tissue absorption metric and the tissue scattering metric. This determination may be accomplished by the execution of decoupling instructions which are held in the memory component. All of these metrics can change over time and these changes can be monitored by this invention.

[0068] Non-limiting examples of dynamic perfusion metrics include speckle flow index (SFI), blood flow index (BFI), Brownian diffusion coefficient (Db), and directed flow speed (vc). Non-limiting examples of tissue absorption metrics include tissue absorption coefficient at different wavelengths, oxygenated hemoglobin concentration, deoxygenated hemoglobin concentration, total hemoglobin concentration water content, and lipid content. Non-limiting examples of tissue scattering metrics include tissue scattering coefficient and reduced scattering coefficient at different wavelengths, scattering amplitude, and scattering slope.

[0069] In some embodiments, the device may comprise two or more detectors with different source-detector separations, and the different source-detector separations may allow the device to distinguish between signals from different depths, for example, between signals from the scalp/skull and signals from the brain. Alternatively, the device may distinguish between signals from the skull and signals from the brain by using (either individually or in combination) different wavelengths, different modulation frequencies, different angles of the source and detector fibers, or time-gating approaches. Location of the source and detector on particular areas of the skull may also enable optimization of signal detection. For example, the temporal bone is thinner. Further, the orbital socket, including the eye and, in particular the retina, can also enable better detection of an optical signal that can provide direct or indirect data about the brain.

[0070] In one embodiment, the brain metric may be indicative of brain perfusion, oxygenation, metabolism, or cerebral edema. As a non-limiting example, the brain metric may comprise cerebral metabolic rate of oxygen (CMRO.sub.2), cerebral blood flow (CBF), tissue concentration of deoxy-hemoglobin (ctHb), tissue concentration of oxygenated hemoglobin (ctHbO.sub.2), tissue oxygenation (StO.sub.2) or a combination thereof. In further embodiments, the dynamic perfusion metric, a tissue absorption metric, or a tissue scattering metric may provide information on neuronal injury, edema, sickle cell disease, depolarization, seizure activity, pharmacologic changes, ischemia, hypoxia, metabolic injury, impaired autoregulation, or a combination thereof.

[0071] The light sources of the present invention may emit one or more light signals. As non-limiting examples, the light signals may include coherent light, modulated light, light at multiple wavelengths, light delivered at different spatial locations, light delivered at different angles relative to the tissue, light delivered to the tissue with different spatial or temporal gates relative to its detection, or a combination thereof. The coherent light signal and the modulated light signal may comprise a single coherent modulated light signal or may comprise separate light signals. In one embodiment, the device may additionally comprise an electroencephalography (EEG) electrode for co-localized EEG monitoring. In some embodiments, concurrent EEG monitoring and optical monitoring may allow for measurement of a time offset between corresponding neurologic and hemodynamic activity. This offset may indicate inefficiencies and/or delays in neurovascular coupling and could potentially be used to quantify disruption in cerebral autoregulation. This may also help in determination of cause and effect. For example, a neuronal dysfunction (e.g. stroke or seizure) may lead to subsequent vascular changes (e.g. hyperperfusion), or a vascular change (e.g. thrombosis) can lead to neuronal dysfunction (e.g. stroke). Thus, offset time can be very helpful for a clinician in the process of a diagnostic workup so that the correct treatment can be provided.

[0072] The device may additionally comprise a therapeutic probe for to administration of a colocalized therapy. As a non-limiting example, the therapy may comprise colocalized transcranial light therapy, photobiomodulation, magnetic stimulation, electrical stimulation, or another therapy. In some embodiments, determination of the absolute value of the brain metric, or a change in the value of the brain metric, may guide the therapy in real time.

[0073] In some embodiments, the device may use diffuse correlation spectroscopy (DCS), frequency-domain diffuse optical spectroscopy (FD-DOS), or a combination thereof, at one or more wavelengths in a range comprising the visible, near-infrared (NIR), and short-wave infrared (SWIR) regimes. DCS may be used to quantify blood flow. FD-DOS may be used to separate tissue absorption and scattering coefficients by using modulated light.

[0074] According to one embodiment, the sources and detectors may be attached to the device body via a plurality of legs or extendable support fibers which may be retracted into the device body. Legs may also be more flexible (e.g. made of a pliable material, such as rubber) or contain flexible optical fibers made of a pliable material that can be extended outward to reach multiple different positions on a patient and then retracted back into the body of the device following completion of the measurement. Additionally, the legs may simply be wires that are pulled out (e.g. unrolled from a rotating device) and extended from the hinge region, locked in length to attach to the target tissue site by sticky material to obtain the signal and then retracted back to its original position when data acquisition and analysis is complete. The latter would also preclude having to hold the device by hand to ensure stability of the data acquisition and thereby allowing the possibility that the device can be clipped onto a patient's clothes or rested next to them to allow for more prolonged data acquisition if desired. The extendable support fibers may be flexible or rigid, and may be designed to set a wide range of different spacings (using retractable flexible fibers or flexible legs) between the various sources and detectors and an angle between each source and detector and the surface of the head. In one embodiment, two or more of the detectors may have different source-detector separations. Without wishing to limit the present invention to any particular theory or mechanism, this may allow for differentiation of signals due to different tissues or materials at different depths. In some embodiments, the various source-detector separations may be measured by a distance sensor, a camera, a ruler, or by a patch which guides placement of the sources and detectors.

[0075] In an alternative embodiment, the present invention may feature a device for quantitative subdermal measurement of a tissue metric. As a non-limiting example, the device may comprise: a device body; one or more light sources; and one or more detectors. The device may additionally comprise a microprocessor, operatively connected to the light sources and the detectors, and a memory component, operatively connected to the microprocessor such that the microprocessor is capable of executing instructions held in the memory component. The light sources and detectors may extend from the device body so as to be positioned in proximity to a body surface of a subject, for example, on a limb or on a torso, or within an endoscope, over the position of an organ of interest. This may also include assessment of the circulatory system (e.g. arteries or veins) in various regions of the body to assess for blood flow and regional metabolism. In an alternative embodiment, the device may be an endoscopic device which is configured to position the sources and detectors within the body. In one preferred embodiment the light sources may be configured to emit a coherent light signal and a modulated light signal. In another preferred embodiment, the detectors may be configured to detect one or more backscattered light signals.

[0076] In some embodiments, the modulated light signal may allow for decoupling of a plurality of components of the backscattered light signals. This decoupling may be accomplished via the microprocessor's execution decoupling instructions which are held in the memory component. As a non-limiting example, the modulated light may enable decoupling of tissue absorption and scattering coefficients. In one embodiment, the components may comprise: a tissue absorption component; a tissue scattering component; and a dynamic scattering flow component. In further embodiments, the decoupled components of the backscattered light signals may allow for determination of an absolute value of the tissue metric. The absolute value of the tissue metric may provide information on the perfusion or metabolism of an organ, or may be indicative of tissue autoregulation. In one embodiment, the device may allow for comparative analysis of the autoregulation of two or more body parts. For example, on COVID-19 patients, the device could measure the autoregulation of the lung versus the brain to quantify the degree to which respiratory impairment is affecting the brain. This is important because in COVID-19, there can be silent hypoxia, wherein the respiratory impairment may lead to low system oxygen levels while the brain is compensating adequately. This may argue against maximal oxygen therapy in the lung by mechanical ventilation, which can cause ventilator-induced lung injury. Thus, understanding the coupling and decoupling of flow-metabolism dynamics between organs can help a clinician tailor the medical therapy for a particular patient's condition in a more precision guided manner without using a suboptimal one-size-fits-all approach.

[0077] In another embodiment, the present invention may feature a portable device for measuring EEG and blood flow simultaneously. As a non-limiting example, the device may include: a body; a signal processing component; a plurality of legs comprising: at least one reference electrode leg, at least one measurement electrode leg, at least one leg comprising an optical source fiber, and at least one leg comprising an optical detection fiber, wherein each leg is pivotably attached to the body; a light source disposed within the body attached to each leg comprising an optical source fiber; a photodetector disposed within the body attached to each leg comprising an optical detection fiber; and a hardware or software correlator disposed within the body communicatively coupled to the photodetector.

[0078] In some embodiments, the device may include a signal processing component that analyzes a signal received by a measurement electrode. In other embodiments, the device may include a signal processing component that analyzes a signal generated by the light source and received by the photodetector. In still other embodiments, the device may include one or more components are wirelessly, operatively connected to a display, where the display shows the data obtained from the components. In some embodiments, the optical detection fiber may comprise an optical or electronic filter that allows for separating coherent light from incoherent light. In some other embodiments, the device may include an on-board or off-board correlator for signal processing and analysis of pulsatile components in cerebral blood flow or cardiovascular blood flow.

[0079] According to one embodiment, the device may additionally comprise: one or more oximetry optical source fibers; one or more oximetry optical detection fibers; one or more light sources connected to the oximetry optical source fibers; and one or more oximetry photodetectors connected to one or more oximetry optical detection fibers. In some embodiments, the first light source may deliver light at a lower-energy range and the other light sources deliver light at higher-energy ranges. In one embodiment, the optical source fibers may be disposed in a single leg. In another embodiment, the optical detection fibers may be capable of collecting electrical signals and converting them to optical signals using an electrical-optical transducer. In still another embodiment, the oximetry optical detection fibers may be connected to one or more optical-electrical transducers. In yet another embodiment, the device may feature an adjustable clamping mechanism for retracting and extracting optical or electronic fibers.

[0080] In one embodiment, the present invention may feature a portable device for therapeutic photobiomodulation and blood flow, oximetry, and/or electrical activity (e.g. from brain or other tissue) measurement. As a non-limiting example, the device may comprise: a body; a plurality of legs comprising: one or more legs having one or more optical source fibers, wherein each leg is pivotably connected to the body by a hinge. one or more legs having one or more optical detection fibers, one or more photobiomodulation light sources connected to one or more optical source fibers; one or more non-photobiomodulation light sources connected to one or more optical source fibers; and one or more photodetectors connected to one or more detection fibers. In some embodiments, the photobiomodulation light sources may produce light capable of therapeutic photobiomodulation of biological tissue. In other embodiments, the non-photobiomodulation light sources may produce light for blood flow, oximetry, and/or electrical activity (e.g. of brain or other tissue) measurements. In still other embodiments, the detection fibers may receive backscattered light produced by the optical source fibers.

[0081] In one embodiment, the photobiomodulation light sources may produce a wide range or specific range of light including but not limited to visible light, near-infrared light, short-wave infrared light, or infrared light. In another embodiment, the photobiomodulation light sources may promote a variety of perfusion and metabolism changes including but not limited to increased cellular metabolism of oxygen, generation of ATP, vasodilation, and enhanced perfusion of tissue. In still another embodiment, the changes produced by the photobiomodulation light sources may be monitored using methods including but not limited to real time detection or delayed detection by the photodetector.

[0082] In some embodiments, the device may have additional legs for delivering an electrical current and/or a magnetic field produced by a plurality of microelectronics, including but not limited to MEMS, disposed within the body for therapeutic electromagnetic cell and tissue modulation. The electrical current and/or magnetic field may promote changes including but not limited to perfusion, metabolism, and electrical activity. In some embodiments, the changes produced by the electrical current and/or magnetic field may be monitored using methods including but not limited to real time analysis by the photodetector and associated electronics.

[0083] According to some embodiments, the present invention features a device with capabilities for optical and electrical (i.e. EEG) monitoring. In other embodiments, the present invention features a device with only optical monitoring capabilities. In still other embodiments, the present invention features a device with only electrical monitoring capabilities. While for some situations, it may be desirable to have both optical and electrical data, for many situations, electrical data alone may be sufficient. For example, the electrical monitoring herein is capable of measuring an evoked potential wherein the stimulus used to evoke a potential may also be provided by the said device or another device. Here, an evoked potential may include, but is not limited to, a visual evoked potential, somatosensory evoked potential, auditory evoked potential, steady-state evoked potential, laser evoked potential, motor evoked potential, evoked compound motor action potential or sensory nerve action potential as seen in nerve conduction studies. As such, this optical EEG pen may be used to rapidly conduct an evoked potential without the need for a bulky equipment and significant preparation time. For example, for each type of evoked potential measurement, the optical EEG pen itself may be able to provide the stimulus (e.g. light pattern for a visual evoked potential, specific frequency and type of sound for an auditory evoked potential, an electrical impulse or shock to a peripheral nerve for a somatosensory evoked potential, etc.). These evoked potentials may be induced by an extension (wired or wirelessly) of the optical EEG pen so that the main body of the optical EEG pen itself may still be touching or attached to the head so that the EEG signal and evoked potential change in the brain can be measured. If necessary, the evoked potential can utilize the optical portions of the optical EEG pen also (e.g. light to conduct a visual evoked potential). However, only the electrical components of the optical EEG pen may be utilized to either record the evoked potential induced by a separate device or the optical EEG pen itself (e.g. an electrical impulse). Similarly, either by itself or with additional components, the optical EEG pen may also allow for an event-related potential to be measured. As such, by allowing for evoked potentials or event-related potentials to be tested, the optical EEG pen's electrical data alone may enable significant diagnostic workup and interrogation of the nervous system, including lesions along the sensory or motor pathways, including peripheral nerves, the spinal cord, brainstem, subcortex, or the cortex. Such lesions may be induced through trauma, surgical procedures, drug-induced dysfunction, or insults. This can also be used to assess for neuroplasticity and monitoring of response to treatments geared towards improvement of a pathway. This may also be used for intraoperative monitoring of the neuroaxis. Current systems that may allow for monitoring and testing are bulky and take time whereas the optical EEG pen will enable rapid analysis of the neuroaxis. Moreover, the optical, electromagnetic, and photobiomodulation aspects of the optical EEG pen can enable for far more diagnostic and therapeutic measures to be conducted. For example, if a lesion is localized in the neuroaxis, treatment may be implemented using these other features of the optical EEG pen described elsewhere here.

[0084] Referring now to FIG. 1A, the present invention may feature a portable device for measuring EEG and blood flow simultaneously. The device may comprise a body and a signal processing component. The device may further comprise a plurality of legs comprising one or more reference electrode legs, one or more measurement electrode legs, at least one leg comprising an optical source fiber, and at least one leg comprising an optical detection fiber. Each leg may be pivotably attached to the body by a hinge or another type of pivotable attachment. The device may further comprise a light source, which may include but is not limited to visible light (e.g. 400-750 nm), near-infrared light (750-1000 nm), short-wave infrared light (approximately 1000-2000 nm) disposed within the body at the hinge of each optical source fiber leg, a photodetector disposed within the body at the hinge of each optical detection fiber leg, and a hardware or software correlator disposed within the body communicatively coupled to the photodetector. In some embodiments, the near-infrared light source may additionally be a coherent light source. The signal processing component may analyze a signal received by a measurement electrode and may analyze a signal generated by the coherent near-infrared light source and received by the photodetector. FIG. 1A shows an embodiment of the Optical EEG pen invention that includes technology for simultaneous measurement of EEG activity and cerebral blood flow. The pen body may contain electronic components for receiving the detected electrical signal from the electrodes. The pen body may also include a coherent light source in the near-infrared range (or visual or short-wave infrared range) that is directed onto the tissue and one or more photodetectors for detecting the optical signal that returns to the surface to measure blood flow. In some embodiments, five or more legs are connected to the pen body via hinges (or another attachment) to allow flexibility in their positioning. Two or more legs may contain measurement (recording) electrodes, one or more legs may contain a reference electrode, one or more legs may contain an optical fiber for delivering coherent light to the tissue, and one or more legs may contain an optical fiber for detecting the light that returned to the surface after delivery of the light by the fiber and propagation of the light through the tissue. In other embodiments, the pen body may be replaced with a clip version with retractable fibers or electrodes that do not need to be held by the hand but rather attached to a surface. In yet other embodiments, the device can be used to obtain measurements from multiple different regions of the brain or bodily tissue to obtain a map of the endpoints measured, including but not limited to a 3-dimensional representation of the brain, heart, lungs, other internal organs, limbs, circulatory system (e.g. arteries or veins), and/or other body parts. As such, this invention need not be limited to a pen body but rather other embodiments of a portable device that are wireless or wired.

[0085] Referring now to FIG. 18, the optical detection fiber may comprise an optical filter (e.g., a spectral filter) for separating coherent light from incoherent light and an on-board or off-board software correlator for signal processing and analysis of pulsatile components in cerebral blood flow. FIG. 1B shows an embodiment of the invention that enables faster on-board or off-board data processing by using a custom-built hardware correlator or software correlator for faster signal processing, enabling real-time on-board or off-board acquisition and analysis of data related to pulsatile components of cerebral hemodynamics.

[0086] FIG. 1C shows an additional embodiment where the correlator could be either a hardware or a software correlator and the data processing can either be on-board or off-board.

[0087] FIG. 1D shows an additional embodiment where multiple detector fibers can be used to detect optical signals from the tissue. The use of multiple detector fibers allows for multiple source-detector separations (e.g., r1 and r2), which allow the device to probe different depths (z) beneath the surface of the tissue to interrogate different regions of the brain or other bodily tissue and/or better separate the signal from the brain or other bodily tissue from that of the skin, subcutaneous tissue, and/or bone. This embodiment can also be applied to all other embodiments described in this invention.

[0088] Referring now to FIG. 2 the reference electrode leg and each measurement electrode leg may further comprise an electrode socket for connecting to recording sites. FIG. 2 shows an embodiment of the invention that shows how the hinges or a related component connect the legs to the body, as well as showing how the tips of three of the legs contain electrode sockets that contact recording sites on the head to record EEG signals from the brain. Of note, signal detection with this invention need not be limited to detection of brain signals. Rather, they can apply to the detection and analysis of cardiac signals (e.g. electrocardiogram; ECG) or muscle (e.g. electromyogram) also if applied to different regions of the body with appropriate signal filtering, calibration, referencing, and analysis. This can be extremely helpful during a medical emergency such as a cardiac arrest wherein very fast deployment of the optical EEG pen onto the chest or other body part to enable ultra-rapid diagnosis of the cardiac rhythm much faster than conventional methods can be life-saving (i.e. shockable versus non-shockable rhythm).

[0089] The electrode socket may comprise a base having a first side and a second side, a joint disposed on the base at the first side and pivotably connected to the leg, a plurality of extrusions disposed on the second side, and a recording site disposed on each extrusion of the plurality of extrusions. The recording site may be coated in Ag or AgCl or another type of material that can maximize detection of electrical signals while minimizing electrical noise. In one embodiment of the invention, the electrode sockets at the tips of the three EEG legs may be composed of small extrusions that serve as recording sites.

[0090] The recording sites are coated with Ag/AgCl or other type of surface to enhance conductivity of detected electrical signals and allow the possibility of AC and DC signals to be measured. The electrodes can be of the wet or dry type depending on the situation and the target organ (e.g. brain EEG signal or heart ECG signal) being measured that will optimize usage, signal detection, indication of need, and cost based on advantages and disadvantages of different types of electrodes.

[0091] FIG. 3 shows an embodiment of the signal processing component of the EEG portion of the invention. At each time point, the detected signals from the two measurement electrodes and the reference electrode are filtered to remove noise or to narrow the signal down to the desired range (e.g. for the clinical indication), calibrated by using the signal from the reference electrode (via common-average referencing or a similar technique). Next, the processed signals are analyzed on-board or off-board via microelectronics, including but not limited to MEMS or nanoelectronics, via methods including but not limited to whole-band EEG, frequency-based EEG (e.g., frequency sub-band raw signal or RMS and amplitude-based measurements by frequency sub-band), entropy-based metrics (e.g., information quantity), synchrony-based metrics (e.g., coherence, phase-amplitude coupling), and seizure detection algorithms. Then, the EEG signal and parameters from the data analysis are read out in real time. Of note, these electrodes are not limited to detection of EEG signals. Rather, they may also be adopted for other purposes, including but not limited to electrocardiogram (ECG) signals and electromyogram (EMG) signals.

[0092] FIG. 4 shows an additional embodiment where the data analysis algorithm can be on-board or off-board and can include a seizure detection algorithm.

[0093] FIG. 5A shows an embodiment of the signal processing component of the optical portion of the invention. At each time point, the detected signal is obtained from light from the coherent source entering the tissue, propagating through the tissue, and eventually returning to the surface and entering the photodetector. The detected signal may be obtained from backscattered light. This signal is then postprocessed via filtering and calibration methods to separate the coherent portion of the detected light and correct the signal for the effect of the skull on signal data to separate out the signal component from the brain. Then, data analysis is performed to extract the measured autocorrelation function from brain tissue and fit that function to a mathematical model to extract the measured blood flow in the tissue, which is read out in real time. With frequency modulation of the light amplitude, the detected light can be analyzed to separately determine the absorption and scattering properties of the interrogated tissue. This approach can be used to improve accuracy of optical oximetry. Scattering measurements provide an additional modality to study bulk changes in the interrogated tissue that may occur due to, for example, neuronal changes associated with but not limited to epileptic seizures, spreading depolarization, ischemia, or neurodegenerative changes.

[0094] FIG. 5B shows an additional embodiment where the data analysis algorithm can be on-board or off-board.

[0095] Referring now to FIG. 6A, the device may further comprise one or more oximetry optical source fibers and one or more oximetry optical detection fibers, including but not limited to a first LED connected to the first oximetry optical source fiber, a second LED connected to the second oximetry optical source fiber, and an oximetry photodetector connected to the oximetry optical detection fiber. The first LED or other type of light source may deliver light at a lower-energy range of a visible spectrum, near-infrared spectrum, or short-wave infrared spectrum, and the second LED or other type of light source may deliver light at a higher-energy range of a visible spectrum, near-infrared spectrum, or short-wave infrared spectrum. FIG. 6A shows an embodiment of the invention where three additional legs have been added to the instrument for optical oximetry measurements. Each leg contains an optical fiber; two of these fibers deliver light (e.g., from LEDs) to the tissue. One fiber delivers light at the higher-energy range of the visible spectrum (e.g., 655 nm) for sensitivity to deoxyhemoglobin, and the other fiber delivers light at the lower-energy range of the near-infrared spectrum (e.g., 850 nm) for sensitivity to oxyhemoglobin. These specific wavelengths serve as examples and can be modified depending on the specific component of interest being measured at a specified target organ, underscoring the versatility and flexibility of this invention to enable customization of the optical EEG pen for different users and purposes. Moreover, light from both sources propagates through the tissue until a portion of this light is detected by the detection fiber. The portion of light may be backscattered. The detected light is sent to a photodetector (e.g., APD, PMT).

[0096] Referring now to FIG. 6B, the first and second oximetry optical source fibers may be disposed in a single leg. FIG. 6B shows an embodiment of the invention where the two non-coherent light sources (e.g., 655 nm, 850 nm) are delivered to the tissue via a single leg of the instrument, by using a miniature multiplexer for rapidly switching between the two light sources (e.g., LEDs or laser diodes) or modulating the two light sources at different frequencies to encode them separately.

[0097] Referring now to FIG. 6C, an oximetry optical source fiber may be capable of delivering electrical signals and the oximetry optical detection fiber may be capable of collecting electrical signals. The oximetry optical detection fiber may be connected to an electrical-optical transducer and/or an optical-electrical transducer. FIG. 6C shows an embodiment of the invention where the electrical and optical signals are delivered and collected via a pair of common fibers, with the source fiber connected to an electrical-optical transducer for converting electrical inputs into optical inputs, and with the detection fiber connected to an optical-electrical transducer for converting detected optical signals into electrical signals. Note that the number of fibers is not limited to two, as additional fibers can be added to increase the number of measurement/recording EEG electrodes as well as for additional optical signals as needed.

[0098] FIG. 6D shows an embodiment of the invention where the optical fibers can be retracted into the pen body and extracted out of the pen body via an adjustable clamping mechanism.

[0099] FIG. 7A shows an embodiment of the signal processing component of the additional optical portion of the embodiment of the invention shown in FIG. 6. At each time point, the detected signal is obtained from light from the two additional sources entering the tissue, propagating through the tissue, and eventually backscattering to the surface and entering the additional photodetector. This signal is then postprocessed via filtering and calibration methods specific to the wavelengths and source-detection separations, and the signal is corrected for the effect of the scalp/skull to separate out the component from the brain. Then, on-board data analysis is performed to fit the calibrated data to a light propagation model to extract the relative values of tissue oxygenation and hemoglobin concentration, which are read out in real time. In some embodiments, the invention may further comprise a frequency domain photon migration (FDPM) component for extracting the tissue absorption and scattering coefficients from the data using a light propagation model. These measured parameters can be combined with the blood flow parameter to measure cerebral oxygen metabolism, concomitant with cerebral blood flow, in real time.

[0100] The blood flow and oxygenation parameters can also be combined with the measured EEG parameters to measure neurovascular coupling in real time. In addition to tissue scattering metrics obtained from FDPM, these metrics can also be used to measure cortical spreading depolarization during diagnostic, prognostic, or therapeutic workup of brain injury. Furthermore, these metrics can be used to quantify the metabolic changes and/or neurovascular coupling changes during a variety of metabolic states, including the normal fed state, calorically restricted state, high versus low metabolic state (e.g. during exercise or stress tests conducted for assessment of particular organs), or pathological states such as hypoxia, ischemia, post-ischemia, reperfusion, and acute or delayed injury, as well as degenerative states such as neurodegenerative conditions (e.g. dementia). Application of the device during these states can be for the purpose of diagnostic, prognostic, or therapeutic purposes. Additionally, application of the device can include optimization of performance for the brain and heart to improve health and fitness. Interpretation of data can be made based on relative or absolute changes in the aforementioned metrics for a single person or a multitude of people using population databases.

[0101] FIG. 7B shows an additional embodiment where the data analysis algorithm can be on-board or off-board. Referring now to FIG. 8, the present invention features a portable device for therapeutic photobiomodulation and blood flow measurement. The device may comprise a body. The device may further comprise a plurality of legs comprising a first leg having a first optical source fiber, a second leg having a second optical source fiber, and a third leg having a detection fiber. Additional legs can be added as needed. Each leg may be pivotably connected to the body by a hinge or related type of attachment. The device may further comprise a photobiomodulation light source connected to the first optical source fiber, a non-photobiomodulation light source connected to the second optical source fiber, and a photodetector connected to the detection fiber. The photobiomodulation light source may produce light capable of therapeutic photobiomodulation of cells, the non-photobiomodulation light source may produce light for blood flow, oximetry, and metabolism measurement, and the detection fiber may receive backscattered light produced by the first and second optical source fibers. The photobiomodulation light source may produce visible light, near-infrared light, or short-wave infrared light to promote increased cellular metabolism of oxygen, generation of ATP, vasodilation, and/or enhanced perfusion of tissue or other biological change desired to, for example, promote healing or rapid medical intervention. Changes produced by the photobiomodulation light sources may be monitored in real time or delayed fashion by the photodetector. FIG. 8 shows an embodiment of the invention where one of the optical light sources (e.g., an LED, laser diode, or other light source) is appropriately optimized (tuned) via modifying a number of parameters including but not limited to wavelength, frequency, and power to deliver visual light, near-infrared (e.g., 850 nm) light, or short-wave infrared light to the tissue for therapeutic photobiomodulation. The light may also have additional features, such as being pulsed or dynamically changing in space and time in the said parameters to achieve the desired goal. A portion of this light that is not absorbed by the tissue is scattered back to the surface (as is the case with the non-photobiomodulating light source), so the backscattered light from both sources is still detected by the detection fibers and sent to the photodetector to perform real-time measurements of absolute and/or relative values of tissue perfusion, oxygenation, and metabolism.

[0102] FIG. 9 shows a flow chart of an embodiment of the invention where visual light, near-infrared light (e.g., 850 nm), or short-wave infrared light is sent to the tissue for photobiomodulation. The input characteristics of this light (e.g., wavelength, power, frequency, etc.) are tuned for optimal uptake by cellular mitochondria or other biological components. A portion of the light energy is absorbed by the mitochondria, promoting increased cellular metabolism of oxygen and generation of ATP, or other biological components in addition to vasodilation and enhanced perfusion of the tissue via release of nitric oxide (NO). On the other hand, if vasoconstriction or reduced perfusion is desired as in the case of pathological hyperperfusion (e.g. reperfusion injury), such modulation can be implemented. These and potentially other factors can contribute to neuronal repair/healing. These healing processes enhanced by the photobiomodulation are manifested in changes in the blood flow, oxygenation, metabolism, tissue scattering, and electrical activity of the end organ (e.g. EEG in the brain, ECG in the heart, or EMG in the muscle). Changes in these parameters can be monitored in real time by the optical and electrical measurement capabilities of the instrument (via simultaneous real-time detection of backscattered coherent light, detection of backscattered incoherent light, and detection of EEG activity). These measurements can be used to quantify the effect of the photobiomodulation treatment on cerebral hemodynamic, metabolic, and electrical activity, and they can be used to provide feedback to modify the characteristics of the photobiomodulation source (e.g., wavelength, power, frequency) in real time to facilitate improved treatment.

[0103] Referring now to FIG. 10, the device may further comprise a fourth leg for delivering an electrical current and/or a magnetic field produced by a plurality of microelectronics, including but not limited to MEMS, disposed within the body for therapeutic electromagnetic biomodulation. The electrical current and/or magnetic field may promote perfusion, metabolism, and electrical activity and changes produced by the electrical current and/or magnetic field may be monitored in real time by the photodetector. FIG. 10 shows an embodiment of the invention where one or more additional legs are added to deliver an electric current and/or magnetic field to the tissue to stimulate repair/healing. This electromagnetic (EM) stimulation may be performed using miniaturized electronics in the pen body to generate electric currents and/or magnetic fields whose characteristics (e.g., power, frequency) are optimized to promote uptake by the tissue for therapeutic purposes, and transmitting the currents and/or fields through the additional legs of the device for delivery to the brain tissue. The additional legs may also have a hinge or related component like the other legs for adjustable positioning.

[0104] FIG. 11 shows a flow chart of an embodiment of the invention where an electric current and/or magnetic field is sent to the tissue for electromagnetic (EM) stimulation (e.g., direct current stimulation, transcranial magnetic stimulation). The input characteristics of this current and/or field (e.g., power, frequency, etc.) are tuned for optimal uptake by the brain tissue. A portion of this EM energy is absorbed by the brain or other biological tissue, promoting increased perfusion, metabolism; and electrical activity. These and potentially other factors contribute to neuronal repair/healing. These healing processes enhanced by the EM stimulation are manifested in changes in the blood flow, oxygenation, metabolism, tissue scattering, and electrical activity in the brain or other tissue. Changes in these parameters can be monitored in real time by the optical and electrical measurement capabilities of the instrument (via simultaneous real-time detection of backscattered coherent light, detection of backscattered incoherent light, and detection of electrical activity). These measurements can be used to quantify the effect of the EM stimulation treatment on hemodynamic, metabolic, and electrical activity, and they can be used to provide feedback to modify the characteristics of the EM source (e.g., power, frequency) in real time to facilitate improved treatment.

EXAMPLES

[0105] The following are non-limiting examples of the present invention. It is to be understood that said examples are not intended to limit the present invention in any way. Equivalents or substitutes are within the scope of the present invention.

Example 1: Use of a Device of the Present Invention to Measure an Absolute Value of CMRO.SUB.2

[0106] Step 1: Extend a light source leg and a detector leg from the portable device.

[0107] Step 2: Illuminate the tissue with the coherent and incoherent light source(s), and detect the backscattered light with the photodetector(s). One or more of these light sources may be modulated at different frequencies, and there may be multiple sources, detectors, and wavelengths of light.

[0108] Step 3: Perform calibration and data processing procedures to obtain a blood flow index from the coherent light measurement and obtain a tissue absorption and scattering coefficient from the incoherent light measurement. The tissue absorption coefficient can then be fit to a linear combination of the different endogenous chromophores in the tissue (e.g., oxygenated hemoglobin, deoxygenated hemoglobin, water, lipid) to obtain the concentration or fractional contribution of each chromophore.

[0109] Step 4: Use the measured tissue absorption and scattering coefficient along with the blood flow index to obtain a corrected value of blood flow in physiological units.

[0110] Step 5: Combine the corrected blood flow data with the hemodynamic data to obtain an absolute value of tissue metabolic rate of oxygen in physiological units.

Example 2: Use of a Device of the Present Invention to Quantify Cerebral Autoregulation

[0111] Step 1: Perform the process described in Example 1 to obtain a quantitative value of blood flow and a quantitative value of tissue metabolic rate of oxygen.

[0112] Step 2: Calculate an autoregulation metric by dividing the quantitative value of the absolute perfusion metric by the quantitative value of the absolute metabolic metric. Alternatively, the quantitative value of cerebral autoregulation may be calculated by dividing the rate of change of an absolute or relative perfusion metric by the rate of change of an absolute or relative metabolic metric, or by other methods described above (in claim [0043]).

Example 3: Use of a Device of the Present Invention to Monitor a COVID-19 Patient

[0113] Step 1: Measure blood flow, metabolism, and/or an autoregulation metric for the brain with the device, using the method described in Example 2.

[0114] Step 2: Measure blood flow, metabolism, and/or an autoregulation metric for a different part of the body (e.g., lung or limb) with the device, using the method described in Example 2.

[0115] Step 3: Calculate a ratio of, or a difference between, the measured blood flow, metabolism, and/or autoregulation metrics for the two different parts of the body.

[0116] Step 4: Compare this metric with a previously measured range for healthy subjects and/or a previously-obtained baseline value of the metric for the same patient (e.g., a value measured upon admission to the hospital) to assess whether the patient exhibits a deficit in blood flow, metabolism, and/or autoregulation between two different parts of the body.

[0117] Step 5: Monitor this metric periodically over time, including in response to clinical treatment, to assess the recovery of the patient and the effectiveness of the clinical intervention. As a simple example, if the COVID-19 patient has thromboses in certain body parts or end organs demonstrated by a reduction in blood flow, a treatment such as anticoagulation can be initiated to treat the potential ischemia or risk of embolization that can lead to a pulmonary embolism or stroke. As a more advanced example, if the blood flow, metabolism, or autoregulation metric between the brain and lung are disparate, including for example a hypoxic lung but normoxic brain, then the ventilator settings can be adjusted so as to avoid excessive oxygen and air pressure that could damage the lung while avoiding a potential stroke or long term hypoxic brain damage. At the best case scenario, such a patient might avoid needing endotracheal intubation and mechanical ventilation altogether and instead may tolerate non-invasive ventilation such as a high-flow nasal cannula, CPAP or BiPAP device because the brain, the organ that is most sensitive to hypoxia and ischemia, might not benefit from additional oxygen delivered by invasive mechanical ventilation. Avoidance of mechanical ventilation can help to save the life of another patient due to a potential shortage of ventilators during a pandemic such as COVID-19.

Example 4: Use of a Device of the Present Invention to Diagnose, Treat, and Monitor an Acute Brain Injury

[0118] Step 1: When encountering a patient with altered mental status with the differential diagnosis of seizure, stroke, or drug overdose, the device would be deployed to the brain to assess cerebral electrical activity (i.e. EEG), blood flow, metabolism, and/or an absolute CMRO.sub.2 as stated in Example 1.

[0119] Step 2: The information provided by the device would help with the diagnosis of the patient. For example, epileptiform discharges seen on the EEG data would suggest possible ongoing seizure activity, diminished cerebral blood flow in a particular part of the brain that corroborates with contralateral motor weakness in a limb might suggest an ongoing ischemic stroke, or a particular EEG signal (e.g. triphasic waves) or a diminished EEG signal or metabolic activity throughout the whole brain might support a drug overdose (e.g. benzodiazepine or opioid overdose).

[0120] Step 3: In each of these cases, a treatment can be administered. The treatment might be pharmacologic in nature or the device itself can be utilized to induce a rapid optical wave (e.g. photobiomodulation) or electromagnetic wave to treat the condition. As a hypothetical example, the said treatment might result in treatment of the seizure (e.g. by rapid induction of cortical spreading depolarization), stroke (e.g. optical or electromagnetic thrombolysis), or drug overdose (e.g. optical or electromagnetic stimulation of breathing centers to ovoid respiratory depression due to drug overdose).

Example 5: Use of a Device of the Present Invention to Diagnose, Treat, and Monitor an Acute Cardiac Event

[0121] Step 1: When encountering a patient who is suspected of a potential heart attack and CPR is being initiated by other bystanders (e.g. basic life support) or health care providers (e.g. advanced cardiac life support), the device is deployed to the chest area to assess the cardiac electrical activity (i.e. ECG). This will help determine the cardiac rhythm (e.g. non-shockable rhythm versus a shockable rhythm) to help with diagnosis of the condition.

[0122] Step 2: Treatment may be guided by the device and potentially provided by the device. For example, if the device's ECG data indicates a shockable rhythm such as ventricular fibrillation or ventricular tachycardia, a shock needs to be administered. The device itself can have the potential to deliver a shock sufficient enough to convert the cardiac rhythm to normal sinus rhythm similar to an Automated External Defibrillator (AED) except that the said device would be handheld and more rapidly deployed than an AED.

[0123] Step 3: The device can additionally be used to monitor the brain during ongoing CPR to ensure that chest compressions and administered breaths are providing sufficient blood flow and oxygen to the brain to avoid hypoxic-ischemic brain damage. Additionally, after successful resuscitation of the patient, the brain can be monitored using the device to monitor brain electrical activity (e.g. spreading depolarizations and repolarizations), blood flow, and metabolism to help guide ongoing medical therapy to optimize a favorable outcome for the patient.

Example 6: High-Speed Quantitative Optical Imaging of Absolute Metabolism in the Rat Cortex

Abstract:

[0124] Quantitative measures of blood flow and metabolism are essential for improved assessment of brain health and response to ischemic injury. This example demonstrates a multimodal technique for measuring the cerebral metabolic rate of oxygen (CMRO.sub.2) in the rodent brain on an absolute scale (M O.sub.2/min). This example uses laser speckle imaging (LSI) at 809 nm and spatial frequency domain imaging (SFDI) at 655 nm, 730 nm, and 850 nm to obtain spatiotemporal maps of cerebral blood flow (CBF), tissue absorption (.sub.a), and tissue scattering (.sub.s). Knowledge of these three values enables calculation of a characteristic blood flow speed, which in turn is input to a mathematical model with a zero-flow boundary condition to calculate absolute CMRO.sub.2. This method is applied to a rat model of cardiac arrest (CA) and cardiopulmonary resuscitation. With this model, the zero-flow condition occurs during entry into CA. The CMRO.sub.2 values calculated with this method are in good agreement with those measured with magnetic resonance (MR) and positron emission tomography (PET). Baseline absolute CMRO.sub.2 values were statistically significant for distinguishing rats that exhibited good vs. poor short-term neurological recovery, as measured by electrocorticography. This technique provides a quantitative metric of cerebral metabolism that can potentially be used for comparison between animals and longitudinal monitoring of a single animal over multiple days, to assess differences in baseline metabolism and track recovery of metabolism in survival studies following ischemia and reperfusion.

INTRODUCTION

[0125] Measurements of cerebral metabolic rate of oxygen (CMRO.sub.2) may provide insight into the viability of brain tissue following ischemia. In particular, knowledge of CMRO.sub.2 in absolute units would obviate the need for baseline measurements and facilitate longitudinal measurements to track longer-term cerebral recovery following ischemia and reperfusion. Also, absolute CMRO.sub.2 measurements would enable quantitative comparisons between the values of different subjects at baseline and at subsequent time points in preclinical or clinical studies.

[0126] Unfortunately, conventional clinical monitoring techniques (e.g., arterial blood pressure, jugular bulb oximetry, pulse oximetry, laser Doppler flowmetry) typically cannot separate alterations in cerebral metabolism from changes in blood flow. Established techniques to measure CMRO.sub.2 include medical imaging modalities such as Positron Emission Tomography (PET) and functional Magnetic Resonance Imaging (fMRI). PET can measure absolute CMRO.sub.2, but it is expensive and non-portable, and requires use of exogenous contrast agents containing radioactive tracers. fMRI measures CMRO.sub.2 changes via the blood oxygen level dependent (BOLD) signal, but it provides absolute CMRO.sub.2 only with extensive calibration, as the BOLD signal only serves as a surrogate for cerebral blood flow and hemoglobin content and not as a direct measurement of these quantities. Additionally, both PET and BOLD fMRI typically have limited temporal resolution and cannot be performed repeatedly on a patient over a short period to monitor hyper-dynamic changes caused by acute insults.

[0127] Diffuse optical spectroscopy (DOS) and diffuse optical imaging (DOI) techniques are an attractive alternative, as they are noncontact and use measurements of low-irradiance visible and near-infrared light to extract endogenous tissue absorption and scattering coefficients. Further analysis of the absorption coefficient yields measurements of relative changes in hemoglobin content and oxygen saturation. Frequency-modulated and time-resolved techniques enable absolute measurements of these two parameters. Coherent light techniques, such as diffuse correlation spectroscopy, enable measurements of blood flow. Combined use of these techniques can yield measurements of CMRO.sub.2. However, the majority of these techniques have limited spatial resolution, oftentimes serving as point measurements.

[0128] It has been demonstrated that the combination of Spatial Frequency Domain Imaging (SFDI) and Laser Speckle Imaging (LSI)) can quantify tissue metabolic changes with both high spatial and temporal resolution. High-speed LSI and SFDI has been used to measure perfusion, oxygenation, and tissue scattering in the brain in a cardiac arrest (CA) model of global cerebral ischemia. The capability of this rapid multimodal SFDI+LSI system to image blood flow and hemoglobin concentration simultaneously enables high-speed measurement of CMRO.sub.2. Furthermore, it is possible to account for the effects of time-varying tissue scattering at multiple wavelengths and the contribution of venous regions versus parenchyma when calculating CMRO.sub.2. The present example features an approach to analyze multimodal optical imaging data with a mathematical model of CMRO.sub.2 that incorporates a zero-flow boundary condition during the onset of ischemia in a CA model, to obtain the parameters necessary for absolute CMRO.sub.2 measurements.

Methods

[0129] Animal Preparation:

[0130] Ten male Wistar rats (weight 300-400 g) were imaged. Before the experiment, all subjects were endotracheally intubated under isoflurane anesthesia. Each subject had epidural screw electrodes implanted for electrocorticography (ECoG) and a hemicraniectomy (4 mm right-to-left6 mm anterior-to-posterior) was performed to enable imaging of a portion of the right sensory and visual cortices. Cannulation of the femoral artery allowed the delivery of drugs, sampling of blood, and monitoring of blood pressure.

[0131] Cardiac Arrest (CA) and Cardiopulmonary Resuscitation (CPR):

[0132] FIG. 15 shows the multimodal setup employed in the experiments. At the onset of each experiment, the level of isoflurane was decreased from 2% to 0.5-1.0%. Concurrently, the mixture of inhaled gases was altered from 50% O.sub.2+50% N.sub.2 to 100% O.sub.2. Two minutes later, to reduce confounding effects of isoflurane on cerebral perfusion and metabolism, the anesthesia was turned off, at which time the subject breathed room air (21% O.sub.2). During this same period, 1 mL of 2 mg/kg Vecuronium (a neuromuscular blocker) and 1 mL of heparinized saline were administered intravenously, which led to respiration controlled solely by the ventilator. This stage of the experiment lasted for 3 min, after which the ventilator was turned off to induce asphyxia, leading to progressive hypoxic hypercarbic hypotension. CA was defined as the period over which the pulse pressure was below 10 mmHg and systolic pressure below 30 mmHg. The conditions of these experiments induced pulseless electrical activity, which is common in CA patients in a hospital setting.

[0133] Forty-five seconds before the end of the CA period, the ventilator was turned on (respiratory rate=85 breaths/min, PIP=17.5-18.5 cmH.sub.2O, PEEP=3 cmH.sub.2O at 2.5 LPM), and 100% oxygen was delivered. Immediately before the onset of CPR, 0.01 mg/kg epinephrine, 1 mmol/kg sodium bicarbonate, and 2 mL of heparinized saline were administered intravenously. Then, CPR was performed via external cardiac massage and terminated upon return of spontaneous circulation (ROSC), as identified from arterial blood pressure measurements. Subsequently, the animal was monitored continuously with arterial blood pressure, optical imaging, and ECoG for an additional 2 hr, after which the animal was euthanized with pentobarbital. Recovery of ECoG signal following ROSC was quantified by (1) time to initial resumption (burst) of ECoG activity, and (2) ECoG Information Quantity (IQ) 90 min post-ROSC.

[0134] Laser Speckle Imaging (LSI):

[0135] For LSI, an 809 nm laser with long coherence length served as the light source. To increase uniformity of illumination over the imaged region of interest (ROI), a ground-glass diffuser was placed between the laser and the brain. A CCD camera detected the backscattered light with a 10 ms exposure time, resulting in image acquisition at a frame rate of 60 Hz. Using a 55 sliding spatial window filter, the equation K=/<I> was employed to calculate the local speckle contrast K at each pixel, where <I> was the mean intensity within the filter and a the standard deviation within the filter. Then, the speckle flow index (SFI) was determined from the values of K and the exposure time T via a simplified speckle imaging equation SFI=1/(2TK.sup.2). Time-resolved SFI curves were generated by taking the mean of the SFI over a selected ROI at each time point.

[0136] Spatial Frequency Domain Imaging (SFDI):

[0137] For SFDI, light-emitting diodes (LEDs) of three different wavelengths (655 nm, 730 nm, 850 nm) were used as light sources. The light was directed to a spatial light modulator that projected square-wave patterns onto the brain. Backscattered light was captured using a scientific complementary metal-oxide semiconductor (sCMOS) camera. An Arduino Due microcontroller board was used to synchronize the camera acquisition, spatial light modulator, and LEDs. For each wavelength, four patterns were projected onto the tissue in sequence. The first pattern was non-modulated (i.e., DC illumination), and the three subsequent patterns were modulated at spatial frequency 0.3 mm.sup.1 with three distinct spatial phases to enable demodulation. Thus, there were a total of (3 wavelengths4 frames)=12 frames of SFDI data for each measurement time point. The detected square wave pattern could be approximated as a sinusoid, thus allowing demodulation. With this acquisition scheme, it was possible to reconstruct tissue hemodynamics and CMRO.sub.2, at an effective imaging rate of 14 Hz.

[0138] After demodulating the spatially-modulated data, the diffuse reflectance at each time point and wavelength was calculated from the raw data via calibration against a tissue-simulating phantom with known optical properties. The diffuse reflectance maps were then fit with a Monte Carlo model to extract the tissue absorption coefficient .sub.a and reduced scattering coefficient .sub.s at each wavelength. Next, the average .sub.s was determined for a selected ROI and a new .sub.a determined using diffuse reflectance with the non-modulated pattern and this average .sub.s. To calculate the concentrations of oxygenated and deoxygenated hemoglobin (ctHbO.sub.2 and ctHb, respectively) within the tissue, this new .sub.a() spectrum was fit with the model spectrum .sub.a()=2.303(ctHbO.sub.2.sub.HbO2+ctHb.sub.Hb), where .sub.HbO2 and .sub.Hb were the molar extinction coefficients of oxy- and deoxy-hemoglobin, respectively. The total tissue hemoglobin concentration (ctHb.sub.tot) was calculated by summing ctHb and ctHbO.sub.2. The tissue oxygen saturation was determined using the equation StO.sub.2=ctHbO.sub.2/(ctHbO.sub.2+ctHb).

[0139] Correction of Speckle Flow Index for Tissue Absorption and Scattering:

[0140] K, and hence SFI, depends on local optical properties. To correct the measured SFI for dynamic optical properties, the measured K values were converted to a characteristic flow speed (v.sub.c) by using the following equation:

[00001]$\begin{array}{cc}{K}^2=\frac{\left(\frac{2}{T}\right).Math.{}_0^T.Math..Math..Math.G_1^2\left(\right).Math.\left(1-/T\right).Math.d.Math..Math.}{G_1^2\left(=0\right)}& \left(1\right)\end{array}$

[0141] is a constant (typically set to 1) related to polarization and coherence properties of the LSI instrumentation. From the Siegert relationship, the intensity autocorrelation function G.sub.2() is related to G.sub.1(), which, in turn, is described by the correlation diffusion equation:

.sup.2G.sub.1().sub.effG.sub.1()=q(2)

[0142] In Eq. (2), q is the source term: and .sub.eff=(3.sub.a,dyn.sub.tr).sup.1/2, where .sub.tr=(.sub.a+.sub.s) is the tissue transport coefficient and .sub.a,dyn=(.sub.a+.sub.sk.sub.o.sup.2<r.sup.2()>/3) the dynamic tissue absorption coefficient. In the equation for .sub.a,dyn, <r.sup.2()> is the mean square displacement of the moving scatterers (i.e., the red blood cells) and k.sub.o is the photon wavenumber. Solving Eq. (2), G.sub.1() can be written as:

[00002]$\begin{array}{cc}{G}_1\left(\right)=\frac{3.Math.P_o.Math.A.Math..Math.\frac{{}_s^{}}{{}_{\mathrm{tr}}}}{\left(\frac{{}_{\mathrm{eff}}}{{}_{\mathrm{tr}}}+1\right).Math.\left(\frac{{}_{\mathrm{eff}}}{{}_{\mathrm{tr}}}+3.Math.A\right)}& \left(3\right)\end{array}$

[0143] In Eq. (3), P.sub.o is the incident optical power, and A is a function of the tissue refractive index. All other terms in Eq. (3) are exclusively functions of the static and dynamic tissue absorption and scattering coefficients (.sub.a, .sub.s, .sub.a,dyn). .sub.a,dyn is a function of <r.sup.2()>, and <r.sup.2()> is related to the characteristic flow speed v.sub.c via the equation <r.sup.2()>=v.sub.c.sup.2 (for directional flow). Using this framework and inputting the measured value of K from LSI and the measured .sub.a and .sub.s from SFDI at each time point, Eq. (1) was solved for v.sub.c at each time point and each pixel by iterating over a pre-defined grid of potential v.sub.c values and minimizing a least-squares cost function. The resulting spatiotemporal values of v.sub.c were used in place of SFI in the subsequent steps to achieve an optical property-corrected calculation of CMRO.sub.2.

[0144] Absolute Cerebral Metabolic Rate of Oxygen (CMRO.sub.2) Calculation:

[0145] To calculate absolute CMRO.sub.2, the following equation is used as a starting point:

CMRO.sub.2=(CBF)(OEF)([O.sub.2].sub.a)(4)

[0146] In Eq. (4), CBF is the cerebral blood flow, [O.sub.2].sub.a is the arterial concentration of oxygen, and OEF is the oxygen extraction fraction, equal to ([O.sub.2].sub.a[O.sub.2].sub.v/[O.sub.2].sub.a, where [O.sub.2].sub.v is the venous concentration of oxygen. For a single arteriole, (OEF)([O.sub.2].sub.a) represents the molar concentration of oxygen that was extracted from that arteriole and used by the brain for metabolic processes related to the synthesis of ATP. This quantity is equivalent to the molar concentration of deoxygenated hemoglobin that arrives in a nearby venule following oxygen extraction by the brain. Therefore, within our measurement paradigm, Eq. (4) is re-written as:

CMRO.sub.2=4(v.sub.c)(ctHb.sub.v)(Hb.sub.bl/<ctHb.sub.tot>.sub.p)(5)

[0147] In Eq. (5), ctHb.sub.v is the tissue concentration of deoxygenated hemoglobin in a region of interest atop a large vein in the ctHb maps obtained from SFDI. The factor of 4 accounts for the fact that the hemoglobin molecule has four binding sites for oxygen. Since v.sub.e is a characteristic flow parameter and not an absolute value of blood flow, it is necessary to include the proportionality constant in the equation to convert v.sub.c into a quantity with units of absolute flow speed.

[0148] The factor (Hb.sub.bl/<ctHb.sub.tot>.sub.p) accounts for partial-volume effects caused by the diffuse nature of light propagation in the brain. Eq. (4) requires an intra-vascular oxygen concentration, but SFDI measures a bulk tissue deoxyhemoglobin concentration. Hence, a blood-volume fraction term is required to convert between these two quantities. The numerator, Hb.sub.bl, is the concentration of hemoglobin in the blood sampled from the femoral artery of the animal during the arterial blood gas measurement (ABG). The denominator, <ctHb.sub.tot>.sub.p, is the mean total tissue hemoglobin concentration in the parenchyma during the period that the ABG was acquired. The factor (Hb.sub.bl/<ctHb.sub.tot>.sub.p) enables the required conversion of optical ctHb measurements from the scale of a tissue hemoglobin concentration to the scale of a vascular hemoglobin concentration, mitigating the partial volume effect and allowing us to measure CMRO.sub.2 on an absolute scale.

[0149] The parameter is typically unknown; thus, the quantity reported in optical brain imaging studies is usually the relative CMRO.sub.2 (rCMRO.sub.2). However, in this example, absolute CMRO.sub.2 was measured using a zero-flow boundary condition, which is provided by the onset of global cerebral ischemia in the animal model:

4(v.sub.c)|.sub.t->tasph(ctHb.sub.v)|.sub.t->tasph=4(dctHb.sub.v/dt)|.sub.t->tasph+(6)

[0150] This procedure was performed for each of the 10 subjects in this study, using the values of SFI and ctHb.sub.v during the period immediately before asphyxia (t->t.sub.asph.sup.) and the mean rate of change dctHb.sub.v/dt immediately after the onset of asphyxia (t->t.sub.asph.sup.+). The value of dctHb.sub.v/dt was measured by fitting a sigmoid function to the ctHb curve during the beginning of the zero-flow period, finding the t.sub.50 value of the sigmoid, linearizing the sigmoid within a 30 sec window centered on the t.sub.50 point, and calculating the slope of the resulting line segment. The values of and absolute CMRO.sub.2 then were calculated over the entire craniectomy region at each measurement time point.

[0151] Results:

[0152] Immediately following cardiac arrest, cerebral hemodynamics are spatially heterogeneous (FIGS. 16A-C). For example, the rate of change in cerebral ctHb following the start of asphyxia was 58.332.3 M/min in a ROI selected over a large vein (FIG. 16B), but only 33.613.6 M/min in a ROI selected over the parenchyma.

[0153] Maps of absolute CMRO.sub.2 throughout a representative CA/CPR experiment are shown in FIG. 17. At baseline, the animal is under anesthesia (2% isoflurane). After 2 min of anesthesia washout, the CMRO.sub.2 increased by a factor of 2 as the subject woke up. Following the onset of ischemia, the CMRO.sub.2 rapidly decreased as the subject entered CA. After resuscitation, the CMRO.sub.2 rapidly increased until reaching a maximum value at 8 min post CPR (during hyperemia). Subsequently, the CMRO.sub.2 decreased towards baseline as cerebral electrical activity resumed (12 min post-CPR).

[0154] Calculation of changes in CMRO.sub.2 are affected by optical properties (FIG. 18). Optical properties measured with SFDI, along with Eq. (1), enable calculation of a characteristic flow speed (v.sub.c) that can be used in place of SFI for the calculation of relative CMRO.sub.2 (rCMRO.sub.2). During each of the experimental phases, a comparison of rCMRO.sub.2 trends suggests that the metabolic activity is at times greater (i.e., during the hyperemic phase) and lower (i.e., during cardiac arrest) with use of SFI instead of v.sub.c in the calculation of CMRO.sub.2. This result demonstrates the effects that hyperdynamic changes in cerebral optical properties can have on calculations of CMRO.sub.2 dynamics using SFI alone to measure flow.

[0155] Measurements of baseline CMRO.sub.2 alone are associated with ECoG IQ (FIGS. 19A-C). For the 10 subjects included in this study, it was observed that higher absolute CMRO2 values at baseline (with 5 min and 7 min CA durations pooled) were associated with increased cerebral electrical activity (ECoG IQ) 90 min post-CPR (Spearman; r=0.63; p=0.076) (FIG. 19B). Further analysis was performed by dividing the subjects into High IQ and Low IQ groups with ECoG values greater than or less than 0.75, respectively. The baseline absolute CMRO.sub.2 was significantly different between the High IQ and Low IQ groups (FIG. 19C; p=0.016). These data collectively suggest that the absolute CMRO.sub.2 measurement at baseline may help to predict cerebral electrical recovery following CA and CPR.

[0156] FIG. 20 shows distributions of CMRO.sub.2 values measured with the imaging setup of the present example, as compared to values reported in the literature using various medical imaging approaches, including magnetic resonance methods (MRI/MRS) and positron emission tomography (PET). The LSI+SFDI method reported here measures absolute CMRO2 values that are within the range measured with these approaches, suggesting the accuracy of this optical imaging approach to determine CMRO.sub.2.

[0157] Discussion:

[0158] The present example is believed to be the first demonstration of dynamic imaging of absolute cerebral metabolic rate of oxygen (CMRO.sub.2) in the living brain using a combination of LSI and SFDI techniques. Tissue optical properties were measured with SFDI to account for their effects on interpretation of the LSI information. The zero-flow condition inherent in the CA experimental paradigm were then used to solve for the coefficient in the CMRO.sub.2 equation by using a continuity condition at the boundary between normal flow and zero-flow states. Using this technique, quantitative spatial mapping of absolute CMRO.sub.2 was performed continuously throughout the different stages of the CA+CPR experiment. The CMRO.sub.2 obtained from this optical system agreed well with established brain imaging techniques (PET, MRI/MRS). As an application, from ongoing research on optical and electrocerebral measurements during CA and CPR, there is a significant difference between baseline absolute CMRO.sub.2 values of rats with high ECoG IQ (>0.75) 90 min post-ROSC and rats with low ECoG IQ (<0.75) 90 min post-CPR, pooled across both 5 and 7 min CA durations.

[0159] This paradigm for measuring absolute CMRO.sub.2, in units of M O.sub.2/min, enables direct comparison of metabolic activity among subjects, across separate imaging sessions, and on different days for a single subject. This approach potentially enables longitudinal monitoring of cerebral recovery for days or weeks following ischemia and reperfusion. The methods of the present example may be applied to quantitative measurement of metabolic recovery and flow-metabolic coupling and uncoupling in preclinical models of ischemic conditions such as CA and stroke.

[0160] Optical Imaging Segments Venous Regions to Better Quantify Cerebral Oxygen Extraction:

[0161] The imaging capability of the device of the present example allows the segmentation of a ROI atop a prominent vein, which enables more accurate measurements of the quantity of deoxygenated venous blood and, hence, the quantity of oxygen consumed by the brain. With the use of a larger ROI, the local CMRO.sub.2 would be systematically underestimated due to inclusion of the parenchyma in the ROI, as oxygen extraction in the parenchyma is lower than in individual vessels. CMRO.sub.2 models of diffuse light transport implicitly assume that the concentration of deoxygenated hemoglobin is that within the veins specifically, and not the bulk tissue. However, most diffuse optics-based CMRO.sub.2 measurements are unable to satisfy this condition, as they typically use fiber-based spectroscopic techniques that sample the bulk tissue and thus cannot distinguish between venous and mixed arterial-venous parenchymal regions. In this example, the use of diffuse optical imaging allows the use of deoxyhemoglobin concentrations measured in a venous ROI to overcome this limitation and thus obtain more accurate quantitative values of CMRO.sub.2.

[0162] Correction of CMRO.sub.2 Data for Partial-Volume Effects:

[0163] Simply selecting a ROI that is coincident with a venule is not enough. To further refine the measurement of deoxygenated hemoglobin into an intravascular hemoglobin concentration, a partial-volume correction to the CMRO.sub.2 equation was employed. To accurately incorporate this scaling term, it is necessary to know the concentration of total hemoglobin (Hb.sub.bl) within the blood of each animal. In this example, these values were acquired via arterial blood gas (ABG) measurement before CA. A coefficient of variation of 13% in Hb.sub.bl was determined from the measurements. If the variation in Hb.sub.bl among the different subjects were not considered, an additional error of 12-25% in the measured CMRO.sub.2 would be achieved due to this within-group variability in Hb.sub.bl values.

[0164] Contributions of Directed Flow Versus Diffuse Flow:

[0165] Here, it was assumed that the corrected flow speed could be attributed solely to a directed-flow term (Eq. 2). Previous studies have used a Brownian diffusion term as the free parameter when fitting for flow speed or constrained the fit in a model system such that one could choose to fit for either diffuse or directed flow, but not both simultaneously. Some have used high-speed LSI to map the autocorrelation function pixel-by-pixel in the rodent brain, identifying the dominant type of particle motion at each pixel. They observed that the directed flow term was dominant in large vessels, while the diffuse flow term was dominant in the parenchyma.

[0166] In the present example, it was not possible to rigorously solve for the autocorrelation function because the sampling frequency of the LSI data acquisition was too low. Instead, a two-step approach was used, which consisted of (1) using SFDI data to account for the effects of optical properties on interpreting the LSI data, and (2) fitting the resulting corrected data to a model of directed flow to extract a characteristic flow speed. This method provided characteristic flow speeds that were similar to previously-reported values.

[0167] Limitations of Zero-Flow Condition:

[0168] One approach for measuring absolute CMRO.sub.2 requires temporary induction of a zero-flow condition in the brain. In this example, this condition was met by using a CA model in rats. However, there is a clear need for alternative approaches for interrogating absolute CMRO.sub.2 without creating harmful perturbations. Theses alternative approaches may incorporate techniques such as temporarily clamping the middle cerebral artery or administering sub-lethal doses of potassium chloride to temporarily induce a zero-flow condition that can be quickly reversed without long-term harm to the animal.

Conclusion:

[0169] This is believed to be the first example of absolute cerebral metabolic rate of oxygen (CMRO.sub.2) mapping in the rat brain using diffuse optical imaging. CMRO.sub.2 allows for quantitative assessment of cerebral metabolism without the need for baseline measurements, enabling longitudinal comparison between animals and among multiple days of measurement on an absolute scale. The CMRO2 measurements provided by this multimodal system were in good agreement with those previously measured in the brain of anesthetized rats using PET and MRI. This method shows significant potential for assessing and monitoring cerebral metabolism and predicting cerebral response to ischemic injury.

[0170] Although there has been shown and described the preferred embodiment of the present invention, it will be readily apparent to those skilled in the art that modifications may be made thereto which do not exceed the scope of the appended claims. Therefore, the scope of the invention is only to be limited by the following claims. In some embodiments, the figures presented in this patent application are drawn to scale, including the angles, ratios of dimensions, etc. In some embodiments, the figures are representative only and the claims are not limited by the dimensions of the figures. In some embodiments, descriptions of the inventions described herein using the phrase comprising includes embodiments that could be described as consisting essentially of or consisting of, and as such the written description requirement for claiming one or more embodiments of the present invention using the phrase consisting essentially of or consisting of is met.