Phantoms for diffuse optical imaging based on totally absorbing objects, part 2: experimental implementation

Criteria for the design of tissue-mimicking phantoms for the standardization of biophotonic instrumentation

A lack of accepted standards and standardized phantoms suitable for the technical validation of biophotonic instrumentation hinders the reliability and reproducibility of its experimental outputs. In this Perspective, we discuss general criteria for the design of tissue-mimicking biophotonic phantoms, and use these criteria and state-of-the-art developments to critically review the literature on phantom materials and on the fabrication of phantoms. By focusing on representative examples of standardization in diffuse optical imaging and spectroscopy, fluorescence-guided surgery and photoacoustic imaging, we identify unmet needs in the development of phantoms and a set of criteria (leveraging characterization, collaboration, communication and commitment) for the standardization of biophotonic instrumentation.

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.

Probe-hosted large area silicon photomultiplier and high-throughput timing electronics for enhanced performance time-domain functional near-infrared spectroscopy

Two main bottlenecks prevent time-domain diffuse optics instruments to reach their maximum performances, namely the limited light harvesting capability of the detection chain and the bounded data throughput of the timing electronics. In this work, for the first time to our knowledge, we overcome both those limitations using a probe-hosted large area silicon photomultiplier detector coupled to high-throughput timing electronics. The system performances were assessed based on international protocols for diffuse optical imagers showing better figures with respect to a state-of-the-art device. As a first step towards applications, proof-of-principle in-vivo brain activation measurements demonstrated superior signal-to-noise ratio as compared to current technologies.

Space-enhanced time-domain diffuse optics for determination of tissue optical properties in two-layered structures

A novel methodology for solving the inverse problem of diffuse optics for two-layered structures is proposed to retrieve the absolute quantities of optical absorption and reduced scattering coefficients of the layers simultaneously. A liquid phantom with various optical absorption properties in the deep layer is prepared and experimentally investigated using the space-enhanced time-domain method. Monte-Carlo simulations are applied to analyze the different measurements in time domain, space domain, and by the new methodology. The deviations of retrieved values from nominal values of both layers' optical properties are simultaneously reduced to a very low extent compared to the single-domain methods. The reliability and uncertainty of the retrieval performance are also considerably improved by the new methodology. It is observed in time-domain analyses that for the deep layer the retrieval of absorption coefficient is almost not affected by the scattering properties and this kind of "deep scattering neutrality" is investigated and overcome as well.

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.

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

Time-domain optical brain imaging techniques introduce a number of different measurands for analyzing absorption changes located deep in the tissue, complicated by superficial absorption changes and noise. We implement a method that allows analysis, quantitative comparison and performance ranking of measurands under various conditions - including different values of reduced scattering coefficient, thickness of the superficial layer, and source-detector separation. Liquid phantom measurements and Monte Carlo simulations were carried out in two-layered geometry to acquire distributions of times of flight of photons and to calculate the total photon count, mean time of flight, variance, photon counts in time windows and ratios of photon counts in different time windows. Quantitative comparison of performance was based on objective metrics: relative contrast, contrast-to-noise ratio (CNR) and depth selectivity. Moreover, the product of CNR and depth selectivity was used to rank the overall performance and to determine the optimal source-detector separation for each measurand. Variance ranks the highest under all considered conditions.

Time-Gated Single-Photon Detection in Time-Domain Diffuse Optics: A Review

This work reviews physical concepts, technologies and applications of time-domain diffuse optics based on time-gated single-photon detection. This particular photon detection strategy is of the utmost importance in the diffuse optics field as it unleashes the full power of the time-domain approach by maximizing performances in terms of contrast produced by a localized perturbation inside the scattering medium, signal-to-noise ratio, measurement time and dynamic range, penetration depth and spatial resolution. The review covers 15 years of theoretical studies, technological progresses, proof of concepts and design of laboratory systems based on time-gated single-photon detection with also few hints on other fields where the time-gated detection strategy produced and will produce further impact.

Implementation of the extended Kalman filter for determining the optical and geometrical properties of turbid layered media by time-resolved single distance measurements

In this article we propose an implementation of the extended Kalman filter (EKF) for the retrieval of optical and geometrical properties in two-layered turbid media assuming a dynamic setting, where absorption of each layer was changed in different steps. Prior works implemented the EKF in frequency-domain with several pairs of light sources and detectors and for static parameters estimation problems. Here we explore the use of the EKF in single distance, time-domain measurements, together with a corresponding forward model. Results show good agreement between retrieved and nominal values, with rather narrow analytical credibility intervals, indicating that the recovery process has low uncertainty, especially for the absorption coefficients.

Broadband Time Domain Diffuse Optical Reflectance Spectroscopy: A Review of Systems, Methods, and Applications

This review presents recent developments and a wide overview of broadband time domain diffuse optical spectroscopy (TD-DOS). Various topics including physics of photon migration, advanced instrumentation, methods of analysis, applications covering multiple domains (tissue chromophore, in vivo studies, food, wood, pharmaceutical industry) are elaborated. The key role of standardization and recent studies in that direction are discussed. Towards the end, a brief outlook is presented on the current status and future trends in broadband TD-DOS.

Cerebral oximetry performance testing with a 3D-printed vascular array phantom

Cerebral oximetry based on near-infrared spectroscopy represents a unique noninvasive tool for real-time surgical monitoring, yet studies have shown a significant discrepancy in accuracy among commercial systems. Towards the establishment of a standardized method for performance testing, we have studied a solid phantom approach - based on a 3D-printed cerebrovascular module (CVM) incorporating an array of 148 cylindrical channels - that has several advantages over liquid phantoms. Development and characterization of a CVM prototype are described, including high-resolution imaging and spectrophotometry measurements. The CVM was filled with whole bovine blood tuned over an oxygen saturation range of 30-90% and molded-silicone layers simulating extracerebral tissues were used to evaluate penetration depth. Saturation measurement accuracy was assessed in two commercially-available clinical cerebral oximeters. For one oximeter, both neonatal and pediatric sensors showed a high degree of precision, whereas accuracy was strongly dependent on saturation level and extracerebral geometry. The second oximeter showed worse precision, yet greater robustness to variations in extracerebral layers. These results indicate that 3D-printed channel array phantoms represent a promising new approach for standardized testing of clinical oximeters.

MAESTROS: A Multiwavelength Time-Domain NIRS System to Monitor Changes in Oxygenation and Oxidation State of Cytochrome-C-Oxidase

We present a multiwavelength, multichannel, time-domain near-infrared spectroscopy system named MAESTROS. This instrument can measure absorption and scattering coefficients and can quantify the concentrations of oxy- and deoxy-haemoglobin ([HbO₂], [HHb]), and oxidation state of cytochrome-c-oxidase ([oxCCO]). This system is composed of a supercontinuum laser source coupled with two acousto-optic tuneable filters. The light is collected by four photomultipliers tubes, connected to a router to redirect the signal to a single time-correlated single-photon counting card. The interface between the system and the tissue is based on optical fibres. This arrangement allows us to resolve up to 16 wavelengths, within the range of 650-900 nm, at a sampling rate compatible with the physiology (from 0.5 to 2 Hz). In this paper, we describe the system and assess its performance based on two specifically designed protocols for photon migration instruments, the basic instrument protocol and nEUROPt protocols, and on a well characterized liquid phantom based on Intralipid and water. Then, the ability to resolve [HbO₂ ], [HHb], and [oxCCO] is demonstrated on a homogeneous liquid phantom, based on blood for [HbO₂], [HHb], and yeast for [oxCCO]. In the future, the system could be used to monitor brain tissue physiology.

Instrumental, optical and geometrical parameters affecting time-gated diffuse optical measurements: a systematic study

In time-domain diffuse optics the sensitivity to localized absorption changes buried inside a diffusive medium depends strongly on the interplay between instrumental, optical and geometrical parameters, which can hinder the theoretical advantages of novel measurement strategies like the short source-detector distance approach. Here, we present a study based on experimental measurements and simulations to comprehensively evaluate the effect of all different parameters. Results are evaluated exploiting standardized figures of merit, like contrast and contrast-to-noise ratio, to quantify the system sensitivity to deep localized absorption perturbations. Key findings show that the most critical hardware parameter is the memory effect which ultimately limits the dynamic range. Further, a choice of the source-detector distance around 10 mm seems to be a good compromise to compensate non-idealities in practical systems still preserving the advantages of short distances. This work provides both indications for users about the best measurement conditions and strategies, and for technology developers to identify the most crucial hardware features in view of next generation diffuse optics systems.

Liquid phantoms for near-infrared and diffuse correlation spectroscopies with tunable optical and dynamic properties

We present the recipe and characterization for preparing liquid phantoms that are suitable for both near-infrared spectroscopy and diffuse correlation spectroscopy. The phantoms have well-defined and tunable optical and dynamic properties, and consist of a solution of water and glycerol with fat emulsion as the scattering element. The recipe takes into account the effect of bulk refractive index changes due to the addition of glycerol, which is commonly used to alter the sample viscosity.

Time Domain Near Infrared Spectroscopy Device for Monitoring Muscle Oxidative Metabolism: Custom Probe and In Vivo Applications

Measurement of muscle oxidative metabolism is of interest for monitoring the training status in athletes and the rehabilitation process in patients. Time domain near infrared spectroscopy (TD NIRS) is an optical technique that allows the non-invasive measurement of the hemodynamic parameters in muscular tissue: concentrations of oxy- and deoxy-hemoglobin, total hemoglobin content, and tissue oxygen saturation. In this paper, we present a novel TD NIRS medical device for muscle oxidative metabolism. A custom-printed 3D probe, able to host optical elements for signal acquisition from muscle, was develop for TD NIRS in vivo measurements. The system was widely characterized on solid phantoms and during in vivo protocols on healthy subjects. In particular, we tested the in vivo repeatability of the measurements to quantify the error that we can have by repositioning the probe. Furthermore, we considered a series of acquisitions on different muscles that were not yet previously performed with this custom probe: a venous-arterial cuff occlusion of the arm muscle, a cycling exercise, and an isometric contraction of the vastus lateralis.

Characterization of the Optical Properties of Turbid Media by Supervised Learning of Scattering Patterns

Fabricated tissue phantoms are instrumental in optical in-vitro investigations concerning cancer diagnosis, therapeutic applications, and drug efficacy tests. We present a simple non-invasive computational technique that, when coupled with experiments, has the potential for characterization of a wide range of biological tissues. The fundamental idea of our approach is to find a supervised learner that links the scattering pattern of a turbid sample to its thickness and scattering parameters. Once found, this supervised learner is employed in an inverse optimization problem for estimating the scattering parameters of a sample given its thickness and scattering pattern. Multi-response Gaussian processes are used for the supervised learning task and a simple setup is introduced to obtain the scattering pattern of a tissue sample. To increase the predictive power of the supervised learner, the scattering patterns are filtered, enriched by a regressor, and finally characterized with two parameters, namely, transmitted power and scaled Gaussian width. We computationally illustrate that our approach achieves errors of roughly 5% in predicting the scattering properties of many biological tissues. Our method has the potential to facilitate the characterization of tissues and fabrication of phantoms used for diagnostic and therapeutic purposes over a wide range of optical spectrum.

Chromophore decomposition in multispectral time-resolved diffuse optical tomography

Multicomponent phantom measurements are carried out to evaluate the ability of multispectral time domain diffuse optical tomography in reflectance geometry to quantify the position and the composition of small heterogeneities at depths of 1-1.5 cm in turbid media. Time-resolved data were analyzed with the Mellin-Laplace transform. Results show good localization and correct composition gradation of objects but still a lack of absolute material composition accuracy when no a priori geometry information is known.

Time-Resolved Diffuse Optical Spectroscopy and Imaging Using Solid-State Detectors: Characteristics, Present Status, and Research Challenges

Diffuse optical spectroscopy (DOS) and diffuse optical imaging (DOI) are emerging non-invasive imaging modalities that have wide spread potential applications in many fields, particularly for structural and functional imaging in medicine. In this article, we review time-resolved diffuse optical imaging (TR-DOI) systems using solid-state detectors with a special focus on Single-Photon Avalanche Diodes (SPADs) and Silicon Photomultipliers (SiPMs). These TR-DOI systems can be categorized into two types based on the operation mode of the detector (free-running or time-gated). For the TR-DOI prototypes, the physical concepts, main components, figures-of-merit of detectors, and evaluation parameters are described. The performance of TR-DOI prototypes is evaluated according to the parameters used in common protocols to test DOI systems particularly basic instrumental performance (BIP). In addition, the potential features of SPADs and SiPMs to improve TR-DOI systems and expand their applications in the foreseeable future are discussed. Lastly, research challenges and future developments for TR-DOI are discussed for each component in the prototype separately and also for the entire system.

Depth sensitivity of frequency domain optical measurements in diffusive media

The depth sensitivity functions for AC amplitude, phase (PH) and DC intensity signals have been obtained in the frequency domain (where the source amplitude is modulated at radio-frequencies) by making use of analytical solutions of the photon diffusion equation in an infinite slab geometry. Furthermore, solutions for the relative contrast of AC, PH and DC signals when a totally absorbing plane is placed at a fixed depth of the slab have also been obtained. The solutions have been validated by comparisons with gold standard Monte Carlo simulations. The obtained results show that the AC signal, for modulation frequencies < 200 MHz, has a depth sensitivity with similar characteristics to that of the continuous-wave (CW) domain (source modulation frequency of zero). Thus, the depth probed by such a signal can be estimated by using the formula of penetration depth for the CW domain (Sci. Rep.6, 27057 (2016)). However, the PH signal has a different behavior compared to the CW domain, showing a larger depth sensitivity at shallow depths and a less steep relative contrast as a function of depth. These results mark a clear difference in term of depth sensitivity between AC and PH signals, and highlight the complexity of the estimation of the actual depth probed in tissue spectroscopy.

Time-domain diffuse optical tomography using silicon photomultipliers: feasibility study

Silicon photomultipliers (SiPMs) have been very recently introduced as the most promising detectors in the field of diffuse optics, in particular due to the inherent low cost and large active area. We also demonstrate the suitability of SiPMs for time-domain diffuse optical tomography (DOT). The study is based on both simulations and experimental measurements. Results clearly show excellent performances in terms of spatial localization of an absorbing perturbation, thus opening the way to the use of SiPMs for DOT, with the possibility to conceive a new generation of low-cost and reliable multichannel tomographic systems.

Characterization of a time-resolved non-contact scanning diffuse optical imaging system exploiting fast-gated single-photon avalanche diode detection

We present a system for non-contact time-resolved diffuse reflectance imaging, based on small source-detector distance and high dynamic range measurements utilizing a fast-gated single-photon avalanche diode. The system is suitable for imaging of diffusive media without any contact with the sample and with a spatial resolution of about 1 cm at 1 cm depth. In order to objectively assess its performances, we adopted two standardized protocols developed for time-domain brain imagers. The related tests included the recording of the instrument response function of the setup and the responsivity of its detection system. Moreover, by using liquid turbid phantoms with absorbing inclusions, depth-dependent contrast and contrast-to-noise ratio as well as lateral spatial resolution were measured. To illustrate the potentialities of the novel approach, the characteristics of the non-contact system are discussed and compared to those of a fiber-based brain imager.

Mechanically switchable solid inhomogeneous phantom for performance tests in diffuse imaging and spectroscopy

A mechanically switchable solid inhomogeneous phantom simulating localized absorption changes was developed and characterized. The homogeneous host phantom was made of epoxy resin with black toner and titanium dioxide particles added as absorbing and scattering components, respectively. A cylindrical rod, movable along a hole in the block and made of the same material, has a black polyvinyl chloride cylinder embedded in its center. By varying the volume and position of the black inclusion, absorption perturbations can be generated over a large range of magnitudes. The phantom has been characterized by various time-domain diffuse optics instruments in terms of absorption and scattering spectra, transmittance images, and reflectance contrast. Addressing a major application of the phantom for performance characterization for functional near-infrared spectroscopy of the brain, the contrast was measured in reflectance mode while black cylinders of volumes from ≈20  mm3 to ≈270  mm3 were moved in lateral and depth directions, respectively. The new type of solid inhomogeneous phantom is expected to become a useful tool for routine quality check of clinical instruments or implementation of industrial standards provided an experimental characterization of the phantom is performed in advance.

Optimal estimation reconstruction of the optical properties of a two-layered tissue phantom from time-resolved single-distance measurements

In this work, we have tested the optimal estimation (OE) algorithm for the reconstruction of the optical properties of a two-layered liquid tissue phantom from time-resolved single-distance measurements. The OE allows a priori information, in particular on the range of variation of fit parameters, to be included. The purpose of the present investigations was to compare the performance of OE with the Levenberg–Marquardt method for a geometry and real experimental conditions typically used to reconstruct the optical properties of biological tissues such as muscle and brain. The absorption coefficient of the layers was varied in a range of values typical for biological tissues. The reconstructions performed demonstrate the substantial improvements achievable with the OE provided a priori information is available. We note the extreme reliability, robustness, and accuracy of the retrieved absorption coefficient of the second layer obtained with the OE that was found for up to six fit parameters, with an error in the retrieved values of less than 10%. A priori information on fit parameters and fixed forward model parameters clearly improves robustness and accuracy of the inversion procedure.

Spatial resolution in depth for time-resolved diffuse optical tomography using short source-detector separations

Diffuse optical tomography for medical applications can require probes with small dimensions involving short source-detector separations. Even though this configuration is seen at first as a constraint due to the challenge of depth sensitivity, we show here that it can potentially be an asset for spatial resolution in depth. By comparing two fiber optic probes on a test object, we first show with simulations that short source-detector separations improve the spatial resolution down to a limit depth. We then confirm these results in an experimental study with a state-of-the-art setup involving a fast-gated single-photon avalanche diode allowing maximum depth sensitivity. We conclude that short source-detector separations are an option to consider for the design of probes so as to improve image quality for diffuse optical tomography in reflectance.

Performance assessment of time-domain optical brain imagers, part 2: nEUROPt protocol

The nEUROPt protocol is one of two new protocols developed within the European project nEUROPt to characterize the performances of time-domain systems for optical imaging of the brain. It was applied in joint measurement campaigns to compare the various instruments and to assess the impact of technical improvements. This protocol addresses the characteristic of optical brain imaging to detect, localize, and quantify absorption changes in the brain. It was implemented with two types of inhomogeneous liquid phantoms based on Intralipid and India ink with well-defined optical properties. First, small black inclusions were used to mimic localized changes of the absorption coefficient. The position of the inclusions was varied in depth and lateral direction to investigate contrast and spatial resolution. Second, two-layered liquid phantoms with variable absorption coefficients were employed to study the quantification of layer-wide changes and, in particular, to determine depth selectivity, i.e., the ratio of sensitivities for deep and superficial absorption changes. We introduce the tests of the nEUROPt protocol and present examples of results obtained with different instruments and methods of data analysis. This protocol could be a useful step toward performance tests for future standards in diffuse optical imaging.