Parcellisation of MRI using glasser atlas 360

23 views
Skip to first unread message

Lou Natacha Fourneaux

unread,
Jul 11, 2024, 9:14:31 AM (5 days ago) Jul 11
to HCP-Users
Hello, 

I am trying to parcellate my MRI using the glasser atlas but I am not sure about the steps... 

First, I perform the FreeSurfer recon-all to extract the surfaces, brain volume, aseg....

Once I have them I just need to create a dense scalar with a wb-command -cifti-desne-scalar and then I can apply my parcellation. 

However, my wb-command -cifti-desne-scalar doesn't seem to be working and I can't find any information on how it works specifically... I input the file from the left surface and right surface so surf/lh.white.gii and surf/rh.white.gii  but I get the error "All data arrays (columns) in the file must have the same number of rows.  The first array (column) contains 134814 rows.  Array 2 contains 269624 rows." If you have any idea on what information I need to condense and how it would be great. 

Best, 
Lou

Tim Coalson

unread,
Jul 11, 2024, 4:53:18 PM (4 days ago) Jul 11
to hcp-...@humanconnectome.org
That sounds like you are trying to put the coordinate and triangle arrays together for some reason - .surf.gii files are not metric files, but all gifti formats share the base parsing code, so you might get that error when you give a .surf.gii file to a command that expects a metric file.  Cifti files do not store the surface coordinates or the triangle (vertex neighbor) information, they are intended only for values such as thickness, fMRI activation, etc.  You could put the surface vertex coordinates into one, but they would be interpreted as data to display using a palette, not as locations.  We continue to use .surf.gii to store surface coordinate and neighbor relation data, which is why spatial operations like -cifti-gradient require surfaces as separate arguments.

See also these explanations, though they don't explicitly make the point about surface coordinates/neighbor info:

That also sounds like a native mesh surface, which is difficult to compare across subjects directly.  We register and resample our surface data to fs_LR standard meshes before creating cifti files, so that we can do comparisons and math between subjects directly by indices.  The 32k resolution is matched to human 2mm data, and is what the Glasser MMP 1.0 is defined on, so that is probably the one you want to use.  You might consider using the HCP pipelines, possibly via the QuNex container to handle the software dependencies, in order to do our recommended alignment and cleanup preprocessing, and it will give you cifti versions of the data that you can then parcellate.

Tim


--
You received this message because you are subscribed to the Google Groups "HCP-Users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to hcp-users+...@humanconnectome.org.
To view this discussion on the web visit https://groups.google.com/a/humanconnectome.org/d/msgid/hcp-users/80d1a692-b9bc-4948-bcd9-77a51ea35d9dn%40humanconnectome.org.
Reply all
Reply to author
Forward
0 new messages