Detailed soil surveys involve costly and time-consuming work and require expert knowledge. Since soil surveys provide information to meet a wide range of needs, new methods are necessary to map soils quickly and accurately. In this study, multilayer perceptron artificial neural networks (ANNs) were developed to map soil units using digital elevation model (DEM) attributes. Several optimal ANNs were produced based on a number of input data and hidden units. The approach used test and validation areas to calculate the accuracy of interpolated and extrapolated data. The results showed that the system and level of soil classification employed had a direct effect on the accuracy of the results. At the lowest level, smaller errors were observed with the World Reference Base (WRB) classification criteria than the Soil Taxonomy (ST) system, but more soil classes could be predicted when using ST (7 soils in the case of ST vs. 5 with WRB). Training errors were below 11% for all the ANN models applied, while the test error (interpolation error) and validation error (extrapolation error) were as high as 50% and 70%, respectively. As expected, soil prediction using a higher level of classification presented a better overall level of accuracy. To obtain better predictions, in addition to DEM attributes, data related to landforms and/or lithology as soil-forming factors, should be used as ANN input data.