Magnetoencephalography (MEG) and electroencephalography (EEG) use non-invasive sensors to detect neural currents. Since the contribution of superficial neural sources to the measured M/EEG signals are orders-of-magnitude stronger than the contribution of subcortical sources, most MEG and EEG studies have focused on cortical activity. Subcortical structures, however, are centrally involved in both healthy brain function as well as in many neurological disorders such as Alzheimer's disease and Parkinson's disease. In this paper, we present a method that can separate and suppress the cortical signals while preserving the subcortical contributions to the M/EEG data. The resulting signal subspace of the data mainly originates from subcortical structures. Our method works by utilizing short-baseline planar gradiometers with short-sighted sensitivity distributions as reference sensors for cortical activity. Since the method is completely data-driven, forward and inverse modeling are not required. In this study, we use simulations and auditory steady state response experiments in a human subject to demonstrate that the method can remove the cortical signals while sparing the subcortical signals. We also test our method on MEG data recorded in an essential tremor patient with a deep brain stimulation implant and show how it can be used to reduce the DBS artifact in the MEG data by ~ 99.9% without affecting low frequency brain rhythms.
Keywords: Electroencephalography; Magnetoencephalography; Signal processing; Spatial filtering; Subcortical imaging; Temporal subspace projection.