“…Compared with geostatistical and spatial interpolation methods (e.g., Kriging procedures, fuzzy clustering algorithms) based on comprehensive, time-consuming, and costly field surveys and soil sampling, non-invasive remote and proximal sensors, combined with empirical-or physical-based data analysis, offer potentially more effective, quick, and cost-efficient continuous direct or indirect data on physiochemical soil characteristics, which are determined by spatial, temporal or spectral sensor resolutions [1,6,[11][12][13][14][15][16][17][18]. Due to the advantages of existing, easily accessible data archives, relatively low-cost and high-temporal, high-spatial resolution multispectral imagery and time series are available for qualitative and partly-quantitative soil information extraction, deduction of soil patterns, and mapping of SSM zones and soil surface units [7,15].…”