Hi Bruce,
Sorry, I did a poor job of describing the problem.
We are trying to do an analysis of subject specific abnormalities on the surface. Instead of a traditional analysis assuming that patients have homogenous abnormalities (for example cortical thinning in anterior cingulate) relative to a control group, we want to count the number of abnormalities over the entire surface in each subject. On the surface, abnormalities can be defined as average thickness in a small region that is >2 SD away from mean thickness in a control group in that same region. Following that, we can compare the average number of subject specific abnormalities between groups.
This method has been used to identify white matter abnormalities in patient populations in which heterogenous abnormalities are expected, such as TBI patients. However, it requires an a priori decision about what constitutes an abnormality. In this case, we need to specify it in number of vertices.
We have previously used 128 ul as our defined abnormality size in the volume. Our first thought was to use that as a starting point and translate the volume to surface area or number of vertices.
Is there a reasonable way to approximate how many voxels get translated into surface space? And is there a reasonable way to approximate how much area each data point would have on the surface?
One approach we are considering is to calculate the average gray matter volume, surface area, and number of vertices from the aparc.stats file. With this information we'd have an approximate relationship between volume, surface area, and # vertices. We are curious if this seems like a reasonable choice, or if you have any thoughts on how to go about this.
Thanks a lot,
Tim