Performance of measurands in time-domain optical brain imaging: depth selectivity versus contrast-to-noise ratio

Quantitative evaluation of the impact of variation of optical parameters on the estimation of blood hematocrit and oxygen saturation for dual-wavelength photoacoustics

Photoacoustic (PA) spectroscopy is considered to be one of the most effective ways to measure the levels of hematocrit (H) and oxygenation saturation (S O 2) of blood, which are essential for diagnosing blood-related illnesses. This simulation study aims to investigate the impact of individual optical parameters, i.e., optical absorption coefficient (μ a), scattering coefficient (μ s), and anisotropy factor (g), on the accuracy of this technique in estimating the blood properties. We first performed the Monte Carlo simulations, using realistic optical parameters, to obtain the fluence maps for various samples. The wavelengths of the incident light were chosen to be 532, 700, 1000, and 1064 nm. Thereafter, the k-Wave simulations were executed, incorporating those fluence maps to generate the PA signals. The blood properties were obtained using the PA signals. We introduced variations in μ a, μ s, and g ranging from -10% to +10%, -10% to +10%, and -5% to +1%, respectively, at 700 and 1000 nm wavelengths. One parameter, at both wavelengths, was changed at a time, keeping others fixed. Subsequently, we examined how accurately the blood parameters could be determined at physiological hematocrit levels. A 10% variation in μ a induces a 10% change in H estimation but no change in S O 2 determination. Almost no change has been seen for μ s variation. However, a 5% (-5% to 0%) variation in the g factor resulted in approximately 160% and 115% changes in the PA signal amplitudes at 700 and 1000 nm, respectively, leading to ≈125% error in hematocrit estimation and ≈14% deviation in S O 2 assessment when nominal S O 2=70%. It is clear from this study that the scattering anisotropy factor is a very sensitive parameter and a small change in its value can result in large errors in the PA estimation of blood properties. In the future, in vitro experiments with pathological blood (inducing variation in the g parameter) will be performed, and accordingly, the accuracy of the PA technique in quantifying blood H and S O 2 will be evaluated.

Depth-selective method for time-domain diffuse reflectance measurements: validation study of the dual subtraction technique

Research on the spatial distribution of sensitivity of time-domain near infrared diffuse reflectance measurement is reported in this paper. The main objective of the investigation is to validate theoretically calculated sensitivity profiles for a measurement geometry with two detectors and two sources in which sensitivity profiles of statistical moments of distributions of time of flight of photons (DTOFs) are spatially restricted to a region underneath the detectors. For this dual subtraction method, smaller sensitivities to changes appearing in the superficial layer of the medium were observed compared to the single distance and single subtraction methods. Experimental validation of this approach is based on evaluation of changes in the statistical moments of DTOFs measured on a liquid phantom with local absorption perturbations. The spatial distributions of sensitivities, depth-related sensitivity and depth selectivities were obtained from the dual subtraction method and compared with those from single distance and single subtraction approaches. Also, the contrast to noise ratio (CNR) was calculated for the dual subtraction technique and combined with depth selectivity in order to assess the overall performance (product of CNR and depth selectivity) of the method. Spatial sensitivity profiles from phantom experiments are in a good agreement with the results of theoretical studies and feature more locally restricted sensitivity volume with the point of maximal sensitivity located deeper. The highest value of overall performance was obtained experimentally for the second statistical moment in the dual subtraction method (∼10.8) surpassing that of the single distance method (∼8.7). This confirms the advantage of dual subtraction measurement geometries in the suppression of optical signals originated in the superficial layer of the medium.

Two-layered blood-lipid phantom and method to determine absorption and oxygenation employing changes in moments of DTOFs

Near-infrared spectroscopy (NIRS) is an established technique for measuring tissue oxygen saturation (StO₂), which is of high clinical value. For tissues that have layered structures, it is challenging but clinically relevant to obtain StO₂ of the different layers, e.g. brain and scalp. For this aim, we present a new method of data analysis for time-domain NIRS (TD-NIRS) and a new two-layered blood-lipid phantom. The new analysis method enables accurate determination of even large changes of the absorption coefficient (Δµ_a) in multiple layers. By adding Δµ_a to the baseline µ_a, this method provides absolute µ_a and hence StO₂ in multiple layers. The method utilizes (i) changes in statistical moments of the distributions of times of flight of photons (DTOFs), (ii) an analytical solution of the diffusion equation for an N-layered medium, (iii) and the Levenberg-Marquardt algorithm (LMA) to determine Δµ_a in multiple layers from the changes in moments. The method is suitable for NIRS tissue oximetry (relying on µ_a) as well as functional NIRS (fNIRS) applications (relying on Δµ_a). Experiments were conducted on a new phantom, which enabled us to simulate dynamic StO₂ changes in two layers for the first time. Two separate compartments, which mimic superficial and deep layers, hold blood-lipid mixtures that can be deoxygenated (using yeast) and oxygenated (by bubbling oxygen) independently. Simultaneous NIRS measurements can be performed on the two-layered medium (variable superficial layer thickness, L), the deep (homogeneous), and/or the superficial (homogeneous). In two experiments involving ink, we increased the nominal µ_a in one of two compartments from 0.05 to 0.25 cm^-1, L set to 14.5 mm. In three experiments involving blood (L set to 12, 15, or 17 mm), we used a protocol consisting of six deoxygenation cycles. A state-of-the-art multi-wavelength TD-NIRS system measured simultaneously on the two-layered medium, as well as on the deep compartment for a reference. The new method accurately determined µ_a (and hence StO₂) in both compartments. The method is a significant progress in overcoming the contamination from the superficial layer, which is beneficial for NIRS and fNIRS applications, and may improve the determination of StO₂ in the brain from measurements on the head. The advanced phantom may assist in the ongoing effort towards more realistic standardized performance tests in NIRS tissue oximetry. Data and MATLAB codes used in this study were made publicly available.

Lock-in functional near-infrared spectroscopy for measurement of the haemodynamic brain response

Here we show a method of the lock-in amplifying near-infrared signals originating within a human brain. It implies using two 90-degree rotated source-detector pairs fixed on a head surface. Both pairs have a joint sensitivity region located towards the brain. A direct application of the lock-in technique on both signals results in amplifying common frequency components, e.g. related to brain cortex stimulation and attenuating the rest, including all components not related to the stimulation: e.g. pulse, instrumental and biological noise or movement artefacts. This is a self-driven method as no prior assumptions are needed and the noise model is provided by the interfering signals themselves. We show the theory (classical modified Beer-Lambert law and diffuse optical tomography approaches), the algorithm implementation and tests on a finite element mathematical model and in-vivo on healthy volunteers during visual cortex stimulation. The proposed hardware and algorithm complexity suit the entire spectrum of (continuous wave, frequency domain, time-resolved) near-infrared spectroscopy systems featuring real-time, direct, robust and low-noise brain activity registration tool. As such, this can be of special interest in optical brain computer interfaces and high reliability/stability monitors of tissue oxygenation.

Time-domain NIRS system based on supercontinuum light source and multi-wavelength detection: validation for tissue oxygenation studies

We present and validate a multi-wavelength time-domain near-infrared spectroscopy (TD-NIRS) system that avoids switching wavelengths and instead exploits the full capability of a supercontinuum light source by emitting and acquiring signals for the whole chosen range of wavelengths. The system was designed for muscle and brain oxygenation monitoring in a clinical environment. A pulsed supercontinuum laser emits broadband light and each of two detection modules acquires the distributions of times of flight of photons (DTOFs) for 16 spectral channels (used width 12.5 nm / channel), providing a total of 32 DTOFs at up to 3 Hz. Two emitting fibers and two detection fiber bundles allow simultaneous measurements at two positions on the tissue or at two source-detector separations. Three established protocols (BIP, MEDPHOT, and nEUROPt) were used to quantitatively assess the system's performance, including linearity, coupling, accuracy, and depth sensitivity. Measurements were performed on 32 homogeneous phantoms and two inhomogeneous phantoms (solid and liquid). Furthermore, measurements on two blood-lipid phantoms with a varied amount of blood and Intralipid provide the strongest validation for accurate tissue oximetry. The retrieved hemoglobin concentrations and oxygen saturation match well with the reference values that were obtained using a commercially available NIRS system (OxiplexTS) and a blood gas analyzer (ABL90 FLEX), except a discrepancy occurs for the lowest amount of Intralipid. In-vivo measurements on the forearm of three healthy volunteers during arterial (250 mmHg) and venous (60 mmHg) cuff occlusions provide an example of tissue monitoring during the expected hemodynamic changes that follow previously well-described physiologies. All results, including quantitative parameters, can be compared to other systems that report similar tests. Overall, the presented TD-NIRS system has an exemplary performance evaluated with state-of-the-art performance assessment methods.

The Use of Supercontinuum Laser Sources in Biomedical Diffuse Optics: Unlocking the Power of Multispectral Imaging

Optical techniques based on diffuse optics have been around for decades now and are making their way into the day-to-day medical applications. Even though the physics foundations of these techniques have been known for many years, practical implementation of these technique were hindered by technological limitations, mainly from the light sources and/or detection electronics. In the past 20 years, the developments of supercontinuum laser (SCL) enabled to unlock some of these limitations, enabling the development of system and methodologies relevant for medical use, notably in terms of spectral monitoring. In this review, we focus on the use of SCL in biomedical diffuse optics, from instrumentation and methods developments to their use for medical applications. A total of 95 publications were identified, from 1993 to 2021. We discuss the advantages of the SCL to cover a large spectral bandwidth with a high spectral power and fast switching against the disadvantages of cost, bulkiness, and long warm up times. Finally, we summarize the utility of using such light sources in the development and application of diffuse optics in biomedical sciences and clinical applications.

Depth-selective data analysis for time-domain fNIRS: moments vs. time windows

Time-domain measurements facilitate the elimination of the influence of extracerebral, systemic effects, a key problem in functional near-infrared spectroscopy (fNIRS) of the adult human brain. The analysis of measured time-of-flight distributions of photons often relies on moments or time windows. However, a systematic and quantitative characterization of the performance of these measurands is still lacking. Based on perturbation simulations for small localized absorption changes, we compared spatial sensitivity profiles and depth selectivity for moments (integral, mean time of flight and variance), photon counts in time windows and their ratios for different time windows. The influence of the instrument response function (IRF) was investigated for all measurands and for various source-detector separations. Variance exhibits the highest depth selectivity among the moments. Ratios of photon counts in different late time windows can achieve even higher selectivity. An advantage of moments is their robustness against the shape of the IRF and instrumental drifts.