Improving the diffuse optical imaging spatial resolution of the cerebral hemodynamic response to brain activation in humans

Somatosensory activation of two fingers can be discriminated with ultrahigh-density diffuse optical tomography

Topographic non-invasive near infrared spectroscopy (NIRS) has become a well-established tool for functional brain imaging. Applying up to 100 optodes over the head of a subject, allows achieving a spatial resolution in the centimeter range. This resolution is poor compared to other functional imaging tools. However, recently it was shown that diffuse optical tomography (DOT) as an extension of NIRS based on high-density (HD) probe arrays and supplemented by an advanced image reconstruction procedure allows describing activation patterns with a spatial resolution in the millimeter range. Building on these findings, we hypothesize that HD-DOT may render very focal activations accessible which would be missed by the traditionally used sparse arrays. We examined activation patterns in the primary somatosensory cortex, since its somatotopic organization is very fine-grained. We performed a vibrotactile stimulation study of the first and fifth finger in eight human subjects, using a 900-channel continuous-wave DOT imaging system for achieving a higher resolution than conventional topographic NIRS. To compare the results to a well-established high-resolution imaging technique, the same paradigm was investigated in the same subjects by means of functional magnetic resonance imaging (fMRI). In this work, we tested the advantage of ultrahigh-density probe arrays and show that highly focal activations would be missed by classical next-nearest neighbor NIRS approach, but also by DOT, when using a sparse probe array. Distinct activation patterns for both fingers correlated well with the expected neuroanatomy in five of eight subjects. Additionally we show that activation for different fingers is projected to different tissue depths in the DOT image. Comparison to the fMRI data yielded similar activation foci in seven out of ten finger representations in these five subjects when comparing the lateral localization of DOT and fMRI results.

Implicit and explicit prior information in near-infrared spectral imaging: accuracy, quantification and diagnostic value

Near-infrared spectroscopy (NIRS) of tissue provides quantification of absorbers, scattering and luminescent agents in bulk tissue through the use of measurement data and assumptions. Prior knowledge can be critical about things such as (i) the tissue shape and/or structure, (ii) spectral constituents, (iii) limits on parameters, (iv) demographic or biomarker data, and (v) biophysical models of the temporal signal shapes. A general framework of NIRS imaging with prior information is presented, showing that prior information datasets could be incorporated at any step in the NIRS process, with the general workflow being: (i) data acquisition, (ii) pre-processing, (iii) forward model, (iv) inversion/reconstruction, (v) post-processing, and (vi) interpretation/diagnosis. Most of the development in NIRS has used ad hoc or empirical implementations of prior information such as pre-measured absorber or fluorophore spectra, or tissue shapes as estimated by additional imaging tools. A comprehensive analysis would examine what prior information maximizes the accuracy in recovery and value for medical diagnosis, when implemented at separate stages of the NIRS sequence. Individual applications of prior information can show increases in accuracy or improved ability to estimate biochemical features of tissue, while other approaches may not. Most beneficial inclusion of prior information has been in the inversion/reconstruction process, because it solves the mathematical intractability. However, it is not clear that this is always the most beneficial stage.

Optical detection of brain function: simultaneous imaging of cerebral vascular response, tissue metabolism, and cellular activity in vivo

It is known that a remaining challenge for functional brain imaging is to distinguish the coupling and decoupling effects among neuronal activity, cerebral metabolism, and vascular hemodynamics, which highlights the need for new tools to enable simultaneous measures of these three properties in vivo. Here, we review current neuroimaging techniques and their prospects and potential limitations for tackling this challenge. We then report a novel dual-wavelength laser speckle imaging (DW-LSI) tool developed in our labs that enables simultaneous imaging of cerebral blood flow (CBF), cerebral blood volume, and tissue hemoglobin oxygenation, which allows us to monitor neurovascular and tissue metabolic activities at high spatiotemporal resolutions over a relatively large field of view. Moreover, we report digital frequency ramping Doppler optical coherence tomography (DFR-OCT) that allows for quantitative 3D imaging of the CBF network in vivo. In parallel, we review calcium imaging techniques to track neuronal activity, including intracellular calcium approach using Rhod2 fluorescence technique that we develop to detect neuronal activity in vivo. We report a new multimodality imaging platform that combines DW-LSI, DFR-OCT, and calcium fluorescence imaging for simultaneous detection of cortical hemodynamics, cerebral metabolism, and neuronal activities of the animal brain in vivo, as well as its integration with microprobes for imaging neuronal function in deep brain regions in vivo. Promising results of in vivo animal brain functional studies suggest the potential of this multimodality approach for future awake animal and behavioral studies.

Functional near infrared spectroscopy study of age-related difference in cortical activation patterns during cycling with speed feedback

Functional decline of lower-limb affects the ability of locomotion and the age-related brain differences have been elucidated among the elderly. Cycling exercise is a common training program for restoring motor function in the deconditioned elderly or stroke patients. The provision of speed feedback has been commonly suggested to clinical therapists for facilitating learning of controlled cycling performance and maintaining motivation in training programs with elderly participants. However, the cortical control of pedaling movements and the effect of external feedback remain poorly understanding. This study investigated the regional cortical activities detected by functional near infrared spectroscopy (fNIRS) in 12 healthy young and 13 healthy elderly subjects under conditions of cycling without-(free cycling) and with feedback (target cycling). The elderly exhibited predominant activation of the sensorimotor cortex during free cycling similar to young subjects but with poorer cycling performance. The cycling performance improved in both groups, and the elderly showed increased brain activities of the supplementary motor area and premotor cortex under target cycling condition. These findings demonstrated age-related changes in the cortical control in processing external feedback and pedaling movements. Use of fNIRS to evaluate brain activation patterns after training may facilitate brain-based design of tailored therapeutic rehabilitation strategies.

Contrast enhanced high-resolution diffuse optical tomography of the human brain using ICG

Non-invasive diffuse optical tomography (DOT) of the adult brain has recently been shown to improve the spatial resolution for functional brain imaging applications. Here we show that high-resolution (HR) DOT is also advantageous for clinical perfusion imaging using an optical contrast agent. We present the first HR-DOT results with a continuous wave near infrared spectroscopy setup using a dense grid of optical fibers and indocyanine green (ICG) as an exogenic contrast agent. We find an early arrival of the ICG bolus in the intracerebral tissue and a delayed arrival of the bolus in the extracerebral tissue, achieving the separation of both layers. This demonstrates the method's potential for brain perfusion monitoring in neurointensive care patients.

Quantitative photoacoustic imaging: correcting for heterogeneous light fluence distributions using diffuse optical tomography

The specificity of molecular and functional photoacoustic (PA) images depends on the accuracy of the photoacoustic absorption spectroscopy. The PA signal is proportional to the product of the optical absorption coefficient and local light fluence; quantitative PA measurements of the optical absorption coefficient therefore require an accurate estimation of optical fluence. Light-modeling aided by diffuse optical tomography (DOT) can be used to map the required fluence and to reduce errors in traditional PA spectroscopic analysis. As a proof-of-concept, we designed a tissue-mimicking phantom to demonstrate how fluence-related artifacts in PA images can lead to misrepresentations of tissue properties. To correct for these inaccuracies, the internal fluence in the tissue phantom was estimated by using DOT to reconstruct spatial distributions of the absorption and reduced scattering coefficients of multiple targets within the phantom. The derived fluence map, which only consisted of low spatial frequency components, was used to correct PA images of the phantom. Once calibrated to a known absorber, this method reduced errors in estimated absorption coefficients from 33% to 6%. These results experimentally demonstrate that combining DOT with PA imaging can significantly reduce fluence-related errors in PA images, while producing quantitatively accurate, high-resolution images of the optical absorption coefficient.

Enhanced functional brain imaging by using adaptive filtering and a depth compensation algorithm in diffuse optical tomography

Reflectance diffuse optical tomography (rDOT) of brain function is limited by its high sensitivity to the superficial tissues (i.e., the scalp and skull) and by its severe decrease in measurement sensitivity with increased depth. Significant interference in rDOT results from spontaneous fluctuations that are embedded in both the superficial tissues and brain, such as arterial pulsation and vasomotion. In this study, first we investigate coherence and phase shift of the spontaneous fluctuations in the resting state, within the superficial tissues and at various depths of the brain, respectively. We demonstrate that the spontaneous fluctuations originating from arterial pulsations ( ∼ 1 Hz) are spatially global and temporally coherent, while the fluctuations originating from vasomotion ( ∼ 0.1 Hz) tend to have less coherence with increased depth. Second, adaptive cancellation of spontaneous fluctuations with a frequency-specific strategy is utilized and validated in both resting and activation (evoked by a finger-tapping task) states. Third, improved depth localization of motor activation in reconstructed rDOT images is achieved by combining adaptive cancellation with a depth compensation algorithm that we recently reported.

Quantification of functional near infrared spectroscopy to assess cortical reorganization in children with cerebral palsy

Cerebral palsy (CP) is the most common motor disorder in children. Currently available neuroimaging techniques require complete body confinement and steadiness and thus are extremely difficult for pediatric patients. Here, we report the use and quantification of functional near infrared spectroscopy (fNIRS) to investigate the functional reorganization of the sensorimotor cortex in children with hemiparetic CP. Ten of sixteen children with congenital hemiparesis were measured during finger tapping tasks and compared with eight of sixteen age-matched healthy children, with an overall measurement success rate of 60%. Spatiotemporal analysis was introduced to quantify the motor activation and brain laterality. Such a quantitative approach reveals a consistent, contralateral motor activation in healthy children at 7 years of age or older. In sharp contrast, children with congenital hemiparesis exhibit all three of contralateral, bilateral and ipsilateral motor activations, depending on specific ages of the pediatric subjects. This study clearly demonstrates the feasibility of fNIRS to be utilized for investigating cortical reorganization in children with CP or other cortical disorders.

Comprehensive investigation of three-dimensional diffuse optical tomography with depth compensation algorithm

A depth compensation algorithm (DCA) can effectively improve the depth localization of diffuse optical tomography (DOT) by compensating the exponentially decreased sensitivity in the deep tissue. In this study, DCA is investigated based on computer simulations, tissue phantom experiments, and human brain imaging. The simulations show that DCA can largely improve the spatial resolution of DOT in addition to the depth localization, and DCA is also effective for multispectral DOT with a wide range of optical properties in the background tissue. The laboratory phantom experiment demonstrates that DCA can effectively differentiate two embedded objects at different depths in the medium. DCA is further validated by human brain imaging using a finger-tapping task. To our knowledge, this is the first demonstration to show that DCA is capable of accurately localizing cortical activations in the human brain in three dimensions.

High-resolution optical functional mapping of the human somatosensory cortex

Non-invasive optical imaging of brain function has been promoted in a number of fields in which functional magnetic resonance imaging (fMRI) is limited due to constraints induced by the scanning environment. Beyond physiological and psychological research, bedside monitoring and neurorehabilitation may be relevant clinical applications that are yet little explored. A major obstacle to advocate the tool in clinical research is insufficient spatial resolution. Based on a multi-distance high-density optical imaging setup, we here demonstrate a dramatic increase in sensitivity of the method. We show that optical imaging allows for the differentiation between activations of single finger representations in the primary somatosensory cortex (SI). Methodologically our findings confirm results in a pioneering study by Zeff et al. (2007) and extend them to the homuncular organization of SI. After performing a motor task, eight subjects underwent vibrotactile stimulation of the little finger and the thumb. We used a high-density diffuse-optical sensing array in conjunction with optical tomographic reconstruction. Optical imaging disclosed three discrete activation foci one for motor and two discrete foci for vibrotactile stimulation of the first and fifth finger, respectively. The results were co-registered to the individual anatomical brain anatomy (MRI) which confirmed the localization in the expected cortical gyri in four subjects. This advance in spatial resolution opens new perspectives to apply optical imaging in the research on plasticity notably in patients undergoing neurorehabilitation.

Quantitative evaluation of high-density diffuse optical tomography: in vivo resolution and mapping performance

Despite the unique brain imaging capabilities and advantages of functional near-infrared spectroscopy (fNIRS), including portability and comprehensive hemodynamic measurement, widespread acceptance in the neuroimaging community has been hampered by low spatial resolution and image localization errors. While recent technical developments such as high-density diffuse optical tomography (HD-DOT) have, in principle, been shown to have superior in silico image quality, the majority of optical imaging studies are still conducted with sparse fNIRS arrays, perhaps partially because the performance increases of HD-DOT appear incremental. Without a quantitative comparative analysis between HD-DOT and fNIRS, using both simulation and in vivo neuroimaging, the implications of the new HD-DOT technology have been difficult to judge. We present a quantitative comparison of HD-DOT and two commonly used fNIRS geometries using (1) standard metrics of image quality, (2) simulated brain mapping tasks, and (3) in vivo visual cortex mapping results in adult humans. The results show that better resolution and lower positional errors are achieved with HD-DOT and that these improvements provide a substantial advancement in neuroimaging capability. In particular, we demonstrate that HD-DOT enables detailed phase-encoded retinotopic mapping, while sparse arrays are limited to imaging individual block-design visual stimuli.

Image reconstruction in diffuse optical tomography based on simplified spherical harmonics approximation

The use of higher order approximations to the Radiative transport equation, through simplified spherical harmonics expansion (SP(N)) in optical tomography are presented. It is shown that, although the anisotropy factor can be modeled in the forward problem, its sensitivity to the measured boundary data is limited to superficial regions and more importantly, due to uniqueness of the inverse problem it cannot be determined using frequency domain data. Image reconstruction through the use of higher ordered models is presented. It is demonstrated that at higher orders (for example SP7) the image reconstruction becomes highly under-determined due to the large increase in the number of unknowns which cannot be adequately recovered. However, reconstruction of diffuse parameters, namely optical absorption and reduced scatter have shown to be more accurate where only the sensitivity matrix used in the inverse problem is based on SP(N) method and image reconstruction is limited to these two diffuse parameters.

Topographic localization of brain activation in diffuse optical imaging using spherical wavelets

Diffuse optical imaging is a non-invasive technique that uses near-infrared light to measure changes in brain activity through an array of sensors placed on the surface of the head. Compared to functional MRI, optical imaging has the advantage of being portable while offering the ability to record functional changes in both oxy- and deoxy-hemoglobin within the brain at a high temporal resolution. However, the reconstruction of accurate spatial images of brain activity from optical measurements represents an ill-posed and underdetermined problem that requires regularization. These reconstructions benefit from incorporating prior information about the underlying spatial structure and function of the brain. In this work, we describe a novel image reconstruction approach which uses surface-based wavelets derived from structural MRI to incorporate high-resolution anatomical and structural prior information about the brain. This surface-based approach is used to approximate brain activation patterns through the reconstruction and presentation of topographical (two-dimensional) maps of brain activation directly onto the folded surface of the cortex. The set of wavelet coefficients is directly estimated by a truncated singular-value decomposition based pseudo-inversion of the wavelet projection of the optical forward model. We use a reconstruction metric based on Shannon entropy which quantifies the sparse loading of the wavelet coefficients and is used to determine the optimal truncation and regularization of this inverse model. In this work, examples of the performance of this model are illustrated for several cases of numerical simulation and experimental data with comparison to functional magnetic resonance imaging.

A quantitative assessment of the depth sensitivity of an optical topography system using a solid dynamic tissue-phantom

A solid dynamic phantom with tissue-like optical properties is presented, which contains seven discrete targets impregnated with thermochromic pigment located at different depths from the surface. Changes in absorption are obtained in response to localized heating of the targets, simulating haemodynamic changes occurring in the brain and other tissues. The depth sensitivity of a continuous wave optical topography system was assessed successfully using the phantom. Images of the targets have been reconstructed using a spatially variant regularization, and the determined spatial localization in the depth direction is shown to be accurate within an uncertainty of about 3 mm down to a depth of about 30 mm.

Comparison of principal and independent component analysis in removing extracerebral interference from near-infrared spectroscopy signals

Near-infrared spectroscopy (NIRS) is a method for noninvasive estimation of cerebral hemodynamic changes. Principal component analysis (PCA) and independent component analysis (ICA) can be used for decomposing a set of signals to underlying components. Our objective is to determine whether PCA or ICA is more efficient in identifying and removing scalp blood flow interference from multichannel NIRS signals. Concentration changes of oxygenated (HbO(2)) and deoxygenated (HbR) hemoglobin are measured on the forehead with multichannel NIRS during hyper- and hypocapnia. PCA and ICA are used separately to identify and remove signal contribution from extracerebral tissue, and the resulting estimates of cerebral responses are compared to the expected cerebral responses. Both methods were able to reduce extracerebral contribution to the signals, but PCA typically performs equal to or better than ICA. The improvement in 3-cm signal quality achieved with both methods is comparable to increasing the source-detector separation from 3 to 5 cm. Especially PCA appears to be well suited for use in NIRS applications where the cerebral activation is diffuse, such as monitoring of global cerebral oxygenation and hemodynamics. Performance differences between PCA and ICA could be attributed primarily to different criteria for identifying the surface effect.

Investigation of the prefrontal cortex in response to duration-variable anagram tasks using functional near-infrared spectroscopy

We hypothesize that nonlinearity between short-term anagram tasks and corresponding hemodynamic responses can be observed by functional near-infrared spectroscopy (fNIRS) in the prefrontal cortex (PFC). The PFC of six human subjects in response to anagram tasks is investigated using multichannel fNIRS. Concentration changes of oxyhemoglobin and deoxyhemoglobin in the PFC are measured with variable anagram durations and at two difficulty levels (four- and six-letter anagrams). The durations to perform the selected anagram tasks range from several seconds to more than one minute. The dorsolateral PFC areas exhibit consistent and strong hemodynamic deactivation during and shortly after task execution. The superposition principle of a linear system is employed to investigate nonlinear hemodynamic features among three task duration subgroups: D1 = 2.0 sec, D2 = 4.0 sec, and D3 = 8.0 sec. Such analysis shows clear nonlinearity in hemodynamic responses on the PFC with task durations shorter than 4 sec. Our observation of significant deactivation in early hemodynamic responses in the PFC is consistent with multiple fNIRS studies and several reports given in the field of functional magnetic resonance imaging. A better understanding of nonlinearity in fNIRS signals will have potential for us to investigate brain adaptation and to extrapolate neuronal activities from hemodynamic signals.

Diffuse optical imaging

Diffuse optical imaging is a medical imaging technique that is beginning to move from the laboratory to the hospital. It is a natural extension of near-infrared spectroscopy (NIRS), which is now used in certain niche applications clinically and particularly for physiological and psychological research. Optical imaging uses sophisticated image reconstruction techniques to generate images from multiple NIRS measurements. The two main clinical applications--functional brain imaging and imaging for breast cancer--are reviewed in some detail, followed by a discussion of other issues such as imaging small animals and multimodality imaging. We aim to review the state of the art of optical imaging.

Numerical modelling and image reconstruction in diffuse optical tomography

The development of diffuse optical tomography as a functional imaging modality has relied largely on the use of model-based image reconstruction. The recovery of optical parameters from boundary measurements of light propagation within tissue is inherently a difficult one, because the problem is nonlinear, ill-posed and ill-conditioned. Additionally, although the measured near-infrared signals of light transmission through tissue provide high imaging contrast, the reconstructed images suffer from poor spatial resolution due to the diffuse propagation of light in biological tissue. The application of model-based image reconstruction is reviewed in this paper, together with a numerical modelling approach to light propagation in tissue as well as generalized image reconstruction using boundary data. A comprehensive review and details of the basis for using spatial and structural prior information are also discussed, whereby the use of spectral and dual-modality systems can improve contrast and spatial resolution.

Phase-encoded retinotopy as an evaluation of diffuse optical neuroimaging

Optical techniques enable portable, non-invasive functional neuroimaging. However, low lateral resolution and poor discrimination between brain hemodynamics and systemic contaminants have hampered the translation of near infrared spectroscopy from research instrument to widespread neuroscience tool. In this paper, we demonstrate that improvements in spatial resolution and signal-to-noise, afforded by recently developed high-density diffuse optical tomography approaches, now permit detailed phase-encoded mapping of the visual cortex's retinotopic organization. Due to its highly organized structure, the visual cortex has long served as a benchmark for judging neuroimaging techniques, including the original development of functional magnetic resonance imaging (fMRI) and positron emission tomography. Using phase-encoded visual stimuli that create traveling waves of cortical activations, we are able to discriminate the representations of multiple visual angles and eccentricities within an individual hemisphere, reproducing classic fMRI results. High contrast-to-noise and repeatable imaging allow the detection of inter-subject differences. These results represent a significant advancement in the level of detail that can be obtained from non-invasive optical imaging of functional brain responses. In addition, these phase-encoded paradigms and the maps they generate form a standardized model with which to judge new developments in optical algorithms and systems, such as new image reconstruction techniques and registration with anatomic imaging. With these advances in techniques and validation paradigms, optical neuroimaging can be extended into studies of higher-order brain function and of clinical utility with greater performance and confidence.

Near infrared optical tomography using NIRFAST: Algorithm for numerical model and image reconstruction

Diffuse optical tomography, also known as near infrared tomography, has been under investigation, for non-invasive functional imaging of tissue, specifically for the detection and characterization of breast cancer or other soft tissue lesions. Much work has been carried out for accurate modeling and image reconstruction from clinical data. NIRFAST, a modeling and image reconstruction package has been developed, which is capable of single wavelength and multi-wavelength optical or functional imaging from measured data. The theory behind the modeling techniques as well as the image reconstruction algorithms is presented here, and 2D and 3D examples are presented to demonstrate its capabilities. The results show that 3D modeling can be combined with measured data from multiple wavelengths to reconstruct chromophore concentrations within the tissue. Additionally it is possible to recover scattering spectra, resulting from the dominant Mie-type scatter present in tissue. Overall, this paper gives a comprehensive over view of the modeling techniques used in diffuse optical tomographic imaging, in the context of NIRFAST software package.

Real-time imaging of human brain function by near-infrared spectroscopy using an adaptive general linear model

Near-infrared spectroscopy is a non-invasive neuroimaging method which uses light to measure changes in cerebral blood oxygenation associated with brain activity. In this work, we demonstrate the ability to record and analyze images of brain activity in real-time using a 16-channel continuous wave optical NIRS system. We propose a novel real-time analysis framework using an adaptive Kalman filter and a state-space model based on a canonical general linear model of brain activity. We show that our adaptive model has the ability to estimate single-trial brain activity events as we apply this method to track and classify experimental data acquired during an alternating bilateral self-paced finger tapping task.

Optimization of probe geometry for diffuse optical brain imaging based on measurement density and distribution

Optode geometry plays an important role in achieving both good spatial resolution and spatial uniformity of detection in diffuse-optical-imaging-based brain activation studies. The quality of reconstructed images for six optode geometries were studied and compared using a laboratory tissue phantom model that contained an embedded object at two separate locations. The number of overlapping measurements per pixel (i.e., the measurement density) and their spatial distributions were quantified for all six geometries and were correlated with the quality of the resulting reconstructed images. The latter were expressed by the area ratio (AR) and contrast-to-noise ratio (CNR) between reconstructed and actual objects. Our results revealed clearly that AR and CNR depended on the measurement density asymptotically, having an optimal point for measurement density beyond which more overlapping measurements would not significantly improve the quality of reconstructed images. Optimization of probe geometry based on our method demonstrated that a practical compromise can be attained between DOI spatial resolution and measurement density.

A broadband continuous-wave multichannel near-infrared system for measuring regional cerebral blood flow and oxygen consumption in newborn piglets

Near-infrared spectroscopy (NIRS) is a promising technique for assessing brain function in newborns, particularly due to its portability and sensitivity to cerebral hemodynamics and oxygenation. Methods for measuring cerebral blood flow (CBF) and cerebral metabolic rate of oxygen (CMRO(2)) have been developed based on broadband continuous-wave NIRS. However, broadband NIRS apparatus typically have only one detection channel, which limits their applicability to measuring regional CBF and CMRO(2). In this study, a relatively simple multiplexing approach based on electronically controlled mechanical shutters is proposed to expand the detection capabilities from one to eight channels. The tradeoff is an increase in the sampling interval; however, this has negligible effects on CBF measurements for intervals less than or equal to 1 s. The ability of the system to detect focal brain injury was demonstrated in piglets by injecting endothelin-1 (ET-1) into the cerebral cortex. For validation, CBF was independently measured by computed tomography (CT) perfusion. The average reduction in CBF from the source-detector pair that interrogated the injured region was 51%+/-9%, which was in good agreement with the CBF reduction measured by CT perfusion (55%+/-5%). No significant changes in regional CMRO(2) were observed. The average regional differential pathlength prior to ET-1 injection was 8.4+/-0.2 cm (range of 7.1-9.6 cm) and did not significantly change after the injury.

Selection of regularization parameter for optical topography

The choice of the regularization parameter has a profound effect on the solution of ill-posed inverse problems such as optical topography. We review 11 different methods for selecting the Tikhonov regularization parameter that have been described previously in the literature. We test them on two trial problems, deblurring and optical topography, and conclude that the L-curve method is the method of choice, though in particularly ill-posed problems, generalized cross-validation may provide an alternative.

HomER: a review of time-series analysis methods for near-infrared spectroscopy of the brain

Near-infrared spectroscopy (NIRS) is a noninvasive neuroimaging tool for studying evoked hemodynamic changes within the brain. By this technique, changes in the optical absorption of light are recorded over time and are used to estimate the functionally evoked changes in cerebral oxyhemoglobin and deoxyhemoglobin concentrations that result from local cerebral vascular and oxygen metabolic effects during brain activity. Over the past three decades this technology has continued to grow, and today NIRS studies have found many niche applications in the fields of psychology, physiology, and cerebral pathology. The growing popularity of this technique is in part associated with a lower cost and increased portability of NIRS equipment when compared with other imaging modalities, such as functional magnetic resonance imaging and positron emission tomography. With this increasing number of applications, new techniques for the processing, analysis, and interpretation of NIRS data are continually being developed. We review some of the time-series and functional analysis techniques that are currently used in NIRS studies, we describe the practical implementation of various signal processing techniques for removing physiological, instrumental, and motion-artifact noise from optical data, and we discuss the unique aspects of NIRS analysis in comparison with other brain imaging modalities. These methods are described within the context of the MATLAB-based graphical user interface program, HomER, which we have developed and distributed to facilitate the processing of optical functional brain data.

Significance of background optical properties, time-resolved information and optode arrangement in diffuse optical imaging of term neonates

The significance of accurate knowledge of background optical properties and time-resolved information in reconstructing images of hemodynamic changes in the neonatal brain from diffuse optical imaging data was studied using Monte Carlo (MC) simulation. A segmented anatomical magnetic resonance (MR) image and literature-derived optical properties for each tissue type were used to create a voxel-based anatomical model. Small absorbing perturbations were introduced into the anatomical model to simulate localized hemodynamic responses related to brain activation. Perturbation MC (pMC) was used as the primary method of image reconstruction. For comparison, reconstructions were also performed using the finite element method (FEM) to solve the diffusion approximation (DA) to the radiative transfer equation (RTE). The effect of optode layout was investigated using three different grids. Of the factors studied, the density of the optode grid was found to have the greatest effect on image quality. The use of time-resolved information significantly improved the spatial accuracy with all optode grids. Adequate knowledge and modeling of the optical properties of the background was found to significantly improve the spatial accuracy of the reconstructed images and make the recovery of contrast of absorption changes more consistent over simplified modeling. Localization accuracy of small perturbations was found to be 2-3 mm with accurate a priori knowledge of the background optical properties, when a grid with high optode density (>1 optode cm(-2)) was used.

Wavelength dependence of sensitivity in spectral diffuse optical imaging: effect of normalization on image reconstruction

Near Infrared Diffuse Optical Tomography has the potential to be used as a non-invasive imaging tool for biological tissue specifically for the diagnosis and characterization of breast cancer. Most model based reconstruction algorithms rely on calculating and inverting a large Jacobian matrix. Although this method is flexible for a wide range of complex problems, it usually results in large image artifacts from hypersensitivity around the detectors. In this work a Jacobian normalization technique is presented which takes into account the varying magnitude of different optical parameters creating a more uniform update within a spectral image reconstruction model. Using simulated data the Jacobian normalization method is used to reconstructed images of absolute chromophore and scattering parameters which are qualitatively and quantitatively as compared to conventional methods. The hypersensitivity resulting in boundary artifacts are shown to be minimized with only a small additional computational cost.

Improving image quality of diffuse optical tomography with a projection-error-based adaptive regularization method

Diffuse optical tomography (DOT) reconstructs the images of internal optical parameter distribution using noninvasive boundary measurements. The image reconstruction procedure is known to be an ill-posed problem. In order to solve such a problem, a regularization technique is needed to constrain the solution space. In this study, a projection-error-based adaptive regularization (PAR) technique is proposed to improve the reconstructed image quality. Simulations are performed using a diffusion approximation model and the simulated results demonstrate that the PAR technique can improve reconstruction precision of object more effectively. The method is demonstrated to have low sensitivity to noise at various noise levels. Moreover, with the PAR method, the detectability of an object located both at the center and near the peripheral regions has been increased largely.

Auditory stimulus repetition effects on cortical hemoglobin oxygenation: a near-infrared spectroscopy investigation

The cortical response to repeated sensory stimuli plateaus (or declines) as repetition frequencies increase beyond 2-8 Hz. This study examined the underlying changes in cortical oxygenated and deoxygenated hemoglobin associated with this phenomenon using near-infrared spectroscopy. The optical signal was measured from 11 healthy volunteers listening to noise-burst trains presented at 2, 10, and 35 Hz. In a bilateral region consistent with the posterior superior temporal gyrus there was an inverse relationship between deoxyhemoglobin concentration change and stimulus frequency: greatest at 2 Hz, intermediate at 10 Hz, and smallest at 35 Hz. These findings provide preliminary support for a relationship between the perceptual characteristics of auditory stimuli and modulation of cortical oxygenation as measured via an emerging neuromonitoring technique.

Event-related functional near-infrared spectroscopy (fNIRS) based on craniocerebral correlations: reproducibility of activation?

The purpose of the present study was to assess the retest reliability of cortical activation detected by event-related functional near-infrared spectroscopy (fNIRS) based on craniocerebral correlations. Isolated functional activation was evoked in the motor cortex by a periodically performed finger-tapping task. During 44-channel fNIRS recording, 12 subjects performed 30 trials of right and left index finger tapping in two sessions. The retest interval was set to 3 weeks. Simple correlations of the contrast t-values supplemented by scatterplots, channel-wise intraclass correlation coefficients (ICC), as well as reproducibility indices for the size and the location of the detected activation were calculated. The results at the group level showed sufficient single measure ICCs (up to 0.80) and excellent reproducibility of the size and the location (up to 89% were reproducible). Comparisons of the intersession group amplitudes demonstrate that the fNIRS signals were stable across time in a retest study design: the number of significant differences was less than randomly occurring false-positive activated channels if an alpha level of 5% is chosen. Effect size analyses indicated that the intersession amplitude differences are small (mean < 0.25). For deoxyhemoglobin and oxyhemoglobin distinct statistical power profiles were revealed regarding the activation vs. baseline contrast as well as the intersession amplitude differences, indicating a higher sensitivity of deoxyhemoglobin for local hemodynamic changes. The results suggest that sensorimotor activation assessed by event-related fNIRS based on craniocerebral correlations is sufficiently reproducible at the group level.

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.

Retinotopic mapping of adult human visual cortex with high-density diffuse optical tomography

Functional neuroimaging is a vital element of neuroscience and cognitive research and, increasingly, is an important clinical tool. Diffuse optical imaging is an emerging, noninvasive technique with unique portability and hemodynamic contrast capabilities for mapping brain function in young subjects and subjects in enriched or clinical environments. We have developed a high-performance, high-density diffuse optical tomography (DOT) system that overcomes previous limitations and enables superior image quality. We show herein the utility of the DOT system by presenting functional hemodynamic maps of the adult human visual cortex. The functional brain images have a high contrast-to-noise ratio, allowing visualization of individual activations and highly repeatable mapping within and across subjects. With the improved spatial resolution and localization, we were able to image functional responses of 1.7 cm in extent and shifts of <1 cm. Cortical maps of angle and eccentricity in the visual field are consistent with retinotopic studies using functional MRI and positron-emission tomography. These results demonstrate that high-density DOT is a practical and powerful tool for mapping function in the human cortex.

Time-resolved optical imager for assessment of cerebral oxygenation

A time-resolved optical instrument allowing for noninvasive assessment of cerebral oxygenation is presented. The instrument is equipped with picosecond diode lasers, fast photodetectors, and time-correlated single photon counting electronics. This technology enables depth-resolved estimation of changes in absorption and, in consequence, assessment of changes in hemoglobin concentrations in the brain cortex. Changes in oxyhemoglobin (HbO(2)) and deoxyhemoglobin (Hb) can be evaluated selectively in extra- and intracerebral tissue compartments using the moments of distributions of times of flight of photons measured at two wavelengths in the near-infrared region. The combination of the data acquired from multiple sources and detectors located on the surface of the head with the depth-resolved analysis, based on the moments, enables imaging of cortex oxygenation. Results of the tests on physical phantoms as well as in vivo validation of the instrument during the motor stimulation experiment are presented.

Effect of probe arrangement on reproducibility of images by near-infrared topography evaluated by a virtual head phantom

The effect of the probe arrangement on the reproducibility of topographic images of the concentration changes in oxygenated hemoglobin and deoxygenated hemoglobin is evaluated by a virtual head phantom. A virtual head phantom consists of five types of tissue the 3D structure of which is based on a magnetic resonance imaging scan of an adult head. Localized and broadened brain activation is assumed in a virtual head phantom. The topographic images are obtained from the reflectance detected by the standard probe arrangement and the double-density probe arrangement. The uneven thickness of the superficial layer, which cannot be evaluated by the previous slab model, affects the distribution of measured activation in the topographic image, and this reduces the position reproducibility of near-infrared (NIR) topography with the standard probe arrangement. The overlapping measurements by the double-density probe arrangement can improve the reproducibility of the image obtained by NIR topography.

A spatial and temporal comparison of hemodynamic signals measured using optical and functional magnetic resonance imaging during activation in the human primary visual cortex

Functional near infrared spectro-imaging (fNIRSI) is potentially a very useful technique for obtaining information about the underlying physiology of the blood oxygenation level dependent (BOLD) signal used in functional magnetic resonance imaging (fMRI). In this paper the temporal and spatial statistical characteristics of fNIRSI data are compared to those of simultaneously acquired fMRI data in the human visual cortex during a variable-frequency reversing checkerboard activation paradigm. Changes in the size of activated volume caused by changes in checkerboard reversal frequency allowed a comparison of the behavior of the spatial responses measured by the two imaging methods. fNIRSI and fMRI data were each analyzed using standard correlation and fixed-effect group analyses of variance pathways. The statistical significance of fNIRSI data was found to be much lower than that of the fMRI data, due mainly to the low signal-to-noise of the measurements. Reconstructed images also showed that, while the time-course of changes in the oxy-, deoxy-, and total hemoglobin concentrations all exhibit high correlation with that of the BOLD response, the changes in the volume of tissue measured as "activated" by the BOLD response demonstrate a closer similarity to the corresponding changes in the oxy- and total hemoglobin concentrations than to that of the deoxyhemoglobin.

An integrated measurement system for simultaneous functional magnetic resonance imaging and diffuse optical tomography in human brain mapping

An integrated measurement system is described for performing simultaneous functional magnetic resonance imaging (fMRI) and diffuse optical tomography (DOT) for human brain mapping experiments. The components of this system consist of an MRI-compatible multi-overlapping-channel optical probe, methods for co-registration of optical and fMRI measurements, and DOT reconstruction algorithms with structural and physiological constraints derived from the MRI data. The optical probe is fully MRI-compatible in the sense that it produces negligible MR image distortion and does not require any modification to the MRI scanner or data acquisition protocol. The probe can be attached to any part of the head without posing any limitation on optical data acquisition. Co-registration of images from fMRI and optical measurements was achieved by localizing the positions of the optical fibers using MRI markers. Human studies show successful implementation of the entire system.

Improving performance of reflectance diffuse optical imaging using a multicentered mode

We propose a novel multicentered mode for arrangement of optical fibers to improve the imaging performance of reflectance diffuse optical imaging (rDOI). Simulations performed using a semi-infinite model show that the proposed multicentered geometries can achieve a maximum of 42 overlapping measurements. The contrast-to-noise ratio (CNR) analysis indicates that the best spatial resolution is 1 mm in radius and the contrast resolution is less than 1.05 for the multicentered geometries. The results from simulations indicate significant improvement in image quality compared to the single-centered mode and previous geometries. Additional experimental results on a single human subject lead to the conclusion that the proposed multicentered geometries are appropriate for exploring activations in the human brain. From the results of this research, we conclude that the proposed multicentered mode could advance the performance of rDOI both in image quality and practical convenience.

Diffuse optical tomography system to image brain activation with improved spatial resolution and validation with functional magnetic resonance imaging

Although most current diffuse optical brain imaging systems use only nearest- neighbor measurement geometry, the spatial resolution and quantitative accuracy of the imaging can be improved through the collection of overlapping sets of measurements. A continuous-wave diffuse optical imaging system that combines frequency encoding with time-division multiplexing to facilitate overlapping measurements of brain activation is described. Phantom measurements to confirm the expected improvement in spatial resolution and quantitative accuracy are presented. Experimental results showing the application of this instrument for imaging human brain activation are also presented. The observed improvement in spatial resolution is confirmed by functional magnetic resonance imaging.

Diffuse optical imaging of the whole head

Near-Infrared Spectroscopy (NIRS) and diffuse optical imaging (DOI) are increasingly used to detect hemodynamic changes in the cerebral cortex induced by brain activity. Until recently, the small number of optodes in NIRS instruments has hampered measurement of optical signals from diverse brain regions. Our new DOI system has 32 detectors and 32 sources; by arranging them in a specific pattern, we can cover most of the adult head. With the increased number of optodes, we can collect optical data from prefrontal, sensorimotor, and visual cortices in both hemispheres simultaneously. We describe the system and report system characterization measurements on phantoms as well as on human subjects at rest and during visual, motor, and cognitive stimulation. Taking advantage of the system's larger number of sources and detectors, we explored the spatiotemporal patterns of physiological signals during rest. These physiological signals, arising from cardiac, respiratory, and blood-pressure modulations, interfere with measurement of the hemodynamic response to brain stimulation. Whole-head optical measurements, in addition to providing maps of multiple brain regions' responses to brain activation, will enable better understandings of the physiological signals, ultimately leading to better signal processing algorithms to distinguish physiological signal clutter from brain activation signals.

Event-related functional near-infrared spectroscopy (fNIRS): Are the measurements reliable?

The purpose of the present study was to investigate the retest reliability of event-related functional near-infrared spectroscopy (fNIRS). Therefore, isolated functional activation was evoked in the occipital cortex by a periodic checkerboard stimulation. During a 52-channel fNIRS recording, 12 subjects underwent 60 trials of visual stimulation in two sessions. The retest interval was set to 3 weeks. Linear correlations of the contrast t values supplemented by scatter plots, channel-wise intraclass correlation coefficients (ICC) as well as reproducibility indices for the quantity of activated channels (RQUANTITY) and the location (ROVERLAP) of the detected activation were calculated. The results at the group level showed good reliability in terms of the single measure ICCs (up to 0.84) and excellent reproducibility quantified by RQUANTITY and ROVERLAP (up to 96% of the quantity and the location were reproducible), whereas the results at the single subjects' level were mediocre. Furthermore, the reliability assessed by single measurement ICCs improved if regarded at a cluster level.

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.

Evidence that cerebral blood volume can provide brain activation maps with better spatial resolution than deoxygenated hemoglobin

With the aim of evaluating the relative performance of hemodynamic contrasts for mapping brain activity, the spatio-temporal response of oxy-, deoxy-, and total-hemoglobin concentrations were imaged with diffuse optical tomography during electrical stimulation of the rat somatosensory cortex. For both 6-s and 30-s stimulus durations, total hemoglobin images provided smaller activation areas than oxy- or deoxy-hemoglobin images. In addition, analysis of regions of interest near the sagittal sinus vein show significantly greater contrast in both oxy- and deoxy-relative to total hemoglobin, suggesting that oximetric contrasts have larger draining vein contributions compared to total hemoglobin contrasts under the given stimulus conditions. These results indicate that total hemoglobin and cerebral blood volume may have advantages as hemodynamic mapping contrasts, particularly for large amplitude, longer duration stimulus paradigms.

Combining near-infrared tomography and magnetic resonance imaging to study in vivo breast tissue: implementation of a Laplacian-type regularization to incorporate magnetic resonance structure

An imaging system that simultaneously performs near infrared (NIR) tomography and magnetic resonance imaging (MRI) is used to study breast tissue phantoms and a healthy woman in vivo. An NIR image reconstruction that exploits the combined data set is presented that implements the MR structure as a soft-constraint in the NIR property estimation. The algorithm incorporates the MR spatially segmented regions into a regularization matrix that links locations with similar MR properties, and applies a Laplacian-type filter to minimize variation within each region. When prior knowledge of the structure of phantoms is used to guide NIR property estimation, root mean square (rms) image error decreases from 26 to 58%. For a representative in vivo case, images of hemoglobin concentration, oxygen saturation, water fraction, scattering power, and scattering amplitude are derived and the properties of adipose and fibroglandular breast tissue types, identified from MRI, are quantified. Fibroglandular tissue is observed to have more than four times as much water content as adipose tissue, almost twice as much blood volume, and slightly reduced oxygen saturation. This approach is expected to improve recovery of abnormalities within the breast, as the inclusion of structural information increases the accuracy of recovery of embedded heterogeneities, at least in phantom studies.

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.

Molecular imaging perspectives

Molecular imaging is an emerging technology at the life science/physical science interface which is set to revolutionize our understanding and treatment of disease. The tools of molecular imaging are the imaging modalities and their corresponding contrast agents. These facilitate interaction with a biological target at a molecular level in a number of ways. The diverse nature of molecular imaging requires knowledge from both the life and physical sciences for its successful development and implementation. The aim of this review is to introduce the subject of molecular imaging from both life science and physical science perspectives. However, we will restrict our coverage to the prominent in vivo molecular imaging modalities of magnetic resonance imaging, optical imaging and nuclear imaging. The physical basis of these imaging modalities, the use of contrast agents and the imaging parameters of sensitivity, temporal resolution and spatial resolution are described. Then, the specificity of contrast agents for targeting and sensing molecular events, and some applications of molecular imaging in biology and medicine are given. Finally, the diverse nature of molecular imaging and its reliance on interdisciplinary collaboration is discussed.

Simulation study of magnetic resonance imaging-guided cortically constrained diffuse optical tomography of human brain function

Diffuse optical imaging can measure brain activity noninvasively in humans through the scalp and skull by measuring the light intensity modulation arising from localized-activity-induced absorption changes within the cortex. Spatial resolution and localization accuracy are currently limited by measurement geometry to approximately 3 cm in the plane parallel to the scalp. Depth resolution is a more significant challenge owing to the limited angle tomography permitted by reflectance-only measurements. We combine previously established concepts for improving image quality and demonstrate, through simulation studies, their application for improving the image quality of adult human brain function. We show in a three-dimensional human head model that localization accuracy is significantly improved by the addition of measurements that provide overlapping samples of brain tissue. However, the reconstructed absorption contrast is significantly underestimated because its depth is underestimated. We show that the absorption contrast amplitude accuracy can be significantly improved by providing a cortical spatial constraint in the image reconstruction to obtain a better depth localization. The cortical constraint makes physiological sense since the brain-activity-induced absorption changes are occurring in the cortex and not in the scalp, skull, and cerebral spinal fluid. This spatial constraint is provided by segmentation of coregistered structural magnetic resonance imaging (MRI). However, the absorption contrast deep within the cortex is reconstructed superficially, resulting in an underestimation of the absorption contrast. The synthesis of techniques described here indicates that multimodality imaging of brain function with diffuse optical imaging and MRI has the potential to provide more quantitative estimates of the total and deoxyhemoglobin response to brain activation, which is currently not provided by either method independently. However, issues of depth resolution within the cortex remain to be resolved.