[pymvpa] Antwort: Re: Antwort: Re: Loading normalized data with 'get_model_bold_dataset' fails

Yaroslav Halchenko debian at onerussian.com
Fri Jun 2 15:21:08 UTC 2017


On Fri, 02 Jun 2017, lydia.hellrung at econ.uzh.ch wrote:

>    No, they all have 195 volumes, but the geometry information is apparently
>    saved differently.
>    The output says:

>    ./sub002/BOLD/task001_run001/bold_irbadvols.nii.gz    float32 [112, 112,
>     27, 195] 1.96x1.96x4.06x2.00  #exts: 1 sform
>    ./sub002/BOLD/task001_run001/bold_iwrbadvols.nii.gz   float32 [105, 127,
>    105, 195] 1.50x1.50x1.50x0.00  #exts: 1

>    while the first file could be loaded, the second one not. Could that be an
>    issue? And in case, what would be needed?

well -- you are saying   flavor='iwrbadvols'  so only that one should be
loaded, and I guess the 'brain.nii.gz'  mask is of the 'compatible'
dimensionality I guess.  otherwise I would have expected it to fail
earlier.

But the above shows the "culprit" --   1.50x1.50x1.50x0.00   so your TR
there is set to 0 sec...  so time indices are degenerate and I guess
that played out somewhere (may be compressed everything into a single
sample or smth like that)

You could fix it  via a quick python code using nibabel or using
nifti_tool  ... --help gives

       e.g. modifying nifti_1_header fields:
            -mod_field descrip 'toga, toga, toga'
            -mod_field qoffset_x 19.4 -mod_field qoffset_z -11
            -mod_field pixdim '1 0.9375 0.9375 1.2 1 1 1 1'

-- 
Yaroslav O. Halchenko
Center for Open Neuroscience     http://centerforopenneuroscience.org
Dartmouth College, 419 Moore Hall, Hinman Box 6207, Hanover, NH 03755
Phone: +1 (603) 646-9834                       Fax: +1 (603) 646-1419
WWW:   http://www.linkedin.com/in/yarik        



More information about the Pkg-ExpPsy-PyMVPA mailing list