m/eeg forward problem & solutions brussels 2011 spm-m/eeg course january 2011 c. phillips,...
TRANSCRIPT
![Page 1: M/EEG forward problem & solutions Brussels 2011 SPM-M/EEG course January 2011 C. Phillips, Cyclotron Research Centre, ULg, Belgium](https://reader036.vdocument.in/reader036/viewer/2022062519/56649cb75503460f9497c7ca/html5/thumbnails/1.jpg)
M/EEG forward problem& solutions
Brussels 2011 SPM-M/EEG courseJanuary 2011
C. Phillips, Cyclotron Research Centre, ULg, Belgium
![Page 2: M/EEG forward problem & solutions Brussels 2011 SPM-M/EEG course January 2011 C. Phillips, Cyclotron Research Centre, ULg, Belgium](https://reader036.vdocument.in/reader036/viewer/2022062519/56649cb75503460f9497c7ca/html5/thumbnails/2.jpg)
Forward Problem
Inverse Problem
Source localisation in M/EEG
![Page 3: M/EEG forward problem & solutions Brussels 2011 SPM-M/EEG course January 2011 C. Phillips, Cyclotron Research Centre, ULg, Belgium](https://reader036.vdocument.in/reader036/viewer/2022062519/56649cb75503460f9497c7ca/html5/thumbnails/3.jpg)
Head anatomy
Solution by Maxwell’s equations
Head model : conductivity layoutSource model : current dipole
Forward problem
![Page 4: M/EEG forward problem & solutions Brussels 2011 SPM-M/EEG course January 2011 C. Phillips, Cyclotron Research Centre, ULg, Belgium](https://reader036.vdocument.in/reader036/viewer/2022062519/56649cb75503460f9497c7ca/html5/thumbnails/4.jpg)
tE
jB
tB
E
B
E
0Ohm’s law :
Continuity equation :
Ej
tj
Maxwell’s equations (1873)
![Page 5: M/EEG forward problem & solutions Brussels 2011 SPM-M/EEG course January 2011 C. Phillips, Cyclotron Research Centre, ULg, Belgium](https://reader036.vdocument.in/reader036/viewer/2022062519/56649cb75503460f9497c7ca/html5/thumbnails/5.jpg)
From Maxwell’s equations find:
where M are the measurements and f(.) depends on:•signal recorded, EEG or MEG•head model, i.e. conductivity layout adopted•source location•source orientation &
amplitude
Solving the forward problem
![Page 6: M/EEG forward problem & solutions Brussels 2011 SPM-M/EEG course January 2011 C. Phillips, Cyclotron Research Centre, ULg, Belgium](https://reader036.vdocument.in/reader036/viewer/2022062519/56649cb75503460f9497c7ca/html5/thumbnails/6.jpg)
From Maxwell’s equations find:
with f(.) as•an analytical solution
- highly symmetrical geometry, e.g. spheres, concentric spheres, etc.
- homogeneous isotropic conductivity•a numerical solution
- more general (but still limited!) head model
Solving the forward problem
![Page 7: M/EEG forward problem & solutions Brussels 2011 SPM-M/EEG course January 2011 C. Phillips, Cyclotron Research Centre, ULg, Belgium](https://reader036.vdocument.in/reader036/viewer/2022062519/56649cb75503460f9497c7ca/html5/thumbnails/7.jpg)
Con’s:•Human head is not spherical •Conductivity is not homogeneous and isotropic.
Analytical solution
Example: 3 concentric spheres
Pro’s:•Simple model•Exact mathematical solution•Fast calculation
![Page 8: M/EEG forward problem & solutions Brussels 2011 SPM-M/EEG course January 2011 C. Phillips, Cyclotron Research Centre, ULg, Belgium](https://reader036.vdocument.in/reader036/viewer/2022062519/56649cb75503460f9497c7ca/html5/thumbnails/8.jpg)
Overlapping spheres in MEG
Better but…
![Page 9: M/EEG forward problem & solutions Brussels 2011 SPM-M/EEG course January 2011 C. Phillips, Cyclotron Research Centre, ULg, Belgium](https://reader036.vdocument.in/reader036/viewer/2022062519/56649cb75503460f9497c7ca/html5/thumbnails/9.jpg)
Usually “Boundary Element Method” (BEM) :• Concentric sub-volumes of homogeneous and
isotropic conductivity,• Estimate values on the interfaces.
Numerical solution
Pro’s:•More exact head shape
modellingCon’s:•Mathematical approximations of solution numerical errors•Slow and intensive calculation
![Page 10: M/EEG forward problem & solutions Brussels 2011 SPM-M/EEG course January 2011 C. Phillips, Cyclotron Research Centre, ULg, Belgium](https://reader036.vdocument.in/reader036/viewer/2022062519/56649cb75503460f9497c7ca/html5/thumbnails/10.jpg)
Features of forward solution
Find forward solution (any):
with f(.)•linear in
•non-linear in
If source location is fixed, then
![Page 11: M/EEG forward problem & solutions Brussels 2011 SPM-M/EEG course January 2011 C. Phillips, Cyclotron Research Centre, ULg, Belgium](https://reader036.vdocument.in/reader036/viewer/2022062519/56649cb75503460f9497c7ca/html5/thumbnails/11.jpg)
SPM solutions
Source space & head model•Template Cortical Surface (TCS), in MNI space.•Canonical Cortical Surface (CCS) = TCS warped to subject’s anatomy•Subject’s Cortical Surface (SCS) = extracted from subject’s own structural image (BrainVisa/FreeSurfer)
![Page 12: M/EEG forward problem & solutions Brussels 2011 SPM-M/EEG course January 2011 C. Phillips, Cyclotron Research Centre, ULg, Belgium](https://reader036.vdocument.in/reader036/viewer/2022062519/56649cb75503460f9497c7ca/html5/thumbnails/12.jpg)
SPM solutions
Jérémie Mattout, Richard N. Henson, and Karl J. Friston, 2007, Canonical Source Reconstruction for MEG
![Page 13: M/EEG forward problem & solutions Brussels 2011 SPM-M/EEG course January 2011 C. Phillips, Cyclotron Research Centre, ULg, Belgium](https://reader036.vdocument.in/reader036/viewer/2022062519/56649cb75503460f9497c7ca/html5/thumbnails/13.jpg)
SPM solutions
Source space & head model•Template Cortical Surface (TCS), in MNI space.•Canonical Cortical Surface (CCS) = TCS warped to subject’s anatomy•Subject’s Cortical Surface (SCS) = extracted from subject’s own structural image (BrainVisa/FreeSurfer)
Forward solutions:•Single sphere•Overlapping spheres•Concentric spheres•BEM •… (new things get added to SPM & FieldTrip)
![Page 14: M/EEG forward problem & solutions Brussels 2011 SPM-M/EEG course January 2011 C. Phillips, Cyclotron Research Centre, ULg, Belgium](https://reader036.vdocument.in/reader036/viewer/2022062519/56649cb75503460f9497c7ca/html5/thumbnails/14.jpg)
Solving the Forward Problem is not exciting but necessary…
Solution depends on data at hand:•MEG vs. EEG:
single sphere or overlapping spheres vs. 3 concentric spheres or BEM
•Structural image available?CCS (good enough) or SCS (tedious work) vs. TCS (best proxy)
•Sensor location/head position available?
Conclusions
![Page 15: M/EEG forward problem & solutions Brussels 2011 SPM-M/EEG course January 2011 C. Phillips, Cyclotron Research Centre, ULg, Belgium](https://reader036.vdocument.in/reader036/viewer/2022062519/56649cb75503460f9497c7ca/html5/thumbnails/15.jpg)
Thank you for your attention