This study aimed to investigate the spatial variability of soil characteristics, digital map and assess the best method of interpolation was performed. In this study, 62 points sampled from topsoil Karun East region were examined. These include the percentage of sand, silt, clay, EC, phosphorus, organic carbon, total nitrogen, calcium carbonate, cation exchange capacity, exchangeable sodium percentage is the sodium adsorption ratio. After the normalization of variables with different methods, first of all they draw the most appropriate theoretical model variograms fitted to experimental data gestated. Choose a suitable model for geostatistics studied characteristics, estimators kriging and cokriging in the GS⁺, was evaluated. Choose a suitable model based on the least error occurred MAE and RMSE. The draw Varyvgramhay kriging parameters of sand, silt, clay, calcium carbonate, cation exchange capacity, sodium absorption ratio of Gaussian, exponential organic carbon, total neutron, phosphorus spherical models were used. Kriging was used for other variables. Studies show that the Gaussian kriging model most used in the preparation of a digital map of soil variables. The silt highest (94,552 meters) and total nitrogen, the highest (2750 m) the impact in the region. Arc GIS software developed digital maps of soil variables and + GS represents the limestone in the area of the region. Soils and sodium and are deficient in organic matter.