Multiple-slice imaging of a tissue-equivalent phantom by use of time-resolved optical tomography

Nonlinear iterative perturbation scheme with simplified spherical harmonics (SP3 ) light propagation model for quantitative photoacoustic tomography

When using quantitative photoacoustic tomography (q-PAT) reconstruction to recover the optical absorption coefficients of tissue, the commonly used diffusion equation has several limitations in the case of the objects that have small geometries and high-absorption or low-scattering areas. Furthermore, the conventional perturbation reconstruction strategy is unsatisfactory when the target tissue containing large heterogeneous features. We herein present a modified q-PAT implementation that employs the higher-order photon migration model achieving the tradeoff between mathematical rigidity and computational efficiency. Besides, a nonlinear iterative method is proposed to obtain the perturbations of optical absorption considering the updating of the sensitivity matrix in calculating the fluence perturbations. Consequently, the distribution of tissue optical properties can be recovered in a robust way even if the targets with high absorption are included. The proposed approach has been validated by simulation, phantom and in vivo experiments, exhibiting promising performances in image fidelity and quantitative feasibility for practical applications.

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.

Fabrication and Characterization of Optical Tissue Phantoms Containing Macrostructure

The rapid development of new optical imaging techniques is dependent on the availability of low-cost, customizable, and easily reproducible standards. By replicating the imaging environment, costly animal experiments to validate a technique may be circumvented. Predicting and optimizing the performance of in vivo and ex vivo imaging techniques requires testing on samples that are optically similar to tissues of interest. Tissue-mimicking optical phantoms provide a standard for evaluation, characterization, or calibration of an optical system. Homogenous polymer optical tissue phantoms are widely used to mimic the optical properties of a specific tissue type within a narrow spectral range. Layered tissues, such as the epidermis and dermis, can be mimicked by simply stacking these homogenous slab phantoms. However, many in vivo imaging techniques are applied to more spatially complex tissue where three dimensional structures, such as blood vessels, airways, or tissue defects, can affect the performance of the imaging system. This protocol describes the fabrication of a tissue-mimicking phantom that incorporates three-dimensional structural complexity using material with optical properties of tissue. Look-up tables provide India ink and titanium dioxide recipes for optical absorption and scattering targets. Methods to characterize and tune the material optical properties are described. The phantom fabrication detailed in this article has an internal branching mock airway void; however, the technique can be broadly applied to other tissue or organ structures.

Toward whole-body quantitative photoacoustic tomography of small-animals with multi-angle light-sheet illuminations

Several attempts to achieve the quantitative photoacoustic tomography (q-PAT) have been investigated using point sources or a single-angle wide-field illumination. However, these schemes normally suffer from low signal-to-noise ratio (SNR) or poor quantification in imaging applications on large-size domains, due to the limitation of ANSI-safety incidence and incompleteness in the data acquisition. We herein present a q-PAT implementation that uses multi-angle light-sheet illuminations and calibrated recovering-and-averaging iterations. The scheme can obtain more complete information on the intrinsic absorption from the multi-angle illumination mode, and collect SNR-boosted photoacoustic signals in the selected planes from the wide-field light-sheet excitation. Therefore, the sliced absorption maps over whole body of small-animals can be recovered in a measurement-flexible, noise-robust and computation-economic way. The proposed approach is validated by phantom, ex vivo and in vivo experiments, exhibiting promising performances in image fidelity and quantitative accuracy for practical applications.

Quantitative multi-spectral oxygen saturation measurements independent of tissue optical properties

Imaging of tissue oxygenation is important in several applications associated with patient care. Optical sensing is commonly applied for assessing oxygen saturation but is often restricted to local measurements or else it requires spectral and spatial information at the expense of time. Many methods proposed so far require assumptions on the properties of measured tissue. In this study we investigated a computational method that uses only multispectral information and quantitatively computes tissue oxygen saturation independently of tissue optical properties. The method is based on linear transformations of measurements in three isosbestic points. We investigated the ideal isosbestic point combination out of six isosbestic points available for measurement in the visible and near-infrared region that enable accurate oxygen saturation computation. We demonstrate this method on controlled tissue mimicking phantoms having different optical properties and validated the measurements using a gas analyzer. A mean error of 2.9 ± 2.8% O2 Sat was achieved. Finally, we performed pilot studies in tissues in-vivo by measuring dynamic changes in fingers subjected to vascular occlusion, the vasculature of mouse ears and exposed mouse organs.

Diffuse optical tomography in the presence of a chest wall

Diffuse optical tomography (DOT) has been employed to derive spatial maps of physiologically important chromophores in the human breast, but the fidelity of these images is often compromised by boundary effects such as those due to the chest wall. We explore the image quality in fast, data-intensive analytic and algebraic linear DOT reconstructions of phantoms with subcentimeter target features and large absorptive regions mimicking the chest wall. Experiments demonstrate that the chest wall phantom can introduce severe image artifacts. We then show how these artifacts can be mitigated by exclusion of data affected by the chest wall. We also introduce and demonstrate a linear algebraic reconstruction method well suited for very large data sets in the presence of a chest wall.

Three-dimensional phantoms for curvature correction in spatial frequency domain imaging

The sensitivity to surface profile of non-contact optical imaging, such as spatial frequency domain imaging, may lead to incorrect measurements of optical properties and consequently erroneous extrapolation of physiological parameters of interest. Previous correction methods have focused on calibration-based, model-based, and computation-based approached. We propose an experimental method to correct the effect of surface profile on spectral images. Three-dimensional (3D) phantoms were built with acrylonitrile butadiene styrene (ABS) plastic using an accurate 3D imaging and an emergent 3D printing technique. In this study, our method was utilized for the correction of optical properties (absorption coefficient μ(a) and reduced scattering coefficient μ(s)') of objects obtained with a spatial frequency domain imaging system. The correction method was verified on three objects with simple to complex shapes. Incorrect optical properties due to surface with minimum 4 mm variation in height and 80 degree in slope were detected and improved, particularly for the absorption coefficients. The 3D phantom-based correction method is applicable for a wide range of purposes. The advantages and drawbacks of the 3D phantom-based correction methods are discussed in details.

A hematoma detector-a practical application of instrumental motion as signal in near infra-red imaging

In this paper we discuss results based on using instrumental motion as a signal rather than treating it as noise in Near Infra-Red (NIR) imaging. As a practical application to demonstrate this approach we show the design of a novel NIR hematoma detection device. The proposed device is based on a simplified single source configuration with a dual separation detector array and uses motion as a signal for detecting changes in blood volume in the dural regions of the head. The rapid triage of hematomas in the emergency room will lead to improved use of more sophisticated/expensive imaging facilities such as CT/MRI units. We present simulation results demonstrating the viability of such a device and initial phantom results from a proof of principle device. The results demonstrate excellent localization of inclusions as well as good quantitative comparisons.

Wavelength optimization for rapid chromophore mapping using spatial frequency domain imaging

Spatial frequency-domain imaging (SFDI) utilizes multiple-frequency structured illumination and model-based computation to generate two-dimensional maps of tissue absorption and scattering properties. SFDI absorption data are measured at multiple wavelengths and used to fit for the tissue concentration of intrinsic chromophores in each pixel. This is done with a priori knowledge of the basis spectra of common tissue chromophores, such as oxyhemoglobin (ctO(2)Hb), deoxyhemoglobin (ctHHb), water (ctH(2)O), and bulk lipid. The quality of in vivo SFDI fits for the hemoglobin parameters ctO(2)Hb and ctHHb is dependent on wavelength selection, fitting parameters, and acquisition rate. The latter is critical because SFDI acquisition time is up to six times longer than planar two-wavelength multispectral imaging due to projection of multiple-frequency spatial patterns. Thus, motion artifact during in vivo measurements compromises the quality of the reconstruction. Optimal wavelength selection is examined through matrix decomposition of basis spectra, simulation of data, and dynamic in vivo measurements of a human forearm during cuff occlusion. Fitting parameters that minimize cross-talk from additional tissue chromophores, such as water and lipid, are determined. On the basis of this work, a wavelength pair of 670 nm∕850 nm is determined to be the optimal two-wavelength combination for in vivo hemodynamic tissue measurements provided that assumptions for water and lipid fractions are made in the fitting process. In our SFDI case study, wavelength optimization reduces acquisition time over 30-fold to 1.5s compared to 50s for a full 34-wavelength acquisition. The wavelength optimization enables dynamic imaging of arterial occlusions with improved spatial resolution due to reduction of motion artifacts.

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.

Linear and nonlinear reconstruction for optical tomography of phantoms with nonscattering regions

Most research in optical imaging incorrectly assumes that light transport in nonscattering regions in the head may be modeled by use of the diffusion approximation. The effect of this assumption is examined in a series of experiments on tissue-equivalent phantoms. Images from cylindrical and head-shaped phantoms with and without clear regions [simulating the cerebrospinal fluid (CSF) filled ventricles] and a clear layer (simulating the CSF layer surrounding the brain) are reconstructed with linear and nonlinear reconstruction techniques. The results suggest that absorbing and scattering perturbations can be identified reliably with nonlinear reconstruction methods when the clear regions are also present in the reference data but that the quality of the image degrades considerably if the reference data does not contain these features. Linear reconstruction performs similarly to nonlinear reconstruction, provided the clear regions are present in the reference data, but otherwise linear reconstruction fails. This study supports the use of linear reconstruction for dynamic imaging but suggests that, in all cases, image quality is likely to improve if the clear regions are modeled correctly.

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.

Optical tomography of a realistic neonatal head phantom

We have begun clinical trials of optical tomography of the neonatal brain. To validate this research, we have built and imaged an anatomically realistic, tissue-equivalent neonatal head phantom that is hollow, allowing contrasting objects to be placed inside it. Images were reconstructed by use of two finite-element meshes, one generated from a computed tomography image of the phantom and the other spherical. The phantom was filled with a liquid of the same optical properties as the outer region, and two perturbations were placed inside. These were successfully imaged with good separation between the absorption and scatter coefficients. The phantom was then refilled with a liquid of increased absorption compared with the background to simulate the brain, and the absolute properties of the two regions were found. These were used as a priori information for the complete reconstruction. Both perturbations were visible, superimposed on the increased absorption of the central region. The head-shaped mesh performed slightly better than the spherical mesh, particularly when the absorption of the central region of the phantom was increased.

State-estimation approach to the nonstationary optical tomography problem

We propose a new numerical approach to the nonstationary optical (diffusion) tomography (OT) problem. The assumption in the method is that the absorption and/or diffusion coefficients are nonstationary in the sense that they may exhibit significant changes during the time that is needed to measure data for one traditional image frame. In the proposed method, the OT problem is formulated as a state-estimation problem. Within the state-estimation formulation, the absorption and/or diffusion coefficients are considered a stochastic process. The objective is to estimate a sequence of states for the process when the state evolution model for the process, the observation model for OT experiments, and data on the exterior boundary are given. In the proposed method, the state estimates are computed by using Kalman filtering techniques. The performance of the proposed method is evaluated on the basis of synthetic data. The simulations also illustrate that further improvements to the results in nonstationary applications can be obtained by adjustment of the measurement protocol.

Assessment of an in situ temporal calibration method for time-resolved optical tomography

A 32-channel time-resolved optical imaging device is developed at University College London to produce functional images of the neonatal brain and the female breast. Reconstruction of images using time-resolved measurements of transmitted light requires careful calibration of the temporal characteristics of the measurement system. Since they can often vary over a period of time, it is desirable to evaluate these characteristics immediately after, or prior to, the acquisition of image data. A calibration technique is investigated that is based on the measurement of light back-reflected from the surface of the object being imaged. This is facilitated by coupling each detector channel with an individual source fiber. A Monte Carlo model is employed to investigate the influence of the optical properties of the object on the back-reflected signal. The results of simulations indicate that their influence may be small enough to be ignored in some cases, or could be largely accounted for by a small adjustment to the calibrated data. The effectiveness of the method is briefly demonstrated by imaging a solid object with tissue-equivalent optical properties.

Three-dimensional optical tomography of the premature infant brain

For the first time, three-dimensional images of the newborn infant brain have been generated using measurements of transmitted light. A 32-channel time-resolved imaging system was employed, and data were acquired using custom-made helmets which couple source fibres and detector bundles to the infant head. Images have been reconstructed using measurements of mean flight time relative to those acquired on a homogeneous reference phantom, and using a head-shaped 3D finite-element-based forward model with an external boundary constrained to match the measured positions of the sources and detectors. Results are presented for a premature infant with a cerebral haemorrhage predominantly located within the left ventricle. Images representing the distribution of absorption at 780 nm and 815 nm reveal an asymmetry consistent with the haemorrhage, and corresponding maps of blood volume and fractional oxygen saturation are generally within expected physiological values.

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.

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 Bayesian optical diffusion tomography with experimental data

Reconstructions of a three-dimensional absorber embedded in a scattering medium by use of frequency domain measurements of the transmitted light in a single source-detector plane are presented. The reconstruction algorithm uses Bayesian regularization and iterative coordinate descent optimization, and it incorporates estimation of the detector noise level, the source-detector coupling coefficient, and the background diffusion coefficient in addition to the absorption image. The use of multiple modulation frequencies is also investigated. The results demonstrate the utility of this algorithm, the importance of a three-dimensional model, and that out-of-plane scattering permits recovery of three-dimensional features from measurements in a single plane.

Three-dimensional imaging of a tissuelike phantom by diffusion optical tomography

Three-dimensional imaging of three small absorbers embedded in a tissuelike cylindrical solid phantom was conducted by diffusion optical tomography. Each absorber, which was 10 mm in diameter and 10 mm high, was located on the same plane in a phantom, which was 80 mm in diameter and 140 mm high. The optical properties of the phantom were similar to those of the human breast; that is, one absorber had lower absorption and the other two absorbers had higher absorption than that of the phantom. Reemission from the phantom was measured with a multichannel photon counting system. Image reconstruction was performed by our average value method. We were able to distinguish lower and higher absorbers quantitatively. This result shows that our method can diagnose not only the existence of but also a morbid state of breast cancer.

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.

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.