The soil profile and its spatial distribution are two essential aspects for promoting sustainable agriculture, with precise inputs in quantity, space, and time. This work’s objective was to elaborate a digital map of soil fertility considering the complete profile for the accurate management of amendments and fertilizers. For the preparation of the soil fertility map, the following inputs were used: a digital elevation model, information from 44 soil profiles, the conversion of the properties of the soil profiles into surface units, geostatistical analysis of the soil properties, and the preparation of the final map with a geographic information system. The best spatial models were achieved with CEC, pH, Ca, Mg, Na, and K. The map of the soil fertility classes was produced considering CEC and the pH value. The soil fertility classes presented the following sequence of occupied surfaces: very low, medium, very high, high, and low. A process was generated to elaborate digital maps (geostatistics) of soil fertility using taxonomic information from soil profiles and considering the complete profile. The process converts soil classification into geographic and soil fertility information from basic science to application.