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
|
# emacs: -*- mode: python; py-indent-offset: 4; indent-tabs-mode: nil -*-
# vi: set ft=python sts=4 ts=4 sw=4 et:
import numpy as np
import pylab as pl
from nipy.neurospin.utils.simul_multisubject_fmri_dataset import \
surrogate_2d_dataset
pos = np.array([[10, 10],
[14, 20],
[23, 18]])
ampli = np.array([4, 5, 2])
# First generate some noiseless data
noiseless_data = surrogate_2d_dataset(nbsubj=1, noise_level=0, spatial_jitter=0,
signal_jitter=0, pos=pos, ampli=ampli)
pl.figure(figsize=(10, 3))
pl.subplot(1, 4, 1)
pl.imshow(noiseless_data[0])
pl.title('Noise-less data')
# Second, generate some group data, with default noise parameters
group_data = surrogate_2d_dataset(nbsubj=3, pos=pos, ampli=ampli)
pl.subplot(1, 4, 2)
pl.imshow(group_data[0])
pl.title('Subject 1')
pl.subplot(1, 4, 3)
pl.title('Subject 2')
pl.imshow(group_data[1])
pl.subplot(1, 4, 4)
pl.title('Subject 3')
pl.imshow(group_data[2])
|
