tensor decompositions for modelling epileptic seizures in eeg borbála hunyadi daan campsmaarten de...

19
Tensor decompositions for modelling epileptic seizures in EEG Borbála Hunyadi Daan Camps Maarten De Vos Laurent Sorber Sabine Van Huffel Wim Van Paesschen Lieven De Lathauwer

Upload: timothy-washington

Post on 15-Jan-2016

221 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Tensor decompositions for modelling epileptic seizures in EEG Borbála Hunyadi Daan CampsMaarten De Vos Laurent SorberSabine Van Huffel Wim Van Paesschen

Tensor decompositions for modelling epileptic seizures in EEGBorbála Hunyadi

Daan Camps Maarten De Vos

Laurent Sorber Sabine Van HuffelWim Van Paesschen Lieven De Lathauwer

Page 2: Tensor decompositions for modelling epileptic seizures in EEG Borbála Hunyadi Daan CampsMaarten De Vos Laurent SorberSabine Van Huffel Wim Van Paesschen

Outline

• Introductiono Epileptic seizureso EEG

• Tensor decompositionso CPDo BTD

• Signal modelo Oscillatory behaviour o Sum of exponentially damped sinusoids

• Simulation study

• Real EEG examples

• Conclusions

Page 3: Tensor decompositions for modelling epileptic seizures in EEG Borbála Hunyadi Daan CampsMaarten De Vos Laurent SorberSabine Van Huffel Wim Van Paesschen

Epilepsy

• Manifestation: o epileptic seizureso severe clinical symptoms

• Epileptic seizure:o abnormal, synchronous

activity of a large group of neurons

o Can be recorded in the EEG

Page 4: Tensor decompositions for modelling epileptic seizures in EEG Borbála Hunyadi Daan CampsMaarten De Vos Laurent SorberSabine Van Huffel Wim Van Paesschen

Seizures and EEG

• Repetitive, oscillatory pattern

• Evolution in o Amplitueo Frequencyo Topography

• Expert visual analysiso Determinte seizure type,

epilepsy syndromeo Important for proper

treatment

Page 5: Tensor decompositions for modelling epileptic seizures in EEG Borbála Hunyadi Daan CampsMaarten De Vos Laurent SorberSabine Van Huffel Wim Van Paesschen

Seizures and EEG

• Repetitive, oscillatory pattern

• Evolution in o Amplitueo Frequencyo Topography

• Expert visual analysiso Determinte seizure type,

epilepsy syndromeo Important for proper

treatment

• BUT! Artefacts...

Page 6: Tensor decompositions for modelling epileptic seizures in EEG Borbála Hunyadi Daan CampsMaarten De Vos Laurent SorberSabine Van Huffel Wim Van Paesschen

6

Nature of EEG

Mixture and

indirect measurement

EEG

Key considerations:Low SNR

Retrieve patterns of interest relying on a structured signal model

Appropriate representation and decomposition

s1

s2

sn

x1

xm

X = AS

Page 7: Tensor decompositions for modelling epileptic seizures in EEG Borbála Hunyadi Daan CampsMaarten De Vos Laurent SorberSabine Van Huffel Wim Van Paesschen

7

Tensor decompositions

= + + ... +

Ta

R

bR

cR

a2

b2

c2

a1

b1

c1

= + + ... +

TI1A1

c1

I2

I3B1T

A

c2

I2B2T

A2

I3 L2I1AR

cR

I2

I3BRTI1I1

I2I3 L

1LR

CPD:

BTD-(L,L,1):

Page 8: Tensor decompositions for modelling epileptic seizures in EEG Borbála Hunyadi Daan CampsMaarten De Vos Laurent SorberSabine Van Huffel Wim Van Paesschen

8

Signal model: oscillatory behaviourBTD of wavelet expanded EEG tensors

freq

uenc

ych

anne

l

time

CWT-CPD (Acar 2007, De Vos 2007)

CWT-BTD

Page 9: Tensor decompositions for modelling epileptic seizures in EEG Borbála Hunyadi Daan CampsMaarten De Vos Laurent SorberSabine Van Huffel Wim Van Paesschen

9

Signal model: sum of exp. damped sinusoidsBTD of Hankel expanded tensors

chan

nel

hankel

H-BTD (De Lathauwer, 2011)

Page 10: Tensor decompositions for modelling epileptic seizures in EEG Borbála Hunyadi Daan CampsMaarten De Vos Laurent SorberSabine Van Huffel Wim Van Paesschen

10

Simulation study

• 3 scenarioso Stationary ictal patterno Ictal pattern with evolving frequencyo Ictal pattern propagating towards remote brain regions

• Ictal pattern superimposed ono background EEG patterno muscle artefact (extracted from healthy EEG)

• Increasing noise levels (SNR: 1-0.1)

Page 11: Tensor decompositions for modelling epileptic seizures in EEG Borbála Hunyadi Daan CampsMaarten De Vos Laurent SorberSabine Van Huffel Wim Van Paesschen

11

Simulation studyStationary ictal pattern

• sinusoidal CWT-CPD or H-BTD-(1,2,2) is optimal

• CWT-BTD can be useful to model artefact sources

• H-BTD performs best to reconstruct time course

• All models equally good for retrieving the spatial map

Page 12: Tensor decompositions for modelling epileptic seizures in EEG Borbála Hunyadi Daan CampsMaarten De Vos Laurent SorberSabine Van Huffel Wim Van Paesschen

12

Simulation studyIctal pattern with evolving frequency

• CWT-BTD or H-BTD is the optimal model (L=?), while CPD cannot capture the frequency evolution

• CWT-BTD retrieves the TF matrices better than CPD (ICWT problem!)

• All models equally good in retrieving the localisation

Page 13: Tensor decompositions for modelling epileptic seizures in EEG Borbála Hunyadi Daan CampsMaarten De Vos Laurent SorberSabine Van Huffel Wim Van Paesschen

13

Simulation studyPropagating ictal pattern

• Fit a dipole on the reconstructed EEG

• CWT-BTD-(2,1,2) can reveal both sourceso fit 2 dipoleso fit 1 moving dipole

• CPD retrieves 1 source located in between the 2 simulated sources

Page 14: Tensor decompositions for modelling epileptic seizures in EEG Borbála Hunyadi Daan CampsMaarten De Vos Laurent SorberSabine Van Huffel Wim Van Paesschen

14

Clinical examplesSevere artefact

Page 15: Tensor decompositions for modelling epileptic seizures in EEG Borbála Hunyadi Daan CampsMaarten De Vos Laurent SorberSabine Van Huffel Wim Van Paesschen

15

Clinical examplesEvolution in frequency

Page 16: Tensor decompositions for modelling epileptic seizures in EEG Borbála Hunyadi Daan CampsMaarten De Vos Laurent SorberSabine Van Huffel Wim Van Paesschen

16

Clinical examplesSpatial evolution

Page 17: Tensor decompositions for modelling epileptic seizures in EEG Borbála Hunyadi Daan CampsMaarten De Vos Laurent SorberSabine Van Huffel Wim Van Paesschen

17

Conclusion

• CWT-CPDo Model stationary sourceso Onset localisation

• CWT-BTDo Sources with evolving frequency or spatial distributiono High power, complex artefacts

• H-BTDo Seizure with fixed topography with arbitrary time courseo Precise reconstruction of time course

Page 18: Tensor decompositions for modelling epileptic seizures in EEG Borbála Hunyadi Daan CampsMaarten De Vos Laurent SorberSabine Van Huffel Wim Van Paesschen

18

Future work

• Automatic model selection

• Applications:o Onset localisation:

• automatic model selection is needed• Test on large real EEG dataset

o Seizure detection: • find optimal model with trial-error and use the model to detect

subsequent seizures

Page 19: Tensor decompositions for modelling epileptic seizures in EEG Borbála Hunyadi Daan CampsMaarten De Vos Laurent SorberSabine Van Huffel Wim Van Paesschen

19

Thank you!

• Any questions?