[FieldTrip] possible bug with either ft_volumelookup or AAL atlas

Tyler Grummett tyler.grummett at flinders.edu.au
Thu May 29 09:26:39 CEST 2014


Thank you Eelke,

I also noticed that with the AFNI atlas, the labels in the file dont match up to the labels
that come out of ft_prepare_sourcemodel.

when I say label I mean data.( label)

I tested it with the AAL atlas and the labels in the file and the data.( label) match up perfectly.

Tyler

*************************

Tyler Grummett ( BBSc, BSc(Hons I))
PhD Candidate
Brain Signals Laboratory
Flinders University
Rm 5A301
Ext 66124

________________________________________
From: fieldtrip-bounces at science.ru.nl <fieldtrip-bounces at science.ru.nl> on behalf of Eelke Spaak <eelke.spaak at donders.ru.nl>
Sent: Thursday, 29 May 2014 4:30 PM
To: FieldTrip discussion list
Subject: Re: [FieldTrip] possible bug with either ft_volumelookup or AAL        atlas

Hi Tyler,

That sounds like a bug to me. I've filed it on our bugzilla:
http://bugzilla.fcdonders.nl/show_bug.cgi?id=2594

Best,
Eelke

On 29 May 2014 03:56, Tyler Grummett <tyler.grummett at flinders.edu.au> wrote:
> Hello fieldtrippers,
>
>
> I just wanted to check if this is a bug or not.
>
>
> When running the following code:    atlas = ft_read_atlas( fullfile(
> matlabrootpath, 'Matlab/fieldtrip/template/atlas/aal/ROI_MNI_V4.nii'))
>
>
> It gives you the following:
>
>
>             dim: [91 109 91]
>             hdr: [1x1 struct]
>       transform: [4x4 double]
>            unit: 'mm'
>          tissue: [91x109x91 double]
>     tissuelabel: {1x116 cell}
>        coordsys: 'mni'
>
> Take note that tissuelabel is 1x116
>
> When running the following code through ft_volumelookup:
>
>             cfg = [];
>             cfg.atlas = atlas;
>             cfg.inputcoord = 'mni';
>             cfg.maskparameter = 'inside';
>             labels = ft_volumelookup( cfg, source);
>
> An error occurs:
>
> Error using cat
> Dimensions of matrices being concatenated are not consistent.
>
> Error in ft_volumelookup (line 313)
>     labels.name = cat(1, labels.name, atlas.([fn{k},'label']));
>
> If you invert tissuelabel, so that it becomes 116x1, then the code works.
>
> For some reason ft_volumelookup doesnt like a 1xY cell.
>
> Thank you,
>
> Tyler
>
>
>
> *************************
>
> Tyler Grummett ( BBSc, BSc(Hons I))
> PhD Candidate
> Brain Signals Laboratory
> Flinders University
> Rm 5A301
> Ext 66124
>
> _______________________________________________
> fieldtrip mailing list
> fieldtrip at donders.ru.nl
> http://mailman.science.ru.nl/mailman/listinfo/fieldtrip
_______________________________________________
fieldtrip mailing list
fieldtrip at donders.ru.nl
http://mailman.science.ru.nl/mailman/listinfo/fieldtrip




More information about the fieldtrip mailing list