There is a consensus about the necessity to achieve a quick soil spatial information with few human resources. Remote/proximal sensing and pedotransference are methods that can be integrated into this approach. On the other hand, there is still a lack of strategies indicating on how to put this in practice, especially in the tropics. Thus, the objective of this work was to suggest a strategy for the spatial prediction of soil classes by using soil spectroscopy from ground laboratory spectra to space images platform, as associated with terrain attributes and spectral libraries. The study area is located in São Paulo State, Brazil, which was covered by a regular grid (one per ha), with 473 boreholes collected at top and undersurface. All soil samples were analyzed in laboratory (granulometry and chemical), and scanned with a VIS-NIR-SWIR (400-2500 nm) spectroradiometer. We developed two traditional pedological maps with different detail levels for comparison: TFS-1 regarding orders and subgroups; and TFS-2 with additional information such as color, iron and fertility. Afterwards, we performed a digital soil map, generated by models, which used the following information: (i) predicted soil attributes from undersurface layer (diagnostic horizon), obtained by using a local spectral library; (ii) spectral reflectance of a bare soil surface obtained by Landsat image; and (iii) derivative of terrain attributes. Thus, the digital map was generated by a combination of three variables: remote sensing (Landsat data), proximal sensing (laboratory spectroscopy) and relief. Landsat image with bare soil was used as a first observation of surface. This strategy assisted on the location of topossequences to achieve soil variation in the area. Afterwards, spectral undersurface information from these locations was used to modelize soil attributes quantification (156 samples). The model was used to quantify samples in the entire area by spectra (other 317 samples). Since the surface was bare soil, it was sampled by image spectroscopy. Indeed, topsoil spectral laboratory information presented great similarity with satellite spectra. We observed angle variation on spectra from clayey to sandy soils as differentiated by intensity. Soil lines between bands 3/4 and 5/7 were helpful on the link between laboratory and satellite data. The spectral models of soil attributes (i.e., clay, sand, and iron) presented a high predictive performance (R 2 0.71 to 0.90) with low error. The spatial prediction of these attributes also presented a high performance (validations with R 2 > 0.78). The models increased spatial resolution of soil weathering information using a known spectral library. Elevation (altitude) improved mapping due to correlation with soil attributes (i.e., clay, iron and chemistry). We observed a close relationship between soil weathering index map and laboratory spectra + image spectra + relief parameters. The color composite of the 5R, 4G and 3B had great performance on the detection of soils along topossequences, since...