Spike cluster analysis in neocortical localization related epilepsy yields clinically significant equivalent source localization results in magnetoencephalogram (MEG).
OBJECTIVE
In magnetoencephalogram (MEG) recordings of patients with epilepsy several types of sharp transients with different spatiotemporal distributions are commonly present. Our objective was to develop a computer based method to identify and classify groups of epileptiform spikes, as well as other transients, in order to improve the characterization of irritative areas in the brain of epileptic patients.


METHODS
MEG data centered on selected spikes were stored in signal matrices of C channels by T time samples. The matrices were normalized and euclidean distances between spike representations in vector space R(CxT) were input to a Ward's hierarchical clustering algorithm.


RESULTS
The method was applied to MEG data from 4 patients with localization-related epilepsy. For each patient, distinct spike subpopulations were found with clearly different topographical field maps. Inverse computations to selected spike subaverages yielded source solutions in agreement with seizure classification and location of structural lesions, if present, on magnetic resonance images.


CONCLUSIONS
With the proposed method a reliable categorization of epileptiform spikes is obtained, that can be applied in an automatic way. Computation of subaverages of similar spikes enhances the signal-to-noise ratio of spike field maps and allows for more accurate reconstruction of sources generating the epileptiform discharges.
D000328:Adult; D001931:Brain Mapping; D016000:Cluster Analysis; D003625:Data Collection; D004574:Electromagnetic Fields; D004827:Epilepsy / Q000145:classification / Q000473:pathology / Q000503:physiopathology*; D005071:Evoked Potentials; D005260:Female; D006801:Humans; D007206:Individuality; D008279:Magnetic Resonance Imaging / Q000379:methods; D015225:Magnetoencephalography; D008297:Male; D019579:Neocortex / Q000503:physiopathology*; D015203:Reproducibility of Results; D012815:Signal Processing, Computer-Assisted; D013997:Time Factors