[pymvpa] opening and closing a dataset

[Wolfgang Pauli]

Hi,

Finally found some time to work with mvpa again. However, despite my
impression that I am now decently familiar with pymvpa2, I realized that I
can not open and close a file.

This is what I do:

dhandle = OpenFMRIDataset(data_path)
ds = dhandle.get_model_bold_dataset(model_id=1,
                                        subj_id=subj,
                                        noinfolabel='rest',
                                        mask=atlas_mask_file)
tmp = map2nifti(ds, data=ds.samples[0])
tmp.to_filename("one_ds_sample.nii.gz")

Shouldn't this work??? The output is complete garbage. The image as the
correct dimension, but the data is just frantically smeared all across
space.

ds.a.mapper looks good to me:
ChainMapper(nodes=[FlattenMapper(shape=(80, 91, 40), auto_train=True,
space='voxel_indices'), StaticFeatureSelection(dshape=(291200,),
slicearg=array([False, False, False, ..., False, False, False],
dtype=bool))])

Thanks!

Wolfgang
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.alioth.debian.org/pipermail/pkg-exppsy-pymvpa/attachments/20160811/a97bd1ac/attachment.html>