In this paper, we present a variational method that recovers both the shape and the reﬂectance of the Lambertian scene using multiple images. Although we consider only Lambertian surfaces in this paper, he proposed method, which is global and completely model based, is the ﬁrst and unavoidable stage for reaching a shape and reﬂectance estimation method for non-Lambertian surfaces. Basically, our method is a multiview stereo/shape from shading algorithm which allows to recover 3D shapes from Lambertian shading with known illumination conditions. Contrary to previous works that deal with a single material object of the constant albedo, our method works for surfaces with non-constant reﬂectance parameters, in particular with non-constant albedo. In addition, our algorithm is not based on two or more separate steps – shape and reﬂectance are jointly recovered in a same process. We veriﬁed the proposed method using synthetic images. We will extend our method for non-Lambertian surfaces to improve the robustness to non-Lambertian effects.