Spatially explicit mapping depicting the distribution of an interest attribute is approached in a model-assisted framework using the inverse distance weighting (IDW) interpolator. The procedure also allows for the inclusion of wall-to-wall auxiliary variables. Moreover, a pseudo-population bootstrap procedure to estimate the precision of the resulting map is adopted. This methodology is applied to the Acquerino-Cantagallo Nature Reserve in central Italy to produce the wild boar rooting impact map and the corresponding precision map.