An Invariant Shape Representation Using the Anisotropic Helmholtz Equation

Analyzing geometry of sulcal curves on the human cortical surface requires a shape representation invariant to Euclidean motion. We present a novel shape representation that characterizes the shape of a curve in terms of a coordinate system based on the eigensystem of the anisotropic Helmholtz equation. This representation has many desirable properties: stability, uniqueness and invariance to scaling and isometric transformation. Under this representation, we can find a point-wise shape distance between curves as well as a bijective smooth point-to-point correspondence. When the curves are sampled irregularly, we also present a fast and accurate computational method for solving the eigensystem using a finite element formulation. This shape representation is used to find symmetries between corresponding sulcal shapes between cortical hemispheres. For this purpose, we automatically generate 26 sulcal curves for 24 subject brains and then compute their invariant shape representation. Left-right sulcal shape symmetry as measured by the shape representation's metric demonstrates the utility of the presented invariant representation for shape analysis of the cortical folding pattern.