In order for soil resources to be sustainably managed, it is necessary to have reliable, valid data on the spatial distribution of their environmental impact. However, in practice, one often has to cope with spatial interpolation achieved from few data that show a skewed distribution and uncertain information about soil contamination. We present a case study with 76 soil samples taken from a site of 15 square km in order to assess the usability of information gleaned from sparse data. The soil was contaminated with cadmium predominantly as a result of airborne emissions from a metal smelter. The spatial interpolation applies lognormal anisotropic kriging and conditional simulation for log-transformed data. The uncertainty of cadmium concentration acquired through data sampling, sample preparation, analytical measurement, and interpolation is factor 2 within 68.3 % confidence. Uncertainty predominantly results from the spatial interpolation necessitated by low sampling density and spatial heterogeneity. The interpolation data are shown in maps presenting likelihoods of exceeding threshold values as a result of a lognormal probability distribution. Although the results are not deterministic, this procedure yields a quantified and transparent estimation of the contamination, which can be used to delineate areas for soil improvement, remediation, or restricted area use, based on the decision-makers' probability safety requirement.