In the work, the problem of presenting the spatial distribution of contamination with heavy metals in small-scale areas has been described. Two study areas have been chosen, located near industrial areas in south Poland – The “Miasteczko Śląskie” Zinc Smelter and the “Częstochowa” Steel Mill. A network of measurement points was planned, where from the top layer of soil, total of 108 samples have been taken. The content of three heavy metals was determined and included: Zn, Pb, Ba. Six methods of generating spatial models were selected: IDW, RBF, GPI, Ordinary Kriging, Simple Kriging and Universal Kriging. Two main criteria of evaluation were adopted: accuracy and visual quality of interpolation methods. In terms of accuracy analysis, no significant differences between the different methods have been observed. In the case of visual evaluation, it was found that the most suitable methods are IDW and RBF