Liquid phantom for investigating light propagation through layered diffusive media

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.

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.

A Compact Multi-Distance DCS and Time Domain NIRS Hybrid System for Hemodynamic and Metabolic Measurements

In this work, we present a new multi-distance diffuse correlation spectroscopy (DCS) device integrated with a compact state-of-the-art time domain near infrared spectroscopy (TD-NIRS) device. The hybrid DCS and TD-NIRS system allows to retrieve information on blood flow, tissue oxygenation, and oxygen metabolic rate. The DCS device performances were estimated in terms of stability, repeatability, ability in retrieving variations of diffusion coefficient, influence of the tissue optical properties, effect of varying count rates and depth sensitivity. Crosstalk between DCS and TD-NIRS optical signals was also evaluated. Finally, in vivo experiments (venous and arterial cuff occlusions on the arm) were conducted to test the ability of the hybrid system in measuring blood flow variations.

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.

Systematic study of the effect of ultrasound gel on the performances of time-domain diffuse optics and diffuse correlation spectroscopy

Recently, multimodal imaging has gained an increasing interest in medical applications thanks to the inherent combination of strengths of the different techniques. For example, diffuse optics is used to probe both the composition and the microstructure of highly diffusive media down to a depth of few centimeters, but its spatial resolution is intrinsically low. On the other hand, ultrasound imaging exhibits the higher spatial resolution of morphological imaging, but without providing solid constitutional information. Thus, the combination of diffuse optical imaging and ultrasound may improve the effectiveness of medical examinations, e.g. for screening or diagnosis of tumors. However, the presence of an ultrasound coupling gel between probe and tissue can impair diffuse optical measurements like diffuse optical spectroscopy and diffuse correlation spectroscopy, since it may provide a direct path for photons between source and detector. A systematic study on the effect of different ultrasound coupling fluids was performed on tissue-mimicking phantoms, confirming that a water-clear gel can produce detrimental effects on optical measurements when recovering absorption/reduced scattering coefficients from time-domain spectroscopy acquisitions as well as particle Brownian diffusion coefficient from diffuse correlation spectroscopy ones. On the other hand, we show the suitability for optical measurements of other types of diffusive fluids, also compatible with ultrasound imaging.

Heuristic model for ballistic photon detection in collimated transmittance measurements

An heuristic model for ballistic photon detection in continuous-wave measurements of collimated transmittance through a slab is presented. The model is based on the small angle approximation and the diffusion equation and covers all the ranges of optical thicknesses of the slab from the ballistic to the diffusive regime. The performances of the model have been studied by means of comparisons with the results of gold standard Monte Carlo simulations for a wide range of optical thicknesses and two types of scattering functions. For a non-absorbing slab and field of view of the receiver less than 3° the model shows errors less than 15% for any value of the optical thickness. Even for an albedo value of 0.9, and field of view of the receiver less than 3° the model shows errors less than 20%. These results have been verified for a large set of scattering functions based on the Henyey-Greenstein model and Mie theory for spherical scatterers. The latter has also been used to simulate the scattering function of Intralipid, a diffusive material widely used as reference standard for tissue simulating phantoms. The proposed model represents an effective improvement compared to the existing literature.

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.

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.

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.

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

We present the experimental implementation and validation of a phantom for diffuse optical imaging based on totally absorbing objects for which, in the previous paper [J. Biomed. Opt.18(6), 066014, (2013)], we have provided the basic theory. Totally absorbing objects have been manufactured as black polyvinyl chloride (PVC) cylinders and the phantom is a water dilution of intralipid-20% as the diffusive medium and India ink as the absorber, filled into a black scattering cell made of PVC. By means of time-domain measurements and of Monte Carlo simulations, we have shown the reliability, the accuracy, and the robustness of such a phantom in mimicking typical absorbing perturbations of diffuse optical imaging. In particular, we show that such a phantom can be used to generate any absorption perturbation by changing the volume and position of the totally absorbing inclusion.

Multi-channel medical device for time domain functional near infrared spectroscopy based on wavelength space multiplexing

We have designed a compact dual wavelength (687 nm, 826 nm) multi-channel (16 sources, 8 detectors) medical device for muscle and brain imaging based on time domain functional near infrared spectroscopy. The system employs the wavelength space multiplexing approach to reduce wavelength cross-talk and increase signal-to-noise ratio. System performances have been tested on homogeneous and heterogeneous tissue phantoms following specifically designed protocols for photon migration instruments. Preliminary in vivo measurements have been performed to validate the instrument capability to monitor hemodynamic parameters changes in the arm muscle during arterial occlusion and in the adult head during a motor task experiment.

Optical characterization of two-layered turbid media for non-invasive, absolute oximetry in cerebral and extracerebral tissue

We introduce a multi-distance, frequency-domain, near-infrared spectroscopy (NIRS) method to measure the optical coefficients of two-layered media and the thickness of the top layer from diffuse reflectance measurements. This method features a direct solution based on diffusion theory and an inversion procedure based on the Levenberg-Marquardt algorithm. We have validated our method through Monte Carlo simulations, experiments on tissue-like phantoms, and measurements on the forehead of three human subjects. The Monte Carlo simulations and phantom measurements have shown that, in ideal two-layered samples, our method accurately recovers the top layer thickness (L), the absorption coefficient (µ a ) and the reduced scattering coefficient (µ' s ) of both layers with deviations that are typically less than 10% for all parameters. Our method is aimed at absolute measurements of hemoglobin concentration and saturation in cerebral and extracerebral tissue of adult human subjects, where the top layer (layer 1) represents extracerebral tissue (scalp, skull, dura mater, subarachnoid space, etc.) and the bottom layer (layer 2) represents cerebral tissue. Human subject measurements have shown a significantly greater total hemoglobin concentration in cerebral tissue (82±14 µM) with respect to extracerebral tissue (30±7 µM). By contrast, there was no significant difference between the hemoglobin saturation measured in cerebral tissue (56%±10%) and extracerebral tissue (62%±6%). To our knowledge, this is the first time that an inversion procedure in the frequency domain with six unknown parameters with no other prior knowledge is used for the retrieval of the optical coefficients and top layer thickness with high accuracy on two-layered media. Our absolute measurements of cerebral hemoglobin concentration and saturation are based on the discrimination of extracerebral and cerebral tissue layers, and they can enhance the impact of NIRS for cerebral hemodynamics and oxygenation assessment both in the research arena and clinical practice.

Effect of dependent scattering on the optical properties of Intralipid tissue phantoms

The calibration of optical tissue-simulating phantoms remains an open question in spite of the many techniques proposed for accurate measurements of optical properties. As a consequence, a reference phantom with well known optical properties is still missing. As a first step towards a reference phantom we have recently proposed to use dilutions of Intralipid 20%. In this paper we discuss a matter that is commonly ignored when dilutions are prepared, i.e., the possibility of deviations from the simple linear relationships between the optical properties of the dilution and the Intralipid concentration due to the effects of dependent scattering. The results of an experimental investigation showed that dependent scattering does not affect absorption. As for the reduced scattering coefficient the effect can be described adding a term proportional to the square of the concentration. However, for concentrations of interest for tissue optics deviations from linearity remain within about 2%. The experimental investigation also showed that the microphysical properties of Intralipid are not affected by dilution. These results show the possibility to easily obtain a liquid diffusive phantom whose optical properties are known with error smaller than about 1%. Due to the intrinsic limitations of the different techniques proposed for measuring the optical properties it seems difficult to obtain a similar accuracy for solid phantoms.

Modeling the performance characteristics of computed radiography (CR) systems

Computed radiography (CR) using storage phosphors is widely used in digital radiography and mammography. A cascaded linear systems approach wherein several parameter values were estimated using Monte Carlo methods was used to model the image formation process of a single-side read ;;flying spot'' CR system using a granular phosphor. Objective image quality metrics such as modulation transfer function and detective quantum efficiency were determined using this model and show good agreement with published empirical data. A model such as that addressed in this work could allow for improved understanding of the effect of storage phosphor physical properties and CR reader parameters on objective image quality metrics for existing and evolving CR systems.

Performance assessment of photon migration instruments: the MEDPHOT protocol

We propose a comprehensive protocol for the performance assessment of photon migration instruments. The protocol has been developed within the European Thematic Network MEDPHOT (optical methods for medical diagnosis and monitoring of diseases) and is based on five criteria: accuracy, linearity, noise, stability, and reproducibility. This protocol was applied to a total of 8 instruments with a set of 32 phantoms, covering a wide range of optical properties.