The drawing of maps of soil quality at a large scale is increasingly being more useful to land planners and stakeholders. Nevertheless, it involves different methodological steps from the description of soil profiles in the field until the regional mapping of integrative soil quality index (IQI) values. The development of proper scoring functions is a paramount task for the calculation of these IQI values since every parameter needs to be standardized accordingly and weighting factors are usually estimated by multivariate techniques. The main goal of this study was to map soil quality in the Spanish region of Extremadura (commonly known by its rangelands called dehesas). To do that, i) we gathered information from 194 soil profiles described throughout the region, ii) we calculated the weighting factors of ten meaningful parameters used as indicators by using multivariate techniques (Principal Component Analysis, PCA; and Analytic Hierarchy Process, AHP), and iii) we developed standard scoring functions (SSFs) that represent the singularity of every variable (less is better, more is better). We established upper and lower limits for standardizing the values of each indicator properly. Regarding weighting factors, soil texture was highlighted by the PCA and nutrients by the AHP. Once IQI values were calculated, two regional maps of soil quality were drawn by using interpolation methods (ordinary kriging). The IQI maps showed remarkable spatial differences in soil quality presumably induced by land management. We conclude this methodology could be useful and we encourage other colleagues to test its effectiveness in places where soil data are available.