In this paper, we present a novel method to estimate curvature of iso gray-level surfaces in gray-value images. Our method succeeds where standard isophote curvature estimation methods fail. There is neither a segmentation of the surface needed nor a parametric model assumed. Our estimator works on the orientation field of the surface. This orientation field and a description of local structure is obtained by the Gradient Structure Tensor. The estimated orientation field has discontinuities mod pi. It is mapped via the Knutsson mapping to a continuous representation. The principal curvatures of the surface, a coordinate invariant property, are computed in this mapped representation. From these curvatures, locally the bending energy is computed to describe the surface shape. An extensive evaluation shows that our curvature estimation is robust even in the presence of noise, independent of the scale of the object and furthermore the relative error stays small.