[med-svn] [python-mne] 343/353: ENH : plot something in plot_lcmv_beamformer_volume.py

[Yaroslav Halchenko]

This is an automated email from the git hooks/post-receive script.

yoh pushed a commit to tag 0.4
in repository python-mne.

commit 52a67fb538633adae682fdb63170954d139782e3
Author: Alexandre Gramfort <alexandre.gramfort at inria.fr>
Date:   Sun Aug 5 14:32:38 2012 +0200

    ENH : plot something in plot_lcmv_beamformer_volume.py
---
 examples/inverse/plot_lcmv_beamformer_volume.py | 21 +++++++++++++++++++++
 1 file changed, 21 insertions(+)

diff --git a/examples/inverse/plot_lcmv_beamformer_volume.py b/examples/inverse/plot_lcmv_beamformer_volume.py
index 053ac1f..9d7fee7 100644
--- a/examples/inverse/plot_lcmv_beamformer_volume.py
+++ b/examples/inverse/plot_lcmv_beamformer_volume.py
@@ -15,6 +15,8 @@ Freeview.
 
 print __doc__
 
+import numpy as np
+import pylab as pl
 import mne
 from mne.datasets import sample
 from mne.fiff import Raw, pick_types
@@ -64,3 +66,22 @@ stc.crop(0.0, 0.2)
 # Save result in a 4D nifti file
 img = mne.save_stc_as_volume('lcmv_inverse.nii.gz', stc,
         forward['src'], mri_resolution=False)  # True for full MRI resolution
+
+# plot result (one slice)
+pl.close('all')
+data = img.get_data()
+coronal_slice = data[:, 10, :, 60]
+pl.figure()
+pl.imshow(np.ma.masked_less(coronal_slice, 1), cmap=pl.cm.Reds,
+          interpolation='nearest')
+pl.colorbar()
+pl.contour(coronal_slice != 0, 1, colors=['black'])
+pl.xticks([])
+pl.yticks([])
+
+# plot source time courses with the maximum peak amplitudes
+pl.figure()
+pl.plot(stc.times, stc.data[np.argsort(np.max(stc.data, axis=1))[-40:]].T)
+pl.xlabel('Time (ms)')
+pl.ylabel('LCMV value')
+pl.show()

-- 
Alioth's /usr/local/bin/git-commit-notice on /srv/git.debian.org/git/debian-med/python-mne.git