Evaluation of Three Cortical Surface Flattening Methods
Lili Ju1 , Monica Hurdal2 , Kelly Rehm3 , Bill Wood2 , Kirt Schaper4 , Josh Stern4 , David Rottenberg3,4
1Institute for Mathematics and its Applications, University of Minnesota, Minneapolis, MN 55455, USA, 2Department of Mathematics, Florida State University, Tallahassee, FL 32306, USA, 3Department of Radiology, University of Minnesota, Minneapolis, MN 55455, USA, 4Department of Neurology, University of Minnesota, Minneapolis, MN 55455, USA

Modeling & Analysis

Abstract
Intoduction
FreeSurfer (FS) [1], a popular software package for cortical surface flattening, explicitly minimizes metric distortion of the flattened surface. In contrast, LSCM [3] and CirclePack (CP) [2] are quasi-conformal flattening methods that preserve angles locally; however, spherical maps generated by these methods can be normalized by a Moebius transformation that minimizes metric distortion among the auomorphism group. All three methods can flatten user-defined patches as well as an entire cortical surface.

Methods
Topologically-correct triangulated surfaces of the left hemisphere were extracted from two high-resolution T1 MRI brain volumes (MNI, PENN) that had been parcellated by an expert neuroanatomist; parcellation of the extracted cortical sheet was guided by the volume parcellation using in-house software (ParcelMan [4]). For each brain, the left hemispheral cortex (equivalent to a topological sphere) and four lobar cortical patches (equivalent to topological discs) were flattened using FS, LSCM, and CP to create the corresponding spherical or planar maps. Measurements of angular and metric distortion of the flattened surfaces [3], which are invariant under the similarity transformations were used to describe the quality of resulting maps.

Results and Conclusions
Measurements of angular and metric distortion for the three methods, two brains, and five surfaces (the left hemispheral cortex and four lobar patches) are presented in Table 1, and frequency histograms of angular and metric distortion for the MNI left hemispheral cortex are illustrated in Figure 1. For the lobar patches FS clearly outperforms both conformal methods with regard to the preservation of metric information; however, for the MNI left hemispheral cortex FS and LSCM produce similar results for mean metric distortion: 31.72% and 36.07%, respectively. For all five surfaces LSCM is superior to the other methods with regard to the preservation of angular information, and CP and LSCM perform similarly with regard to metric distortion. Thus, by preserving angular (shape) information and adequately preserving metric information, LSCM may offer advantages to methods such as FS, which preserve only metric information, for flattening hemispheral cortical surfaces.

References
1. Fischl B., et al., Neuroimage 9:195-207, 1999.
2. Hurdal M.K., et al., LNCompSci, 1679:279-286, 1999.
3. Ju L.L., et al., NeuroImage 19:S864, 2003.
4. Rehm K., et al., HBM2004, submitted.

Acknowledgements
This work is supported in part by NIH grant MH57180 and NSF grant DMS101339.

Table 1. Angular (degrees) and metric (%) distortion of flat maps generated by FS, LSCM and CP.

MNI surface PENN Surface
Flat Maps S F O P T S F O P T
FS MAD 18.76 11.38 11.82 7.27 27.49 18.99 15.12 14.02 13.32 23.14

MMD 31.72 18.14 18.81 15.24 25.79 24.47 22.60 22.29 20.73 30.46
LSCM MAD 4.65 1.80 2.06 1.57 2.35 7.21 3.02 3.25 2.55 2.70

MMD 36.07 34.94 32.38 29.14 36.19 40.70 36.10 36.59 39.45 41.93
CP MAD 16.39 11.27 11.25 11.28 11.17 15.73 11.93 10.67 10.95 11.29

MMD 40.18 31.51 31.69 24.23 39.98 41.42 32.81 35.32 30.74 45.20
MAD, mean angular distortion; MMD, mean metric distortion; S, spherical map of the left hemispheral cortex; F,O,P,T, planar maps of frontal, occipital, parietal and temporal patches, respectively.




Figure 1. Frequency histograms illustrating the angular and metric distortion of spherical maps of the MNI left hemispheral cortex generated by FS (left), LSCM (middle) and CP(right). Top row, angular distortion; bottom row, metric distortion.