Jim,

Attached are two scripts for getting the union or intersect of two label
files.

Nick

On Thu, 2009-08-27 at 20:32 -0500, James Porter wrote:
> Hello-
> 
> Is there a simple way to obtain the intersection of two labels, output
> into a third label? Or the intersection of a label and an .mgh file 
> output by mri_surfcluster?
> 
> For instance, I have a cluster defined by mri_glmfit/mri_surfcluster 
> that spans multiple anatomical regions, and I want to investigate the
> values from subsections of the cluster. That is to say, I want to run 
> mris_anatomical_stats on only the section of my cluster that covers, 
> say, precentral cortex but not caudal middle frontal or pars opercularis.
> 
> I've tried several strategies without success. The most promising seemed 
> to be just having mri_surfcluster output an aparc file, and then pass 
> that aparc file into mris_anatomical_stats. However, that aparc just 
> puts the very first region definition to every vertex. The same happens 
> if I run mri_label2aparc on the labels output by mri_surfcluster. This 
> happens no matter which LUT or ctab I use, i.e., with aparc.annot.ctab 
> the whole volume gets labeled "bankssts" and with FreeSurferColorLUT.txt 
> every vertex in the cluster volume gets labeled "Left-Cerebral-Exterior."
> 
> I guess there are a few issues in this email, but the most important one 
> to me is being able to run mris_anatomical_stats on the union of a 
> cluster file and the individual regions from the standard aparc atlas.
> 
> Any suggestions?
> 
> Thanks in advance,
> ---------
> Jim Porter
> Graduate Student
> Clinical Science & Psychopathology Research
> University of Minnesota
> 
> _______________________________________________
> Freesurfer mailing list
> Freesurfer@nmr.mgh.harvard.edu
> https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
> 
> 

Attachment: labels_union
Description: C-Shell script

Attachment: labels_intersect
Description: C-Shell script

_______________________________________________
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer

Reply via email to