3D shape reconstruction from images is an active topic in computer vision. Shape-from-Shading is an important approach which requires the surface properties and light source position to infer the 3D shape. A L2 regularizer is typically used to penalize the irradiance equation. In this article, anisotropic diffusion (AD) is introduced as a regularizer to solve the image irradiance equation. The method is then compared with L1 and L2 regularization methods, where all of the three techniques are formulated using gradient descent. Results shows that with AD, edges can be better preserved. AD shows lower depth error and higher correlation when compared with L1 and L2 regularization methods.