Three-dimensional optical tomography of hemodynamics in the human head

Dynamics of hemoglobin states in the sensorimotor cortex during motor tasks: a functional near infrared spectroscopy study

In this study, we used functional near infrared spectroscopy (fNIRS) imaging to investigate the independent levels of oxygenated, deoxygenated, and total hemoglobin (oxy-Hb, deoxy-Hb, and total-Hb, respectively) at the sensorimotor cortex during hand-grasping motor tasks. Our results showed that the activation of contralateral primary motor cortex (M1) exhibited increased oxy-Hb and reduced dexoy-Hb after hand grasping began. Meanwhile, the contralateral primary somatosensory cortex (S1) was deactivated with reductions of both oxy-Hb and deoxy-Hb concentration. The Hb circulation patterns indicated that the hand grasping demanded rapid and sufficient O2 supply at contralateral M1, which was achieved by the local vasodilation. The contralateral S1 presented decreased total-Hb via the mechanism of vasoconstriction, and maintained the local oxygenation level in a relatively stable state (mostly with O2 debt) to compensate the blood demanding at nearby M1. This study presented that fNIRS data can efficiently differentiate the activation of M1 from the deactivation of S1 during motor tasks, which can provide full interpretations of hemodynamic response to the neuronal activation in comparison with the Blood Oxygenation Level Dependent signal of functional magnetic resonance imaging.

Measuring cerebral hemodynamics with a modified magnetoencephalography system

Magnetoencephalography (MEG) systems are designed to noninvasively measure magnetic fields produced by neural electrical currents. This project examines the possibility of measuring hemodynamics with an MEG system that has been modified with dc electromagnets to measure magnetic susceptibility while maintaining the capability of measuring neural dynamics. A forward model is presented that simulates the interaction of an applied magnetic field with changes in magnetic susceptibility in the brain associated with hemodynamics. Model predictions are compared with an experiment where deionized water was pumped into an inverted flask under the MEG sensor array of superconducting quantum interference device (SQUID) gradiometers (R(2) = 0.98, p < 0.001). The forward model was used to simulate the SQUID readouts from hemodynamics in the scalp and brain induced by performing the Valsalva maneuver. Experimental human subject recordings (N = 10) were made from the prefrontal region during Valsalva using concurrent measurement with the modified MEG system and near-infrared spectroscopy (NIRS). The NIRS deoxyhemoglobin signal was found to correlate significantly with the SQUID readouts (R(2) = 0.84, p < 0.01). SQUID noise was found to increase with the applied field, which will need to be mitigated in future work. These results demonstrate the potential and technical challenges of measuring cerebral hemodynamics with a modified MEG system.

Singular value decomposition based regularization prior to spectral mixing improves crosstalk in dynamic imaging using spectral diffuse optical tomography

The spectrally constrained diffuse optical tomography (DOT) method relies on incorporating spectral prior information directly into the image reconstruction algorithm, thereby correlating the underlying optical properties across multiple wavelengths. Although this method has been shown to provide a solution that is stable, the use of conventional Tikhonov-type regularization techniques can lead to additional crosstalk between parameters, particularly in linear, single-step dynamic imaging applications. This is due mainly to the suboptimal regularization of the spectral Jacobian matrix, which smoothes not only the image-data space, but also the spectral mapping space. In this work a novel regularization technique based on the singular value decomposition (SVD) is presented that preserves the spectral prior information while regularizing the Jacobian matrix, leading to dramatically reduced crosstalk between the recovered parameters. Using simulated data, images of changes in oxygenated and deoxygenated hemoglobin concentrations are reconstructed via the SVD-based approach and compared with images reconstructed by using non-spectral and conventional spectral methods. In a 2D, two wavelength example, it is shown that the proposed approach provides a 98% reduction in crosstalk between recovered parameters as compared with conventional spectral reconstruction algorithms, and 60% as compared with non-spectrally constrained algorithms. Using a subject specific multilayered model of the human head, a noiseless dynamic simulation of cortical activation is performed to further demonstrate such improvement in crosstalk. However, with the addition of realistic noise in the data, both non-spectral and proposed algorithms perform similarly, indicating that the use of spectrally constrained reconstruction algorithms in dynamic DOT may be limited by the contrast of the signal as well as the noise characteristics of the system.

Optical imaging modalities for biomedical applications

Optical photographic imaging is a well known imaging method that has been successfully translated into biomedical applications such as microscopy and endoscopy. Although several advanced medical imaging modalities are used today to acquire anatomical, physiological, metabolic, and functional information from the human body, optical imaging modalities including optical coherence tomography, confocal microscopy, multiphoton microscopy, multispectral endoscopy, and diffuse reflectance imaging have recently emerged with significant potential for non-invasive, portable, and cost-effective imaging for biomedical applications spanning tissue, cellular, and molecular levels. This paper reviews methods for modeling the propagation of light photons in a biological medium, as well as optical imaging from organ to cellular levels using visible and near-infrared wavelengths for biomedical and clinical applications.

High-density diffuse optical tomography of term infant visual cortex in the nursery

Advancements in antenatal and neonatal medicine over the last few decades have led to significant improvement in the survival rates of sick newborn infants. However, this improvement in survival has not been matched by a reduction in neurodevelopmental morbidities with increasing recognition of the diverse cognitive and behavioral challenges that preterm infants face in childhood. Conventional neuroimaging modalities, such as cranial ultrasound and magnetic resonance imaging, provide an important definition of neuroanatomy with recognition of brain injury. However, they fail to define the functional integrity of the immature brain, particularly during this critical developmental period. Diffuse optical tomography methods have established success in imaging adult brain function; however, few studies exist to demonstrate their feasibility in the neonatal population. We demonstrate the feasibility of using recently developed high-density diffuse optical tomography (HD-DOT) to map functional activation of the visual cortex in healthy term-born infants. The functional images show high contrast-to-noise ratio obtained in seven neonates. These results illustrate the potential for HD-DOT and provide a foundation for investigations of brain function in more vulnerable newborns, such as preterm infants.

Image quality analysis of high-density diffuse optical tomography incorporating a subject-specific head model

High-density diffuse optical tomography (HD-DOT) methods have shown significant improvement in localization accuracy and image resolution compared to traditional topographic near infrared spectroscopy of the human brain. In this work we provide a comprehensive evaluation of image quality in visual cortex mapping via a simulation study with the use of an anatomical head model derived from MRI data of a human subject. A model of individual head anatomy provides the surface shape and internal structure that allow for the construction of a more realistic physical model for the forward problem, as well as the use of a structural constraint in the inverse problem. The HD-DOT model utilized here incorporates multiple source-detector separations with continuous-wave data with added noise based on experimental results. To evaluate image quality we quantify the localization error and localized volume at half maximum (LVHM) throughout a region of interest within the visual cortex and systematically analyze the use of whole-brain tissue spatial constraint within image reconstruction. Our results demonstrate that an image quality with less than 10 mm in localization error and 1000 m³ in LVHM can be obtained up to 13 mm below the scalp surface with a typical unconstrained reconstruction and up to 18 mm deep when a whole-brain spatial constraint based on the brain tissue is utilized.

Validating atlas-guided DOT: a comparison of diffuse optical tomography informed by atlas and subject-specific anatomies

We describe the validation of an anatomical brain atlas approach to the analysis of diffuse optical tomography (DOT). Using MRI data from 32 subjects, we compare the diffuse optical images of simulated cortical activation reconstructed using a registered atlas with those obtained using a subject's true anatomy. The error in localization of the simulated cortical activations when using a registered atlas is due to a combination of imperfect registration, anatomical differences between atlas and subject anatomies and the localization error associated with diffuse optical image reconstruction. When using a subject-specific MRI, any localization error is due to diffuse optical image reconstruction only. In this study we determine that using a registered anatomical brain atlas results in an average localization error of approximately 18 mm in Euclidean space. The corresponding error when the subject's own MRI is employed is 9.1 mm. In general, the cost of using atlas-guided DOT in place of subject-specific MRI-guided DOT is a doubling of the localization error. Our results show that despite this increase in error, reasonable anatomical localization is achievable even in cases where the subject-specific anatomy is unavailable.

Shedding light on words and sentences: near-infrared spectroscopy in language research

Investigating the neuronal network underlying language processing may contribute to a better understanding of how the brain masters this complex cognitive function with surprising ease and how language is acquired at a fast pace in infancy. Modern neuroimaging methods permit to visualize the evolvement and the function of the language network. The present paper focuses on a specific methodology, functional near-infrared spectroscopy (fNIRS), providing an overview over studies on auditory language processing and acquisition. The methodology detects oxygenation changes elicited by functional activation of the cerebral cortex. The main advantages for research on auditory language processing and its development during infancy are an undemanding application, the lack of instrumental noise, and its potential to simultaneously register electrophysiological responses. Also it constitutes an innovative approach for studying developmental issues in infants and children. The review will focus on studies on word and sentence processing including research in infants and adults.

A brief review on the history of human functional near-infrared spectroscopy (fNIRS) development and fields of application

This review is aimed at celebrating the upcoming 20th anniversary of the birth of human functional near-infrared spectroscopy (fNIRS). After the discovery in 1992 that the functional activation of the human cerebral cortex (due to oxygenation and hemodynamic changes) can be explored by NIRS, human functional brain mapping research has gained a new dimension. fNIRS or optical topography, or near-infrared imaging or diffuse optical imaging is used mainly to detect simultaneous changes in optical properties of the human cortex from multiple measurement sites and displays the results in the form of a map or image over a specific area. In order to place current fNIRS research in its proper context, this paper presents a brief historical overview of the events that have shaped the present status of fNIRS. In particular, technological progresses of fNIRS are highlighted (i.e., from single-site to multi-site functional cortical measurements (images)), introduction of the commercial multi-channel systems, recent commercial wireless instrumentation and more advanced prototypes.

Source localization approach for functional DOT using MUSIC and FDR control

In this paper, we formulate diffuse optical tomography (DOT) problems as a source localization problem and propose a MUltiple SIgnal Classification (MUSIC) algorithm for functional brain imaging application. By providing MUSIC spectra for major chromophores such as oxy-hemoglobin (HbO) and deoxy-hemoglobin (HbR), we are able to investigate the spatial distribution of brain activities. Moreover, the false discovery rate (FDR) algorithm can be applied to control the family-wise error in the MUSIC spectra. The minimum distance between the center of mass in DOT and the Montreal Neurological Institute (MNI) coordinates of target regions in experiments was between approximately 6 and 18 mm, and the displacement of the center of mass in DOT and fMRI ranged between 12 and 28 mm, which demonstrate the legitimacy of the DOT-based imaging. The proposed brain mapping method revealed its potential as an alternative algorithm to monitor the brain activation.

Three-dimensional optical topography of brain activity in infants watching videos of human movement

We present 3D optical topography images reconstructed from data obtained previously while infants observed videos of adults making natural movements of their eyes and hands. The optical topography probe was placed over the temporal cortex, which in adults is responsible for cognitive processing of similar stimuli. Increases in oxyhaemoglobin were measured and reconstructed using a multispectral imaging algorithm with spatially variant regularization to optimize depth discrimination. The 3D optical topography images suggest that similar brain regions are activated in infants and adults. Images were presented showing the distribution of activation in a plane parallel to the surface, as well as changes in activation with depth. The time-course of activation was followed in the pixel which demonstrated the largest change, showing that changes could be measured with high temporal resolution. These results suggest that infants a few months old have regions which are specialized for reacting to human activity, and that these subtle changes can be effectively analysed using 3D optical topography.

A quantitative spatial comparison of high-density diffuse optical tomography and fMRI cortical mapping

Functional neuroimaging commands a dominant role in current neuroscience research. However its use in bedside clinical and certain neuro-scientific studies has been limited because the current tools lack the combination of being non-invasive, non-ionizing and portable while maintaining moderate resolution and localization accuracy. Optical neuroimaging satisfies many of these requirements, but, until recent advances in high-density diffuse optical tomography (HD-DOT), has been hampered by limited resolution. While early results of HD-DOT have been promising, a quantitative voxel-wise comparison and validation of HD-DOT against the gold standard of functional magnetic resonance imaging (fMRI) has been lacking. Herein, we provide such an analysis within the visual cortex using matched visual stimulation protocols in a single group of subjects (n=5) during separate HD-DOT and fMRI scanning sessions. To attain the needed voxel-to-voxel co-registration between HD-DOT and fMRI image spaces, we implemented subject-specific head modeling that incorporated MRI anatomy, detailed segmentation, and alignment of source and detector positions. Comparisons of the visual responses found an average localization error between HD-DOT and fMRI of 4.4+/-1mm, significantly less than the average distance between cortical gyri. This specificity demonstrates that HD-DOT has sufficient image quality to be useful as a surrogate for fMRI.

Optomechanical imaging system for breast cancer detection

Imaging studies of the breast comprise three principal sensing domains: structural, mechanical, and functional. Combinations of these domains can yield either additive or wholly new information, depending on whether one domain interacts with the other. In this report, we describe a new approach to breast imaging based on the interaction between controlled applied mechanical force and tissue hemodynamics. Presented is a description of the system design, performance characteristics, and representative clinical findings for a second-generation dynamic near-infrared optical tomographic breast imager that examines both breasts simultaneously, under conditions of rest and controlled mechanical provocation. The expected capabilities and limitations of the developed system are described in relation to the various sensing domains for breast imaging.

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.

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.

Bedside optical imaging of occipital resting-state functional connectivity in neonates

Resting-state networks derived from temporal correlations of spontaneous hemodynamic fluctuations have been extensively used to elucidate the functional organization of the brain in adults and infants. We have previously developed functional connectivity diffuse optical tomography methods in adults, and we now apply these techniques to study functional connectivity in newborn infants at the bedside. We present functional connectivity maps in the occipital cortices obtained from healthy term-born infants and premature infants, including one infant with an occipital stroke. Our results suggest that functional connectivity diffuse optical tomography has potential as a valuable clinical tool for the early detection of functional deficits and for providing prognostic information on future development.

Altered prefrontal cortical function during processing of fear-relevant stimuli in pregnancy

In non-pregnant individuals, the prefrontal cortex (PFC) is involved in the regulation of emotion, and appears to play a role in anxiety. Near-infrared spectroscopy (NIRS) detects cortical neural activation without harmful radiation making it safe for use in pregnancy. The aims of this study were to assess neural circuitry involved in processing fear-relevant stimuli during pregnancy using NIRS, and to determine associations between activation of this circuitry, distress and anxiety symptoms, attention to threat, cortisol, estrogen, progesterone and testosterone levels. There was significant activation of the PFC in response to fearful faces compared to rest in both pregnant and control groups. Within pregnancy, the activation was most pronounced at trimester 2, compared to the other trimesters. In pregnant women only (all trimesters), PFC activation was significantly associated with increased distress and anxiety, but with decreased selective attention to masked fear. PFC activation was also significantly associated with increased levels of cortisol and testosterone in pregnancy. PFC function appears to be altered during processing of fear-relevant stimuli in pregnancy. Changes in hormone levels may lead to changes in PFC function, and in turn to changes in cognitive-affective processing and anxiety. Further work is needed, however, to explore precisely how PFC function is altered in pregnancy; it is possible that certain changes reflect altered processing of threat stimuli, while others reflect attempts to compensate for distressing and anxious symptoms that emerge during pregnancy.

The generation of tetrahedral mesh models for neuroanatomical MRI

In this article, we describe a detailed method for automatically generating tetrahedral meshes from 3D images having multiple region labels. An adaptively sized tetrahedral mesh modeling approach is described that is capable of producing meshes conforming precisely to the voxelized regions in the image. Efficient tetrahedral mesh improvement is then performed minimizing an energy function containing three terms: a smoothing term to remove the voxelization, a fidelity term to maintain continuity with the image data, and a novel elasticity term to prevent the tetrahedra from becoming flattened or inverted as the mesh deforms while allowing the voxelization to be removed entirely. The meshing algorithm is applied to structural MR image data that has been automatically segmented into 56 neuroanatomical sub-divisions as well as on two other examples. The resulting tetrahedral representation has several desirable properties such as tetrahedra with dihedral angles away from 0 and 180 degrees, smoothness, and a high resolution. Tetrahedral modeling via the approach described here has applications in modeling brain structure in normal as well as diseased brain in human and non-human data and facilitates examination of 3D object deformations resulting from neurological illness (e.g. Alzheimer's disease), development, and/or aging.

Multispectral quantitative photoacoustic imaging of osteoarthritis in finger joints

We present in vivo experimental evidence that multispectral quantitative photoacoustic tomography (qPAT) has the potential to detect osteoarthritis (OA) in the finger joints. In this pilot study, two OA patients and three healthy volunteers were enrolled, and their distal interphalangeal (DIP) joints were examined photoacoustically by a multispectral PAT scanner. Images of tissue physiological/functional parameters including oxyhemoglobin, deoxyhemoglobin, oxygen saturation, and water content, along with the tissue acoustic velocity of all the examined joints, were simultaneously recovered using a finite element reconstruction algorithm for multispectral photoacoustic measurements. The recovered multispectral photoacoustic images show that the OA joints have significantly elevated water content, decreased oxygen saturation, and increased acoustic velocity compared to the normal joints.

Brain specificity of diffuse optical imaging: improvements from superficial signal regression and tomography

Functional near infrared spectroscopy (fNIRS) is a portable monitor of cerebral hemodynamics with wide clinical potential. However, in fNIRS, the vascular signal from the brain is often obscured by vascular signals present in the scalp and skull. In this paper, we evaluate two methods for improving in vivo data from adult human subjects through the use of high-density diffuse optical tomography (DOT). First, we test whether we can extend superficial regression methods (which utilize the multiple source–detector pair separations) from sparse optode arrays to application with DOT imaging arrays. In order to accomplish this goal, we modify the method to remove physiological artifacts from deeper sampling channels using an average of shallow measurements. Second, DOT provides three-dimensional image reconstructions and should explicitly separate different tissue layers. We test whether DOT's depth-sectioning can completely remove superficial physiological artifacts. Herein, we assess improvements in signal quality and reproducibility due to these methods using a well-characterized visual paradigm and our high-density DOT system. Both approaches remove noise from the data, resulting in cleaner imaging and more consistent hemodynamic responses. Additionally, the two methods act synergistically, with greater improvements when the approaches are used together.

Quantitative two-dimensional photoacoustic tomography of osteoarthritis in the finger joints

We present for the first time in vivo experimental evidence that quantitative photoacoustic tomography (qPAT) has the potential to detect osteoarthritis (OA) in the finger joints. In this pilot study, 2 OA patients and 4 healthy volunteers were enrolled, and their distal interphalangeal (DIP) joints were examined photoacoustically by a PAT scanner. Tissue absorption coefficient images of all the examined joints, which were unavailable from the conventional PAT, were recovered using our finite element based qPAT approach. The recovered quantitative photoacoustic images revealed significant differences in the tissue absorption coefficient of the joint cavity (cartilage and synovial fluid) between the OA and normal joints.

Perturbation theory for the diffusion equation by use of the moments of the generalized temporal point-spread function. III. Frequency-domain and time-domain results

We study the performance of a previously proposed perturbation theory for the diffusion equation in frequency and time domains as they are known in the field of near infrared spectroscopy and diffuse optical tomography. We have derived approximate formulas for calculating higher order self- and mixed path length moments, up to the fourth order, which can be used in general diffusive media regardless of geometry and initial distribution of the optical properties, for studying the effect of absorbing defects. The method of Padé approximants is used to extend the validity of the theory to a wider range of absorption contrasts between defects and background. By using Monte Carlo simulations, we have tested these formulas in the semi-infinite and slab geometries for the cases of single and multiple absorbing defects having sizes of interest (d=4-10 mm, where d is the diameter of the defect). In frequency domain, the discrepancy between the two methods of calculation (Padé approximants and Monte Carlo simulations) was within 10% for absorption contrasts Deltamu(a)<or=0.2 mm(-1) for alternating current data, and usually to within 1 degrees for Deltamu(a)<or=0.1 mm(-1) for phase data. In time domain, the average discrepancy in the temporal range of interest (a few nanoseconds) was 2%-3% for Deltamu(a)<or=0.06 mm(-1). The proposed method is an effective fast forward problem solver: all the time-domain results presented in this work were obtained with a computational time of less than about 15 s with a Pentium IV 1.66 GHz personal computer.

Neonatal hemodynamic response to visual cortex activity: high-density near-infrared spectroscopy study

The neurodevelopmental outcome of neonatal intensive care unit (NICU) infants is a major clinical concern with many infants displaying neurobehavioral deficits in childhood. Functional neuroimaging may provide early recognition of neural deficits in high-risk infants. Near-infrared spectroscopy (NIRS) has the advantage of providing functional neuroimaging in infants at the bedside. However, limitations in traditional NIRS have included contamination from superficial vascular dynamics in the scalp. Furthermore, controversy exists over the nature of normal vascular, responses in infants. To address these issues, we extend the use of novel high-density NIRS arrays with multiple source-detector distances and a superficial signal regression technique to infants. Evaluations of healthy term-born infants within the first three days of life are performed without sedation using a visual stimulus. We find that the regression technique significantly improves brain activation signal quality. Furthermore, in six out of eight infants, both oxy- and total hemoglobin increases while deoxyhemoglobin decreases, suggesting that, at term, the neurovascular coupling in the visual cortex is similar to that found in healthy adults. These results demonstrate the feasibility of using high-density NIRS arrays in infants to improve signal quality through superficial signal regression, and provide a foundation for further development of high-density NIRS as a clinical tool.

Anatomical atlas-guided diffuse optical tomography of brain activation

We describe a neuroimaging protocol that utilizes an anatomical atlas of the human head to guide diffuse optical tomography of human brain activation. The protocol is demonstrated by imaging the hemodynamic response to median-nerve stimulation in three healthy subjects, and comparing the images obtained using a head atlas with the images obtained using the subject-specific head anatomy. The results indicate that using the head atlas anatomy it is possible to reconstruct the location of the brain activation to the expected gyrus of the brain, in agreement with the results obtained with the subject-specific head anatomy. The benefits of this novel method derive from eliminating the need for subject-specific head anatomy and thus obviating the need for a subject-specific MRI to improve the anatomical interpretation of diffuse optical tomography images of brain activation.

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.

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.

Resting-state functional connectivity in the human brain revealed with diffuse optical tomography

Mapping resting-state networks allows insight into the brain's functional architecture and physiology and has rapidly become important in contemporary neuroscience research. Diffuse optical tomography (DOT) is an emerging functional neuroimaging technique with the advantages, relative to functional magnetic resonance imaging (fMRI), of portability and the ability to simultaneously measure both oxy- and deoxyhemoglobin. Previous optical studies have evaluated the temporal features of spontaneous resting brain signals. Herein, we develop techniques for spatially mapping functional connectivity with DOT (fc-DOT). Simultaneous imaging over the motor and visual cortices yielded robust correlation maps reproducing the expected functional neural architecture. The localization of the maps was confirmed with task-response studies and with subject-matched fc-MRI. These fc-DOT methods provide a task-less approach to mapping brain function in populations that were previously difficult to research. Our advances may permit new studies of early childhood development and of unconscious patients. In addition, the comprehensive hemoglobin contrasts of fc-DOT enable innovative studies of the biophysical origin of the functional connectivity signal.

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.

Light transport in biological tissue using three-dimensional frequency-domain simplified spherical harmonics equations

The accuracy of the commonly used diffusion approximation as used in diffuse optical tomography is known to be limited in cases involving strong absorption and in these situations a higher ordered approximation is necessary. In this study, a light transport model has been developed based upon the three-dimensional frequency-domain simplified spherical harmonics (SP(N)) approximation for orders up to N = 7. The SP(N) data are tested against a semi-infinite multi-layered Monte Carlo model. It has been shown that the SP(N) approximation for higher orders (N >1) provides an increase in accuracy over the diffusion equation specifically near sources and at boundaries of regions with increased optical absorption. It is demonstrated that the error of fluence calculated near the sources between the diffusion approximation and the SP(N) model (N = 7) can be as large as 60%, therefore limiting the use of the diffusion approximation for small animal imaging and in situations where optical changes near sources are critical for tomographic reconstructions.

Dynamic functional and mechanical response of breast tissue to compression

Physiological tissue dynamics following breast compression offer new contrast mechanisms for evaluating breast health and disease with near infrared spectroscopy. We monitored the total hemoglobin concentration and hemoglobin oxygen saturation in 28 healthy female volunteers subject to repeated fractional mammographic compression. The compression induces a reduction in blood flow, in turn causing a reduction in hemoglobin oxygen saturation. At the same time, a two phase tissue viscoelastic relaxation results in a reduction and redistribution of pressure within the tissue and correspondingly modulates the tissue total hemoglobin concentration and oxygen saturation. We observed a strong correlation between the relaxing pressure and changes in the total hemoglobin concentration bearing evidence of the involvement of different vascular compartments. Consequently, we have developed a model that enables us to disentangle these effects and obtain robust estimates of the tissue oxygen consumption and blood flow. We obtain estimates of 1.9+/-1.3 micromol/100 mL/min for OC and 2.8+/-1.7 mL/100 mL/min for blood flow, consistent with other published values.

A femtosecond laser pacemaker for heart muscle cells

The intracellular effects of focused near-infrared femtosecond laser irradiation are shown to cause contraction in cultured neonatal rat cardiomyocytes. By periodic exposure to femtosecond laser pulse-trains, periodic contraction cycles in cardiomyocytes could be triggered, depleted, and synchronized with the laser periodicity. This was observed in isolated cells, and in small groups of cardiomyocytes with the laser acting as pacemaker for the entire group. A window for this effect was found to occur between 15 and 30 mW average power for an 80 fs, 82 MHz pulse train of 780 nm, using 8 ms exposures applied periodically at 1 to 2 Hz. At power levels below this power window, laser-induced cardiomyocyte contraction was not observed, while above this power window, cells typically responded by a high calcium elevation and contracted without subsequent relaxation. This laser-cell interaction allows the laser irradiation to act as a pacemaker, and can be used to trigger contraction in dormant cells as well as synchronize or destabilize contraction in spontaneously contracting cardiomyocytes. By increasing laser power above the window available for laser-cell synchronization, we also demonstrate the use of cardiomyocytes as optically-triggered actuators. To our knowledge, this is the first demonstration of remote optical control of cardiomyocytes without requiring exogenous photosensitive compounds.

A new planar left-handed metamaterial composed of metal-dielectric-metal structure

An improved planar structure of left-handed (LH) metamaterial is presented, and then designed and analyzed in microwave regime. In the anticipated LH frequency regime, the LH property is validated from the phenomena of backward wave propagation and negative refraction. To characterize the electromagnetic property of the planar metamaterial, we introduce the wedge method by constructing a wedge-shaped bulk LH metamaterial by stacking the planar LH metamaterials. The effective refractive index estimated by the wedge method is in excellent agreement with that retrieved by the inversion method from the transmission and reflection spectra.

UV-induced modification of stress distribution in optical fibers and its contribution to Bragg grating birefringence

This paper discusses the importance of stress-induced contributions to the photo-induced birefringence observed in fiber Bragg gratings. Optical tomography measurements are performed in exposed and unexposed fibers to extract the stress profiles induced by UV-writing of fiber Bragg gratings for various exposure levels. A photoelastic analysis and a high-order isoparametric finite elements method are then used to calculate the birefringence caused by stress profile modifications. The results are compared to the birefringence directly measured by spectral analysis of a chirped fiber grating with multiple phase-shifts. We can therefore estimate the fraction of the photo-induced birefringence due to stress-induced anisotropy following UV exposure.

Calibration of the second-order nonlinear optical susceptibility of surface and bulk of glass

A two-beam second-harmonic generation technique is developed to calibrate the magnitude of the second-order nonlinear optical susceptibility components of surface and bulk (multipolar origin) of isotropic materials. The values obtained for fused silica calibrated against ChiXXX of crystalline quartz are chi parallel parallel perpendicular = 7.9(4), chi perpendicular parallel parallel (+)gamma = 3.8(4), parallel perpendicular perpendicular perpendicular(+)gamma = 59(4), and delta' = 7.8(4) in units of 10(-22) m(2)/V. Similar values are obtained for BK7 glass. An alternative way of calibration against ChiXYZ of quartz is demonstrated. The technique could also be extended to characterize the susceptibility tensor of crystals as a convenient alternative to the Maker-fringe technique.

Polarization dependence of non-linear gain compression factor in semiconductor optical amplifier

We investigate the power and the polarization dependence of the intraband dynamics in a bulk semiconductor optical amplifier using both a 2.5-ps pump-probe experimental set-up in contra-propagation and a theoretical model. Our model is based on the rate equations and takes into account the polarization dependence of the gain. By comparing experimental and computational results we are able to highlight the dependences of the intraband dynamics and to extract the non-linear gain compression factor as a function of both pulse energy and polarization of the injected pulses.

Simple and robust image-based autofocusing for digital microscopy

A simple image-based autofocusing scheme for digital microscopy is demonstrated that uses as few as two intermediate images to bring the sample into focus. The algorithm is adapted to a commercial inverted microscope and used to automate brightfield and fluorescence imaging of histopathology tissue sections.

Silicon cross-connect filters using microring resonator coupled multimode-interference-based waveguide crossings

We report silicon cross-connect filters using microring resonator coupled multimode-interference (MMI) based waveguide crossings. Our experiments reveal that the MMI-based cross-connect filters impose lower crosstalk at the crossing than the conventional cross-connect filters using plain crossings, while offering a nearly symmetric resonance line shape in the drop-port transmission. As a proof-of-concept for cross-connection applications, we demonstrate on a silicon-on-insulator substrate (i) a 4-channel 1 x 4 linear-cascaded MMI-based cross-connect filter, and (ii) a 2-channel 2 x 2 array-cascaded MMI-based cross-connect filter.

Blazed high-efficiency x-ray diffraction via transmission through arrays of nanometer-scale mirrors

Diffraction gratings are ubiquitous wavelength dispersive elements for photons as well as for subatomic particles, atoms, and large molecules. They serve as enabling devices for spectroscopy, microscopy, and interferometry in numerous applications across the physical sciences. Transmission gratings are required in applications that demand high alignment and figure error tolerances, low weight and size, or a straight-through zero-order beam. However, photons or particles are often strongly absorbed upon transmission, e.g., in the increasingly important extreme ultraviolet (EUV) and soft x-ray band, leading to low diffraction efficiency. We demonstrate the performance of a critical-angle transmission (CAT) grating in the EUV and soft x-ray band that for the first time combines the advantages of transmission gratings with the superior broadband efficiency of blazed reflection gratings via reflection from nanofabricated periodic arrays of atomically smooth nanometer-thin silicon mirrors at angles below the critical angle for total external reflection. The efficiency of the CAT grating design is not limited to photons, but also opens the door to new, sensitive, and compact experiments and applications in atom and neutron optics, as well as for the efficient diffraction of electrons, ions, or molecules.

Birefringence and optical power confinement in horizontal multi-slot waveguides made of Si and SiO2

Through simulations and measurements, we show that in multi-slot thin film waveguides, the TM polarized modes can be confined mostly in the low refractive index layers of the waveguide. The structure consisted of alternating layers of a-Si and SiO(2), in the thickness range between 3 and 40 nm, for which the slots were the SiO(2) layers. Simulations were performed using the transfer matrix method and experiments using the m-line technique at 1.55 mum. The dependence of the birefringence and of the power confinement in the slots was studied as a function of the waveguide thickness, the Si and SiO(2) layer thicknesses, and the SiO(2) / Si layer thickness ratio. We find a large birefringence-a refractive index difference between TE and TM modes-as large as 0.8. For TM polarized modes, up to ~ 85% of the total power in the fundamental mode can be confined in the slots.

High-power, femtosecond, thermal-lens-shaped Yb:KGW oscillator

Thermal lens shaping for astigmatism compensation is extended to a high-power, diode-pumped, Yb:KGW laser by employing a gain crystal geometry designed for efficient polarized pumping. The 63MHz oscillator is soliton mode-locked with the aid of a saturable Bragg reflector to yield 250fs (347fs) pulses at an output power of 3.5W (5W). Frequency doubling of the 250fs pulses with an intrinsic efficiency >60% provides 1.65W of average green power.

Interaction of transverse modes in a single-frequency few-mode fiber amplifier caused by local gain saturation

We report on the behavior of modal polarization states in a single-frequency, ytterbium-doped, few-mode fiber amplifier. Experimental data show that the polarization of the individual transverse modes depends on the pump power and that the modes tend towards orthogonally polarized states with increasing gain. The observations can be explained by local gain saturation that favors the amplification of differently polarized modes.

Measurement of internal tissue optical properties at ultraviolet and visible wavelengths: Development and implementation of a fiberoptic-based system

A novel, multi-wavelength, fiberoptic system was constructed, evaluated and implemented to determine internal tissue optical properties at ultraviolet A (UVA) and visible (VIS) wavelengths. Inverse modeling was performed with a neural network to estimate absorption and reduced scattering coefficients based on spatially-resolved reflectance distributions. The model was calibrated with simulated reflectance datasets generated using a condensed Monte Carlo approach with absorption coefficients up to 85 cm(-1) and reduced scattering coefficients up to 118 cm(-1). After theoretical and experimental evaluations of the system, optical properties of porcine bladder, colon, esophagus, oral mucosa, and liver were measured at 325, 375, 405, 445 and 532 nm. These data provide evidence that as wavelengths decrease into the UVA, the dominant tissue chromophore shifts from hemoglobin to structural proteins such as collagen. This system provides a high level of accuracy over a wide range of optical properties, and should be particularly useful for in situ characterization of highly attenuating biological tissues in the UVA-VIS.

Low jitter Q-switched fiber laser using optically driven surface-normal saturable absorber modulator

A technique for stabilizing the repetition frequency of a passively Q-switched laser is presented using an optically driven surface-normal semiconductor modulator. A method is capable of significant reduction of the timing jitter in a passively Q-switched laser by optical triggering the saturable absorber semiconductor reflector. The experimental demonstration using passively Q-switched ytterbium-doped fiber laser shows the jitter reduction by factor of 1.66??10(3) from 50 mus down to 30 ns.

Fluorescence tomography characterization for sub-surface imaging with protoporphyrin IX

Optical imaging of fluorescent objects embedded in a tissue simulating medium was characterized using non-contact based approaches to fluorescence remittance imaging (FRI) and sub-surface fluorescence diffuse optical tomography (FDOT). Using Protoporphyrin IX as a fluorescent agent, experiments were performed on tissue phantoms comprised of typical in-vivo tumor to normal tissue contrast ratios, ranging from 3.5:1 up to 10:1. It was found that tomographic imaging was able to recover interior inclusions with high contrast relative to the background; however, simple planar fluorescence imaging provided a superior contrast to noise ratio. Overall, FRI performed optimally when the object was located on or close to the surface and, perhaps most importantly, FDOT was able to recover specific depth information about the location of embedded regions. The results indicate that an optimal system for localizing embedded fluorescent regions should combine fluorescence reflectance imaging for high sensitivity and sub-surface tomography for depth detection, thereby allowing more accurate localization in all three directions within the tissue.

Aberration reduction and unique light focusing in a photonic crystal negative refractive lens

Light focusing characteristics of a negative refractive lens fabricated out of a silicon-on-insulator photonic crystal (PC) slab are investigated theoretically and experimentally. It focuses in the near infrared, but the focal spot is degraded by a lens aberration. To reduce the aberration, we designed a composite PC that gives rise to a narrower focal spot. In addition, two unique functions of this lens are demonstrated: refocusing outside of the PC and parallel focusing, enabling image transfer and real image formation, respectively. These results prove the feasibility of an in-plane free space optical network based on negative refraction.

Measurement of layer-like hemodynamic trends in scalp and cortex: implications for physiological baseline suppression in functional near-infrared spectroscopy

A multidetector, continuous wave, near-infrared spectroscopy (NIRS) system is developed to examine whether the hemodynamics of the scalp and brain in adults contain significant layer-like hemodynamic trends. NIRS measurements are made using contrasting geometries, one with four detectors equidistant from a source 33 mm away, and one with detectors collinear with the source (5 to 33 mm away). When NIRS time series are acquired over the prefrontal cortex from resting adults using both geometries, variations among the time series are consistent with a substantially homogeneous two-layer model (p<0.001) and inconsistent with one dominated by heterogeneities. Additionally, when time series measured 5 mm from the source are subtracted from corresponding 33-mm signals via a least-squares algorithm, 60% of the hemoglobin changes are on average removed. These results suggest that hemodynamic trends present in the scalp can contribute significantly to NIRS measurements, and that attempts to reduce this noise by subtracting a simultaneous near-channel measurement using a two-layer model are justified. Such subtractions are then performed on NIRS measurements from two stimulus protocols. For systemic stimulations (Valsalva maneuver), the subtraction cancels the hemodynamic response, as desired. For localized stimulation of the occipital lobe (viewing a flickering pattern), the subtraction isolated a stimulus-correlated hemodynamic feature from background noise.

Visualizing localized dynamic changes during epileptic seizure onset in vivo with diffuse optical tomography

Diffuse optical tomography (DOT) is a promising functional imaging modality due to its ability to provide quantitative and dynamic tomographic imaging of brain functions. This pilot study was conducted to demonstrate that DOT can be used to visualize the changes in local hemodynamics during seizures. The focal seizure was induced by microinjection of 10 microl of 1.9 mM GABAA antagonist bicuculline methiodide (BMI) into the left parietal neocortex of male Harlen Sprague-Dawley rats, which was imaged by a multispectral continuous-wave DOT system. Functional images were obtained by our finite element based image reconstruction algorithm. A series of dynamic 2D images were obtained to delineate the time course of concentration changes of oxyhaemoglobin, deoxyhaemoglobin, and total hemoglobin in the rat brain during seizure onset. The BMI induced epileptic foci were localized and observed over time from the images obtained. Our results suggest that diffuse optical tomography may be a promising modality for epilepsy imaging due to its ability to localize epileptic foci as well as its potential to map the functional activity in the area of human cerebral cortex in planning of epilepsy surgery.