9th International Conference on Human Brain Mapping, New York City, U.S.A., June 18-22, 2003
Cortical Surface Flattening Using Discrete Conformal Mapping with Minimal Metric Distortion
Lili Ju1, Yongwu Shao2, Josh Stern3, Kelly Rehm4, Kirt Schaper3, Monica Hurdal5, and David Rottenberg3,4

1Institute for Mathematics and its Applications, University of Minnesota, Minneapolis, MN 55455, USA
2Department of Statistics, University of Minnesota, Minneapolis, MN 55455, USA
3Department of Neurology, University of Minnesota, Minneapolis, MN 55455, USA
4Department of Radiology, University of Minnesota, Minneapolis, MN 55455, USA
5Department of Mathematics, Florida State University, Tallahassee, FL 32306, USA

Abstract
Although flattening a cortical surface necessarily introduces metric distortion due to the non-constant Gaussian curvature of the surface, the Riemann Mapping Theorem states that continuously differentiable surfaces can be mapped without angular distortion. Several techniques have been proposed for flattening polygonal representations of surfaces while substantially minimizing metric distortion [1,2], and methods for conformal flattening of polygonal surfaces have also been proposed [3,4]. We describe an efficient method for generating conformal flat maps of triangulated surfaces while minimizing metric distortion within the class of conformal maps. Our method, which controls both angular and metric distortion, involves the solution of a linear system [5] and a small scale nonlinear minimization. It can be applied to user-defined "patches" or to an entire cortical surface.

Methods
Given any pair of surfaces, P and Q, that are both topological disks, and a fixed homeomorphism between their boundaries, there exists a unique harmonic map from P to Q that minimizes the Dirichlet energy and is conformal. We utilize a discretization scheme and optimization-based algorithm for computing the harmonic mapping described in [5]. Since there is always a class of conformal mappings from P to Q related to each other by an automorphism of Q, our algorithm then uses numerical methods to search for that conformal mapping which minimizes the metric distortion of the resulting mesh. Conformal mappings between topological spheres proceed via the mapping for disks using the same technique as [3]. A vertex v and all edges containing it are removed from the input mesh, after that we generate a discrete conformal map from the pruned mesh to the unit disk D; this map is then stereographically projected to the unit sphere S while v is mapped to the "north pole" of S, and the resulting map is again optimized by choosing the S automorphism that minimizes metric distortion.

Results and Conclusions
We created discrete conformal flat maps on the unit sphere from topologically-correct human cerebellar and cerebral cortices; we also flattened selected cortical patches from the cerebellum (Lobes IV and V) and cerebrum (occipital lobe). Mappings obtained for the cerebellar cortex are illustrated in Figure 1. Values for mean angular and normalized metric distortion are compared to those obtained using FreeSurfer [1] in Table 1. With our conformal method, metric distortion is no worse than that produced by FreeSurfer; moreover, angular information is preserved, and computation time greatly reduced (20 minutes for the spherical flattening of cerebral cortex vs. 10 hours using FreeSurfer).

Table 1.  Computational results for discrete flat maps of human cerebellar and cerebral cortices produced by our method and by FreeSurfer. Flat maps of the [entire] cerebellar and cerebral cortices on S (A); Flat maps of cortical patches (lobes IV and V for cerebellum, occipital lobe for cerebrum)in a predefined planar region (B) and on D (C).
Cerebellar Cortex Cerebral Cortex
Our MethodFreeSurfer Our MethodFreeSurfer
Flat maps ABCA ABCA
Mean angular distortion 4.44°6.12°11.51°23.37° 6.50°2.66°8.15°18.76°
Mean normalized metric distortion 40.44%23.47%35.53%46.16% 36.52%25.37%21.79%30.81%


Figure 1.  Conformal flat maps of the human cerebellar cortex. Left panel : Parcellated surface of the cerebellum and its flat map on S; Right panel : Patch from cerebellar Lobes IV and V and its flat maps in a predefined planar region and on D.

References
1. Fischl B., et al., Neuroimage 9:195-207, 1999.
2. Drury H.M. et al., JCN, 8:1-28, 1996.
3. Hurdal M.K., et al., LNCompSci, 1679:279-286, 1999.
4. Angenent S., et al., TMI 18:700-711, 1999.
5. Eck M., et al., SIGGRAPH95, 173-182.

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

NeuroImage, Volume 19, Number 2, Supplement 1, Page S44, CD-Rom Abstract 864, 2003