iSDAsoil: soil pH for Africa predicted at 30 m resolution at 0-20 and 20-50 cm depths

<p>iSDAsoil dataset soil pH (1:1 Soil-Water Suspension) predicted at 30 m resolution for 0&ndash;20 and 20&ndash;50 cm depth intervals. Data has been projected in WGS84 coordinate system and compiled as <a href="https://gdal.org/drivers/raster/cog.html">COG</a>.&nbsp;Predictions have been generated using multi-scale Ensemble Machine Learning with 250 m (MODIS, PROBA-V, climatic variables and similar) and 30 m (DTM derivatives, Landsat, Sentinel-2 and similar) resolution covariates. For model training we use a pan-African compilations of soil samples and profiles (<a href="https://www.isda-africa.com/national-soil-services/">iSDA points</a>, <a href="https://www.isric.org/projects/africa-soil-profiles-database-afsp">AfSPDB</a>, and other national and regional soil datasets). Cite as:</p> <p>Hengl, T., Miller, M.A.E., Križan, J.&nbsp;<em>et al.</em>&nbsp;African soil properties and nutrients mapped at 30&nbsp;m spatial resolution using two-scale ensemble machine learning.&nbsp;<em>Sci Rep</em>&nbsp;<strong>11,&nbsp;</strong>6130 (2021). <a href="https://doi.org/10.1038/s41598-021-85639-y">https://doi.org/10.1038/s41598-021-85639-y</a></p> <p>To open the maps in QGIS and/or directly compute with them, please use the <a href="https://gitlab.com/openlandmap/africa-soil-and-agronomy-data-cube"><strong>Cloud-Optimized GeoTIFF version</strong></a>.</p> <p>Layer description:</p> <ul> <li>sol_ph_h2o_m_30m_*..*cm_2001..2017_v0.13_wgs84.tif = predicted soil pH mean value,</li> <li>sol_ph_h2o_md_30m_*..*cm_2001..2017_v0.13_wgs84.tif = predicted soil pH model (prediction) errors,</li> </ul> <p>Model errors were derived using bootstrapping: md is derived as standard deviation of individual learners from 5-fold cross-validation (using spatial blocking). The model 5-fold cross-validation (<a href="https://mlr.mlr-org.com/reference/makeStackedLearner.html">mlr::makeStackedLearner</a>) for this variable indicates:</p> <pre><code>Variable: ph_h2o R-square: 0.818 Fitted values sd: 0.972 RMSE: 0.459 Random forest model: Call: stats::lm(formula = f, data = d) Residuals: Min 1Q Median 3Q Max -5.5939 -0.2328 -0.0066 0.2222 4.7477 Coefficients: Estimate Std. Error t value Pr(&gt;|t|) (Intercept) 1.113440 1.164473 0.956 0.338986 regr.ranger 1.032918 0.003138 329.116 &lt; 2e-16 *** regr.xgboost -0.014201 0.004185 -3.393 0.000691 *** regr.cubist 0.049667 0.003709 13.392 &lt; 2e-16 *** regr.nnet -0.188570 0.188214 -1.002 0.316398 regr.cvglmnet -0.059763 0.003636 -16.438 &lt; 2e-16 *** --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 0.4591 on 133378 degrees of freedom Multiple R-squared: 0.8176, Adjusted R-squared: 0.8176 F-statistic: 1.195e+05 on 5 and 133378 DF, p-value: &lt; 2.2e-16</code></pre> <p>To back-transform values (y) to index use:</p> <pre><code>index = y / 10</code></pre> <p>To submit an issue or request support please visit <a href="https://isda-africa.com/isdasoil"><strong>https://isda-africa.com/isdasoil</strong></a></p>