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
|
# emacs: -*- mode: python; py-indent-offset: 4; indent-tabs-mode: nil -*-
# vi: set ft=python sts=4 ts=4 sw=4 et:
"""
Example of a script that uses the BSA (Bayesian Structural Analysis)
-- nipy.neurospin.spatial_models.bayesian_structural_analysis --
module
Please adapt the image paths to make it work on your own data
Author : Bertrand Thirion, 2008-2009
"""
#autoindent
import numpy as np
import scipy.stats as st
import os.path as op
import tempfile
from nipy.neurospin.spatial_models.bsa_io import make_bsa_image
import get_data_light
# Get the data
#get_data_light.getIt()
nbsubj = 12
nbeta = 29
data_dir = op.expanduser(op.join('~', '.nipy', 'tests', 'data',
'group_t_images'))
mask_images = [op.join(data_dir,'mask_subj%02d.nii'%n)
for n in range(nbsubj)]
betas =[ op.join(data_dir,'spmT_%04d_subj_%02d.nii'%(nbeta,n))
for n in range(nbsubj)]
# set various parameters
subj_id = ['%04d' %i for i in range(12)]
theta = float(st.t.isf(0.01,100))
dmax = 5.
ths = 0 #2# or nbsubj/4
thq = 0.95
verbose = 1
smin = 5
swd = tempfile.mkdtemp()
method='simple'
# call the function
AF, BF = make_bsa_image(mask_images, betas, theta, dmax,
ths, thq, smin, swd, method, subj_id, '%04d'%nbeta)
# Write the result. OK, this is only a temporary solution
import pickle
picname = op.join(swd,"AF_%04d.pic" %nbeta)
pickle.dump(AF, open(picname, 'w'), 2)
picname = op.join(swd,"BF_%04d.pic" %nbeta)
pickle.dump(BF, open(picname, 'w'), 2)
print "Wrote all the results in directory %s"%swd
|
