Quantitative Mapping of Soil Property Based on Laboratory and Airborne Hyperspectral Data Using Machine Learning

Soil visible and near-infrared spectroscopy provides a non-destructive, rapid and low-cost approach to quantify various soil physical and chemical properties based on their reflectance in the spectral range of 400–2500 nm. With an increasing number of large-scale soil spectral libraries established across the world and new space-borne hyperspectral sensors, there is a need to explore methods to extract informative features from reflectance spectra and produce accurate soil spectroscopic models using machine learning. Features generated from regional or large-scale soil spectral data play a key role in the quantitative spectroscopic model for soil properties. The Land Use/Land Cover Area Frame Survey (LUCAS) soil library was used to explore PLS-derived components and fractal features generated from soil spectra in this study. The gradient-boosting method performed well when coupled with extracted features on the estimation of several soil properties. Transfer learning based on convolutional neural networks (CNNs) was proposed to make the model developed from laboratory data transferable for airborne hyperspectral data. The soil clay map was successfully derived using HyMap imagery and the fine-tuned CNN model developed from LUCAS mineral soils, as deep learning has the potential to learn transferable features that generalise from the source domain to target domain. The external environmental factors like the presence of vegetation restrain the application of imaging spectroscopy. The reflectance data can be transformed into a vegetation suppressed domain with a force invariance approach, the performance of which was evaluated in an agricultural area using CASI airborne hyperspectral data. However, the relationship between vegetation and acquired spectra is complicated, and more efforts should put on removing the effects of external factors to make the model transferable from one sensor to another.:Abstract I Kurzfassung III Table of Contents V List of Figures IX List of Tables XIII List of Abbreviations XV 1 Introduction 1 1.1 Motivation 1 1.2 Soil spectra from different platforms 2 1.3 Soil property quantification using spectral data 4 1.4 Feature representation of soil spectra 5 1.5 Objectives 6 1.6 Thesis structure 7 2 Combining Partial Least Squares and the Gradient-Boosting Method for Soil Property Retrieval Using Visible Near-Infrared Shortwave Infrared Spectra 9 2.1 Abstract 10 2.2 Introduction 10 2.3 Materials and methods 13 2.3.1 The LUCAS soil spectral library 13 2.3.2 Partial least squares algorithm 15 2.3.3 Gradient-Boosted Decision Trees 15 2.3.4 Calculation of relative variable importance 16 2.3.5 Assessment 17 2.4 Results 17 2.4.1 Overview of the spectral measurement 17 2.4.2 Results of PLS regression for the estimation of soil properties 19 2.4.3 Results of PLS-GBDT for the estimation of soil properties 21 2.4.4 Relative important variables derived from PLS regression and the gradient-boosting method 24 2.5 Discussion 28 2.5.1 Dimension reduction for high-dimensional soil spectra 28 2.5.2 GBDT for quantitative soil spectroscopic modelling 29 2.6 Conclusions 30 3 Quantitative Retrieval of Organic Soil Properties from Visible Near-Infrared Shortwave Infrared Spectroscopy Using Fractal-Based Feature Extraction 31 3.1 Abstract 32 3.2 Introduction 32 3.3 Materials and Methods 35 3.3.1 The LUCAS topsoil dataset 35 3.3.2 Fractal feature extraction method 37 3.3.3 Gradient-boosting regression model 37 3.3.4 Evaluation 41 3.4 Results 42 3.4.1 Fractal features for soil spectroscopy 42 3.4.2 Effects of different step and window size on extracted fractal features 45 3.4.3 Modelling soil properties with fractal features 47 3.4.3 Comparison with PLS regression 49 3.5 Discussion 51 3.5.1 The importance of fractal dimension for soil spectra 51 3.5.2 Modelling soil properties with fractal features 52 3.6 Conclusions 53 4 Transfer Learning for Soil Spectroscopy Based on Convolutional Neural Networks and Its Application in Soil Clay Content Mapping Using Hyperspectral Imagery 55 4.1 Abstract 55 4.2 Introduction 56 4.3 Materials and Methods 59 4.3.1 Datasets 59 4.3.2 Methods 62 4.3.3 Assessment 67 4.4 Results and Discussion 67 4.4.1 Interpretation of mineral and organic soils from LUCAS dataset 67 4.4.2 1D-CNN and spectral index for LUCAS soil clay content estimation 69 4.4.3 Application of transfer learning for soil clay content mapping using the pre-trained 1D-CNN model 72 4.4.4 Comparison between spectral index and transfer learning 74 4.4.5 Large-scale soil spectral library for digital soil mapping at the local scale using hyperspectral imagery 75 4.5 Conclusions 75 5 A Case Study of Forced Invariance Approach for Soil Salinity Estimation in Vegetation-Covered Terrain Using Airborne Hyperspectral Imagery 77 5.1 Abstract 78 5.2 Introduction 78 5.3 Materials and Methods 81 5.3.1 Study area of Zhangye Oasis 81 5.3.2 Data description 82 5.3.3 Methods 83 5.3.3 Model performance assessment 85 5.4 Results and Discussion 86 5.4.1 The correlation between NDVI and soil salinity 86 5.4.2 Vegetation suppression performance using the Forced Invariance Approach 86 5.4.3 Estimation of soil properties using airborne hyperspectral data 88 5.5 Conclusions 90 6 Conclusions and Outlook 93 Bibliography 97 Acknowledgements 11

Remote sensing approaches and mapping methods for monitoring soil salinity under different climate regimes

Soil salinization is one of the severe land-degradation problems due to its adverse effects on land productivity. Each year several hectares of lands are degraded due to primary or secondary soil salinization, and as a result, it is becoming a major economic and environmental concern in different countries. Spatio-temporal mapping of soil salinity is therefore important to support decisionmaking procedures for lessening adverse effects of land degradation due to the salinization. In that sense, satellite-based technologies provide cost effective, fast, qualitative and quantitative spatial information on saline soils. The main objective of this work is to highlight the recent remote sensing (RS) data and methods to assess soil salinity that is a worldwide problem. In addition, this study indicates potential linkages between salt-affected land and the prevailing climatic conditions of the case study areas being examined. Web of science engine is used for selecting relevant articles. "Soil salinity" is used as the main keyword for finding "articles" that are published from January 1, 2007 up to April 30, 2018. Then, 3 keywords; "remote sensing", "satellite" and "aerial" were used to filter the articles. After that, 100 case studies from 27 different countries were selected. Remote sensing based researches were further overviewed regarding to their location, spatial extent, climate regime, remotely sensed data type, mapping methods, sensing approaches together with the reason of salinity for each case study. In addition, soil salinity mapping methods were examined to present the development of different RS based methods with time. Studies are shown on the Köppen-Geiger climate classification map. Analysis of the map illustrates that 63% of the selected case study areas belong to arid and semi-arid regions. This finding corresponds to soil characteristics of arid regions that are more susceptible to salinization due to extreme temperature, high evaporation rates and low precipitation

A Case Study of the Forced Invariance Approach for Soil Salinity Estimation in Vegetation-Covered Terrain Using Airborne Hyperspectral Imagery

Soil spectroscopy is a promising technique for soil analysis, and has been successfully utilized in the laboratory. When it comes to space, the presence of vegetation significantly affects the performance of imaging spectroscopy or hyperspectral imaging on the retrieval of topsoil properties. The Forced Invariance Approach has been proven able to effectively suppress the vegetation contribution to the mixed image pixel. It takes advantage of scene statistics and requires no specific a priori knowledge of the referenced spectra. However, the approach is still mainly limited to lithological mapping. In this case study, the objective was to test the performance of the Forced Invariance Approach to improve the estimation accuracy of soil salinity for an agricultural area located in the semi-arid region of Northwest China using airborne hyperspectral data. The ground truth data was obtained from an eco-hydrological wireless sensing network. The relationship between Normalized Difference Vegetation Index (NDVI) and soil salinity is discussed. The results demonstrate that the Forced Invariance Approach is able to improve the retrieval accuracy of soil salinity at a depth of 10 cm, as indicated by a higher value for the coefficient of determination (R2). Consequently, the vegetation suppression method has the potential to improve quantitative estimation of soil properties with multivariate statistical methods