In this paper we study a modern, perspective model for shape from shading and its numerical approximation. We show that a new form of the classic concave/convex ambiguity is still present, although the model has been shown to be well-posed under particular assumptions. This analytical result is confirmed by various numerical tests. Moreover, we present convergence results for two iterative approximation schemes recently introduced in the literature. The first one is based on a finite difference discretization, whereas the second one is based on a semi-Lagrangian discretization. The convergence results are obtained in the general framework of viscosity solutions of the underlying partial differential equation. We also show that it is possible to obtain even in complex scenes results of reasonable quality. To this end we solve the constituting equation on a previously-segmented input image, where we use state constraints boundary conditions at the segment borders.