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

Data-driven approach to optimum wavelength selection for diffuse optical imaging

The production of accurate and independent images of the changes in concentration of oxyhemoglobin and deoxyhemoglobin by diffuse optical imaging is heavily dependent on which wavelengths of near-infrared light are chosen to interrogate the target tissue. Although wavelengths can be selected by theoretical methods, in practice the accuracy of reconstructed images will be affected by wavelength-specific and system-specific factors such as laser source power and detector sensitivity. We describe the application of a data-driven approach to optimum wavelength selection for the second generation of University College London's multichannel, time-domain optical tomography system (MONSTIR II). By performing a functional activation experiment using 12 different wavelengths between 690 and 870 nm, we were able to identify the combinations of 2, 3, and 4 wavelengths which most accurately reproduced the results obtained using all 12 wavelengths via an imaging approach. Our results show that the set of 2, 3, and 4 wavelengths which produce the most accurate images of functional activation are [770, 810], [770, 790, 850], and [730, 770, 810, 850] respectively, but also that the system is relatively robust to wavelength selection within certain limits. Although these results are specific to MONSTIR II, the approach we developed can be applied to other multispectral near-infrared spectroscopy and optical imaging systems.

MONSTIR II: a 32-channel, multispectral, time-resolved optical tomography system for neonatal brain imaging

We detail the design, construction and performance of the second generation UCL time-resolved optical tomography system, known as MONSTIR II. Intended primarily for the study of the newborn brain, the system employs 32 source fibres that sequentially transmit picosecond pulses of light at any four wavelengths between 650 and 900 nm. The 32 detector channels each contain an independent photo-multiplier tube and temporally correlated photon-counting electronics that allow the photon transit time between each source and each detector position to be measured with high temporal resolution. The system's response time, temporal stability, cross-talk, and spectral characteristics are reported. The efficacy of MONSTIR II is demonstrated by performing multi-spectral imaging of a simple phantom.

Multimodality characterization of microstructure by the combination of diffusion NMR and time-domain diffuse optical data

Combining datasets with a model of the underlying physics prior to mapping of tissue provides a novel approach improving the estimation of parameters. We demonstrate this approach by merging near infrared diffuse optical signal data with diffusion NMR data to inform a model describing the microstructure of a sample. The study is conducted on a homogeneous emulsion of oil in a dispersion medium of water and proteins. The use of a protein based background, rich in collagen, introduces a similarity to real tissues compared to other models such as intralipids. The sample is investigated with the two modalities separately. Then, the two datasets are used to inform a combined model, and to estimate the size of the microstructural elements and the volume fraction. The combined model fits the microstructural properties by minimizing the difference between experimental and modelled data. The experimental results are validated with confocal laser scanning microscopy. The final results demonstrate that the combined model provides improved estimates of microstructural parameters compared to either individual model alone.

The reproducibility of optical mammography in healthy volunteers

This study was designed to determine the reproducibility of optical mammography. Eight healthy pre-menopausal volunteers were scanned at different time intervals (minutes, weeks and months apart) to investigate the effects of within-subject variation, between-subject variation and systematic variations on both the raw data and images. The study shows that the greatest source of variation in optical mammography raw data and images is between different subjects, and scans of the same subject are very reproducible. The averaged total haemoglobin concentration from the eight volunteers was (24 ± 10) µM, and the average tissue oxygen saturation was (70 ± 10)%, which is comparable with other data in the literature. The average absorption coefficient at 780 nm was (0.0048 ± 0.0017) mm(-1) and the average reduced scatter coefficient at 780 nm was (0.80 ± 0.12) mm(-1). Again, this is comparable with published values. When our data are combined with the published values, the weighted average total haemoglobin concentration and tissue oxygen saturation for pre-menopausal breasts are (29 ± 8) µM and (73 ± 3)%, respectively. The results of our study show that we can be reassured that any changes within the tumour region seen during neoadjuvant therapy are likely to be due to a real physiological response to treatment, as the physiological properties of the breast remain relatively constant. However, in this study, we cannot distinguish between a tumour response to treatment and systemic changes in the healthy breast.

Time-domain reflectance diffuse optical tomography with Mellin-Laplace transform for experimental detection and depth localization of a single absorbing inclusion

We show how to apply the Mellin-Laplace transform to process time-resolved reflectance measurements for diffuse optical tomography. We illustrate this method on simulated signals incorporating the main sources of experimental noise and suggest how to fine-tune the method in order to detect the deepest absorbing inclusions and optimize their localization in depth, depending on the dynamic range of the measurement. To finish, we apply this method to measurements acquired with a setup including a femtosecond laser, photomultipliers and a time-correlated single photon counting board. Simulations and experiments are illustrated for a probe featuring the interfiber distance of 1.5 cm and show the potential of time-resolved techniques for imaging absorption contrast in depth with this geometry.

Effect of task-related extracerebral circulation on diffuse optical tomography: experimental data and simulations on the forehead

The effect of task-related extracerebral circulatory changes on diffuse optical tomography (DOT) of brain activation was evaluated using experimental data from 14 healthy human subjects and computer simulations. Total hemoglobin responses to weekday-recitation, verbal-fluency, and hand-motor tasks were measured with a high-density optode grid placed on the forehead. The tasks caused varying levels of mental and physical stress, eliciting extracerebral circulatory changes that the reconstruction algorithm was unable to fully distinguish from cerebral hemodynamic changes, resulting in artifacts in the brain activation images. Crosstalk between intra- and extracranial layers was confirmed by the simulations. The extracerebral effects were attenuated by superficial signal regression and depended to some extent on the heart rate, thus allowing identification of hemodynamic changes related to brain activation during the verbal-fluency task. During the hand-motor task, the extracerebral component was stronger, making the separation less clear. DOT provides a tool for distinguishing extracerebral components from signals of cerebral origin. Especially in the case of strong task-related extracerebral circulatory changes, however, sophisticated reconstruction methods are needed to eliminate crosstalk artifacts.

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.

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.

Time-resolved optical mammography using a liquid coupled interface

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

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.

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.

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.

Determination of the concentration scaling law of the scattering coefficient of water solutions of Intralipid at 832 nm by comparison between collimated detection measurements and Monte Carlo simulations

BACKGROUND AND OBJECTIVES

Intralipid (IP) is a scatterer extensively used in the building of phantoms for Biomedical Optics measurements. Recently, deviations from the linearity have been shown for the concentration scaling law of the scattering coefficient of IP water solutions at visible wavelengths. In this work this scaling law was determined at 832 nm.

STUDY DESIGN/MATERIALS AND METHODS

Space resolved transmittance measurements of a laser beam at 832 nm through water solutions of IP and ink were performed and compared with the corresponding results of Monte Carlo simulations.

RESULTS

The comparison provides a quadratic dependence of mu'(s) on the volume-to-volume scatterer concentration, C(IP), in the range of C(IP) values (0.0024<C(IP)<0.0075). These deviations from the linear behavior are related to the failure of the independent scatterer approximation.

CONCLUSION

The quadratic dependence of mu'(s) on C(IP) is in agreement with recent results obtained by other groups with different experimental techniques and is validated by a recent theoretical work.

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

Induced haemodynamic and blood oxygenation changes occurring within the brain of a ventilated newborn infant have been imaged in three dimensions using optical tomography. Noninvasive measurements of the flight times of transmitted light were acquired during illumination of the brain by laser pulses at wavelengths of 780 nm and 815 nm. The oxygen and carbon dioxide partial pressures were adjusted through alterations to the ventilator settings, resulting in changes to the cerebral blood volume and oxygenation. Three-dimensional images were generated using the physiologically associated differences in the measured data, obviating the need for data calibration using a separate reference measurement. The results exhibit large changes in absorption coefficient at both wavelengths. Images corresponding to differences in concentrations of oxy- and deoxyhaemoglobin are in qualitative agreement with known physiological data.

The effects of internal refractive index variation in near-infrared optical tomography: a finite element modelling approach

Near-infrared (NIR) tomography is a technique used to measure light propagation through tissue and generate images of internal optical property distributions from boundary measurements. Most popular applications have concentrated on female breast imaging, neonatal and adult head imaging, as well as muscle and small animal studies. In most instances a highly scattering medium with a homogeneous refractive index is assumed throughout the imaging domain. Using these assumptions, it is possible to simplify the model to the diffusion approximation. However, biological tissue contains regions of varying optical absorption and scatter, as well as varying refractive index. In this work, we introduce an internal boundary constraint in the finite element method approach to modelling light propagation through tissue that accounts for regions of different refractive indices. We have compared the results to data from a Monte Carlo simulation and show that for a simple two-layered slab model of varying refractive index, the phase of the measured reflectance data is significantly altered by the variation in internal refractive index, whereas the amplitude data are affected only slightly.

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.

A method for generating patient-specific finite element meshes for head modelling

Finite element modelling of fields within the body, whether electrical or optical, requires knowledge of the geometry of the object being examined. It can be clinically impractical to obtain accurate surface information for individual patients, although a limited set of measurements such as the locations of sensors attached to the body, can be acquired more readily. In this paper, we describe how a generic surface taken from an adult head is warped to fit points measured on a neonatal head surface to provide a new, individual surface from which a finite element mesh was generated. Simulations show that data generated from this mesh and from the original neonatal head surface are similar to within experimental errors. However, data generated from a mesh of the best fit sphere were significantly different from data generated from the original neonatal head surface.