Curvature estimation for enhancement of crossing curves

In this paper we describe a method for estimating curvature of elongated structures in images. The curvature estimation is performed on an invertible orientation score, which is a 3D entity obtained from a 2D image by convolution with a rotating kernel. By considering the group structure we can define left-invariant derivatives, which are essential to construct operations on the orientation score that amount to rotationally invariant operations on the corresponding image. The problem of estimating curvature of an oriented structure is stated as a minimization problem, which can be solved by eigenvector analysis of a matrix constructed from the non-symmetric Hessian matrix. The experiments show the method performs well for a wide range of curvatures and noise levels. The method clearly outperforms a related curvature estimation method by Van Ginkel et al. that tends to give estimates that are too small. We show how we can incorporate the curvature estimate in our method for coherence-enhancing diffusion in orientation scores. This method has superior performance in enhancing crossing contours, which is demonstrated on medical images.