Optical tomography of a realistic neonatal head phantom

AAPM Task Group Report 311: Guidance for performance evaluation of fluorescence-guided surgery systems

The last decade has seen a large growth in fluorescence-guided surgery (FGS) imaging and interventions. With the increasing number of clinical specialties implementing FGS, the range of systems with radically different physical designs, image processing approaches, and performance requirements is expanding. This variety of systems makes it nearly impossible to specify uniform performance goals, yet at the same time, utilization of different devices in new clinical procedures and trials indicates some need for common knowledge bases and a quality assessment paradigm to ensure that effective translation and use occurs. It is feasible to identify key fundamental image quality characteristics and corresponding objective test methods that should be determined such that there are consistent conventions across a variety of FGS devices. This report outlines test methods, tissue simulating phantoms and suggested guidelines, as well as personnel needs and professional knowledge bases that can be established. This report frames the issues with guidance and feedback from related societies and agencies having vested interest in the outcome, coming from an independent scientific group formed from academics and international federal agencies for the establishment of these professional guidelines.

Optical Epidermal Mimicry from Ultraviolet to Infrared Wavelengths

Optical tissue phantoms present substantial value for medical imaging and therapeutic applications. We have developed an epidermal tissue phantom to mimic the optical properties of human skin from the ultraviolet to the infrared region, exceeding the breadth of existing studies. An epoxy matrix is combined with melanin-mimicking polydopamine via a cost-effective fabrication strategy. Reflectance and transmittance measurements enable calculation of the wavelength-dependent complex refractive index and absorption coefficient. Results are compared with literature data to establish agreement with a real human epidermis. By analyzing emissive power at a typical skin temperature, the epidermal tissue phantom is shown to accurately mimic the radiative properties of human skin. This simple, multifunctional material represents a promising substitute for human tissue for a variety of medical and bioengineering applications.

Mechanically switchable solid inhomogeneous phantom for performance tests in diffuse imaging and spectroscopy

A mechanically switchable solid inhomogeneous phantom simulating localized absorption changes was developed and characterized. The homogeneous host phantom was made of epoxy resin with black toner and titanium dioxide particles added as absorbing and scattering components, respectively. A cylindrical rod, movable along a hole in the block and made of the same material, has a black polyvinyl chloride cylinder embedded in its center. By varying the volume and position of the black inclusion, absorption perturbations can be generated over a large range of magnitudes. The phantom has been characterized by various time-domain diffuse optics instruments in terms of absorption and scattering spectra, transmittance images, and reflectance contrast. Addressing a major application of the phantom for performance characterization for functional near-infrared spectroscopy of the brain, the contrast was measured in reflectance mode while black cylinders of volumes from ≈20  mm3 to ≈270  mm3 were moved in lateral and depth directions, respectively. The new type of solid inhomogeneous phantom is expected to become a useful tool for routine quality check of clinical instruments or implementation of industrial standards provided an experimental characterization of the phantom is performed in advance.

Nonlinear approach to difference imaging in diffuse optical tomography

Difference imaging aims at recovery of the change in the optical properties of a body based on measurements before and after the change. Conventionally, the image reconstruction is based on using difference of the measurements and a linear approximation of the observation model. One of the main benefits of the linearized difference reconstruction is that the approach has a good tolerance to modeling errors, which cancel out partially in the subtraction of the measurements. However, a drawback of the approach is that the difference images are usually only qualitative in nature and their spatial resolution can be weak because they rely on the global linearization of the nonlinear observation model. To overcome the limitations of the linear approach, we investigate a nonlinear approach for difference imaging where the images of the optical parameters before and after the change are reconstructed simultaneously based on the two datasets. We tested the feasibility of the method with simulations and experimental data from a phantom and studied how the approach tolerates modeling errors like domain truncation, optode coupling errors, and domain shape errors.

3D Rapid Prototyping for Otolaryngology-Head and Neck Surgery: Applications in Image-Guidance, Surgical Simulation and Patient-Specific Modeling

The aim of this study was to demonstrate the role of advanced fabrication technology across a broad spectrum of head and neck surgical procedures, including applications in endoscopic sinus surgery, skull base surgery, and maxillofacial reconstruction. The initial case studies demonstrated three applications of rapid prototyping technology are in head and neck surgery: i) a mono-material paranasal sinus phantom for endoscopy training ii) a multi-material skull base simulator and iii) 3D patient-specific mandible templates. Digital processing of these phantoms is based on real patient or cadaveric 3D images such as CT or MRI data. Three endoscopic sinus surgeons examined the realism of the endoscopist training phantom. One experienced endoscopic skull base surgeon conducted advanced sinus procedures on the high-fidelity multi-material skull base simulator. Ten patients participated in a prospective clinical study examining patient-specific modeling for mandibular reconstructive surgery. Qualitative feedback to assess the realism of the endoscopy training phantom and high-fidelity multi-material phantom was acquired. Conformance comparisons using assessments from the blinded reconstructive surgeons measured the geometric performance between intra-operative and pre-operative reconstruction mandible plates. Both the endoscopy training phantom and the high-fidelity multi-material phantom received positive feedback on the realistic structure of the phantom models. Results suggested further improvement on the soft tissue structure of the phantom models is necessary. In the patient-specific mandible template study, the pre-operative plates were judged by two blinded surgeons as providing optimal conformance in 7 out of 10 cases. No statistical differences were found in plate fabrication time and conformance, with pre-operative plating providing the advantage of reducing time spent in the operation room. The applicability of common model design and fabrication techniques across a variety of otolaryngological sub-specialties suggests an emerging role for rapid prototyping technology in surgical education, procedure simulation, and clinical practice.

Compensation of modeling errors due to unknown domain boundary in diffuse optical tomography

Diffuse optical tomography is a highly unstable problem with respect to modeling and measurement errors. During clinical measurements, the body shape is not always known, and an approximate model domain has to be employed. The use of an incorrect model domain can, however, lead to significant artifacts in the reconstructed images. Recently, the Bayesian approximation error theory has been proposed to handle model-based errors. In this work, the feasibility of the Bayesian approximation error approach to compensate for modeling errors due to unknown body shape is investigated. The approach is tested with simulations. The results show that the Bayesian approximation error method can be used to reduce artifacts in reconstructed images due to unknown domain shape.

Phantoms for diffuse optical imaging based on totally absorbing objects, part 2: experimental implementation

We present the experimental implementation and validation of a phantom for diffuse optical imaging based on totally absorbing objects for which, in the previous paper [J. Biomed. Opt.18(6), 066014, (2013)], we have provided the basic theory. Totally absorbing objects have been manufactured as black polyvinyl chloride (PVC) cylinders and the phantom is a water dilution of intralipid-20% as the diffusive medium and India ink as the absorber, filled into a black scattering cell made of PVC. By means of time-domain measurements and of Monte Carlo simulations, we have shown the reliability, the accuracy, and the robustness of such a phantom in mimicking typical absorbing perturbations of diffuse optical imaging. In particular, we show that such a phantom can be used to generate any absorption perturbation by changing the volume and position of the totally absorbing inclusion.

Detectability of absorption and reduced scattering coefficients in frequency-domain measurements using a realistic head phantom

Detection limits of the changes in absorption and reduced scattering coefficients were investigated using a frequency-domain near-infrared system in a realistic head phantom. The results were quantified in terms of the maximum detectable depth for different activation volumes in the range of 0.8-20 microliters. The non-linear relation between the maximum detectable depth and the magnitude of changes in the absorption coefficient conform well with the Born approximation to the diffusion equation. The minimal detectable changes in the reduced scattering coefficient measured in terms of the phase signal were found to be approximately twice as large as that of the absorption coefficient using the AC signal for the same volume and at the same depth. The phase delay, which can be used to quantify the fast neuronal optical response in the human brain, showed a linear dependence on the reciprocal of the reduced scattering coefficient, as predicted by the Rytov approximation.

Time-domain diffuse optical tomography processing by using the Mellin-Laplace transform

We investigate the use of the Mellin-Laplace transform for reconstructing optical parameters from time-resolved optical tomography in diffusive media. We present here its definition, its mathematical properties, and its sensitivity to variations of optical properties. The method was validated on two-dimensional reconstructions from simulation in the reflection geometry. We conclude that reconstructions based on the Mellin-Laplace transform are more robust to noise than the methods using first moments.

High-density diffuse optical tomography of term infant visual cortex in the nursery

Advancements in antenatal and neonatal medicine over the last few decades have led to significant improvement in the survival rates of sick newborn infants. However, this improvement in survival has not been matched by a reduction in neurodevelopmental morbidities with increasing recognition of the diverse cognitive and behavioral challenges that preterm infants face in childhood. Conventional neuroimaging modalities, such as cranial ultrasound and magnetic resonance imaging, provide an important definition of neuroanatomy with recognition of brain injury. However, they fail to define the functional integrity of the immature brain, particularly during this critical developmental period. Diffuse optical tomography methods have established success in imaging adult brain function; however, few studies exist to demonstrate their feasibility in the neonatal population. We demonstrate the feasibility of using recently developed high-density diffuse optical tomography (HD-DOT) to map functional activation of the visual cortex in healthy term-born infants. The functional images show high contrast-to-noise ratio obtained in seven neonates. These results illustrate the potential for HD-DOT and provide a foundation for investigations of brain function in more vulnerable newborns, such as preterm infants.

A quantitative spatial comparison of high-density diffuse optical tomography and fMRI cortical mapping

Functional neuroimaging commands a dominant role in current neuroscience research. However its use in bedside clinical and certain neuro-scientific studies has been limited because the current tools lack the combination of being non-invasive, non-ionizing and portable while maintaining moderate resolution and localization accuracy. Optical neuroimaging satisfies many of these requirements, but, until recent advances in high-density diffuse optical tomography (HD-DOT), has been hampered by limited resolution. While early results of HD-DOT have been promising, a quantitative voxel-wise comparison and validation of HD-DOT against the gold standard of functional magnetic resonance imaging (fMRI) has been lacking. Herein, we provide such an analysis within the visual cortex using matched visual stimulation protocols in a single group of subjects (n=5) during separate HD-DOT and fMRI scanning sessions. To attain the needed voxel-to-voxel co-registration between HD-DOT and fMRI image spaces, we implemented subject-specific head modeling that incorporated MRI anatomy, detailed segmentation, and alignment of source and detector positions. Comparisons of the visual responses found an average localization error between HD-DOT and fMRI of 4.4+/-1mm, significantly less than the average distance between cortical gyri. This specificity demonstrates that HD-DOT has sufficient image quality to be useful as a surrogate for fMRI.

Biomedical tissue phantoms with controlled geometric and optical properties for Raman spectroscopy and tomography

To support the translation of Raman spectroscopy into clinical applications, synthetic models are needed to accurately test, optimize and validate prototype fiber optic instrumentation. Synthetic models (also called tissue phantoms) are widely used for developing and testing optical instrumentation for diffuse reflectance, fluorescence, and Raman spectroscopies. While existing tissue phantoms accurately model tissue optical scattering and absorption, they do not typically model the anatomic shapes and chemical composition of tissue. Because Raman spectroscopy is sensitive to molecular composition, Raman tissue phantoms should also approximate the bulk tissue composition. We describe the fabrication and characterization of tissue phantoms for Raman tomography and spectroscopy. These phantoms have controlled chemical and optical properties, and also multilayer morphologies which approximate the appropriate anatomic shapes. Tissue phantoms were fabricated to support on-going Raman studies by simulating the human wrist and rat leg. Surface meshes (triangle patch models) were generated from computed tomography (CT) images of a human arm and rat leg. Rapid prototyping was used to print mold templates with complex geometric patterns. Plastic casting techniques used for movie special effects were adapted to fabricate molds from the rapid prototypes, and finally to cast multilayer gelatin tissue phantoms. The gelatin base was enriched with additives to model the approximate chemistry and optical properties of individual tissue layers. Additional studies were performed to determine optimal casting conditions, phantom stability, layer delamination and chemical diffusion between layers. Recovery of diffuse reflectance and Raman spectra in tissue phantoms varied with probe placement. These phantoms enable optimization of probe placement for human or rat studies. These multilayer tissue phantoms with complex geometries are shown to be stable, with minimal layer delamination and chemical diffusion.

Bedside optical imaging of occipital resting-state functional connectivity in neonates

Resting-state networks derived from temporal correlations of spontaneous hemodynamic fluctuations have been extensively used to elucidate the functional organization of the brain in adults and infants. We have previously developed functional connectivity diffuse optical tomography methods in adults, and we now apply these techniques to study functional connectivity in newborn infants at the bedside. We present functional connectivity maps in the occipital cortices obtained from healthy term-born infants and premature infants, including one infant with an occipital stroke. Our results suggest that functional connectivity diffuse optical tomography has potential as a valuable clinical tool for the early detection of functional deficits and for providing prognostic information on future development.

Implementation of a computationally efficient least-squares algorithm for highly under-determined three-dimensional diffuse optical tomography problems

Three-dimensional (3D) diffuse optical tomography is known to be a nonlinear, ill-posed and sometimes under-determined problem, where regularization is added to the minimization to allow convergence to a unique solution. In this work, a generalized least-squares (GLS) minimization method was implemented, which employs weight matrices for both data-model misfit and optical properties to include their variances and covariances, using a computationally efficient scheme. This allows inversion of a matrix that is of a dimension dictated by the number of measurements, instead of by the number of imaging parameters. This increases the computation speed up to four times per iteration in most of the under-determined 3D imaging problems. An analytic derivation, using the Sherman-Morrison-Woodbury identity, is shown for this efficient alternative form and it is proven to be equivalent, not only analytically, but also numerically. Equivalent alternative forms for other minimization methods, like Levenberg-Marquardt (LM) and Tikhonov, are also derived. Three-dimensional reconstruction results indicate that the poor recovery of quantitatively accurate values in 3D optical images can also be a characteristic of the reconstruction algorithm, along with the target size. Interestingly, usage of GLS reconstruction methods reduces error in the periphery of the image, as expected, and improves by 20% the ability to quantify local interior regions in terms of the recovered optical contrast, as compared to LM methods. Characterization of detector photo-multiplier tubes noise has enabled the use of the GLS method for reconstructing experimental data and showed a promise for better quantification of target in 3D optical imaging. Use of these new alternative forms becomes effective when the ratio of the number of imaging property parameters exceeds the number of measurements by a factor greater than 2.

An electrically-activated dynamic tissue-equivalent phantom for assessment of diffuse optical imaging systems

A novel design of solid dynamic phantom with tissue-like optical properties is presented, which contains variable regions of contrast which are activated electrically. Reversible changes in absorption are produced by localized heating of targets impregnated with thermochromic pigment. A portable, battery-operated prototype has been constructed, and its optical and temporal characteristics have been investigated. The phantom has been developed as a means of assessing the performance of diffuse optical imaging systems, such as those used to monitor haemodynamic changes in the brain and other tissues. Images of the phantom have been reconstructed using data acquired with a continuous wave optical topography system.

Optical brain imaging in vivo: techniques and applications from animal to man

Optical brain imaging has seen 30 years of intense development, and has grown into a rich and diverse field. In-vivo imaging using light provides unprecedented sensitivity to functional changes through intrinsic contrast, and is rapidly exploiting the growing availability of exogenous optical contrast agents. Light can be used to image microscopic structure and function in vivo in exposed animal brain, while also allowing noninvasive imaging of hemodynamics and metabolism in a clinical setting. This work presents an overview of the wide range of approaches currently being applied to in-vivo optical brain imaging, from animal to man. Techniques include multispectral optical imaging, voltage sensitive dye imaging and speckle-flow imaging of exposed cortex, in-vivo two-photon microscopy of the living brain, and the broad range of noninvasive topography and tomography approaches to near-infrared imaging of the human brain. The basic principles of each technique are described, followed by examples of current applications to cutting-edge neuroscience research. In summary, it is shown that optical brain imaging continues to grow and evolve, embracing new technologies and advancing to address ever more complex and important neuroscience questions.

Image reconstruction in optical tomography in the presence of coupling errors

Image reconstruction in optical tomography is a nonlinear and generally ill- posed inverse problem. Noise in the measured surface data can give rise to substantial artifacts in the recovered volume images of optical coefficients. Apart from random shot noise caused by the limited number of photons detected at the measurement site, another class of systematic noise is associated with losses specific to individual source and detector locations. A common cause for such losses in data acquisition systems based on fiber-optic light delivery is the imperfect coupling between the fiber tips and the skin of the patient because of air gaps or surface moisture. Thus the term coupling errors was coined for this type of data noise. However, source and detector specific errors can also occur in noncontact measurement systems not using fiber-optic delivery, for example, owing to local skin pigmentation, hair and hair follicles, or instrumentation calibration errors. Often it is not possible to quantify coupling effects in a way that allows us to remove them from the data or incorporate them into the light transport model. We present an alternative method of eliminating coupling errors by regarding the complex-valued coupling factors for each source and detector as unknowns in the reconstruction process and recovering them simultaneously with the images of absorption and scattering. Our method takes into account the possibility that coupling effects have an influence on both the amplitude and the phase shift of the measurements. Reconstructions from simulated and experimental phantom data are presented, which show that including the coupling coefficients in the reconstruction greatly improves the recovery of absorption and scattering images.

Optical tomography of the neonatal brain

A new method of assessing neurological function and pathology in the newborn infant is being developed based on the transmission of near-infrared light across the brain. Absorption by blood over a range of wavelengths reveals a strong dependency on oxygenation status, and measurements of transmitted light enable the spatial variation in the concentrations of the oxygenated and de-oxygenated forms of hemoglobin to be derived. Optical tomography has so far provided static three-dimensional maps of blood volume and oxygenation as well as dynamic images revealing the brain's response to sensory stimulation and global hemodynamic changes. The imaging modality is being developed as a safe and non-invasive tool that can be utilized at the cotside in intensive care. Optical tomography of the healthy infant brain is also providing a means of studying neurophysiological processes during early development and the potential consequences of prematurity.

Multiplicative optical tomography of cardiac electrical activity

Cardiac electrical activity can be mapped today through the response of voltage-sensitive dyes; but poor transparency of muscle tissue has enforced shallow-depth imaging. We present a three-dimensional (3D) reconstruction method for electrical activity deep inside the myocardial wall. Our approach is nonlinear and differs substantially from standard diffusive optical tomography. It does not require matrix inversion, data regularization or a priori information concerning the original object. Opposite sides of a slab-shaped preparation are scanned in parallel by detection and illumination points with a constant vector offset between illumination and detection axes (biaxial scanning). Scanning is performed in two perpendicular directions. In each direction, a pair of 2D images is obtained under offsets of opposite signs. These two pairs are the input for a multiplicative reconstruction algorithm, whose output is a 3D image. The overall procedure was successfully tested on computer-generated sources that include points, lines and hemispheres, patterned after actual electrophysiological excitations. The algorithm is computationally efficient and stable with respect to varying noise levels in the raw data.

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.

Effect of errors in baseline optical properties on accuracy of transabdominal near-infrared spectroscopy in fetal sheep brain during hypoxic stress

A continuous-wave (cw) near-infrared spectroscopy (NIRS) instrument has been developed to noninvasively quantify fetal cerebral blood oxygen saturation (StO2). A linear Green's function formulism was used to analytically solve the photon diffusion equation and extract the time-varying fetal tissue oxy- and deoxy-hemoglobin concentrations from the NIR measurements. Here we explored the accuracy with which this instrument can be expected to perform over a range of fetal hypoxic states. We investigated the dependence of this accuracy on the accuracy of the reference optical properties chosen based on the literature. The fetal oxygenation of a pregnant ewe model was altered via maternal aortic occlusion. The NIR cw instrument was placed on the maternal abdomen directly above the fetal head, continuously acquiring diffuse optical measurements. Blood was sampled periodically from the fetus to obtain fetal arterial saturation (SaO2) measurements from blood gas analysis. The NIR StO2 values were compared with the fetal SaO2 measurements. Variations in the NIR results due to uncertainty in the reference optical properties were relatively small within the fetal SaO2 range of 30 to 80%. Under hypoxic conditions, however, the variability of the NIR StO2 calculations with changes in the assumed reference properties became more significant.

Three-dimensional whole-head optical tomography of passive motor evoked responses in the neonate

Optical tomography has been used to reconstruct three-dimensional images of the entire neonatal head during motor evoked responses. Data were successfully acquired during passive movement of each arm on four out of six infants examined, from which eight sets of bilateral images of hemodynamic parameters were reconstructed. Six out of the eight images showed the largest change in total hemoglobin in the region of the contralateral motor cortex. The mean distance between the peak response in the image and the estimated position of the contralateral motor cortex was 10.8 mm. These results suggest that optical tomography may provide an appropriate technique for non-invasive cot-side imaging of brain function.

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.

Magnetic-resonance-imaging-coupled broadband near-infrared tomography system for small animal brain studies

A novel magnetic-resonance-coupled broadband near-infrared (NIR) tomography system for small animal brain studies is described. Several features of the image formation approach are new in NIR tomography and represent major advances in the path to recovering high-resolution hemoglobin and oxygen saturation images of tissue. The NIR data were broadband and continuous wave and were used along with a second-derivative-based estimation of the path length from water absorption. The path length estimation from water was then used along with the attenuation spectrum to recover absorption and reduced scattering coefficient images at multiple wavelengths and then to recover images of total hemoglobin and oxygen saturation. Going beyond these basics of NIR tomography, software has been developed to allow inclusion of structures derived from MR imaging (MRI) for the external and internal tissue boundaries, thereby improving the accuracy and spatial resolution of the properties in each tissue type. The system has been validated in both tissue-simulating phantoms, with 10% accuracy observed, and in a rat cranium imaging experiment. The latter experiment used variation in inspired oxygen (FiO2) to vary the observed hemoglobin and oxygen saturation images. Quantitative agreement was observed between the changes in deoxyhemoglobin values derived from NIR and the changes predicted with blood-oxygen-level-dependent (BOLD) MRI. This system represents the initial stage in what will likely be a larger role for NIR tomography, coupled to MRI, and illustrates that the technological challenges of using continuous-wave broadband data and inclusion of a priori structural information can be met with careful phantom studies.

Excitation patterns in three-dimensional electrical impedance tomography

Electrical impedance tomography (EIT) is a non-invasive technique that aims to reconstruct images of internal electrical properties of a domain, based on electrical measurements on the periphery. Improvements in instrumentation and numerical modeling have led to three-dimensional (3D) imaging. The availability of 3D modeling and imaging raises the question of identifying the best possible excitation patterns that will yield to data, which can be used to produce the best image reconstruction of internal properties. In this work, we describe our 3D finite element model of EIT. Through singular value decomposition as well as examples of reconstructed images, we show that for a homogenous female breast model with four layers of electrodes, a driving pattern where each excitation plane is a sinusoidal pattern out-of-phase with its neighboring plane produces better qualitative images. However, in terms of quantitative imaging an excitation pattern where all electrode layers are in phase produces better results.

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.

Theoretical evaluation of accuracy in position and size of brain activity obtained by near-infrared topography

Near-infrared (NIR) topography can obtain a topographical distribution of the activated region in the brain cortex. Near-infrared light is strongly scattered in the head, and the volume of tissue sampled by a source-detector pair on the head surface is broadly distributed in the brain. This scattering effect results in poor resolution and contrast in the topographic image of the brain activity. In this study, a one-dimensional distribution of absorption change in a head model is calculated by mapping and reconstruction methods to evaluate the effect of the image reconstruction algorithm and the interval of measurement points for topographic imaging on the accuracy of the topographic image. The light propagation in the head model is predicted by Monte Carlo simulation to obtain the spatial sensitivity profile for a source-detector pair. The measurement points are one-dimensionally arranged on the surface of the model, and the distance between adjacent measurement points is varied from 4 mm to 28 mm. Small intervals of the measurement points improve the topographic image calculated by both the mapping and reconstruction methods. In the conventional mapping method, the limit of the spatial resolution depends upon the interval of the measurement points and spatial sensitivity profile for source-detector pairs. The reconstruction method has advantages over the mapping method which improve the results of one-dimensional analysis when the interval of measurement points is less than 12 mm. The effect of overlapping of spatial sensitivity profiles indicates that the reconstruction method may be effective to improve the spatial resolution of a two-dimensional reconstruction of topographic image obtained with larger interval of measurement points. Near-infrared topography with the reconstruction method potentially obtains an accurate distribution of absorption change in the brain even if the size of absorption change is less than 10 mm.

Three-dimensional optical tomography: resolution in small-object imaging

Near-infrared (NIR) optical tomography provide estimates of the internal distribution of optical absorption and transport scattering from boundary measurement of light propagation within biological tissue. Although this is a truly three-dimensional (3D) imaging problem, most research to date has concentrated on two-dimensional modeling and image reconstruction. More recently, 3D imaging algorithms are demonstrating better estimation of the light propagation within the imaging region and are providing the basis of more accurate image construction algorithms. As 3D methods emerge, it will become increasingly important to evaluate their resolution, contrast, and localization of optical property heterogeneity. We present a concise study of 3D reconstructed resolution of a small, low-contrast, absorbing and scattering anomaly as it is placed in different locations within a cylindrical phantom. The object is an 8-mm-diameter cylinder, which represents a typical small target that needs to be resolved in NIR mammographic imaging. The best resolution and contrast is observed when the object is located near the periphery of the imaging region (12-22 mm from the edge) and is also positioned within the multiple measurement planes, with the most accurate results seen for the scatter image when the anomaly is at 17 mm from the edge. Furthermore, the accuracy of quantitative imaging is increased to almost 100% of the target values when a priori information regarding the internal structure of imaging domain is utilized.