Absorption spectroscopy in tissue-simulating materials: a theoretical and experimental study of photon paths

Spatially resolved reflectance spectroscopy for diagnosis of cervical precancer: Monte Carlo modeling and comparison to clinical measurements

We present Monte Carlo modeling studies to provide a quantitative understanding of contrast observed in spatially resolved reflectance spectra of normal and highly dysplastic cervical tissue. Simulations have been carried out to analyze the sensitivity of spectral measurements to a range of changes in epithelial and stromal optical properties that are reported to occur as dysplasia develops and to predict reflectance spectra of normal and highly dysplastic tissue at six different source-detector separations. Simulation results provide important insights into specific contributions of different optical parameters to the overall spectral response. Predictions from simulations agree well with in vivo measurements from cervical tissue and successfully describe spectral differences observed in reflectance measurements from normal and precancerous tissue sites. Penetration depth statistics of photons detected at the six source-detector separations are also presented to reveal the sampling depth profile of the fiber-optic probe geometry simulated. The modeling studies presented provide a framework to meaningfully interpret optical signals obtained from epithelial tissues and to optimize design of optical sensors for in vivo reflectance measurements for precancer detection. Results from this study can facilitate development of analytical photon propagation models that enable inverse estimation of diagnostically relevant optical parameters from in vivo reflectance measurements.

Noninvasive diffuse optical measurement of blood flow and blood oxygenation for monitoring radiation therapy in patients with head and neck tumors: a pilot study

This pilot study explores the potential of noninvasive diffuse correlation spectroscopy (DCS) and diffuse reflectance spectroscopy (DRS) for monitoring early relative blood flow (rBF), tissue oxygen saturation (StO(2)), and total hemoglobin concentration (THC) responses to chemo-radiation therapy in patients with head and neck tumors. rBF, StO(2), and THC in superficial neck tumor nodes of eight patients are measured before and during the chemo-radiation therapy period. The weekly rBF, StO(2), and THC kinetics exhibit different patterns for different individuals, including significant early blood flow changes during the first two weeks. Averaged blood flow increases (52.7+/-9.7)% in the first week and decreases (42.4+/-7.0)% in the second week. Averaged StO(2) increases from (62.9+/-3.4)% baseline value to (70.4+/-3.2)% at the end of the second week, and averaged THC exhibits a continuous decrease from pretreatment value of (80.7+/-7.0) [microM] to (73.3+/-8.3) [microM] at the end of the second week and to (63.0+/-8.1) [microM] at the end of the fourth week of therapy. These preliminary results suggest daily diffuse-optics-based therapy monitoring is feasible during the first two weeks and may have clinical promise.

Characterization of optical parameters with a human forearm at the region from 1.15 to 1.52 microm using diffuse reflectance measurements

Time- and space-resolved diffuse reflectance measurements were used to identify the optical parameters, the reduced scattering and absorption coefficients, of bulk living tissue in the region from 1.15 to 1.52 microm. Although in this region the detector was limited in its temporal resolution, we applied a peak-time shift analysis successfully to determine these coefficients in a human forearm, and then determined the absorption spectrum by space-resolved diffuse reflectance measurements. The absorption spectrum of a water content of 52% determined by magnetic resonance imaging experiments is in good agreement with the absorption coefficient obtained by optical measurements. Moreover, magnetic resonance imaging measurements suggest that the deviation of the absorption coefficients from the water spectrum in the strong water absorption band is caused by the heterogeneity of water distribution in tissue: the low content of water in the skin. These findings indicate that this optical method is potentially applicable to the non-invasive measurement of water in tissue, especially in a region lower than about 1.3-1.35 microm, which may be useful in monitoring oedema and tissue swelling.

Image analysis methods for diffuse optical tomography

Three major analytical tools in imaging science are summarized and demonstrated relative to optical imaging in vivo. Standard resolution testing is optimal when infinite contrast is used and hardware evaluation is the goal. However, deep tissue imaging of absorption or fluorescent contrast agents in vivo often presents a different problem, which requires contrast-detail analysis. This analysis shows that the minimum detectable sizes are in the range of 1/10 the outer diameter, whereas minimum detectable contrast values are in the range of 10 to 20% relative to the continuous background values. This is estimated for objects being in the center of the domain being imaged, and as the heterogeneous region becomes closer to the surface, the lower limit on size and contrast can become arbitrarily low and more dictated by hardware specifications. Finally, if human observer detection of abnormalities in the images is the goal, as is standard in most radiological practice, receiver operating characteristic (ROC) curve and location receiver operating characteristic curve (LROC) are used. Each of these three major areas of image interpretation and analysis are reviewed in the context of medical imaging as well as how they are used to quantify the performance of diffuse optical imaging of tissue.

Modulation of an optical needle's reflectivity alters the average photon path through scattering media

We introduce the concept of deliberate placement of absorbers to alter the average path of photons through tissue for a biomedical optical device. By changing the reflectivity of a needle that separates a source and detector, the average photon path through a turbid medium can be changed. Totally reflective needles have photon scattering density functions similar to a point source and detector in an infinite medium. An absorbing needle moves the average photon path of photons that reach the detector away from the needle. Thus, by modulating the reflectivity of the needle, it is possible to modify the sensitive volume, and simple tomography data should be possible. These results are confirmed by Monte Carlo simulations and experiment. Experiments include moving a black target relative to an optical "needle" and measuring the resulting intensity and phase lag of light reaching a detector at the distal end of the needle.

Polarization-correlation mapping of biological tissue coherent images

We investigate the statistical polarization parameters of biological tissue histological section images with different morphological structure. First we outline the results of polarization coordinate mapping and analysis of the statistics of the first to fourth orders of biological tissue image polarization azimuth and ellipticities. Second, we study the statistics of the first to fourth orders of coordinate distributions of the complex degree of mutual polarization (CDMP) of biological tissue images. Finally, we consider the diagnostic possibilities of investigating 2-D distributions of CDMP of images that correspond to physiologically normal and degeneratively and/or dystrophycally changed biological tissues that are being analyzed.

Reflectance spectroscopy for diagnosis of epithelial precancer: model-based analysis of fiber-optic probe designs to resolve spectral information from epithelium and stroma

Reflectance spectroscopy is a promising technology for detection of epithelial precancer. Fiber-optic probes that selectively collect scattered light from both the epithelium and the underlying stroma are likely to improve diagnostic performance of in vivo reflectance spectroscopy by revealing diagnostic features unique to each layer. We present Monte Carlo models with which to evaluate fiber-optic probe geometries with respect to sampling depth and depth resolution. We propose a probe design that utilizes half-ball lens coupled source and detector fibers to isolate epithelial scattering from stromal scattering and hence to resolve spectral information from the two layers. The probe is extremely compact and can provide easy access to different organ sites.

Characterization of normal breast tissue heterogeneity using time-resolved near-infrared spectroscopy

In recent years, extensive efforts have been made in developing near-infrared optical techniques to be used in detection and diagnosis of breast cancer. Variations in optical properties of normal breast tissue set limits to the performance of such techniques and must therefore be thoroughly examined. In this paper, we present intra- and intersubject as well as contralateral variations of optical and physiological properties in breast tissue as measured by using four-wavelength time-resolved spectroscopy (at 660, 786, 916 and 974 nm). In total, 36 volunteers were examined at five regions at each breast. Optical properties (absorption, mu(a), and reduced scattering, mu'(s)) are derived by employing diffusion theory. The use of four wavelengths enables determination of main tissue chromophores (haemoglobin, water and lipids) as well as haemoglobin oxygenation. Variations in all evaluated properties seen over the entire breast are approximately twice those for small-scale heterogeneity (millimetre scale). Intrasubject variations in optical properties are almost in all cases below 20% for mu'(s), and 40% for mu(a). Overall variations in water, lipid and haemoglobin concentrations are all in the order of 20%. Oxygenation is the least variable of the quantities evaluated, overall intrasubject variations being 6% on average. Extracted physiological properties confirm differences between pre- and post-menopausal breast tissue. Results do not indicate systematic differences between left and right breasts.

Sampling tissue volumes using frequency-domain photon migration

We demonstrate the use of Monte Carlo simulations to generate photon scattering density functions (PSDFs) that represent the tissue volume sampled by steady-state and frequency-domain photon migration. We use these results to illustrate how scaling laws can be developed to determine the mean sampling depth of the multiply scattered photons detected by photon migration methods that remain valid outside the bounds of the standard diffusion approximation, i.e., at small source-detector separations and in media where the optical absorption is significant relative to scattering. Using both the PSDF computation and the newly formulated scaling laws, we focus on a comprehensive description of the effects of source modulation frequency, optical absorption, and source-detector separation on the depth of the sampled tissue volume as well as the sensitivity of frequency-domain photon migration measurements to the presence of a localized absorption heterogeneity.

Determination of the optical properties of a two-layer tissue model by detecting photons migrating at progressively increasing depths

We have investigated a method for solving the inverse problem of determining the optical properties of a two-layer turbid model. The method is based on deducing the optical properties (OPs) of the top layer from the absolute spatially resolved reflectance that results from photon migration within only the top layer by use of a multivariate calibration model. Then the OPs of the bottom layer are deduced from relative frequency-domain (FD) reflectance measurements by use of the two-layer FD diffusion model. The method was validated with Monte Carlo FD reflectance profiles and experimental measurements of two-layer phantoms. The results showed that the method is useful for two-layer models with interface depths of >5 mm; the OPs were estimated, within a relatively short time (<1 min), with a mean error of <10% for the Monte Carlo reflectance profiles and with errors of <25% for the phantom measurements.

Comparison of spatially and temporally resolved diffuse-reflectance measurement systems for determination of biomedical optical properties

Time-resolved and spatially resolved measurements of the diffuse reflectance from biological tissue are two well-established techniques for extracting the reduced scattering and absorption coefficients. We have performed a comparison study of the performance of a spatially resolved and a time-resolved instrument at wavelengths 660 and 786 nm and also of an integrating-sphere setup at 550-800 nm. The first system records the diffuse reflectance from a diode laser by means of a fiber bundle probe in contact with the sample. The time-resolved system utilizes picosecond laser pulses and a single-photon-counting detection scheme. We extracted the optical properties by calibration using known standards for the spatially resolved system, by fitting to the diffusion equation for the time-resolved system, and by using an inverse Monte Carlo model for the integrating sphere. The measurements were performed on a set of solid epoxy tissue phantoms. The results showed less than 10% difference in the evaluation of the reduced scattering coefficient among the systems for the phantoms in the range 9-20 cm(-1), and absolute differences of less than 0.05 cm(-1) for the absorption coefficient in the interval 0.05-0.30 cm(-1).

Coregistration of diffuse optical spectroscopy and magnetic resonance imaging in a rat tumor model

We report coregistration of near-infrared diffuse optical spectroscopy (DOS) and magnetic resonance imaging (MRI) for the study of animal model tumors. A combined broadband steady-state and frequency-domain apparatus was used to determine tissue oxyhemoglobin, deoxyhemoglobin, and water concentration locally in tumors. Simultaneous MRI coregistration provided structural (T2-weighted) and contrast-enhanced images of the tumor that were correlated with the optical measurements. By use of Monte Carlo simulations, the optically sampled volume was superimposed on the MR images, showing precisely which tissue structure was probed optically. DOS and MRI coregistration measurements were performed on seven rats over 20 days and were separated into three tumor tissue classifications: viable, edematous, and necrotic. A ratio of water concentration to total hemoglobin concentration, as measured optically, was performed for each tissue type and showed values for edematous tissue to be greater than viable tissue (1.2 +/- 0.49 M/microM versus 0.48 +/- 0.15 M/microM). Tissue hemoglobin oxygen saturation (StO2) also showed a large variation between tissue types: viable tissue had an optically measured StO2 value of 61 +/- 5%, whereas StO2 determined for necrotic tissue was 43 +/- 6%.

Accelerated Monte Carlo models to simulate fluorescence spectra from layered tissues

Two efficient Monte Carlo models are described, facilitating predictions of complete time-resolved fluorescence spectra from a light-scattering and light-absorbing medium. These are compared with a third, conventional fluorescence Monte Carlo model in terms of accuracy, signal-to-noise statistics, and simulation time. The improved computation efficiency is achieved by means of a convolution technique, justified by the symmetry of the problem. Furthermore, the reciprocity principle for photon paths, employed in one of the accelerated models, is shown to simplify the computations of the distribution of the emitted fluorescence drastically. A so-called white Monte Carlo approach is finally suggested for efficient simulations of one excitation wavelength combined with a wide range of emission wavelengths. The fluorescence is simulated in a purely scattering medium, and the absorption properties are instead taken into account analytically afterward. This approach is applicable to the conventional model as well as to the two accelerated models. Essentially the same absolute values for the fluorescence integrated over the emitting surface and time are obtained for the three models within the accuracy of the simulations. The time-resolved and spatially resolved fluorescence exhibits a slight overestimation at short delay times close to the source corresponding to approximately two grid elements for the accelerated models, as a result of the discretization and the convolution. The improved efficiency is most prominent for the reverse-emission accelerated model, for which the simulation time can be reduced by up to two orders of magnitude.

Quantitative oximetry of breast tumors: a near-infrared method that identifies two optimal wavelengths for each tumor

We present a noninvasive optical method to measure the oxygen saturation of hemoglobin in breast lesions. This method introduces the novel concept that the best choice of near-infrared wavelengths for noninvasive tumor oximetry consists of a wavelength pair (lambda1, lambda2) within the range 680-880 nm, where the specific values of lambda1 and lambda2 depend on the optical properties of the specific tumor under examination. Our method involves two steps: (1) identify the optimal wavelength pair for each tumor and (2) measure the tumor oxygenation using the optical data at the two selected wavelengths. We have tested our method by acquiring experimental optical data from turbid media containing cylindrical or irregularly shaped inhomogeneities and by computing theoretical data for the case of spherical lesions embedded in a highly scattering medium. We have found that our optical method can provide accurate and quantitative measurements of the oxygenation of embedded lesions without requiring knowledge of their size, shape, and depth.

Non-invasive optical monitoring of the newborn piglet brain using continuous-wave and frequency-domain spectroscopy

We have used continuous-wave (CW) and frequency-domain spectroscopy to investigate the optical properties of the newborn piglet brain in vivo and non-invasively. Three anaesthetized, intubated, ventilated and instrumented newborn piglets were placed into a stereotaxic instrument for optimal experimental stability, reproducible probe-to-scalp optical contact and 3D adjustment of the optical probe. By measuring the absolute values of the brain absorption and reduced scattering coefficients at two wavelengths (758 and 830 nm), frequency-domain spectroscopy provided absolute readings (in contrast to the relative readings of CW spectroscopy) of cerebral haemoglobin concentration and saturation during experimentally induced perturbations in cerebral haemodynamics and oxygenation. Such perturbations included a modulation of the inspired oxygen concentration, transient brain asphyxia, carotid artery occlusion and terminal brain asphyxia. The baseline cerebral haemoglobin saturation and concentration, measured with frequency-domain spectroscopy, were about 60% and 42 microM respectively. The cerebral saturation values ranged from a minimum of 17% (during transient brain asphyxia) to a maximum of 80% (during recovery from transient brain asphyxia). To analyse the CW optical data, we have (a) derived a mathematical relationship between the cerebral optical properties and the differential pathlength factor and (b) introduced a method based on the spatial dependence of the detected intensity (dc slope method). The analysis of the cerebral optical signals associated with the arterial pulse and with respiration demonstrates that motion artefacts can significantly affect the intensity recorded from a single optode pair. Motion artefacts can be strongly reduced by combining data from multiple optodes to provide relative readings in the dc slope method. We also report significant biphasic changes (initial decrease and successive increase) in the reduced scattering coefficient measured in the brain after the piglet had been sacrificed.

Influence of a superficial layer in the quantitative spectroscopic study of strongly scattering media

We have experimentally investigated the meaning of the effective optical absorption [mu(a)((eff))] and the reduced scattering [mu(s)?((eff))] coefficients measured on the surfaces of two-layered turbid media, using the diffusion equation for homogeneous, semi-infinite media. We performed frequency-domain spectroscopy in a reflectance geometry, using source-detector distances in the range 1.5-4.5 cm. We measured 100 samples, each made of one layer (thickness in the range 0.08-1.6 cm) on top of one semi-infinite block. The optical properties of the samples were similar to those of soft tissues in the near infrared. We found that the measured effective optical coefficients are representative of the underlying block if the superficial layer is less than ~0.4 cm thick, whereas they are representative of the superficial layer if it is more than ~1.3 cm thick.

Absorbed photodynamic dose from pulsed versus continuous wave light examined with tissue-simulating dosimeters

A dosimetric system has been developed to measure the spatially resolved light dose absorbed by a photosensitizer in a tissue-simulating medium. These gelatin-based dosimeters had macroscopic optical scattering and absorption properties that are typical for homogeneous tissue and contained the photosensitizer benzoporphyrin derivative monoacid (BPD-MA). A reporter molecule, 2?7?-dichlorofluorescin diacetate (DCF-DA), served as an actinometer, which could be photosensitized by BPD-MA to generate a highly fluorescent photoproduct. The relative photosensitizing efficiencies of high-intensity pulsed and cw laser light were compared in these tissue-simulating dosimeters. These measurements demonstrate an increase in penetration for pulsed light as compared with cw light in the dosimeters. A numerical simulation of the light propagation based on optical diffusion theory was used along with the energy levels of the photosensitizer molecule to examine the mechanisms involved in the absorbed dose. The increased penetration of high-intensity pulsed light was due to a transient decrease in the absorption of the photosensitizer, resulting from saturation of the photosensitizer optical transitions. This study provides the first direct comparison of the photodynamic dose absorbed by a photosensitizer using both high-intensity pulsed and cw laser light in a tissue-simulating medium. These measurements demonstrate that a small increase in depth of treatment is possible with pulsed laser light as compared with cw laser light simply on the basis of the unique photochemistry of the photosensitizer. However, this effect still needs to be examined carefully in tumor tissue, where other biological or chemical effects may become significant.

Forward-adjoint fluorescence model: Monte Carlo integration and experimental validation

The adjoint form of the photon transport equation is applied to a generalized fluorescence detection problem, and its accuracy is empirically tested. This approach can be interpreted as mathematically reversing the temporal flow of fluorescent photons; that is, they are tracked from the detector back to potential sites of origin in the scattering medium. The result is a distribution of potential fluorescing sites that, when properly normalized, gives a probability field of the relative importance of the photon starting position and direction to the resulting signal. This adjoint solution can be combined with the temporally forward-derived distribution of absorbed excitation photons to evaluate the fluorescence excitation detection scheme. This bypasses the normal, temporal derivation wherein the fluorescence transport solution is dependent on the result of the excitation transport solution.

Accuracy of noninvasive in vivo measurements of photosensitizer uptake based on a diffusion model of reflectance spectroscopy

This study compares the photosensitizer concentration measured noninvasively in vivo by diffuse reflectance spectroscopy with the results of postmortem tissue solubilization and fluorometric assay. The reflectance spectrometer consists of a fiber optic surface probe, spectrometer and charge-coupled device (CCD) array detector. The surface probe has eight detection fibers separated from the light source fiber by distances ranging from 0.85 to 10 mm. The imaging spectrometer disperses the light from each detector fiber onto the two-dimensional CCD array, while maintaining spatial separation of each individual spectrum. A single exposure of the CCD therefore captures the reflectance spectrum ar eight distances and over a range of 300 nm. From the spectra, the tissue's optical scattering and absorption coefficients are determined using a diffusion model of light propagation. Changes in the tissue absorption are used to estimate the photosensitizer concentration. Normal New Zealand White rabbits were injected with aluminum phthalocyanine tetrasulfonate (AlPcS4) and probe measurements made 24 h after injection on the dorsal skin, on muscle after surgically turning the skin back and on liver. For skin, the noninvasive estimate to proportional to the true concentration but low by a factor of 3. Based on Monte Carlo modeling of multilayered systems, this underestimate is attributed to the layered structure of the skin and nonuniform AlPcS4 distribution. A comparison of the noninvasive concentration estimates to the postmortem assay results finds good agreement for liver tissue even though application of the diffusion model is not strictly justified.

A solubilization technique for photosensitizer quantification in ex vivo tissue samples

The determination of the photosensitizer concentration in ex vivo tissue samples is commonly used for pharmacokinetic and dosimetric studies of photodynamic therapy, both clinically and pre-clinically. In this report, a new method is presented based on tissue solubilization and subsequent fluorometry. This method has the advantages of good sensitivity, accuracy and reproducibility, as well as low cost and ease of handling of the tissue samples. The method was tested for six different photosensitizers in a variety of tissues. The accuracy and concentration detection limits are compared with those of other published extraction methods.

The finite-element method for the propagation of light in scattering media: frequency domain case

A frequency domain light transport model to simulate the transillumination of a scattering object with radio frequency intensity modulated light is presented. The model is based on the diffusion approximation to the radiative transfer equation and uses a finite-element model to allow for complex geometries and an inhomogeneous distribution of absorption and scattering. It calculates the complex photon density within the object and the complex exitance on the boundary of the object. The model is validated against an analytic Green's function model for a circular geometry in the homogeneous case, and its accuracy is investigated for a range of mesh resolutions, optical parameters, and modulation frequencies.

Predictions and measurements of scattering and absorption over broad wavelength ranges in tissue phantoms

Predictions from Mie theory regarding the wavelength dependence of scattering in tissue from the near UV to the near IR are discussed and compared with experiments on tissue phantoms. For large fiber separations it is shown that rapid, simultaneous measurements of the elastic scatter signal for several fiber separations can yield the absorption coefficient and reduced scattering coefficient. With this information, the size of the scattering particles can be estimated, and this is done for Intralipid. Measurements made at smaller source detector separations support Mie theory calculations, demonstrating that the sensitivity of elastic scatter measurements to morphological features, such as scatterer size, is enhanced when the distance between the source and detector fibers is small.

Optical polarization imaging

The temporal profiles of the parallel and perpendicular polarization components of a light pulse backscattered from a scattering medium are different. The depth of penetration into the tissue and depolarization of the backscattered light depend on the scattering and absorption characteristics of the tissue. Based on these facts, a novel technique is demonstrated for noninvasive surface and beneath-the-surface imaging of biological systems.

Photon diffusion on the velocity sphere

The concept of photon diffusion on the velocity sphere and the non-Euclidean diffusion equation (NED) are introduced to describe photon migration in highly forward-scattering random media. The NED covers the ballistic, transient, and developed diffusion modes of photon migration in random media. An approximate analytic solution to the NED is presented.

Fermat photons in turbid media: an exact analytic solution for most favorable paths-a step toward optical tomography

Photon propagation in a turbid medium characterized by highly forward-directed light scattering can be described by a generalized Fermat principle, a counterpart of the conventional Fermat principle for transparent media. An exact analytic solution for the most favorable photon paths (Fermat paths) has been found and is investigated. An analysis intended to encourage and facilitate experimental observation of Fermat photons is presented. The feasibility of using the theory in time-resolved optical tomography for the recognition of objects hidden in turbid media is discussed.