1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
|
#!/usr/bin/env python
# emacs: -*- mode: python; py-indent-offset: 4; indent-tabs-mode: nil -*-
# vi: set ft=python sts=4 ts=4 sw=4 et:
### ### ### ### ### ### ### ### ### ### ### ### ### ### ### ### ### ### ### ##
#
# See COPYING file distributed along with the PyMVPA package for the
# copyright and license terms.
#
### ### ### ### ### ### ### ### ### ### ### ### ### ### ### ### ### ### ### ##
"""
Generating Topography plots
===========================
Example demonstrating a topography plot."""
from mvpa.suite import *
# Sanity check if we have griddata available
externals.exists("griddata", raiseException=True)
# EEG example splot
P.subplot(1, 2, 1)
# load the sensor information from their definition file.
# This file has sensor names, as well as their 3D coordinates
sensors=XAVRSensorLocations(os.path.join(pymvpa_dataroot, 'xavr1010.dat'))
# make up some artifical topography
# 'enable' to channels, all others set to off ;-)
topo = N.zeros(len(sensors.names))
topo[sensors.names.index('O1')] = 1
topo[sensors.names.index('F4')] = 1
# plot with sensor locations shown
plotHeadTopography(topo, sensors.locations(), plotsensors=True)
# MEG example plot
P.subplot(1, 2, 2)
# load MEG sensor locations
sensors=TuebingenMEGSensorLocations(
os.path.join(pymvpa_dataroot, 'tueb_meg_coord.xyz'))
# random values this time
topo = N.random.randn(len(sensors.names))
# plot without additional interpolation
plotHeadTopography(topo, sensors.locations(),
interpolation='nearest')
if cfg.getboolean('examples', 'interactive', True):
# show all the cool figures
P.show()
"""
The ouput of the provided example should look like
.. image:: ../pics/ex_topo_plot.*
:align: center
:alt: Topography plot of MEG data
"""
|
