The invention relates to a method for estimating a set of variables associated with sites of an area, the method comprising the following steps: a) a set of observation data is supplied, b) a variable associated with a site of a first sub-area is estimated by krigeing, using the precision matrix of the first sub-area, c) a second adjacent sub-area is determined, d) a variable associated with a site of the second sub-area is estimated by krigeing, on the basis of the first sub-space, the second sub-space, and the precision matrix determined for the first sub-space.