Search

Many approaches to fluorescence tomography utilize some form of regularized nonlinear least-squares algorithm for data inversion, thus requiring repeated computation of the Jacobian sensitivity matrix relating changes in observable quantities, such as emission fluence, to changes in underlying optical parameters, such as fluorescence absorption. An exact adjoint formulation of these sensitivities comprises three terms, reflecting the individual contributions of 1) sensitivities of diffusion and decay coefficients at the emission wavelength, 2) sensitivities of diffusion and decay coefficients at the excitation wavelength, and 3) sensitivity of the emission source term. Simplifying linearity assumptions are computationally attractive in that they cause the first and second terms to drop out of the formulation. The relative importance of the three terms is thus explored in order to determine the extent to which these approximations introduce error. Computational experiments show that, while the third term of the sensitivity matrix has the largest magnitude, the second term becomes increasingly significant as target fluorophore concentration or volume increases. Image reconstructions from experimental data confirm that neglecting the second term results in overestimation of sensitivities and consequently overestimation of the value and volume of the fluorescent target, whereas contributions of the first term are so low that they are probably not worth the additional computational costs.

Fluorescence lifetime spectroscopy in the presence of tissuelike scattering is demonstrated from measurements of phase and modulation ratio as a function of modulation frequency using a pH-sensitive dye, Carboxy Seminaphthofluorescein-1 (C-SNAFL-1). From the optical diffusion equation describing the propagation and generation of fluorescence within solutions of 0.5 microM C-SNAFL-1 containing 2.0% (by volume) of Intralipid as a scatterer, the values of the average lifetime of C-SNAFL-1 were determined as the solution pH varied between 5 and 9. Average lifetime values were found to match those measured using traditional phase-modulation measurement in nonscattering media. Furthermore, the robustness of the spectroscopic technique was demonstrated by conducting lifetime measurements at varying scatterer concentrations (1.5-3.0 vol % Intralipid). These results confirm the approach for analytical sensing in scattering media via fluorescence lifetime kinetics in order to track changes in analyte concentrations.

Fluorescence frequency-domain photon migration measurements were acquired from tissue phantoms, each containing a fluorescent target, by means of area illumination and area detection on the same surface and for the first time, to our knowledge, compared with predictions computed with a numerical solution to the coupled photon diffusion equations. We accomplished area illumination and area detection using a planar, intensity-modulated excitation light source and a gain-modulated intensified charge-coupled device camera, respectively. A 1-ml vessel containing 1-microm solution of Indocyanine Green in 1% Liposyn was immersed 1 cm deep in each 512-ml tissue phantom. For most tissue phantoms, the background surrounding the 1-ml target was composed of Liposyn solution containing Indocyanine Green or 3,3'-Diethylthiatricarbocyanine Iodide such that the target-to-background ratio of fluorescence yield was > or = 10:1. Measurements of fluorescence modulation amplitude and phase were predicted with a mean error ranging from 10.1% to 13.6% and 0.56 degrees to 1.72 degrees, respectively. These numbers are similar to those obtained by use of single-pixel frequency-domain photon migration techniques and validate the potential use of area illumination and area detection for biomedical imaging of tissues. Results also demonstrate that target-to-background ratios of fluorescence yield and fluorescence lifetime significantly affect target detectability.

A three-dimensional image reconstruction for fluorescence-enhanced frequency-domain photon migration (FDPM) measurements in turbid media is developed and investigated for three different simulated measurement types: 1) absolute emission measurement, or emission measurements of phase and amplitude attenuation made for a given incident point source of excitation light; 2) referenced emission measurements made relative to an excitation measurement conducted at a single reference point away from the incident source; and 3) referenced emission measurements made relative to the excitation measurement conducted at identical points of detection. The image reconstruction algorithm employs a gradient-based constrained truncated Newton (CONTN) method which implements a bounding parameter, which can be used to govern the level of contrast used to discriminate tissue volumes from heterogeneous background tissues. Reverse differentiation technique is used to calculate the gradients. Using simulated data with superimposed noise to achieve a signal-to-noise ratio of 55 and 35 dB to mimic experimental excitation and emission FDPM measurements, respectively, we show the robustness of emission measurements referenced to excitation light. We investigate the performance of algorithm CONTN using these measurement techniques and show that the absorption coefficients due to fluorophore are reconstructed by CONTN accurately and efficiently. Furthermore, we demonstrate the performance of the bounding parameter for rejection of background artifacts owing to background tissue heterogeneity.

A novel image-intensified charge-coupled device (ICCD) imaging system has been developed to perform 3D fluorescence tomographic imaging in the frequency-domain using near-infrared contrast agents. The imager is unique since it (i) employs a large tissue-mimicking phantom, which is shaped and sized to resemble a female breast and part of the extended chest-wall region, and (ii) enables rapid data acquisition in the frequency-domain by using a gain-modulated ICCD camera. Diffusion model predictions are compared to experimental measurements using two different referencing schemes under two different experimental conditions of perfect and imperfect uptake of fluorescent agent into a target. From these experimental measurements, three-dimensional images of fluorescent absorption were reconstructed using a computationally efficient variant of the approximate extended Kalman filter algorithm. The current work represents the first time that 3D fluorescence-enhanced optical tomographic reconstructions have been achieved from experimental measurements of the time-dependent light propagation on a clinically relevant breast-shaped tissue phantom using a gain-modulated ICCD camera.

The development of near-infrared fluorescent contrast agents and imaging techniques depends on the deep penetration of excitation light through several centimeters of tissue and the sensitive collection of the re-emitted fluorescence. In this contribution, the sensitivity and depth penetration of various fluorescence-enhanced imaging studies is surveyed and compared with current studies using continuous wave (CW) and frequency-domain photon migration (FDPM) measurements with planar wave illumination of modulated excitation light at 100 MHz and area collection of reemitted fluorescent light using a previously developed modulated intensified charge-coupled device camera system. Fluorescence was generated from nanomolar to micromolar solutions of indocyanine green (ICG) in a 100 microL volume submerged at 1-4 cm depths in a 1% Liposyn solution to mimic tissue scattering properties. Enhanced depth penetration and sensitivity are achieved with optimal filter rejection of excitation light, and FDPM rejection of background light is not achieved using CW methods. We show the ability to detect as few as 100 fmol of ICG from area illumination of 785 nm light (5.5 mW/cm2) and FDPM area collection of 830 nm fluorescent light generated from 3 cm below the phantom surface. The lowered noise floor of FDPM measurements enables greater sensitivity and penetration depth than comparable CW measurements.

The frequency domain photon migration (FDPM) technique was employed to determine mean particle size of pharmaceutical powders. Results show that the FDPM-measured scattering coefficient increases linearly with reciprocal mean particle size of powdered samples. In contrast to near-infrared spectroscopy techniques, FDPM technique enables determination of scattering and absorption separately so that it does not require data pretreatment and chemometric calibration models. In addition, this unique advantage provides more detailed information about powder samples, which can be used as a potential tool for on-line monitoring of not only variation of active pharmaceutical ingredient concentrations from changes in the absorption coefficient but also variation of particle sizes from changes in the scattering coefficient.

Fluorescence frequency-domain photon migration (FDPM) through tissue refers to the propagation of intensity-modulated fluorescent light that originates from tissue-laden fluorophores following illumination with an intensity-modulated excitation light source. FDPM measurements of modulation amplitude and phase are ultimately employed in an inversion algorithm for tomographic reconstruction of interior optical and fluorescent property maps that delineate disease enhanced with fluorescent contrast agent. Because the inverse problem is underdetermined, measurement precision and accuracy crucially impact its solution. Reported here are the precision and accuracy of FDPM measurements acquired using an intensified CCD homodyne detection system. By introducing 32 phase delays between the oscillators used to modulate the intensifier gain and light source intensity at 100 MHz, mean precision is maximized at +/-0.46% and +/-0.26 deg for measurements of modulation amplitude and phase, respectively. Measurement precision improves when the number of phase delays increases. Measurements of fluorescence modulation amplitude and phase, acquired from the surface of a tissue phantom at distances ranging between 0.71 and 3.6 cm from an incident excitation point source, exhibit a mean accuracy of 17% and 1.9 deg, respectively. Measurement accuracy deteriorates with increasing distance from the point source, but for distances up to 1.0 cm from the point source, measurements of fluorescence modulation amplitude and phase exhibit a mean accuracy of 5.4% and 0.30 deg, respectively., (Copyright 2003 Society of Photo-Optical Instrumentation Engineers)

The advent of recent advances in near-infrared laser diodes and fast electro-optic detection has spawned a new research field of diagnostic spectroscopy and imaging based on targeting and reporting exogenous fluorescent agents. This review seeks to concisely address the physics, instrumentation, advancements in tomography, and near-infrared fluorescent contrast agent development that promises selective and specific molecular targeting of diseased tissues. As an example of one area of the field, recent work focusing on pharmacokinetic analysis of fluorophores targeting the epidermal growth factor receptor (EGFR) is presented in a human breast cancer xenograft mouse model to demonstrate specificity of molecularly targeted contrast agents. Finally, a critical evaluation of the limitations and the opportunities for future translation of fluorescence-enhanced optical imaging of deep tissues is presented.

The deep tissue propagation of near-infrared (NIR) light between 700-900 nm offers new opportunities for diagnostic imaging when employing sensitive detection techniques and NIR excitable fluorescent agents that target and report disease and metabolism. Herein, we highlight approaches for illuminating tissues and monitoring the re-emitted fluorescence for tomographic reconstruction, strategies for developing fluorescent dye constructs, and clinical opportunities for fluorescence-enhanced NIR optical imaging.

Using probability distribution analysis to describe the propagation of multiply scattered light between a point source and point detector located a known distance apart, mathematical expressions predicting the volume of sampled powder were determined for infinite and semi-infinite powder beds and then compared to experimental measurements of frequency-domain photon migration (FDPM). Our results show that the volume of powder sampled varies with optical properties and, when using FDPM techniques, with modulation frequency. For a typical measurement in lactose powder, the volume of powder sampled by multiply scattered light propagating between a 1000-microm-diameter point source and point detector pair separated by 8 mm is predicted to be 1.8 and 1.5 cm3 at modulated frequencies of 50 and 100 MHz, respectively. Experimental measurements in lactose powder confirm the predictions.

To investigate fluorescence lifetime spectroscopy in tissue-like scattering, measurements of phase modulation as a function of modulation frequency were made using two fluorescent dyes exhibiting single exponential decay kinetics in a 2% intralipid solution. To experimentally simulate fluorescence multiexponential decay kinetics, we varied the concentration ratios of the two dyes, 3,3-diethylthiatricarbocyanine iodide and indocynanine green (ICG), which exhibit distinctly different lifetimes of 1.33 and 0.57 ns, respectively. The experimental results were then compared with values predicted using the optical diffusion equation incorporating 1) biexponential decay, 2) average of the biexponential decay, as well as 3) stretched exponential decay kinetic models to describe kinetics owing to independent and quenched relaxation of the two dyes. Our results show that while all kinetic models could describe phase-modulation data in nonscattering solution, when incorporated into the diffusion equation, the kinetic parameters failed to likewise predict phase-modulation data in scattering solutions. We attribute the results to the insensitivity of phase-modulation measurements in nonscattering solutions and the inaccuracy of the derived kinetic parameters. Our results suggest the high sensitivity of phase-modulation measurements in scattering solutions may provide greater opportunities for fluorescence lifetime spectroscopy.

A method for inverting measurements made on the surfaces of tissues for recovery of interior optical property maps is demonstrated for sparse near-infrared (NIR) fluorescence measurement sets on large tissue-simulating volumes with highly variable signal-to-noise ratio. A Bayesian minimum-variance reconstruction algorithm compensates for the spatial variability in signal-to-noise ratio that must be expected to occur in actual NIR contrast-enhanced diagnostic medical imaging. Image reconstruction is demonstrated by using frequency-domain photon migration measurements on 256-cm(3) tissue-mimicking phantoms containing none, one, or two 1-cm(3) heterogeneities with 50- to 100-fold greater concentration of Indocyanine Green dye over background levels. The spatial parameter estimate of absorption owing to the dye was reconstructed from only 160 to 296 surface measurements of emission light at 830 nm in response to incident 785-nm excitation light modulated at 100 MHz. Measurement error of acquired fluence at fluorescent emission wavelengths is shown to be highly variable. Convergence and quality of image reconstructions are improved by Bayesian conditioning incorporating (i) experimentally determined measurement error variance, (ii) recursively updated estimates of parameter uncertainty, and (iii) dynamic zonation. The results demonstrate that, to employ NIR fluorescence-enhanced optical imaging for large volumes, reconstruction approaches must account for the large range of signal-to-noise ratio associated with the measurements.

The isotropic scattering coefficients of 143-nm diameter polystyrene latex suspensions were measured using frequency-domain photon migration (FDPM) at 687 and 828 nm as a function of volume fraction (0.05-0.3) and ionic strength (1.0 to 120 mM NaCl equivalents) in order to derive the angle-integrated structure factor, S(q), and structure factor at zero wave vector, S(0). The effective surface charges of the dispersions were estimated by fitting the measured isotropic scattering coefficients at each wavelength as a function of volume fraction to the solution of the Orstein-Zernike integral equation using the hard sphere Yukawa potential model and mean spherical approximation as a closure relation. The estimates of surface charges were comparable at both wavelengths, but decreased with ionic strength. At 120 mM NaCl equivalents, the values of S(0) obtained from FDPM matched those predicted by the Percus-Yevick model, and decreased with volume fraction, consistent with prediction by the Carnahan-Starling equation.

The ultimate success of near-infrared optical tomography rests on the precise measurement of light propagation within tissues or random media, the accurate prediction of these measurements from a light propagation model, and an efficient three-dimensional solution of the inverse imaging problem. To date, optical tomography algorithms have focused on frequency-domain photon migration (FDPM) measurements of phase-delay and amplitude attenuation, which are reported relative to the incident light, even though phase-delay and amplitude of incident light are nearly impossible to measure directly. In this contribution, we examine referenced, fluorescence-enhanced frequency-domain photon migration measured at excitation and/or emission wavelengths and report on a measurement strategy to minimize measurement and calibration error for efficient coupling of data to a distorted Born iterative imaging algorithm. We examine three referencing approaches and develop associated inversion algorithms for (1) normalizing detected emission FDPM data to the predicted emission wave arising from a homogeneous medium, (2) referencing detected emission FDPM data to that detected at a reference point, and (3) referencing detected emission FDPM data to detected excitation FDPM data detected at a reference point. Our results show the latter approach to be practical while reducing the nonlinearity of the inverse problem. Finally, in light of our results, we demonstrate the method for eliminating the influence of source strength and instrument functions for effective fluorescence-enhanced optical tomography using FDPM.

Frequency domain photon migration (FDPM) measurements were conducted to assess particle interactions of concentrated, monodisperse, polystyrene samples obtained directly from industry by using multiple scattering light. The angle-integrated static structure factor, S(q), and static structure factor at small wave vector q, S(0), were obtained from FDPM measurements at high volume fractions ranging from 0.05 to 0.3, and were compared with those obtained from the monodisperse hard sphere Percus-Yevick (HSPY) model. Effects of different colloid sizes on structure factor evaluated at two different wavelengths were also investigated. Results show that the monodisperse HSPY model is suitable for accounting for particle interactions and local microstructures in these colloidal suspensions. Upon using the HSPY model, particle sizes of polystyrene suspensions were recovered from FDPM measurements at high volume fractions (up to 0.3), which agree well with the DLS measurement of diluted sample ( approximately 0.001). The study of polydispersity effect shows that the FDPM method can be successfully used for recovering the mean particle size of polydisperse colloidal suspension with low polydispersity (<16%) under the assumption of monodisperse hard sphere systems.

We present three-dimensional tomographic images of the absorption coefficient that is due to the presence of a fluorophore reconstructed from frequency domain fluence measurements of a tissue phantom containing a single, fluorescence contrast-enhanced inclusion. We show that such a reconstruction may be improved when the importance of measurement error correlations between relative phase shift and amplitude is assessed and when measurements are preprocessed to reduce the magnitude and the bias of system error.

The development of near-infrared (NIR) optical imaging for biomedical optical imaging is hampered by the computational intensiveness of large-scale three-dimensional (3-D) image reconstruction and the potential lack of endogenous contrast for detection of relevant tissue features. In this contribution the inverse optical imaging problem is formulated in three dimensions in a noncompressive geometry as a simple-bound constrained minimization problem in order to recover the interior fluorescence properties of exogenous contrast agent from frequency-domain photon migration measurements at the boundary. The solution of the forward optical diffusion problem for the frustum shape containing fluorescence inclusions of 10:1 contrast is accomplished by use of the Galerkin finite-element formulation. The inverse approach employs the truncated Newton method with trust region and a modification of automatic reverse differentiation to speed the computation of the optimization problem. The image-reconstruction results confirm that the constrained minimization may offer a more logical approach for the 3-D optical imaging problem than unconstrained optimization.

Most current efforts in near-infrared optical tomography are effectively limited to two-dimensional reconstructions due to the computationally intensive nature of full three-dimensional (3-D) data inversion. Previously, we described a new computationally efficient and statistically powerful inversion method APPRIZE (automatic progressive parameter-reducing inverse zonation and estimation). The APPRIZE method computes minimum-variance estimates of parameter values (here, spatially variant absorption due to a fluorescent contrast agent) and covariance, while simultaneously estimating the number of parameters needed as well as the size, shape, and location of the spatial regions that correspond to those parameters. Estimates of measurement and model error are explicitly incorporated into the procedure and implicitly regularize the inversion in a physically based manner. The optimal estimation of parameters is bounds-constrained, precluding infeasible values. In this paper, the APPRIZE method for optical imaging is extended for application to arbitrarily large 3-D domains through the use of domain decomposition. The effect of subdomain size on the performance of the method is examined by assessing the sensitivity for identifying 112 randomly located single-voxel heterogeneities in 58 3-D domains. Also investigated are the effects of unmodeled heterogeneity in background optical properties. The method is tested on simulated frequency-domain photon migration measurements at 100 MHz in order to recover absorption maps owing to fluorescent contrast agent. This study provides a new approach for computationally tractable 3-D optical tomography.

The use of near-infrared (NIR) light to interrogate deep tissues has enormous potential for molecular-based imaging when coupled with NIR excitable dyes. More than a decade has now passed since the initial proposals for NIR optical tomography for breast cancer screening using time-dependent measurements of light propagation in the breast. Much accomplishment in the development of optical mammography has been demonstrated, most recently in the application of time-domain, frequency-domain, and continuous-wave measurements that depend on endogenous contrast owing to angiogenesis and increased hemoglobin absorbance for contrast. Although exciting and promising, the necessity of angiogenesis-mediated absorption contrast for diagnostic optical mammography minimizes the potential for using NIR techniques to assess sentinel lymph node staging, metastatic spread, and multifocality of breast disease, among other applications. In this review, we summarize the progress made in the development of optical mammography, and focus on the emerging work underway in the use of diagnostic contrast agents for the molecular-based, diagnostic imaging of breast.

In the past, nonlinear unconstrained optimization of the optical imaging problem has focused on Newton-Raphson techniques. Besides requiring expensive computation of the Jacobian, the unconstrained minimization with Tikhonov regularization can pose significant storage problems for large-scale reconstructions, involving a large number of unknowns necessary for realization of optical imaging. We formulate the inverse optical imaging problem as both simple-bound constrained and unconstrained minimization problems in order to illustrate the reduction in computational time and storage associated with constrained image reconstructions. The forward simulator of excitation and generated fluorescence, consisting of the Galerkin finite-element formulation, is used in an inverse algorithm to find the spatial distribution of absorption and lifetime that minimizes the difference between predicted and synthetic frequency-domain measurements. The inverse approach employs the truncated Newton method with trust region and a modification of automatic reverse differentiation to speed the computation of the optimization problem. The reconstruction results confirm that the physically based, constrained minimization with efficient optimization schemes may offer a more logical approach to the large-scale optical imaging problem than unconstrained minimization with regularization.

We present in vivo fluorescent, near-infrared (NIR), reflectance images of indocyanine green (ICG) and carotene-conjugated 2-devinyl-2-(1-hexyloxyethyl) pyropheophorbide (HPPH-car) to discriminate spontaneous canine adenocarcinoma from normal mammary tissue. Following intravenous administration of 1.0 mg kg-1 ICG or 0.3 mg kg-1 HPPH-car into the canine, a 25 mW, 778 nm or 70 mW, 660 nm laser diode beam, expanded by a diverging lens to approximately 4 cm in diameter, illuminated the surface of the mammary tissue. Successfully propagating to the tissue surface, ICG or HPPH-car fluorescence generated from within the tissue was collected by an image-intensified, charge-coupled device camera fitted with an 830 or 710 nm bandpass interference filter. Upon collecting time-dependent fluorescence images at the tissue surface overlying both normal and diseased tissue volumes, and fitting these images to a pharmacokinetic model describing the uptake (wash-in) and release (wash-out) of fluorescent dye, the pharmacokinetics of fluorescent dye was spatially determined. Mapping the fluorescence intensity owing to ICG indicates that the dye acts as a blood pool or blood persistent agent, for the model parameters show no difference in the ICG uptake rates between normal and diseased tissue regions. The wash-out of ICG was delayed for up to 72 h after intravenous injection in tissue volumes associated with disease, because ICG fluorescence was still detected in the diseased tissue 72 h after injection. In contrast, HPPH-car pharmacokinetics illustrated active uptake into diseased tissues, perhaps owing to the overexpression of LDL receptors associated with the malignant cells. HPPH-car fluorescence was not discernable after 24 h. This work illustrates the ability to monitor the pharmacokinetic delivery of NIR fluorescent dyes within tissue volumes as great as 0.5-1 cm from the tissue surface in order to differentiate normal from diseased tissue volumes on the basis of parameters obtained from the pharmacokinetic models.

Dual wavelength frequency-domain measurements of photon migration (FDPM) are conducted on filtrate samples obtained from an industrial centrifugation process designed to separate Escherichia coli cell debris from the inclusion bodies. FDPM measurements consist of detecting phase delay of intensity-modulated light at 670 and 820 (or 830) nm. Optical properties of isotropic scattering and absorption are obtained from the regression of phase delay data to the optical diffusion equation. We show that the corresponding intensity-based measurements alone cannot provide accurate and independent estimates for these optical properties. However, FDPM-derived scattering coefficients of filtrate solutions (primarily consisting of 0.1-0.2 micrometer E. coli cell debris) are sensitive to approximately 1 vol % of added inclusion bodies (of 1-2 micrometer size). The technique, theory, and future adaptation of FDPM as an on-line monitor to detect the loss of inclusion bodies in centrifugation following homogenization are presented and contrasted to conventional, intensity-based measurements.

The measurement and analysis of frequency-domain photon migration (FDPM) measurements of powder absorbance in pharmaceutical powders is described in the context of other optical techniques. FDPM consists of launching intensity-modulated light into a powder and detecting the phase delay and amplitude modulation of the re-emitted light as a function of the modulation frequency. From analysis of the data using the diffusion approximation to the radiative transport equation, the absorption coefficient can be obtained. Absorption coefficient measurements of riboflavin in lactose mixtures are presented at concentrations of 0.1 to 1% (w/w) at near-infrared wavelengths where solution absorption cross sections are difficult to accurately measure using traditional transmission measurements in nonscattering solutions. FDPM measurements in powders enabled determinations of absorption coefficients that increase linearly with concentration (w/w) according to Beer-Lambert relationship. The extension of FDPM for monitoring absorbance of low-dose and ultralow-dose powder blending operations is presented.

A method is presented to determine fluorescence decay lifetimes within tissuelike scattering media. Fluorescence lifetimes are determined for micromolar concentrations of the dyes 3,3'-Diethylthiatricarbocyanine Iodide and Indocyanine Green by frequency-domain investigations of light propagating in turbid media. Dual-wavelength photon-migration measurements that use intensity-modulated sources at excitation and emission wavelengths of the fluorophores provide optical parameters of the media as well as fluorescence properties of the dyes. The deduction of fluorescence lifetimes requires no calibration with reference fluorophores, and the results are shown to be independent of dye concentration.

We present near-infrared frequency-domain photon migration imaging for the lifetime sensitive detection and localization of exogenous fluorescent contrast agents within tissue-simulating phantoms and actual tissues. We employ intensity-modulated excitation light that is expanded and delivered to the surface of a tissue or tissue-simulating phantom. The intensity-modulated fluorescence generated from within the volume propagates to the surface and is collected using a gain-modulated image-intensified charge-coupled device camera. From the spatial values of modulation amplitude and phase of the detected fluorescent light, micromolar volumes of diethylthiatricarbocyanine iodide (tau = 1.17 ns) and indocyanine green (ICG) (tau = 0.58 ns) embedded 1.0 cm deep in a tissue phantom are localized and discriminated on the basis of their lifetime differences. To demonstrate the utility of frequency-domain fluorescent measurements for imaging disease, we image the fluorescence emitted from the surface of in vivo and ex vivo canine mammary gland tissues containing lesions with preferential uptake of ICG. Pathology confirms the ability to detect spontaneous mammary tumors and regional lymph nodes amidst normal mammary tissue and fat as deep as 1.5 cm from the tissue surface.

Stochastic reconstruction techniques are developed for mapping the interior optical properties of tissues from exterior frequency-domain photon migration measurements at the air-tissue interface. Parameter fields of absorption cross section, fluorescence lifetime, and quantum efficiency are accurately reconstructed from simulated noisy measurements of phase shift and amplitude modulation by use of a recursive, Bayesian, minimum-variance estimator known as the approximate extended Kalman filter. Parameter field updates are followed by data-driven zonation to improve the accuracy, stability, and computational efficiency of the method by moving the system from an underdetermined toward an overdetermined set of equations. These methods were originally developed by Eppstein and Dougherty [Water Resources Res. 32, 3321 (1996)] for applications in geohydrology. Estimates are constrained to within feasible ranges by modeling of parameters as beta-distributed random variables. No arbitrary smoothing, regularization, or interpolation is required. Results are compared with those determined by use of Newton-Raphson-based inversions. The speed and accuracy of these preliminary Bayesian reconstructions suggest the near-future application of this inversion technology to three-dimensional biomedical imaging with frequency-domain photon migration.

Time-dependent measurements of light propagation were conducted in aqueous dispersions of 523 nm diameter polystyrene at concentrations between 0.1 and 0.4 solids volume fraction in order to assess how particle correlation is influenced by depletion interactions arising from the addition of soluble polyethyleneoxide (PEO). In the absence of polymer, the transport scattering length can be predicted from Mie scattering theory and the Percus-Yevick (P-Y) model for static structure of a dense hard-sphere colloidal solution. Depletion forces arising from the addition of PEO of varying molecular weights influenced the spatial ordering of the dispersion and caused a further increase in the transport scattering length beyond that predicted by hard-sphere static structure factor but similar to that predicted by the mean sphere approximation (MSA) to the P-Y model described by Ye et al. (1996). Onset of flocculation occurred with increased PEO addition and correlated with PEO molecular weight. Phase separation was noted by no further change in the transport scattering length, except when flocculation was induced by the highest molecular weight PEO. The use of time-dependent measurements of light propagation in dense systems provides an alternative to small-angle light, neutron, and X-ray scattering characterization of interaction potentials in dense, multiply scattering samples and promises further fruitful investigation of colloidal particle interactions in suspensions. Copyright 1999 Academic Press.

Frequency-domain measurements of photon migration are coupled with a model of fluorescence generation and propagation in order to develop a method for reconstructing maps of fluorescent properties within interior tissue volumes from exterior measurements at the air-tissue interface. Simulation results confirm the feasibility of optical imaging through the use of exogenously administered contrast agents on the basis of fluorophore decay kinetics and yield. Experimental measurements using single-pixel and multipixel devices illustrate that the contrast owing to exogenous fluorescence exceeds that owing to absorption or scattering caused by endogenous chromophores or tissue structure and owing to absorption caused by exogenous contrast agents.

From analytical and numerical solutions that predict the scattering of diffuse photon density waves and from experimental measurements of changes in phase shift theta and ac amplitude demodulation M caused by the presence of single and double cylindrical heterogeneities, we show that second- and higher-order perturbations can affect the prediction of the propagation characteristics of diffuse photon density waves. Our experimental results for perfect absorbers in a lossless medium suggest that the performance of fast inverse-imaging algorithms that use first-order Born or Rytov approximations might have inherent limitations compared with inverse solutions that use iterative solutions of a linear perturbation equation or numerical solutions of the diffusion equation.

The ability to map interior optical properties of a highly scattering medium from exterior measurements of light propagation is afforded by optical tomography. In this communication, we describe the problem of optical tomography, the techniques of photon migration measurements necessary to accomplish it, and the development of multipixel measurements for rapid collection of optical signals. These multipixel measurements are shown to provide detection of contrast owing to the optical properties of absorption and fluorescence associated with dye-laden heterogeneities embedded in a tissue-like scattering medium. From these rapid measurements, successful reconstruction of an interior optical property map may now be possible with clinically realistic data acquisition times. Applications for the technology arise for biomedical optical imaging for the in vivo detection of disease and the diagnosis of tissue (bio-) chemistry.

The ability to optically image or detect diseased tissue volumes located deep within tissues depends upon the degree of contrast provided by differences in local optical properties. In this report, we show that the exogenous contrast offered by fluorescent compounds is superior to that provided by nonfluorescing, light-absorbing compounds when time-dependent measurements are employed. In addition, we show that the induced contrast is not only moderated by the preferential uptake of fluorescent agents into diseased tissue volumes of interest but also by the fluorescent optical properties and the fluorescence dynamics in the specific tissue volume. Using tissue phantom studies, we demonstrated experimentally that near-infrared-absorbing and fluorescent dyes such as indocyanine green can provide detection of diseased tissue volumes from fluorescence measurements made at the periphery of tissue when there is perfect, 100-fold and 10-fold partitioning in diseased tissues over that in surrounding normal tissues. Experimental results of common laser dyes show the contrast is also mediated by the quantum yield and lifetime parameters that may be dependent upon the local tissue environment.

The feasibility of employing fluorescent contrast agents to perform optical imaging in tissues and other scattering media has been examined through computational studies. Fluorescence lifetime and yield can give crucial information about local metabolite concentrations or environmental conditions within tissues. This information can be employed toward disease detection, diagnosis, and treatment if noninvasively quantitated from reemitted optical signals. However, the problem of inverse image reconstruction of fluorescence yield and lifetime is complicated because of the highly scattering nature of the tissue. Here a light propagation model employing the diffusion equation is used to account for the scattering of both the excitation and fluorescent light. Simulated measurements of frequency-domain parameters of fluorescent modulated ac amplitude and phase lag are used as inputs to an inverse image-reconstruction algorithm, which employs the diffusion model to predict frequency-domain measurements resulting from a modulated input at the phantom periphery. In the inverse image-reconstruction algorithm, a Newton-Raphson technique combined with a Marquardt algorithm is employed to converge on the fluorescent properties within the medium. The successful reconstruction of both the fluorescence yield and lifetime in the case of a heterogeneous fluorophore distribution within a scattering medium has been demonstrated without a priori information or without the necessity of obtaining absence images.

Measurements of nanosecond and subnanosecond fluorescence lifetimes are restricted to dilute, nonscattering systems since excitation and emission photon times of flight significantly affect measured fluorescent decay kinetics. We provide the theoretical rationale for frequency-domain measurements of phase-shift and amplitude demodulation made at excitation and emission wavelengths for direct determination of lifetimes in tissues and other scattering media. We confirm our analytical expressions using standard laser dyes such as 3,3'-diethylthiatricarbocyanine iodide, IR- 125, and IR- 140 in polystyrene suspensions with similar scattering properties as tissues. Our results have significant implication for lifetime-based spectroscopy in tissues and other scattering media.

We have numerically solved the photon diffusion equation to predict the distribution of light in a tissue model system with a uniform concentration of fluorophore. Our results show that time-dependent measurements of light propagation can be used to monitor the fluorescent lifetimes of a uniformly distributed fluorophore in tissues. With proper referencing, frequency-domain measurements of phase-shift, theta, may allow quantitation of fluorescent lifetimes, tau, independent of changes in the local absorption and scattering properties. These results point to a new approach for noninvasive diagnostic monitoring through quantitation of fluorescent lifetime, tau, when the lifetime of the fluorophore is comparable with photon migration times.

We investigate the origin of fluorescent or phosphorescent signals reemitted from highly scattering media (such as tissues), using diffusion theory and probability analysis. Results show that the lifetime of a uniformly distributed phosphorescent or fluorescent optical probe will profoundly affect the volume interrogated by noninvasive reflectance measurements. When the lifetime is greater than photon migration times, the origin of the reemitted signal is confined closest to the surface. Our computations suggest that noninvasive measurements of tissue oxygen concentration may not necessarily interrogate deep tissues when systemically administered phosphorescent dyes are used.

In summary, Table II is a listing of the pitfalls and advantages of using Monte Carlo simulations and numerical solution of the diffusion equation to describe photon migration in tissues. Judicious use of these techniques to describe the solution to the forward imaging problem may allow determination of the best theoretical resolution and the smallest detectable volume for the range of optical property differences expected in situ or imposed by contrast agent administration. Furthermore, an understanding of the forward imaging problem through these numerical techniques also contributes to our understanding of the most efficient solution to the inverse imaging problem.