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:

1. Read the lh.area file directly [A, f]=read_curv('lh.area'); Total surface area: sum(A) = 86905.15mm^2
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