Initial studies of in vivo absorbing and scattering heterogeneity in near-infrared tomographic breast imaging

Time-Gated Single-Photon Detection in Time-Domain Diffuse Optics: A Review

This work reviews physical concepts, technologies and applications of time-domain diffuse optics based on time-gated single-photon detection. This particular photon detection strategy is of the utmost importance in the diffuse optics field as it unleashes the full power of the time-domain approach by maximizing performances in terms of contrast produced by a localized perturbation inside the scattering medium, signal-to-noise ratio, measurement time and dynamic range, penetration depth and spatial resolution. The review covers 15 years of theoretical studies, technological progresses, proof of concepts and design of laboratory systems based on time-gated single-photon detection with also few hints on other fields where the time-gated detection strategy produced and will produce further impact.

High-density diffuse optical tomography for imaging human brain function

This review describes the unique opportunities and challenges for noninvasive optical mapping of human brain function. Diffuse optical methods offer safe, portable, and radiation free alternatives to traditional technologies like positron emission tomography or functional magnetic resonance imaging (fMRI). Recent developments in high-density diffuse optical tomography (HD-DOT) have demonstrated capabilities for mapping human cortical brain function over an extended field of view with image quality approaching that of fMRI. In this review, we cover fundamental principles of the diffusion of near infrared light in biological tissue. We discuss the challenges involved in the HD-DOT system design and implementation that must be overcome to acquire the signal-to-noise necessary to measure and locate brain function at the depth of the cortex. We discuss strategies for validation of the sensitivity, specificity, and reliability of HD-DOT acquired maps of cortical brain function. We then provide a brief overview of some clinical applications of HD-DOT. Though diffuse optical measurements of neurophysiology have existed for several decades, tremendous opportunity remains to advance optical imaging of brain function to address a crucial niche in basic and clinical neuroscience: that of bedside and minimally constrained high fidelity imaging of brain function.

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

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

Study of a simple model for the transition between the ballistic and the diffusive regimes in diffusive media

A Monte Carlo simulation was utilized to investigate a simple model for the transition between the ballistic and the diffusive regimes in diffusive media. The simulation focuses on the propagation of visible and near-infrared light in biological tissues. This research has mainly two findings: (1) the transition can be described, as was found experimentally, with good accuracy by only two terms (ballistic and diffusive). (2) The model can be utilized for cases where the absorption coefficient is not negligible compared to the scattering coefficient by adding a power-law prefactor to the diffusive term.

Region-of-interest diffuse optical tomography system

Diffuse optical tomography (DOT) using near-infrared light is a promising tool for non-invasive imaging of deep tissue. This technique is capable of quantitative reconstruction of absorption (μa) and scattering coefficient (μs) inhomogeneities in the tissue. The rationale for reconstructing the optical property map is that the absorption coefficient variation provides diagnostic information about metabolic and disease states of the tissue. The aim of DOT is to reconstruct the internal tissue cross section with good spatial resolution and contrast from noisy measurements non-invasively. We develop a region-of-interest scanning system based on DOT principles. Modulated light is injected into the phantom/tissue through one of the four light emitting diode sources. The light traversing through the tissue gets partially absorbed and scattered multiple times. The intensity and phase of the exiting light are measured using a set of photodetectors. The light transport through a tissue is diffusive in nature and is modeled using radiative transfer equation. However, a simplified model based on diffusion equation (DE) can be used if the system satisfies following conditions: (a) the optical parameter of the inhomogeneity is close to the optical property of the background, and (b) μs of the medium is much greater than μa (μs > > μa). The light transport through a highly scattering tissue satisfies both of these conditions. A discrete version of DE based on finite element method is used for solving the inverse problem. The depth of probing light inside the tissue depends on the wavelength of light, absorption, and scattering coefficients of the medium and the separation between the source and detector locations. Extensive simulation studies have been carried out and the results are validated using two sets of experimental measurements. The utility of the system can be further improved by using multiple wavelength light sources. In such a scheme, the spectroscopic variation of absorption coefficient in the tissue can be used to arrive at the oxygenation changes in the tissue.

Joint sparsity-driven non-iterative simultaneous reconstruction of absorption and scattering in diffuse optical tomography

Some optical properties of a highly scattering medium, such as tissue, can be reconstructed non-invasively by diffuse optical tomography (DOT). Since the inverse problem of DOT is severely ill-posed and nonlinear, iterative methods that update Green's function have been widely used to recover accurate optical parameters. However, recent research has shown that the joint sparse recovery principle can provide an important clue in achieving reconstructions without an iterative update of Green's function. One of the main limitations of the previous work is that it can only be applied to absorption parameter reconstruction. In this paper, we extended this theory to estimate the absorption and scattering parameters simultaneously when the background optical properties are known. The main idea for such an extension is that a joint sparse recovery step gives us unknown fluence on the estimated support set, which eliminates the nonlinearity in an integral equation for the simultaneous estimation of the optical parameters. Our numerical results show that the proposed algorithm reduces the cross-talk artifacts between the parameters and provides improved reconstruction results compared to existing methods.

In vivo quantitative imaging of normal and cancerous breast tissue using broadband diffuse optical tomography

PURPOSE

A NIR tomography system that combines frequency domain (FD) and continuous wave (CW) measurements was used to image normal and malignant breast tissues.

METHODS

FD acquisitions were confined to wavelengths less than 850 nm because of detector limitations, whereas light from longer wavelengths (up to 948 nm) was measured in CW mode with CCD-coupled spectrometer detection. The two data sets were combined and processed in a single spectrally constrained reconstruction to map concentrations of hemoglobin, water, and lipid, as well as scattering parameters in the breast.

RESULTS

Chromophore concentrations were imaged in the breasts of nine asymptomatic volunteers to evaluate their intrasubject and intersubject variability. Normal subject data showed physiologically expected trends. Images from three cancer patients indicate that the added CW data is critical to recovering the expected increases in water and decreases in lipid content within malignancies. Contrasts of 1.5 to twofold in hemoglobin and water values were found in cancers.

CONCLUSIONS

In vivo breast imaging with instrumentation that combines FD and CW NIR data acquisition in a single spectral reconstruction produces more accurate hemoglobin, water, and lipid results relative to FD data alone.

Globally accelerated reconstruction algorithm for diffusion tomography with continuous-wave source in an arbitrary convex shape domain

A new numerical imaging algorithm is presented for reconstruction of optical absorption coefficients from near-infrared light data with a continuous-wave source. As a continuation of our earlier efforts in developing a series of methods called "globally convergent reconstruction methods" [J. Opt. Soc. Am. A23, 2388 (2006)], this numerical algorithm solves the inverse problem through solution of a boundary-value problem for a Volterra-type integral partial differential equation. We deal here with the particular issues in solving the inverse problems in an arbitrary convex shape domain. It is demonstrated in numerical studies that this reconstruction technique is highly efficient and stable with respect to the complex distribution of actual unknown absorption coefficients. The method is particularly useful for reconstruction from a large data set obtained from a tissue or organ of particular shape, such as the prostate. Numerical reconstructions of a simulated prostate-shaped phantom with three different settings of absorption-inclusions are presented.

In vivo fluorescence imaging: a personal perspective

In vivo fluorescence imaging with near-infrared (NIR) light holds enormous potential for a wide variety of molecular diagnostic and therapeutic applications. Because of its quantitative sensitivity, inherent biological safety, and relative ease of use (i.e., with respect to cost, time, mobility, and its familiarity to a diverse population of investigators), fluorescence-based imaging techniques are being increasingly utilized in small-animal research. Moreover, there is substantial interest in the translation of novel optical techniques into the clinic, where they will prospectively aid in noninvasive and quantitative screening, disease diagnosis, and post-treatment monitoring of patients. Effective deep-tissue fluorescence imaging requires the application of exogenous NIR-emissive contrast agents. Currently, available probes fall into two major categories: organic and inorganic NIR fluorophores (NIRFs). In the studies reviewed herein, we utilized polymersomes (50 nm to 50 microm diameter polymer vesicles) for the incorporation and delivery of large numbers of highly emissive oligo (porphyrin)-based, organic NIRFs.

Implementation of a phase array diffuse optical tomographic imager

Diffuse optical tomography (DOT) using near-infrared (NIR) light is a promising tool for noninvasive imaging of deep tissue. This technique is capable of quantitative reconstructions of absorption coefficient inhomogeneities of tissue. The motivation for reconstructing the optical property variation is that it, and, in particular, the absorption coefficient variation, can be used to diagnose different metabolic and disease states of tissue. In DOT, like any other medical imaging modality, the aim is to produce a reconstruction with good spatial resolution and accuracy from noisy measurements. We study the performance of a phase array system for detection of optical inhomogeneities in tissue. The light transport through a tissue is diffusive in nature and can be modeled using diffusion equation if the optical parameters of the inhomogeneity are close to the optical properties of the background. The amplitude cancellation method that uses dual out-of-phase sources (phase array) can detect and locate small objects in turbid medium. The inverse problem is solved using model based iterative image reconstruction. Diffusion equation is solved using finite element method for providing the forward model for photon transport. The solution of the forward problem is used for computing the Jacobian and the simultaneous equation is solved using conjugate gradient search. The simulation studies have been carried out and the results show that a phase array system can resolve inhomogeneities with sizes of 5 mm when the absorption coefficient of the inhomogeneity is twice that of the background tissue. To validate this result, a prototype model for performing a dual-source system has been developed. Experiments are carried out by inserting an inhomogeneity of high optical absorption coefficient in an otherwise homogeneous phantom while keeping the scattering coefficient same. The high frequency (100 MHz) modulated dual out-of-phase laser source light is propagated through the phantom. The interference of these sources creates an amplitude null and a phase shift of 180 degrees along a plane between the two sources with a homogeneous object. A solid resin phantom with inhomogeneities simulating the tumor is used in our experiment. The amplitude and phase changes are found to be disturbed by the presence of the inhomogeneity in the object. The experimental data (amplitude and the phase measured at the detector) are used for reconstruction. The results show that the method is able to detect multiple inhomogeneities with sizes of 4 mm. The localization error for a 5 mm inhomogeneity is found to be approximately 1 mm.

A boundary element approach for image-guided near-infrared absorption and scatter estimation

Multimodality NIR spectroscopy systems offer the possibility of region-based vascular and molecular characterization of tissue in vivo. However, computationally efficient 3D image reconstruction algorithms specific to these image-guided systems currently do not exist. Image reconstruction is often based on finite-element methods (FEMs), which require volume discretization. Here, a boundary element method (BEM) is presented using only surface discretization to recover the optical properties in an image-guided setting. The reconstruction of optical properties using BEM was evaluated in a domain containing a 30 mm inclusion embedded in two layer media with different noise levels and initial estimates. For 5% noise in measurements, and background starting values for reconstruction, the optical properties were recovered to within a mean error of 6.8%. When compared with FEM for this case, BEM showed a 28% improvement in computational time. BEM was also applied to experimental data collected from a gelatin phantom with a 25 mm inclusion and could recover the true absorption to within 6% of expected values using less time for computation compared with FEM. When applied to a patient-specific breast mesh generated using MRI, with a 2 cm ductal carcinoma, BEM showed successful recovery of optical properties with less than 5% error in absorption and 1% error in scattering, using measurements with 1% noise. With simpler and faster meshing schemes required for surface grids as compared with volume grids, BEM offers a powerful and potentially more feasible alternative for high-resolution 3D image-guided NIR spectroscopy.

Developments in quantitative oxygen-saturation imaging of breast tissue in vivo using multispectral near-infrared tomography

Imaging of oxygen saturation provides a spatial map of the tissue metabolic activity and has potential in diagnosis and treatment monitoring of breast cancer. Oxygen-saturation imaging is possible through near-infrared (NIR) tomography, but has low signal-to-noise ratio (SNR). This can be augmented by using NIR tomography as an add-on to MRI. Presented are results from a free-standing NIR system and a hybrid MR-guided system for breast imaging. In results from imaging 60 healthy volunteers in the initial NIR system, oxygen saturation was a significant discriminator between the BIRADS classifications of adipose tissue, heterogeneously dense, and extremely dense tissue. By using the MR-guided NIR system, more accurate tissue-specific data were obtained on adipose and fibroglandular volumes, with 11 healthy volunteers. In these data, oxygen saturation in the adipose tissue correlated with percentage of adipose tissue. In two case studies of infiltrating ductal carcinomas, oxygen saturation was reduced at the site of the tumor, as compared with the surrounding healthy tissue, agreeing with conventional thought that hypoxia exists in larger solid tumors. The MRI-guided NIR images of oxygen saturation provide higher resolution and superior SNR and will likely be used in the future to study and characterize specific tissue volumes.

Noise-tolerance analysis for detection and reconstruction of absorbing inhomogeneities with diffuse optical tomography using single- and phase-correlated dual-source schemes

An iterative reconstruction procedure is used to invert intensity data from both single- and phase-correlated dual-source illuminations for absorption inhomogeneities. The Jacobian for the dual source is constructed by an algebraic addition of the Jacobians estimated for the two sources separately. By numerical simulations, it is shown that the dual-source scheme performs superior to the single-source system in regard to (i) noise tolerance in data and (ii) ability to reconstruct smaller and lower contrast objects. The quality of reconstructions from single-source data, as indicated by mean-square error at convergence, is markedly poorer compared to their dual-source counterpart, when noise in data was in excess of 2%. With fixed contrast and decreasing inhomogeneity diameter, our simulations showed that, for diameters below 7 mm, the dual-source scheme has a higher percentage contrast recovery compared to the single-source scheme. Similarly, the dual-source scheme reconstructs to a higher percentage contrast recovery from lower contrast inhomogeneity, in comparison to the single-source scheme.

Reconstruction method with data from a multiple-site continuous-wave source for three-dimensional optical tomography

A method is presented for reconstruction of the optical absorption coefficient from transmission near-infrared data with a cw source. As it is distinct from other available schemes such as optimization or Newton's iterative method, this method resolves the inverse problem by solving a boundary value problem for a Volterra-type integral-differential equation. It is demonstrated in numerical studies that this technique has a better than average stability with respect to the discrepancy between the initial guess and the actual unknown absorption coefficient. The method is particularly useful for reconstruction from a large data set obtained from a CCD camera. Several numerical reconstruction examples are presented.

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

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

Near-infrared characterization of breast tumors in vivo using spectrally-constrained reconstruction

Multi-wavelength Near-Infrared (NIR) Tomography was utilized in this study to non-invasively quantify physiological parameters of breast tumors using direct spectral reconstruction. Frequency domain NIR measurements were incorporated with a new spectrally constrained direct chromophore and scattering image reconstruction algorithm, which was validated in simulations and experimental phantoms. Images of total hemoglobin, oxygen saturation, water, and scatter parameters were obtained with higher accuracy than previously reported. Using this spectral approach, in vivo NIR images are presented and interpreted through a series of case studies (n=6 subjects) having differing abnormalities. The corresponding mammograms and ultrasound images are also evaluated. Three of six cases were malignant (infiltrating ductal carcinomas) and showed higher hemoglobin (34-86% increase), a reduction in oxygen saturation, an increase in water content as well as scatter changes relative to surrounding normal tissue. Three of six cases were benign, two of which were diagnosed with fibrocystic disease and showed a dominant contrast in water, consistent with fluid filled cysts. Scatter amplitude was the main source of contrast in the volunteer with the benign condition fibrosis, which typically contains denser collagen tissue. The changes monitored correspond to physiological changes associated with angiogenesis, hypoxia and cell proliferation anticipated in cancers. These changes represent potential diagnostic indicators, which can be assessed to characterize breast tumors.

Measurement of particle-size distribution and concentration in heterogeneous turbid media with multispectral diffuse optical tomography

We present a method that is capable of extracting particle-size distribution (PSD) and concentration in heterogeneous turbid media by use of multispectral diffuse optical tomography (MSDOT). After the spectroscopic scattering images of the heterogeneous turbid media are obtained with MSDOT, the morphologic information of particles in the heterogeneities is recovered with an iterative regularized reconstruction algorithm based on Mie scattering theory when a particular form of PSD is assumed (Gaussian distribution is used in this study). The method described is tested and evaluated with both simulated and experimental data. The simulations are intended to test the sensitivity of the overall approach to noise effect. A series of phantom experiments are conducted with our newly developed ten-wavelength MSDOT system. Polystyrene microsphere suspensions contain particles of varying size from 2 to 6 microm as targets are embedded in a scattering background medium in these experiments. To achieve optimized results from experimental data, we developed a data preprocessing method for MSDOT as well as a scheme for calibrating scattering spectra. The results from both simulations and experiments show that the particle mean size and concentration can be reconstructed with acceptable accuracy, whereas the recovery of the standard deviation is sensitive to noise effect and can be as large as 86% from the experimental data.

Non-PET functional imaging techniques: optical

The strong and steady development of diffuse optical spectroscopy and tomography as new biomedical optics technologies promises to bring these optical techniques into clinical practice. This article provides a brief review of the light-tissue interaction, the instrumentation, and the theory relevant to this field. This is followed by a survey of the three main applications: brain imaging, muscle imaging, and breast imaging. Lastly, the future outlook of the technology is presented, highlighting the new promises based on recent breakthroughs.

Fast photoacoustic imaging system based on 320-element linear transducer array

A fast photoacoustic (PA) imaging system, based on a 320-transducer linear array, was developed and tested on a tissue phantom. To reconstruct a test tomographic image, 64 time-domain PA signals were acquired from a tissue phantom with embedded light-absorption targets. A signal acquisition was accomplished by utilizing 11 phase-controlled sub-arrays, each consisting of four transducers. The results show that the system can rapidly map the optical absorption of a tissue phantom and effectively detect the embedded light-absorbing target. By utilizing the multi-element linear transducer array and phase-controlled imaging algorithm, we thus can acquire PA tomography more efficiently, compared to other existing technology and algorithms. The methodology and equipment thus provide a rapid and reliable approach to PA imaging that may have potential applications in noninvasive imaging and clinic diagnosis.

Improved quantification of small objects in near-infrared diffuse optical tomography

Diffuse optical tomography allows quantification of hemoglobin, oxygen saturation, and water in tissue, and the fidelity in this quantification is dependent on the accuracy of optical properties determined during image reconstruction. In this study, a three-step algorithm is proposed and validated that uses the standard Newton minimization with Levenberg-Marquardt regularization as the first step. The second step is a modification to the existing algorithm using a two-parameter regularization to allow lower damping in a region of interest as compared to background. This second stage allows the recovery of the actual size of an inclusion. A region-based reconstruction is the final third step, which uses the estimated size and position information from step 2 to yield quantitatively accurate average values for the optical parameters. The algorithm is tested on simulated and experimental data and is found to be insensitive to object contrast and position. The percentage error between the true and the average recovered value for the absorption coefficient in test images is reduced from 47 to 27% for a 10-mm inclusion, from 38 to 13% for a 15-mm anomaly, and from 28 to 5.5% for a 20-mm heterogeneity. Simulated data with absorbing and scattering heterogeneities of 15 mm diam located in different positions show recovery with less than 15% error in absorption and 6% error in reduced scattering coefficients. The algorithm is successfully applied to clinical data from a subject with a breast abnormality to yield quantitatively increased absorption coefficients, which enhances the contrast to 3.8 compared to 1.23 previously.

Endoscopically compatible near-infrared photon migration probe

We have developed a 2.3-mm-diameter fiber-optic probe for near-infrared photon migration spectroscopy that can be inserted into the body through an endoscope or biopsy needle. This probe is specifically designed to be inserted into a core biopsy needle to facilitate optical sampling of lesions during breast needle biopsy. This probe was tested on tissue phantoms containing heterogeneities (to stimulate breast lesions) of various sizes and optical properties. Under the conditions tested, the probe can measure the absorption coefficient to within 30% for heterogeneities with radii as small as 10 mm.

Method for reduction of background artifacts of images in scanning holography with a Fresnel-zone-plate coded aperture

Near-infrared scanning holography with a Fresnel zone plate (FZP) coded aperture has potential applications in imaging through turbid media. However, the nonnegative intensity-distribution function of the FZP coded aperture introduces the background artifacts into the reconstructed images, reducing the contrast and the signal-to-noise ratio (SNR) of the images. A novel method termed as the composite hologram is proposed to reduce the artifacts. The computer simulations showed that the contrast and the SNR of the reconstructed images had improvements of at least 50.2% and 5.58-dB, respectively, compared with the conventional method. The composite hologram of a metal ring with a 6.0-mm diameter made by a wire with a 0.4-mm diameter immersing in 1% intralipid solution was recorded, and the reconstruction was performed numerically. The experimental results demonstrated that the contrast and the SNR of the reconstructed image had improvements of at least 32.3% and 2.51-dB, respectively.

Diffuse optical tomography with physiological and spatiala prioriconstraints

Diffuse optical tomography is a typical inverse problem plagued by ill-condition. To overcome this drawback, regularization or constraining techniques are incorporated in the inverse formulation. In this work, we investigate the enhancement in recovering functional parameters by using physiological and spatial a priori constraints. More accurate recovery of the two main functional parameters that are the blood volume and the relative saturation is demonstrated through simulations by using our method compared to actual techniques.

Metabolism-enhanced tumor localization by fluorescence imaging: in vivo animal studies

We present a high-sensitivity near-infrared optical imaging system for noninvasive cancer detection and localization based on molecularly labeled fluorescent contrast agents. This frequency-domain system utilizes the interferencelike pattern of diffuse photon density waves to achieve high detection sensitivity and localization accuracy for the fluorescent heterogeneity embedded inside the scattering media. A two-dimensional localization map is obtained through reflectance probe geometry and goniometric reconstruction. In vivo measurements with a tumor-bearing mouse model by use of the novel Cypate-mono-2-deoxy-glucose fluorescent contrast agent, which targets the enhanced tumor glycolysis, demonstrate the feasibility of detection of a 2-cm-deep subsurface tumor in the tissuelike medium, with a localization accuracy within 2-3 mm.

Spectral ballistic imaging: a novel technique for viewing through turbid or obstructing media

We propose a new method for viewing through turbid or obstructing media. The medium is illuminated with a modulated cw laser and the amplitude and phase of the transmitted (or reflected) signal is measured. This process takes place for a set of wavelengths in a certain wide band. In this way we acquire the Fourier transform of the temporal output. With this information we can reconstruct the temporal shape of the transmitted signal by computing the inverse transform. The proposed method benefits from the advantages of the first-light technique: high resolution, simple algorithms, insensitivity to boundary condition, etc., without suffering from its main deficiencies: complex and expensive equipment.

In vivo continuous-wave optical breast imaging enhanced with Indocyanine Green

We investigate the uptake of a nontargeted contrast agent by breast tumors using a continuous wave diffuse optical tomography apparatus. The instrument operates in the near-infrared spectral window and employs 16 sources and 16 detectors to collect light in parallel on the surface of the tumor-bearing breast (coronal geometry). In our protocol an extrinsic contrast agent, Indocyanine Green (ICG), was injected by bolus. Three clinical scenarios with three different pathologies were investigated. A two-compartment model was used to analyze the pharmacokinetics of ICG and preprocess the data, and diffuse optical tomography was used for imaging. Localization and delineation of the tumor was achieved in good agreement with a priori information. Moreover, different dynamical features were observed for differing pathologies. The malignant cases exhibited slower rate constants (uptake and outflow) compared to healthy tissue. These results provide further evidence that in vivo pharmacokinetics of ICG in breast tumors may be a useful diagnostic tool for differentiation of benign and malignant pathologies.

Mesh-based enhancement schemes in diffuse optical tomography

Two mesh-based methods including dual meshing and adaptive meshing are developed to improve the finite element-based reconstruction of both absorption and scattering images of heterogeneous turbid media. The idea of dual meshing scheme is to use a fine mesh for the solution of photon propagation and a coarse mesh for the inversion of optical property distributions. The adaptive meshing method is accomplished by the automatic mesh refinement in the region of heterogeneity during reconstruction. These schemes are validated using tissue-like phantom measurements. Our results demonstrate the capabilities of the dual meshing and adaptive meshing in both qualitative and quantitative improvement of optical image reconstruction.

Quantitative analysis of near-infrared tomography: sensitivity to the tissue-simulating precalibration phantom

A near-infrared (NIR) imaging system is evaluated as a diagnostic clinical tool to image total hemoglobin concentration and oxygen saturation within tissue. Calibration of this type of system requires measurement of the response at each detector and source location from a homogeneous tissue-simulating phantom. The effect of using calibration phantoms of varying composition, size, and optical properties is examined to determine how it affects the overall image accuracy. All of the calibration phantoms investigated result in accurate reconstruction of absorbing heterogeneities due to increased blood concentration with less than 4% standard deviation. Images from a patient with a biopsy-confirmed ductal carcinoma are also evaluated and found to be insensitive to the choice of calibration object, with only 1% variation between images generated with different calibration objects. The tumor total hemoglobin contrast is approximately 240% higher than the average total hemoglobin concentration in contralateral breast. Soft calibration phantoms, which mimic the elastic properties of human breast tissue, are also considered and found to diminish positioning errors in the fibers relative to the actual breast exam, thereby reducing the artifacts in the periphery of the reconstructed image.

Poly(vinyl alcohol) gels for use as tissue phantoms in photoacoustic mammography

Materials for solid photoacoustic breast phantoms, based on poly(vinyl alcohol) hydrogels, are presented. Phantoms intended for use in photoacoustics must possess both optical and acoustic properties of tissue. To realize the optical properties of tissue, one approach was to optimize the number of freezing and thawing cycles of aqueous poly(vinyl alcohol) solutions, a procedure which increases the turbidity of the gel while rigidifying it. The second approach concentrated on forming a clear matrix of the rigid poly(vinyl alcohol) gel without any scattering, so that appropriate amounts of optical scatterers could be added at the time of formation, to tune the optical properties as per requirement. The relevant optical and acoustic properties of such samples were measured to be close to the average properties of human breast tissue. Tumour simulating gel samples of suitable absorption coefficient were created by adding appropriate quantities of dye at the time of formation; the samples were then cut into spheres. A breast phantom embedded with such 'tumours' was developed for studying the applicability of photoacoustics in mammography.

Three-dimensional diffuse optical tomography in the parallel plane transmission geometry: evaluation of a hybrid frequency domain/continuous wave clinical system for breast imaging

Three-dimensional diffuse optical tomography (DOT) of breast requires large data sets for even modest resolution (1 cm). We present a hybrid DOT system that combines a limited number of frequency domain (FD) measurements with a large set of continuous wave (cw) measurements. The FD measurements are used to quantitatively determine tissue averaged absorption and scattering coefficients. The larger cw data sets (10(5) measurements) collected with a lens coupled CCD, permit 3D DOT reconstructions of a 1-liter tissue volume. To address the computational complexity of large data sets and 3D volumes we employ finite difference based reconstructions computed in parallel. Tissue phantom measurements evaluate imaging performance. The tests include the following: point spread function measures of resolution, characterization of the size and contrast of single objects, field of view measurements and spectral characterization of constituent concentrations. We also report in vivo measurements. Average tissue optical properties of a healthy breast are used to deduce oxy- and deoxy-hemoglobin concentrations. Differential imaging with a tumor simulating target adhered to the surface of a healthy breast evaluates the influence of physiologic fluctuations on image noise. This tomography system provides robust, quantitative, full 3D image reconstructions with the advantages of high data throughput, single detector-tissue coupling path, and large (1L) imaging domains. In addition, we find that point spread function measurements provide a useful and comprehensive representation of system performance.

Analysis on performance and optimization of frequency-domain near-infrared instruments

Frequency-domain near-infrared techniques have been widely used to detect the optical properties of biological tissues noninvasively. In this paper we propose an analytical model to evaluate the performance of frequency-domain instruments. Based on the diffusion equation and the transfer properties of optoelectronic components, we treat all parts, including the medium, as two-port networks and apply systematic methods to answer questions concerning frequency-domain instruments. Experiments show that this method can reasonably reflect the properties of the instrument within an accuracy of 7%. This kind of method can be used to design suitable instruments for various applications. We also analyze the selection of the instrument parameters to achieve optimal performance at an efficient cost using this analytical model.

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.

Seeing the light in bone metabolism imaging

Conventional in vivo imaging of bone metabolism is dominated by gamma-ray bone scintigraphy: a technique in which gamma-ray emissions from radioactively labeled regions of metabolically active bone are mapped. More recently, however, near-infrared fluorescent probes have been developed that optically emulate these radionuclides. Although still in their infancy, techniques based on the use of such functionally targeted fluorophores might one-day offer improved resolution, sensitivity and speed in bone metabolism imaging -- without any of the health risks posed by the internalization of radioactive sources.

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.

Statistical analysis of nonlinearly reconstructed near-infrared tomographic images: Part II--Experimental interpretation

Image error analysis of a diffuse near-infrared tomography (NIR) system has been carried out on simulated data using a statistical approach described in Part I of this paper (Pogue et al., 2002). The methodology is used here with experimental data acquired on phantoms with a prototype imaging system intended for characterizing breast tissue. Results show that imaging performance is not limited by random measurement error, but rather by calibration issues. The image error over the entire field of view is generally not minimized when an accurate homogeneous estimate of the phantom properties is available; however, local image error over a target region of interest (ROI) is reduced. The image reconstruction process which includes a Levenberg-Marquardt style regularization provides good minimization of the objective function, yet its reduction is not always correlated with an overall image error decrease. Minimization of the bias in an ROI which contains localized changes in the optical properties can be achieved through five to nine iterations of the algorithm. Precalibration of the algorithm through statistical evaluation of phantom studies may provide a better measure of the image accuracy than that implied by minimization of the standard objective function.

Fast analytical approximation for arbitrary geometries in diffuse optical tomography

Diffuse optical tomography is a novel imaging technique that resolves and quantifies the optical properties of objects buried in turbid media. Typically, numerical solutions of the diffusion equation are employed to construct the tomographic problem when media of complex geometries are investigated. Numerical methods offer implementation simplicity but also significant computation burden, especially when large three-dimensional reconstructions are involved. We present an alternative method of performing tomography of diffuse media of arbitrary geometries by means of an analytical approach, the Kirchhoff approximation. We show that the method is extremely efficient in computation times and consider its potential as a real-time three-dimensional imaging tool.

Multispectral near-infrared tomography: a case study in compensating for water and lipid content in hemoglobin imaging of the breast

Images of hemoglobin concentration and oxygen saturation are presented from multispectral near-infrared tomographic measurements in the breast of a woman with an invasive cancer. Images of the absorption coefficient and reduced scattering coefficient are recovered from the measured data using a finite element reconstruction algorithm based on the frequency-domain diffusion equation. Three methods of recovering the hemoglobin concentration and oxygen saturation images are presented which compensate for water and lipid absorption in different ways: (1) an assumed bulk content of water and lipids is used, (2) four chromophores are imaged, and (3) scattering power data are applied to deduce water and lipid images. In all three cases, a large increase in the hemoglobin concentration (3:1) is observed at the location of the cancer while a maximum of 15% difference is observed in the hemoglobin images between each of these methods for water and lipid compensation.

Absorbance of opaque microstructures in optically diffuse media

In this study experimental measurements are used to determine that the observed absorbance of opaque microstructures in optically diffuse media correlates with the total surface area rather than the attenuation as calculated in a nonscattering environment. The data suggest that it may be possible to use remote measurements of optical diffuse transmission to quantify surface areas of microcapillaries that are highly absorbing or larger blood vessels that can have high intrinsic attenuation because of hematocrit alone. Results obtained in a transmission geometry are insensitive to the position of the microstructure along the line between source and detector, whereas those collected in a remission geometry are highly sensitive to the depth at which the structure is located. These types of measurement involving microscopic structures embedded in diffuse media have potential application in quantifying blood vessel surface areas that contain contrast agents or other microparticles within tissue.