Fluorescence optical diffusion tomography using multiple-frequency data

Determining optical material parameters with motion in structured illumination

A set of power measurements as a function of controlled nanopositioner movement of a planar film arrangement in a standing wave field is presented as a means to obtain the thicknesses and the dielectric constants to a precision dictated by noise in an exciting laser beam and the positioning and detector process, all of which can be refined with averaging. From a mutual information perspective, knowing the set of positions at which measurements are performed adds information. While applicable to any arrangement of planar films, the implementation considered involves thin transmissive membranes, as are employed in applications such as optomechanics. We show that measured power data as a function of object position provides sensitivity to the film refractive index and far-subwavelength thickness. Use of a cost function allows iterative retrieval of the film parameters, and a multi-resolution framework is described as a computationally efficient procedure. The approach is complementary to ellipsometry and could play an important role in routine film characterization studies for fields involving solid state material processing, as is common in the semiconductor device field.

Improving Localization of Deep Inclusions in Time-Resolved Diffuse Optical Tomography

Time-resolved diffuse optical tomography is a technique used to recover the optical properties of an unknown diffusive medium by solving an ill-posed inverse problem. In time-domain, reconstructions based on datatypes are used for their computational efficiency. In practice, most used datatypes are temporal windows and Fourier transform. Nevertheless, neither theoretical nor numerical studies assessing different datatypes have been clearly expressed. In this paper, we propose an overview and a new process to compute efficiently a long set of temporal windows in order to perform diffuse optical tomography. We did a theoretical comparison of these large set of temporal windows. We also did simulations in a reflectance geometry with a spherical inclusion at different depths. The results are presented in terms of inclusion localization and its absorption coefficient recovery. We show that (1) the new windows computed with the developed method improve inclusion localization for inclusions at deep layers, (2) inclusion absorption quantification is improved at all depths and, (3) in some cases these windows can be equivalent to frequency based reconstruction at GHz order.

Recent methodology advances in fluorescence molecular tomography

Molecular imaging (MI) is a novel imaging discipline that has been continuously developed in recent years. It combines biochemistry, multimodal imaging, biomathematics, bioinformatics, cell & molecular physiology, biophysics, and pharmacology, and it provides a new technology platform for the early diagnosis and quantitative analysis of diseases, treatment monitoring and evaluation, and the development of comprehensive physiology. Fluorescence Molecular Tomography (FMT) is a type of optical imaging modality in MI that captures the three-dimensional distribution of fluorescence within a biological tissue generated by a specific molecule of fluorescent material within a biological tissue. Compared with other optical molecular imaging methods, FMT has the characteristics of high sensitivity, low cost, and safety and reliability. It has become the research frontier and research hotspot of optical molecular imaging technology. This paper took an overview of the recent methodology advances in FMT, mainly focused on the photon propagation model of FMT based on the radiative transfer equation (RTE), and the reconstruction problem solution consist of forward problem and inverse problem. We introduce the detailed technologies utilized in reconstruction of FMT. Finally, the challenges in FMT were discussed. This survey aims at summarizing current research hotspots in methodology of FMT, from which future research may benefit.

Optical property uncertainty estimates for spatial frequency domain imaging

Spatial frequency domain imaging (SFDI) is a wide-field diffuse optical imaging modality that has attracted considerable interest in recent years. Typically, diffuse reflectance measurements of spatially modulated light are used to quantify the optical absorption and reduced scattering coefficients of tissue, and with these, chromophore concentrations are extracted. However, uncertainties in estimated absorption and reduced scattering coefficients are rarely reported, and we know of no method capable of providing these when look-up table (LUT) algorithms are used to recover the optical properties. We present a method to generate optical property uncertainty estimates from knowledge of diffuse reflectance measurement errors. By employing the Cramér-Rao bound, we can quickly and efficiently explore theoretical SFDI performance as a function of spatial frequencies and sample optical properties, allowing us to optimize spatial frequency selection for a given application. In practice, we can also obtain useful uncertainty estimates for optical properties recovered with a two-frequency LUT algorithm, as we demonstrate with tissue-simulating phantom and in vivo experiments. Finally, we illustrate how absorption coefficient uncertainties can be propagated forward to yield uncertainties for chromophore concentrations, which could significantly impact the interpretation of experimental results.

Compactly Supported Radial Basis Function-Based Meshless Method for Photon Propagation Model of Fluorescence Molecular Tomography

Fluorescence Molecular Tomography (FMT) is a powerful imaging modality for the research of cancer diagnosis, disease treatment and drug discovery. Via three-dimensional (3-D) imaging reconstruction, it can quantitatively and noninvasively obtain the distribution of fluorescent probes in biological tissues. Currently, photon propagation of FMT is conventionally described by the Finite Element Method (FEM), and it can obtain acceptable image quality. However, there are still some inherent inadequacies in FEM, such as time consuming, discretization error and inflexibility in mesh generation, which partly limit its imaging accuracy. To further improve the solving accuracy of photon propagation model (PPM), we propose a novel compactly supported radial basis functions (CSRBFs)-based meshless method (MM) to implement the PPM of FMT. We introduced a series of independent nodes and continuous CSRBFs to interpolate the PPM, which can avoid complicated mesh generation. To analyze the performance of the proposed MM, we carried out numerical heterogeneous mouse simulation to validate the simulated surface fluorescent measurement. Then we performed an in vivo experiment to observe the tomographic reconstruction. The experimental results confirmed that our proposed MM could obtain more similar surface fluorescence measurement with the golden standard (Monte-Carlo method), and more accurate reconstruction result was achieved via MM in in vivo application.

Novel l 2,1-norm optimization method for fluorescence molecular tomography reconstruction

Fluorescence molecular tomography (FMT) is a promising tomographic method in preclinical research, which enables noninvasive real-time three-dimensional (3-D) visualization for in vivo studies. The ill-posedness of the FMT reconstruction problem is one of the many challenges in the studies of FMT. In this paper, we propose a l 2,1-norm optimization method using a priori information, mainly the structured sparsity of the fluorescent regions for FMT reconstruction. Compared to standard sparsity methods, the structured sparsity methods are often superior in reconstruction accuracy since the structured sparsity utilizes correlations or structures of the reconstructed image. To solve the problem effectively, the Nesterov's method was used to accelerate the computation. To evaluate the performance of the proposed l 2,1-norm method, numerical phantom experiments and in vivo mouse experiments are conducted. The results show that the proposed method not only achieves accurate and desirable fluorescent source reconstruction, but also demonstrates enhanced robustness to noise.

A Bayesian approach for characterization of soft tissue viscoelasticity in acoustic radiation force imaging

Biomechanical imaging techniques based on acoustic radiation force (ARF) have been developed to characterize the viscoelasticity of soft tissue by measuring the motion excited by ARF non-invasively. The unknown stress distribution in the region of excitation limits an accurate inverse characterization of soft tissue viscoelasticity, and single degree-of-freedom simplified models have been applied to solve the inverse problem approximately. In this study, the ARF-induced creep imaging is employed to estimate the time constant of a Voigt viscoelastic tissue model, and an inverse finite element (FE) characterization procedure based on a Bayesian formulation is presented. The Bayesian approach aims to estimate a reasonable quantification of the probability distributions of soft tissue mechanical properties in the presence of measurement noise and model parameter uncertainty. Gaussian process metamodeling is applied to provide a fast statistical approximation based on a small number of computationally expensive FE model runs. Numerical simulation results demonstrate that the Bayesian approach provides an efficient and practical estimation of the probability distributions of time constant in the ARF-induced creep imaging. In a comparison study with the single degree of freedom models, the Bayesian approach with FE models improves the estimation results even in the presence of large uncertainty levels of the model parameters.

Reconstruction of fluorescence molecular tomography via a nonmonotone spectral projected gradient pursuit method

Fluorescence molecular tomography (FMT) is a promising imaging technique in preclinical research, enabling three-dimensional location of the specific tumor position for small animal imaging. However, FMT presents a challenging inverse problem that is quite ill-posed and ill-conditioned. Thus, the reconstruction of FMT faces various challenges in its robustness and efficiency. We present an FMT reconstruction method based on nonmonotone spectral projected gradient pursuit (NSPGP) with /₁-norm optimization. At each iteration, a spectral gradient-projection method approximately minimizes a least-squares problem with an explicit one-norm constraint. A nonmonotone line search strategy is utilized to get the appropriate updating direction, which guarantees global convergence. Additionally, the Barzilai-Borwein step length is applied to build the optimal step length, further improving the convergence speed of the proposed method. Several numerical simulation studies, including multisource cases as well as comparative analyses, have been performed to evaluate the performance of the proposed method. The results indicate that the proposed NSPGP method is able to ensure the accuracy, robustness, and efficiency of FMT reconstruction. Furthermore, an in vivo experiment based on a heterogeneous mouse model was conducted, and the results demonstrated that the proposed method held the potential for practical applications of FMT.

On the use of the Cramér-Rao lower bound for diffuse optical imaging system design

We evaluated the potential of the Cramér-Rao lower bound (CRLB) to serve as a design metric for diffuse optical imaging systems. The CRLB defines the best achievable precision of any estimator for a given data model; it is often used in the statistical signal processing community for feasibility studies and system design. Computing the CRLB requires inverting the Fisher information matrix (FIM), however, which is usually ill-conditioned (and often underdetermined) in the case of diffuse optical tomography (DOT). We regularized the FIM by assuming that the inhomogeneity to be imaged was a point target and assessed the ability of point-target CRLBs to predict system performance in a typical DOT setting in silico. Our reconstructions, obtained with a common iterative algebraic technique, revealed that these bounds are not good predictors of imaging performance across different system configurations, even in a relative sense. This study demonstrates that agreement between the trends predicted by the CRLBs and imaging performance obtained with reconstruction algorithms that rely on a different regularization approach cannot be assumed a priori. Moreover, it underscores the importance of taking into account the intended regularization method when attempting to optimize source-detector configurations.

Fast and robust reconstruction for fluorescence molecular tomography via a sparsity adaptive subspace pursuit method

Fluorescence molecular tomography (FMT), as a promising imaging modality, can three-dimensionally locate the specific tumor position in small animals. However, it remains challenging for effective and robust reconstruction of fluorescent probe distribution in animals. In this paper, we present a novel method based on sparsity adaptive subspace pursuit (SASP) for FMT reconstruction. Some innovative strategies including subspace projection, the bottom-up sparsity adaptive approach, and backtracking technique are associated with the SASP method, which guarantees the accuracy, efficiency, and robustness for FMT reconstruction. Three numerical experiments based on a mouse-mimicking heterogeneous phantom have been performed to validate the feasibility of the SASP method. The results show that the proposed SASP method can achieve satisfactory source localization with a bias less than 1mm; the efficiency of the method is much faster than mainstream reconstruction methods; and this approach is robust even under quite ill-posed condition. Furthermore, we have applied this method to an in vivo mouse model, and the results demonstrate the feasibility of the practical FMT application with the SASP method.

Small animal fluorescence and bioluminescence tomography: a review of approaches, algorithms and technology update

Emerging fluorescence and bioluminescence tomography approaches have several common, yet several distinct features from established emission tomographies of PET and SPECT. Although both nuclear and optical imaging modalities involve counting of photons, nuclear imaging techniques collect the emitted high energy (100-511 keV) photons after radioactive decay of radionuclides while optical techniques count low-energy (1.5-4.1 eV) photons that are scattered and absorbed by tissues requiring models of light transport for quantitative image reconstruction. Fluorescence imaging has been recently translated into clinic demonstrating high sensitivity, modest tissue penetration depth, and fast, millisecond image acquisition times. As a consequence, the promise of quantitative optical tomography as a complement of small animal PET and SPECT remains high. In this review, we summarize the different instrumentation, methodological approaches and schema for inverse image reconstructions for optical tomography, including luminescence and fluorescence modalities, and comment on limitations and key technological advances needed for further discovery research and translation.

Small animal optical diffusion tomography with targeted fluorescence

Despite the broad impact in medicine that optics can bring, thus far practical approaches are limited to weak scatter or near-surface monitoring. We show a method that utilizes a laser topography scan and a diffusion equation model to describe the photon transport, together with a multiresolution unstructured grid solution to the nonlinear optimization measurement functional, that overcomes these limitations. We conclude that it is possible to achieve whole body optical imaging with a resolution suitable for finding cancer nodules within an organ during surgery, with the aid of a targeted imaging agent.

Low-cost three-dimensional imaging system combining fluorescence and ultrasound

In this paper, we present a dual-modality imaging system combining three-dimensional (3D) continuous-wave transillumination fluorescence tomography with 3D ultrasound (US) imaging. We validated the system with two phantoms, one containing fluorescent inclusions (Cy5.5) at different depths, and another varying-thickness semicylindrical phantom. Using raster scanning, the combined fluorescence/US system was used to collect the boundary fluorescent emission in the X-Y plane, as well as recovered the 3D surface and position of the inclusions from US signals. US images were segmented to provide soft priors for the fluorescence image reconstruction. Phantom results demonstrated that with priors derived from the US images, the fluorescent reconstruction quality was significantly improved. As further evaluation, we show pilot in vivo results using an Apo-E mouse to assess the feasibility and performance of this system in animal studies. Limitations and potential to be used in artherosclerosis studies are then discussed.

A photo-multiplier tube-based hybrid MRI and frequency domain fluorescence tomography system for small animal imaging

Fluorescence tomography (FT) is a promising molecular imaging technique that can spatially resolve both fluorophore concentration and lifetime parameters. However, recovered fluorophore parameters highly depend on the size and depth of the object due to the ill-posedness of the FT inverse problem. Structural a priori information from another high spatial resolution imaging modality has been demonstrated to significantly improve FT reconstruction accuracy. In this study, we have constructed a combined magnetic resonance imaging (MRI) and FT system for small animal imaging. A photo-multiplier tube is used as the detector to acquire frequency domain FT measurements. This is the first MR-compatible time-resolved FT system that can reconstruct both fluorescence concentration and lifetime maps simultaneously. The performance of the hybrid system is evaluated with phantom studies. Two different fluorophores, indocyanine green and 3-3' diethylthiatricarbocyanine iodide, which have similar excitation and emission spectra but different lifetimes, are utilized. The fluorescence concentration and lifetime maps are both reconstructed with and without the structural a priori information obtained from MRI for comparison. We show that the hybrid system can accurately recover both fluorescence intensity and lifetime within 10% error for two 4.2 mm-diameter cylindrical objects embedded in a 38 mm-diameter cylindrical phantom when MRI structural a priori information is utilized.

Cramer-Rao analysis of steady-state and time-domain fluorescence diffuse optical imaging

Using a Cramer-Rao analysis, we study the theoretical performances of a time and spatially resolved fDOT imaging system for jointly estimating the position and the concentration of a point-wide fluorescent volume in a diffusive sample. We show that the fluorescence lifetime is a critical parameter for the precision of the technique. A time resolved fDOT system that does not use spatial information is also considered. In certain cases, a simple steady-state configuration may be as efficient as this time resolved fDOT system.

Diffuse Optics for Tissue Monitoring and Tomography

This review describes the diffusion model for light transport in tissues and the medical applications of diffuse light. Diffuse optics is particularly useful for measurement of tissue hemodynamics, wherein quantitative assessment of oxy- and deoxy-hemoglobin concentrations and blood flow are desired. The theoretical basis for near-infrared or diffuse optical spectroscopy (NIRS or DOS, respectively) is developed, and the basic elements of diffuse optical tomography (DOT) are outlined. We also discuss diffuse correlation spectroscopy (DCS), a technique whereby temporal correlation functions of diffusing light are transported through tissue and are used to measure blood flow. Essential instrumentation is described, and representative brain and breast functional imaging and monitoring results illustrate the workings of these new tissue diagnostics.

Reconstruction of fluorescence/bioluminescence sources in biological medium with spatial filter

We propose a new method for reconstruction of emitting source distributions by use of a spatial filter and a successive updating process of the forward model for fluorescence/bioluminescence diffuse optical tomography. The spatial filter transforms a set of the measurement data to a single source strength at a position of interest, and the forward model is updated by use of the estimated source strengths. This updating process ignores the dispensable source positions from reconstruction according to the reconstructed source distribution, and the spatial resolution of the reconstructed image is improved. The estimated sources are also used for the reduction of artifacts induced by noises based on the singular value decomposition. Some numerical experiments show the advantages of the proposed method by comparing the present results with those obtained by the conventional methods of the least squares method nd Algebraic Reconstruction Technique. Finally the criteria for practical use of the method are quantitatively presented by the simulations for 2D and 3D geometries.

Illumination pattern optimization for fluorescence tomography: theory and simulation studies

Fluorescence molecular tomography is a powerful tool for 3D visualization of molecular targets and pathways in vivo in small animals. Owing to the high degrees of absorption and scattering of light through tissue, the fluorescence tomographic inverse problem is inherently ill-posed. In order to improve source localization and the conditioning of the light propagation model, multiple sets of data are acquired by illuminating the animal surface with different spatial patterns of near-infrared light. However, the choice of these patterns in most experimental setups is ad hoc and suboptimal. This paper presents a systematic approach for designing efficient illumination patterns for fluorescence tomography. Our objective here is to determine how to optimally illuminate the animal surface so as to maximize the information content in the acquired data. We achieve this by improving the conditioning of the Fisher information matrix. We parameterize the spatial illumination patterns and formulate our problem as a constrained optimization problem that, for a fixed number of illumination patterns, yields the optimal set of patterns. For geometric insight, we used our method to generate a set of three optimal patterns for an optically homogeneous, regular geometrical shape and observed expected symmetries in the result. We also generated a set of six optimal patterns for an optically homogeneous cuboidal phantom set up in the transillumination mode. Finally, we computed optimal illumination patterns for an optically inhomogeneous realistically shaped mouse atlas for different given numbers of patterns. The regularized pseudoinverse matrix, generated using the singular value decomposition, was employed to reconstruct the point spread function for each set of patterns in the presence of a sample fluorescent point source deep inside the mouse atlas. We have evaluated the performance of our method by examining the singular value spectra as well as plots of average spatial resolution versus estimator variance corresponding to different illumination schemes.

Multimodal optical molecular image reconstruction with frequency domain measurements

Multimodality molecular imaging is becoming more and more important to understand both the structural and the functional characteristics of tissue, organs and tumors. So far, invasive nuclear methods utilizing ionizing radiation have been the "gold standard" of molecular imaging. We investigate non-contact, non-invasive, patient-tolerant and inexpensive near infrared (NIR) frequency domain optical tomography (FDOT) as a functional complement to structural X-ray computed tomography (CT) data. We show a novel multifrequency NIR FDOT approach both in transmission and reflectance mode and employ radiative transport equation (RTE) for 3D reconstruction of a target with novel fluorescent gold nanoshell indocyanine green (NS ICG) in an ex vivo nude mouse. The results demonstrate that gold NS ICG with multifrequency NIR FDOT is a promising fluorophore for multimodal optical molecular image reconstruction.

Pre-clinical whole-body fluorescence imaging: Review of instruments, methods and applications

Fluorescence sampling of cellular function is widely used in all aspects of biology, allowing the visualization of cellular and sub-cellular biological processes with spatial resolutions in the range from nanometers up to centimeters. Imaging of fluorescence in vivo has become the most commonly used radiological tool in all pre-clinical work. In the last decade, full-body pre-clinical imaging systems have emerged with a wide range of utilities and niche application areas. The range of fluorescent probes that can be excited in the visible to near-infrared part of the electromagnetic spectrum continues to expand, with the most value for in vivo use being beyond the 630 nm wavelength, because the absorption of light sharply decreases. Whole-body in vivo fluorescence imaging has not yet reached a state of maturity that allows its routine use in the scope of large-scale pre-clinical studies. This is in part due to an incomplete understanding of what the actual fundamental capabilities and limitations of this imaging modality are. However, progress is continuously being made in research laboratories pushing the limits of the approach to consistently improve its performance in terms of spatial resolution, sensitivity and quantification. This paper reviews this imaging technology with a particular emphasis on its potential uses and limitations, the required instrumentation, and the possible imaging geometries and applications. A detailed account of the main commercially available systems is provided as well as some perspective relating to the future of the technology development. Although the vast majority of applications of in vivo small animal imaging are based on epi-illumination planar imaging, the future success of the method relies heavily on the design of novel imaging systems based on state-of-the-art optical technology used in conjunction with high spatial resolution structural modalities such as MRI, CT or ultrasound.

Three-dimensional reconstruction in free-space whole-body fluorescence tomography of mice using optically reconstructed surface and atlas anatomy

We present a 3-D image reconstruction method for free-space fluorescence tomography of mice using hybrid anatomical prior information. Specifically, we use an optically reconstructed surface of the experimental animal and a digital mouse atlas to approximate the anatomy of the animal as structural priors to assist image reconstruction. Experiments are carried out on a cadaver of a nude mouse with a fluorescent inclusion (2.4-mm-diam cylinder) implanted in the chest cavity. Tomographic fluorescence images are reconstructed using an iterative algorithm based on a finite element method. Coregistration of the fluorescence reconstruction and micro-CT (computed tomography) data acquired afterward show good localization accuracy (localization error 1.2+/-0.6 mm). Using the optically reconstructed surface, but without the atlas anatomy, image reconstruction fails to show the fluorescent inclusion correctly. The method demonstrates the utility of anatomical priors in support of free-space fluorescence tomography.

Direct reconstruction of pharmacokinetic-rate images of optical fluorophores from NIR measurements

In this paper, we present a new method to form pharmacokinetic-rate images of optical fluorophores directly from near infra-red (NIR) boundary measurements. We first derive a mapping from spatially resolved pharmacokinetic rates to NIR boundary measurements by combining compartmental modeling with a diffusion based NIR photon propagation model. We express this mapping as a state-space equation. Next, we introduce a spatio-temporal prior model for the pharmacokinetic-rate images and combine it with the state-space equation. We address the image formation problem using the extended Kalman filtering framework. We analyze the computational complexity of the resulting algorithms and evaluate their performance in numerical simulations. An important feature of our approach is that the reconstruction of fluorescence concentrations and compartmental modeling are combined into a single step 1) to take advantage of the inherent temporal correlations in dynamic NIR measurements, and 2) to incorporate spatio-temporal a priori information on pharmacokinetic-rate images. Simulation results show that the resulting algorithms are more robust and lead to higher signal-to-noise ratio as compared to existing approaches where the reconstruction of concentrations and compartmental modeling are treated separately. Additionally, we reconstructed pharmacokinetic-rate images using in vivo data obtained from three patients with breast tumors. The reconstruction results show that the pharmacokinetic rates of indocyanine green are higher inside the tumor region as compared to the surrounding tissue.

Towards in vivo imaging of intramolecular fluorescence resonance energy transfer parameters

Fluorescence resonance energy transfer (FRET) is a nonradiative energy transfer process based on dipole-dipole interaction between donor and acceptor fluorophores that are spatially separated by a distance of a few nanometers. FRET has proved to be of immense value in the study of cellular function and the underlying cause of disease due to, for example, protein misfolding (of consequence in Alzheimer's disease). The standard parameterization in intramolecular FRET is the lifetime and yield, which can be related to the donor-acceptor (DA) distance. FRET imaging has thus far been limited to in vitro or near-surface microscopy because of the deleterious effects of substantial scatter. We show that it is possible to extract the microscopic FRET parameters in a highly scattering environment by incorporating the FRET kinetics of an ensemble of DA molecules connected by a flexible or rigid linker into an optical diffusion tomography (ODT) framework. We demonstrate the efficacy of our approach for extracting the microscopic DA distance through simulations and an experiment using a phantom with scattering properties similar to tissue. Our method will allow the in vivo imaging of FRET parameters in deep tissue, and hence provide a new vehicle for the fundamental study of disease.

Fluorescence lifetime imaging from time resolved measurements using a shape-based approach

We present a novel fluorescent tomography algorithm to estimate the spatial distribution of fluorophores and the fluorescence lifetimes from surface time resolved measurements. The algorithm is a hybridization of the level set technique for recovering the distributions of distinct fluorescent markers with a gradient method for estimating their lifetimes. This imaging method offers several advantages compared to more traditional pixel-based techniques as, for example, well defined boundaries and a better resolution of the images. The numerical experiments show that our imaging method gives rise to accurate reconstructions in the presence of data noise and fluorescence background even for complicated fluorophore distributions in several-centimiter-thick biological tissue.

Performance dependence of hybrid x-ray computed tomography/fluorescence molecular tomography on the optical forward problem

Hybrid imaging systems combining x-ray computed tomography (CT) and fluorescence tomography can improve fluorescence imaging performance by incorporating anatomical x-ray CT information into the optical inversion problem. While the use of image priors has been investigated in the past, little is known about the optimal use of forward photon propagation models in hybrid optical systems. In this paper, we explore the impact on reconstruction accuracy of the use of propagation models of varying complexity, specifically in the context of these hybrid imaging systems where significant structural information is known a priori. Our results demonstrate that the use of generically known parameters provides near optimal performance, even when parameter mismatch remains.

Effects of sampling strategy on image quality in noncontact panoramic fluorescence diffuse optical tomography for small animal imaging

Fluorescence diffuse optical tomography is an emerging technology for molecular imaging with recent technological advances in biomarkers and photonics. The introduction of noncontact imaging methods enables very large-scale data acquisition that is orders of magnitude larger than that from earlier systems. In this study, the effects of sampling strategy on image quality were investigated using an imaging phantom mimicking small animals and further analyzed using singular value analysis (SVA). The sampling strategy was represented in terms of a number of key acquisition parameters, namely the numbers of sources, detectors, and imaging angles. A number of metrics were defined to quantitatively evaluate image quality. The effects of acquisition parameters on image quality were subsequently studied by varying each of the parameters within a reasonable range while maintaining the other parameters constant, a method analogue to partial derivative in mathematical analysis. It was found that image quality improves at a much slower rate if the acquisition parameters are above certain critical values (approximately 5 sources, approximately 15 detectors, and approximately 20 angles for our system). These critical values remain virtually the same even if other acquisition parameters are doubled. It was also found that increasing different acquisition parameters improves image quality with different efficiencies in terms of the number of measurements: for a system characterized by a smaller threshold in SVA (less than 10(-5) in our study), the number of sources is the most efficient, followed by the number of detectors and subsequently the number of imaging angles. However, for systems characterized by a larger threshold, the numbers of sources and angles are equally more efficient than the number of detectors.

In vivo fluorescence lifetime tomography

Local molecular and physiological processes can be imaged in vivo through perturbations in the fluorescence lifetime (FLT) of optical imaging agents. In addition to providing functional information, FLT methods can quantify specific molecular events and multiplex diagnostic and prognostic information. We have developed a fluorescence lifetime diffuse optical tomography (DOT) system for in vivo preclinical imaging. Data is captured using a time-resolved intensified charge coupled device (ICCD) system to measure fluorescence excitation and emission in the time domain. Data is then converted to the frequency domain, and we simultaneously reconstruct images of yield and lifetime using an extension to the normalized Born approach. By using differential phase measurements, we demonstrate DOT imaging of short lifetimes (from 350 ps) with high precision (+/-5 ps). Furthermore, this system retains the efficiency, speed, and flexibility of transmission geometry DOT. We demonstrate feasibility of FLT-DOT through a progressive series of experiments. Lifetime range and repeatability are first measured in phantoms. Imaging of subcutaneous implants then verifies the FLT-DOT approach in vivo in the presence of inhomogeneous optical properties. Use in a common research scenario is ultimately demonstrated by imaging accumulation of a targeted near-infrared (NIR) fluorescent-labeled peptide probe (cypate-RGD) in a mouse with a subcutaneous tumor.

Total light approach of time-domain fluorescence diffuse optical tomography

In this study, time-domain fluorescence diffuse optical tomography in biological tissue is numerically investigated using a total light approach. Total light is a summation of excitation light and zero-lifetime emission light divided by quantum yield. The zero-lifetime emission light is an emitted fluorescence light calculated by assuming that the fluorescence lifetime is zero. The zero-lifetime emission light is calculated by deconvolving the actually measured emission light with a lifetime function, an exponential function for fluorescence decay. The object for numerical simulation is a 2-D 10 mm-radius circle with the optical properties simulating biological tissues for near infrared light, and contains regions with fluorophore. The inverse problem of fluorescence diffuse optical tomography is solved using time-resolved simulated measurement data of the excitation and total lights for reconstructing the bsorption coefficient and fluorophore concentration simultaneously. The mean time of flight is used as the featured data-type extracted from the time-resolved data. The reconstructed images of fluorophore concentration show good quantitativeness and spatial reproducibility. By use of the total light approach, computation is performed much faster than the conventional ones.

Three-dimensional fluorescence-enhanced optical tomography using a hand-held probe based imaging system

Hand-held based optical imaging systems are a recent development towards diagnostic imaging of breast cancer. To date, all the hand-held based optical imagers are used to perform only surface mapping and target localization, but are not capable of demonstrating tomographic imaging. Herein, a novel hand-held probe based optical imager is developed towards three-dimensional (3-D) optical tomography studies. The unique features of this optical imager, which primarily consists of a hand-held probe and an intensified charge coupled device detector, are its ability to; (i) image large tissue areas (5 x 10 sq. cm) in a single scan, (ii) perform simultaneous multiple point illumination and collection, thus reducing the overall imaging time; and (iii) adapt to varying tissue curvatures, from a flexible probe head design. Experimental studies are performed in the frequency domain on large slab phantoms (approximately 650 ml) using fluorescence target(s) under perfect uptake (1:0) contrast ratios, and varying target depths (1-2 cm) and X-Y locations. The effect of implementing simultaneous over sequential multiple point illumination towards 3-D tomography is experimentally demonstrated. The feasibility of 3-D optical tomography studies has been demonstrated for the first time using a hand-held based optical imager. Preliminary fluorescence-enhanced optical tomography studies are able to reconstruct 0.45 ml target(s) located at different target depths (1-2 cm). However, the depth recovery was limited as the actual target depth increased, since only reflectance measurements were acquired. Extensive tomography studies are currently carried out to determine the resolution and performance limits of the imager on flat and curved phantoms.

Choice of data types in time resolved fluorescence enhanced diffuse optical tomography

In this paper we examine possible data types for time resolved fluorescence enhanced diffuse optical tomography (FDOT). FDOT is a particular case of diffuse optical tomography, where our goal is to analyze fluorophores deeply embedded in a turbid medium. We focus on the relative robustness of the different sets of data types to noise. We use an analytical model to generate the expected temporal point spread function (TPSF) and generate the data types from this. Varying levels of noise are applied to the TPSF before generating the data types. We show that local data types are more robust to noise than global data types, and should provide enhanced information to the inverse problem. We go on to show that with a simple reconstruction algorithm, depth and lifetime (the parameters of interest) of the fluorophore are better reconstructed using the local data types. Further we show that the relationship between depth and lifetime is better preserved for the local data types, suggesting they are in some way not only more robust, but also self-regularizing. We conclude that while the local data types may be more expensive to generate in the general case, they do offer clear advantages over the standard global data types.

A time domain fluorescence tomography system for small animal imaging

We describe the application of a time domain diffuse fluorescence tomography system for whole body small animal imaging. The key features of the system are the use of point excitation in free space using ultrashort laser pulses and noncontact detection using a gated, intensified charge-coupled device (CCD) camera. Mouse shaped epoxy phantoms, with embedded fluorescent inclusions, were used to verify the performance of a recently developed asymptotic lifetime-based tomography algorithm. The asymptotic algorithm is based on a multiexponential analysis of the decay portion of the data. The multiexponential model is shown to enable the use of a global analysis approach for a robust recovery of the lifetime components present within the imaging medium. The surface boundaries of the imaging volume were acquired using a photogrammetric camera integrated with the imaging system, and implemented in a Monte-Carlo model of photon propagation in tissue. The tomography results show that the asymptotic approach is able to separate axially located fluorescent inclusions centered at depths of 4 and 10 mm from the surface of the mouse phantom. The fluorescent inclusions had distinct lifetimes of 0.5 and 0.95 ns. The inclusions were nearly overlapping along the measurement axis and shown to be not resolvable using continuous wave (CW) methods. These results suggest the practical feasibility and advantages of a time domain approach for whole body small animal fluorescence molecular imaging, particularly with the use of lifetime as a contrast mechanism.

Three-dimensional scheme for time-domain fluorescence molecular tomography based on Laplace transforms with noise-robust factors

As a visualizing and quantitative method, Fluorescence Molecular Tomography (FMT) has many potential applications in biomedical field and its three-dimensional (3D) implementation is needed in both theory and practice. In this paper, we propose a 3D scheme for time-domain FMT within the normalized Born-ratio formulation. A finite element method solution to the Laplace transformed time-domain coupled diffusion equations is employed as the forward model, and the resultant linear inversions at two distinct transform-factors are solved with an algebraic reconstruction technique to separate fluorescent yield and lifetime images. The algorithm is validated using simulated data for 3D cylinder phantoms, and the spatial resolution and quantitativeness of the reconstruction assessed. We demonstrate that the proposed approach can accurately retrieve the positions and shapes of the targets with high spatial resolution and quantitative accuracy, and tolerate a signal-to-noise ratio down to 25dB by appropriately choosing the transform factors.

Fluorescence diffuse optical tomography with functional and anatomicala prioriinformation: feasibility study

Fluorescence diffuse optical tomography (FT) is an emerging molecular imaging technique that can spatially resolve both fluorophore concentration and lifetime parameters. In this study, we investigate the performance of a frequency-domain FT system for small inclusions that are embedded in a heterogeneous background. The results demonstrate that functional and structural a priori information is crucial to be able to recover both parameters with high accuracy. The functional a priori information is defined by the absorption and scattering maps at both excitation and emission wavelengths. Similarly, the boundaries of the small inclusion and different regions in the background are utilized as the structural a priori information. Without a priori information, the fluorophore concentration of a 5 mm inclusion in a 40 mm medium is recovered with 50% error, while the lifetime cannot be recovered at all. On the other hand, when both functional and structural information are available, the true lifetime can be recovered and the fluorophore concentration can be estimated only with 5% error. This study shows that a hybrid system that can acquire diffuse optical absorption tomography (DOT), FT and anatomical images in the same setting is essential to be able to recover the fluorophore concentration and lifetime accurately in vivo.

Design and development of a hand-held optical probe toward fluorescence diagnostic imaging

Near-infrared optical imaging is an emerging noninvasive technology toward breast cancer diagnosis. The optical imaging systems available to date are limited either by flexibility to image any given breast volume, patient comfort, or instrument portability. Here, a hand-held optical probe is designed and developed, 1. employing a unique measurement scheme of simultaneous multiple point illumination and collection for rapid data acquisition and minimal patient discomfort, and 2. employing a curved probe head such that it allows flexible imaging of tissue curvatures. Simulation studies are carried out on homogeneous slab phantoms (5x10x8 cc) to determine an appropriate source-detector configuration for the probe head. These design features are implemented in the development of the probe, which consisted of six simultaneous illuminating and 165 simultaneous collecting fibers, spaced 0.5 cm apart on a 5x10 sq-cm probe head. Simulation studies on 3-D slab and curved phantoms demonstrate an increase in the total area of predicted fluorescence amplitude and overall signal strength on using simultaneous multiple point sources over a single point source. The probe is designed and developed such that on coupling with a detection system in the future, the hand-held probe based imager can be clinically assessed toward cancer diagnostic imaging.

Optimal probing of optical contrast of breast lesions of different size located at different depths by US localization

We report a frequency domain optical tomography system utilizing three RF modulation frequencies, which are optimized for probing breast lesions of different size located at different depths. A real-time co-registered ultrasound scanner is used to provide on-site estimation of lesion size and location. Based on the lesion information, an optimal light modulation frequency can be selected, which may yield more accurate estimates of lesion angiogenesis and hypoxia. Phantom experiments have demonstrated that a high modulation frequency, such as 350Mhz, is preferable for probing small lesions closer to the surface while a low modulation frequency, such as 50Mhz, is desirable for imaging deeper and larger lesions. A clinical example of a large invasive carcinoma is presented to demonstrate the application of this novel technique.

Diffuse optical tomographic reconstruction using multifrequency data

We investigated the use of multifrequency diffuse optical tomography (MF-DOT) data for the reconstruction of the optical parameters. The experiments were performed in a 63 mm diameter cylindrical phantom containing a 15 mm diameter cylindrical object. Modulation frequencies ranging from 110 MHz to 280 MHz were used in the phantom experiments changing the contrast in absorption of the object with respect to the phantom while keeping the scattering value the same. The diffusion equation was solved using the finite element method. The sensitivity information from each frequency was combined to form a single Jacobian. The inverse problem was solved iteratively by minimizing the difference between the measurements and forward problem using single and multiple modulation frequency data. A multiparameter Tikhonov scheme was used for regularization. The phantom results show that the peak absorption coefficient in a region of interest was obtained with an error less then 5% using two-frequency reconstruction for absorption contrast values up to 2.2 times higher than background and 10% for the absorption contrast values larger than 2.2. The use of two-frequency data is sufficient to improve the quantitative accuracy compared with the single frequency reconstruction with appropriate selection of these frequencies.

Fluorescence molecular tomography in the presence of background fluorescence

Fluorescence molecular tomography is an emerging imaging technique that resolves the bio-distribution of engineered fluorescent probes developed for in vivo reporting of specific cellular and sub-cellular targets. The method can detect fluorochromes in picomole amounts or less, imaged through entire animals, but the detection sensitivity and imaging performance drop in the presence of background, non-specific fluorescence. In this study, we carried out a theoretical and an experimental investigation on the effect of background fluorescence on the measured signal and on the tomographic reconstruction. We further examined the performance of three subtraction methods based on physical models of photon propagation, using experimental data on phantoms and small animals. We show that the data pre-processing with subtraction schemes can improve image quality and quantification when non-specific background florescence is present.

Design and implementation of a multifrequency near-infrared diffuse optical tomography system

The design and implementation of a multifrequency and multispectral diffuse optical tomography system is described. Four wavelengths are utilized: 665, 785, 808, and 830 nm. The system is based on a network analyzer, which provides rf modulation signals for the laser diodes, as well as measures the amplitude and the phase of the detected signals. Six different modulation frequencies ranging from 110 to 280 MHz are used. The details of instrumentation, calibration, data acquisition, and performance of the system are given. A finite element algorithm is used to solve the diffusion equation, and an inverse solver based on this forward solver is implemented to calculate the absorption and scattering maps from the acquired data. Data acquisition for one wavelength is completed in less than 2.5 min for a single modulation frequency. The measurement repeatability is 0.5% in ac intensity and 0.2 deg in phase. The performance of the system is evaluated with phantom studies. A multifrequency reconstruction algorithm is used, in which a single absorption and scattering image pair is obtained using the whole dataset obtained at different modulation frequencies. It is shown that the multifrequency reconstruction approach provides superior image quality compared to the single frequency counterpart.

Estimation of kinetic model parameters in fluorescence optical diffusion tomography

We present a technique for reconstructing the spatially dependent dynamics of a fluorescent contrast agent in turbid media. The dynamic behavior is described by linear and nonlinear parameters of a compartmental model or some other model with a deterministic functional form. The method extends our previous work in fluorescence optical diffusion tomography by parametrically reconstructing the time-dependent fluorescent yield. The reconstruction uses a Bayesian framework and parametric iterative coordinate descent optimization, which is closely related to Gauss-Seidel methods. We demonstrate the method with a simulation study.

Validation of in vivo fluorochrome concentrations measured using fluorescence molecular tomography

Fluorescence molecular tomography (FMT) has emerged as a means of quantitatively imaging fluorescent molecular probes in three dimensions in living systems. To assess the accuracy of FMT in vivo, translucent plastic tubes containing a turbid solution with a known concentration of Cy 5.5 fluorescent dye are constructed and implanted subcutaneously in nude mice, simulating the presence of a tumor accumulating a fluorescent molecular reporter. Comparisons between measurements of fluorescent tubes made before and after implantation demonstrate that the accuracy of FMT reported for homogeneous phantoms extends to the in vivo situation. The sensitivity of FMT to background fluorescence is tested by imaging fluorescent tubes in mice injected with Cy 5.5-labeled Annexin V. For small tube fluorochrome concentrations, the presence of background fluorescence results in increases in the reconstructed concentration. This phenomenon is counteracted by applying a simple subtraction correction to the measured fluorescence data. The effects of varying tumor photon absorption are simulated by imaging fluorescent tubes with varying ink concentrations, and are found to be minor. These findings demonstrate the in vivo quantitative accuracy of fluorescence tomography, and encourage further development of this imaging modality as well as application of FMT in molecular imaging studies using fluorescent reporters.

Statistical approach for detection and localization of a fluorescing mouse tumor in Intralipid

We present a method for detecting and localizing a fluorescing tumor obscured underneath several millimeters of a multiply scattering, homogeneous medium from fluorescence measurements made above the surface. Using a statistical model of the measurement system, we develop approaches for detection by use of a binary hypothesis testing approach and localization by use of maximum-likelihood estimation. We also compute the probability of tumor detection and the Cramér-Rao lower bound for the localization estimate error, which are performance metrics that could potentially be optimized in an experimental design. We validate the methods in an experimental study involving an excised mouse tumor tagged with a new folate-indocyanine dye and obscured under a tissue-simulating lipid suspension.

Optical tomographic imaging of small animals

Diffuse optical tomography is emerging as a viable new biomedical imaging modality. Using visible and near-infrared light this technique can probe the absorption and scattering properties of biological tissues. The main applications are currently in brain, breast, limb and joint imaging; however, optical tomographic imaging of small animals is attracting increasing attention. This interest is fuelled by recent advances in the transgenic manipulation of small animals that has led to many models of human disease. In addition, an ever increasing number of optically reactive biochemical markers has become available, which allow diseases to be detected at the molecular level long before macroscopic symptoms appear. The past three years have seen an array of novel technological developments that have led to the first optical tomographic studies of small animals in the areas of cerebral ischemia and cancer.