- Open Access
Kalman estimator- and general linear model-based on-line brain activation mapping by near-infrared spectroscopy
© Hu et al; licensee BioMed Central Ltd. 2010
- Received: 5 July 2010
- Accepted: 8 December 2010
- Published: 8 December 2010
Near-infrared spectroscopy (NIRS) is a non-invasive neuroimaging technique that recently has been developed to measure the changes of cerebral blood oxygenation associated with brain activities. To date, for functional brain mapping applications, there is no standard on-line method for analysing NIRS data.
In this paper, a novel on-line NIRS data analysis framework taking advantages of both the general linear model (GLM) and the Kalman estimator is devised. The Kalman estimator is used to update the GLM coefficients recursively, and one critical coefficient regarding brain activities is then passed to a t-statistical test. The t-statistical test result is used to update a topographic brain activation map. Meanwhile, a set of high-pass filters is plugged into the GLM to prevent very low-frequency noises, and an autoregressive (AR) model is used to prevent the temporal correlation caused by physiological noises in NIRS time series. A set of data recorded in finger tapping experiments is studied using the proposed framework.
The obtained results suggest that the method can effectively track the task related brain activation areas, and prevent the noise distortion in the estimation while the experiment is running. Thereby, the potential of the proposed method for real-time NIRS-based brain imaging was demonstrated.
This paper presents a novel on-line approach for analysing NIRS data for functional brain mapping applications. This approach demonstrates the potential of a real-time-updating topographic brain activation map.
- Brain Activation
- Temporal Correlation
- Hemodynamic Response
- NIRS Data
- Kalman Estimator
Near-infrared spectroscopy (NIRS), an emerging brain imaging technique, measures the hemodynamic changes that effectively reflect the brain activity occurring while people perform a wide range of mental tasks [1–5]. It can provide both topographic [2, 4, 6] and tomographic [1, 7] brain images. Specifically, NIRS monitors the regional cerebral blood flow (rCBF) variation by measuring, through the skull, the absorption changes of near-infrared light at wavelengths between 650 nm and 950 nm . These changes are caused by the concentration variations of oxy-hemoglobin (HbO) and deoxy-hemoglobin (HbR), two primary absorbing chromophores in brain capillary blood.
NIRS, compared with other prevalent brain imaging and activity measurement techniques such as electroencephalography (EEG) and functional magnetic resonance imaging (fMRI), offers itself as a trade-off between spatial and temporal resolutions. The usability and drawbacks of NIRS methods, as discussed in a detailed review and comparison with other neuroimaging methods, was provided by Perrey . fMRI has been used over the past decade in a growing number of applications. The critical drawbacks of the fMRI-based approaches, however, are the cost and the non-portability of the fMRI scanner. In fact, another comprehensive review , in comparing the respective features of NIRS and fMRI, concluded that NIRS has great potentials for neurological and psychiatric applications, due to its simplicity, portability, and insensitivity to motion artifacts. Meanwhile, the EEG technique is limited, due to its poor spatial resolution and low signal-to-noise ratios in many applications; NIRS can provide comparatively better quality in these aspects . Indeed NIRS, in its wide applicability, might help to bring functional imaging to the patient's bedside .
There is currently no standard method of topographic NIRS data analysis for brain mapping. In NIRS detection of hemodynamic responses, the light attenuation measured by the equipment needs to be converted to HbO and HbR concentration changes via the modified Beer-Lambert law (MBLL) . Hence, a differential path length factor (DPF, intra- and inter-subject varying) in the MBLL has to be assumed to account for the increase of the path length between a source and a detector .
The classical approach in topographic NIRS data analysis is a paired t-test to determine if a concentration change between two states (for instance, "rest" vs. "task") is statistically significant. Many researchers nowadays use this approach [12–14], because it is simple and, thus, can provide a quick assessment to the task. One of the most popular tools in this regard is a Matlab-based program known as HomER ( http://www.nmr.mgh.harvard.edu/PMI/).
However, there are limitations to the classical t-test. First, a maximum activation period needs to be predefined. The remaining temporal information not included in the defined activation period therefore is ignored, leading to underestimation of brain activation. Another problem is the DPF assumption. Since the DPF is intra- and inter-subject variant and impossible to be measured for every measurement location with commonly obtainable continuous wave NIRS equipment, use of a constant DPF leads to biased estimation of concentration changes .
To overcome such problems of the classical t-test, a number of research groups have used various general linear model (GLM)-based methods for analysis of NIRS data [16–20]. The GLM-based methods were initially developed for fMRI-based functional brain mapping . The GLM is a statistical linear concept that explains measured data in the form of a linear combination of several explanatory variables plus an error term. The explanatory variables, modelled according to the time course, separately account for the brain-activity-evoked signals and noises. Therefore, the estimation of brain activity is reduced to the problem of estimating the relevant coefficients with proper statistics.
The GLM-based methods negate the need for user-defined rest and task periods, because the response is modelled according to the entire time course. The temporal information over the entire time course, thus, is examined. On the other hand, the GLM-based methods investigate the temporal variation pattern of the signal, and estimate the coefficients with statistics at different measurement locations, separately. Therefore, these methods are robust in cases where an assumed constant DPF is used.
GLM-based methods, however, cannot offer on-line analysis. This constrains their use in applications where real-time information or feedback is required. Real-time brain imaging data analysis in comparison with off-line methods, significantly, improves the information acquisition rate and the feedback speed. Furthermore, it may potentially benefit the development of brain-computer interfaces (BCI).
The critical task in modifying a GLM-based method on-line is the recursive estimation of GLM coefficients. Previous researchers have used different methods to achieve recursive estimation, including recursive least square , Cholesky-decomposition-based recursive least square , and Kalman filtering [24, 25]. All these approaches can effectively and recursively estimate GLM coefficients.
It is not sufficient to draw an updated brain activation map only by estimating GLM coefficients on-line. It is known that NIRS time series contain noises from different sources. Very-low-frequency noises caused by optodes shifts or slow cardiac/vascular artifacts , for example, might lead to biased estimation. The temporal correlation caused by physiological (cardiac, respiratory, blood pressure) noises might lead to an inflated t-value and, thereby, overestimation of brain activation. Furthermore, for functional brain mapping applications, a statistical test is very important, since it will provide significance verification of the derived estimation. All these issues with regard to on-line versions of GLM-based methods need to be addressed.
On-line estimation of GLM coefficients is generally discussed in . In , the whole framework provided on-line estimation of GLM coefficients, and a relevant statistical test analysed the fMRI data. However, this work did not consider temporal correlation in fMRI data, which might also exist in NIRS data. In , a framework for NIRS-based BCI applications was developed that can estimate GLM coefficients without statistical information in classifying different hand tasks in real-time. In , the feasibility of estimating GLM coefficients on-line using a Kalman filter by studying a set of fMRI data was examined.
In the present study, we develop an on-line Kalman-estimator- and GLM-based NIRS data processing framework for brain activation mapping. We aim to answer two questions. (i) Is it possible to covert the off-line GLM-based method to an on-line version for NIRS-based brain activation mapping? The proposed method can provide updated brain activation maps on-line as well as track task-related brain-active areas from an early stage while the experiment is running. (ii) Can the proposed method prevent noises that distort the estimation while data is sequentially incorporated?
NIRS measurement system and experimental procedure
In the experiment, the subjects were asked to perform a finger-tapping task. The experiment consisted of a 42 sec preparation period and 10 sessions. Each session included a 21 sec finger-tapping period and a 30 sec rest period. Accordingly, the total duration of the experiment was 552 sec.
Analysis framework of NIRS data
NIRS measurement model
where , λ is the wavelength of the laser source, i indicates channel number, α HbO [μM-1mm-1] and α HbR [μM-1mm-1] are the extinction coefficients of the HbO and HbR, L is the distance between the source and the detector, and DPF is the differential path length factor. In the present study, the optical density variation was derived by dividing the light intensity measured at each time step by the light intensity measured at the first time step.
General linear model (GLM)
In NIRS-based studies, both HbO and HbR concentration changes can reflect changes in the rCBF. However, it has been suggested that HbO is a more sensitive indicator of such changes . Therefore, only the HbO concentration change data was considered in the present study.
The GLM design process is described in detail in ; we provide only a brief description here. A design matrix H including a set of explanatory variables is predefined in order to model the observed NIRS time series. Five explanatory variables are considered. The first variable models the HbO concentration changes (the brain activity signals) using a stimulus vector convolved with the basis function (BF, a double-gamma model; ). The second one models the baseline level, and the remaining variables represent a set of high-pass filters (discrete cosine transform, DCT)  with a cut-off frequency of 0.0006 Hz.
where β i (k) is the coefficient vector quantifying the magnitude of the explanatory variables. In vector β i (k), we are interested in the component β i 1(k), which reflects the magnitude of the task-evoked brain response: By statistically determining if it is greater than zero, the existence of brain activity at the area covered by channel i can be confirmed (it is greater than zero) or ruled out (it is less than zero).
It is worth noting that we make an assumption, |β(k) - β(k- 1)| < ζ, for the AR(1) model used on-line in the current study, where ζ > 0 is an arbitrary small number. This assumption compromises the model's robustness. However, as the result shows, the temporal correlation can be effectively reduced. In the present study, β i (k) was updated with a Kalman estimator.
where E(k) = H(k)P-(k)H T +R, , K(k) is the Kalman gain, and P(k) is the updated error covariance matrix. In this notation, the superscript (-) refers to the intermediate state and covariance predictions provided by the state update model, which are then modified by the measured data to produce the next state value. In the present study, the state vector was initialized to zero. The a priori estimates of the process and observation noise covariances (Q and R respectively) were (1%/sec)2 and (0.5 μM/sec)2, according to a restricted maximum likelihood (ReML) estimation and an empirical-experimental performance check based on a set of training data. We collected the training data during 3 sessions of finger tapping for each subject. We estimated the Q and R values in two steps: They were estimated separately from each of the subjects by ReML, averaged, and then adjusted according to the performance in practically estimation based on the training data.
where L is the number of regressors. Therefore, the null hypothesis c T β i (k) = 0 was assessed by comparing t i (k) with a t-distribution with k-L degrees of freedom. By setting proper p-values with Bonferroni correction, a statistical activation map of the detected area could be displayed.
To simulate a real-time process, each measured data was incorporated sequentially in the analysis and updated at each time step. The entire procedure was simulated with Matlab at this stage. The computation time for one processing step was approximately 0.015 ± 0.0025 sec (mean ± standard deviation, 4 subjects averaged).
Time and location of the brain activities tracked
(p in < 0.05)
(p b , with Bonferroni correction)
Changes in cerebral oxygenation reflect cerebral functional activity. In the present work, an updated version of the GLM was devised and used for on-line brain activation mapping. To demonstrate this framework, a finger-tapping task activation mapping study was carried out. This study allowed us to highlight several important features of our framework.
An obvious advantage of our method is its real-time applicability. The Kalman estimator has an acceptable computational overload, which allowed us to implement this method using Matlab in real-time (1.81 Hz in this study). The proposed method displays an updated brain activation map while the experiment is running. The finger-tapping-related brain-active area can be identified and tracked at the early stage of an experiment. The data analysis stage in the classical GLM framework can be conducted only after an experiment, and thus needs extra time. By contrast, with the proposed method, the model coefficients are recursively estimated, and the brain activation map can be updated at each time step. Accordingly, the method is able to track the brain-active area in an ongoing experiment, and provide an early warning to the experimenter when the subject is not responding appropriately or the system is not working properly or well. Thus, both the subject/patient and researcher can receive feedback in real-time.
There are some limitations of our method. First, the design matrix in the GLM was designed before the experiment. Different subjects might have individual hemodynamic responses. However, we assumed a general hemodynamic response model for different subjects in our framework. Our method might not work properly if a notable difference exists between the hemodynamic response of a subject and our modelled one. In , a set of linear combined Gaussian-based temporal functions was used to model the hemodynamic response. This model can approximate complex hemodynamic responses among different subjects, by means of multiple overlapping Gaussian functions. However, a recent work  suggested that event-evoked hemodynamic responses were similar among different subjects. In the current study, we also found a similarity in the hemodynamic responses, among all of our subjects. Therefore, we consider that it was sufficient to use the general hemodynamic model for estimation in current study.
The second limitation lies in the fact that we made only a preliminary conclusion about the feasibility of the proposed on-line framework, based on the results from 5 subjects. In the current study, brain activity at the primary motor cortex was found in 4 out of the 5 subjects investigated using the proposed framework. There was no brain activation area identified for subject 5, because the measurement noises (electric voltage fluctuations) induced by the equipment were very large in the case of this subject. Thus, no brain-active area could be identified based on the measured data. In order to produce a complete picture of the feasibility of the proposed framework, data needs to be measured for a greater number of subjects.
The t-statistics from the current study are based on the estimation of mean square error (MSE) via Equation (14). Therefore, while the experiment is running, the cumulative type I error increases as the test number grows. Thus, a post-hoc analysis is also needed to guarantee the statistical significance. We used Bonferroni correction for post-hoc analysis at each step in the current study. The p in for an individual test was set to 0.05, and the p b after Bonferroni correction was calculated as p b = (1 - (1 - p in ) j ]/j, where j represents the number of tests. As indicated in Table 1 column 4, the post-hoc analysis reduced the performance of the proposed framework, though it still yields promising results. Nevertheless, it might not be suitable for long-duration experiments (e.g. 6 hours' measurement for sleep studies). When the number of tests increases to a large value, the p b might be difficult to satisfy, resulting in an underestimation of brain activation. This is one of the limitations of the proposed framework. The problem probably can be avoided, either by reducing the sampling frequency over a long-duration measurement, thereby reducing the number of tests, or by controlling the experimentation time.
The time needed for locating different brain-active areas differed in the current study. Because Kalman filtering is a recursive process, the new information will be added as it arrives at each sampling time. Therefore, the estimation results from early stages are less reliable than those from later ones. In another aspect, the brain is a heterogeneous and dynamic organ that can process parallel work (in different brain areas), and so the same estimation process might perform differently for different brain areas. In our study, channel 4 of subject 4 was identified as a brain-active area during the finger-tapping task. The area covered by channel 4 is Brodmann area 8, responsible for voluntary eye movement, which is not expected to be active during finger tapping. This particular result might be accounted for simply as an indication that subject 4 was moving his eyes frequently while tapping his finger.
The DYNOT equipment in the present study employs a brush-less DC servomotor to provide for light source beam positioning . The motor is driven by a controller containing freely programmable microprocessor. In this experiment, only a simple open-loop control scheme was designed for controlling the motor moving. Some complex control algorithms can be applied to achieve a more precise motor controlling [41–46].
One further issue in NIRS topographical applications is the fact that only a small portion of a detected area is measurable with sparsely distributed channels. Either tomographical NIRS detection or application of interpolation techniques might help increase it . Both solutions would require extra processing time. Therefore, the issue of the development of a proper means of achieving higher spatial resolution in real-time applications also will have to be resolved in future research.
A new Kalman-estimator- and GLM-based NIRS data analysis framework was demonstrated for real-time imaging of brain activity. In an experiment, this framework allowed for updating of the topographic activation map of the detected brain area, and, additionally, it could locate the activated brain area at an early stage by analyzing the noise-containing raw data.
This research was supported by the World Class University program through the National Research Foundation of Korea funded by the Ministry of Education, Science and Technology, Korea (grant no. R31-2008-000-20004-0).
- Boas DA, Dale AM, Franceschini MA: Diffuse optical imaging of brain activation: approaches to optimizing image sensitivity, resolution, and accuracy. Neuroimage 2004, 23: S275-S288. 10.1016/j.neuroimage.2004.07.011View ArticleGoogle Scholar
- Koizumi H, Yamamoto T, Maki A, Yamashita Y, Sato H, Kawaguchi H, Ichikawa N: Optical topography: practical problems and new applications. Appl Opt 2003, 42: 3054–3062. 10.1364/AO.42.003054View ArticleGoogle Scholar
- Obrig H, Villringer A: Beyond the visible--imaging the human brain with light. J Cereb Blood Flow Metab 2003, 23: 1–18. 10.1097/00004647-200301000-00001View ArticleGoogle Scholar
- Wolf M, Ferrari M, Quaresima V: Progress of near-infrared spectroscopy and topography for brain and muscle clinical applications. J Biomed Opt 2007, 12: 062104. 10.1117/1.2804899View ArticleGoogle Scholar
- Perrey S: Non-invasive NIR spectroscopy of human brain function during exercise. Methods 2008, 45: 289–299. 10.1016/j.ymeth.2008.04.005View ArticleGoogle Scholar
- Toronov VY, Zhang X, Webb AG: A spatial and temporal comparison of hemodynamic signals measured using optical and functional magnetic resonance imaging during activation in the human primary visual cortex. Neuroimage 2007, 34: 1136–1148. 10.1016/j.neuroimage.2006.08.048View ArticleGoogle Scholar
- Barbour RL, Graber HL, Pei YL, Zhong S, Schmitz CH: Optical tomographic imaging of dynamic features of dense-scattering media. Journal of the Optical Society of America a-Optics Image Science and Vision 2001, 18: 3018–3036. 10.1364/JOSAA.18.003018View ArticleGoogle Scholar
- Irani F, Platek SM, Bunce S, Ruocco AC, Chute D: Functional near infrared spectroscopy (fNIRS): An emerging neuroimaging technology with important applications for the study of brain disorders. Clinical Neuropsychologist 2007, 21: 9–37. 10.1080/13854040600910018View ArticleGoogle Scholar
- Gratton G, Fabiani M: The event-related optical signal: a new tool for studying brain function. International Journal of Psychophysiology 2001, 42: 109–121. 10.1016/S0167-8760(01)00161-1View ArticleGoogle Scholar
- Cope M, Delpy DT, Reynolds EO, Wray S, Wyatt J, van der Zee P: Methods of quantitating cerebral near infrared spectroscopy data. Adv Exp Med Biol 1988, 222: 183–189.View ArticleGoogle Scholar
- Zhao H, Tanikawa Y, Gao F, Onodera Y, Sassaroli A, Tanaka K, Yamada Y: Maps of optical differential pathlength factor of human adult forehead, somatosensory motor and occipital regions at multi-wavelengths in NIR. Phys Med Biol 2002, 47: 2075–2093. 10.1088/0031-9155/47/12/306View ArticleGoogle Scholar
- Hashimoto T, Minagawa-Kawai Y, Kojima S: Motion illusion activates the visual motion area of the brain: A near-infrared spectroscopy (NIRS) study. Brain Research 2006, 1077: 116–122. 10.1016/j.brainres.2005.11.050View ArticleGoogle Scholar
- Shibuya K, Kuboyama N: Human motor cortex oxygenation during exhaustive pinching task. Brain Research 2007, 1156: 120–124. 10.1016/j.brainres.2007.05.009View ArticleGoogle Scholar
- Watanabe H, Homae F, Nakano T, Taga G: Functional activation in diverse regions of the developing brain of human infants. Neuroimage 2008, 43: 346–357. 10.1016/j.neuroimage.2008.07.014View ArticleGoogle Scholar
- Huppert TJ, Diamond SG, Franceschini MA, Boas DA: HomER: a review of time-series analysis methods for near-infrared spectroscopy of the brain. Applied Optics 2009, 48: D280-D298. 10.1364/AO.48.00D280View ArticleGoogle Scholar
- Schroeter ML, Bucheler MM, Muller K, Uludag K, Obrig H, Lohmann G, Tittgemeyer M, Villringer A, von Cramon DY: Towards a standard analysis for functional near-infrared imaging. Neuroimage 2004, 21: 283–290. 10.1016/j.neuroimage.2003.09.054View ArticleGoogle Scholar
- Koh PH, Glaser DE, Flandin G, Kiebel S, Butterworth B, Maki A, Delpy DT, Elwell CE: Functional optical signal analysis: a software tool for near-infrared spectroscopy data processing incorporating statistical parametric mapping. J Biomed Opt 2007, 12: 064010. 10.1117/1.2804092View ArticleGoogle Scholar
- Lu CM, Zhang YJ, Biswal BB, Zang YF, Peng DL, Zhu CZ: Use of fNIRS to assess resting state functional connectivity. J Neurosci Methods 186: 242–249. 10.1016/j.jneumeth.2009.11.010Google Scholar
- Plichta MM, Heinzel S, Ehlis AC, Pauli P, Fallgatter AJ: Model-based analysis of rapid event-related functional near-infrared spectroscopy (NIRS) data: a parametric validation study. Neuroimage 2007, 35: 625–634. 10.1016/j.neuroimage.2006.11.028View ArticleGoogle Scholar
- Ye JC, Tak S, Jang KE, Jung J, Jang J: NIRS-SPM: statistical parametric mapping for near-infrared spectroscopy. Neuroimage 2009, 44: 428–447. 10.1016/j.neuroimage.2008.08.036View ArticleGoogle Scholar
- Worsley KJ, Friston KJ: Analysis of fMRI time-series revisited--again. Neuroimage 1995, 2: 173–181. 10.1006/nimg.1995.1023View ArticleGoogle Scholar
- McGilchrist CA, Sandland RL: Recursive Estimation of the General Linear Model with Dependent Errors. Journal of the Royal Statistical Society 1979, 41: 65–68.MathSciNetGoogle Scholar
- Bagarinao E, Matsuo K, Nakai T, Sato S: Estimation of general linear model coefficients for real-time application. Neuroimage 2003, 19: 422–429. 10.1016/S1053-8119(03)00081-8View ArticleGoogle Scholar
- Abdelnour AF, Huppert T: Real-time imaging of human brain function by near-infrared spectroscopy using an adaptive general linear model. Neuroimage 2009, 46: 133–143. 10.1016/j.neuroimage.2009.01.033View ArticleGoogle Scholar
- Buchel C, Friston KJ: Dynamic changes in effective connectivity characterized by variable parameter regression and Kalman filtering. Hum Brain Mapp 1998, 6: 403–408. 10.1002/(SICI)1097-0193(1998)6:5/6<403::AID-HBM14>3.0.CO;2-9View ArticleGoogle Scholar
- Obrig H, Neufang M, Wenzel R, Kohl M, Steinbrink J, Einhaupl K, Villringer A: Spontaneous low frequency oscillations of cerebral hemodynamics and metabolism in human adults. Neuroimage 2000, 12: 623–639. 10.1006/nimg.2000.0657View ArticleGoogle Scholar
- Matthews F, Pearlmutter BA, Ward TE, Soraghan C, Markham C: Hemodynamics for braincomputer interfaces. IEEE Signal Processing Magazine 2008, 25: 87–94. 10.1109/MSP.2008.4408445View ArticleGoogle Scholar
- Hoshi Y, Kobayashi N, Tamura M: Interpretation of near-infrared spectroscopy signals: a study with a newly developed perfused rat brain model. J Appl Physiol 2001, 90: 1657–1662.Google Scholar
- Friston KJ: Statistical parametric mapping: the analysis of funtional brain images. 1st edition. Amsterdam; Boston: Elsevier/Academic Press; 2007.Google Scholar
- Glover GH: Deconvolution of impulse response in event-related BOLD fMRI. Neuroimage 1999, 9: 416–429. 10.1006/nimg.1998.0419View ArticleGoogle Scholar
- Kalman RE: A new approach to linear filtering and prediction problems. Journal of Basic Engineering 1960, 82: 35–45.View ArticleGoogle Scholar
- Kim YS, Hong KS: An IMM algorithm with federated information mode-matched filters for AGV. International Journal of Adaptive Control and Signal Processing 2007, 21: 533–555. 10.1002/acs.939View ArticleGoogle Scholar
- Kim YS, Hong KS: Federated information mode-matched filters in ACC environment. International Journal of Control Automation and Systems 2005, 3: 173–182.Google Scholar
- Turnip A, Hong KS, Park S: Modeling of a hydraulic engine mount for active pneumatic engine vibration control using the extended Kalman filter. Journal of Mechanical Science and Technology 2009, 23: 229–236. 10.1007/s12206-008-1105-2View ArticleGoogle Scholar
- Diamond SG, Huppert TJ, Kolehmainen V, Franceschini MA, Kaipio JP, Arridge SR, Boas DA: Physiological system identification with the Kalman filter in diffuse optical tomography. Medical Image Computing and Computer-Assisted Intervention - Miccai 2005, Pt 2 2005, 3750: 649–656. full_textView ArticleGoogle Scholar
- Izzetoglu M, Chitrapu P, Bunce S, Onaral B: Motion artifact cancellation in NIR spectroscopy using discrete Kalman filtering. Biomed Eng Online 9: 16. 10.1186/1475-925X-9-16Google Scholar
- Prince S, Kolehmainen V, Kaipio JP, Franceschini MA, Boas D, Arridge SR: Time-series estimation of biological factors in optical diffusion tomography. Physics in Medicine and Biology 2003, 48: 1491–1504. 10.1088/0031-9155/48/11/301View ArticleGoogle Scholar
- Diamond SG, Huppert TJ, Kolehmainen V, Franceschini MA, Kaipio JP, Arridge SR, Boas DA: Dynamic physiological modeling for functional diffuse optical tomography. Neuroimage 2006, 30: 88–101. 10.1016/j.neuroimage.2005.09.016View ArticleGoogle Scholar
- Huppert TJ, Hoge RD, Diamond SG, Franceschini MA, Boas DA: A temporal comparison of BOLD, ASL, and NIRS hemodynamic responses to motor stimuli in adult humans. Neuroimage 2006, 29: 368–382. 10.1016/j.neuroimage.2005.08.065View ArticleGoogle Scholar
- Schmitz CH, Locker M, Lasker JM, Hielscher AH, Barbour RL: Instrumentation for fast functional optical tomography. Review of Scientific Instruments 2002, 73: 429–439. 10.1063/1.1427768View ArticleGoogle Scholar
- Bentsman J, Hong KS: Vibrational Stabilization of Nonlinear Parabolic-Systems with Neumann Boundary-Conditions. IEEE Transactions on Automatic Control 1991, 36: 501–507. 10.1109/9.75112MathSciNetView ArticleGoogle Scholar
- Bentsman J, Hong KS, Fakhfakh J: Vibrational Control of Nonlinear Time-Lag Systems - Vibrational Stabilization and Transient-Behavior. Automatica 1991, 27: 491–500. 10.1016/0005-1098(91)90106-CMathSciNetView ArticleGoogle Scholar
- Hong KS: An open-loop control for underactuated manipulators using oscillatory inputs: Steering capability of an unactuated joint. IEEE Transactions on Control Systems Technology 2002, 10: 469–480. 10.1109/87.998037View ArticleGoogle Scholar
- Sohn HC, Hong KT, Hong KS, Yoo WS: An adaptive LQG control for semi-active suspension systems. International Journal of Vehicle Design 2004, 34: 309–326. 10.1504/IJVD.2004.004060View ArticleGoogle Scholar
- Choi JY, Hong KS, Yang KJ: Exponential stabilization of an axially moving tensioned strip by passive damping and boundary control. Journal of Vibration and Control 2004, 10: 661–682. 10.1177/1077546304038103MathSciNetView ArticleGoogle Scholar
- Yang KJ, Hong KS, Matsuno F: Energy-based control of axially translating beams: Varying tension, varying speed, and disturbance adaptation. IEEE Transactions on Control Systems Technology 2005, 13: 1045–1054. 10.1109/TCST.2005.854368View ArticleGoogle Scholar
- Singh AK, Okamoto M, Dan H, Jurcak V, Dan I: Spatial registration of multichannel multi-subject fNIRS data to MNI space without MRI. Neuroimage 2005, 27: 842–851. 10.1016/j.neuroimage.2005.05.019View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<url>http://creativecommons.org/licenses/by/2.0</url>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.