There are a couple of things happening.
The lh.area has the vertices in the medial wall zeroed so #1 below does not include all of the unknown area. Here the "medial wall" is defined as anything outside of lh.cortex.label
The mris_anatomical_stats (#3 ) command below does include the medial wall.
I'm not sure about #2.
Note that there is a #4 which is from ?.aparc.stats:
# Measure Cortex, WhiteSurfArea, White Surface Total Area, 85048.6, mm^2
It will not give the same answer as #1 because #4 computes the total over all all vertices that have an annotation (implicitly defining the medial wall as the unknown in the aparc which does not quite line up with lh.cortex.label)
#4 is what I would use
doug
On 09/30/2013 04:38 PM, XJ Kang wrote:
Hi Experts,
I am trying to figure out the total surface area of a hemisphere. I found out 3 ways to read the area data:
- Read the lh.area file directly [A, f]=read_curv('lh.area'); Total surface area: sum(A) = 86905.15mm^2
- Read from the ~/stats/lh.curv.stats (generated by
mris_curvature_stats in recon-all): Raw Total Surface Area: 87347.86719 mm^2 3. Read from the results: % mris_anatomical_stats sub lh white total surface area = 91046 mm^2
First, the total area includes the cutting area (unknown region?), right? Secondly, which total area is more reasonable? Thanks.
XJ Kang
Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer