Accurate interpretation of functional MRI (fMRI) signals requires knowledge of the relationship between the hemodynamic response and the neuronal activity that underlies it. Here we address the question of coupling between pre-and postsynaptic neuronal activity and the hemodynamic response in rodent somatosensory (Barrel) cortex in response to single-whisker deflection. Using full-field multiwavelength optical imaging of hemoglobin oxygenation and electrophysiological recordings of spiking activity and local field potentials, we demonstrate that a point hemodynamic measure is influenced by neuronal activity across multiple cortical columns. We demonstrate that the hemodynamic response is a spatiotemporal convolution of the neuronal activation. Therefore, positive hemodynamic response in one cortical column might be explained by neuronal activity not only in that column but also in the neighboring columns. Thus, attempts at characterizing the neurovascular relationship based on point measurements of electrophysiology and hemodynamics may yield inconsistent results, depending on the spatial extent of neuronal activation. The finding that the hemodynamic signal observed at a given location is a function of electrophysiological activity over a broad spatial region helps explain a previously observed increase of local vascular response beyond the saturation of local neuronal activity. We also demonstrate that the oxy-and total-hemoglobin hemodynamic responses can be well approximated by space-time separable functions with an antagonistic center-surround spatial pattern extending over several millimeters. The surround ''negative'' hemodynamic activity did not correspond to observable changes in neuronal activity. The complex spatial integration of the hemodynamic response should be considered when interpreting fMRI data.Barrel cortex ͉ blood oxygenation ͉ intrinsic signals ͉ optical imaging T he advent of noninvasive imaging methods such as functional MRI (fMRI) has made it possible to obtain spatial maps of hemodynamic ''activation'' in the human brain under a variety of conditions (1, 2). However, the indirect and poorly understood nature of the coupling between these hemodynamic signals and the underlying neuronal activity has greatly limited the interpretability of neuroimaging results. Recently, several groups have attempted to characterize this coupling in the form of a linear or nonlinear neurovascular ''transfer function '' (3-8). In principle, if such a function could be defined, it would provide a basis for inferring time-averaged local neuronal activity based on hemodynamic measurements. Furthermore, it would permit more accurate integration of hemodynamic imaging methods with noninvasive electrophysiological recordings such as electroencephalography and magnetoencephalography (9, 10).In a previous publication (3), using simultaneous spectroscopic optical imaging and electrophysiological measurements in rodent somatosensory cortex during brief and spatially localized stimuli, we found a strongly nonlinear r...