The relationship between electrophysiological and functional magnetic resonance imaging (fMRI) signals remains poorly understood. To date, studies have required invasive methods and have been limited to single functional regions and thus cannot account for possible variations across brain regions. Here, we present a method that uses functional magnetic resonance imaging (fMRI) data and singe-trial electroencephalography (EEG) analyses to assess the spatial and spectral dependencies between the blood oxygenation level dependent (BOLD) responses and the non-invasively estimated local field potentials (eLFPs) over a wide range of frequencies (0-256Hz) throughout the entire brain volume. This method was applied in a study where human subjects completed separate fMRI and EEG sessions while performing a passive visual task. Intracranial LFPs were estimated from the scalp-recorded data using the ELECTRA source model (Grave de Peralta Menendez et al. 2000). We compared statistical images from BOLD signals with statistical images of each frequency of the eLFPs. In agreement with previous studies in animals, we found a significant correspondence between LFP and BOLD statistical images in the gamma band (44-78Hz) within primary visual cortices. In addition, significant correspondence was observed at low frequencies (<14Hz) and also at very high frequencies (>100Hz). Effects within extrastriate visual areas showed a different correspondence that not only included those frequency ranges observed in primary cortices, but also additional frequencies. Results, therefore, suggest that the relationship between electrophysiological and hemodynamic signals thus might vary both as a function of frequency and anatomical region.