The spatiotemporal dynamics of cortical oscillations across human brain regions remain poorly understood because of a lack of adequately validated methods for reconstructing such activity from noninvasive electrophysiological data. We present a novel adaptive spatial filtering algorithm optimized for robust source time-frequency reconstruction from magnetoencephalography (MEG) and electroencephalography (EEG) data. The efficacy of the method is demonstrated with real MEG data from a self-paced finger movement task. The algorithm reliably reveals modulations both in the beta band (12-30 Hz) and a high gamma band (65-90 Hz) in sensorimotor cortex. The performance is validated by both across-subjects statistical comparisons and by intracranial electrocorticography (ECoG) data from two epilepsy patients. We also revealed observed high gamma activity in the cerebellum. The proposed algorithm is highly parallelizable and runs efficiently on modern high performance computing clusters. This method enables non-invasive five-dimensional imaging of space, time, and frequency activity in the brain and renders it applicable for widespread studies of human cortical dynamics.