Hey there! In this tutorial I want to show and explain how to spatially estimate the distribution of clay content across an area in Python, using jupyter notebook and some packages such as geopandas and shap. Because the geopandas package is difficult to install, I created a new environment in the conda prompt and installed there the geopandas, following the instructions given in this link.
The follwing datasets were used:
dataset.csv
: This file contains four columnsX
andY
: Longitude and Latitude coordinates in EPSG: 4326 (WGS84).Clay.gkg
: Clay content in the soil in grams per kilogram.OM.gkg
: Organic matter content in the soil in grams per kilogram. This attribute is not used in predictions
SYSI.tif
: Bare soil image with six spectral bands used as predictors of the behavior of the clay content in the modeling process.