Mapping Particle Size and Soil Organic Matter in Tropical Soil Based on Hyperspectral Imaging and Non-Imaging Sensors

: We evaluated the use of airborne hyperspectral imaging and non-imaging sensors in the Vis—NIR—SWIR spectral region to assess particle size and soil organic matter in the surface layer of tropical soils (Oxisols, Ultisols, Entisols). The study area is near Piracicaba municipality, S ã o Paulo state, Brazil, in a sugarcane cultivation area of 135 hectares. The study area, with bare soil, was imaged in April 2016 by the AisaFENIX aerotransported hyperspectral sensor, with spectral resolution of 3.5 nm between 380 and 970 nm, and 12 nm between 970 and 2500 nm. We collected 66 surface soil samples. The samples were analyzed for particle size and soil organic matter content. Laboratory spectral measurements were performed using a non-imaging spectroradiometer (ASD FieldSpec 3 Jr). Partial Least Square Regression (PLSR) was used to predict clay, silt, sand and soil organic matter (SOM). The PLSR functions developed were applied to the hyperspectral image of the study area, allowing development of a prediction map of clay, sand, and SOM. The developed PLSR models demonstrated the relationship between the predictor variables at the cross-validation step, both for the non-imaging and imaging sensors, when the highest r and R 2 values were obtained for clay, sand, and SOM, with R 2 over 0.67. We did not obtain a satisfactory model for silt content. For the non-imaging sensor at the prediction step, R 2 values for clay and SOM were over 0.7 and sand was lower than 0.54. The imaging sensor yielded models for clay, sand, and SOM with R 2 values of 0.62, 0.66, and 0.67, respectively. Pearson correlation between sensors was greater than 0.849 for the prediction of clay, sand, and SOM. Our study successfully generated, from the imaging sensor, a large-scale and detailed predicted soil maps for particle size and SOM, which are important in the management of tropical soils. validation of prediction models. Even with the increase, these values can still suggest good performance for the content estimate of clay, sand and SOM [64]. These results demonstrate that Vis—NIR—SWIR spectroscopy is suitable for the multidimensional evaluation of soil attributes, just as reported by Paz-Kagan et al. and et al.


Introduction
Brazil has a great soil diversity due to a range of topography, biota, climate, and parent materials among the soil-forming factors. Agriculture is the main economic activity in Brazil and, therefore, the knowledge of pedogenetic processes and the capacity to differentiate soil classes are essential to improve soil management and productivity [1]. However, detailed, or even semi-detailed, soil use maps are rare in Brazil and, according to Mendonça-Santos [2] and Santos et al. [3], only 0.25-1% of available soil maps are plotted on an adequate scale. Furthermore, soil map production requires numerous on-theground observations and sample collection. The use of mathematical models obtained via Oxisols predominate on gentle, planar slopes and summits in the higher portions of the landscape. Ultisols occur on steeper slopes and Entisols occupy floodplains or low stream terraces ( Figure 2). Soil texture varies from sandy to clayey. Both the clayey and sandy soils are derived from parent materials from members of the Permian age Irati, Tatuí, and Itararé formations [27][28][29][30][31][32][33]. A total of 66 soil samples (0-20 cm depth- Figure 2) was uniformly distributed over the area and collected in July 2017. The soil samples were analyzed for particle size by the pipette method and wet sieving using 0.1 mol L −1 NaOH as a dispersing agent. Soil organic matter was obtained by determining the organic carbon content (Walkley-Black method) and subsequent multiplication by the 1.724 factor according to Teixeira et al. [34]. The The soil landscape is characterized by rolling hills with rounded summits and smooth slopes with lengths varying from 200 to 1500 m. Narrow valleys are cut into and separate the hillslopes. Elevation ranges between 532 to 570 m. Three soil orders Oxisols, Entisols, and Ultisols occur in the area.
Oxisols predominate on gentle, planar slopes and summits in the higher portions of the landscape. Ultisols occur on steeper slopes and Entisols occupy floodplains or low stream terraces ( Figure 2). Soil texture varies from sandy to clayey. Both the clayey and sandy soils are derived from parent materials from members of the Permian age Irati, Tatuí, and Itararé formations [27][28][29][30][31][32][33]. Oxisols predominate on gentle, planar slopes and summits in the higher portions of the landscape. Ultisols occur on steeper slopes and Entisols occupy floodplains or low stream terraces ( Figure 2). Soil texture varies from sandy to clayey. Both the clayey and sandy soils are derived from parent materials from members of the Permian age Irati, Tatuí, and Itararé formations [27][28][29][30][31][32][33]. A total of 66 soil samples (0-20 cm depth- Figure 2) was uniformly distributed over the area and collected in July 2017. The soil samples were analyzed for particle size by the pipette method and wet sieving using 0.1 mol L −1 NaOH as a dispersing agent. Soil organic matter was obtained by determining the organic carbon content (Walkley-Black method) and subsequent multiplication by the 1.724 factor according to Teixeira et al. [34]. The A total of 66 soil samples (0-20 cm depth- Figure 2) was uniformly distributed over the area and collected in July 2017. The soil samples were analyzed for particle size by the pipette method and wet sieving using 0.1 mol L −1 NaOH as a dispersing agent. Soil organic matter was obtained by determining the organic carbon content (Walkley-Black method) and subsequent multiplication by the 1.724 factor according to Teixeira et al. [34]. The descriptive analysis was carried out with the STATISTICA v. 12 software (StatSoft, Tulsa, OK, USA).

Non-Imaging Sensor Data Acquisition and Processing
For the spectral analysis in the laboratory, the samples were previously dried at 45 • C for 48 h, sieved in mesh with 2 mm to homogenize particles and to reduce humidity and sampling effects [35]. Afterward each sample was put on a Petri dish (with a diameter of 9 cm and 1.5 cm in height) and samples were subjected to spectral measurements. Reflectance curves were obtained by a ASD Fieldspec 3 Jr spectroradiometer (ASD, Boulder, CO, USA), with a spectral resolution of 3 nm between 350 and 1400 nm and 30 nm between 1400 and 2500 nm. Before obtaining spectral readings, the calibration in the non-imaging sensor (standard equipment) was performed using the Spectralon standard white plate as the maximum reflectance reference, according to Labsphere Reflectance Calibration Laboratory [36]. The optical fiber was attached to the pistol grip and placed on a support platform in a vertical position, at a distance of 8 cm from the samples. A 650 W halogen lamp with a non-collimated beam, located 35 cm from the platform and an angle of 30 • to the horizontal plane, was used as light source. The spectroradiometer was programmed to take 30 readings for each sample, thus leading to an average spectral curve, according to the Cezar et al. [13].
The spectral readings were carried out with three replicates for each sample while the Petri dish was successively replaced at 120 • clockwise, therefore analyzing the entire sample. Afterward a simple arithmetic mean of the three measurements was calculated, as described by Nanni and Demattê [37].

Imaging Sensor Data Acquisition and Processing
To compare the results of the data obtained by the spectroradiometer, considered standard equipment for obtaining spectral information when working with remote sensing [13,38,39], hyperspectral aerial images were obtained by the AisaFENIX hyperspectral image sensor, with the Vis-NIR-SWIR spectral range, and spectral resolution of 3.5 nm between 380 and 970 nm and 12 nm between 970 and 2500 nm.
A PA-34-200T Seneca II model twin-engined plane flew over the study area, on 30 April 2016, with the AisaFENIX imaging sensor and the GPS/IMU unit (OXTS-RT3050 model) installed on it, to collect the hyperspectral images and data related to the geometric correction of the image. The flight was performed at 14:35 p.m. (GMT) under open-sky conditions, fully sunny at an altitude of 660 m, resulting in a spatial resolution of about 1 m. Even though it is a sugar-growing region, the study area had bare soil, previously plowed and prepared for planting a new crop cycle.
The collected images by the sensor during the aerial assessment were pre-processed for geometric and radiometric correction through the CaliGeo PRO ® software [40]. The input data files were the raw data (.raw) generated by the AisaFENIX sensor, the sensor calibration file (.cal and the LUT), and the dark reference file (.hdr) to optimize spectral radiation data quality. This resulted in radiance and radiometrically corrected images.
After atmospheric correction, the images were transformed from radiance to reflectance, standardizing the differences of lighting and sensors. To perform image processing (atmospheric/topographic correction), the ATCOR-4 software, based in MODTRAM algorithm [41] was used. This computational packet receives the input data, which is the images in the radiance unit and converts it into a reflectance unit.
Then, from the AisaFENIX hyperspectral image, spectral curves were collected at the correspondent location of each soil sample through the Z-profile tool, software ENVI 5.1 [42]. The file was exported in txt format, then opened by the Excel software and the spectral curves were plotted, as well as the reflectance values of the points. Subsequently, the data was imported into The Unscrambler ® software (CAMO Software, Oslo, Norway) for statistical analysis.
Wavelengths between 1350 and 1400 nm and between 1800 and 2000 nm were excluded due to the interference in their reflectance values (greater absorption) caused by atmospheric water vapor [43]. In addition, the Savitzky-Golay filtering and smoothing algorithm (2nd order polynomial) [44] was applied to the data using the ENVI 5.1 software.

Multivariate Statistical Analysis and Mapping
The statistical analyses were carried out by both sensors. Initially, the marginal wavelengths (350-399 and 2361-2500 nm for the FieldSpec) and (380-467 and 2128-2500 nm for the AisaFENIX), have been removed due to noise [18,19].
For the non-imaging sensor, each spectral curve was submitted to baseline and light scattering correction by first and second derivatives (Savitzky-Golay), standard normal variate (SNV) and multiplicative scatter correction (MSC). Based on the best results, it was chosen the multiplicative scatter correction (MSC) method, according to Buddenbaum and Steffens [45]. This process was necessary to mathematically reduce the noise and to prevent the negative influence on the quality of the adjusted models [46].
The Unscrambler ® software (CAMO Software) was used to analyze major components to group samples and to identify possible atypical data for curves called outliers. Four samples that presented results not representative soils types of the area of interest and spectrally vastly different were excluded from the dataset to achieve better results not only in the calibration of the models but also in their validation.
Then prediction models of clay, silt, sand and organic matter contents were adjusted through partial least squares regression (PLSR) of data from AisaFENIX hyperspectral imaging sensor and from FieldSpec. To develop the statistical models, the 66 samples were randomly divided in two sets, independent of the soil orders. One set, composed of 46 samples (the calibration set), was used to calibrate the models and cross-validation testing, while the other set, composed of 20 samples (prediction set), was used to validate the calibrated models, according to Nanni et al. [47].
In the calibration step, matrices X and Y are modeled simultaneously to find the latent variables (PLS factors) in the wavelengths (X) that best predict the measured parameter (Y), thus, the data set is reduced to some factors, which are used in data prediction and construction of the statistical model [48]. The factor is characterized by the main changes observed in X, and, when this factor cannot explain any variation in the sample set, a second factor perpendicular to the first one will be used, and so on [49].
After data calibration, a cross-validation process was carried out, in which each sample used to build the model is tested. A leave-one-out cross-validation procedure was adopted to verify the prediction capability of the PLSR model for the calibration set. After the cross-validation step, the model was tested by assessing the ability to predict new samples from data that did not participate in the construction of the PLSR model.
The R 2 was used for prior assessment of calibration (c) models. RMSEc (calibration root mean square error) was the subsequent parameter for analyzing the model adjustment. However, the best estimate of the calibration model future performance is the RMSEcv (which is obtained during the cross-validation phase), producing more realistic and independent results for the data used in the construction of models [50]. Systematic error (BIAS) was also used as an evaluation parameter. Finally, the PLSR regression coefficients were used to study the significant wavelengths used in PLSR.
In order to obtain qualitative visual information, maps of the clay, sand and SOM content estimated by PLSR were also drawn up over the study area. For this purpose, the mathematical model generated by PLSR in the calibration step from the hyperspectral image, was inserted in the raster calculator of the QGIS software (QGIS Development Team), which from the original image, calculated the estimate of each attribute over the entire area of the study. This process was carried out individually for each attribute mentioned. Table 1 presents the results of the granulometric analysis and the SOM content of the 66 samples collected in the study area. The results confirm the contrasting aspect of the parent materials (siltstone and sandstone) observed in Figure 2. The measured data of particle-size of the samples were plotted in a USDA textural triangle ( Figure 3).   Table 2 presents the results from descriptive statistics of clay, silt, sand, and ganic matter attributes. The coefficients of variation (CV) showed high values for a uated attributes. According to Harr [51], the CV below 10% are considered "low", b 15% and 30% "moderate", and greater than 30% "high". Thus, the CVs were high fo clay and silt and moderate to sand. The high variability for particle size fractio plained by the influence of different parent materials on soils of the region, such a stone and siltstone [52], which presents varied mineralogical formation, consisting of Fe/Al oxides and Kaolinite, as reported by Vidal-Torrado and Lepsch [33].

Descriptive Analysis of Soil Attributes
For SOM, the high CV value is due to the high variation of the particle size f in the study area, since SOM is directly linked to the clay content for tropical soils observed in Table 3. Moreover, the agricultural use and management in the study also a determining factor [13], as there were places with and without the practic moving crop residues from sugarcane cultivation.   Table 2 presents the results from descriptive statistics of clay, silt, sand, and soil organic matter attributes. The coefficients of variation (CV) showed high values for all evaluated attributes. According to Harr [51], the CV below 10% are considered "low", between 15% and 30% "moderate", and greater than 30% "high". Thus, the CVs were high for SOM, clay and silt and moderate to sand. The high variability for particle size fraction is explained by the influence of different parent materials on soils of the region, such as sandstone and siltstone [52], which presents varied mineralogical formation, consisting mainly of Fe/Al oxides and Kaolinite, as reported by Vidal-Torrado and Lepsch [33].  Figure 4a,b demonstrate great variability in the reflectance of the evaluated soils, due to the difference in reflectance between soils developed from sandstone and siltstone formations, which can be seen in Figure 4c,d. In fact, this result was expected because sandy soils tend to have much higher reflectance due to their mineralogical formation (rich in quartz) and due to their low levels of SOM [54] and opaque minerals [55]. In the case of the curves obtained by the airborne imaging sensor (Figure 4b), a higher concentration is observed among the curves, with a lower amplitude of the reflectance factor. This fact is probably due to the natural moisture of the soil (not collected at the time of imaging) which ends up reducing the albedo of sandy and silty soils, making them closer to the clayey portions of the area. When the samples are dried in the laboratory and the standardized humidity is maintained, as in the case of the non-imager sensor (Figure 4a), there is a gap between the curves of the sandy spots (higher reflectance across spectral range) and the clayey ones (lower reflectance). In addition, for this sensor, concavities can be observed at 450 and 2200 nm due to iron oxides and kaolinite, respectively [55].  Figure 4a,b demonstrate great variability in the reflectance of the evaluated soils, due to the difference in reflectance between soils developed from sandstone and siltstone formations, which can be seen in Figure 4c,d. In fact, this result was expected because sandy soils tend to have much higher reflectance due to their mineralogical formation (rich in quartz) and due to their low levels of SOM [54] and opaque minerals [55]. In the case of the curves obtained by the airborne imaging sensor (Figure 4b), a higher concentration is observed among the curves, with a lower amplitude of the reflectance factor. This fact is probably due to the natural moisture of the soil (not collected at the time of imaging) which ends up reducing the albedo of sandy and silty soils, making them closer to the clayey portions of the area. When the samples are dried in the laboratory and the standardized humidity is maintained, as in the case of the non-imager sensor (Figure 4a), there is a gap between the curves of the sandy spots (higher reflectance across spectral range) and the clayey ones (lower reflectance). In addition, for this sensor, concavities can be observed at 450 and 2200 nm due to iron oxides and kaolinite, respectively [55].     Table 4 presents model parameters and performance indicators for soil attributes obtained through data modeling using the PLSR technique for both sensors. During the process for obtaining the PLSR models (Table 4), the relationship between the predictor variables (reflectance) and the predicted (attributes) were better explained for the nonimaging sensor by using the models containing: 10 factors for clay, five factors for silt, four factors for sand and eight factors for SOM. For the imaging sensor, the model presented five factors for clay, four factors for sand and three factors for SOM. For silt, it was not possible to establish a satisfactory model due to the low relation between reflectance obtained for this attribute, as demonstrated by Demattê et al. [56], Nanni and Demattê [57] and Nanni et al. [47].

Prediction Models Generated for Particle-Size Fraction and SOM from Spectral Data Obtained by ASD FieldSpec 3 Jr and AisaFENIX Sensors
It was observed that, for cross-validation, both in the imaging sensor and non imaging sensor, the values of R 2 varied from 0.70 to 0.81 for clay, from 0.69 to 0.72 for sand, and 0.67 to 0.80 for SOM (Table 4), which, according to Minasny and McBratney [58], are values considered satisfactory. Furthermore, these results corroborate the findings obtained by Paz-Kagan et al. [18], who studied the ability of an imager sensor to estimate attributes of an Israeli soil and found R 2 values for clay, sand, silt above 0.61. Nanni et al. [47], using a non-imaging sensor, found for clay and sand an R 2 of 0.89 and 0.92, respectively. But the authors found R 2 for SOM equal to 0.44, a very low value compared to the results obtained in the present study.
For silt, as expected, the cross-validation models presented low values possibly because of this fraction, in weathered soils, has many different minerals, which confused with clay, Fe and Al oxides, and that can make aggregates with organic matters [47,59]. Furthermore, the difficulty to estimate the silt content occurs because the fraction does not have linear relation with the soil reflectance, but a curvilinear one (hyperbolic), as described by Cezar et al. [60]. For both sensors, the bias of cross-validation for clay, sand and SOM was close to zero, indicating that there was a random point distribution on the regression line, showing no tendency of the values obtained [13]. The standard error (SE) close to the RMSE indicates the good fit of the obtained models [61,62].
The error in the prediction value (RMSEp), for both sensors, in most cases, had a small increase in relation to the cross-validation step (RMSEcv) for all studied attributes ( Figure 5). This is due to the reduction in precision between processes due to the decrease in the number of samples used in each step. Even so, the obtained results are satisfactory since no significant error variations were observed, indicating the good prediction capacity of the models [63]. However, the bias value was low, showing that there was no bias in the validation of prediction models. Even with the increase, these values can still suggest good performance for the content estimate of clay, sand and SOM [64]. These results demonstrate that Vis-NIR-SWIR spectroscopy is suitable for the multidimensional evaluation of soil attributes, just as reported by Paz-Kagan et al. [18] and Nanni et al. [47].
Remote Sens. 2021, 13, 1782 10 of 19 in the validation of prediction models. Even with the increase, these values can still suggest good performance for the content estimate of clay, sand and SOM [64]. These results demonstrate that Vis-NIR-SWIR spectroscopy is suitable for the multidimensional evaluation of soil attributes, just as reported by Paz-Kagan et al. [18] and Nanni et al. [47]. The statistical metrics show that the imaging sensor presented (in the prediction step), except for the silt variable, results similar to that of the non-imaging sensor. This is very promising, since the use of an imaging sensor allows a greater possibility of collecting information from an area, since with this type of sensor it is possible to obtain a spectral curve in each pixel of the image.
Furthermore, the generated models demonstrate the potential to be applied in areas with similar geological and edaphoclimatic characteristics, and further investigations are necessary to evaluate the applicability of this method in areas with different characteristics. During the prediction step (non-imaging sensor), values of r and R 2 for clay and SOM remained high (above 0.70). For clay, silt, sand and SOM, the results of R 2 found in the present work had different results in relation to those practiced by Martins et al. [65], For those same attributes, the results obtained were lower than those found by Paz-Kagan et al. [18], who, studying the use of image spectroscopy to evaluate the soil attributes, found R 2 values of 0.87, 0.73, 0.67 and 0.95 for clay, silt, sand and SOM determination models, respectively. Although the model has presented very consistent metrics, mainly for the silt attribute, the low number of samples evaluated to create the PLSR model (prediction set = 20 samples) may have limited the creation of models with higher R 2 p values, for example. The statistical metrics show that the imaging sensor presented (in the prediction step), except for the silt variable, results similar to that of the non-imaging sensor. This is very promising, since the use of an imaging sensor allows a greater possibility of collecting information from an area, since with this type of sensor it is possible to obtain a spectral curve in each pixel of the image.
Furthermore, the generated models demonstrate the potential to be applied in areas with similar geological and edaphoclimatic characteristics, and further investigations are necessary to evaluate the applicability of this method in areas with different characteristics.
During the prediction step (non-imaging sensor), values of r and R 2 for clay and SOM remained high (above 0.70). For clay, silt, sand and SOM, the results of R 2 found in the present work had different results in relation to those practiced by Martins et al. [65], For those same attributes, the results obtained were lower than those found by Paz-Kagan et al. [18], who, studying the use of image spectroscopy to evaluate the soil attributes, found R 2 values of 0.87, 0.73, 0.67 and 0.95 for clay, silt, sand and SOM determination models, respectively. Although the model has presented very consistent metrics, mainly for the silt attribute, the low number of samples evaluated to create the PLSR model (prediction set = 20 samples) may have limited the creation of models with higher R 2 p-values, for example.
According to Hair Jr et al. [66], for multivariate analysis, the recommended sample size is at least 20 samples for each group. The authors stated that larger number of samples generally reduces the sampling error and increases the sensitivity (power) of the test. The set of cross-validation (n = 46) and prediction (n = 20) of this work are within the aforementioned guideline (although at the limit) and the values of R 2 in general were satisfactory. However, it is recommended, for future work, to carry out tests with a larger sample size and to verify the real influence.
For the imaging sensor, in accordance to the error estimates, the models for clay, sand and SOM presented relevant capacity of prediction, with R 2 p of 0.62, 0.66 and 0.67 in the prediction respectively. These results were slightly lower than those obtained by Gomez et al. [67] (clay-R 2 p = 0.75), Selige et al. [68] (sand-R 2 p = 0.87) and Schwanghart and Jarmer [69] (SOM-R 2 p = 0.77).
The RMSEp, standard error (SEp) and bias p parameters, in all evaluated attributes in the prediction, presented low values. The observed low RMSEp values indicate that the model presents high prediction capacity [70]. For bias p , the low values indicate that the model is unbiased [71].
The results of the statistical parameters of the models (R 2 , RMSE, bias etc.), are similar between the two sensors. These results demonstrate that the use of the imaging sensor in a tropical environment for clay, sand and SOM determinations, can be as useful as the use of the non-imaging sensor (considered as a standard). In addition, the ability of the imaging sensor to carry out the simultaneous analysis of large areas is an important factor, which can greatly contribute to the generation of soil attribute maps, assisting the decision making of the production process, such as the choice of more suitable crops to the area in question, and even the rational use of correctives and fertilizers. Furthermore, these results are corroborated by other studies using image spectroscopy to determine the particle size and SOM, such as those by Hbirkou et al. [72] and Franceschini et al. [73].
By PLSR definition, all wavelengths participate, in different proportions or weights, in the construction of the model. In this regard, Figure 6 illustrates the regression coefficients for PLSR models of clay, silt, sand and SOM for both sensors, which are considered important to the prediction of these attributes through its determination coefficients (R 2 ).
For clay and sand, Vis (400 to 700 nm) and SWIR (2100 to 2400 nm) regions presented high values of regression coefficients for the both sensors ( Figure 6), indicating that these regions of the electromagnetic spectrum presents characteristics to predict these attributes. These results were expected for clay, since the Vis wavelengths (around to 600 nm) are mainly associated to minerals that have iron, e.g., hematite and goethite, which are frequent in Brazilian highly weathered soils [74][75][76], and the absorption in the SWIR range is due to soil clay minerals [23,77]. More specifically, in the SWIR region, at 2189 and 2260 nm, the highest values of regression coefficients were due to kaolinite and gibbsite, respectively, which are typical minerals of clay fraction in the tropical soils [78,79]. Around 1400 and 1900 nm, the highest values of regression coefficients were possibly obtained as a result from OH-groups of silicate clays, which are present in the samples [80]. For sand, at 1453 and 1910 nm, may be associated with minerals 2:1 covering over the quartz particles [79]. For silt attribute, the regression coefficients for the non-imaging sensor ( Figure 6) and imaging sensor do not explain the PLSR model due to low r and R 2 values (Table 4).
Pearson's correlation was performed between the estimated contents of clay, sand, and SOM by the imaging and non-imaging sensors through PLSR models, as shown in Figure 7. Significant positive correlations were observed, with an alpha significance level of 5% and p value lower than 0.0001, among the contents estimated of clay (0.858), sand (0.849), and SOM (0.897) by the PLSR models for the imaging and non-imaging sensors. This high correlation is due to the similarity between the curves, which provided similarity between the PLSR models of the sensors. These results are corroborated by the study of Franceschini et al. [73], who, evaluating the ability to predict SOM, clay and sand from an imaging and another non-imaging sensor, found similar relationships between these sensors.  The highest values of regression coefficients for SOM for the imaging sensor ( Figure 6) were found in the wavelengths between 400 and 600 nm, and between 1700 and 1900 nm, in agreement with Nawar et al. [81]. For the non-imaging sensor (SOM), observed high values of regression coefficients in wavelengths around 580, 1401, 1900-1940, 2180 and 2220 nm, agreeing with Chen et al. [8], Viscarra Rossel and Hicks [82] and Yang and Li [83]. These results were possibly found because of the interaction between organic matter and iron oxides (580 nm) and kaolinite (2220 nm) [84]. On the order hand, for both sensors, the high values of regression coefficients in wavelengths around 1401 and 1900-1940 nm come from water molecules vibrations and OH-associated with the SOM [52].
Pearson's correlation was performed between the estimated contents of clay, sand, and SOM by the imaging and non-imaging sensors through PLSR models, as shown in Figure 7. Significant positive correlations were observed, with an alpha significance level of 5% and p value lower than 0.0001, among the contents estimated of clay (0.858), sand (0.849), and SOM (0.897) by the PLSR models for the imaging and non-imaging sensors. This high correlation is due to the similarity between the curves, which provided similarity between the PLSR models of the sensors. These results are corroborated by the study of Franceschini et al. [73], who, evaluating the ability to predict SOM, clay and sand from an imaging and another non-imaging sensor, found similar relationships between these sensors.

Attributes Mapping Using PLSR Models
Pixel based spatial distribution maps for each soil attribute are shown in Figure 8 using the best cross validated regression models from the airborne imaging sensor.

Attributes Mapping Using PLSR Models
Pixel based spatial distribution maps for each soil attribute are shown in Figure 8 using the best cross validated regression models from the airborne imaging sensor. The mathematical model generated by PLSR in the calibration step for clay, sand and SOM were applied to all pixels from the image. The result was a quantification map for each pixel, divided into classes for better representation and visualization. As the soil at- The mathematical model generated by PLSR in the calibration step for clay, sand and SOM were applied to all pixels from the image. The result was a quantification map for each pixel, divided into classes for better representation and visualization. As the soil attributes presented a strong correlation with the image sensor data (Table 4), it was possible to quantify the estimates of interest. Similarly, Pullanagari et al. [85], mapped nutrients in pasture using classification algorithms.
It can be seen in Figure 8 that the highest levels of clay were found in the highest (and flat) portions of the landscape and the areas of the lowest (and steepest) portions have higher levels of sand. These results are in agreement with the results found in the samples analyzed in the laboratory and are also in agreement with Vidal-Torrado and Lepsch [33], working in the same area, who found clay contents in horizon A ranging from 16 to 58%. In addition, it is noted that it was possible to identify, through the map, areas with higher levels of clay on the terrace (marked with waypoints). These results were due to the exposure of the clayey horizon, typical of the Ultisols present in the area.
The highest levels of organic matter were found in portions of the landscape with the highest levels of clay. In fact, it is expected that the levels of organic matter are related to the texture of the soil. Lepsch [86] found a high correlation between the carbon and clay contents of the superficial horizon. Indeed, due to soil management for sugarcane, the soil is plowed and harrowed. With soil mobilization, the decomposition of SOM increases through its greater exposure to microbial attack. In portions with higher clay levels, the process is minimized, due to the higher humidity and lower temperature. On the other hand, as with more sandy portions, the opposite occurs. Lower water absorption increases the temperature of the soil and, therefore, greater decomposition of sound in these areas. In addition, higher levels of SOM can be observed within the terrace, possibly caused by the accumulation of this attribute through an erosive process.
The spatial distribution maps of content of soil attributes are information that becomes important for decision making at the time of planting, since physical attributes and organic matter are essential for a relevant crop. Moreover, this study contributes to understand the optimum analysis methods when mapping soil attributes at a regional, national or global scale using spaceborne imaging spectrometers.

Conclusions
The prediction of particle size and soil organic matter through hyperspectral imaging and non-imaging sensors yielded comparable results to traditional analytical techniques. It was possible to carry out the quantification of the clay (R 2 > 0.62), sand (R 2 > 0.54) and SOM (R 2 > 0.67) attributes through the use of data obtained by the hyperspectral imaging and non-imaging sensors.
Furthermore, highly significant positive correlations were observed among the estimated contents of clay (0.858), sand (0.849), and SOM (0.897) by the PLSR models for the imaging and non-imaging sensors.
Our study successfully generated large-scale and detailed soil prediction maps by the imaging sensor for particle size and SOM. Considering that Brazil is mainly occupied by large agricultural areas, the possibility of generating large-scale soil attributes prediction maps allows better knowledge about soil attributes (particle size and SOM) important to the management of tropical agricultural production.