Performance assessment of photon migration instruments: the MEDPHOT protocol

Cerebral time domain-NIRS: reproducibility analysis, optical properties, hemoglobin species and tissue oxygen saturation in a cohort of adult subjects

The reproducibility of cerebral time-domain near-infrared spectroscopy (TD-NIRS) has not been investigated so far. Besides, reference intervals of cerebral optical properties, of absolute concentrations of deoxygenated-hemoglobin (HbR), oxygenated-hemoglobin (HbO), total hemoglobin (HbT) and tissue oxygen saturation (StO₂) and their variability have not been reported. We have addressed these issues on a sample of 88 adult healthy subjects. TD-NIRS measurements at 690, 785, 830 nm were fitted with the diffusion model for semi-infinite homogenous media. Reproducibility, performed on 3 measurements at 5 minutes intervals, ranges from 1.8 to 6.9% for each of the hemoglobin species. The mean ± SD global values of HbR, HbO, HbT, StO₂ are respectively 24 ± 7 μM, 33.3 ± 9.5 μM, 57.4 ± 15.8 μM, 58 ± 4.2%. StO₂ displays the narrowest range of variability across brain regions.

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.

Miniaturized pulsed laser source for time-domain diffuse optics routes to wearable devices

We validate a miniaturized pulsed laser source for use in time-domain (TD) diffuse optics, following rigorous and shared protocols for performance assessment of this class of devices. This compact source (12×6  mm2) has been previously developed for range finding applications and is able to provide short, high energy (∼100  ps, ∼0.5  nJ) optical pulses at up to 1 MHz repetition rate. Here, we start with a basic level laser characterization with an analysis of suitability of this laser for the diffuse optics application. Then, we present a TD optical system using this source and its performances in both recovering optical properties of tissue-mimicking homogeneous phantoms and in detecting localized absorption perturbations. Finally, as a proof of concept of in vivo application, we demonstrate that the system is able to detect hemodynamic changes occurring in the arm of healthy volunteers during a venous occlusion. Squeezing the laser source in a small footprint removes a key technological bottleneck that has hampered so far the realization of a miniaturized TD diffuse optics system, able to compete with already assessed continuous-wave devices in terms of size and cost, but with wider performance potentialities, as demonstrated by research over the last two decades.

Broadband diffuse optical characterization of elastin for biomedical applications

Elastin is a key structural protein of dynamic connective tissues widely found in the extracellular matrix of skin, arteries, lungs and ligaments. It is responsible for a range of diseases related to aging of biological tissues. The optical characterization of elastin can open new opportunities for its investigation in biomedical studies. In this work, we present the absorption spectra of elastin using a broadband (550-1350nm) diffuse optical spectrometer. Distortions caused by fluorescence and finite bandwidth of the laser source on estimated absorption were effectively accounted for in measurements and data analysis and compensated. A comprehensive summary and comparison between collagen and elastin is presented, highlighting distinct features for its accurate quantification in biological applications.

Probe-hosted silicon photomultipliers for time-domain functional near-infrared spectroscopy: phantom and in vivo tests

We report the development of a compact probe for time-domain (TD) functional near-infrared spectroscopy (fNIRS) based on a fast silicon photomultiplier (SiPM) that can be put directly in contact with the sample without the need of optical fibers for light collection. We directly integrated an avalanche signal amplification stage close to the SiPM, thus reducing the size of the detection channel and optimizing the signal immunity to electromagnetic interferences. The whole detection electronics was placed in a plastic screw holder compatible with the electroencephalography standard cap for measurement on brain or with custom probe holders. The SiPM is inserted into a transparent and insulating resin to avoid the direct contact of the scalp with the 100-V bias voltage. The probe was integrated in an instrument for TD fNIRS spectroscopy. The system was characterized on tissue phantoms in terms of temporal resolution, responsivity, linearity, and capability to detect deep absorption changes. Preliminary in vivo tests on adult volunteers were performed to monitor hemodynamic changes in the arm during a cuff occlusion and in the brain cortex during a motor task.

Quantification in time-domain diffuse optical tomography using Mellin-Laplace transforms

Simulations and phantom measurements are used to evaluate the ability of time-domain diffuse optical tomography using Mellin-Laplace transforms to quantify the absorption perturbation of centimetric objects immersed at depth 1-2 cm in turbid media. We find that the estimated absorption coefficient varies almost linearly with the absorption change in the range of 0-0.15 cm^-1 but is underestimated by a factor that depends on the inclusion depth (~2, 3 and 6 for depths of 1.0, 1.5 and 2.0 cm respectively). For larger absorption changes, the variation is sublinear with ~20% decrease for δμ_a = 0.37 cm^-1. By contrast, constraining the absorption change to the actual volume of the inclusion may considerably improve the accuracy and linearity of the reconstructed absorption.

New frontiers in time-domain diffuse optics, a review

The recent developments in time-domain diffuse optics that rely on physical concepts (e.g., time-gating and null distance) and advanced photonic components (e.g., vertical cavity source-emitting laser as light sources, single photon avalanche diode, and silicon photomultipliers as detectors, fast-gating circuits, and time-to-digital converters for acquisition) are focused. This study shows how these tools could lead on one hand to compact and wearable time-domain devices for point-of-care diagnostics down to the consumer level and on the other hand to powerful systems with exceptional depth penetration and sensitivity.

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.

Optical phantoms with adjustable subdiffusive scattering parameters

A new epoxy-resin-based optical phantom system with adjustable subdiffusive scattering parameters is presented along with measurements of the intrinsic absorption, scattering, fluorescence, and refractive index of the matrix material. Both an aluminium oxide powder and a titanium dioxide dispersion were used as scattering agents and we present measurements of their scattering and reduced scattering coefficients. A method is theoretically described for a mixture of both scattering agents to obtain continuously adjustable anisotropy values g between 0.65 and 0.9 and values of the phase function parameter γ in the range of 1.4 to 2.2. Furthermore, we show absorption spectra for a set of pigments that can be added to achieve particular absorption characteristics. By additional analysis of the aging, a fully characterized phantom system is obtained with the novelty of g and γ parameter adjustment.

Development of traceable measurement of the diffuse optical properties of solid reference standards for biomedical optics at National Institute of Standards and Technology

The development of a national reference instrument dedicated to the measurement of the scattering and absorption properties of solid tissue-mimicking phantoms used as reference standards is presented. The optical properties of the phantoms are measured with a double-integrating sphere setup in the steady-state domain, coupled with an inversion routine of the adding-doubling procedure that allows for the computation of the uncertainty budget for the measurements. The results are compared to the phantom manufacturer's values obtained by a time-resolved approach. The results suggest that the agreement between these two independent methods is within the estimated uncertainties. This new reference instrument will provide optical biomedical research laboratories with reference values for absolute diffuse optical properties of phantom materials.

Fast silicon photomultiplier improves signal harvesting and reduces complexity in time-domain diffuse optics

We present a proof of concept prototype of a time-domain diffuse optics probe exploiting a fast Silicon PhotoMultiplier (SiPM), featuring a timing resolution better than 80 ps, a fast tail with just 90 ps decay time-constant and a wide active area of 1 mm². The detector is hosted into the probe and used in direct contact with the sample under investigation, thus providing high harvesting efficiency by exploiting the whole SiPM numerical aperture and also reducing complexity by avoiding the use of cumbersome fiber bundles. Our tests also demonstrate high accuracy and linearity in retrieving the optical properties and suitable contrast and depth sensitivity for detecting localized inhomogeneities. In addition to a strong improvement in both instrumentation cost and size with respect to legacy solutions, the setup performances are comparable to those of state-of-the-art time-domain instrumentation, thus opening a new way to compact, low-cost and high-performance time-resolved devices for diffuse optical imaging and spectroscopy.

Optical properties of developing pip and stone fruit reveal underlying structural changes

Analyzing the optical properties of fruits represents a powerful approach for non-destructive observations of fruit development. With classical spectroscopy in the visible and near-infrared wavelength ranges, the apparent attenuation of light results from its absorption or scattering. In horticultural applications, frequently, the normalized difference vegetation index (NDVI) is employed to reduce the effects of varying scattering properties on the apparent signal. However, this simple approach appears to be limited. In the laboratory, with time-resolved reflectance spectroscopy, the absorption coefficient, μa , and the reduced scattering coefficient, μs ', can be analyzed separately. In this study, these differentiated optical properties were recorded (540-940 nm), probing fruit tissue from the skin up to 2 cm depth in apple (Malus × domestica 'Elstar') and plum (Prunus domestica 'Tophit plus') harvested four times (65-145 days after full bloom). The μa spectra showed typical peak at 670 nm of the chlorophyll absorption. The μs ' at 670 nm in apple changed by 14.7% (18.2-15.5 cm(-1) ), while in plum differences of 41.5% (8.5-5.0 cm(-1) ) were found. The scattering power, the relative change of μs ', was zero in apple, but enhanced in plum over the fruit development period. This mirrors more isotropic and constant structures in apple compared with plum. For horticultural applications, the larger variability in scattering properties of plum explains the discrepancy between commercially assessed NDVI values or similar indices and the absolute μa values in plum (R < 0.05), while the NDVI approach appeared reasonable in apple (R ≥ 0.80).

Macroscopic optical imaging technique for wide-field estimation of fluorescence depth in optically turbid media for application in brain tumor surgical guidance

A diffuse imaging method is presented that enables wide-field estimation of the depth of fluorescent molecular markers in turbid media by quantifying the deformation of the detected fluorescence spectra due to the wavelength-dependent light attenuation by overlying tissue. This is achieved by measuring the ratio of the fluorescence at two wavelengths in combination with normalization techniques based on diffuse reflectance measurements to evaluate tissue attenuation variations for different depths. It is demonstrated that fluorescence topography can be achieved up to a 5 mm depth using a near-infrared dye with millimeter depth accuracy in turbid media having optical properties representative of normal brain tissue. Wide-field depth estimates are made using optical technology integrated onto a commercial surgical microscope, making this approach feasible for real-world applications.

Microfabrication of polydimethylsiloxane phantoms to simulate tumor hypoxia and vascular anomaly

We introduce a microfluidic approach to simulate tumor hypoxia and vascular anomaly. Polydimethylsiloxane (PDMS) phantoms with embedded microchannel networks were fabricated by a soft lithography process. A dialysis membrane was sandwiched between two PDMS slabs to simulate the controlled mass transport and oxygen metabolism. A tortuous microchannel network was fabricated to simulate tumor microvasculature. A dual-modal multispectral and laser speckle imaging system was used for oxygen and blood flow imaging in the tumor-simulating phantom. The imaging results were compared with those of the normal vasculature. Our experiments demonstrated the technical feasibility of simulating tumor hypoxia and vascular anomalies using the proposed PDMS phantom. Such a phantom fabrication technique may be potentially used to calibrate optical imaging devices, to study the mechanisms for tumor hypoxia and angiogenesis, and to optimize the drug delivery strategies.

Diffuse optics using a dual window fast-gated counter

In this paper we demonstrate the advantages of a fast-gated counter in achieving high count-rate and reducing costs of timing equipment in a time-resolved diffuse optical spectroscopy setup. We experimentally prove the equivalence between the fast-gated counter we developed and a traditional time-correlated single-photon counting setup in terms of depth sensitivity and signal-to-noise ratio. Additionally, we show the suitability of this device for bilayer analysis and to estimate the absorption coefficient of homogeneous diffusing media. Finally, we present a proof-of-principle arterial occlusion measurement on a healthy volunteer to validate the proposed approach in a real application. Fast-gated counters can dramatically reduce both costs and complexity in time-resolved multichannel systems, while achieving high count-rate, thus offering a great advantage in applications like brain and muscle functional imaging.

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.

Performance assessment of time-domain optical brain imagers, part 1: basic instrumental performance protocol

Performance assessment of instruments devised for clinical applications is of key importance for validation and quality assurance. Two new protocols were developed and applied to facilitate the design and optimization of instruments for time-domain optical brain imaging within the European project nEUROPt. Here, we present the "Basic Instrumental Performance" protocol for direct measurement of relevant characteristics. Two tests are discussed in detail. First, the responsivity of the detection system is a measure of the overall efficiency to detect light emerging from tissue. For the related test, dedicated solid slab phantoms were developed and quantitatively spectrally characterized to provide sources of known radiance with nearly Lambertian angular characteristics. The responsivity of four time-domain optical brain imagers was found to be of the order of 0.1 m² sr. The relevance of the responsivity measure is demonstrated by simulations of diffuse reflectance as a function of source-detector separation and optical properties. Second, the temporal instrument response function (IRF) is a critically important factor in determining the performance of time-domain systems. Measurements of the IRF for various instruments were combined with simulations to illustrate the impact of the width and shape of the IRF on contrast for a deep absorption change mimicking brain activation.

Determination of reference values for optical properties of liquid phantoms based on Intralipid and India ink

A multi-center study has been set up to accurately characterize the optical properties of diffusive liquid phantoms based on Intralipid and India ink at near-infrared (NIR) wavelengths. Nine research laboratories from six countries adopting different measurement techniques, instrumental set-ups, and data analysis methods determined at their best the optical properties and relative uncertainties of diffusive dilutions prepared with common samples of the two compounds. By exploiting a suitable statistical model, comprehensive reference values at three NIR wavelengths for the intrinsic absorption coefficient of India ink and the intrinsic reduced scattering coefficient of Intralipid-20% were determined with an uncertainty of about 2% or better, depending on the wavelength considered, and 1%, respectively. Even if in this study we focused on particular batches of India ink and Intralipid, the reference values determined here represent a solid and useful starting point for preparing diffusive liquid phantoms with accurately defined optical properties. Furthermore, due to the ready availability, low cost, long-term stability and batch-to-batch reproducibility of these compounds, they provide a unique fundamental tool for the calibration and performance assessment of diffuse optical spectroscopy instrumentation intended to be used in laboratory or clinical environment. Finally, the collaborative work presented here demonstrates that the accuracy level attained in this work for optical properties of diffusive phantoms is reliable.

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.

Metamodeling approach for efficient estimation of optical properties of turbid media from spatially resolved diffuse reflectance measurements

A metamodeling approach is introduced and applied to efficiently estimate the bulk optical properties of turbid media from spatially resolved spectroscopy (SRS) measurements. The model has been trained on a set of liquid phantoms covering a wide range of optical properties representative for food and agricultural products and was successfully validated in forward and inverse mode on phantoms not used for training the model. With relative prediction errors of 10% for the estimated bulk optical properties the potential of this metamodeling approach for the estimation of the optical properties of turbid media from spatially resolved spectroscopy measurements has been demonstrated.

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.

Broadband photon time-of-flight spectroscopy of pharmaceuticals and highly scattering plastics in the VIS and close NIR spectral ranges

We present extended spectroscopic analysis of pharmaceutical tablets in the close near infrared spectral range performed using broadband photon time-of-flight (PTOF) absorption and scattering spectra measurements. We show that the absorption spectra can be used to perform evaluation of the chemical composition of pharmaceutical tablets without need for chemo-metric calibration. The spectroscopic analysis was performed using an advanced PTOF spectrometer operating in the 650 to 1400 nm spectral range. By employing temporal stabilization of the system we achieve the high precision of 0.5% required to evaluate the concentration of tablet ingredients. In order to further illustrate the performance of the system, we present the first ever reported broadband evaluation of absorption and scattering spectra from pure and doped Spectralon®.

Phantoms for diffuse optical imaging based on totally absorbing objects, part 1: Basic concepts

The design of inhomogeneous phantoms for diffuse optical imaging purposes using totally absorbing objects embedded in a diffusive medium is proposed and validated. From time-resolved and continuous-wave Monte Carlo simulations, it is shown that a given or desired perturbation strength caused by a realistic absorbing inhomogeneity of a certain absorption and volume can be approximately mimicked by a small totally absorbing object of a so-called equivalent black volume (equivalence relation). This concept can be useful in two ways. First, it can be exploited to design realistic inhomogeneous phantoms with different perturbation strengths simply using a set of black objects with different volumes. Further, it permits one to grade physiological or pathological changes on a reproducible scale of perturbation strengths given as equivalent black volumes, thus facilitating the performance assessment of clinical instruments. A set of plots and interpolating functions to derive the equivalent black volume corresponding to a given absorption change is provided. The application of the equivalent black volume concept for grading different optical perturbations is demonstrated for some examples.

Tissue phantoms in multicenter clinical trials for diffuse optical technologies

Tissue simulating phantoms are an important part of instrumentation validation, standardization/training and clinical translation. Properly used, phantoms form the backbone of sound quality control procedures. We describe the development and testing of a series of optically turbid phantoms used in a multi-center American College of Radiology Imaging Network (ACRIN) clinical trial of Diffuse Optical Spectroscopic Imaging (DOSI). The ACRIN trial is designed to measure the response of breast tumors to neoadjuvant chemotherapy. Phantom measurements are used to determine absolute instrument response functions during each measurement session and assess both long and short-term operator and instrument reliability.

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.

Assessment of the frequency-domain multi-distance method to evaluate the brain optical properties: Monte Carlo simulations from neonate to adult

The near infrared spectroscopy (NIRS) frequency-domain multi-distance (FD-MD) method allows for the estimation of optical properties in biological tissue using the phase and intensity of radiofrequency modulated light at different source-detector separations. In this study, we evaluated the accuracy of this method to retrieve the absorption coefficient of the brain at different ages. Synthetic measurements were generated with Monte Carlo simulations in magnetic resonance imaging (MRI)-based heterogeneous head models for four ages: newborn, 6 and 12 month old infants, and adult. For each age, we determined the optimal set of source-detector separations and estimated the corresponding errors. Errors arise from different origins: methodological (FD-MD) and anatomical (curvature, head size and contamination by extra-cerebral tissues). We found that the brain optical absorption could be retrieved with an error between 8-24% in neonates and infants, while the error increased to 19-44% in adults over all source-detector distances. The dominant contribution to the error was found to be the head curvature in neonates and infants, and the extra-cerebral tissues in adults.

Intralipid: towards a diffusive reference standard for optical tissue phantoms

Measurements of optical properties carried out at visible and NIR wavelengths on many samples of Intralipid 20% showed a high stability and surprisingly small batch-to-batch variations. Measurements have been carried out in a short time interval using samples from nine different batches with expiry dates spreading over ten years. For the specific reduced scattering coefficient, the values we have obtained, averaged over the nine batches, are 25.9, 21.2, and 18.4 mm(-1) at λ = 632.8, 751, and 833 nm, respectively, and the corresponding maximum deviations from the average were 2.2%, 1.1%, and 1.4%. For the absorption coefficient, we obtained values slightly smaller with respect to the absorption coefficient of pure water at 751 and 833 nm, and slightly larger at 632.8 nm. These results suggest that Intralipid 20% can be the first step towards a diffusive reference standard for tissue-simulating phantoms.

A compact time-resolved system for near infrared spectroscopy based on wavelength space multiplexing

We designed and developed a compact dual-wavelength and dual-channel time-resolved system for near-infrared spectroscopy studies of muscle and brain. The system employs pulsed diode lasers as sources, compact photomultipliers, and time-correlated single photon counting boards for detection. To exploit the full temporal and dynamic range of the acquisition technique, we implemented an approach based on wavelength space multiplexing: laser pulses at the two wavelengths are alternatively injected into the two channels by means of an optical 2×2 switch. In each detection line (i.e., in each temporal window), the distribution of photon time-of-flights at one wavelength is acquired. The proposed approach increases the signal-to-noise ratio and avoids wavelength cross-talk with respect to the typical approach based on time multiplexing. The instrument was characterized on tissue phantoms to assess its properties in terms of linearity, stability, noise, and reproducibility. Finally, it was successfully tested in preliminary in vivo measurements on muscle during standard cuff occlusion and on the brain during a motor cortex response due to hand movements.

Fluorescence lifetime optical tomography with Discontinuous Galerkin discretisation scheme

We develop discontinuous Galerkin framework for solving direct and inverse problems in fluorescence diffusion optical tomography in turbid media. We show the advantages and the disadvantages of this method by comparing it with previously developed framework based on the finite volume discretization. The reconstruction algorithm was used with time-gated experimental dataset acquired by imaging a highly scattering cylindrical phantom concealing small fluorescent tubes. Optical parameters, quantum yield and lifetime were simultaneously reconstructed. Reconstruction results are presented and discussed.

Time-resolved diffusing wave spectroscopy with a CCD camera

We show how time-resolved measurements of the diffuse light transmitted through a thick scattering slab can be performed with a standard CCD camera, thanks to an interferometric protocol. Time-resolved correlations measured at a fixed photon transit time are also presented. The high number of pixels of the camera allows us to attain a quite good sensitivity for a reasonably low acquisition time.

Optimization of the hyperspectral imaging-based spatially-resolved system for measuring the optical properties of biological materials

This paper reports on the optimization and assessment of a hyperspectral imaging-based spatially-resolved system for determination of the optical properties of biological materials over the wavelengths of 500-1,000 nm. Twelve model samples covering a wide range of absorption and reduced scattering coefficients were created to validate the hyperspectral imaging system, and their true values of absorption and reduced scattering coefficients were determined and then cross-validated using three commonly used methods (i.e., transmittance, integrating sphere, and empirical equation). Light beam and source-detector distance were optimized through Monte Carlo simulations and experiments for the model samples. The optimal light beam should be of Gaussian type with the diameter of less than 1 mm, and the optimal minimum and maximum source-detector distance should be 1.5 mm and 10-20 mean free paths, respectively. The optimized hyperspectral imaging-based spatially-resolved system achieved good estimation of the optical parameters.

Reference optical phantoms for diffuse optical spectroscopy. Part 1--Error analysis of a time resolved transmittance characterization method

Development, production quality control and calibration of optical tissue-mimicking phantoms require a convenient and robust characterization method with known absolute accuracy. We present a solid phantom characterization technique based on time resolved transmittance measurement of light through a relatively small phantom sample. The small size of the sample enables characterization of every material batch produced in a routine phantoms production. Time resolved transmittance data are pre-processed to correct for dark noise, sample thickness and instrument response function. Pre-processed data are then compared to a forward model based on the radiative transfer equation solved through Monte Carlo simulations accurately taking into account the finite geometry of the sample. The computational burden of the Monte-Carlo technique was alleviated by building a lookup table of pre-computed results and using interpolation to obtain modeled transmittance traces at intermediate values of the optical properties. Near perfect fit residuals are obtained with a fit window using all data above 1% of the maximum value of the time resolved transmittance trace. Absolute accuracy of the method is estimated through a thorough error analysis which takes into account the following contributions: measurement noise, system repeatability, instrument response function stability, sample thickness variation refractive index inaccuracy, time correlated single photon counting system time based inaccuracy and forward model inaccuracy. Two sigma absolute error estimates of 0.01 cm(-1) (11.3%) and 0.67 cm(-1) (6.8%) are obtained for the absorption coefficient and reduced scattering coefficient respectively.

Chi2 analysis for estimating the accuracy of optical properties derived from time resolved diffuse-reflectance

Weighted residuals and the reduced chi(2) (chi(R) (2)) value are investigated with regard to their relevance for assessing optical property estimates using the diffusion equation for time-resolved measurements in turbid media. It is shown and explained, for all photon counting experiments including lifetime estimation, why chi(R) (2) increases linearly with the number of photons when there is a model bias. Only when a sufficient number of photons has been acquired, chi(R) (2) is a pertinent value for assessing the accuracy of mu(a) and mu(s)' estimates. It was concluded that chi(R) (2) is of particular interest for cases of small interfiber separation, low-level scattering, strong absorption and incorrect measurement of instrument response function. It was also found that chi(R) (2) is less pertinent for judging mu(a) in case of air boundary effects.

Near-infrared photon time-of-flight spectroscopy of turbid materials up to 1400 nm

Photon time-of-flight spectroscopy (PTOFS) is a powerful tool for analysis of turbid materials. We have constructed a time-of-flight spectrometer based on a supercontinuum fiber laser, acousto-optical tunable filtering, and an InP/InGaAsP microchannel plate photomultiplier tube. The system is capable of performing PTOFS up to 1400 nm, and thus covers an important region for vibrational spectroscopy of solid samples. The development significantly increases the applicability of PTOFS for analysis of chemical content and physical properties of turbid media. The great value of the proposed approach is illustrated by revealing the distinct absorption features of turbid epoxy resin. Promising future applications of the approach are discussed, including quantitative assessment of pharmaceuticals, powder analysis, and calibration-free near-infrared spectroscopy.

Characterizing liquid turbid media by frequency-domain photon-migration spectroscopy

We present a wavelength-tunable frequency-domain instrument for the characterization of liquid turbid media. The instrument employs a tunable titanium-sapphire laser modulated by an acousto-optic modulator. The absorption and reduced scattering coefficient of Intralipid(R) 20%, diluted to concentrations of 0.94 to 4.00%, are measured over the wavelength range 710 to 850 nm at 10-nm intervals. The standard measurement errors for the absorption and reduced scattering coefficients are 1 and 2.5%, respectively. Extrapolation to 0% Intralipid(R) concentration gives an absorption coefficient that closely follows that of water, overestimating the absorption of pure water by less than 10%. The reduced scattering coefficient is compared at 750 nm with published results and is found consistent within the experimental error. We compare the reduced scattering coefficient to an estimate based on Mie theory and find the reduced scattering coefficient underestimated the Mie theory result by about 9%.

Combined optical imaging and mammography of the healthy breast: optical contrast derived from breast structure and compression

In this paper, we report new progress in developing the instrument and software platform of a combined X-ray mammography/diffuse optical breast imaging system. Particularly, we focus on system validation using a series of balloon phantom experiments and the optical image analysis of 49 healthy patients. Using the finite-element method for forward modeling and a regularized Gauss-Newton method for parameter reconstruction, we recovered the inclusions inside the phantom and the hemoglobin images of the human breasts. An enhanced coupling coefficient estimation scheme was also incorporated to improve the accuracy and robustness of the reconstructions. The recovered average total hemoglobin concentration (HbT) and oxygen saturation (SO2) from 68 breast measurements are 16.2 microm and 71%, respectively, where the HbT presents a linear trend with breast density. The low HbT value compared to literature is likely due to the associated mammographic compression. From the spatially co-registered optical/X-ray images, we can identify the chest-wall muscle, fatty tissue, and fibroglandular regions with an average HbT of 20.1+/-6.1 microm for fibroglandular tissue, 15.4+/-5.0 microm for adipose, and 22.2+/-7.3 microm for muscle tissue. The differences between fibroglandular tissue and the corresponding adipose tissue are significant (p < 0.0001). At the same time, we recognize that the optical images are influenced, to a certain extent, by mammographical compression. The optical images from a subset of patients show composite features from both tissue structure and pressure distribution. We present mechanical simulations which further confirm this hypothesis.

Combined reconstruction of fluorescent and optical parameters using time-resolved data

We present an algorithm for simultaneous reconstruction of optical parameters, quantum yield, and lifetime in turbid media with embedded fluorescent inclusions. This algorithm is designed in the Fourier domain as an iterative solution of a system of differential equations of the Helmholtz type and does not involve full ill-conditioned matrix computations. The approach is based on allowing the unknown optical parameters, quantum yield, and lifetime to depend on the Fourier spectral parameter. The algorithm was applied to a time-gated experimental data set acquired by imaging a highly scattering cylindrical phantom concealing small fluorescent tubes. Relatively accurate reconstruction demonstrates the potential of the method.

White Monte Carlo for time-resolved photon migration

A novel scheme for fully scalable White Monte Carlo (WMC) has been developed and is used as a forward solver in the evaluation of experimental time-resolved spectroscopy. Previously reported scaling problems are avoided by storing detection events individually, turning spatial and temporal binning into post-simulation activities. The approach is suitable for modeling of both interstitial and noninvasive settings (i.e., infinite and semi-infinite geometries). Motivated by an interest in in vivo optical properties of human prostate tissue, we utilize WMC to explore the low albedo regime of time-domain photon migration--a regime where the diffusion approximation of radiative transport theory breaks down, leading to the risk of overestimating both reduced scattering (mu(s)') and absorption (mu(a)). Experimental work supports our findings and establishes the advantages of Monte Carlo-based evaluation.

Time-resolved optical spectroscopy of wood

We have proposed and experimentally demonstrated that picosecond time-resolved optical spectroscopy in the visible/near-infrared (NIR) region (700-1040 nm) is a useful technique for noninvasive characterization of wood. This technique has been demonstrated on both softwood and hardwood samples treated in different ways simulating the aging process suffered by waterlogged woods. In all the cases, alterations of absorption and scattering spectra were observed, revealing changes of chemical and structural composition.

Fully automated time domain spectrometer for the absorption and scattering characterization of diffusive media

We describe a system for absorption and scattering spectroscopy of diffusive media based on time-resolved reflectance and transmittance measurements. The system is operated with mode-locked lasers tunable in the 550-1050 nm spectral range and on a detection chain based on time-correlated single-photon counting. All measurement procedures such as laser tuning and optimization, signal conditioning, data acquisition, and analysis are completely automated, permitting spectral measurements over the whole range in a few minutes. The criticalities of the system are discussed together with the strategies to compensate them. The Medphot protocol devised for the characterization of photon migration instruments was applied to assess the system performances in terms of accuracy, linearity, noise, stability, and reproducibility. Finally, an example of application of the instrument to the spectroscopy of powders is presented.

Time-resolved optical imager for assessment of cerebral oxygenation

A time-resolved optical instrument allowing for noninvasive assessment of cerebral oxygenation is presented. The instrument is equipped with picosecond diode lasers, fast photodetectors, and time-correlated single photon counting electronics. This technology enables depth-resolved estimation of changes in absorption and, in consequence, assessment of changes in hemoglobin concentrations in the brain cortex. Changes in oxyhemoglobin (HbO(2)) and deoxyhemoglobin (Hb) can be evaluated selectively in extra- and intracerebral tissue compartments using the moments of distributions of times of flight of photons measured at two wavelengths in the near-infrared region. The combination of the data acquired from multiple sources and detectors located on the surface of the head with the depth-resolved analysis, based on the moments, enables imaging of cortex oxygenation. Results of the tests on physical phantoms as well as in vivo validation of the instrument during the motor stimulation experiment are presented.

In vivo time-resolved reflectance spectroscopy of the human forehead

We present an in vivo broadband spectroscopic characterization of the human forehead. Absorption and scattering properties are measured on five healthy volunteers at five different interfiber distances, using time-resolved diffuse spectroscopy and interpreting data with a model of the diffusion equation for a homogeneous semi-infinite medium. A wavelength-tunable mode-locked laser and time-correlated single-photon counting detection are employed, enabling fully spectroscopic measurements in the range of 700-1000 nm. The results show a large variation in the absorption and scattering properties of the head depending on the subject, whereas intrasubject variations, assessed at different interfiber distances, appear less relevant, particularly for what concerns the absorption coefficient. The high intersubject variability observed indicates that a unique set of optical properties for modeling the human head cannot be used correctly. To better interpret the results of the analysis of in vivo measurements, we performed a set of four-layer model Monte Carlo simulations based on different data sets for the optical properties of the human head, partially derived from the literature. The analysis indicated that, when simulated time-resolved curves are fitted with a homogeneous model for the photon migration, the retrieved absorption and reduced scattering coefficients are much closer to superficial layer values (i.e., scalp and skull) than to deeper layer ones (white and gray matter). In particular, for the shorter interfiber distances, the recovered values can be assumed as a good estimate of the optical properties of the first layer.

Calibration of scattering and absorption properties of a liquid diffusive medium at NIR wavelengths. CW method

In spite of many progresses achieved both with theories and with experiments in studying light propagation through diffusive media, a reliable method for accurate measurements of the optical properties of diffusive media at NIR wavelengths is, in our opinion, still missing. It is therefore difficult to create a diffusive medium with well known optical properties to be used as a reference. In this paper we describe a method to calibrate the reduced scattering coefficient, mu'(s) , of a liquid diffusive medium and the absorption coefficient, mu(a), of an absorbing medium with a standard error smaller than 2% both on mu'(s) and on mu(a). The method is based on multidistance measurements of fluence into an infinite medium illuminated by a CW source. The optical properties are retrieved with simple inversion procedures (linear fits) exploiting the knowledge of the absorption coefficient of the liquid into which the diffuser and the absorber are dispersed. In this study Intralipid diluted in water has been used as diffusive medium and Indian ink as absorber. For a full characterization of these media measurements of collimated transmittance have also been carried out, from which the asymmetry factor of the scattering function of Intralipid and the single scattering albedo of Indian ink have been determined.

In vivo optical characterization of human prostate tissue using near-infrared time-resolved spectroscopy

The development of photodynamic therapy into a modality for treatment of prostate cancer calls for reliable optical dosimetry. We employ, for the first time, interstitial time-resolved spectroscopy to determine in vivo optical properties of human prostate tissue. Nine patients are included in the study, and measurements are conducted prior to primary brachytherapy treatment of prostate cancer. Intrasubject variability is examined by measuring across three tissue volumes within each prostate. The time-resolved instrumentation proves its usefulness by producing good signal levels in all measurements. We are able to present consistent values on reduced scattering coefficients (mu(s)'), absorption coefficients (mu(a)), and effective attenuation (mu(eff)) at the wavelengths 660, 786, and 916 nm. At 660 nm, mu(s)' is found to be 9+/-2 cm(-1), and mu(a) is 0.5+/-0.1 cm(-1). Derived values of mu(eff) are in the range of 3 to 4 cm(-1) at 660 nm, a result in good agreement with previously published steady state data. Total hemoglobin concentration (THC) and oxygen saturation are spectroscopically determined using derived absorption coefficients. Derived THC values are fairly variable (215+/-65 microM), while derived values of oxygen saturation are gathered around 75% (76+/-4%). Intrasubject variations in derived parameters correlate (qualitatively) with the heterogeneity exhibited in acquired ultrasound images.

Review of tissue simulating phantoms for optical spectroscopy, imaging and dosimetry

Optical spectroscopy, imaging, and therapy tissue phantoms must have the scattering and absorption properties that are characteristic of human tissues, and over the past few decades, many useful models have been created. In this work, an overview of their composition and properties is outlined, by separating matrix, scattering, and absorbing materials, and discussing the benefits and weaknesses in each category. Matrix materials typically are water, gelatin, agar, polyester or epoxy and polyurethane resin, room-temperature vulcanizing (RTV) silicone, or polyvinyl alcohol gels. The water and hydrogel materials provide a soft medium that is biologically and biochemically compatible with addition of organic molecules, and are optimal for scientific laboratory studies. Polyester, polyurethane, and silicone phantoms are essentially permanent matrix compositions that are suitable for routine calibration and testing of established systems. The most common three choices for scatters have been: (1.) lipid based emulsions, (2.) titanium or aluminum oxide powders, and (3.) polymer microspheres. The choice of absorbers varies widely from hemoglobin and cells for biological simulation, to molecular dyes and ink as less biological but more stable absorbers. This review is an attempt to indicate which sets of phantoms are optimal for specific applications, and provide links to studies that characterize main phantom material properties and recipes.

Scatter correction of transmission near-infrared spectra by photon migration data: quantitative analysis of solids

The scope of this work is a new methodology to correct conventional near-infrared (NIR) data for scattering effects. The technique aims at measuring the absorption coefficient of the samples rather than the total attenuation measured in conventional NIR spectroscopy. The main advantage of this is that the absorption coefficient is independent of the path length of the light inside the sample and therefore independent of the scattering effects. The method is based on time-resolved spectroscopy and modeling of light transport by diffusion theory. This provides an independent measure of the scattering properties of the samples and therefore of the path length of light. This yields a clear advantage over other preprocessing techniques, where scattering effects are estimated and corrected for by using the shape of the measured spectrum only. Partial least squares (PLS) calibration models show that, by using the proposed evaluation scheme, the predictive ability is improved by 50% as compared to a model based on conventional NIR data alone. The method also makes it possible to predict the concentration of active substance in samples with other physical properties than the samples included in the calibration model.