CN114065643A - Plant soil water content estimation method and system based on SAR and polarization decomposition - Google Patents
Plant soil water content estimation method and system based on SAR and polarization decomposition Download PDFInfo
- Publication number
- CN114065643A CN114065643A CN202111405759.4A CN202111405759A CN114065643A CN 114065643 A CN114065643 A CN 114065643A CN 202111405759 A CN202111405759 A CN 202111405759A CN 114065643 A CN114065643 A CN 114065643A
- Authority
- CN
- China
- Prior art keywords
- water content
- polarization
- soil water
- soil
- sar
- Prior art date
- Legal status (The legal status 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 status listed.)
- Granted
Links
- 239000002689 soil Substances 0.000 title claims abstract description 119
- 230000010287 polarization Effects 0.000 title claims abstract description 89
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 81
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 39
- 239000011159 matrix material Substances 0.000 claims abstract description 54
- 238000005070 sampling Methods 0.000 claims abstract description 48
- 230000001427 coherent effect Effects 0.000 claims abstract description 17
- 238000012549 training Methods 0.000 claims abstract description 14
- 238000007781 pre-processing Methods 0.000 claims abstract description 5
- 238000012360 testing method Methods 0.000 claims description 17
- 230000006870 function Effects 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000013507 mapping Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 238000012795 verification Methods 0.000 claims description 5
- 238000011156 evaluation Methods 0.000 claims description 3
- 238000004590 computer program Methods 0.000 claims description 2
- 238000005259 measurement Methods 0.000 abstract description 3
- 241000196324 Embryophyta Species 0.000 description 18
- 241000209140 Triticum Species 0.000 description 15
- 235000021307 Triticum Nutrition 0.000 description 15
- 238000011160 research Methods 0.000 description 6
- 238000012544 monitoring process Methods 0.000 description 5
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 244000061456 Solanum tuberosum Species 0.000 description 1
- 235000002595 Solanum tuberosum Nutrition 0.000 description 1
- 238000012271 agricultural production Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 235000013339 cereals Nutrition 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000001704 evaporation Methods 0.000 description 1
- 230000008020 evaporation Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000035558 fertility Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000005251 gamma ray Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 235000012015 potatoes Nutrition 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000007637 random forest analysis Methods 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 238000009331 sowing Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
- 239000002352 surface water Substances 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Systems 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/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Software Systems (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Business, Economics & Management (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Marine Sciences & Fisheries (AREA)
- Algebra (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- General Health & Medical Sciences (AREA)
- Computing Systems (AREA)
- Economics (AREA)
- Health & Medical Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Animal Husbandry (AREA)
- Marketing (AREA)
- Agronomy & Crop Science (AREA)
- Databases & Information Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computer Networks & Wireless Communication (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Electromagnetism (AREA)
- Medical Informatics (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
Abstract
The invention discloses a method and a system for estimating the water content of plant soil based on SAR and polarization decomposition, which comprises the steps of firstly obtaining SAR image data and preprocessing the SAR image data; extracting a coherent matrix T3; calculating a surface scattering matrix TG and a backscattering coefficient which represent the surface scattering components; calculating a surface backscattering coefficient by utilizing a soil moisture estimation data set through a training set; simulating a surface backscattering coefficient; obtaining the soil water content of each sampling point by using a lookup table and a minimum cost function strategy; and an effective roughness parameter; and finally estimating the soil water content of the plants. According to the method, the water content of the plant soil is estimated through SAR and polarization decomposition, the effective roughness parameter is used for representing the soil roughness of the plant growing area, the solution of the water content of the soil is simplified, the actual measurement roughness parameter is not relied on, meanwhile, the backscattering component of the earth surface is obtained through polarization decomposition, the influence of vegetation scattering contribution is avoided, and the accuracy of estimating the water content of the soil is improved.
Description
Technical Field
The invention relates to the technical field of remote sensing image data processing, in particular to a method and a system for estimating the water content of plant soil based on SAR and polarization decomposition.
Background
Soil moisture is an important component of surface ecological water circulation and controls surface runoff and surface water evaporation. In agricultural production, soil moisture is an important source of moisture absorbed by crops, is an important factor influencing soil fertility, and is also an important condition for crop survival. The volume water content of the soil refers to the ratio of the volume of water in the soil to the total volume of the soil. The traditional soil water content monitoring method is time-consuming and labor-consuming, only information of a limited point position can be obtained, large-area soil water monitoring cannot be achieved, and an effective way is provided for soil water content estimation due to the appearance of a remote sensing technology. The quality of optical remote sensing data is often influenced by weather conditions such as cloud and fog, and in contrast, microwave remote sensing has strong penetrating power and can realize all-weather monitoring of ground objects.
In addition, the surface microwave radiation is very sensitive to the change of soil moisture, so the microwave remote sensing is gradually the main means for monitoring the soil moisture, wherein the Synthetic Aperture Radar (SAR) has higher spatial resolution, can provide multi-angle and multi-polarization data, and has become the main means for actively monitoring the soil moisture by the microwave remote sensing. Wheat is the grain crop with the largest sowing area, the largest output and the widest distribution in the world, so that the accurate and reliable acquisition of the soil moisture information of the wheat planting area has important guiding significance for agricultural management.
Currently, many scholars construct soil water content inversion models based on SAR, and the research mainly focuses on the relationship between the soil water content and SAR backscattering coefficients. Common models are empirical or semi-empirical models, such as the Oh model, the Dubois model; physical models, such as IEM models; machine learning models such as support vector machines, random forests, artificial neural networks, and the like. For most physical models or empirical/semi-empirical models, the model construction usually requires many additional input parameters, such as root mean square height and correlation length, representing roughness parameters, but these roughness parameters are difficult to measure and are subject to human interference, especially in heavily vegetated areas, leading to uncertainty in soil moisture estimation.
Furthermore, in vegetation covered areas, such as wheat growing areas, the presence of wheat will produce backscattering and attenuate backscattering from the surface, so that the observed backscattering coefficient contains multiple scattering components, making estimation of soil moisture under vegetation coverage difficult.
Disclosure of Invention
In view of the above, the present invention provides a method and a system for estimating water content in plant soil based on SAR and polarization decomposition, in which the method utilizes polarization information provided by SAR and adopts polarization decomposition technique to decompose a full polarization signal into scattering components representing different scattering components, thereby eliminating scattering contribution of vegetation and further estimating water content in soil.
In order to achieve the purpose, the invention provides the following technical scheme:
the invention provides a plant soil water content estimation method based on SAR and polarization decomposition, which comprises the following steps:
acquiring and preprocessing SAR image data, wherein the image data comprises a plant and soil image;
extracting a coherent matrix T3 according to the image data;
calculating a surface scattering matrix TG representing the surface scattering components according to the coherent matrix T3;
calculating the backscattering coefficients of VV polarization and HH polarization of the earth surface according to the surface scattering matrix TG;
constructing a soil moisture estimation data set, and determining a training set and a testing set;
calculating the surface backscattering coefficient through a training set;
constructing a surface scattering model, and simulating a surface backscattering coefficient through the surface scattering model;
constructing a lookup table of the earth surface scattering model;
based on a minimum cost function strategy, inverting from the lookup table to obtain the soil water content of each sampling point;
calculating an effective roughness parameter at the current sampling date;
and constructing a lookup table of the earth surface scattering model according to the effective roughness parameters, and estimating the soil water content of the earth surface covered by the plant by combining the observed backscattering coefficient and the polarization decomposition method.
Further, the training set is determined by obtaining measured soil moisture content data.
Further, the calculation of the effective roughness parameter at the current sampling date is performed according to the following steps:
calculating the estimation precision of the current sampling date under the roughness condition;
and traversing all the given roughness parameter values, repeating the cycle, and taking the roughness parameter with the estimation precision meeting the preset condition as the effective roughness parameter at the current sampling date.
Further, the method also comprises the following steps:
the evaluation and the precision verification of the soil water content of the test set are as follows:
and for each sampling date, constructing a lookup table of the earth surface scattering model by using the effective roughness parameters, obtaining the estimated soil water content on the test set on the basis of the surface backscattering coefficient obtained by polarization decomposition and the backscattering coefficient simulated by the earth surface scattering model on the test set in combination with a minimum cost function strategy, comparing the estimated soil water content with the actually measured soil water content of the test set, carrying out precision verification and judging the estimation precision.
Further, the method also comprises the following steps:
regional estimation and mapping of the water content of plant soil are specifically as follows:
and for each sampling date, constructing a backscattering coefficient lookup table of a CIEM (common information element) model or a Dubois model by using the effective roughness parameters, traversing all pixels of the plant area in the image data, estimating the soil water content information of each pixel by using a soil water content estimation method, and finally completing the soil water content estimation of the plant area and carrying out soil water content area mapping.
Further, the calculation of the backscattering coefficient is performed according to the following steps:
extracting a coherent matrix T3 of each sampling point from the image data according to the longitude and latitude information of the sampling point;
calculating a polarization coherence matrix TG corresponding to surface scattering from the coherence matrix T3 by using different volume scattering matrices;
based on a non-negative eigenvalue decomposition method and a conversion relation between a polarization coherent matrix and a backscattering coefficient, calculating the backscattering coefficient corresponding to the surface scattering component according to the following formula:
wherein TG represents a polarization coherence matrix;
TGijis the ith row and jth column element of TG;
backscattering coefficients corresponding to surface scattering components in the HH polarization mode and the VV polarization mode are represented by Intensity (Intensity);
further, the determination of the effective roughness parameter is performed according to the following steps:
for different sampling dates, constructing a surface scattering model based on different roughness parameters to obtain a backscattering coefficient lookup table of VV polarization and HH polarization:
under the conditions of different roughness parameters and different soil water contents, simulating the backscattering coefficients of VV polarization and HH polarization of the bare earth surface;
on the training set, obtaining an effective surface backscattering coefficient and a surface backscattering coefficient simulated by a surface scattering model according to polarization decomposition, and obtaining the soil water content of each sampling point from a lookup table by inversion according to a minimum cost function strategy shown by the following formula:
wherein,backscattering coefficients respectively representing VV polarization and HH polarization simulated by the earth surface scattering model;
when the F value is the minimum value, taking the soil water content corresponding to the current lookup table element as the estimated soil water content of the sampling point;
calculating the root mean square error between the actually measured soil water content and the estimated soil water content under the roughness condition, and taking the root mean square error as the estimation precision of the current sampling date;
and traversing all the given roughness parameter values, repeating the cycle, and taking the roughness parameter with the precision meeting the requirement as the effective roughness parameter at the sampling date.
The invention provides a plant soil water content estimation system based on SAR and polarization decomposition, which comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the processor executes the program to realize the steps of the method of any one of claims 1-7.
The invention has the beneficial effects that:
compared with the traditional SAR-based soil moisture estimation model, the SAR-and polarization-decomposition-based plant soil moisture estimation method uses the effective roughness parameter to represent the soil roughness of the plant growing region, reduces the unknown parameters of the earth surface scattering model, simplifies the solution of the soil moisture content, does not depend on the actually measured roughness parameter, and uses the polarization decomposition method to decompose the observed backscattering signals into scattering components representing different components, thereby obtaining the backscattering components of the earth surface, avoiding the influence of vegetation scattering contribution and improving the accuracy of soil moisture estimation.
Additional advantages, objects, and features of the invention will be set forth in part in the description which follows and in part will become apparent to those having ordinary skill in the art upon examination of the following or may be learned from practice of the invention. The objectives and other advantages of the invention may be realized and attained by the means of the instrumentalities and combinations particularly pointed out hereinafter.
Drawings
In order to make the object, technical scheme and beneficial effect of the invention more clear, the invention provides the following drawings for explanation:
FIG. 1 is a schematic flow diagram of the overall process of the present invention.
FIG. 2 is a schematic diagram of a research area of a Middling farm in an embodiment of the present invention, wherein the selected research area is a wheat growing area in the southeast of Ontario, Canada.
FIG. 3 is a graph showing the accuracy of estimation on a test set using two different surface scattering models in an embodiment of the present invention
Fig. 4 is a regional plot of soil water content for wheat growing areas at different sampling dates obtained in an embodiment of the present invention.
Detailed Description
The present invention is further described with reference to the following drawings and specific examples so that those skilled in the art can better understand the present invention and can practice the present invention, but the examples are not intended to limit the present invention.
Example 1
As shown in fig. 1, the method for estimating the water content of plant soil based on SAR and polarization decomposition provided by this embodiment includes the following steps:
step 1: preprocessing an SAR image, and extracting a T3 matrix; t3 denotes a coherence matrix.
For radar targets, the scattering properties can be fully characterized by a polarization scattering matrix S,
wherein VV, VH, HH represent the backscattering coefficients of VV polarization, VH polarization, HH polarization observed by radar, respectively, and the unit is Intensity (Intensity).
If the polarization scattering matrix is vectorized on the basis of Pauli basis, its scattering matrix vector kp is expressed as,
where T denotes a matrix transpose. The polarization coherence matrix T3 corresponding to the n-view processed scattering vector kp can be expressed as:
wherein denotes a matrix conjugate;
step 2: a surface scattering matrix TG characterizing the earth surface scattering components is calculated from the T3 matrix based on the vertical volume scattering matrix (TVv) and the horizontal volume scattering matrix (TVh), and the backscattering coefficients of VV polarization and HH polarization of the earth surface are further calculated based on the TG.
Based on the polarization decomposition method of the three-component mechanism (Freeman decomposition), the T3 matrix can be expressed in the form,
wherein T isijRow i and column j elements representing the coherence matrix T3; TG, TD and TV are polarization coherent matrixes corresponding to surface scattering, secondary scattering and bulk scattering respectively; beta is the normalized difference of H-polarized and V-polarized Bragg scattering; alpha is the normalized difference of Fresnel scattering of H-polarized and V-polarized earth-vegetation stems; f. ofs、fdThe intensity coefficients of the surface scattering component and the secondary scattering component are respectively; "+" denotes the conjugate transpose.
The vertical volume scattering matrix (TVv) and the horizontal volume scattering matrix (TVh) proposed by Yamaguchi are selected as the volume scattering matrix,
wherein, VijRow i and column j elements representing a volume scattering matrix; f. ofvRepresenting the intensity coefficient of the bulk scatter component.
Through the analysis of the scattering mechanism components, the proportion of the secondary scattering component TD is negligible, and T3 can be expressed as T3 ═ TG + TV, so TG can be expressed as: TG is expressed as a band unknown parameter f, T3-TVvBased on a non-negative eigenvalue decomposition method (NNED), the minimum f when the minimum characteristic of TG is not less than 0 is obtainedvAs coefficient of bulk scattering intensity fv。
After solving a coherent matrix TG corresponding to the surface scattering, calculating a backscattering coefficient corresponding to the surface scattering component based on a conversion relation between the polarized coherent matrix and the backscattering coefficient as shown in the following formula:
wherein TG represents a polarization coherence matrix; TG (gamma-ray) in a single phaseijIs the ith row and jth column element of TG;calculating a backscattering coefficient corresponding to a surface scattering component of VV polarization and an HH polarization mode based on a polarization decomposition method, wherein the unit is Intensity (Intensity); selecting one obtained by calculation based on vertical body scattering matrixSelecting the back scattering coefficient obtained by calculation based on a horizontal volume scattering matrix as the back scattering coefficient of the surface in the VV polarization modeThe surface backscatter coefficient is expressed as the HH polarization mode.
And step 3: calculated surface backscattering coefficientAnd constructing a soil moisture estimation data set together with the actually measured soil moisture. The data of each day is divided into training data and testing data according to a certain proportion.
And 4, step 4: for a certain sampling date, on training data, a look-up table of a surface scattering model (CIEM or Dubois) is constructed based on different roughness parameters (root mean square height): and the simulated VV polarization of the bare earth surface and the backscattering coefficient of HH polarization under different root mean square heights under different soil water contents.
And then calculating according to polarization decomposition to obtain a surface backscattering coefficient and a surface backscattering coefficient obtained by simulating a scattering model, and inverting from the lookup table to obtain the soil water content of each sampling point based on a minimum cost function strategy.
Wherein,backscattering coefficient values respectively representing VV polarization and HH polarization simulated by the earth surface scattering model; and when the F is minimum, taking the soil moisture content corresponding to the simulated backscattering coefficient as the estimated soil moisture content value of the sampling point. Then, under the roughness condition, the Root Mean Square Error (RMSE) between the estimated soil moisture content and the actually measured soil moisture content on the current sampling date is calculated as the estimation accuracy on the current date.
Wherein,and yiRespectively representing the estimated soil moisture content and the measured soil moisture content of the ith sampling point.
And traversing all the given roughness parameter values, repeating the step, and taking the roughness parameter with the highest estimation precision as the effective roughness parameter at the sampling date.
And 5: for a certain sampling date, a lookup table of a surface scattering model (CIEM or Dubois) is constructed by using effective roughness parameters, all pixels of a plant coverage area of the SAR image are traversed to obtain an observed total backscattering coefficient, then surface backscattering coefficients of all pixels in the area are obtained based on the steps 1 and 2, and finally soil water content of each pixel in the plant area is obtained through inversion according to a minimum cost function strategy.
The plant provided by the embodiment can be a farmland for planting wheat, and can also be plants for planting rice, potatoes, flowers and plants, trees and the like.
Example 2
The method for estimating the soil water content of the wheat farmland based on SAR and polarization decomposition provided by the embodiment is further described as follows:
data acquisition: the remote sensing data adopted by the invention is RADARSAT-2 full-polarization SAR backscattering coefficient data, the actual measurement data comprises soil water content information of a wheat planting research area, the remote sensing data and the actual measurement data are acquired on the same date, eight dates in the growth period of the wheat are covered, and the specific date is as shown in figure 4.
The experimental procedure is shown in FIG. 1 and is described in detail below.
The method comprises the following steps: preprocessing the remote sensing image: the method mainly comprises radiometric calibration, coherent matrix extraction T3, polarization filtering, terrain correction and projection conversion.
Step two: solving the surface backscattering coefficient based on a polarization decomposition method:
extracting a coherent matrix T3 of each sampling point from a T3 image according to longitude and latitude information of the sampling point, then calculating a polarized coherent matrix (TG) corresponding to surface scattering from the coherent matrix T3 by using different volume scattering matrixes (a vertical volume scattering matrix and a horizontal volume scattering matrix), and then calculating a backscattering coefficient corresponding to a surface scattering component based on a non-negative eigenvalue decomposition method (NNED) and a conversion relation between the polarized coherent matrix and the backscattering coefficient.
Wherein TG represents a polarization coherence matrix;
TGijis the ith row and jth column element of TG;
backscattering coefficients corresponding to surface scattering components in the HH polarization mode and the VV polarization mode are represented by Intensity (Intensity);
in step two, the matrix is obtained by using a scattering matrix based on a vertical bodyAs effective VV surface scattering component of wheat research area, obtained based on horizontal volume scattering matrixAs an effective HH surface scattering component in wheat research.
Step three: and (3) constructing a soil moisture estimation data set:
and (3) constructing a soil moisture estimation data set by actually measuring soil moisture of the sampling points and the effective surface scattering components (VV, HH) extracted in the second step, and dividing the data into 7 for each sampling date: 3, training set and testing set;
step four: estimating soil water content and determining effective roughness parameters in a training set:
for different sampling dates, constructing a surface scattering model (a CIEM model and a Dubois model) based on different roughness parameters (root mean square height) to obtain a backscattering coefficient lookup table of VV polarization and HH polarization: and simulating the backscattering coefficients of VV polarization and HH polarization of the bare earth surface under different roughness parameters and different soil water contents. On the training set, an effective surface backscattering coefficient and a surface backscattering coefficient simulated by a surface scattering model are obtained according to polarization decomposition, and the soil water content of each sampling point is obtained through inversion from a lookup table based on a minimum cost function strategy shown in the specification.
Wherein,representing the VV polarization backscattering coefficient value of the earth surface scattering model simulation;values of HH polarization backscattering coefficients representing a surface scattering model simulation;
and when the F value is the minimum value, taking the soil water content corresponding to the current lookup table element as the estimated soil water content of the sampling point. Then, the Root Mean Square Error (RMSE) between the measured soil moisture content and the estimated soil moisture content under the roughness condition is calculated as the estimation accuracy of the current sampling date. And traversing all given roughness parameter values, repeating the step, and taking the roughness parameter with the highest precision as an effective roughness parameter at the sampling date.
Step five: the evaluation and the precision verification of the soil water content of the test set are as follows:
and for each sampling date, constructing a lookup table of a surface scattering model (CIEM model or Dubois model) by using the effective roughness parameters obtained in the fourth step, obtaining the estimated soil moisture content on the test set on the basis of the surface backscattering coefficient obtained by polarization decomposition and the backscattering coefficient simulated by the surface scattering model on the test set in combination with a minimum cost function strategy, comparing the obtained soil moisture content with the actually measured soil moisture content of the test set, and verifying the precision (as shown in figure 3) so as to judge the effectiveness and the estimation precision of the method.
Step six: regional estimation and mapping of water content of wheat farmland soil
For each sampling date, a backscattering coefficient lookup table of a surface scattering model (a CIEM model or a Dubois model) is constructed by using the effective roughness parameters, then all pixels of a wheat planting area in the RADARSAT-2 image are traversed, soil water content information of each pixel is estimated by using the soil water content estimation method in the steps, and finally, soil water content estimation of a wheat farmland area is completed and soil water content area mapping is carried out (as shown in figure 4).
The above-mentioned embodiments are merely preferred embodiments for fully illustrating the present invention, and the scope of the present invention is not limited thereto. The equivalent substitution or change made by the technical personnel in the technical field on the basis of the invention is all within the protection scope of the invention. The protection scope of the invention is subject to the claims.
Claims (8)
1. The plant soil water content estimation method based on SAR and polarization decomposition is characterized by comprising the following steps: the method comprises the following steps:
acquiring and preprocessing SAR image data, wherein the image data comprises a plant and soil image;
extracting a coherent matrix T3 according to the image data;
calculating a surface scattering matrix TG representing the surface scattering components according to the coherent matrix T3;
calculating the backscattering coefficients of VV polarization and HH polarization of the earth surface according to the surface scattering matrix TG;
constructing a soil moisture estimation data set, and determining a training set and a testing set;
calculating the surface backscattering coefficient through a training set;
constructing a surface scattering model, and simulating a surface backscattering coefficient through the surface scattering model;
constructing a lookup table of the earth surface scattering model;
based on a minimum cost function strategy, inverting from the lookup table to obtain the soil water content of each sampling point;
calculating an effective roughness parameter at the current sampling date;
and constructing a lookup table of the earth surface scattering model according to the effective roughness parameters, and estimating the soil water content of the earth surface covered by the plant by combining the observed backscattering coefficient and the polarization decomposition method.
2. The SAR and polarization decomposition based plant soil water content estimation method of claim 1, characterized in that: the training set is determined by obtaining measured soil moisture content data.
3. The SAR and polarization decomposition based plant soil water content estimation method of claim 1, characterized in that: the calculation of the effective roughness parameter at the current sampling date is carried out according to the following steps:
calculating the estimation precision of the current sampling date under the roughness condition;
and traversing all the given roughness parameter values, repeating the cycle, and taking the roughness parameter with the estimation precision meeting the preset condition as the effective roughness parameter at the current sampling date.
4. The SAR and polarization decomposition based plant soil water content estimation method of claim 1, characterized in that: further comprising the steps of:
the evaluation and the precision verification of the soil water content of the test set are as follows:
and for each sampling date, constructing a lookup table of the earth surface scattering model by using the effective roughness parameters, obtaining the estimated soil water content on the test set on the basis of the surface backscattering coefficient obtained by polarization decomposition and the backscattering coefficient simulated by the earth surface scattering model on the test set in combination with a minimum cost function strategy, comparing the estimated soil water content with the actually measured soil water content of the test set, carrying out precision verification and judging the estimation precision.
5. The SAR and polarization decomposition based plant soil water content estimation method of claim 1, characterized in that: further comprising the steps of:
regional estimation and mapping of the water content of plant soil are specifically as follows:
and for each sampling date, constructing a backscattering coefficient lookup table of a CIEM (common information element) model or a Dubois model by using the effective roughness parameters, traversing all pixels of the plant area in the image data, estimating the soil water content information of each pixel by using a soil water content estimation method, and finally completing the soil water content estimation of the plant area and carrying out soil water content area mapping.
6. The SAR and polarization decomposition based plant soil water content estimation method of claim 1, characterized in that: the calculation of the backscattering coefficient is carried out according to the following steps:
extracting a coherent matrix T3 of each sampling point from the image data according to the longitude and latitude information of the sampling point;
calculating a polarization coherence matrix TG corresponding to surface scattering from the coherence matrix T3 by using different volume scattering matrices;
based on a non-negative eigenvalue decomposition method and a conversion relation between a polarization coherent matrix and a backscattering coefficient, calculating the backscattering coefficient corresponding to the surface scattering component according to the following formula:
wherein TG represents a polarization coherence matrix;
TGijis the ith row and jth column element of TG;
7. The SAR and polarization decomposition based plant soil water content estimation method of claim 1, characterized in that: the effective roughness parameter is determined according to the following steps:
for different sampling dates, constructing a surface scattering model based on different roughness parameters to obtain a backscattering coefficient lookup table of VV polarization and HH polarization:
under the conditions of different roughness parameters and different soil water contents, simulating the backscattering coefficients of VV polarization and HH polarization of the bare earth surface;
on the training set, obtaining an effective surface backscattering coefficient and a surface backscattering coefficient simulated by a surface scattering model according to polarization decomposition, and obtaining the soil water content of each sampling point from a lookup table by inversion according to a minimum cost function strategy shown by the following formula:
wherein,backscattering coefficients respectively representing VV polarization and HH polarization simulated by the earth surface scattering model;
when the F value is the minimum value, taking the soil water content corresponding to the current lookup table element as the estimated soil water content of the sampling point;
calculating the root mean square error between the actually measured soil water content and the estimated soil water content under the roughness condition, and taking the root mean square error as the estimation precision of the current sampling date;
and traversing all the given roughness parameter values, repeating the cycle, and taking the roughness parameter with the precision meeting the requirement as the effective roughness parameter at the sampling date.
8. A plant soil water content estimation system based on SAR and polarization decomposition, comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the steps of the method of any one of claims 1 to 7 when executing the program.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111405759.4A CN114065643B (en) | 2021-11-24 | 2021-11-24 | SAR and polarization decomposition-based plant soil water content estimation method and system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111405759.4A CN114065643B (en) | 2021-11-24 | 2021-11-24 | SAR and polarization decomposition-based plant soil water content estimation method and system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114065643A true CN114065643A (en) | 2022-02-18 |
CN114065643B CN114065643B (en) | 2024-06-28 |
Family
ID=80277630
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111405759.4A Active CN114065643B (en) | 2021-11-24 | 2021-11-24 | SAR and polarization decomposition-based plant soil water content estimation method and system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114065643B (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114548277A (en) * | 2022-02-22 | 2022-05-27 | 电子科技大学 | Method and system for fitting ground points and extracting crop height based on point cloud data |
CN115061112A (en) * | 2022-08-09 | 2022-09-16 | 中国科学院地理科学与资源研究所 | Planar snow water equivalent obtaining method and device and electronic equipment |
CN115618174A (en) * | 2022-10-24 | 2023-01-17 | 安徽师范大学 | Soil humidity inversion method based on pixel scale surface roughness spectrum parameters |
CN117826112A (en) * | 2024-03-05 | 2024-04-05 | 天津智云水务科技有限公司 | Soil water content inversion method based on sar |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108681652A (en) * | 2018-06-04 | 2018-10-19 | 中国科学院遥感与数字地球研究所 | A kind of vegetation-covered area Soil Moisture Inversion method based on three number of high score |
CN109521182A (en) * | 2018-10-30 | 2019-03-26 | 武汉大学 | A kind of PolSAR soil moisture content inversion method based on two component decomposition models |
CN111751286A (en) * | 2020-06-22 | 2020-10-09 | 电子科技大学 | Soil moisture extraction method based on change detection algorithm |
-
2021
- 2021-11-24 CN CN202111405759.4A patent/CN114065643B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108681652A (en) * | 2018-06-04 | 2018-10-19 | 中国科学院遥感与数字地球研究所 | A kind of vegetation-covered area Soil Moisture Inversion method based on three number of high score |
CN109521182A (en) * | 2018-10-30 | 2019-03-26 | 武汉大学 | A kind of PolSAR soil moisture content inversion method based on two component decomposition models |
CN111751286A (en) * | 2020-06-22 | 2020-10-09 | 电子科技大学 | Soil moisture extraction method based on change detection algorithm |
Non-Patent Citations (1)
Title |
---|
刘伟;施建成;余琴;: "基于机载极化雷达技术的农作物覆盖区土壤水分估算", 干旱区地理, no. 06, 30 December 2005 (2005-12-30) * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114548277A (en) * | 2022-02-22 | 2022-05-27 | 电子科技大学 | Method and system for fitting ground points and extracting crop height based on point cloud data |
CN114548277B (en) * | 2022-02-22 | 2023-09-08 | 电子科技大学 | Method and system for ground point fitting and crop height extraction based on point cloud data |
CN115061112A (en) * | 2022-08-09 | 2022-09-16 | 中国科学院地理科学与资源研究所 | Planar snow water equivalent obtaining method and device and electronic equipment |
CN115618174A (en) * | 2022-10-24 | 2023-01-17 | 安徽师范大学 | Soil humidity inversion method based on pixel scale surface roughness spectrum parameters |
CN117826112A (en) * | 2024-03-05 | 2024-04-05 | 天津智云水务科技有限公司 | Soil water content inversion method based on sar |
CN117826112B (en) * | 2024-03-05 | 2024-05-31 | 天津智云水务科技有限公司 | Soil water content inversion method based on sar |
Also Published As
Publication number | Publication date |
---|---|
CN114065643B (en) | 2024-06-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Betbeder et al. | Assimilation of LAI and dry biomass data from optical and SAR images into an agro-meteorological model to estimate soybean yield | |
CN114065643B (en) | SAR and polarization decomposition-based plant soil water content estimation method and system | |
Huang et al. | Jointly assimilating MODIS LAI and ET products into the SWAP model for winter wheat yield estimation | |
El Hajj et al. | Soil moisture retrieval over irrigated grassland using X-band SAR data | |
Jongschaap | Run-time calibration of simulation models by integrating remote sensing estimates of leaf area index and canopy nitrogen | |
Olioso et al. | Estimation of evapotranspiration and photosynthesis by assimilation of remote sensing data into SVAT models | |
Dente et al. | Assimilation of leaf area index derived from ASAR and MERIS data into CERES-Wheat model to map wheat yield | |
Leroux et al. | Crop monitoring using vegetation and thermal indices for yield estimates: case study of a rainfed cereal in semi-arid West Africa | |
Huang et al. | BESS-Rice: A remote sensing derived and biophysical process-based rice productivity simulation model | |
WO2011102520A1 (en) | Method of generating paddy rice crop yield forecasting model, and method of forecasting crop yield of paddy rice | |
Wenzhi et al. | Estimation of sunflower seed yield using partial least squares regression and artificial neural network models | |
Yang et al. | In-season biomass estimation of oilseed rape (Brassica napus L.) using fully polarimetric SAR imagery | |
Lee et al. | A combined passive/active microwave remote sensing approach for surface variable retrieval using Tropical Rainfall Measuring Mission observations | |
CN115479960B (en) | Method for monitoring soil humidity in crop growth process by combining SAR (synthetic aperture radar) with optical remote sensing data | |
CN114782834A (en) | Farmland soil moisture retrieval method and system based on ratio method | |
Huang et al. | Assimilating remotely sensed information with the wheatgrow model based on the ensemble square root filter for improving regional wheat yield forecasts | |
CN109615551A (en) | The wheat crops inversion method of leaf area index simulated based on microwave scattering and canopy | |
Fieuzal et al. | Estimation of sunflower yield using a simplified agrometeorological model controlled by optical and SAR satellite data | |
Feng et al. | Performance of the improved two-source energy balance model for estimating evapotranspiration over the heterogeneous surface | |
Zhang et al. | Estimation of rice yield from a C-band radar remote sensing image by integrating a physical scattering model and an optimization algorithm | |
Abimbola et al. | Improving crop modeling to better simulate maize yield variability under different irrigation managements | |
CN115931917A (en) | Soil moisture inversion method and system based on fully-polarized SAR data | |
CN114611699A (en) | Soil moisture downscaling method and device, electronic equipment and storage medium | |
Seidl et al. | Integrating remotely sensed images with a soybean model to improve spatial yield simulation | |
Lewis | 3D canopy modelling as a tool in remote-sensing research |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |