Imaging changes in blood volume and oxygenation in the newborn infant brain using three-dimensional optical tomography

Using a priori structural information from magnetic resonance imaging to investigate the feasibility of prostate diffuse optical tomography and spectroscopy: a simulation study

Implementation of diffuse optical tomography (DOT) for prostate cancer is challenging because the prostate is a deep-seated organ. We investigated whether diffuse optical tomography (DOT) and spectroscopy could be applied to monitor the physiology of prostate cancer using a small probe that could be placed endorectally. We manually segmented the prostate, the intraprostatic tumor, and the rectum using data from endorectal magnetic resonance imaging. These structures were reconstructed and meshed with tetrahedral finite elements in three dimensions. A 2 x 4 cm probe that has ten sources and 52 detectors were placed to face the anterior wall of the rectum in our simulation. Optical properties of the organs were obtained from the literature in the near infrared regime. Diffusion approximation was used to simulate photon migration with finite element method. Five wavelengths were used to simulate tissue absorption with realistic water, oxy- and deoxyhaemoglobin concentrations in the prostate. We combined a global search based on genetic algorithm with gradient-driven local search methods to fit the simulated data. Our results suggest that the optical properties and the concentrations of the chromophores of the prostate and the prostate cancer can be reliably recovered from the measurements using an endorectal probe. Prostate DOT is worth further investigation for clinical application.

Automated classification schemes for optical tomographic arthritis scans

We have recently developed a sagittal laser optical tomographic (SLOT) imaging system for the diagnosis and monitoring of inflammatory processes in proximal interphalangeal (PIP) joints of patients with rheumatoid arthritis (RA). While cross sectional images of distribution of optical properties can now be generated easily, clinical interpretation of these images remains a challenge. In this paper, we apply and analyse two machine learning methods for optimal identification and severity classification of RA in a data set of 78 joints. The methods surveyed include Fisher Face with Support Vector Machines (SVMs), and Transformed Mixtures of Gausians (TMG). It appears that TMG methods outperform the approach using Fisher Face with SVMs; however, the results need to be further validated in studies involving larger patient populations.

Comparison between a time-domain and a frequency-domain system for optical tomography

The quality of phase and amplitude data from two medical optical tomography systems were compared. The two systems are a 32-channel time-domain system developed at University College London (UCL) and a 16-channel frequency-domain system developed at Helsinki University of Technology (HUT). Difference data measured from an inhomogeneous and a homogeneous phantom were compared with a finite-element method (diffusion equation) and images of scattering and absorption were reconstructed based on it. The measurements were performed at measurement times between 1 and 30 s per source. The mean rms errors in the data measured by the HUT system were 3.4% for amplitude and 0.51 deg for phase, while the corresponding values for the UCL data were 6.0% and 0.46 deg, respectively. The reproducibility of the data measured with the two systems was tested with a measurement time of 5 s per source. It was 0.4% in amplitude for the HUT system and 4% for the UCL system, and 0.08 deg in phase for both systems. The image quality of the reconstructions from the data measured with the two systems were compared with several quantitative criteria. In general a higher contrast was observed in the images calculated from the HUT data.

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.

Future of brain support: multimodality monitoring

Multimodality monitoring of cerebral physiology encompasses the application of different monitoring techniques and integration of several measured physiological and biochemical variables into the assessment of brain metabolism, structure, perfusion and oxygenation status, in addition to clinical evaluation. Novel monitoring techniques include transcranial Doppler ultrasonography, neuroimaging, intracranial pressure, cerebral perfusion and cerebral blood flow monitors, brain tissue oxygen tension monitoring, microdialysis, evoked potentials and continuous electroencephalography. Multimodality monitoring enables the immediate detection and prevention of acute neurological events, as well as appropriate intervention based on a patient’s individual disease state in the neurocritical care unit. Simultaneous real-time analysis of cerebral physiological, metabolic and cardiovascular parameters has broadened knowledge regarding complex brain pathophysiology and cerebral hemodynamics. Integration of this information allows for a more precise diagnosis and optimization of management of patients with brain injury.

Three dimensional optical imaging of blood volume and oxygenation in the neonatal brain

Optical methods provide a means of monitoring cerebral oxygenation in newborn infants at risk of brain injury. A 32-channel optical imaging system has been developed with the aim of reconstructing three-dimensional images of regional blood volume and oxygenation. Full image data sets were acquired from 14 out of 24 infants studied; successful images have been reconstructed in 8 of these infants. Regional variations in cerebral blood volume and tissue oxygen saturation are present in healthy preterm infants. In an infant with a large unilateral intraventricular haemorrhage, a corresponding region of low oxygen saturation was detected. These results suggest that optical tomography may provide an appropriate technique for investigating regional cerebral haemodynamics and oxygenation at the cotside.

Time-gated optical system for depth-resolved functional brain imaging

We present a time-domain optical system for functional imaging of the adult head. We first describe the instrument, which is based on a Ti:Sapphire pulsed laser (wavelength 750-850 nm) and an intensified CCD camera enabling parallel detection of multiple fibers. We characterize the system in terms of sensitivity and signal-to-noise ratio, instrument response function, cross-talk, stability, and reproducibility. We then describe two applications of the instrument: the characterization of baseline optical properties of homogeneous scattering media, and functional brain imaging. For the second application, we developed a two-part probe consisting in two squares of 4 x 4 sources and 3 x 3 detectors. The laser source is time-multiplexed to define 4 states of 8 sources that can be turned on during the same camera frame while minimizing cross-talk. On the detection side, we use for each detector 7 fibers of different lengths creating an optical delay, and enabling simultaneous detection in 7 windows (by steps of 500 ps) for each detector. This multiple window detection allows depth sensitivity. The imaging probe was tested on dynamic phantoms and a preliminary result on an adult performing a motor task shows discrimination between superficial and cortical responses to the stimulus on both hemispheres.

Method for the three-dimensional localization of intramyocardial excitation centers using optical imaging

This study explores the possibility of localizing the excitation centers of electrical waves inside the heart wall using voltage-sensitive dyes (fluorescent or absorptive). In the present study, we propose a method for the 3-D localization of excitation centers from pairs of 2-D images obtained in two modes of observation: reflection and transillumination. Such images can be obtained using high-speed charge-coupled device (CCD) cameras and photodiode arrays with time resolution up to 0.5 ms. To test the method, we simulate optical signals produced by point sources and propagating ellipsoidal waves in 1-cm-thick slabs of myocardial tissue. Solutions of the optical diffusion equation are constructed by employing the method of images with Robin boundary conditions. The coordinates of point sources as well as of the centers of expanding waves can be accurately determined using the proposed algorithm. The method can be extended to depth estimations of the outer boundaries of the expanding wave. The depth estimates are based on ratios of spatially integrated images. The method shows high tolerance to noise and can give accurate results even at relatively low signal-to-noise ratios. In conclusion, we propose a novel and efficient algorithm for the localization of excitation centers in 3-D cardiac tissue.

Diffuse photon propagation in multilayered geometries

Diffuse optical tomography (DOT) is an emerging functional medical imaging modality which aims to recover the optical properties of biological tissue. The forward problem of the light propagation of DOT can be modelled in the frequency domain as a diffusion equation with Robin boundary conditions. In the case of multilayered geometries with piecewise constant parameters, the forward problem is equivalent to a set of coupled Helmholtz equations. In this paper, we present solutions for the multilayered diffuse light propagation for a three-layer concentric sphere model using a series expansion method and for a general layered geometry using the boundary element method (BEM). Results are presented comparing these solutions to an independent Monte Carlo model, and for an example three layered head model.

Integrated semiconductor optical sensors for chronic, minimally-invasive imaging of brain function

Intrinsic optical signal (IOS) imaging is a widely accepted technique for imaging brain activity. We propose an integrated device consisting of interleaved arrays of gallium arsenide (GaAs) based semiconductor light sources and detectors operating at telecommunications wavelengths in the near-infrared. Such a device will allow for long-term, minimally invasive monitoring of neural activity in freely behaving subjects, and will enable the use of structured illumination patterns to improve system performance. In this work we describe the proposed system and show that near-infrared IOS imaging at wavelengths compatible with semiconductor devices can produce physiologically significant images in mice, even through skull.

Dialysis-assisted fiber optic spectroscopy for in situ biomedical sensing

A miniature fiber optic spectrometer enclosed within a semipermeable (dialysis) membrane is proposed for in vivo interstitial sensing applications. The semipermeable membrane acts as a molecular filter, allowing only small molecules to pass through to the sampling volume. This filtering, in principle, should enable continuous in vivo drug sensing, removing the necessity for complex microdialysis systems. We use a biological phantom to examine the reliable detection of a fluorescence signal from small dye molecules in the presence of large fluorophores and scatterers. We find that spectral artefacts arising from scatterers and large fluorophores are substantially suppressed, simplifying the spectral analysis. In addition, the measured sampling rate of 157 s is superior to existing in vivo tissue assaying techniques such as microdialysis, which can take tens of minutes.

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.

Dynamic physiological modeling for functional diffuse optical tomography

Diffuse optical tomography (DOT) is a noninvasive imaging technology that is sensitive to local concentration changes in oxy- and deoxyhemoglobin. When applied to functional neuroimaging, DOT measures hemodynamics in the scalp and brain that reflect competing metabolic demands and cardiovascular dynamics. The diffuse nature of near-infrared photon migration in tissue and the multitude of physiological systems that affect hemodynamics motivate the use of anatomical and physiological models to improve estimates of the functional hemodynamic response. In this paper, we present a linear state-space model for DOT analysis that models the physiological fluctuations present in the data with either static or dynamic estimation. We demonstrate the approach by using auxiliary measurements of blood pressure variability and heart rate variability as inputs to model the background physiology in DOT data. We evaluate the improvements accorded by modeling this physiology on ten human subjects with simulated functional hemodynamic responses added to the baseline physiology. Adding physiological modeling with a static estimator significantly improved estimates of the simulated functional response, and further significant improvements were achieved with a dynamic Kalman filter estimator (paired t tests, n=10, P<0.05). These results suggest that physiological modeling can improve DOT analysis. The further improvement with the Kalman filter encourages continued research into dynamic linear modeling of the physiology present in DOT. Cardiovascular dynamics also affect the blood-oxygen-dependent (BOLD) signal in functional magnetic resonance imaging (fMRI). This state-space approach to DOT analysis could be extended to BOLD fMRI analysis, multimodal studies and real-time analysis.

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.

Recent advances in diffuse optical imaging

We review the current state-of-the-art of diffuse optical imaging, which is an emerging technique for functional imaging of biological tissue. It involves generating images using measurements of visible or near-infrared light scattered across large (greater than several centimetres) thicknesses of tissue. We discuss recent advances in experimental methods and instrumentation, and examine new theoretical techniques applied to modelling and image reconstruction. We review recent work on in vivo applications including imaging the breast and brain, and examine future challenges.

Practical and adequate approach to modeling light propagation in an adult head with low-scattering regions by use of diffusion theory

A practical and adequate approach to modeling light propagation in an adult head with a low-scattering cerebrospinal fluid (CSF) region by use of diffusion theory was investigated. The diffusion approximation does not hold in a nonscattering or low-scattering regions. The hybrid radiosity-diffusion method was adopted to model the light propagation in the head with a nonscattering region. In the hybrid method the geometry of the nonscattering region is acquired as a priori information. In reality, low-level scattering occurs in the CSF region and may reduce the error caused by the diffusion approximation. The partial optical path length and the spatial sensitivity profile calculated by the finite-element method agree well with those calculated by the Monte Carlo method in the case in which the transport scattering coefficient of the CSF layer is greater than 0.3 mm(-1). Because it is feasible to assume that the transport scattering coefficient of a CSF layer is 0.3 mm(-1), it is practical to adopt diffusion theory to the modeling of light propagation in an adult head as an alternative to the hybrid method.

Detection of cortical activation with time-resolved diffuse optical methods

Simulations based on diffusion theory that use a finite-element method and rely on an magnetic resonance imaging head model suggest that time-resolved diffuse optical techniques could provide information about the depth at which variations in perfusion take place and improve the detection of cortical activation. Experimental investigations were performed with sequentially driven picosecond laser diodes and an eight-channel time-correlated single-photon-counting detection system. The experimental results obtained for activation in the motor cortex, and for the Valsalva maneuver, confirm our assumptions and are in good agreement with the simulated data.

Effects of refractive index on near-infrared tomography of the breast

Near infrared (NIR) optical tomography is an imaging technique in which internal images of optical properties are reconstructed with the boundary measurements of light propagation through the medium. Recent advances in instrumentation and theory have led to the use of this method for the detection and characterization of tumors within the female breast tissue. Most image reconstruction approaches have used the diffusion approximation and have assumed that the refractive index of the breast is constant, with a bulk value of approximately 1.4. We have applied a previously reported modified diffusion approximation, in which the refractive index for different tissues can be modeled. The model was used to generate NIR data from a realistic breast geometry containing a localized anomaly. Using this simulated data, we have reconstructed optical images, both with and without correct knowledge of the refractive-index distribution to show that the modified diffusion approximation can accurately recover the anomaly given a priori knowledge of refractive index. But using a reconstruction algorithm without the use of correct a priori information regarding the refractive-index distribution is shown as recovering the anomaly but with a degraded quality, depending on the degree of refractive index mismatch. The results suggest that provided the refractive index of breast tissue is approximately 1.3-1.4, their exclusion will have minimal effect on the reconstructed images.

Time-resolved diffuse optical tomographic imaging for the provision of both anatomical and functional information about biological tissue

We present in vivo images of near-infrared (NIR) diffuse optical tomography (DOT) of human lower legs and forearm to validate the dual functions of a time-resolved (TR) NIR DOT in clinical diagnosis, i.e., to provide anatomical and functional information simultaneously. The NIR DOT system is composed of time-correlated single-photon-counting channels, and the image reconstruction algorithm is based on the modified generalized pulsed spectral technique, which effectively incorporates the TR data with reasonable computation time. The reconstructed scattering images of both the lower legs and the forearm revealed their anatomies, in which the bones were clearly distinguished from the muscles. In the absorption images, some of the blood vessels were observable. In the functional imaging, a subject was requested to do handgripping exercise to stimulate physiological changes in the forearm tissue. The images of oxyhemoglobin, deoxyhemoglobin, and total hemoglobin concentration changes in the forearm were obtained from the differential images of the absorption at three wavelengths between the exercise and the rest states, which were reconstructed with a differential imaging scheme. These images showed increases in both blood volume and oxyhemoglobin concentration in the arteries and simultaneously showed hypoxia in the corresponding muscles. All the results have demonstrated the capability of TR NIR DOT by reconstruction of the absolute images of the scattering and the absorption with a high spatial resolution that finally provided both the anatomical and functional information inside bulky biological tissues.

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

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

Grundlagen optischer und fluoreszenzgestützter Tomographie in diffusen Medien

Three-dimensional, tomographic imaging of biological tissues by means of visible light is becoming increasingly important. Current progress in the mathematical-physical modelling of photon propagation in scattering media allows spatially-resolved reconstructions of optical parameters with common computing hardware. Especially in the field of molecular imaging, optical tomography promises a transfer of knowledge from successful in vitro assays (as evaluated by fluorescence microscopy) to in vivo imaging of living animals. In the latter case, spatial resolution is not as critical as the ability to quantify the concentration of fluorescence-labelled probes, a task not solvable by the use of common planar imaging techniques. In this article, the theoretical foundations of optical tomography are introduced along with some examples of applications.

Near infrared spectroscopy in brain injury: today’s perspective

The technique of near infrared spectroscopy (NIRS) is based on the principle of light attenuation by the chromophores oxyhaemoglobin (HbO2), deoxyhaemoglobin (Hb) and cytochrome oxidase. Changes in the detected light levels can therefore represent changes in concentrations of these chromophores. Clinical use of NIRS in the brain has been well established in neonates where transillumination is possible. While it has become a useful research tool for monitoring the adult brain, clinical application has been hampered by the fact that it must be applied in reflectance mode. This has resulted in a number of concerns, most significantly the issue of signal contamination by the extracranial tissue layers. Algorithms have been applied to try to overcome this problem, and techniques such as time resolved, phase resolved and spatially resolved spectroscopy have been developed. There has been renewed interest in NIRS as an easy to use, non-invasive technique for measuring tissue oxygenation in the adult brain. Recent technical advances have led to the development of compact, portable instruments that detect changes in optical attenuation of several wavelengths of light. Near infrared spectroscopy is an evolving technology that holds significant potential for technical advancement. In particular, NIRS shows future promise as a clinical tool for bedside cerebral blood flow measurements and as a cerebral imaging modality for mapping structure and function.

A modelling error approach for the estimation of optical absorption in the presence of anisotropies

Optical tomography is an emerging method for non-invasive imaging of human tissues using near-infrared light. Generally, the tissue is assumed isotropic, but this may not always be true. In this paper, we present a method for the estimation of optical absorption coefficient allowing the background to be anisotropic. To solve the forward problem, we model the light propagation in tissue using an anisotropic diffusion equation. The inverse problem consists of the estimation of the absorption coefficient based on boundary measurements. Generally, the background anisotropy cannot be assumed to be known. We treat the uncertainties in the background anisotropy parameter values as modelling error, and include this in our model and reconstruction. We present numerical examples based on simulated data. For reference, examples using an isotropic inversion scheme are also included. The estimates are qualitatively different for the two methods.