This paper describes a method for 3-D surface estimation based on reflection measurements in different orientationsof the light source and the paper. To facilitate the estimation, surface facets that probably support specular reflection weremarked and isolated, so as to enable the application of sample Lambertian reflectance model over those non-specular reflectionpoints. A computer controlled imaging device was used to capture images in a large number of orientations. The image capturewas performed using several different exposures in order to present larger dynamic irradiance range. Images together with theircounterpoints under antithetical illuminations were used to deal with the possible situation that light failed illumination arrivalson some points locate behind bulges.Iteration using difference approximation as well as enforcing integrability algorithm wereused to calculate the surface height based on Shaping from Shading algorithm.