In forest management, it is of interest to obtain detailed inventories such that the local prediction errors on forest attributes are less than a prespecified threshold, while keeping the number of ground samples as low as possible. Given an initial sampling design, we propose an algorithm to determine the additional sample locations. The algorithm relies on two tools: geostatistical simulation, which allows measuring the uncertainty in the values of the attribute of interest, and simulated annealing, which allows finding an infill design that minimizes a given objective function. The proposed approach is applied to a data set from a Prosopis spp. plantation located in the Atacama Desert, in which the measured attribute is the rate of tree survival.Résumé : En gestion des forêts, il est intéressant d'obtenir des inventaires détaillés tels que les erreurs de prédiction locale sur les variables forestières soient inférieures à un seuil prédéfini, tout en réduisant le plus possible le nombre d'échan-tillons de terrain. É tant donné un plan d'échantillonnage initial, nous proposons un algorithme pour déterminer l'emplacement d'échantillons supplémentaires. L'algorithme repose sur deux outils : la simulation géostatistique, qui permet de mesurer l'incertitude sur les valeurs de la variable étudiée, et le recuit simulé, qui permet de trouver un plan d'échantil-lonnage minimisant une fonction objectif donnée. L'approche proposée est appliquée à un jeu de données d'une plantation de Prosopis située dans le désert d'Atacama, où la variable mesurée est le taux de survie des arbres.