In this report, we present a new method to recover an approximation of the bidirectional reflectance distribution function (BRDF) of the surfaces present in a real scene. This is done from a single photograph and a 3D geometric model of the scene. The result is a full model of the reflectance properties of all surfaces, which can be rendered under novel illumination conditions, with, for example, viewpoint modification and the addition of new synthetic objects. Our technique produces a reflectance model using a small number of parameters. Using a new hierarchical algorithm assumptions are done for each surface about their BRDF, which becomes more and more complex during the rerendering iterations. We treat each type of BRDF separately, applying a specific algorithm to take into account its own characteristics and its parameters. The number of parameters varies with the complexity of the model. When an approximation of the BRDF of all surfaces has been found, high quality rerenderings are produced using a global illumination software. We present several synthetic images that are compared to the original ones, as well as some possible applications in augmented reality.