An algorithm for automated, noninvasive detection of cortical spreading depolarizations based on EEG simulations.

We present a novel signal processing algorithm for automated, noninvasive detection of Cortical Spreading Depolarizations (CSDs) using electroencephalography (EEG) signals and validate the algorithm on simulated EEG signals. CSDs are waves of neurochemical changes that suppress neuronal activity as they propagate across the brain's cortical surface. CSDs are believed to mediate secondary brain damage after brain trauma and cerebrovascular diseases like stroke. We address key challenges in detecting CSDs from EEG signals: (i) decay of high spatial frequencies as they travel from the cortical surface to the scalp surface; and (ii) presence of sulci and gyri, which makes it difficult to track the CSD waves as they travel across the cortex. Our algorithm detects and tracks "wavefronts" of the CSD wave, and stitches together data across space and time to decide on the presence of a CSD wave. To test our algorithm, we provide different models and complex patterns of CSD waves, including different widths of CSD suppressions, and use these models to simulate scalp EEG signals using head models of 4 subjects from the OASIS dataset. Our results suggest that the average width of suppression that a low-density EEG grid of 40 electrodes can detect is 1.1 cm, which includes a vast majority of CSD suppressions, but not all. A higher density EEG grid having 340 electrodes can detect complex CSD patterns as thin as 0.43 cm (less than minimum widths reported in prior works), among which single-gyrus propagation is the hardest to detect because of its small suppression area.