IDNE 701 ยท Week 6 Lecture

Optical Neuroimaging & fNIRS

Diffuse Optical Tomography, Hemodynamic Inversion & Bedside Neuromonitoring

๐Ÿ“„ Lecture Notes ๐Ÿ”ฌ Demo Lab ๐Ÿ“… Course Syllabus

The Near-Infrared Diagnostic Window

Why NIR light (650โ€“900 nm) uniquely penetrates deep into human brain tissue.

ฮป < 650 nm Electronic Absorption

Visible light is heavily absorbed by hemoglobin, melanin, and cytochromes. Penetration depth is limited to < 1 mm (skin surface only).

650โ€“900 nm Diagnostic Window

Absorption drops by orders of magnitude ($\mu_a \sim 0.01 - 0.1\,\text{mm}^{-1}$). Scattering dominates, allowing photons to diffuse several centimeters into tissue.

ฮป > 950 nm Water Absorption

Vibrational overtones of water ($H_2O$) and lipids cause strong infrared absorption, closing the deep penetration window.

760 nm
HbR Absorption Peak
805 nm
Isosbestic Point (Total Hb)
850 nm
HbO Dominant Absorption
1.40
Tissue Refractive Index (n)

Photon Transport in Turbid Media

Tissue is multiple-scattering: transport mean free path $l_{\text{tr}} = 1/\mu_s' \approx 1\,\text{mm}$.

$$\frac{1}{v} \frac{\partial \Phi(\vec{r}, t)}{\partial t} - \nabla \cdot \left[ D(\vec{r}) \nabla \Phi(\vec{r}, t) \right] + \mu_a(\vec{r}) \Phi(\vec{r}, t) = S(\vec{r}, t)$$
Diffusion Equation for Fluence Rate ฮฆ(r,t) where D โ‰ˆ 1 / (3ยทฮผs')

๐ŸŒ Banana-Shaped Adjoint Profile

Photons entering at the source diffuse outwards and upwards to the detector, forming a curved banana-shaped sensitivity kernel.

Mean Penetration Depth: $\bar{z}_{\text{max}} \approx \frac{1}{2} d_{\text{SD}}$ (at $d_{\text{SD}} = 30\,\text{mm}$, depth $\approx 15\,\text{mm}$).

๐Ÿ“ Scattering Anisotropy ($g \approx 0.9$)

Mitochondria and cell membranes cause forward-peaked Mie scattering. Reduced scattering coefficient:

$$\mu_s' = \mu_s (1 - g) \approx 1.0\,\text{mm}^{-1}$$

The Modified Beer-Lambert Law (MBLL)

Accounting for multiple scattering via Differential Pathlength Factor (DPF).

$$\Delta \text{OD}(\lambda) = -\ln\left(\frac{I(t)}{I_0}\right) = \left[ \varepsilon_{\text{HbO}}(\lambda) \cdot \Delta[\text{HbO}] + \varepsilon_{\text{HbR}}(\lambda) \cdot \Delta[\text{HbR}] \right] \cdot d_{\text{SD}} \cdot \text{DPF}(\lambda)$$
Effective Optical Pathlength L_opt = d_SD ยท DPF (Adult DPF โ‰ˆ 6.0; Neonatal DPF โ‰ˆ 4.6)

Spectroscopic Inversion Matrix

Evaluating at $\lambda_1 = 760\,\text{nm}$ and $\lambda_2 = 850\,\text{nm}$ allows exact linear separation:

$$\begin{bmatrix} \Delta[\text{HbO}] \\ \Delta[\text{HbR}] \end{bmatrix} = \frac{1}{d_{\text{SD}}} \mathbf{E}^{-1} \begin{bmatrix} \Delta\text{OD}_{760} / \text{DPF}_{760} \\ \Delta\text{OD}_{850} / \text{DPF}_{850} \end{bmatrix}$$

Total Hemoglobin & Blood Volume

Summing both chromophores gives total hemoglobin change:

$$\Delta[\text{HbT}] = \Delta[\text{HbO}] + \Delta[\text{HbR}] \propto \Delta\text{CBV}$$

Serves as direct surrogate for regional cerebral blood volume changes.

Optical Instrumentation: CW vs FD vs TD

CW Continuous-Wave

Mechanism: Steady-state DC light injection.

Measures: Relative concentration changes ($\Delta\text{HbO}, \Delta\text{HbR}$).

Pros/Cons: Extremely low cost, lightweight portable caps, robust. Cannot measure absolute absorption ($\mu_a$).

FD Frequency-Domain

Mechanism: Intensity-modulated RF lasers ($50 - 200\,\text{MHz}$).

Measures: Amplitude attenuation $A$ and phase shift $\theta$.

Pros/Cons: Uncouples $\mu_a$ and $\mu_s'$ directly. Yields absolute baseline $\text{StO}_2$. Moderate hardware cost.

TD Time-Domain (TCSPC)

Mechanism: Picosecond laser pulses + SPAD detectors.

Measures: Photon time-of-flight distribution (TPSF).

Pros/Cons: Gold-standard depth discrimination via late photon time-gating. Bulky and expensive.

Grand Challenge 06

Neonatal & Pediatric Bedside Neuromonitoring

๐Ÿ‘ถ Hypoxic-Ischemic Encephalopathy

1โ€“3 per 1,000 live births. Standard care: 72 hours of whole-body Therapeutic Hypothermia ($33.5^\circ\text{C}$).

Continuous bedside fNIRS assesses cerebral metabolic extraction and autoregulation without MRI transport risks.

๐Ÿฉธ FTOE Biomarker

Fractional Tissue Oxygen Extraction:

$$\text{FTOE} = \frac{\text{SpO}_2 - \text{rSO}_2}{\text{SpO}_2}$$

Elevated FTOE (> 0.38) flags severe ischemia/hypoperfusion; collapsed FTOE (< 0.15) signals mitochondrial cell death.

โšก Autoregulation Index (COx)

Pearson correlation between Mean Arterial Pressure (MAP) and cerebral $\text{rSO}_2$.

$\text{COx} \approx 0 \implies$ Intact autoregulation.
$\text{COx} > 0.3 \implies$ Pressure-passive circulation (high risk of IVH).

Neonatal Cranial Optics vs Adults

๐Ÿ‘ถ Neonatal Head Characteristics

Scalp + Skull Thickness: $1.5 - 2.5\,\text{mm}$ (vs adult $8 - 14\,\text{mm}$).

Reduced Extracerebral Shunting: Minimal scalp muscle/bone allows photons to enter cortex almost directly.

Optode Separation: $d_{\text{SD}} = 15 - 20\,\text{mm}$ gives deep cortical sensitivity.

DPF: $\approx 4.4 - 4.8$ (calibrated for thin unmyelinated head).

โš ๏ธ NICU Physical Constraints

Skin Fragility: Stratum corneum is immature; interface pressure must remain < 20 mmHg to prevent necrosis.

Phototherapy Interference: High-power blue phototherapy lamps (450 nm) for hyperbilirubinemia require > 60 dB optical filtering.

Incubator Heating: Optode surface must not exceed 41ยฐC (IEC 60601-1 standard).

Signal Processing & Artifact Removal

Isolating true cortical signals from physiological noise and motion.

๐ŸŒŠ Motion Artifact Correction

Head movement causes shear displacement and optode decoupling.

Wavelet Filtering: Strips high-kurtosis detail coefficients.
Spline Interpolation: Reconstructs baseline offsets seamlessly.

๐Ÿ“ Short-Separation Regression

Short channels ($d_{\text{short}} \approx 8\,\text{mm}$) sample scalp/skin exclusively.

Adaptive GLM subtraction eliminates systemic extracerebral cardiac and vasomotor Mayer waves ($0.1\,\text{Hz}$).

Bandpass Filtering

Passband: $0.01 - 0.2\,\text{Hz}$

โ€ข Rejects cardiac pulsation ($1 - 2\,\text{Hz}$ neonate)
โ€ข Rejects respiratory rhythm ($0.3 - 0.7\,\text{Hz}$)
โ€ข Preserves slow neurovascular BOLD coupling

Today's Demo: Worksheet 6.1 Lab

๐Ÿ”ฌ Real-Time Synthetic Pipeline

1. Simulate Raw Photons: Dual-wavelength attenuation curves at $760\,\text{nm}$ and $850\,\text{nm}$.

2. MBLL Transformation: Solve the $2\times 2$ extinction matrix for $\Delta[\text{HbO}]$ and $\Delta[\text{HbR}]$.

3. Calculate Bedside Biomarkers: Compute regional saturation $\text{rSO}_2$ and extraction $\text{FTOE}$.

๐ŸŽฏ Student Objectives

โ€ข Evaluate HIE clinical scenarios (therapeutic hypothermia vs rewarming).

โ€ข Assess how fontanelle proximity alters required optical source-detector spacing.

โ€ข Prepare system block diagrams for next week's Mid-Semester Design Review.

๐Ÿš€ Launch Worksheet 6.1 Demo Lab โ†’

Looking Ahead to Week 7

Next Milestone: Mid-Semester Concept Pin-Up & Preliminary Design Review

๐Ÿ”ฌ Complete Worksheet 6.1 ๐Ÿ“„ Read Lecture Notes ๐Ÿ“… Week 7 Milestone Details
IDNE 701: Introduction to Neuroengineering ยท Department of Biomedical Engineering