Hello Freesurfers,
I was wondering if someone might have information for how I can measure overlapping areas between a cluster that was generated as a results of a cortical thickness analysis, with a portion of a Yeo Network atlas. Specifically, I would like to calculate the percentage that the cortical thickness significant cluster (cache.th13.abs.sig.ocn.annot) overlaps with one region of the "lh.Yeo2011_7Networks_N1000.annot" overlay.
I am not sure if there is a way to trace the two in tksurfer, and measure their overlapping area in that manner.
Any help would be greatly appreciated.
Thank you,
Arkadiy
THere is not an easy way to do this from the command line. You can do it from matlab, something like
[Yvertices Ylabel Ycolortable] = read_annotation('lh.Yeo2011_7Networks_N1000.annot')
[Cvertices Clabel Ccolortable] = read_annotation('cache.th13.abs.sig.ocn.annot');
Then find the vertices in Ylabel that correspond to the label you want (number Ny), then find those in Clabel that correspond to the cluster your want (Nc), then those vertices that are in both (Nyc). You can then compute the overlap in the way you want.
Look at the help for read_annotation to see how to get the label you want
On 08/24/2016 10:07 PM, Arkadiy Maksimovskiy wrote:
Hello Freesurfers,
I was wondering if someone might have information for how I can measure overlapping areas between a cluster that was generated as a results of a cortical thickness analysis, with a portion of a Yeo Network atlas. Specifically, I would like to calculate the percentage that the cortical thickness significant cluster (cache.th13.abs.sig.ocn.annot) overlaps with one region of the "lh.Yeo2011_7Networks_N1000.annot" overlay.
I am not sure if there is a way to trace the two in tksurfer, and measure their overlapping area in that manner.
Any help would be greatly appreciated.
Thank you,
Arkadiy
Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
freesurfer@nmr.mgh.harvard.edu