Samuel Powell

Functional imaging of the human brain using a modular, fibre-less, high-density diffuse optical tomography system

We present the first three-dimensional, functional images of the human brain to be obtained using a fibre-less, high-density diffuse optical tomography system. Our technology consists of independent, miniaturized, silicone-encapsulated DOT modules that can be placed directly on the scalp. Four of these modules were arranged to provide up to 128, dual-wavelength measurement channels over a scalp area of approximately 60 × 65 mm2. Using a series of motor-cortex stimulation experiments, we demonstrate that this system can obtain high-quality, continuous-wave measurements at source-detector separations ranging from 14 to 55 mm in adults, in the presence of hair. We identify robust haemodynamic response functions in 5 out of 5 subjects, and present diffuse optical tomography images that depict functional haemodynamic responses that are well-localized in all three dimensions at both the individual and group levels. This prototype modular system paves the way for a new generation of wearable, wireless, high-density optical neuroimaging technologies.

Highly parallel Monte-Carlo simulations of the acousto-optic effect in heterogeneous turbid media

The development of a highly parallel simulation of the acousto-optic effect is detailed. The simulation supports optically heterogeneous simulation domains under insonification by arbitrary monochromatic ultrasound fields. An adjoint method for acousto-optics is proposed to permit point-source/point-detector simulations. The flexibility and efficiency of this simulation code is demonstrated in the development of spatial absorption sensitivity maps which are in broad agreement with current experimental investigations. The simulation code has the potential to provide guidance in the feasibility and optimization of future studies of the acousto-optic technique, and its speed may permit its use as part of an iterative inversion model.

Towards a wearable near infrared spectroscopic probe for monitoring concentrations of multiple chromophores in biological tissue in vivo

The first wearable multi-wavelength technology for functional near-infrared spectroscopy has been developed, based on a custom-built 8-wavelength light emitting diode (LED) source. A lightweight fibreless probe is designed to monitor changes in the concentrations of multiple absorbers (chromophores) in biological tissue, the most dominant of which at near-infrared wavelengths are oxyhemoglobin and deoxyhemoglobin. The use of multiple wavelengths enables signals due to the less dominant chromophores to be more easily distinguished from those due to hemoglobin and thus provides more complete and accurate information about tissue oxygenation, hemodynamics, and metabolism. The spectroscopic probe employs four photodiode detectors coupled to a four-channel charge-to-digital converter which includes a charge integration amplifier and an analogue-to-digital converter (ADC). Use of two parallel charge integrators per detector enables one to accumulate charge while the other is being read out by the ADC, thus facilitating continuous operation without dead time. The detector system has a dynamic range of about 80 dB. The customized source consists of eight LED dies attached to a 2 mm × 2 mm substrate and encapsulated in UV-cured epoxy resin. Switching between dies is performed every 20 ms, synchronized to the detector integration period to within 100 ns. The spectroscopic probe has been designed to be fully compatible with simultaneous electroencephalography measurements. Results are presented from measurements on a phantom and a functional brain activation study on an adult volunteer, and the performance of the spectroscopic probe is shown to be very similar to that of a benchtop broadband spectroscopy system. The multi-wavelength capabilities and portability of this spectroscopic probe will create significant opportunities for in vivo studies in a range of clinical and life science applications.

Reconstruction-classification method for quantitative photoacoustic tomography

Abstract. We propose a combined reconstruction-classification method for simultaneously recovering absorption and scattering in turbid media from images of absorbed optical energy. This method exploits knowledge that optical parameters are determined by a limited number of classes to iteratively improve their estimate. Numerical experiments show that the proposed approach allows for accurate recovery of absorption and scattering in two and three dimensions, and delivers superior image quality with respect to traditional reconstruction-only approaches.

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

Abstract. 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.

Linear reconstruction of absorption perturbations in coherent ultrasound-modulated optical tomography

Abstract. Discrete form sensitivity functions are derived for a number of measurement types employed in autocorrelation-based ultrasound-modulated optical tomography. The Jacobian for a particular problem is constructed from the set of such sensitivity functions arising from a raster scan of a focused ultrasound field through a turbid medium. A linear reconstruction of an absorption perturbation is performed in a simulated difference data experiment, and the performance of the measurement types is compared under different degrees of added noise.

Quantitative photoacoustic tomography using forward and adjoint Monte Carlo models of radiance

Abstract. Forward and adjoint Monte Carlo (MC) models of radiance are proposed for use in model-based quantitative photoacoustic tomography. A two-dimensional (2-D) radiance MC model using a harmonic angular basis is introduced and validated against analytic solutions for the radiance in heterogeneous media. A gradient-based optimization scheme is then used to recover 2-D absorption and scattering coefficients distributions from simulated photoacoustic measurements. It is shown that the functional gradients, which are a challenge to compute efficiently using MC models, can be calculated directly from the coefficients of the harmonic angular basis used in the forward and adjoint models. This work establishes a framework for transport-based quantitative photoacoustic tomography that can fully exploit emerging highly parallel computing architectures.

Evaluating real-time image reconstruction in diffuse optical tomography using physiologically realistic test data.

In diffuse optical tomography (DOT), real-time image reconstruction of oxy- and deoxy-haemoglobin changes occurring in the brain could give valuable information in clinical care settings. Although non-linear reconstruction techniques could provide more accurate results, their computational burden makes them unsuitable for real-time applications. Linear techniques can be employed under the assumption that the expected change in absorption is small. Several approaches exist, differing primarily in their handling of regularization and the noise statistics. In real experiments, it is impossible to compute the true noise statistics, because of the presence of physiological oscillations in the measured data. This is even more critical in real-time applications, where no off-line filtering and averaging can be performed to reduce the noise level. Therefore, many studies substitute the noise covariance matrix with the identity matrix. In this paper, we examined two questions: does using the noise model with realistic, imperfect data yield an improvement in image quality compared to using the identity matrix; and what is the difference in quality between online and offline reconstructions. Bespoke test data were created using a novel process through which simulated changes in absorption were added to real resting-state DOT data. A realistic multi-layer head model was used as the geometry for the reconstruction. Results validated our assumptions, highlighting the validity of computing the noise statistics from the measured data for online image reconstruction, which was performed at 2 Hz. Our results can be directly extended to a real application where real-time imaging is required.

Optimization of the acousto-optic signal detection in cylindrical geometry

The use of ultrasonic tagging of multiple scattered photons within turbid media for tomographic imaging is typically performed using optical detection in transmission mode. This study aimed to optimize the detection of the acousto-optic (AO) signal in cylindrical geometry, with a view to using the technique to measure blood oxygenation in the internal jugular vein of infants in the future. In our experiments, homogeneous phantoms of multiple transport scattering coefficients were constructed for the described geometry mimicking the infant neck. The optical source was systematically repositioned at different angles relative to the optical detector and the resulting AO signal was measured. The experimental results were also compared to focused ultrasound AO Monte Carlo (MC) simulation results. It was found that the optimal modulation depth and noise variance were highly dependent on the overlap region between the optical path length of the optical source-detector pair and the ultrasound focal zone. Therefore the optimal positions for both the optical and ultrasound probes could be estimated from both experimental and simulation results for a given geometry.

Image reconstruction of oxidized cerebral cytochrome C oxidase changes from broadband near-infrared spectroscopy data

Abstract. In diffuse optical tomography (DOT), overlapping and multidistance measurements are required to reconstruct depth-resolved images of oxy- (HbO2) and deoxy- (HHb) hemoglobin concentration changes occurring in the brain. These can be considered an indirect measure of brain activity, under the assumption of intact neurovascular coupling. Broadband systems also allow changes in the redox state of cytochrome c oxidase (oxCCO) to be measured, which can be an important biomarker when neurovascular coupling is impaired. We used DOT to reconstruct images of Δ[HbO2], Δ[HHb], and Δ[oxCCO] from data acquired with a broadband system. Four healthy volunteers were measured while performing a visual stimulation task (4-Hz inverting checkerboard). The broadband system was configured to allow multidistance and overlapping measurements of the participants’ visual cortex with 32 channels. A multispectral approach was employed to reconstruct changes in concentration of the three chromophores during the visual stimulation. A clear and focused activation was reconstructed in the left occipital cortex of all participants. The difference between the residuals of the three-chromophore model and of the two-chromophore model (recovering only Δ[HbO2] and Δ[HHb]) exhibits a spectrum similar to that of oxCCO. These results form a basis for further studies aimed to further optimize image reconstruction of Δ[oxCCO].