Improved surface feature recognition in CT images is provided by extracting a triangulated mesh representation of the surface of interest. Shape operators are computed at each vertex of the mesh from finite differences of vertex normals. The shape operators at each vertex are smoothed according to an iterative weighted averaging procedure. Principal curvatures at each vertex are computed from the smoothed shape operators. Vertices are marked as maxima and/or minima according to the signs of the principal curvatures. Vertices marked as having the same feature type are clustered together by adjacency on the mesh to provide candidate patches. Feature scores are computed for each candidate patch and the scores are provided as output to a user or for further processing.