A numerical scheme for the curvature equation near the singularities

Abstract

In this paper we propose a modification in the usual numerical method for computing the solutions of the curvature equation in the plane . This modification takes place near the singularities of the image. We propose to use zero as the vertical speed at a saddle point and, at an extremum, the geometric mean of the eigenvalues of the Hessian matrix. This modification is theoretically justified and the preliminary experimental results show that it makes the algorithm more reliable. © 2005 Springer Science + Business Media, Inc.