WO2022038623A1 - System and method for remote quantification of electrical conductivity of soil - Google Patents

System and method for remote quantification of electrical conductivity of soil Download PDF

Info

Publication number
WO2022038623A1
WO2022038623A1 PCT/IN2020/051018 IN2020051018W WO2022038623A1 WO 2022038623 A1 WO2022038623 A1 WO 2022038623A1 IN 2020051018 W IN2020051018 W IN 2020051018W WO 2022038623 A1 WO2022038623 A1 WO 2022038623A1
Authority
WO
WIPO (PCT)
Prior art keywords
soil
module
saline
dielectric constant
electrical conductivity
Prior art date
Application number
PCT/IN2020/051018
Other languages
French (fr)
Inventor
Shoba Periasamy
Original Assignee
Shoba Periasamy
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shoba Periasamy filed Critical Shoba Periasamy
Publication of WO2022038623A1 publication Critical patent/WO2022038623A1/en

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9021SAR image post-processing techniques
    • G01S13/9027Pattern recognition for feature extraction
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/904SAR modes
    • G01S13/9076Polarimetric features in SAR
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/024Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using polarisation effects
    • G01S7/025Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using polarisation effects involving the transmission of linearly polarised waves

Definitions

  • the present invention is generally related to a system and method for remote identification of salination of soil.
  • the present invention is particularly related to a system and method for identifying salination of soil by remotely determining the electrical conductivity of soil in vegetated and bare lands.
  • the present invention is further related to a system and method for identifying salination of soil by remotely determining the imaginary part of the dielectric constant of the soil through C-band Synthetic-Aperture Radar (SAR).
  • SAR C-band Synthetic-Aperture Radar
  • the primary object of the present invention is to provide a system and method for identifying salination of soil by remotely determining the electrical conductivity (EC) of soil in vegetated lands.
  • EC electrical conductivity
  • Another object of the present invention is to simulate the SAR product to act as a proxy for soil textural characteristics.
  • Another object of the present invention is to provide a system and method for identifying salination of soil by remotely determining the imaginary part of the dielectric constant of the soil through C-band Synthetic-Aperture Radar (SAR).
  • SAR C-band Synthetic-Aperture Radar
  • Yet another object of the present invention is to enable retrieving the soil EC with optimum accuracy level even when the field samples are not equally distributed in all ranges of values.
  • Yet another object of the present invention is to quantify the soil Electrical Conductivity of the bare soils and the soils covered with crops of less than 0.5m height.
  • Yet another object of the present invention is to enable carrying out essential and appropriate agricultural processes according to actual salinity level in the soil.
  • the various embodiments of the present invention provide a system and method for identifying salination of soil by remotely determining the electrical conductivity (EC) of soil in vegetated lands.
  • the embodiments also provide a system and method for identifying salination of soil by remotely determining the imaginary part of the dielectric constant of the soil through C-band Synthetic-Aperture Radar (SAR).
  • SAR C-band Synthetic-Aperture Radar
  • a system for identifying salination of soil by remotely determining the electrical conductivity of soil in vegetated and bare lands.
  • the system comprises a Synthetic-Aperture Radar (SAR) electronics module, a Synthetic-Aperture Radar (SAR) antenna module and an application platform module.
  • the SAR electronics module further comprises a C-band signal generator, a processor module, a digitizer module, a Ground Range Detected (GRD) module and a communication module.
  • the SAR antenna module further comprises an antenna module, a transmitter module, an acquisition modes module, a receiver module and a waveguide radiators module.
  • the application platform module provides an architecture for processing and analysis for Earth Observation, and further comprises a pre-processing module and a calibrated backscattering module.
  • the calibrated backscattering module further comprises a soil moisture information database, soil textural index database, a Spatial Simulation module for dielectric constant of non-saline soils, a Spatial Simulation module for dielectric constant of saline soils and a Spatial simulation module of Density Space Dielectric Model for Soil Electrical Conductivity.
  • the communication module is configured to communicatively couple with a plurality of remote computing devices to determine the type of activities that are suitable on a piece of land based on preset data and observed and/or simulated data from the system.
  • the communication module is also configured to render the observations and simulated data on a plurality of computing devices in user-intelligible form.
  • a method for identifying salination of soil by remotely determining the electrical conductivity of soil in vegetated lands.
  • the method provides a Density Space Dielectric Model for Soil Electrical Conductivity (DSDM-SEC), wherein the DSDM-SEC is enabled by three-dimensional density space information obtained as Synthetic-Aperture Radar (SAR) data product in C-band frequency.
  • SAR Synthetic-Aperture Radar
  • the DSDM-SEC comprises the following steps: describing the available total soil moisture content ( ) by investigating the scattering and penetration behavior of SAR signals of VV and VH polarizations with an interacting medium in two-dimensional scatterplot; determining soil moisture content ( ) by decoupling the effect of vegetation moisture ( ) from the using DPSVI as a vegetation descriptor; providing an “Soil Textural Index (STI)” by mathematically investigating the water holding capacity of the soil in terms of available bound and free water in the soil; identifying models for dielectric constant of non-saline soil ( ) and saline soil ( ) from parameters such as field-observed dielectric constant measurements ( ), STI, and as three axial parameters in a three-dimensional density plot; providing an empirical Density Space Dielectric Model for Soil Electrical Conductivity and identifying Imaginary part of dielectric constant ( ) the soil by scaling down the values of from .
  • STI Soil Textural Index
  • the empirical model provided to estimate the soil moisture further comprises the following: constructing a theoretical scatterplot to the shape of a triangle in which the possible positions of a plurality of features such as Water bodies (WB), Dense Vegetation (D-Veg), Sparse Vegetation (S-Veg), Moisture Soil (MS) and Bare Soil without moisture (BS) are located according the de-polarization characteristics of the SAR signal, and wherein the scattering behavior in the scatter plot was analyzed from VH polarization (y-axis) and VV polarization (x-axis); constructing the functional two-dimensional scatterplot between VH and VV polarization to cross verify the theoretical assumptions of the possible positions of the different features; and, representing the soil moisture content as a simulation of diagonal distance measure (from BS to WB and D-Veg) obtained from approximated triangle in the two-dimensional scatter plot.
  • WB Water bodies
  • D-Veg Dense Vegetation
  • S-Veg Sparse Vegetation
  • MS Moisture Soil
  • the empirical model “Soil Textural Index (STI)” is provided for mapping and representing the soil textural characteristics as a relative measure from varying moisture content according to sand and clay fraction in the soil under semi-saturated condition.
  • identifying models for dielectric constant of non-saline soil ( ) and saline soil ( ) further includes: constructing a three-dimensional density space from inputs such as field-observed dielectric constant measurements ( ), STI, and as three axial parameters in which the z-axis ( ) is considered as a conductivity domain, while the x- and y-axis (STI, and ) are considered as a polarization domain; examining the relationship of STI with in representing dielectric constant in three-dimensional density space by keeping the values of the third dimension ( ) as null, and wherein the model of 2-Dimensional Spatial Distance Measure (2D-SDM) is provided by simulating the gradient starting from the pixels having maximum , and minimum STI values to the pixels distributed at the minimum range of , and maximum range of STI to indicate the dielectric constant of the soil according to the polarization property; providing an elevation to the gradient of 2-Dimensional Spatial Distance Measure (2D-SDM) in
  • an empirical model DSDM-SEC Density Space Dielectric Model for Soil Electrical Conductivity
  • DSDM-SEC Density Space Dielectric Model for Soil Electrical Conductivity
  • the model DSDM-SEC Density Space Dielectric Model for Soil Electrical Conductivity
  • DSDM-SEC Density Space Dielectric Model for Soil Electrical Conductivity
  • model DSDM-SEC Density Space Dielectric Model for Soil Electrical Conductivity
  • FIG. 1 illustrates a system for identifying salination of soil by remotely determining the electrical conductivity of soil in vegetated lands, according to one embodiment of the present invention.
  • FIG. 1 illustrates a calibrated backscattering module, according to one embodiment of the present invention.
  • FIG. 1 illustrates a graphical plot that enables the determination of soil moisture content, according to one embodiment of the present invention.
  • FIG. 1 illustrates a three-dimensional graphical plot that enables the determination of dielectric constant value of saline and non-saline soil, and dielectric loss according to one embodiment of the present invention.
  • the various embodiments of the present invention provide a system and method for identifying salination of soil by remotely determining the electrical conductivity (EC) of soil in vegetated and bare lands.
  • the embodiments also provide a system and method for identifying salination of soil by remotely determining the imaginary part of the dielectric constant of the soil through C-band Synthetic-Aperture Radar (SAR).
  • SAR C-band Synthetic-Aperture Radar
  • an empirical simulation for identifying soil textural index is provided by investigating the scattering and penetration behavior of C-band SAR with the interacting surface features in the two-dimensional scatter plot.
  • the spatial model to estimate the Electrical Conductivity of the soil from three-dimensional density space using C-band SAR has not been developed in any of the earlier inventions.
  • the three-dimensional density space is constructed with field-observed dielectric constant measurements of non-saline ( ) and saline soils ( ), STI, and as three axial parameters in which STI, and represent the polarization domain and and represent conductivity domain.
  • the 2D-SDM is given elevation in third-dimension by the angles of optimization to simulate 3D-SDM for saline and non-saline soil.
  • the spatial simulation for the imaginary part of dielectric constant ( ) is provided.
  • a system for identifying salination of soil by remotely determining the electrical conductivity (EC) of soil in vegetated and bare lands.
  • the system comprises a Synthetic-Aperture Radar (SAR) electronics module, a Synthetic-Aperture Radar (SAR) antenna module and an application platform module.
  • the SAR electronics module further comprises a C-band signal generator, a processor module, a digitizer module, a Ground Range Detected (GRD) module and a communication module.
  • the SAR antenna module further comprises an antenna module, a transmitter module, an acquisition modes module, a receiver module and a waveguide radiators module.
  • the application platform module provides an architecture for processing and analysis for Earth Observation, and further comprises a pre-processing module and a calibrated backscattering module.
  • the application platform is operated in C-band frequency (5.36 GHz).
  • the C-band signal is generated inside SAR electronic module and transmitted to the beam-forming network of SAR antenna module.
  • the signal transmission and reception over the earth surface is carried out by SAR antenna module in four different modes such as Strip map, Interferometric Wide swath, Extra-wide swath and Wave mode.
  • the signal transmission and reception of the SAR Antenna Module is controlled by waveguide radiators. i.e. the receiver module is turned off while transmission and transistor module is turn off while reception.
  • the received signals in SAR Antenna Module in terms of phase and amplitude are transmitted again to SAR Electronic module to process like digitization, data compression and formatting.
  • the process is resulted in the generation of various levels of data products among which the current invention concentrates on Ground Range Detected product (GRD) acquired in Interferometric Wide swath mode.
  • GRD Ground Range Detected product
  • an application platform module is provided.
  • the calibrated backscattering products of VV and VH polarization are retrieved from GRD product (delivered in terms of amplitude and phase) by the sequence of steps performed in the application platform module.
  • the steps include: pre-processing for radiometric calibration to retrieve sigma naught of VV and VH polarization; removing speckle from data products of VV and VH polarization; and, Range-Doppler terrain correction for VV and VH polarization to attain geographical orientation.
  • the pixels influenced by effective scatters, and corner reflectors were identified and masked by terrain masking function to reduce the uncertainty in the modeling results due to the incoherent backscattering values induced by terrain variations. The effect of double bounce scattering in the backscattering values is eliminated in this step.
  • an empirical model to arrive soil moisture content ( ) is provided by investigating the scattering characteristics of the SAR signal in two-dimensional scatter plot.
  • the empirical simulation is provided to retrieve “Soil Textural Index (STI)” by considering varying as a crucial indicator.
  • the semi-empirical simulations for (dielectric constant of non-saline soil from SAR data) and (dielectric constant of saline affected soil from SAR data) are invented by simulating the inputs such as field-observed dielectric constant measurements ( ), STI , and as three axial parameters in a three-dimensional density plot. Based on the promising partition observed between the trends of and in third dimension, the model is invented to arrive by scaling down the values of the product from the corresponding values of .
  • soil moisture content is identified.
  • the unique nature of the different surface features to influence the scattering and penetration characteristics of the SAR signal is examined by generating the scatterplot between and .
  • the theoretical scatterplot (Fig. 3a) is approximated to the shape of a triangle in which the possible positions of the different features like Water bodies (WB), Dense Vegetation (D-Veg), Sparse Vegetation (S-Veg), Moisture Soil (MS) and Bare Soil without moisture (BS) are located according to the de-polarization and penetration characteristics of an illuminated SAR signal.
  • the values of for the clay soil is expected to be slightly higher than sandy soil because under semi-saturated (53% moisture content) the water in the clay soil available in the form of bound water which allows the signal to attain penetration and less scattering return.
  • the pixels of bare soil surface would be expected to attain high values in VV polarization and low values in VH polarization because of which they should have been scattered at the lower right corner of the theoretical triangle (Fig. 3a).
  • the pixels representing dense vegetation are identified in the upper left portion of the approximated triangle (Periasamy, 2018) as the rate of depolarization of an illuminated signal is more significant in densely arranged vegetation and the sensitivity of the VV polarization towards its high moisture content.
  • the Diagonal Measure (DM) (represented in grey line in Fig. 3a & 3b) emanating from the point z to the edge xy in the scatterplot was considered as a first derivative to be modeled to derive the soil moisture content.
  • the point where the DM originates represents the bare soil (higher values of and lower values of ) and the edge where the DM ends indicates water bodies and moisture-rich dense vegetation (lower values of and higher values of ).
  • the highest value of DM indicates moisture-rich dense vegetation and water bodies, next to which the saturated sandy and clay soil were found.
  • the lower values of DM represent bare soils with less moisture or no moisture, and when the distance of DM increases towards the lower-left portion of the scatterplot, it starts showing the increment in the moisture content.
  • the DM generally gets saturated at its upper limit or pixels belong water bodies and dense vegetation.
  • the soil moisture content ( ) is estimated by decoupling the effect of vegetation moisture ( ) from the using DPSVI as a vegetation descriptor (Equation 2).
  • a method for identifying soil textural index is provided.
  • the empirical simulation is proposed to retrieve the target parameter “Soil Textural Index (STI )” by considering varying as a crucial indicator.
  • the simulation is done based on the assumption that the surface under investigation has the constant roughness value in the semi-saturated state, and the influence of the system parameter wavelength ( ) is constant and the incidence angle is controlled by adopting optimization procedure.
  • the significance of in representing the soil textural properties is mathematically analyzed and the predicted relationship is shown in equation 3.
  • the equation for STI was proposed in such a way to extract the normalized product ranging from 0 to 1 (Shafian and Maas, 2015) .
  • the stands for the value of soil textural index is soil moisture value in the i th pixel of the corresponding imagery, and stands for the maximum and minimum value of soil moisture.
  • Semi-empirical models are provided to estimate (dielectric constant of non-saline soil from SAR data) and (dielectric constant of saline affected soil from SAR data) by simulating the inputs such as field-observed dielectric constant measurements ( ), STI , and as three axial parameters in a three-dimensional density plot.
  • the inputs such as field-observed dielectric constant measurements ( ), STI , and as three axial parameters in a three-dimensional density plot.
  • the surface under investigation has been assumed to be smooth in a semi-saturated condition.
  • the relationship of STI with in representing dielectric constant has been examined in three-dimensional density space by keeping the values of the third dimension ( ) as null (Fig. 4a).
  • the model of 2-Dimensional Spatial Distance Measure ( 2D-SDM ) is invented by simulating the gradient starting from the pixels having maximum , and minimum STI values to the pixels distributed at the minimum range of , and maximum range of STI (Fig. 4a).
  • the equation for the gradient of 2D-SDM is developed in such a way that it increases with an increasing dielectric constant in the soil (Equation 4).
  • the two-dimensional spatial distance measure could be considered as a preliminary theoretical simulation to indicate the dielectric constant of the soil according to the relationship of polarization property of a signal with the grain-size distribution of the soil. (4)
  • the products of STI and make up the x -axis and y -axis, and the corresponding values of are taken up in the z -axis, respectively.
  • the z -axis is considered as a conductivity domain
  • the x - and y -axis are considered as a polarization domain.
  • the gradient of 2-dimensional spatial distance is given elevation in the z -axis ( ) by two different angles of optimization, namely optimization angle of saline soil ( ) and non-saline soil ( ).
  • the pixels representing non-saline soil are decoupled from saline soil and represented in a separate three-dimensional density space with the corresponding STI and values.
  • the outliers are excluded from the sampling points, and the shape of the three-dimensional density space for non-saline soil is approximated to a shape of a right angle triangle in which pixels oriented along the line indicating 3-Dimensional Spatial Distance Measure ( ) are considered as hypotenuse, are considered as an adjacent side, and are considered as an opposite side to . Because of the constant of proportionality found between the pixels distributed along , and there is no much variation observed in throughout its trend during the upliftment of the pixels from two-dimension to three-dimension. By analyzing the constant rate of the gradient from to , is computed by taking the inverse tangent between the mean of maximum dielectric constant ( readings recorded in the field for non-saline soils and the maximum computed value of the two-dimensional spatial measure ( ).
  • a model for (3-Dimensional Spatial Distance Measure for Non-Saline soil) is proposed by applying an inverse cosine rule from dividing the resultant product of by the cos of (Equ 5).
  • the model is proposed in such a way that it starts increasing the dielectric values of the non-saline pixels distributed in the SAR imagery as the three-dimensional spatial distance measure increases.
  • the same procedure is employed for the pixels recorded as saline affected to find out the , and to retrieve the (3-Dimensional Spatial Distance Measure for Saline soil) (Fig 4d) (Equ 6 & 6a).
  • the system comprises a Synthetic-Aperture Radar (SAR) electronics module 101, a Synthetic-Aperture Radar (SAR) antenna module 102 and an application platform module 103.
  • the SAR electronics module further comprises a C-band signal generator 101a, a processor module 101b, a digitizer module 101c, a Ground Range Detected (GRD) module 101d and a communication module 101e.
  • the SAR antenna module further comprises an antenna module 102a, a transmitter module 102b, an acquisition modes module102e, a receiver module 102c and a waveguide radiators module 102d.
  • the application platform module provides an architecture for processing and analysis for Earth Observation, and further comprises a pre-processing module 104 and a calibrated backscattering module 105.
  • the calibrated backscattering module further comprises a soil moisture information database 105a, soil textural index database 105b, a Spatial Simulation module 105c for dielectric constant of non-saline soils, a Spatial Simulation module 105d for dielectric constant of saline soils and a Spatial simulation module 105e of Density Space Dielectric Model for Soil Electrical Conductivity.
  • FIG. 5 illustrates an exemplary implementation of the present invention.
  • the results of SAR derived have shown optimum agreement with ground measurements (exclusive of water bodies) with a correlation coefficient equal to 0.86 and with the least RMSE (0.04) and Bias (0.005) values.
  • the various embodiments of the present invention provide a system and method for identifying salination of soil by remotely determining the electrical conductivity (EC) of soil in vegetated and bare lands.
  • the embodiments also provide a system and method for identifying salination of soil by remotely determining the imaginary part of the dielectric constant of the soil through C-band Synthetic-Aperture Radar (SAR).
  • SAR C-band Synthetic-Aperture Radar
  • the present invention yields better results irrespective of the variations in the distribution of the samples.
  • the invention is employed to retrieve the soil EC with optimum accuracy level even when the field samples are not equally distributed in all ranges of values.
  • the implementation of the present invention does not depend on the different types of field samples and the corresponding spatial data products, which makes it more convenient to employ in soil salinity studies.
  • the present invention enables quantifying the soil Electrical Conductivity of the bare soils and the soils covered with crops of less than 0.5m height.
  • the invention finds application in research groups working on agriculture related use-cases as the resultant product of the model is essential to carry out appropriate agricultural processes according to actual salinity level in the soil.

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The various embodiments of the present invention provide a system and method for identifying salination of soil by remotely determining the electrical conductivity (EC) of soil in vegetated and bare lands. The embodiments also provide a system and method for identifying salination of soil by remotely determining the imaginary part of the dielectric constant of the soil through C-band Synthetic-Aperture Radar (SAR). In the present invention, the three-dimensional density space is constructed with field-observed dielectric constant measurements of non-saline and saline soils. The present invention enables quantifying the soil Electrical Conductivity of the bare soils and the soils covered with crops of less than 0.5m height. The invention finds application in research groups working on agriculture related use-cases as the resultant product of the model is essential to carry out appropriate agricultural processes according to actual salinity level in the soil.

Description

SYSTEM AND METHOD FOR REMOTE QUANTIFICATION OF ELECTRICAL CONDUCTIVITY OF SOIL
The present invention is generally related to a system and method for remote identification of salination of soil. The present invention is particularly related to a system and method for identifying salination of soil by remotely determining the electrical conductivity of soil in vegetated and bare lands. The present invention is further related to a system and method for identifying salination of soil by remotely determining the imaginary part of the dielectric constant of the soil through C-band Synthetic-Aperture Radar (SAR).
BACKGROUND OF THE INVENTION
The accurate measurement of soil’s electrical conductivity (EC) over vegetated land is challenging in multispectral remote sensing due to the lack of surface penetration capacity. Synthetic Aperture Radar (SAR) imagery has great potential to investigate soil EC due to its penetration capacity and high sensitivity towards the dielectric properties of the medium. However, the theoretical and empirical models of SAR are site-specific and need too many in-situ observations to implement the model in the spatial domain. These practical limitations lead to the development of semi-empirical models over bare soils. Though these semi-empirical models were widely accepted and employed in various studies, some studies reported discrepancies between predicted and measured values.
Furthermore, the models have studied only the dependency of dielectric constant on soil moisture while ignoring the effect of soil EC on the dielectric constant. In Dobson model (1985), which is considered to be the most accurate model so far, the potential of the imaginary part was explored in describing the soil EC as experimental research work. However, the model does not address retrieving spatial products representing soil salinity at a finer resolution.
Hence, there exists a need for an effective system and method for remotely determining the electrical conductivity of soil in vegetated lands. There also exists a need for a method to identify the electrical conductivity of soil from the imaginary part of dielectric constant through remote methodologies such as a C-band Synthetic-Aperture Radar (SAR) with high prediction accuracy level for all ranges.
The above-mentioned shortcomings, disadvantages and problems are addressed herein and which will be understood by reading and studying the following specification.
OBJECT OF THE INVENTION
The primary object of the present invention is to provide a system and method for identifying salination of soil by remotely determining the electrical conductivity (EC) of soil in vegetated lands.
Another object of the present invention is to simulate the SAR product to act as a proxy for soil textural characteristics.
Another object of the present invention is to provide a system and method for identifying salination of soil by remotely determining the imaginary part of the dielectric constant of the soil through C-band Synthetic-Aperture Radar (SAR).
Yet another object of the present invention is to enable retrieving the soil EC with optimum accuracy level even when the field samples are not equally distributed in all ranges of values.
Yet another object of the present invention is to quantify the soil Electrical Conductivity of the bare soils and the soils covered with crops of less than 0.5m height.
Yet another object of the present invention is to enable carrying out essential and appropriate agricultural processes according to actual salinity level in the soil.
These and other objects and advantages of the present invention will become readily apparent from the following detailed description taken in conjunction with the accompanying drawings.
SUMMARY OF THE INVENTION
The various embodiments of the present invention provide a system and method for identifying salination of soil by remotely determining the electrical conductivity (EC) of soil in vegetated lands. The embodiments also provide a system and method for identifying salination of soil by remotely determining the imaginary part of the dielectric constant of the soil through C-band Synthetic-Aperture Radar (SAR).
According to one embodiment of the present invention, a system is provided for identifying salination of soil by remotely determining the electrical conductivity of soil in vegetated and bare lands. The system comprises a Synthetic-Aperture Radar (SAR) electronics module, a Synthetic-Aperture Radar (SAR) antenna module and an application platform module. The SAR electronics module further comprises a C-band signal generator, a processor module, a digitizer module, a Ground Range Detected (GRD) module and a communication module. The SAR antenna module further comprises an antenna module, a transmitter module, an acquisition modes module, a receiver module and a waveguide radiators module. The application platform module provides an architecture for processing and analysis for Earth Observation, and further comprises a pre-processing module and a calibrated backscattering module.
According to one embodiment of the present invention, the calibrated backscattering module further comprises a soil moisture information database, soil textural index database, a Spatial Simulation module for dielectric constant of non-saline soils, a Spatial Simulation module for dielectric constant of saline soils and a Spatial simulation module of Density Space Dielectric Model for Soil Electrical Conductivity.
According to one embodiment of the present invention, the communication module is configured to communicatively couple with a plurality of remote computing devices to determine the type of activities that are suitable on a piece of land based on preset data and observed and/or simulated data from the system. The communication module is also configured to render the observations and simulated data on a plurality of computing devices in user-intelligible form.
According to one embodiment of the present invention, a method is provided for identifying salination of soil by remotely determining the electrical conductivity of soil in vegetated lands. The method provides a Density Space Dielectric Model for Soil Electrical Conductivity (DSDM-SEC), wherein the DSDM-SEC is enabled by three-dimensional density space information obtained as Synthetic-Aperture Radar (SAR) data product in C-band frequency. The DSDM-SEC comprises the following steps: describing the available total soil moisture content (
Figure pctxmlib-appb-M000001
) by investigating the scattering and penetration behavior of SAR signals of VV and VH polarizations with an interacting medium in two-dimensional scatterplot; determining soil moisture content (
Figure pctxmlib-appb-M000002
) by decoupling the effect of vegetation moisture (
Figure pctxmlib-appb-M000003
) from the
Figure pctxmlib-appb-M000004
using DPSVI as a vegetation descriptor; providing an “Soil Textural Index (STI)” by mathematically investigating the water holding capacity of the soil in terms of available bound and free water in the soil; identifying models for dielectric constant of non-saline soil (
Figure pctxmlib-appb-M000005
) and saline soil (
Figure pctxmlib-appb-M000006
) from parameters such as field-observed dielectric constant measurements (
Figure pctxmlib-appb-M000007
), STI, and
Figure pctxmlib-appb-M000008
as three axial parameters in a three-dimensional density plot; providing an empirical Density Space Dielectric Model for Soil Electrical Conductivity and identifying Imaginary part of dielectric constant (
Figure pctxmlib-appb-M000009
) the soil by scaling down the values of
Figure pctxmlib-appb-M000010
from
Figure pctxmlib-appb-M000011
.
According to one embodiment of the present invention, the empirical model provided to estimate the soil moisture further comprises the following: constructing a theoretical scatterplot to the shape of a triangle in which the possible positions of a plurality of features such as Water bodies (WB), Dense Vegetation (D-Veg), Sparse Vegetation (S-Veg), Moisture Soil (MS) and Bare Soil without moisture (BS) are located according the de-polarization characteristics of the SAR signal, and wherein the scattering behavior in the scatter plot was analyzed from VH polarization (y-axis) and VV polarization (x-axis); constructing the functional two-dimensional scatterplot between VH and VV polarization to cross verify the theoretical assumptions of the possible positions of the different features; and, representing the soil moisture content as a simulation of diagonal distance measure (from BS to WB and D-Veg) obtained from approximated triangle in the two-dimensional scatter plot.
According to one embodiment of the present invention, the empirical model “Soil Textural Index (STI)” is provided for mapping and representing the soil textural characteristics as a relative measure from varying moisture content according to sand and clay fraction in the soil under semi-saturated condition.
According to one embodiment of the present invention, identifying models for dielectric constant of non-saline soil (
Figure pctxmlib-appb-M000012
) and saline soil (
Figure pctxmlib-appb-M000013
) further includes: constructing a three-dimensional density space from inputs such as field-observed dielectric constant measurements (
Figure pctxmlib-appb-M000014
), STI, and
Figure pctxmlib-appb-M000015
as three axial parameters in which the z-axis (
Figure pctxmlib-appb-M000016
) is considered as a conductivity domain, while the x- and y-axis (STI, and
Figure pctxmlib-appb-M000017
) are considered as a polarization domain; examining the relationship of STI with
Figure pctxmlib-appb-M000018
in representing dielectric constant in three-dimensional density space by keeping the values of the third dimension (
Figure pctxmlib-appb-M000019
) as null, and wherein the model of 2-Dimensional Spatial Distance Measure (2D-SDM) is provided by simulating the gradient starting from the pixels having maximum
Figure pctxmlib-appb-M000020
, and minimum STI values to the pixels distributed at the minimum range of
Figure pctxmlib-appb-M000021
, and maximum range of STI to indicate the dielectric constant of the soil according to the polarization property; providing an elevation to the gradient of 2-Dimensional Spatial Distance Measure (2D-SDM) in the z-axis (
Figure pctxmlib-appb-M000022
) by two different angles of optimization, namely optimization angle of saline soil (
Figure pctxmlib-appb-M000023
) and non-saline soil (
Figure pctxmlib-appb-M000024
) based on the mean values of maximum 2D-SDM and
Figure pctxmlib-appb-M000025
of saline (
Figure pctxmlib-appb-M000026
) and non-saline soils (
Figure pctxmlib-appb-M000027
); and, identifying the 3D-SDM(NS) (3-Dimensional Spatial Distance Measure for Non-Saline soil) or
Figure pctxmlib-appb-M000028
and 3D-SDM(S) (3-Dimensional Spatial Distance Measure for Saline soil) or
Figure pctxmlib-appb-M000029
by applying an inverse cosine rule from dividing the resultant product of 2D-SDM by the cos of
Figure pctxmlib-appb-M000030
and
Figure pctxmlib-appb-M000031
to represent dielectric constant of non-saline soil (
Figure pctxmlib-appb-M000032
) and saline soil (
Figure pctxmlib-appb-M000033
).
According to one embodiment of the present invention, an empirical model DSDM-SEC (Density Space Dielectric Model for Soil Electrical Conductivity) is provided by subtracting the resultant products of dielectric constant of saline soil (
Figure pctxmlib-appb-M000034
) from the corresponding product of non-saline soil (
Figure pctxmlib-appb-M000035
).
According to one embodiment of the present invention, the model DSDM-SEC (Density Space Dielectric Model for Soil Electrical Conductivity) is dependent on soil textural characteristics under the semi-saturated condition.
According to one embodiment of the present invention, wherein the model DSDM-SEC (Density Space Dielectric Model for Soil Electrical Conductivity) is potential to quantify the soil Electrical Conductivity of the bare soils and soils covered with the crops of less than 0.5m height.
These and other aspects of the embodiments herein will be better appreciated and understood when considered in conjunction with the following description and the accompanying drawings. It should be understood, however, that the following descriptions, while indicating the preferred embodiments and numerous specific details thereof, are given by way of illustration and not of limitation. Many changes and modifications may be made within the scope of the embodiments herein without departing from the spirit thereof, and the embodiments herein include all such modifications.
BRIEF DESCRIPTION OF THE DRAWINGS
The other objects, features and advantages will occur to those skilled in the art from the following description of the preferred embodiment and the accompanying drawings in which:
illustrates a system for identifying salination of soil by remotely determining the electrical conductivity of soil in vegetated lands, according to one embodiment of the present invention.
illustrates a calibrated backscattering module, according to one embodiment of the present invention.
illustrates a graphical plot that enables the determination of soil moisture content, according to one embodiment of the present invention.
illustrates a three-dimensional graphical plot that enables the determination of dielectric constant value of saline and non-saline soil, and dielectric loss according to one embodiment of the present invention.
illustrates an exemplary implementation of the present invention, according to one embodiment of the present invention.
Although the specific features of the present invention are shown in some drawings and not in others. This is done for convenience only as each feature may be combined with any or all of the other features in accordance with the present invention.
DETAILED DESCRIPTION OF THE INVENTION
In the following detailed description, a reference is made to the accompanying drawings that form a part hereof, and in which the specific embodiments that may be practiced is shown by way of illustration. These embodiments are described in sufficient detail to enable those skilled in the art to practice the embodiments and it is to be understood that other changes may be made without departing from the scope of the embodiments. The following detailed description is therefore not to be taken in a limiting sense.
The various embodiments of the present invention provide a system and method for identifying salination of soil by remotely determining the electrical conductivity (EC) of soil in vegetated and bare lands. The embodiments also provide a system and method for identifying salination of soil by remotely determining the imaginary part of the dielectric constant of the soil through C-band Synthetic-Aperture Radar (SAR).
According to one embodiment of the present invention, an empirical simulation for identifying soil textural index (STI) is provided by investigating the scattering and penetration behavior of C-band SAR with the interacting surface features in the two-dimensional scatter plot. The spatial model to estimate the Electrical Conductivity of the soil from three-dimensional density space using C-band SAR has not been developed in any of the earlier inventions. In the present invention, the three-dimensional density space is constructed with field-observed dielectric constant measurements of non-saline (
Figure pctxmlib-appb-M000036
) and saline soils (
Figure pctxmlib-appb-M000037
), STI, and
Figure pctxmlib-appb-M000038
as three axial parameters in which STI, and
Figure pctxmlib-appb-M000039
represent the polarization domain and
Figure pctxmlib-appb-M000040
and
Figure pctxmlib-appb-M000041
represent conductivity domain. The 2D-SDM is given elevation in third-dimension by the angles of optimization to simulate 3D-SDM for saline and non-saline soil. Based on the promising partition found between the 3D-SDM (NS) or (
Figure pctxmlib-appb-M000042
) and 3D-SDM (S) or (
Figure pctxmlib-appb-M000043
) in the conductivity domain, the spatial simulation for the imaginary part of dielectric constant (
Figure pctxmlib-appb-M000044
) is provided.
According to one embodiment of the present invention, a system is provided for identifying salination of soil by remotely determining the electrical conductivity (EC) of soil in vegetated and bare lands. The system comprises a Synthetic-Aperture Radar (SAR) electronics module, a Synthetic-Aperture Radar (SAR) antenna module and an application platform module. The SAR electronics module further comprises a C-band signal generator, a processor module, a digitizer module, a Ground Range Detected (GRD) module and a communication module. The SAR antenna module further comprises an antenna module, a transmitter module, an acquisition modes module, a receiver module and a waveguide radiators module. The application platform module provides an architecture for processing and analysis for Earth Observation, and further comprises a pre-processing module and a calibrated backscattering module. The application platform is operated in C-band frequency (5.36 GHz). The C-band signal is generated inside SAR electronic module and transmitted to the beam-forming network of SAR antenna module. The signal transmission and reception over the earth surface is carried out by SAR antenna module in four different modes such as Strip map, Interferometric Wide swath, Extra-wide swath and Wave mode. The signal transmission and reception of the SAR Antenna Module is controlled by waveguide radiators. i.e. the receiver module is turned off while transmission and transistor module is turn off while reception. The received signals in SAR Antenna Module in terms of phase and amplitude are transmitted again to SAR Electronic module to process like digitization, data compression and formatting. The process is resulted in the generation of various levels of data products among which the current invention concentrates on Ground Range Detected product (GRD) acquired in Interferometric Wide swath mode.
According to one embodiment of the present invention, an application platform module is provided. The calibrated backscattering products of VV and VH polarization are retrieved from GRD product (delivered in terms of amplitude and phase) by the sequence of steps performed in the application platform module. The steps include: pre-processing for radiometric calibration to retrieve sigma naught of VV and VH polarization; removing speckle from data products of VV and VH polarization; and, Range-Doppler terrain correction for VV and VH polarization to attain geographical orientation. The pixels influenced by effective scatters, and corner reflectors were identified and masked by terrain masking function to reduce the uncertainty in the modeling results due to the incoherent backscattering values induced by terrain variations. The effect of double bounce scattering in the backscattering values is eliminated in this step.
According to one embodiment of the present invention, an empirical model to arrive soil moisture content (
Figure pctxmlib-appb-M000045
) is provided by investigating the scattering characteristics of the SAR signal in two-dimensional scatter plot. The empirical simulation is provided to retrieve “Soil Textural Index (STI)” by considering varying
Figure pctxmlib-appb-M000046
as a crucial indicator. The semi-empirical simulations for
Figure pctxmlib-appb-M000047
(dielectric constant of non-saline soil from SAR data) and
Figure pctxmlib-appb-M000048
(dielectric constant of saline affected soil from SAR data) are invented by simulating the inputs such as field-observed dielectric constant measurements (
Figure pctxmlib-appb-M000049
), STI, and
Figure pctxmlib-appb-M000050
as three axial parameters in a three-dimensional density plot. Based on the promising partition observed between the trends of
Figure pctxmlib-appb-M000051
and
Figure pctxmlib-appb-M000052
in third dimension, the model
Figure pctxmlib-appb-M000053
is invented to arrive
Figure pctxmlib-appb-M000054
by scaling down the values of the product
Figure pctxmlib-appb-M000055
from the corresponding values of
Figure pctxmlib-appb-M000056
.
According to one embodiment of the present invention, soil moisture content is identified. The
Figure pctxmlib-appb-M000057
and
Figure pctxmlib-appb-M000058
stands for the optimized backscattering products of VV and VH polarization for incidence angle variations. The unique nature of the different surface features to influence the scattering and penetration characteristics of the SAR signal is examined by generating the scatterplot between
Figure pctxmlib-appb-M000059
and
Figure pctxmlib-appb-M000060
. The theoretical scatterplot (Fig. 3a) is approximated to the shape of a triangle in which the possible positions of the different features like Water bodies (WB), Dense Vegetation (D-Veg), Sparse Vegetation (S-Veg), Moisture Soil (MS) and Bare Soil without moisture (BS) are located according to the de-polarization and penetration characteristics of an illuminated SAR signal.
The pixels representing water bodies attain lower values in both the polarizations as the water surfaces induce additional forward scattering and thus reduces the backscattering return of the SAR signals. The ability of sandy soil to allow free water movement in it, often resulted in the semi-forward scattering, because of which the pixels representing semi-saturated sandy soil attains low
Figure pctxmlib-appb-M000061
values next to water bodies. The values of
Figure pctxmlib-appb-M000062
for the clay soil is expected to be slightly higher than sandy soil because under semi-saturated (53% moisture content) the water in the clay soil available in the form of bound water which allows the signal to attain penetration and less scattering return. The electromagnetic wave experiences scattering while interacting with bound water.
Since the depolarization rate of loose and dry-bare soil is very minimal, the pixels of bare soil surface would be expected to attain high values in VV polarization and low values in VH polarization because of which they should have been scattered at the lower right corner of the theoretical triangle (Fig. 3a). The pixels representing dense vegetation are identified in the upper left portion of the approximated triangle (Periasamy, 2018) as the rate of depolarization of an illuminated signal is more significant in densely arranged vegetation and the sensitivity of the VV polarization towards its high moisture content.
The Diagonal Measure (DM) (represented in grey line in Fig. 3a & 3b) emanating from the point z to the edge xy in the scatterplot was considered as a first derivative to be modeled to derive the soil moisture content. The point where the DM originates represents the bare soil (higher values of
Figure pctxmlib-appb-M000063
and lower values of
Figure pctxmlib-appb-M000064
) and the edge where the DM ends indicates water bodies and moisture-rich dense vegetation (lower values of
Figure pctxmlib-appb-M000065
and higher values of
Figure pctxmlib-appb-M000066
). Hence the highest value of DM indicates moisture-rich dense vegetation and water bodies, next to which the saturated sandy and clay soil were found. The lower values of DM represent bare soils with less moisture or no moisture, and when the distance of DM increases towards the lower-left portion of the scatterplot, it starts showing the increment in the moisture content. The DM generally gets saturated at its upper limit or pixels belong water bodies and dense vegetation. An empirical model was developed for DM in the scatterplot based on the Pythagoras theorem (Equation 1) to exemplify the total moisture content (
Figure pctxmlib-appb-M000067
).
Figure pctxmlib-appb-M000068
(1)
Figure pctxmlib-appb-M000069
(2)
where
Figure pctxmlib-appb-M000070
=
Figure pctxmlib-appb-M000071
* DPSVI
The soil moisture content (
Figure pctxmlib-appb-M000072
) is estimated by decoupling the effect of vegetation moisture (
Figure pctxmlib-appb-M000073
) from the
Figure pctxmlib-appb-M000074
using DPSVI as a vegetation descriptor (Equation 2).
According to one embodiment of the present invention, a method for identifying soil textural index (STI) is provided. The empirical simulation is proposed to retrieve the target parameter “Soil Textural Index ( STI)” by considering varying
Figure pctxmlib-appb-M000075
as a crucial indicator. The simulation is done based on the assumption that the surface under investigation has the constant roughness value in the semi-saturated state, and the influence of the system parameter wavelength (
Figure pctxmlib-appb-M000076
) is constant and the incidence angle is controlled by adopting optimization procedure. Hence, the significance of
Figure pctxmlib-appb-M000077
in representing the soil textural properties is mathematically analyzed and the predicted relationship is shown in equation 3. The equation for STI was proposed in such a way to extract the normalized product ranging from 0 to 1 (Shafian and Maas, 2015) .
Figure pctxmlib-appb-M000078
(3)
Where the
Figure pctxmlib-appb-M000079
stands for the value of soil textural index,
Figure pctxmlib-appb-M000080
is soil moisture value in the i th pixel of the corresponding imagery,
Figure pctxmlib-appb-M000081
and
Figure pctxmlib-appb-M000082
stands for the maximum and minimum value of soil moisture.
Three-dimensional Density space for Dielectric Constant
Semi-empirical models are provided to estimate
Figure pctxmlib-appb-M000083
(dielectric constant of non-saline soil from SAR data) and
Figure pctxmlib-appb-M000084
(dielectric constant of saline affected soil from SAR data) by simulating the inputs such as field-observed dielectric constant measurements (
Figure pctxmlib-appb-M000085
), STI, and
Figure pctxmlib-appb-M000086
as three axial parameters in a three-dimensional density plot. As mentioned earlier, with the average field readings on rms-height (0.27 cm) and correlation length (3.1 cm), the surface under investigation has been assumed to be smooth in a semi-saturated condition. Initially, the relationship of STI with
Figure pctxmlib-appb-M000087
in representing dielectric constant has been examined in three-dimensional density space by keeping the values of the third dimension (
Figure pctxmlib-appb-M000088
) as null (Fig. 4a). Based on the negative relationship that the pixels of STI have followed with corresponding pixels of
Figure pctxmlib-appb-M000089
in representing dielectric constant, the model of 2-Dimensional Spatial Distance Measure ( 2D-SDM) is invented by simulating the gradient starting from the pixels having maximum
Figure pctxmlib-appb-M000090
, and minimum STI values to the pixels distributed at the minimum range of
Figure pctxmlib-appb-M000091
, and maximum range of STI (Fig. 4a). The equation for the gradient of 2D-SDM is developed in such a way that it increases with an increasing dielectric constant in the soil (Equation 4). The two-dimensional spatial distance measure could be considered as a preliminary theoretical simulation to indicate the dielectric constant of the soil according to the relationship of polarization property of a signal with the grain-size distribution of the soil.
Figure pctxmlib-appb-M000092
(4)
In three-dimensional feature space, the products of STI, and
Figure pctxmlib-appb-M000093
make up the x-axis and y-axis, and the corresponding values of
Figure pctxmlib-appb-M000094
are taken up in the z-axis, respectively. From the feature space, the z-axis is considered as a conductivity domain, while the x- and y-axis are considered as a polarization domain. The
Figure pctxmlib-appb-M000095
is shown an increasing trend in two different patterns when the property of conductivity is introduced in third-dimension of the existing polarization domain (Fig. 4b) representing the behavior of saline and non-saline soils.
To simulate the model for
Figure pctxmlib-appb-M000096
(dielectric constant of non-saline soil from SAR data) and
Figure pctxmlib-appb-M000097
(dielectric constant of saline soil from SAR data), the gradient of 2-dimensional spatial distance is given elevation in the z-axis (
Figure pctxmlib-appb-M000098
) by two different angles of optimization, namely optimization angle of saline soil (
Figure pctxmlib-appb-M000099
) and non-saline soil (
Figure pctxmlib-appb-M000100
). The pixels representing non-saline soil are decoupled from saline soil and represented in a separate three-dimensional density space with the corresponding STI and
Figure pctxmlib-appb-M000101
values. The outliers are excluded from the sampling points, and the shape of the three-dimensional density space for non-saline soil is approximated to a shape of a right angle triangle in which pixels oriented along the line indicating 3-Dimensional Spatial Distance Measure (
Figure pctxmlib-appb-M000102
) are considered as hypotenuse,
Figure pctxmlib-appb-M000103
are considered as an adjacent side, and
Figure pctxmlib-appb-M000104
are considered as an opposite side to
Figure pctxmlib-appb-M000105
. Because of the constant of proportionality found between the pixels distributed along
Figure pctxmlib-appb-M000106
, and
Figure pctxmlib-appb-M000107
there is no much variation observed in
Figure pctxmlib-appb-M000108
throughout its trend during the upliftment of the pixels from two-dimension to three-dimension. By analyzing the constant rate of the gradient from
Figure pctxmlib-appb-M000109
to
Figure pctxmlib-appb-M000110
,
Figure pctxmlib-appb-M000111
is computed by taking the inverse tangent between the mean of maximum dielectric constant (
Figure pctxmlib-appb-M000112
readings recorded in the field for non-saline soils and the maximum computed value of the two-dimensional spatial measure (
Figure pctxmlib-appb-M000113
).
Once the angle of optimization is determined, a model for
Figure pctxmlib-appb-M000114
(3-Dimensional Spatial Distance Measure for Non-Saline soil) is proposed by applying an inverse cosine rule from dividing the resultant product of
Figure pctxmlib-appb-M000115
by the cos of
Figure pctxmlib-appb-M000116
(Equ 5). The model
Figure pctxmlib-appb-M000117
is proposed in such a way that it starts increasing the dielectric values of the non-saline pixels distributed in the SAR imagery as the three-dimensional spatial distance measure increases. The same procedure is employed for the pixels recorded as saline affected to find out the
Figure pctxmlib-appb-M000118
, and to retrieve the
Figure pctxmlib-appb-M000119
(3-Dimensional Spatial Distance Measure for Saline soil) (Fig 4d) (Equ 6 & 6a). The modeled products
Figure pctxmlib-appb-M000120
and
Figure pctxmlib-appb-M000121
are resulted in the imagery representing
Figure pctxmlib-appb-M000122
and
Figure pctxmlib-appb-M000123
.
Figure pctxmlib-appb-M000124
(5)
Figure pctxmlib-appb-M000125
(6)
NS in the equation stands for non-saline soil, and S stands for saline soil. An expanded form of equations 5 and 6 is given in 5a and 6a.
Figure pctxmlib-appb-M000126
(5a)
Figure pctxmlib-appb-M000127
(6a)
Though the equations 5a and 6a seem to be the same, the model arrives different dielectric values for saline and non-saline soils due to the change in the angles of optimization
Figure pctxmlib-appb-M000128
and
Figure pctxmlib-appb-M000129
computed from the 3-dimensional density plot. This is due to the variations of the dielectric constant measurements for saline and non-saline soils recorded in the field. The mean value of
Figure pctxmlib-appb-M000130
is 37 for non-saline soils, and 29 for saline affected soils. Note that
Figure pctxmlib-appb-M000131
mentioned in the equations is an optimized product
Figure pctxmlib-appb-M000132
.
DSDM-SEC (Density Space Dielectric Model for Soil Electrical Conductivity)
The
Figure pctxmlib-appb-M000133
is recorded for both saline and non-saline soils, and when plotted those points against STI and
Figure pctxmlib-appb-M000134
in the density space, we recorded two different sets of pixels alignment in the third dimension (Fig 4b). This demonstrated that the difference in the orientation was due to dielectric loss that occurred in saline soil. The simulated products
Figure pctxmlib-appb-M000135
and
Figure pctxmlib-appb-M000136
have shown promising partition in the third dimension that leads to the simple but effective way of extracting soil EC. With the two different paradigm shift observed for saline and non-saline soils in 3-dimensional density space, the model for extracting the imaginary of part of dielectric constant (
Figure pctxmlib-appb-M000137
) is invented by scaling down the values of
Figure pctxmlib-appb-M000138
(simulated from
Figure pctxmlib-appb-M000139
) from the corresponding values of
Figure pctxmlib-appb-M000140
(simulated from
Figure pctxmlib-appb-M000141
) (Fig. 4b).
Since the model for
Figure pctxmlib-appb-M000142
(Equ 6a) is simulated based on the dielectric loss due to ionic conductivity of samples and polarization behavior of SAR (Fig 4b), the resultant imagery could able to describe the relative permittivity (
Figure pctxmlib-appb-M000143
). Though the model for
Figure pctxmlib-appb-M000144
is simulated by investigating the samples of non-saline soils in polarization and conductivity domain of a three-dimensional scatter plot, the resultant product is influenced only by the presence of moisture content (polarization) but not by loss factor and hence is considered as an effective indicator of the real part of dielectric constant (
Figure pctxmlib-appb-M000145
). In this study, we proposed a model DSDM-SEC (Density Space Dielectric Model for Soil Electrical Conductivity) for retrieving
Figure pctxmlib-appb-M000146
based on the conception that the simulated dielectric constant of saline soil is relatively lesser than that of the non-saline soil and the potential difference of which could yield the dielectric loss factor (Fig. 4b). The model DSDM-SEC is simulated to arrive
Figure pctxmlib-appb-M000147
by scaling down the product of
Figure pctxmlib-appb-M000148
from
Figure pctxmlib-appb-M000149
(Equ 7). The resultant product of
Figure pctxmlib-appb-M000150
could act as an effective indicator of soil EC.
Figure pctxmlib-appb-M000151
(7)
The proposed equation 7 is simple but effective method to retrieve an imaginary part of dielectric constant from the SAR imagery.
illustrates a system for identifying salination of soil by remotely determining the electrical conductivity of soil in vegetated and bare lands. The system comprises a Synthetic-Aperture Radar (SAR) electronics module 101, a Synthetic-Aperture Radar (SAR) antenna module 102 and an application platform module 103. The SAR electronics module further comprises a C-band signal generator 101a, a processor module 101b, a digitizer module 101c, a Ground Range Detected (GRD) module 101d and a communication module 101e. The SAR antenna module further comprises an antenna module 102a, a transmitter module 102b, an acquisition modes module102e, a receiver module 102c and a waveguide radiators module 102d. The application platform module provides an architecture for processing and analysis for Earth Observation, and further comprises a pre-processing module 104 and a calibrated backscattering module 105.
illustrates a calibrated backscattering module. The calibrated backscattering module further comprises a soil moisture information database 105a, soil textural index database 105b, a Spatial Simulation module 105c for dielectric constant of non-saline soils, a Spatial Simulation module 105d for dielectric constant of saline soils and a Spatial simulation module 105e of Density Space Dielectric Model for Soil Electrical Conductivity.
illustrates a graphical plot that enables the determination of soil moisture content.
illustrates a three-dimensional graphical plot that enables the determination of dielectric constant value for saline and non-saline soil and dielectric loss.
FIG. 5 illustrates an exemplary implementation of the present invention. The results of SAR derived
Figure pctxmlib-appb-M000152
have shown optimum agreement with ground measurements (exclusive of water bodies) with a correlation coefficient equal to 0.86 and with the least RMSE (0.04) and Bias (0.005) values. When the values of SAR derived STI were plotted against the corresponding values of soil samples, we found that STI values were positively correlated (R 2=0.779) with the soil samples composed with a high percentage of sand like sandy soil, and sandy loam categories. The negative correlation has been recorded between the STI values and the soils with a relatively high percentage of clay fraction (R 2=-0.779). The statistical results showed that the
Figure pctxmlib-appb-M000153
and
Figure pctxmlib-appb-M000154
have maintained optimum correlation with the corresponding samples assigned for training and validation (R 2 (
Figure pctxmlib-appb-M000155
) = 0.839, 0.702; R 2 (
Figure pctxmlib-appb-M000156
) = 0.858, 0.707). The statistical significance of
Figure pctxmlib-appb-M000157
in demonstrating soil EC values was examined for different soil categories (Fig 5). The statistical bias in the relationship has been found decreasing as the percentage of free water increases in the soil. Among four soil categories under investigation, sandy soil has attained good correlation (R2=0.883), least RMSE (0.618), and biased value (-0.112) with ground measurements. The mean values of the correlation coefficient between soil EC values and
Figure pctxmlib-appb-M000158
were found to be higher (>R 2=0.7) with acceptable variance (
Figure pctxmlib-appb-M000159
=0.005) for the soil patches uncovered with vegetation. The accuracy level remains approximately the same for the soil covered with the crops of less than 0.5m height.
Although the embodiments herein are described with various specific embodiments, it will be obvious for a person skilled in the art to practice the embodiments herein with modifications.
ADVANTAGES OF THE INVENTION
The various embodiments of the present invention provide a system and method for identifying salination of soil by remotely determining the electrical conductivity (EC) of soil in vegetated and bare lands. The embodiments also provide a system and method for identifying salination of soil by remotely determining the imaginary part of the dielectric constant of the soil through C-band Synthetic-Aperture Radar (SAR). Unlike multiple regression analysis, the present invention yields better results irrespective of the variations in the distribution of the samples. Hence, the invention is employed to retrieve the soil EC with optimum accuracy level even when the field samples are not equally distributed in all ranges of values. The implementation of the present invention does not depend on the different types of field samples and the corresponding spatial data products, which makes it more convenient to employ in soil salinity studies. The present invention enables quantifying the soil Electrical Conductivity of the bare soils and the soils covered with crops of less than 0.5m height. The invention finds application in research groups working on agriculture related use-cases as the resultant product of the model is essential to carry out appropriate agricultural processes according to actual salinity level in the soil.
The foregoing description of the specific embodiments will so fully reveal the general nature of the embodiments herein that others can, by applying current knowledge, readily modify and/or adapt for various applications such as specific embodiments without departing from the generic concept, and, therefore, such adaptations and modifications should and are intended to be comprehended within the meaning and range of equivalents of the disclosed embodiments.
It is to be understood that the phraseology or terminology employed herein is for the purpose of description and not of limitation. Therefore, while the embodiments herein have been described in terms of preferred embodiments, those skilled in the art will recognize that the embodiments herein can be practiced with modifications. However, all such modifications are deemed to be within the scope of the claims.

Claims (10)

  1. A system for identifying salination of soil by remotely determining the electrical conductivity of soil in vegetated lands, the system comprising:
    a Synthetic-Aperture Radar (SAR) electronics module, wherein the SAR electronics module further comprises a C-band signal generator, a processor module, a digitizer module, a Ground Range Detected (GRD) module and a communication module;
    a Synthetic-Aperture Radar (SAR) antenna module, wherein the SAR antenna module further comprises an antenna module, a transmitter module, an acquisition modes module, a receiver module and a waveguide radiators module; and,
    an application platform module, wherein the application platform module provides an architecture for processing and analysis for Earth Observation, and wherein the application platform module further comprises a pre-processing module and a calibrated backscattering module.
  2. The system according to claim 1, wherein the calibrated backscattering module further comprises a soil moisture information database, soil textural index database, a Spatial Simulation module for dielectric constant of non-saline soils, a Spatial Simulation module for dielectric constant of saline soils and a Spatial simulation module of Density Space Dielectric Model for Soil Electrical Conductivity.
  3. The system according to claim 1, wherein the communication module is configured to communicatively couple with a plurality of remote computing devices to determine the type of activities that are suitable on a piece of land based on preset data and observed and/or simulated data from the system, and wherein the communication module is also configured to render the observations and simulated data on a plurality of computing devices in user-intelligible form.
  4. A method for identifying salination of soil by remotely determining the electrical conductivity of soil in vegetated lands, wherein the method provides a Density Space Dielectric Model for Soil Electrical Conductivity (DSDM-SEC), wherein the DSDM-SEC is enabled by three-dimensional density space information obtained as Synthetic-Aperture Radar (SAR) data product in C-band frequency, the DSDM-SEC comprises:
    describing the available total moisture content (
    Figure pctxmlib-appb-M000160
    ) by investigating the scattering and penetration behavior of SAR signals of VV and VH polarizations with an interacting medium in two-dimensional scatterplot;
    determining soil moisture content (
    Figure pctxmlib-appb-M000161
    ) by decoupling the effect of vegetation moisture (
    Figure pctxmlib-appb-M000162
    ) from the
    Figure pctxmlib-appb-M000163
    using DPSVI as a vegetation descriptor;
    providing an “Soil Textural Index (STI)” by mathematically investigating the water holding capacity of the soil in terms of available bound and free water in the soil;
    identifying models for dielectric constant of non-saline soil (
    Figure pctxmlib-appb-M000164
    ) and saline soil (
    Figure pctxmlib-appb-M000165
    ) from parameters such as field-observed dielectric constant measurements (
    Figure pctxmlib-appb-M000166
    ), STI, and
    Figure pctxmlib-appb-M000167
    as three axial parameters in a three-dimensional density plot; and,
    providing an empirical Density Space Dielectric Model for Soil Electrical Conductivity and identifying Imaginary part of dielectric constant (
    Figure pctxmlib-appb-M000168
    ) the soil by scaling down the values of
    Figure pctxmlib-appb-M000169
    from
    Figure pctxmlib-appb-M000170
    ;
  5. The method according to claim 4, wherein the empirical model to estimate the soil moisture further comprises:
    constructing a theoretical scatterplot to the shape of a triangle in which the possible positions of a plurality of features such as Water bodies (WB), Dense Vegetation (D-Veg), Sparse Vegetation (S-Veg), Moisture Soil (MS) and Bare Soil without moisture (BS) are located according the de-polarization characteristics of the SAR signal, and wherein the scattering behavior in the scatter plot was analyzed from VH polarization (y-axis) and VV polarization (x-axis);
    constructing the functional two-dimensional scatterplot between VH and VV polarization to cross verify the theoretical assumptions of the possible positions of the different features; and,
    representing the total moisture content as a simulation of diagonal distance measure (from BS to WB and D-Veg) obtained from the approximated triangle in the two-dimensional scatter plot.
    representing the soil moisture content by reducing the influence of vegetation moisture from total moisture values by considering DPSVI as a vegetative descriptor.
  6. The method according to claim 4, wherein the empirical model “Soil Textural Index (STI)”is provided for mapping and representing the soil textural characteristics as a relative measure from varying moisture content according to sand and clay fraction in the soil under semi-saturated condition.
  7. The method according to claim 4, wherein identifying models for dielectric constant of non-saline soil (
    Figure pctxmlib-appb-M000171
    ) and saline soil (
    Figure pctxmlib-appb-M000172
    ) further includes:
    constructing a three-dimensional density space from inputs such as field-observed dielectric constant measurements (
    Figure pctxmlib-appb-M000173
    ), STI, and
    Figure pctxmlib-appb-M000174
    as three axial parameters in which the z-axis (
    Figure pctxmlib-appb-M000175
    ) is considered as a conductivity domain, while the x- and y-axis (STI, and
    Figure pctxmlib-appb-M000176
    ) are considered as a polarization domain;
    examining the relationship of STI with
    Figure pctxmlib-appb-M000177
    in representing dielectric constant in three-dimensional density space by keeping the values of the third dimension (
    Figure pctxmlib-appb-M000178
    ) as null, and wherein the model of 2-Dimensional Spatial Distance Measure (2D-SDM) is provided by simulating the gradient starting from the pixels having maximum
    Figure pctxmlib-appb-M000179
    , and minimum STI values to the pixels distributed at the minimum range of
    Figure pctxmlib-appb-M000180
    , and maximum range of STI to indicate the dielectric constant of the soil according to the polarization property;
    providing an elevation to the gradient of 2-Dimensional Spatial Distance Measure (2D-SDM) in the z-axis (
    Figure pctxmlib-appb-M000181
    ) by two different angles of optimization, namely optimization angle of saline soil (
    Figure pctxmlib-appb-M000182
    ) and non-saline soil (
    Figure pctxmlib-appb-M000183
    ) based on the mean values of maximum 2D-SDM and
    Figure pctxmlib-appb-M000184
    of saline (
    Figure pctxmlib-appb-M000185
    ) and non-saline soils (
    Figure pctxmlib-appb-M000186
    ); and,
    identifying the 3D-SDM(NS) (3-Dimensional Spatial Distance Measure for Non-Saline soil) or
    Figure pctxmlib-appb-M000187
    and 3D-SDM(S) (3-Dimensional Spatial Distance Measure for Saline soil) or
    Figure pctxmlib-appb-M000188
    by applying an inverse cosine rule from dividing the resultant product of 2D-SDM by the cos of
    Figure pctxmlib-appb-M000189
    and
    Figure pctxmlib-appb-M000190
    to represent dielectric constant of non-saline soil (
    Figure pctxmlib-appb-M000191
    ) and saline soil (
    Figure pctxmlib-appb-M000192
    ).
  8. The method according to claim 4, wherein an empirical model DSDM-SEC (Density Space Dielectric Model for Soil Electrical Conductivity) is provided by subtracting the resultant products of dielectric constant of saline soil (
    Figure pctxmlib-appb-M000193
    ) from the corresponding product of non-saline soil (
    Figure pctxmlib-appb-M000194
    ).
  9. The method according to claim 4, wherein the model DSDM-SEC (Density Space Dielectric Model for Soil Electrical Conductivity) is dependent on soil textural characteristics under the semi-saturated condition.
  10. The method according to claim 4, wherein the model DSDM-SEC (Density Space Dielectric Model for Soil Electrical Conductivity) is potential to quantify the soil Electrical Conductivity of the bare soils and soils covered with the crops of less than 0.5m height.
PCT/IN2020/051018 2020-08-20 2020-12-11 System and method for remote quantification of electrical conductivity of soil WO2022038623A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
IN202041035999 2020-08-20
IN202041035999 2020-08-20

Publications (1)

Publication Number Publication Date
WO2022038623A1 true WO2022038623A1 (en) 2022-02-24

Family

ID=80323298

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/IN2020/051018 WO2022038623A1 (en) 2020-08-20 2020-12-11 System and method for remote quantification of electrical conductivity of soil

Country Status (1)

Country Link
WO (1) WO2022038623A1 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115203978A (en) * 2022-09-01 2022-10-18 中国科学院西北生态环境资源研究院 Improved soil semi-empirical dielectric model based on Dobson dielectric model
CN116310842A (en) * 2023-05-15 2023-06-23 菏泽市国土综合整治服务中心 Soil saline-alkali area identification and division method based on remote sensing image

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7039523B2 (en) * 2001-07-06 2006-05-02 Gecoz Pty Ltd. Method of determining salinity of an area of soil
CN101614818B (en) * 2009-07-09 2012-01-04 中国科学院遥感应用研究所 Radar remote sensing-based detection method of soil alkalization

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7039523B2 (en) * 2001-07-06 2006-05-02 Gecoz Pty Ltd. Method of determining salinity of an area of soil
CN101614818B (en) * 2009-07-09 2012-01-04 中国科学院遥感应用研究所 Radar remote sensing-based detection method of soil alkalization

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115203978A (en) * 2022-09-01 2022-10-18 中国科学院西北生态环境资源研究院 Improved soil semi-empirical dielectric model based on Dobson dielectric model
CN116310842A (en) * 2023-05-15 2023-06-23 菏泽市国土综合整治服务中心 Soil saline-alkali area identification and division method based on remote sensing image

Similar Documents

Publication Publication Date Title
Minet et al. Mapping shallow soil moisture profiles at the field scale using full-waveform inversion of ground penetrating radar data
Lambot et al. Analysis of air‐launched ground‐penetrating radar techniques to measure the soil surface water content
WO2022038623A1 (en) System and method for remote quantification of electrical conductivity of soil
Yao et al. Comparison of four spatial interpolation methods for estimating soil moisture in a complex terrain catchment
Zhao et al. Mapping cation exchange capacity using a quasi-3d joint inversion of EM38 and EM31 data
Castrignanò et al. Use of EMI, gamma-ray emission and GPS height as multi-sensor data for soil characterisation
Lambot et al. Modeling of ground-penetrating radar for accurate characterization of subsurface electric properties
Shaddad et al. Data fusion techniques for delineation of site-specific management zones in a field in UK
Lambot et al. Effect of soil roughness on the inversion of off‐ground monostatic GPR signal for noninvasive quantification of soil properties
Buttafuoco et al. Geostatistical modelling of within-field soil and yield variability for management zones delineation: a case study in a durum wheat field
Farzamian et al. Mapping soil salinity using electromagnetic conductivity imaging—A comparison of regional and location‐specific calibrations
Valencia et al. Oil slicks detection using GNSS-R
Barca et al. Contribution of EMI and GPR proximal sensing data in soil water content assessment by using linear mixed effects models and geostatistical approaches
Koganti et al. Mapping cation exchange capacity using a Veris-3100 instrument and invVERIS modelling software
Etminan et al. Retrieving root-zone soil moisture profile from P-band radar via hybrid global and local optimization
Lambot et al. Frequency dependence of the soil electromagnetic properties derived from ground-penetrating radar signal inversion
Khongnawang et al. Comparing electromagnetic induction instruments to map soil salinity in two-dimensional cross-sections along the Kham-rean Canal using EM inversion software
Cafarelli et al. A linear mixed effect (LME) model for soil water content estimation based on geophysical sensing: a comparison of an LME model and kriging with external drift
Cavallo et al. Use of geophysical data for assessing 3D soil variation in a durum wheat field and their association with crop yield
Buttafuoco et al. Taking into account change of support when merging heterogeneous spatial data for field partition
Wu et al. A nondestructive conductivity estimating method for saline-alkali land based on ground penetrating radar
Chen et al. An application of Coons patch to generate grid-based digital elevation models
Tyndale-Biscoe et al. A system for collecting spatially variable terrain data
Byndych Using Multispectral Satellite Imagery for Parameterisation of Eroded Chernozem
Attema et al. A radar signature model for partially coherent scattering from irregular surfaces

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 20950210

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 20950210

Country of ref document: EP

Kind code of ref document: A1