[pymvpa] opening and closing a dataset

[Wolfgang Pauli]

Hi again,

I found the mistake. The mask file had incompatible dimensions. The code
does work fine.

Thanks,

Wolfgang

On Thu, Aug 11, 2016 at 3:12 PM, Wolfgang Pauli <wolfgang.m.pauli at gmail.com>
wrote:

> 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/f6aaf24f/attachment.html>