Quantitative evaluation of atlas-based high-density diffuse optical tomography for imaging of the human visual cortex

Effects of atlas-based anatomy on modelled light transport in the neonatal head

Objective.Diffuse optical tomography (DOT) provides a relatively convenient method for imaging haemodynamic changes related to neuronal activity on the cerebral cortex. Due to practical challenges in obtaining anatomical images of neonates, an anatomical framework is often created from an age-appropriate atlas model, which is individualized to the subject based on measurements of the head geometry. This work studies the approximation error arising from using an atlas instead of the neonate's own anatomical model.Approach.We consider numerical simulations of frequency-domain (FD) DOT using two approaches, Monte Carlo simulations and diffusion approximation via finite element method, and observe the variation in (1) the logarithm of amplitude and phase shift measurements, and (2) the corresponding inner head sensitivities (Jacobians), due to varying segmented anatomy. Varying segmentations are sampled by registering 165 atlas models from a neonatal database to the head geometry of one individual selected as the reference model. Prior to the registration, we refine the segmentation of the cerebrospinal fluid (CSF) by separating the CSF into two physiologically plausible layers.Main results.In absolute measurements, a considerable change in the grey matter or extracerebral tissue absorption coefficient was found detectable over the anatomical variation. In difference measurements, a small local 10%-increase in brain absorption was clearly detectable in the simulated measurements over the approximation error in the Jacobians, despite the wide range of brain maturation among the registered models.Significance.Individual-level atlas models could potentially be selected within several weeks in gestational age in DOT difference imaging, if an exactly age-appropriate atlas is not available. The approximation error method could potentially be implemented to improve the accuracy of atlas-based imaging. The presented CSF segmentation algorithm could be useful also in other model-based imaging modalities. The computation of FD Jacobians is now available in the widely-used Monte Carlo eXtreme software.

Intraoperative Optical and Fluorescence Imaging of Blood Flow Distributions in Mastectomy Skin Flaps for Identifying Ischemic Tissues

SUMMARY

Insufficient blood flow causes mastectomy skin flap necrosis in 5 to 30 percent of cases. Fluorescence angiography with the injection of indocyanine green dye has shown high sensitivities (90 to 100 percent) but moderate specificities (72 to 50 percent) in predicting mastectomy skin flap necrosis. However, a number of challenging issues limit its wide acceptance in clinical settings, including allergic reaction, short time-window for observation, and high cost for equipment and supplies. An emerging inexpensive speckle contrast diffuse correlation tomography technology enables noninvasive, noncontact, and continuous three-dimensional imaging of blood flow distributions in deep tissues. This preliminary study tested the hypothesis that speckle contrast diffuse correlation tomography and indocyanine green-fluorescence angiography measurements of blood flow distributions in mastectomy skin flaps are consistent. Eleven female patients undergoing skin-sparing or nipple-sparing mastectomies were imaged sequentially by the dye-free speckle contrast diffuse correlation tomography and dye-based commercial fluorescence angiography (SPY-PHI). Resulting images from these two imaging modalities were co-registered based on the ischemic areas with the lowest blood flow values. Because the ischemic areas have irregular shapes, a novel contour-based algorithm was used to compare three-dimensional images of blood flow distribution and two-dimensional maps of indocyanine green perfusion. Significant correlations were observed between the two measurements in all contours from a selected area of 10 × 10 mm 2 with the lowest blood flow ( r ≥ 0.78; p < 0.004), suggesting that speckle contrast diffuse correlation tomography provides the information for identifying ischemic tissues in mastectomy skin flaps. With further optimization and validation in large populations, speckle contrast diffuse correlation tomography may ultimately be used as a noninvasive and inexpensive imaging tool for intraoperative assessment of skin flap viability to predict mastectomy skin flap necrosis.

CLINICAL QUESTION/LEVEL OF EVIDENCE

Diagnostic, II.

Lawrence K, Tachtsidis I, Tong Y, Torricelli A, Urner T, Wabnitz H, Wolf M, Wolf U, Xu S, Yang C, Yodh AG, Yücel MA, Zhou W. Optical imaging and spectroscopy for the study of the human brain: status report

This report is the second part of a comprehensive two-part series aimed at reviewing an extensive and diverse toolkit of novel methods to explore brain health and function. While the first report focused on neurophotonic tools mostly applicable to animal studies, here, we highlight optical spectroscopy and imaging methods relevant to noninvasive human brain studies. We outline current state-of-the-art technologies and software advances, explore the most recent impact of these technologies on neuroscience and clinical applications, identify the areas where innovation is needed, and provide an outlook for the future directions.

Imaging affective and non-affective touch processing in two-year-old children

Touch is an important component of early parent-child interaction and plays a critical role in the socio-emotional development of children. However, there are limited studies on touch processing amongst children in the age range from one to three years. The present study used frequency-domain diffuse optical tomography (DOT) to investigate the processing of affective and non-affective touch over left frontotemporal brain areas contralateral to the stimulated forearm in two-year-old children. Affective touch was administered by a single stroke with a soft brush over the child's right dorsal forearm at 3 cm/s, while non-affective touch was provided by multiple brush strokes at 30 cm/s. We found that in the insula, the total haemoglobin (HbT) response to slow brushing was significantly greater than the response to fast brushing (slow > fast). Additionally, a region in the postcentral gyrus, Rolandic operculum and superior temporal gyrus exhibited greater response to fast brushing than slow brushing (fast > slow). These findings confirm that an adult-like pattern of haemodynamic responses to affective and non-affective touch can be recorded in two-year-old subjects using DOT. To improve the accuracy of modelling light transport in the two-year-old subjects, we used a published age-appropriate atlas and deformed it to match the exterior shape of each subject's head. We estimated the combined scalp and skull, and grey matter (GM) optical properties by fitting simulated data to calibrated and coupling error corrected phase and amplitude measurements. By utilizing a two-compartment cerebrospinal fluid (CSF) model, the accuracy of estimation of GM optical properties and the localization of activation in the insula was improved. The techniques presented in this paper can be used to study neural development of children at different ages and illustrate that the technology is well-tolerated by most two-year-old children and not excessively sensitive to subject movement. The study points the way towards exciting possibilities in functional imaging of deeper functional areas near sulci in small children.

Variability in atlas registration of optical intrinsic signal imaging and its effect on functional connectivity analysis

To compare neuroimaging data between subjects, images from individual sessions need to be aligned to a common reference or "atlas." Atlas registration of optical intrinsic signal imaging of mice, for example, is commonly performed using affine transforms with parameters determined by manual selection of canonical skull landmarks. Errors introduced by such procedures have not previously been investigated. We quantify the variability that arises from this process and consequent errors from misalignment that affect interpretation of functional neuroimaging data. We propose an improved method, using separately acquired high-resolution images and demonstrate improvements in variability and alignment using this method.

Global motion detection and censoring in high-density diffuse optical tomography

Motion‐induced artifacts can significantly corrupt optical neuroimaging, as in most neuroimaging modalities. For high‐density diffuse optical tomography (HD‐DOT) with hundreds to thousands of source‐detector pair measurements, motion detection methods are underdeveloped relative to both functional magnetic resonance imaging (fMRI) and standard functional near‐infrared spectroscopy (fNIRS). This limitation restricts the application of HD‐DOT in many challenging imaging situations and subject populations (e.g., bedside monitoring and children). Here, we evaluated a new motion detection method for multi‐channel optical imaging systems that leverages spatial patterns across measurement channels. Specifically, we introduced a global variance of temporal derivatives (GVTD) metric as a motion detection index. We showed that GVTD strongly correlates with external measures of motion and has high sensitivity and specificity to instructed motion—with an area under the receiver operator characteristic curve of 0.88, calculated based on five different types of instructed motion. Additionally, we showed that applying GVTD‐based motion censoring on both hearing words task and resting state HD‐DOT data with natural head motion results in an improved spatial similarity to fMRI mapping. We then compared the GVTD similarity scores with several commonly used motion correction methods described in the fNIRS literature, including correlation‐based signal improvement (CBSI), temporal derivative distribution repair (TDDR), wavelet filtering, and targeted principal component analysis (tPCA). We find that GVTD motion censoring on HD‐DOT data outperforms other methods and results in spatial maps more similar to those of matched fMRI data.

Brain segmentation, spatial censoring, and averaging techniques for optical functional connectivity imaging in mice

Resting-state functional connectivity analysis using optical neuroimaging holds the potential to be a powerful bridge between mouse models of disease and clinical neurologic monitoring. However, analysis techniques specific to optical methods are rudimentary, and algorithms from magnetic resonance imaging are not always applicable to optics. We have developed visual processing tools to increase data quality, improve brain segmentation, and average across sessions with better field-of-view. We demonstrate improved performance using resting-state optical intrinsic signal from normal mice. The proposed methods increase the amount of usable data from neuroimaging studies, improve image fidelity, and should be translatable to human optical neuroimaging systems.

Graph- and finite element-based total variation models for the inverse problem in diffuse optical tomography

Total variation (TV) is a powerful regularization method that has been widely applied in different imaging applications, but is difficult to apply to diffuse optical tomography (DOT) image reconstruction (inverse problem) due to unstructured discretization of complex geometries, non-linearity of the data fitting and regularization terms, and non-differentiability of the regularization term. We develop several approaches to overcome these difficulties by: i) defining discrete differential operators for TV regularization using both finite element and graph representations; ii) developing an optimization algorithm based on the alternating direction method of multipliers (ADMM) for the non-differentiable and non-linear minimization problem; iii) investigating isotropic and anisotropic variants of TV regularization, and comparing their finite element (FEM)- and graph-based implementations. These approaches are evaluated on experiments on simulated data and real data acquired from a tissue phantom. Our results show that both FEM and graph-based TV regularization is able to accurately reconstruct both sparse and non-sparse distributions without the over-smoothing effect of Tikhonov regularization and the over-sparsifying effect of L1 regularization. The graph representation was found to out-perform the FEM method for low-resolution meshes, and the FEM method was found to be more accurate for high-resolution meshes.

High-density diffuse optical tomography for imaging human brain function

This review describes the unique opportunities and challenges for noninvasive optical mapping of human brain function. Diffuse optical methods offer safe, portable, and radiation free alternatives to traditional technologies like positron emission tomography or functional magnetic resonance imaging (fMRI). Recent developments in high-density diffuse optical tomography (HD-DOT) have demonstrated capabilities for mapping human cortical brain function over an extended field of view with image quality approaching that of fMRI. In this review, we cover fundamental principles of the diffusion of near infrared light in biological tissue. We discuss the challenges involved in the HD-DOT system design and implementation that must be overcome to acquire the signal-to-noise necessary to measure and locate brain function at the depth of the cortex. We discuss strategies for validation of the sensitivity, specificity, and reliability of HD-DOT acquired maps of cortical brain function. We then provide a brief overview of some clinical applications of HD-DOT. Though diffuse optical measurements of neurophysiology have existed for several decades, tremendous opportunity remains to advance optical imaging of brain function to address a crucial niche in basic and clinical neuroscience: that of bedside and minimally constrained high fidelity imaging of brain function.

A Kalman-based tomographic scheme for directly reconstructing activation levels of brain function

In functional near-infrared spectroscopy (fNIRS), the conventional indirect approaches first separately recover the spatial distribution of the changes in the optical properties at every time point, and then extract the activation levels by a time-course analysis process at every site. In the tomographic implementation of fNIRS, i.e., diffuse optical tomography (DOT), these approaches not only suffer from the ill-posedness of the optical inversions and error propagation between the two successive steps, but also fail to achieve satisfactory temporal resolution due to the requirement for a complete data set. To cope with the above adversities of the indirect approaches, we propose herein a direct approach to tomographically reconstructing the activation levels by incorporating a Kalman scheme. Dynamic simulative and phantom experiments were conducted for the performance validation of the proposed approach, demonstrating its potentials to improve the calculated images and to relax the speed limitation of the instruments.

Lightweight sCMOS-based high-density diffuse optical tomography

Though optical imaging of human brain function is gaining momentum, widespread adoption is restricted in part by a tradeoff among cap wearability, field of view, and resolution. To increase coverage while maintaining functional magnetic resonance imaging (fMRI)-comparable image quality, optical systems require more fibers. However, these modifications drastically reduce the wearability of the imaging cap. The primary obstacle to optimizing wearability is cap weight, which is largely determined by fiber diameter. Smaller fibers collect less light and lead to challenges in obtaining adequate signal-to-noise ratio. Here, we report on a design that leverages the exquisite sensitivity of scientific CMOS cameras to use fibers with ∼30× smaller cross-sectional area than current high-density diffuse optical tomography (HD-DOT) systems. This superpixel sCMOS DOT (SP-DOT) system uses 200-μm -diameter fibers that facilitate a lightweight, wearable cap. We developed a superpixel algorithm with pixel binning and electronic noise subtraction to provide high dynamic range ( >105 ), high frame rate ( >6  Hz ), and a low effective detectivity threshold ( ∼200  fW/Hz1/2-mm2 ), each comparable with previous HD-DOT systems. To assess system performance, we present retinotopic mapping of the visual cortex ( n=5 subjects). SP-DOT offers a practical solution to providing a wearable, large field-of-view, and high-resolution optical neuroimaging system.

L1-norm based nonlinear reconstruction improves quantitative accuracy of spectral diffuse optical tomography

Spectrally constrained diffuse optical tomography (SCDOT) is known to improve reconstruction in diffuse optical imaging; constraining the reconstruction by coupling the optical properties across multiple wavelengths suppresses artefacts in the resulting reconstructed images. In other work, L1-norm regularization has been shown to improve certain types of image reconstruction problems as its sparsity-promoting properties render it robust against noise and enable the preservation of edges in images, but because the L1-norm is non-differentiable, it is not always simple to implement. In this work, we show how to incorporate L1 regularization into SCDOT. Three popular algorithms for L1 regularization are assessed for application in SCDOT: iteratively reweighted least square algorithm (IRLS), alternating directional method of multipliers (ADMM), and fast iterative shrinkage-thresholding algorithm (FISTA). We introduce an objective procedure for determining the regularization parameter in these algorithms and compare their performance in simulated experiments, and in real data acquired from a tissue phantom. Our results show that L1 regularization consistently outperforms Tikhonov regularization in this application, particularly in the presence of noise.

Toward real-time diffuse optical tomography: accelerating light propagation modeling employing parallel computing on GPU and CPU

Parameter recovery in diffuse optical tomography is a computationally expensive algorithm, especially when used for large and complex volumes, as in the case of human brain functional imaging. The modeling of light propagation, also known as the forward problem, is the computational bottleneck of the recovery algorithm, whereby the lack of a real-time solution is impeding practical and clinical applications. The objective of this work is the acceleration of the forward model, within a diffusion approximation-based finite-element modeling framework, employing parallelization to expedite the calculation of light propagation in realistic adult head models. The proposed methodology is applicable for modeling both continuous wave and frequency-domain systems with the results demonstrating a 10-fold speed increase when GPU architectures are available, while maintaining high accuracy. It is shown that, for a very high-resolution finite-element model of the adult human head with ∼600,000 nodes, consisting of heterogeneous layers, light propagation can be calculated at ∼0.25  s/excitation source.

Probabilistic liver atlas construction

BACKGROUND

Anatomical atlases are 3D volumes or shapes representing an organ or structure of the human body. They contain either the prototypical shape of the object of interest together with other shapes representing its statistical variations (statistical atlas) or a probability map of belonging to the object (probabilistic atlas). Probabilistic atlases are mostly built with simple estimations only involving the data at each spatial location.

RESULTS

A new method for probabilistic atlas construction that uses a generalized linear model is proposed. This method aims to improve the estimation of the probability to be covered by the liver. Furthermore, all methods to build an atlas involve previous coregistration of the sample of shapes available. The influence of the geometrical transformation adopted for registration in the quality of the final atlas has not been sufficiently investigated. The ability of an atlas to adapt to a new case is one of the most important quality criteria that should be taken into account. The presented experiments show that some methods for atlas construction are severely affected by the previous coregistration step.

CONCLUSION

We show the good performance of the new approach. Furthermore, results suggest that extremely flexible registration methods are not always beneficial, since they can reduce the variability of the atlas and hence its ability to give sensible values of probability when used as an aid in segmentation of new cases.

Region-based diffuse optical tomography with registered atlas: in vivo acquisition of mouse optical properties

The reconstruction quality in the model-based optical tomography modalities can greatly benefit from a priori information of accurate tissue optical properties, which are difficult to be obtained in vivo with a conventional diffuse optical tomography (DOT) system alone. One of the solutions is to apply a priori anatomical structures obtained with anatomical imaging systems such as X-ray computed tomography (XCT) to constrain the reconstruction process of DOT. However, since X-ray offers low soft-tissue contrast, segmentation of abdominal organs from sole XCT images can be problematic. In order to overcome the challenges, the current study proposes a novel method of recovering a priori organ-oriented tissue optical properties, where anatomical structures of an in vivo mouse are approximately obtained by registering a standard anatomical atlas, i.e., the Digimouse, to the target XCT volume with the non-rigid image registration, and, in turn, employed to guide DOT for extracting the optical properties of inner organs. Simulative investigations have validated the methodological availability of such atlas-registration-based DOT strategy in revealing both a priori anatomical structures and optical properties. Further experiments have demonstrated the feasibility of the proposed method for acquiring the organ-oriented tissue optical properties of in vivo mice, making it as an efficient way of the reconstruction enhancement.

Automated voxel classification used with atlas-guided diffuse optical tomography for assessment of functional brain networks in young and older adults

Atlas-guided diffuse optical tomography (atlas-DOT) is a computational means to image changes in cortical hemodynamic signals during human brain activities. Graph theory analysis (GTA) is a network analysis tool commonly used in functional neuroimaging to study brain networks. Atlas-DOT has not been analyzed with GTA to derive large-scale brain connectivity/networks based on near-infrared spectroscopy (NIRS) measurements. We introduced an automated voxel classification (AVC) method that facilitated the use of GTA with atlas-DOT images by grouping unequal-sized finite element voxels into anatomically meaningful regions of interest within the human brain. The overall approach included volume segmentation, AVC, and cross-correlation. To demonstrate the usefulness of AVC, we applied reproducibility analysis to resting-state functional connectivity measurements conducted from 15 young adults in a two-week period. We also quantified and compared changes in several brain network metrics between young and older adults, which were in agreement with those reported by a previous positron emission tomography study. Overall, this study demonstrated that AVC is a useful means for facilitating integration or combination of atlas-DOT with GTA and thus for quantifying NIRS-based, voxel-wise resting-state functional brain networks.

Improving the quantitative accuracy of cerebral oxygen saturation in monitoring the injured brain using atlas based Near Infrared Spectroscopy models

The application of Near Infrared Spectroscopy (NIRS) for the monitoring of the cerebral oxygen saturation within the brain is well established, albeit using temporal data that can only measure relative changes of oxygenation state of the brain from a baseline. The focus of this investigation is to demonstrate that hybridisation of existing near infrared probe designs and reconstruction techniques can pave the way to produce a system and methods that can be used to monitor the absolute oxygen saturation in the injured brain. Using registered Atlas models in simulation, a novel method is outlined by which the quantitative accuracy and practicality of NIRS for specific use in monitoring the injured brain, can be improved, with cerebral saturation being recovered to within 10.1 ± 1.8% of the expected values.

Fabrication and application of heterogeneous printed mouse phantoms for whole animal optical imaging

This work demonstrates the usefulness of 3D printing for optical imaging applications. Progress in developing optical imaging for biomedical applications requires customizable and often complex objects for testing and evaluation. There is therefore high demand for what have become known as tissue-simulating "phantoms." We present a new optical phantom fabricated using inexpensive 3D printing methods with multiple materials, allowing for the placement of complex inhomogeneities in complex or anatomically realistic geometries, as opposed to previous phantoms, which were limited to simple shapes formed by molds or machining. We use diffuse optical imaging to reconstruct optical parameters in 3D space within a printed mouse to show the applicability of the phantoms for developing whole animal optical imaging methods. This phantom fabrication approach is versatile, can be applied to optical imaging methods besides diffusive imaging, and can be used in the calibration of live animal imaging data.

Experimental validation of a high-resolution diffuse optical imaging modality: photomagnetic imaging

We present experimental results that validate our imaging technique termed photomagnetic imaging (PMI). PMI illuminates the medium under investigation with a near-infrared light and measures the induced temperature increase using magnetic resonance imaging. A multiphysics solver combining light and heat propagation is used to model spatiotemporal distribution of temperature increase. Furthermore, a dedicated PMI reconstruction algorithm has been developed to reveal high-resolution optical absorption maps from temperature measurements. Being able to perform measurements at any point within the medium, PMI overcomes the limitations of conventional diffuse optical imaging. We present experimental results obtained on agarose phantoms mimicking biological tissue with inclusions having either different sizes or absorption contrasts, located at various depths. The reconstructed images show that PMI can successfully resolve these inclusions with high resolution and recover their absorption coefficient with high-quantitative accuracy. Even a 1-mm inclusion located 6-mm deep is recovered successfully and its absorption coefficient is underestimated by only 32%. The improved PMI system presented here successfully operates under the maximum skin exposure limits defined by the American National Standards Institute, which opens up the exciting possibility of its future clinical use for diagnostic purposes.

Fast and efficient image reconstruction for high density diffuse optical imaging of the human brain

Real-time imaging of human brain has become an important technique within neuroimaging. In this study, a fast and efficient sensitivity map generation based on Finite Element Models (FEM) is developed which utilises a reduced sensitivitys matrix taking advantage of sparsity and parallelisation processes. Time and memory efficiency of these processes are evaluated and compared with conventional method showing that for a range of mesh densities from 50000 to 320000 nodes, the required memory is reduced over tenfold and computational time fourfold allowing for near real-time image recovery.

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.

Evaluation of rigid registration methods for whole head imaging in diffuse optical tomography

Functional brain imaging has become an important neuroimaging technique for the study of brain organization and development. Compared to other imaging techniques, diffuse optical tomography (DOT) is a portable and low-cost technique that can be applied to infants and hospitalized patients using an atlas-based light model. For DOT imaging, the accuracy of the forward model has a direct effect on the resulting recovered brain function within a field of view and so the accuracy of the spatially normalized atlas-based forward models must be evaluated. Herein, the accuracy of atlas-based DOT is evaluated on models that are spatially normalized via a number of different rigid registration methods on 24 subjects. A multileveled approach is developed to evaluate the correlation of the geometrical and sensitivity accuracies across the full field of view as well as within specific functional subregions. Results demonstrate that different registration methods are optimal for recovery of different sets of functional brain regions. However, the "nearest point to point" registration method, based on the EEG 19 landmark system, is shown to be the most appropriate registration method for image quality throughout the field of view of the high-density cap that covers the whole of the optically accessible cortex.