. Data processing
The relationships between the individual elements and between some soil parameters and the elements have been
assessed. Maps of isolevels of heavy metal contents have been produced as a first approach to know their spatial
distribution in the studied soils. Geographical information systems are increasingly being used for producing
thematic maps [18-20]. A raster-based GIS has been used to integrate physiographic, lithologic and soil information
and to produce a number of combinations of output data that resulted in several thematic maps. Thus, the spatial
distribution of heavy metals in the soils in and around the industrial complex has been mapped by geostatistical
interpolation. Two-dimensional ordinary block kriging as a most advantageous interpolation technique was applied
to produce regular grids.