Tomographic image reconstruction from optical projections in light-diffusing media

Advanced readout methods for superheated emulsion detectors

Superheated emulsions develop visible vapor bubbles when exposed to ionizing radiation. They consist in droplets of a metastable liquid, emulsified in an inert matrix. The formation of a bubble cavity is accompanied by sound waves. Evaporated bubbles also exhibit a lower refractive index, compared to the inert gel matrix. These two physical phenomena have been exploited to count the number of evaporated bubbles and thus measure the interacting radiation flux. Systems based on piezoelectric transducers have been traditionally used to acquire the acoustic (pressure) signals generated by bubble evaporation. Such systems can operate at ambient noise levels exceeding 100 dB; however, they are affected by a significant dead time (>10 ms). An optical readout technique relying on the scattering of light by neutron-induced bubbles has been recently improved in order to minimize measurement dead time and ambient noise sensitivity. Beams of infra-red light from light-emitting diode (LED) sources cross the active area of the detector and are deflected by evaporated bubbles. The scattered light correlates with bubble density. Planar photodiodes are affixed along the detector length in optimized positions, allowing the detection of scattered light from the bubbles and minimizing the detection of direct light from the LEDs. A low-noise signal-conditioning stage has been designed and realized to amplify the current induced in the photodiodes by scattered light and to subtract the background signal due to intrinsic scattering within the detector matrix. The proposed amplification architecture maximizes the measurement signal-to-noise ratio, yielding a readout uncertainty of 6% (±1 SD), with 1000 evaporated bubbles in a detector active volume of 150 ml (6 cm detector diameter). In this work, we prove that the intensity of scattered light also relates to the bubble size, which can be controlled by applying an external pressure to the detector emulsion. This effect can be exploited during the readout procedure to minimize shadowing effects between bubbles, which become severe when the latter are several thousands. The detector we used in this work is based on superheated C-318 (octafluorocyclobutane), emulsified in 100 μm ± 10% (1 SD) diameter drops in an inert matrix of approximately 150 ml. The detector was operated at room temperature and ambient pressure.

Dimensionality Reduction Based Optimization Algorithm for Sparse 3-D Image Reconstruction in Diffuse Optical Tomography

Diffuse optical tomography (DOT) is a relatively low cost and portable imaging modality for reconstruction of optical properties in a highly scattering medium, such as human tissue. The inverse problem in DOT is highly ill-posed, making reconstruction of high-quality image a critical challenge. Because of the nature of sparsity in DOT, sparsity regularization has been utilized to achieve high-quality DOT reconstruction. However, conventional approaches using sparse optimization are computationally expensive and have no selection criteria to optimize the regularization parameter. In this paper, a novel algorithm, Dimensionality Reduction based Optimization for DOT (DRO-DOT), is proposed. It reduces the dimensionality of the inverse DOT problem by reducing the number of unknowns in two steps and thereby makes the overall process fast. First, it constructs a low resolution voxel basis based on the sensing-matrix properties to find an image support. Second, it reconstructs the sparse image inside this support. To compensate for the reduced sensitivity with increasing depth, depth compensation is incorporated in DRO-DOT. An efficient method to optimally select the regularization parameter is proposed for obtaining a high-quality DOT image. DRO-DOT is also able to reconstruct high-resolution images even with a limited number of optodes in a spatially limited imaging set-up.

Interstitial diffuse radiance spectroscopy of gold nanocages and nanorods in bulk muscle tissues

Radiance spectroscopy was applied to the interstitial detection of localized inclusions containing Au nanocages or nanorods with various concentrations embedded in porcine muscle phantoms. The radiance was quantified using a perturbation approach, which enabled the separation of contributions from the porcine phantom and the localized inclusion, with the inclusion serving as a perturbation probe of photon distributions in the turbid medium. Positioning the inclusion at various places in the phantom allowed for tracking of photons that originated from a light source, passed through the inclusion's location, and reached a detector. The inclusions with high extinction coefficients were able to absorb nearly all photons in the range of 650-900 nm, leading to a spectrally flat radiance signal. This signal could be converted to the relative density of photons incident on the inclusion. Finally, the experimentally measured quantities were expressed via the relative perturbation and arranged into the classical Beer-Lambert law that allowed one to extract the extinction coefficients of various types of Au nanoparticles in both the transmission and back reflection geometries. It was shown that the spatial variation of perturbation could be described as 1/r dependence, where r is the distance between the inclusion and the detector. Due to a larger absorption cross section, Au nanocages produced greater perturbations than Au nanorods of equal particle concentration, indicating a better suitability of Au nanocages as contrast agents for optical measurements in turbid media. Individual measurements from different inclusions were combined into detectability maps.

Spatial information in large-scale neural recordings

To record from a given neuron, a recording technology must be able to separate the activity of that neuron from the activity of its neighbors. Here, we develop a Fisher information based framework to determine the conditions under which this is feasible for a given technology. This framework combines measurable point spread functions with measurable noise distributions to produce theoretical bounds on the precision with which a recording technology can localize neural activities. If there is sufficient information to uniquely localize neural activities, then a technology will, from an information theoretic perspective, be able to record from these neurons. We (1) describe this framework, and (2) demonstrate its application in model experiments. This method generalizes to many recording devices that resolve objects in space and should be useful in the design of next-generation scalable neural recording systems.

Robust and accurate transient light transport decomposition via convolutional sparse coding

Ultrafast sources and detectors have been used to record the time-resolved scattering of light propagating through macroscopic scenes. In the context of computational imaging, decomposition of this transient light transport (TLT) is useful for applications, such as characterizing materials, imaging through diffuser layers, and relighting scenes dynamically. Here, we demonstrate a method of convolutional sparse coding to decompose TLT into direct reflections, inter-reflections, and subsurface scattering. The method relies on the sparsity composition of the time-resolved kernel. We show that it is robust and accurate to noise during the acquisition process.

Diffuse reflectance optical topography: location of inclusions in 3D and detectability limits

In the present contribution we investigate the images of CW diffusely reflected light for a point-like source, registered by a CCD camera imaging a turbid medium containing an absorbing lesion. We show that detection of μa variations (absorption anomalies) is achieved if images are normalized to background intensity. A theoretical analysis based on the diffusion approximation is presented to investigate the sensitivity and the limitations of our proposal and a novel procedure to find the location of the inclusions in 3D is given and tested. An analysis of the noise and its influence on the detection capabilities of our proposal is provided. Experimental results on phantoms are also given, supporting the proposed approach.

Estimating wide-angle, spatially varying reflectance using time-resolved inversion of backscattered light

Imaging through complex media is a well-known challenge, as scattering distorts a signal and invalidates imaging equations. For coherent imaging, the input field can be reconstructed using phase conjugation or knowledge of the complex transmission matrix. However, for incoherent light, wave interference methods are limited to small viewing angles. On the other hand, time-resolved methods do not rely on signal or object phase correlations, making them suitable for reconstructing wide-angle, larger-scale objects. Previously, a time-resolved technique was demonstrated for uniformly reflecting objects. Here, we generalize the technique to reconstruct the spatially varying reflectance of shapes hidden by angle-dependent diffuse layers. The technique is a noninvasive method of imaging three-dimensional objects without relying on coherence. For a given diffuser, ultrafast measurements are used in a convex optimization program to reconstruct a wide-angle, three-dimensional reflectance function. The method has potential use for biological imaging and material characterization.

Tagging photons with gold nanoparticles as localized absorbers in optical measurements in turbid media

We analyze a role of a localized inclusion as a probe for spatial distributions of migrating photons in turbid media. We present new experimental data and two-dimensional analysis of radiance detection of a localized absorptive inclusion formed by gold nanoparticles in Intralipid-1% when the target is translated along the line connecting the light source and detector. Data are analyzed using the novel analytical expression for the relative angular photon distribution function for radiance developed by extending the perturbation approach for fluence. Obtained photon maps allow predicting conditions for detectability of inclusions for which proximity to the detector is essential.

Radiance detection of non-scattering inclusions in turbid media

Detection of non-scattering domains (voids) is an area of active research in biomedical optics. To avoid complexities of image reconstruction algorithms and requirements of a priori knowledge of void locations inherent to diffuse optical tomography (DOT), it would be useful to establish specific experimental signatures of voids that would help identify and detect them by other means. To address this, we present a radiance-based spectro-angular mapping approach that identifies void locations in the angular domain and establishes their spectral features. Using water-filled capillaries in scattering Intralipid as a test platform, we demonstrate perturbations in the directional photon density distribution produced by individual voids.

Diffuser-aided diffuse optical imaging for breast tumor: a feasibility study based on time-resolved three-dimensional Monte Carlo modeling

This study proposed diffuser-aided diffuse optical imaging (DADOI) as a new approach to improve the performance of the conventional diffuse optical tomography (DOT) approach for breast imaging. The 3-D breast model for Monte Carlo simulation is remodeled from clinical MRI image. The modified Beer-Lambert's law is adopted with the DADOI approach to substitute the complex algorithms of inverse problem for mapping of spatial distribution, and the depth information is obtained based on the time-of-flight estimation. The simulation results demonstrate that the time-resolved Monte Carlo method can be capable of performing source-detector separations analysis. The dynamics of photon migration with various source-detector separations are analyzed for the characterization of breast tissue and estimation of optode arrangement. The source-detector separations should be less than 4 cm for breast imaging in DOT system. Meanwhile, the feasibility of DADOI was manifested in this study. In the results, DADOI approach can provide better imaging contrast and faster imaging than conventional DOT measurement. The DADOI approach possesses great potential to detect the breast tumor in early stage and chemotherapy monitoring that implies a good feasibility for clinical application.

Improvement of absorption and scattering discrimination by selection of sensitive points on temporal profile in diffuse optical tomography

We present a new method allowing the reconstruction of 3D time-domain diffuse optical tomography images, based on the time-dependent diffusion equation and an iterative algorithm (ART) using specific points on the temporal profiles. The first advantage of our method versus the full time-resolved scheme consists in considerably reducing the inverse problem resolution time. Secondly, in the presence of scattering heterogeneities, our method provides images of better quality comparatively to classical methods using full-time data or the first moments of the profiles.

Comparative study of tumor hypoxia by diffuse optical spectroscopy and immunohistochemistry in two tumor models

The capabilities of diffuse optical spectroscopy for noninvasive assessing of oxygen status in experimental tumors have been demonstrated. Specific features of the distribution of total hemoglobin, oxygenated hemoglobin, deoxygenated hemoglobin, and blood-oxygen saturation were shown on two tumor models having different histological structure and functional characteristics. The results obtained by the optical technique were verified by immunohistochemical study of tissue samples marked with exogenous marker of hypoxia--pimonidazole.

Fluorescence lifetime optical tomography with Discontinuous Galerkin discretisation scheme

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

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.

Transmission and fluorescence angular domain optical projection tomography of turbid media

When imaging through turbid media, objects are often blurred by scattered light. An optical collimator (i.e., an angular filter array) improves images by accepting only photons propagating within a narrow solid angle about the direction of the incident light. These photons are expected to participate in a limited number of small-angle scattering events, maintaining their original propagation direction and, finally, contributing to the development of a faithful image of an object within a turbid medium. The collimation method, also referred to as angular domain imaging (ADI), applies to a see-through configuration where the incident collimated light beam can be aligned with the collimator in a transillumination mode of operation. In this paper, we present angular domain optical projection tomography (ADOPT), a method that can extract depth information of optical contrast in turbid media with high longitudinal resolution based on ADI technology. The resolution of the ADI system has been tested over various depths in a 5 cm optical cuvette using a resolution target suspended in a homogeneous turbid medium. The ADOPT system reconstructed images from a series of angular domain projections collected at angular intervals. The system was used to measure the attenuation of an absorbing target in transmission mode (t-ADOPT) and to measure the light emitting from a fluorescent target (f-ADOPT). Tissue-mimicking phantoms were used to validate the performance of the method. In the t-ADOPT configuration, a background scattered light estimation and subtraction methodology was introduced to improve the imaging contrast. A target consisting of two graphite rods (0.9 mm diameter) was suspended in the cuvette by a rotation stage. An Indocyanine Green-filled glass rod was used as an imaging target in the f-ADOPT arrangement. The target was placed in a manner such that the line of laser light was perpendicular to the longitudinal axis of the rods. Several projections were collected at increments of 1.8 degrees and compiled into a sinogram. A transverse image was reconstructed from the sinogram by using filtered backprojection and image contrast was improved by experimental scatter measurements using a wedge prism and an image processing algorithm. The submillimeter target embedded in a 2 cm thick scattering medium (reduced scattering coefficient < or = 2.4 cm(-1)) was discernable in both the sinograms and the reconstructed images. In the f-ADOPT system, fluorescent line targets <1 cm in diameter embedded in a 2 cm thick scattering medium (reduced scattering coefficient < or = 0.8 cm(-1)) were discernable in both the sinograms and the reconstructed images. The proposed method could be used as the basis to construct an optical tomographic scanner for simultaneous absorption and fluorescence-based imaging of biological specimens (i.e., up to 7 mm across).

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.

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.

Near-infrared laser tomographic imaging with right-angled scattered coherent light using an optical heterodyne-detection-based confocal scanning system

We demonstrate, what is to the best of our knowledge, a novel optical tomographic method for the visualization of the inner structure of scattering media such as biological tissue in the near-infrared region. We constructed a scanning confocal imaging system with a cross-axes arrangement using optical fibers. This system is based on the optical heterodyne technique and enables the detection of very weak coherence photons that are generated in the spatially restricted confocal region and scattered laterally (90 degrees ) against an incident beam. To evaluate the fundamental imaging capabilities of the system, we assessed measurements from scattering phantoms composed of an Intralipid suspension with varying volume concentrations. The results of this study demonstrate that the right-angled scattered light adheres to the Lambert-Beer law and that the present system can detect light propagating through a distance of approximately 31l of the mean free path. An inclusion as small as 100 microm can be discriminated in a scattering media with an optical thickness of 4. We investigated the potential of the proposed system for imaging biological tissues in preliminary experiments using samples of chicken breast tissue.

Multispectral breast imaging using a ten-wavelength, 64 x 64 source/detector channels silicon photodiode-based diffuse optical tomography system

We describe a compact diffuse optical tomography system specifically designed for breast imaging. The system consists of 64 silicon photodiode detectors, 64 excitation points, and 10 diode lasers in the near-infrared region, allowing multispectral, three-dimensional optical imaging of breast tissue. We also detail the system performance and optimization through a calibration procedure. The system is evaluated using tissue-like phantom experiments and an in vivo clinic experiment. Quantitative two-dimensional (2D) and three-dimensional (3D) images of absorption and reduced scattering coefficients are obtained from these experiments. The ten-wavelength spectra of the extracted reduced scattering coefficient enable quantitative morphological images to be reconstructed with this system. From the in vivo clinic experiment, functional images including deoxyhemoglobin, oxyhemoglobin, and water concentration are recovered and tumors are detected with correct size and position compared with the mammography.

Time-resolved optical mammography using a liquid coupled interface

A method has been devised for generating three-dimensional optical images of the breast using a 32-channel time-resolved system and a liquid-coupled interface. The breast is placed in a hemispherical cup surrounded by sources and detectors, and the remaining space is filled with a fluid with tissue-like optical properties. This approach has three significant benefits. First, cups can accommodate a large range of breast sizes, enabling the entire volume of the breast to be sampled. Second, the coupling of the source and detector optics at the surface is constant and independent of the subject, enabling intensity measurements to be employed in the image reconstruction. Third, the external geometry of the reconstructed volume is known exactly. Images of isolated targets with contrasting absorbing and scattering properties have been acquired, and the performance of the system has been evaluated in terms of the contrast, spatial resolution, and localization accuracy. These parameters were strongly dependent on the location of the targets within the imaged volume. Preliminary images of a healthy human subject are also presented, which reveal subtle heterogeneity, particularly in the distribution of scatter. The ability to detect an absorbing target adjacent to the breast is also demonstrated.

Optical tomography of the breast using a multi-channel time-resolved imager

A time-resolved optical tomography system has been used to generate cross-sectional images of the human breast. Images are reconstructed using an iterative, nonlinear algorithm and measurements of mean photon flight time relative to those acquired on a homogeneous reference phantom. Thirty-eight studies have been performed on three healthy volunteers and 21 patients with a variety of breast lesions including cancer. We have successfully detected 17 out of 19 lesions, and shown that optical images of the healthy breast of the same volunteer display a heterogeneity which is repeatable over a period of months. However, results also indicate that the lack of accurate quantitation of optical parameters and limited morphological information limits the ability to characterize different types of lesions and distinguish benign from malignant tissues. Drawbacks of our current methodology and plans for overcoming them are discussed.

Multiple projection optical diffusion tomography with plane wave illumination

We describe a new data collection scheme for optical diffusion tomography in which plane wave illumination is combined with multiple projections in the slab imaging geometry. Multiple projection measurements are performed by rotating the slab around the sample. The advantage of the proposed method is that the measured data are more compatible with the dynamic range of most commonly used detectors. At the same time, multiple projections improve image quality by mutually interchanging the depth and transverse directions, and the scanned (detection) and integrated (illumination) surfaces. Inversion methods are derived for image reconstructions with extremely large data sets. Numerical simulations are performed for fixed and rotated slabs.

Monitoring recovery after laser surgery of the breast with optical tomography: a case study

Results are presented of a study to monitor the changes in the optical properties of breast tissue over a 12-month period after interstitial laser photocoagulation treatment of a fibroadenoma. The study involved generating cross-sectional images of the breast with a multichannel time-resolved imaging system and a nonlinear image reconstruction algorithm. Images of the internal absorbing and scattering properties revealed the expected initial inflammatory response, followed by the development of low-scattering cysts consistent with corresponding ultrasound examinations. Although results indicate that purely qualitative images can potentially provide clinically valuable data, means of enhancing diagnostic information by overcoming present limitations of the approach are discussed.

Diffuse optical tomography using intensity measurements and the a priori acquired regions of interest: theory and simulations

Light transmission data collected around an object show large variation with source-detector separation owing to the presence of single or multiple inhomogeneous regions in the object. This variation in the measured intensity is made use of to reconstruct regions of the inhomogeneous inclusions. In addition, it is possible to select a set of data from the above which is most likely least affected by the presence of the inhomogeneity, and estimate reasonably accurately the background optical properties from it. The reconstructed region is found to always contain the inhomogeneity and is of size approximately 140% by area of the inhomogeneity. With the regions to be reconstructed a priori known, a model-based iterative reconstruction procedure for reconstructing the optical properties of the region converged five times faster than without such information. It is also shown that whereas for the full object, a view-based propagation-backpropagation reconstruction procedure failed to converge, owing to large underdeterminacy of the problem, a smaller problem attempting to reconstruct a priori specified regions of interest converged and did so faster than a non-view-based approach for similar objects. Reconstruction results are presented from simulated transmitted intensity data from the following objects with regions of inhomogeneity in both absorption and scattering: (i) single centrally located inhomogeneity, (ii) two off-centred inhomogeneous regions of equal size and contrast (iii) two off-centred inhomogeneous regions of unequal size and equal contrast and (iv) two off-centred inhomogeneous regions of unequal size and contrast. Whereas the model-based iterative image reconstruction procedure gave good convergence in the first, second and third cases, in the fourth case the reconstructions failed to recover the exact numerical value of the optical properties in the higher contrast region.

Multislice Tomographic Imaging and Analysis of Human Breast-Equivalent Phantoms and Biological Tissues

A laser transillumination tomographic system, consisting of electrical, optical, mechanical, and software components, to obtain multislice images of tissue-equivalent breast phantoms and biological tissues, is developed. The tissue-equivalent phantoms are prepared from paraffin wax mixed with wax color pigments by matching their surface backscattered profiles as measured by multiprobe laser reflectometer, with that of respective tissues. The optical parameters of these phantoms are determined by matching their reflectance profiles with that as obtained by Monte Carlo simulation of optical scattering. For multislice tomographic analysis conical breast phantoms of height 80.0 mm and 80.0 mm base diameter with inclusions of different optical properties and dimensions are developed. The resolution of the inclusions in the tomograms depends on their sizes and optical parameters. The minimum size of the inclusion as detected by this procedure in a slice of diameter 50.0 mm is 3.0 mm. The structural variation as observed in the tomograms of phantoms of combination of biological tissues indicates its possible applications in detecting the abnormalities developing in human healthy soft tissues.

Dual-projection optical diffusion tomography

We propose a new approach to optical diffusion tomography that incorporates two orthogonal projections. All the data obtained in a double projection measurement are treated simultaneously. The second projection improves image quality due to the fact that the depth and transverse directions are interchanged. An image reconstruction algorithm is derived and illustrated with simulations. It is shown that the spatial resolution of images improves by a factor of 4-5 due to the second projection.

Parallel programming of gradient-based iterative image reconstruction schemes for optical tomography

Optical tomography (OT) is a fast developing novel imaging modality that uses near-infrared (NIR) light to obtain cross-sectional views of optical properties inside the human body. A major challenge remains the time-consuming, computational-intensive image reconstruction problem that converts NIR transmission measurements into cross-sectional images. To increase the speed of iterative image reconstruction schemes that are commonly applied for OT, we have developed and implemented several parallel algorithms on a cluster of workstations. Static process distribution as well as dynamic load balancing schemes suitable for heterogeneous clusters and varying machine performances are introduced and tested. The resulting algorithms are shown to accelerate the reconstruction process to various degrees, substantially reducing the computation times for clinically relevant problems.

Diffuse optical imaging of brain activation: approaches to optimizing image sensitivity, resolution, and accuracy

Near-infrared spectroscopy (NIRS) and diffuse optical imaging (DOI) are finding widespread application in the study of human brain activation, motivating further application-specific development of the technology. NIRS and DOI offer the potential to quantify changes in deoxyhemoglobin (HbR) and total hemoglobin (HbT) concentration, thus enabling distinction of oxygen consumption and blood flow changes during brain activation. While the techniques implemented presently provide important results for cognition and the neurosciences through their relative measures of HbR and HbT concentrations, there is much to be done to improve sensitivity, accuracy, and resolution. In this paper, we review the advances currently being made and issues to consider for improving optical image quality. These include the optimal selection of wavelengths to minimize random and systematic error propagation in the calculation of the hemoglobin concentrations, the filtering of systemic physiological signal clutter to improve sensitivity to the hemodynamic response to brain activation, the implementation of overlapping measurements to improve image spatial resolution and uniformity, and the utilization of spatial prior information from structural and functional MRI to reduce DOI partial volume error and improve image quantitative accuracy.

Time-series estimation of biological factors in optical diffusion tomography

We apply state space estimation techniques to the time-varying reconstruction problem in optical tomography. We develop a stochastic model for describing the evolution of quasi-sinusoidal medical signals such as the heartbeat, assuming these are represented as a known frequency with randomly varying amplitude and phase. We use the extended Kalman filter in combination with spatial regularization techniques to reconstruct images from highly under-determined time-series data. This system also naturally segments activity belonging to different biological processes. We present reconstructions of simulated data and of real data recorded from the human motor cortex (Franceschini et al 2000 Optics Express 6 49-57). It is argued that the application of these time-series techniques improves both the fidelity and temporal resolution of reconstruction in optical tomography.

Semi-three-dimensional algorithm for time-resolved diffuse optical tomography by use of the generalized pulse spectrum technique

Although a foil three-dimensional (3-D) reconstruction with both 3-D forward and inverse models provide, the optimal solution for diffuse optical tomography (DOT), because of the 3-D nature of photon diffusion in tissue, it is computationally costly for both memory requirement and execution time in a conventional computing environment. Thus in practice there is motivation to develop an image reconstruction algorithm with dimensional reduction based on some modeling approximations. Here we have implemented a semi-3-D modified generalized pulse spectrum technique for time-resolved DOT, where a two-dimensional (2-D) distribution of optical properties is approximately assumed, while we retain 3-D distribution of photon migration in tissue. We have validated the proposed algorithm by reconstructing 3-D structural test objects from both numerically simulated and experimental date. We demonstrate our algorithm by comparing it with the calibrated 2-D reconstruction that is in widespread use as a shortcut to 3-D imaging and proving that the semi-3-D algorithm outperforms the calibrated 2-D algorithm.

Absorption and scattering images of heterogeneous scattering media can be simultaneously reconstructed by use of dc data

We present a carefully designed phantom experimental study aimed to provide solid evidence that both absorption and scattering images of heterogeneous scattering media can be reconstructed independently from dc data. We also study the important absorption-scattering cross-talk issue. In this regard, we develop a simple normalizing scheme that is incorporated into our nonlinear finite-element-based reconstruction algorithm. Our results from the controlled phantom experiments show that the cross talk of an absorption object appearing in scattering images can be eliminated and that the cross talk of a scattering object appearing in absorption images can be reduced considerably. In addition, these carefully designed phantom experiments clearly suggest that both absorption and scattering images can be simultaneously recovered and quantitatively separated in highly scattering media by use of dc measurements. Finally, we discuss our results in light of recent theoretical findings on nonuniqueness for dc image reconstruction.

Statistical analysis of nonlinearly reconstructed near-infrared tomographic images: Part I--Theory and simulations

Near-infrared (NIR) diffuse tomography is an emerging method for imaging the interior of tissues to quantify concentrations of hemoglobin and exogenous chromophores non-invasively in vivo. It often exploits an optical diffusion model-based image reconstruction algorithm to estimate spatial property values from measurements of the light flux at the surface of the tissue. In this study, mean-squared error (MSE) over the image is used to evaluate methods for regularizing the ill-posed inverse image reconstruction problem in NIR tomography. Estimates of image bias and image standard deviation were calculated based upon 100 repeated reconstructions of a test image with randomly distributed noise added to the light flux measurements. It was observed that the bias error dominates at high regularization parameter values while variance dominates as the algorithm is allowed to approach the optimal solution. This optimum does not necessarily correspond to the minimum projection error solution, but typically requires further iteration with a decreasing regularization parameter to reach the lowest image error. Increasing measurement noise causes a need to constrain the minimum regularization parameter to higher values in order to achieve a minimum in the overall image MSE.

Application of the photon average trajectories method to real-time reconstruction of tissue inhomogeneities in diffuse optical tomography of strongly scattering media

The possibility of application of the photon average trajectories (PAT) method to real-time reconstruction of tissue inhomogeneities in diffuse optical tomography of strongly scattering media has been substantiated. By this method, the inverse problem is reduced to solution of the integral equation with integration along a conditional PAT. Such an approach allows the standard fast algebraic algorithms commonly used in projection computed tomography to be applied to diffuse optical image reconstruction. To demonstrate the capabilities of the PAT method, a numerical experiment on cross-sectional reconstruction of cylindrical strongly scattering objects with absorbing inhomogeneities has been done. Relative shadows caused by inhomogeneities are simulated via numerical solution of the non-stationary diffusion equation. To solve the inverse problem, the QR-factorization least-squares algorithm and the multiplicative algebraic reconstruction technique are used. The results are compared with those obtained by a well-known software package for temporal optical absorption and scattering tomography based on multiple solution of the diffusion equation. It is shown that the PAT method allows reconstruction of the optical structure of objects with comparable accuracy while saving reconstruction time considerably.

Imaging of in vitro chicken leg using time-resolved near-infrared optical tomography

Near-infrared optical imaging gains much attention because of its noninvasiveness and deep penetration depths into tissue. Here, we report near-infrared optical tomographic imaging of an in vitro chicken leg from time-resolved measurements. The in vitro chicken leg, dipped in a cylindrical container filled with diluted Intralipid-10% solution, was imaged with a multichannel time-resolved imaging system. A two-dimensional reconstruction algorithm based on a modified generalized pulse spectrum technique has been developed to reconstruct the images of both the absorption and reduced scattering coefficients simultaneously and quickly. The results demonstrate that a simultaneous reconstruction of absorption and reduced scattering coefficients from time-resolved measurement has a potential to reveal the changes in the optical properties associated with not only the physiological information but also the anatomical structure of the organ.

Improvement of image quality in diffuse optical tomography by use of full time-resolved data

In the field of diffuse optical tomography (DOT), it is widely accepted that time-resolved (TR) measurement can provide the richest information on photon migration in a turbid medium, such as biological tissue. However, the currently available image reconstruction algorithms for TR DOT are based mostly on the cw component or some featured data types of original temporal profiles, which are related to the solution of a time-independent diffusion equation. Although this methodology can greatly simplify the reconstruction process, it suffers from low spatial resolution and poor quantitativeness owing to the limitation of effectively applicable data types. To improve image quality, it has been argued that exploiting the full TR data is essential. We propose implementation of a DOT algorithm by using full TR data and furthermore a variant algorithm with time slices of TR data to alleviate the computational complexity and enhance noise robustness. Compared with those algorithms where the featured data types are used, our evaluations on the spatial resolution and quantitativeness show that a significant improvement in imaging quality can be achieved when full TR data are used, which convinces the DOT community of the potential advantage of the TR domain over cw and frequency domains.

Three-dimensional diffuse optical tomography of bones and joints

We present for the first time full three-dimensional (3D) volumetric reconstruction of absorption images of in vitro and in vivo bones and joints from near-infrared tomographic measurements. Imaging experiments were conducted on human finger and chicken bones embedded in cylindrical scattering media using a Clemson multichannel diffuse optical imager. The volumetric optical images were recovered with our 3D finite element based reconstruction algorithm. Our results show that 3D imaging methods can provide details of the joint structure/composition that would be impossible from two-dimensional imaging methods.

Three-dimensional optical tomographic imaging of breast in a human subject

We present for the first time a full three-dimensional (3-D) reconstruction of absorption images of breast from continuous-wave (cw) measurements performed on a premenopausal woman. Our 3-D optical images clearly reveal a large primary tumor as well as a small secondary tumor in a separate location of the breast. The multiple tumors identified by our 3-D optical imaging have been confirmed by the subsequent biopsy examination of the breast. Quantitative information of the optical images obtained is provided in terms of the location, size, and absorption coefficient of the tumors.

Photon path distribution and optical responses of turbid media: theoretical analysis based on the microscopic Beer-Lambert law

A concise theoretical treatment has been developed to describe the optical responses of a highly scattering inhomogeneous medium using functions of the photon path distribution (PPD). The treatment is based on the microscopic Beer-Lambert law and has been found to yield a complete set of optical responses by time- and frequency-domain measurements. The PPD is defined for possible photons having a total zigzag pathlength of l between the points of light input and detection. Such a distribution is independent of the absorption properties of the medium and can be uniquely determined for the medium under quantification. Therefore, the PPD can be calculated with an imaginary reference medium having the same optical properties as the medium under quantification except for the absence of absorption. One of the advantages of this method is that the optical responses, the total attenuation, the mean pathlength, etc are expressed by functions of the PPD and the absorption distribution.

Three-dimensional time-resolved optical tomography of a conical breast phantom

A 32-channel time-resolved imaging device for medical optical tomography has been employed to evaluate a scheme for imaging the human female breast. The fully automated instrument and the reconstruction procedure have been tested on a conical phantom with tissue-equivalent optical properties. The imaging protocol has been designed to obviate compression of the breast and the need for coupling fluids. Images are generated from experimental data with an iterative reconstruction algorithm that employs a three-dimensional (3D) finite-element diffusion-based forward model. Embedded regions with twice the background optical properties are revealed in separate 3D absorption and scattering images of the phantom. The implications for 3D time-resolved optical tomography of the breast are discussed.

Time resolved optical tomography of the human forearm

A 32-channel time-resolved optical imaging instrument has been developed principally to study functional parameters of the new-born infant brain. As a prelude to studies on infants, the device and image reconstruction methodology have been evaluated on the adult human forearm. Cross-sectional images were generated using time-resolved measurements of transmitted light at two wavelengths. All data were acquired using a fully automated computer-controlled protocol. Images representing the internal scattering and absorbing properties of the arm are presented, as well as images that reveal physiological changes during a simple finger flexion exercise. The results presented in this paper represent the first simultaneous tomographic reconstruction of the internal scattering and absorbing properties of a clinical subject using purely temporal data, with additional co-registered difference images showing repeatable absorption changes at two wavelengths in response to exercise.

Instrumentation and calibration protocol for imaging dynamic features in dense-scattering media by optical tomography

Instrumentation is described that is suitable for acquiring multisource, multidetector, time-series optical data at high sampling rates (up to 150 Hz) from tissues having arbitrary geometries. The design rationale, calibration protocol, and measured performance features are given for both a currently used, CCD-camera-based instrument and a new silicon-photodiode-based system under construction. Also shown are representative images that we reconstructed from data acquired in laboratory studies using the described CCD-based instrument.

Simultaneous mapping of absorption and scattering coefficients from a three-dimensional model of time-resolved optical tomography

A Newton-Raphson inversion algorithm has been extended for simultaneous absorption and scattering reconstruction of fully three-dimensional (3D) diffuse optical tomographic imaging from time-resolved measurements. The proposed algorithm is derived from the efficient computation of the Jacobian matrix of the forward model and uses either the algebraic reconstruction technique or truncated singular-value decomposition as the linear inversion tool. Its validation was examined with numerically simulated data from 3-D finite-element discretization models of tissuelike phantoms, with several combinations of geometric and optical properties, as well as two commonly used source-detector configurations. Our results show that the fully 3-D image reconstruction of an object can be achieved with reasonable quality when volumetric light propagation in tissues is considered, and temporal information from the measurements can be effectively employed. Also, we investigated the conditions under which 3-D issues could be approximately addressed with two-dimensional reconstruction algorithms and further demonstrated that these conditions are seldom predictable or attainable in practice. Thus the application of 3-D algorithms to realistic situations is necessary.

Quantitative optical image reconstruction of turbid media by use of direct-current measurements

We present a detailed experimental study concerning quantitative optical property reconstruction of heterogeneous turbid media by use of absolute dc data only. We performed experiments by using tissuelike phantoms in both single-target and multitarget configurations in which variations in target size and optical contrast with the background were explored. Our results show that both scattering and absorption images can be reconstructed quantitatively by use of dc data only, whereas it was impossible to obtain such quantitative information in previously reported studies. We believe that this improvement is primarily a result of the realization of a novel data preprocessing/optimization scheme for accurately determining several critical parameters needed for reconstruction. The use of this data preprocessing/optimization scheme also eliminates the calibration reference measurement previously required for reconstruction. Experimental confirmation of this scheme is given in detail.

Dual-spatial integration for longitudinal localization of inclusions in turbid media

We introduce a technique called dual-spatial integration (DSI) that is used to isolate and enhance inclusions that differ only by their longitudinal placement within a scattering medium. DSI uses three different source-detector configurations to section a scattering medium into three longitudinal zones. This sectioning permits the extraction of structures close to surfaces and the enhancement of those structures located in the central part of the medium. Both the simulation and the experimental results indicate that DSI has potential interest for applications in biomedical imaging such as optical mammography.

Developments toward diagnostic breast cancer imaging using near-infrared optical measurements and fluorescent contrast agents

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.