This paper proposes a novel approach for extracting two intrinsic feature curves on hippocampal (HC) surfaces. The hippocampus is a key target of study in medical imaging, as it degenerates in conditions such as epilepsy and Alzheimer's disease (AD), but its structure is complex. To facilitate HC morphometry, we generate two intrinsic feature curves that describe their global geometries. For example, the separation of them captures thickness changes in HC surfaces, which can be used to effectively measure HC atrophy found in patients with AD. They also separate HC surfaces into upper and lower surface patches where intrinsic shape analysis using conformal modules can be carried out. Based on these curves, we further propose a parameterization of HC surfaces called the eigen-harmonic parameterization (EHP). EHP maps each HC surface onto a parameter domain and imposes longitudinal and azimuthal coordinates on each surface, which follow the gradient and level sets of its first nontrivial Laplace-Beltrami eigenfunction, respectively. Each tubular domain is constructed according to the geometry of an individual HC surface. This gives a parameter domain with much less geometric distortion compared to spherical parameterization. With EHP, all HC surfaces are automatically registered with intrinsic feature curves preserved and geometric distortions minimized. This allows shape analysis on any number of HC surfaces to be performed consistently. We studied geometric changes over time in 138 HC surfaces of patients with AD and normal subjects scanned at two different times. We successfully located areas with significantly different shape changes over time between the two groups.

Adolescent Idiopathic Scoliosis (AIS) characterized by the 3D spine deformity affects about 4% schoolchildren worldwide. Several studies have demonstrated the malfunctioning of postural balance, proprioception, and equilibrium control in patients with AIS. Since these functions are closely related to structures in and around the brainstem, the morphometry of the brainstem surface is of utmost importance. In this paper, we propose an effective method to accurately compute the registration between brainstem surfaces. Four consistent features, which describe the global geometry of the brainstem, are automatically extracted to guide the surface registration. Using the discrete Ricci flow method, brainstem surfaces are parameterized conformally onto the quadrilaterally-faced hexahedron, through which the brainstem registration can be obtained. Our registration algorithm can guarantee the exact landmark correspondence between brainstem surfaces. With the obtained registration, a shape energy can be defined to measure the local shape difference between different brainstem surfaces. We have tested our algorithms on 30 real brainstem surfaces extracted from MRIs of 15 normal subjects and 15 AIS patients. Experimental results show the efficacy of the proposed algorithm to register brainstem surfaces, which matches landmark features consistently. The computed registration can be used for the morphometry of brainstems.

We address the problem of surface inpainting, which aims to fill in holes or missing regions on a Riemann surface based on its surface geometry. In practical situation, surfaces obtained from range scanners often have holes where the 3D models are incomplete. In order to analyze the 3D shapes effectively, restoring the incomplete shape by filling in the surface holes is necessary. In this paper, we propose a novel conformal approach to inpaint surface holes on a Riemann surface based on its surface geometry. The basic idea is to represent the Riemann surface using its conformal factor and mean curvature. According to Riemann surface theory, a Riemann surface can be uniquely determined by its conformal factor and mean curvature up to a rigid motion. Given a Riemann surface S, its mean curvature H and conformal factor λ can be computed easily through its conformal parameterization. Conversely, given λ and H, a Riemann surface can be uniquely reconstructed by solving the Gauss-Codazzi equation on the conformal parameter domain. Hence, the conformal factor and the mean curvature are two geometric quantities fully describing the surface. With this λ-H representation of the surface, the problem of surface inpainting can be reduced to the problem of image inpainting of λ and H on the conformal parameter domain. Once λ and H are inpainted, a Riemann surface can be reconstructed which effectively restores the 3D surface with missing holes. Since the inpainting model is based on the geometric quantities λ and H, the restored surface follows the surface geometric pattern. We test the proposed algorithm on synthetic data as well as real surface data. Experimental results show that our proposed method is an effective surface inpainting algorithm to fill in surface holes on an incomplete 3D models based their surface geometry.