Abstract. Brain Cortical surface registration is required for inter-subject studies of functional and anatomical data. Harmonic mapping has been applied for brain mapping, due to its existence, uniqueness, regularity and numerical stability. In order to improve the registration accuracy, sculcal landmarks are usually used as constraints for brain registration. Unfortunately, constrained harmonic mappings may not be diffeomorphic and produces invalid registration. This work conquers this problem by changing the Riemannian metric on the target cortical surface to a hyperbolic metric, so that the harmonic mapping is guaranteed to be a diffeomorphism. With landmarks delineated as boundary condition, it is possible to be integrated with landmark matching. The computational algorithms are based on the Ricci flow method and hyperbolic heat diffusion. Experimental results demonstrate that, by changing the Riemannian metric, the registrations have higher qualities in terms of landmark alignment, curvature matching, area distortion and overlapping of region of interests.