CN115479960B - Method for monitoring soil humidity in crop growth process by combining SAR (synthetic aperture radar) with optical remote sensing data - Google Patents
Method for monitoring soil humidity in crop growth process by combining SAR (synthetic aperture radar) with optical remote sensing data Download PDFInfo
- Publication number
- CN115479960B CN115479960B CN202211111362.9A CN202211111362A CN115479960B CN 115479960 B CN115479960 B CN 115479960B CN 202211111362 A CN202211111362 A CN 202211111362A CN 115479960 B CN115479960 B CN 115479960B
- Authority
- CN
- China
- Prior art keywords
- soil
- sar
- period
- humidity
- time
- 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.)
- Active
Links
- 239000002689 soil Substances 0.000 title claims abstract description 194
- 238000000034 method Methods 0.000 title claims abstract description 65
- 230000003287 optical effect Effects 0.000 title claims abstract description 29
- 238000012544 monitoring process Methods 0.000 title claims abstract description 26
- 230000008569 process Effects 0.000 title claims abstract description 24
- 238000009331 sowing Methods 0.000 claims abstract description 32
- 230000008859 change Effects 0.000 claims abstract description 19
- 230000010287 polarization Effects 0.000 claims description 25
- 241000209140 Triticum Species 0.000 claims description 8
- 235000021307 Triticum Nutrition 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 7
- 235000017060 Arachis glabrata Nutrition 0.000 claims description 3
- 244000105624 Arachis hypogaea Species 0.000 claims description 3
- 235000010777 Arachis hypogaea Nutrition 0.000 claims description 3
- 235000018262 Arachis monticola Nutrition 0.000 claims description 3
- 240000008042 Zea mays Species 0.000 claims description 3
- 235000005824 Zea mays ssp. parviglumis Nutrition 0.000 claims description 3
- 235000002017 Zea mays subsp mays Nutrition 0.000 claims description 3
- 235000005822 corn Nutrition 0.000 claims description 3
- 235000020232 peanut Nutrition 0.000 claims description 3
- 235000014698 Brassica juncea var multisecta Nutrition 0.000 claims 1
- 235000006008 Brassica napus var napus Nutrition 0.000 claims 1
- 240000000385 Brassica napus var. napus Species 0.000 claims 1
- 235000006618 Brassica rapa subsp oleifera Nutrition 0.000 claims 1
- 235000004977 Brassica sinapistrum Nutrition 0.000 claims 1
- 238000004088 simulation Methods 0.000 abstract description 8
- 238000004458 analytical method Methods 0.000 abstract description 5
- 230000001575 pathological effect Effects 0.000 abstract description 3
- 238000004422 calculation algorithm Methods 0.000 description 8
- 238000011160 research Methods 0.000 description 8
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 6
- 230000014509 gene expression Effects 0.000 description 4
- 230000003746 surface roughness Effects 0.000 description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 3
- 241000196324 Embryophyta Species 0.000 description 3
- 238000012271 agricultural production Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000001704 evaporation Methods 0.000 description 2
- 230000008020 evaporation Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 235000015097 nutrients Nutrition 0.000 description 2
- 238000004064 recycling Methods 0.000 description 2
- 238000002310 reflectometry Methods 0.000 description 2
- 238000000611 regression analysis Methods 0.000 description 2
- 208000005156 Dehydration Diseases 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005388 cross polarization Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000004090 dissolution Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 238000003973 irrigation Methods 0.000 description 1
- 230000002262 irrigation Effects 0.000 description 1
- 238000010899 nucleation Methods 0.000 description 1
- 235000016709 nutrition Nutrition 0.000 description 1
- 230000035764 nutrition Effects 0.000 description 1
- 230000029553 photosynthesis Effects 0.000 description 1
- 238000010672 photosynthesis Methods 0.000 description 1
- 238000001556 precipitation Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000002352 surface water Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N22/00—Investigating or analysing materials by the use of microwaves or radio waves, i.e. electromagnetic waves with a wavelength of one millimetre or more
- G01N22/04—Investigating moisture content
-
- 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/86—Combinations of radar systems with non-radar systems, e.g. sonar, direction finder
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A40/00—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production
- Y02A40/10—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production in agriculture
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Electromagnetism (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Databases & Information Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Algebra (AREA)
- Health & Medical Sciences (AREA)
- Software Systems (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
The invention discloses a method for monitoring soil humidity in a crop growth process by combining SAR and optical remote sensing data, which comprises the following steps: s1, based on sowing time T 1 SAR satellite data SAR (T) 1 ,θ 1 ) SAR (T) 1 ,θ 2 ) Calculating a sowing time T 1 Soil surface humidity M of (2) v (T 1 ) The method comprises the steps of carrying out a first treatment on the surface of the S2, based on period T i SAR satellite data SAR (T) i θ), period T i+1 Optical remote sensing data OPT (T) i+1 ) SAR satellite data SAR (T) i+1 θ) determining the period T i And period T i+1 The relative change of the soil humidity in the time phase, and then the time T is obtained i+1 Soil surface humidity M of (2) v (T i+1 ) Wherein i is sequentially from 1 to h; through forward model simulation analysis, an inversion model of the surface humidity of the bare soil is improved; based on a first-order discrete vegetation model, a soil humidity monitoring model based on combined SAR and optical data is constructed, and the problems of pathological inversion and mixed pixels of vegetation coverage soil humidity remote sensing inversion are solved; finally, the time series soil humidity monitoring of the crop growth process independent of ground observation data is realized.
Description
Technical Field
The invention belongs to the technical field of soil environment monitoring, and particularly relates to a method for monitoring soil humidity in a crop growth process by combining SAR (synthetic aperture radar) and optical remote sensing data.
Background
Soil moisture (i.e., soil moisture content) is an important parameter in the study of the earth's system and is one of the core variables affecting the surface process, playing an important role in energy recycling, water recycling, ecosystems, and agricultural production. In the energy circulation of the earth surface and the atmosphere, the change of the soil humidity can influence earth surface parameters such as earth surface albedo, soil heat capacity, earth surface evaporation and the like, change sensible heat, latent heat and long wave radiation flux of earth surface to atmosphere, cause redistribution of earth surface energy, and play a very important role in regulating climate change. Soil moisture is a tie linking surface water and underground water in the process of forming, converting and consuming land water resources, plays a vital role in the hydrologic processes of precipitation, runoff, infiltration, evaporation and the like, and is an important input parameter of various hydrologic models, climate models, ecological models, atmosphere models and the like. Soil moisture is an important component of a Liu Miansheng-state system, is a basic condition for land plants and soil organisms to survive, and moderate soil humidity is beneficial to the dissolution and movement of various nutrients in soil, improves the nutrition condition of plants, and promotes the circulation of the nutrients, thereby influencing the photosynthesis and productivity of vegetation. Meanwhile, soil moisture is a key factor in agricultural production, soil humidity is an important index for researching plant water stress, carrying out drought monitoring, estimating crop yield and the like, is an important parameter of a crop growth state monitoring and estimating model, and can play a guiding role in agricultural production irrigation management. It can be seen that soil moisture is of great importance in the research of the earth science.
According to different sensor types, the current soil humidity remote sensing inversion method is divided into 3 types of optical remote sensing methods, passive microwave remote sensing methods and active microwave remote sensing methods:
1) Optical remote sensing method: soil moisture is estimated using spectral reflectance characteristics of soil surface, emissivity of soil surface, or surface temperature, and there are mainly an index-based method and a soil thermal inertia-based method. The method utilizes the principle that the reflectivity of dry soil is higher, and the reflectivity of similar wet soil in each wave band is correspondingly reduced, and inverts the soil humidity by constructing different spectrum indexes. The indices commonly used herein are mainly a distance-to-flat vegetation index, a temperature vegetation drought index TVDI, a soil humidity index SWI, a vertical drought index PDI, a soil moisture content temperature index SWCTI, a crop water deficiency index CWSI, a vegetation water supply index VSWI, and the like. The soil thermal inertia-based method utilizes good linear relation between the soil thermal inertia and the soil humidity and stronger sensitivity to invert the soil humidity, and has higher precision in the soil drought monitoring with smaller range and single type. The advantage of inverting the soil humidity by the optical remote sensing method is that the spatial resolution is high, and more satellite sensors can be selected; however, the optical wave band can only measure the reflection and emission of 1mm on the earth surface, belongs to a method for indirectly inverting the soil humidity, and cannot penetrate through vegetation canopy and is affected by atmospheric attenuation, so that the method has certain limitation in application.
2) The passive microwave remote sensing method comprises the following steps: soil humidity is inverted by using soil microwave emission or brightness temperature measured by a microwave radiometer. Researches show that for the bare earth surface, a simple linear relation exists between the soil humidity and the microwave emissivity (replaced by the common microwave bright temperature and ground temperature ratio), and the earth humidity of the bare earth surface can be inverted by establishing an empirical model through the linear relation; for vegetation covering the earth surface, the influence of vegetation is quantified by using a vegetation index (which can be a microwave remote sensing vegetation index or an optical remote sensing vegetation index) or a tau-omega model and the like, the soil roughness is described by using an H parameter, and then the soil humidity is inverted by combining the empirical model of bare soil. The multi-band multi-polarization multi-angle soil moisture inversion method is developed, becomes a standard algorithm of an AMSR-E global soil moisture data set and a standard algorithm of an SMOS global soil moisture data set, and is applied in a related way. The passive microwave remote sensing method can penetrate through the vegetation layer to detect soil with a depth of about 5cm below the ground surface, has a strong inversion physical foundation, has large data width and high time resolution, and is suitable for soil humidity inversion and application research on the global scale; in contrast, its low spatial resolution on the order of tens of kilometers greatly affects the effective use of data in small-area scale studies.
3) An active microwave remote sensing method comprises the following steps: soil moisture was inverted using the characteristic that the low frequency band (L, C band) radar backscatter coefficient was highly sensitive to soil moisture. Inversion algorithms can be categorized into 3 classes, empirical, theoretical, and semi-empirical. The empirical model is mostly to invert the soil humidity by establishing a regression relation between the multi-frequency multi-polarization multi-angle backscattering coefficient and the actual measured values such as soil humidity, soil roughness and the like of a certain depth, such as a classical Oh model and a Dubois model. The theoretical model is a forward model for modeling the surface backward radar scattering features. For bare earth surface, the theoretical model mainly comprises a Kirchhoff model and a small disturbance model (SPM), then an Integral Equation Model (IEM) for unifying the Kirchhoff model and the small disturbance model is developed, the real earth surface backscattering condition can be reproduced in a very wide earth surface roughness range, and the real earth surface backscattering condition is continuously improved and perfected to be developed into an Advanced Integral Equation Model (AIEM), so that the simulation result and the simulation precision of the model are continuously improved. The contribution of the vegetation layer to the microwave backscattering is an important factor influencing the inversion of soil humidity, and theoretical models for researching the microwave vegetation scattering characteristics, which are proposed by students, comprise a three-component vegetation scattering model, a first-order discrete vegetation model, a michigan microwave vegetation scattering model (MIMICS) and the like for vegetation covered earth surfaces. The theoretical model has more parameters, is difficult to directly use for inversion, and needs to construct a semi-empirical model on the basis. Aiming at the bare earth surface condition, the semi-empirical model generally utilizes IEM/AIEM to simulate earth surface backscattering characteristics under different earth surface roughness and soil humidity conditions, combines lookup table regression analysis, establishes correlation between backscattering coefficients and earth surface roughness and dielectric constant or soil humidity under different polarization combinations or incidence angle combinations, and further inverts the bare earth surface soil humidity, wherein the Shi model is particularly represented most. Aiming at vegetation covered earth surfaces, attema and the like, a semi-empirical model 'water-cloud' model for estimating soil moisture of the crop covered earth surfaces is provided, and is widely used; the MIMICS model is also commonly used to construct semi-empirical scattering models of vegetation cover surfaces, thereby inverting soil moisture. The active microwave remote sensing method can penetrate through the vegetation layer to detect soil with a depth of about 5cm below the ground surface, and has a strong inversion physical foundation; synthetic aperture radar (Synthetic Aperture Radar, SAR) can also provide high spatial resolution data on the order of meters to hundred meters; and a large number of comparison researches show that the accuracy of the active microwave algorithm for soil humidity inversion is higher than that of the optical algorithm and the passive microwave algorithm. The method for inverting the soil humidity by utilizing SAR data has the advantages of high spatial resolution, large detection depth and high inversion precision, has important application value and prospect in application of regional scale, and becomes one of important research directions in the field. However, in general, the soil moisture inversion method based on SAR has good inversion effect on bare or sparse vegetation ground surface, and the following two problems still exist on the high vegetation coverage ground surface:
firstly, the existing vegetation coverage inversion model depends on the coefficients of ground observation data regression experience and semi-experience inversion model, so that the practical application of the method in the vegetation coverage is limited;
secondly, in the existing vegetation coverage soil humidity inversion research, the vegetation is considered to cover the earth surface completely, and the change process of the earth surface from pure soil pixels to mixed pixels to pure vegetation pixels in the growth process from sowing to maturing of crops is omitted.
Disclosure of Invention
The invention provides a method for monitoring soil humidity in a crop growth process by combining SAR and optical remote sensing data, aiming at solving the problems.
The invention is realized in such a way that a method for monitoring soil humidity in a crop growth process by combining SAR and optical remote sensing data comprises the following steps:
s1, based on sowing time T 1 SAR satellite data SAR (T) 1 ,θ 1 ) SAR (T) 1 ,θ 2 ) Calculating a sowing time T 1 Soil surface humidity M of (2) v (T 1 );
S2, based on period T i SAR satellite data SAR (T) i θ), period T i+1 Optical remote sensing data OPT (T) i+1 ) SAR satellite data SAR (T) i+1 θ) determining the period T i And period T i+1 The relative change of the soil humidity in the time phase, and then the time T is obtained i+1 Soil surface humidity M of (2) v (T i+1 ) Wherein i is sequentially from 1 to h;
SAR(T 1 ,θ 1 ) Represents the sowing time T 1 Radar incident angle theta 1 SAR satellite data, θ 1 For small angles, SAR (T 1 ,θ 2 ) Represents the sowing time T 1 Incidence angle theta of radar 2 SAR satellite data, θ 2 Is at a large angle.
Further, soil surface humidity M v (T 1 ) The acquisition method of (a) is specifically as follows:
s11, SAR (T) based on SAR satellite data 1 ,θ 1 ) SAR (T) 1 ,θ 2 ) Acquiring backscattering coefficientAndWherein θ 1 And theta 2 Respectively representing a small angle and a large angle;
S14, obtaining S2Or->Combined roughness Z of soil s Calculating the sowing time T 1 Soil surface humidity M of (2) v (T 1 )。
Further, the combined roughness Z of bare soil s The calculation formula of (2) is specifically as follows:
wherein, c and d are fitting coefficients,represents the sowing time T 1 Incidence angle theta of radar 1 Soil meter under pq polarization modeA surface backscattering coefficient; />Represents the sowing time T 1 Incidence angle theta of radar 2 Backscattering coefficient in pq polarization mode.
Further, sowing time T 1 Soil surface humidity M of (2) v (T 1 ) The calculation formula is specifically as follows:
A pp (θ)、B pp (θ)、C pp (θ) is a cubic polynomial of the angle of incidence θ of the radar, Z s Indicating the combined roughness of bare soil.
Further, period T i+1 Soil surface humidity M of (2) v (T i+1 ) The acquisition method of (a) is specifically as follows:
s21, based on period T i+1 Optical remote sensing data OPT (T) i+1 ) Obtaining normalized vegetation index NDVI (T i+1 )、NDVI veg (T i+1 ) And NDVI soil (T i+1 );
S22, based on period T i+1 SAR satellite data SAR (Ti+1, theta) acquisition backscattering coefficientAnd
s23, based on NDVI (T) i+1 )、NDVI veg (T i+1 ) And NDVI soil (T i+1 ) Calculating vegetation coverage F (T) i+1 ) Further calculate the time period T i By time period T i+1 Average vegetation coverage of (a)
S24, willZ of combined roughness with bare soil s Calculate period T i And period T i+1 Relative change in phase soil humidity;
s25, period T i Soil surface humidity M of (2) v (T i ) Input period T i And period T i+1 Relative change of soil humidity in time phase, calculate time T i+1 Soil surface humidity M v (T i+1 )。
Further, period T i And period T i+1 The relative change formula of the soil humidity of the time phase is specifically as follows:
further, vegetation coverage F (T i ) The calculation formula is specifically as follows:
F(T i )=(NDVI(T i )-NDVI soil (T i ))/(NDVI veg (T i )-NDVI soil (T i ))
further, the crops are grown in dry farmland, including wheat, corn, rape, peanut, etc.
Further, a large angle θ 2 With a small angle theta 1 The difference in (2) is greater than 10.
According to the invention, through forward model simulation analysis, an inversion model of the humidity of the exposed soil surface is improved; based on a first-order discrete vegetation model, a soil humidity monitoring model based on combined SAR and optical data is constructed, and the problems of pathological inversion and mixed pixels of vegetation coverage soil humidity remote sensing inversion are solved; finally, the time series soil humidity monitoring of the crop growth process independent of ground observation data is realized.
Drawings
Fig. 1 is a flowchart of a method for monitoring soil humidity in a crop growth process by combining SAR and optical remote sensing data according to an embodiment of the present invention.
Detailed Description
The following detailed description of the embodiments of the invention, given by way of example only, is presented in the accompanying drawings to aid in a more complete, accurate, and thorough understanding of the inventive concepts and aspects of the invention by those skilled in the art.
According to the invention, through forward model simulation analysis, an inversion model of the humidity of the exposed soil surface is improved; based on a first-order discrete vegetation model, a soil humidity monitoring model based on combined SAR and optical data is constructed, and the problems of pathological inversion and mixed pixels of vegetation coverage soil humidity remote sensing inversion are solved; finally, the time series soil humidity monitoring of the crop growth process independent of ground observation data is realized, and the crops comprise crops grown in dry farmland, such as wheat, corn, rape, peanut and the like, and can be also suitable for grasslands
(1) Backscattering analog lookup table construction
And simulating the surface backscattering coefficients under different polarization modes, incidence angles, soil humidity and soil surface roughness (related length and root mean square height) of the C wave band by using an AIEM model, constructing a backscattering simulation lookup table, and providing data support for the subsequent combined roughness inversion algorithm improvement and bare soil surface humidity semi-empirical inversion model comparison analysis.
(2) Bare soil surface humidity inversion model improvement
Based on the constructed backscattering simulation lookup table, the inversion algorithm of the combined roughness is improved by combining the analysis of the bare soil surface humidity semi-empirical inversion model expression and utilizing a nonlinear regression analysis method, and then the bare soil surface humidity inversion model is improved. In the prior researches, a plurality of functional forms of fitting the combined roughness by utilizing the difference between the multi-angle backscattering coefficients are adopted, wherein the functional forms comprise polynomial function fitting and logarithmic function fitting, and the functions are respectively shown as a formula (1), a formula (2) and a formula (3):
wherein the subscript s represents soil (soil) surface scattering; z is Z s Representing the combined roughness of the soil surface, wherein the unit cm is the ratio of the square of the root mean square height to the relevant length; sigma (sigma) s The backscattering coefficient of the soil surface is expressed in dB; IS2 and IS7 respectively represent IS2 (19.2-26.7) and IS7 (41-. 5-45.2) incidence angle modes of ASAR; a, b are all coefficients obtained by fitting; the superscripts p, q denote the polarization, pq denotes one of hh, vv, hv or vh, pp denotes one of hh or vv, and therefore,represents the radar incidence angle of 38 DEG, the backscattering coefficient of the soil surface in hh polarization mode,>represents the angle of incidence of 22 DEG radar, the backscattering coefficient of the soil surface in hh polarization mode,>represents the angle of incidence of 38 DEG radar, the backscattering coefficient of the soil surface in the vv polarization mode,>represents the 22 DEG radar incidence angle, the backscattering coefficient of the soil surface in the vv polarization mode,>represents the soil surface backscattering coefficient in the pq polarization mode in IS2 incidence angle mode, < >>The soil surface backscattering coefficient in the pq polarization mode in the IS7 incidence angle mode IS shown.
The complex correlation coefficient of the three fitting expressions reaches more than 0.95, and the fitting effect is good. However, the polynomial models (1), (2) or the logarithmic model (3) have certain problems in practical application: for equations (1), (2), when the difference between the small angle and large angle backscattering coefficients is greater than 8.013dB and 6.272dB, the combined roughness will appear negative, and for equation (3) there will be a negative case. This negative occurrence is not uncommon in SAR images, and thus causes the occurrence of more pixels that fail inversion.
By observing the additive model expression (shown as formula (4)) of the soil surface humidity semi-empirical inversion model, the difference between the multi-angle backscattering coefficients and the combined roughness are found to be more in line with a form of an exponential function, as shown as formula (5):
wherein T is i Represents the growth period of wheat, when i=1, T i =T 1 The sowing time is represented as bare soil; m is M v (T i ) Indicating the time period T i The water content of the soil volume; θ is the angle of incidence of the radar, θ=θ 1 For small angles, θ=θ 2 Is a large angle, a large angle theta 2 With a small angle theta 1 The difference between (2) is more than 10 DEG;indicating the time period T i The soil surface backscattering coefficient under the radar incidence angle theta and pp polarization modes; />Represents the sowing time T 1 Incidence angle theta of radar 1 The soil surface backscattering coefficient under pq polarization mode; />Represents the sowing time T 1 Incidence angle theta of radar 2 Backscattering coefficient in pq polarization mode; c. d is the fitting coefficient, A pp (θ)、B pp (θ)、C pp (θ) is a cubic polynomial of the radar incidence angle θ, and the expression is specifically as follows:
A pp (θ)=a 1 sin 3 θ+b 1 sin 2 θ+c 1 sinθ+d 1 ;
B pp (θ)=a 2 sin 3 θ+b 2 sin 2 θ+c 2 sinθ+d 2 ;
C pp (θ)=a 3 sin 3 θ+b 3 sin 2 θ+c 3 sinθ+d 3 ;
in which a is 1 ,b 1 ,c 1 ,d 1 ,a 2 ,b 2 ,c 2 ,d 2 ,a 3 ,b 3 ,c 3 ,d 3 Is a fitting coefficient.
Fitting a formula (5) by using an AIEM simulation lookup table according to an exponential function relation, wherein under the condition of the incidence angle difference between IS2 and IS7, the fitted complex correlation coefficient reaches 0.996, which IS superior to the polynomial model and the logarithmic model; and the situation that the combined roughness has a negative value can not be generated, the actual situation is more met, and inversion failure is avoided. Therefore, the above method would be more advantageous in soil roughness inversion.
As crops sow and grow, the vegetation coverage of the soil surface gradually changes, and the problem of mixed pixels is brought. The first-order discrete vegetation model is a back scattering conceptual model of mixed pixels of different vegetation coverage earth surfaces, and can not directly realize soil humidity inversion of the vegetation coverage earth surfaces, but can realize relative change monitoring of soil humidity on the premise of acquiring vegetation coverage information and eliminating vegetation double-pass attenuation influence: sub-pel decomposition is performed by using the optical remote sensing data and a pel bipartite model, and pel vegetation coverage information of multiple phases is extracted; the first-order discrete vegetation model and the three-component scattering model are utilized to separate out the surface scattering components of the soil surface, and the influence of a vegetation double-pass attenuation factor is eliminated through the ratio of the multi-time phase surface scattering components; and finally, constructing a multi-time-phase soil humidity relative change monitoring model independent of ground observation data by combining the semi-empirical inversion model of the surface humidity of the bare soil and the multi-time-phase multi-polarization SAR data.
(3) Separating surface scattering components of first-order discrete vegetation model
The first-order discrete vegetation model is a mixed pixel back scattering conceptual model considering different vegetation coverage, and the model is expressed as follows:
wherein F (T) i ) For period T i And θ is the incidence angle of the radar, pp represents v or hh polarization mode, subscripts t, v, s, sv represent total scattering, vegetation body scattering, soil surface scattering, vegetation body and soil surface body-surface interaction scattering items respectively, therefore,indicating the time period T i Back volume scattering coefficient component of vegetation in radar incidence angle theta, pp polarization mode, +.>Indicating the time period T i Backscattering coefficient component of vegetation and soil body-surface interaction under radar incidence angle theta and pp polarization modes, ++>Indicating the time period T i The backward scattering coefficient component of the soil under the radar incidence angle theta and pp polarization mode,indicating the time period T i Total backscattering coefficient in radar incidence angle θ, pp polarization, +.>Indicating the time period T i Is a vegetation double pass attenuation factor.
In the co-polarized pp (hh and vv polarization) state, when the radar is incident at an angle θ>At 20 °, the body-surface cross scatter in the first order discrete vegetation modelNegligible; at the same time, assuming that the crop cover layer is composed of randomly distributed fine cylinders, according to Freeman A three-component scattering model, the bulk scattering of the vegetation layer can be considered to be unaffected by polarization, about 3 times the total scattering of cross polarization, i.e.)>Thus, the surface scattering component->Can be separated into:
(4) Calculating pixel vegetation coverage by optical data pixel bipartite model
In the growth process of crops, the vegetation coverage is continuously changed, and the pixels are also changed from pure bare soil pixels to mixed pixels of soil and vegetation, and then to pure vegetation pixels. The pel bipartite model is the most commonly used mixed pel linear decomposition model, and pel vegetation coverage can be inverted by using the pel bipartite model. The pel bipartite model can be described as:
S t (T i )=F(T i )·S veg (T i )+(1-F(T i ))·S soil (T i ) (8)
wherein F (T) i ) For period T i Is a pel vegetation coverage of (1); s is S t (T i ) Indicating the time period T i Is the total information of the pixels; s is S veg (T i ) Indicating the time period T i Remote sensing information of pure vegetation pixels; s is S soil (T i ) Indicating the time period T i Remote sensing information of pure soil pixels of (2), vegetation coverage of pixels F (T) i ) Can be expressed as formula (9):
F(T i )=(S t (T i )-S soil (T i ))/(S veg (T i )-S soil (T i )) (9)
substituting NDVI as remote sensing information into formula (9) to obtain vegetation coverage F (T) i ) A calculation formula, such as formula (10):
F(T i )=(NDVI(T i )-NDVI soil (T i ))/(NDVI veg (T i )-NDVI soil (T i )) (10)
in the formula, NDVI (T) i ) For period T i An actual NDVI value for the pel of (c); NDVI veg (T i ) For period T i NDVI value of the pure vegetation picture element, NDVI soil (T i ) For period T i NDVI values of bare soil pixels without vegetation.
(5) Construction of soil humidity relative change monitoring model in crop growth process
According to the formula (7), the vegetation coverage F alone cannot calculate the soil surface scattering component, and the vegetation double-pass attenuation factor L needs to be further eliminated 2 pp Is a function of (a) and (b). Considering that the vegetation cover layer changes have longer time scales than the soil humidity changes in the crop growth process, two time phases T with similar time are considered i And T i+1 On the vegetation coverage F and vegetation double-pass attenuation factor L 2 pp The variation is small, namely:
indicating the time period T i By time period T i+1 The average vegetation coverage of (2) can be calculated by comparing the scattering components of soil surface to eliminate L 2 pp A factor as shown in formula (13):
the right side of the equation can be calculated by remote sensing data; the left side of the equation is the ratio of the scattering components of the soil surface, and the substitution of the formula (4) into the equation can be obtained:
and then combining the formula (5) to obtain T i And T i+1 Relative change in phase soil humidity:
[M v (T i+1 )] m /[M v (T i )] n =k (15)
wherein m, n, k can be obtained by simplifying the formula (14).
Fig. 1 is a flowchart of a method for monitoring soil humidity in a crop growth process by combining SAR and optical remote sensing data, provided by an embodiment of the present invention, and the method specifically includes:
s1, SAR (T) based on SAR satellite data 1 ,θ 1 ) SAR (T) 1 ,θ 2 ) Acquiring backscattering coefficientAndWherein θ 1 And theta 2 Respectively represent a small angle and a large angleA degree;
SAR(T 1 ,θ 1 ) Represents the sowing time T 1 Radar incident angle theta 1 SAR satellite data, θ 1 For small angles, SAR (T 1 ,θ 2 ) Represents the sowing time T 1 Incidence angle theta of radar 2 SAR satellite data, θ 2 Is a large angle;represents the sowing time T 1 Incidence angle theta of radar 1 Total backscattering coefficient in pq polarization,>represents the sowing time T 1 Incidence angle theta of radar 2 Total backscattering coefficient in pq polarization mode.
S3, willIs->Substituting the rough soil into the formula (5) to calculate the combined roughness Z of the bare soil s ;
S4, obtaining the product in the step S2Or->Combined roughness Z of soil s Substituting formula (4) to calculate seeding time T 1 Is wet on the soil surfaceDegree M v (T 1 );
S5, based on T 2 Time-phased optical remote sensing data OPT (T 2 ) Obtaining normalized vegetation index NDVI (T) 2 )、NDVI veg (T 2 ) And NDVI soil (T 2 ),T 2 The period is a sowing period T 1 Is the next period of time;
S7, NDVI (T) 2 )、NDVI veg (T 2 ) And NDVI soil (T 2 ) Substituting formula (10) to calculate vegetation coverage F (T) 2 ) F (T) 1 ) And F (T) 2 ) Substituting into formula (11) to obtain
S7, willAnd Z obtained in the step S3 s Substituting formula (14) and simplifying to obtain T 1 And T 2 Phase soil humidity relative change formula (15);
s8, combining the formula (15) and M obtained in the step S4 v (T 1 ) Calculating to obtain M v (T 2 )。
And so on, any period T of the growth process i+1 Soil surface humidity inversion of (2):
s9, based on optical remote sensing data OPT (T) i+1 ) Obtaining normalized vegetation index NDVI (T) i+1 )、NDVI veg (T i+1 ) And NDVI soil (T i+1 ) The method comprises the steps of carrying out a first treatment on the surface of the SAR (T) based i+1 θ) acquisition of backscattering coefficientAnd->
S10, NDVI (T) obtained in the step S9 i+1 )、NDVI veg (T i+1 ) And NDVI soil (T i+1 ) Substituting formula (10) to calculate vegetation coverage F (T) i+1 ) F (T) i ) And F (T) i+1 ) Substituting into formula (11) to obtain
S11, willAnd Z obtained by S3 s Is substituted into formula (14) to be simplified to obtain T i And T i+1 Phase soil humidity relative change formula (15);
s12, combining the formula (15) and the last approaching period T i M of (2) v (T i ) Calculating to obtain M v (T i+1 )。
Winter wheat growth process time sequence soil humidity monitoring
Taking a typical planting area of winter wheat in Mongolian county of Anhui province as a research experiment area, taking a sowing period in the growth process of winter wheat as an initial time phase, and directly calculating by using a formula (4) and a formula (5) according to the prior knowledge of the soil exposure of the sowing period to obtain the soil surface combined roughness Z s Initial soil humidity M v (T 1 ) Then substituting formula (15) to obtain M v (T 2 ). Similarly, M can be calculated v (T 3 ),M v (T 4 ) … …, thereby obtaining time sequence soil humidity and realizing the monitoring of the soil humidity in the winter wheat growth process.
While the present invention has been described by way of example, it should be apparent that the practice of the invention is not limited by the foregoing, but rather is intended to cover various insubstantial modifications of the method concepts and teachings of the invention, either as applied to other applications without modification, or as applied directly to other applications, without departing from the scope of the invention.
Claims (3)
1. The method for monitoring the soil humidity in the crop growth process by combining SAR and optical remote sensing data is characterized by comprising the following steps of:
s1, based on sowing time T 1 SAR satellite data SAR (T) 1 ,θ 1 ) SAR (T) 1 ,θ 2 ) Calculating a sowing time T 1 Soil surface humidity M of (2) v (T 1 );
S2, based on period T i SAR satellite data SAR (T) i θ), period T i+1 Optical remote sensing data OPT (T) i+1 ) SAR satellite data SAR (T) i+1 θ) determining the period T i And period T i+1 The relative change of the soil humidity in the time phase, and then the time T is obtained i+1 Soil surface humidity M of (2) v (T i+1 ) Wherein i is sequentially from 1 to h;
SAR(T 1 ,θ 1 ) Represents the sowing time T 1 Radar incident angle theta 1 SAR satellite data, θ 1 For small angles, SAR (T 1 ,θ 2 ) Represents the sowing time T 1 Incidence angle theta of radar 2 SAR satellite data, θ 2 Is a large angle;
soil surface humidity M v (T 1 ) The acquisition method of (a) is specifically as follows:
s11, SAR (T) based on SAR satellite data 1 ,θ 1 ) SAR (T) 1 ,θ 2 ) Acquiring backscattering coefficientAndWherein θ 1 And theta 2 Respectively representing a small angle and a large angle;
S14, obtaining S2Or->Combined roughness Z of soil s Calculating the sowing time T 1 Soil surface humidity M of (2) v (T 1 );
Combined roughness Z of bare soil s The calculation formula of (2) is specifically as follows:
wherein, c and d are fitting coefficients,represents the sowing time T 1 Incidence angle theta of radar 1 The soil surface backscattering coefficient under pq polarization mode; />Represents the sowing time T 1 Incidence angle theta of radar 2 Under pq polarization modeIs a backscatter coefficient of (2);
sowing time T 1 Soil surface humidity M of (2) v (T 1 ) The calculation formula is specifically as follows:
A pp (θ)、B pp (θ)、C pp (θ) is a cubic polynomial of the angle of incidence θ of the radar, Z s Representing the combined roughness of bare soil;
period T i+1 Soil surface humidity M of (2) v (T i+1 ) The acquisition method of (a) is specifically as follows:
s21, based on period T i+1 Optical remote sensing data OPT (T) i+1 ) Obtaining normalized vegetation index NDVI (T i+1 )、NDVI veg (T i+1 ) And NDVI soil (T i+1 );
S22, based on period T i+1 SAR satellite data SAR (T) i+1 θ) acquisition of backscattering coefficientAnd
s23, based on NDVI (T) i+1 )、NDVI veg (T i+1 ) And NDVI soil (T i+1 ) Calculating vegetation coverage F (T) i+1 ) Further calculate the time period T i By time period T i+1 Average vegetation coverage of (a)
S24, willZ of combined roughness with bare soil s Calculate period T i And period T i+1 Relative change in phase soil humidity; />
S25, period T i Soil surface humidity M of (2) v (T i ) Input period T i And period T i+1 Relative change of soil humidity in time phase, calculate time T i+1 Soil surface humidity M v (T i+1 );
Period T i And period T i+1 The relative change formula of the soil humidity of the time phase is specifically as follows:
vegetation coverage F (T) i ) The calculation formula is specifically as follows:
F(T i )=(NDVI(T i )-NDVI soil (T i ))/(NDVI veg (T i )-NDVI soil (T i ))。
2. the method for monitoring soil moisture during crop growth combining SAR and optical remote sensing data as recited in claim 1, wherein said crop is grown in dry farmland including wheat, corn, canola, peanut.
3. The method for monitoring soil moisture during crop growth combining SAR and optical remote sensing data as recited in claim 1, wherein the large angle θ 2 With a small angle theta 1 The difference in (2) is greater than 10.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211111362.9A CN115479960B (en) | 2022-09-13 | 2022-09-13 | Method for monitoring soil humidity in crop growth process by combining SAR (synthetic aperture radar) with optical remote sensing data |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211111362.9A CN115479960B (en) | 2022-09-13 | 2022-09-13 | Method for monitoring soil humidity in crop growth process by combining SAR (synthetic aperture radar) with optical remote sensing data |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115479960A CN115479960A (en) | 2022-12-16 |
CN115479960B true CN115479960B (en) | 2023-06-13 |
Family
ID=84393060
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211111362.9A Active CN115479960B (en) | 2022-09-13 | 2022-09-13 | Method for monitoring soil humidity in crop growth process by combining SAR (synthetic aperture radar) with optical remote sensing data |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115479960B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116049768B (en) * | 2023-04-03 | 2023-07-11 | 中国科学院空天信息创新研究院 | Active and passive microwave soil moisture fusion algorithm based on generalized linear regression model |
CN117574161B (en) * | 2024-01-17 | 2024-04-16 | 航天宏图信息技术股份有限公司 | Surface parameter estimation method, device and equipment based on generation of countermeasure network |
CN117826112B (en) * | 2024-03-05 | 2024-05-31 | 天津智云水务科技有限公司 | Soil water content inversion method based on sar |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112418016B (en) * | 2020-11-09 | 2024-07-12 | 中国农业大学 | SAR-based irrigation information extraction method and device |
EP4050334A1 (en) * | 2021-02-26 | 2022-08-31 | Tata Consultancy Services Limited | System and method for root zone soil moisture estimation for vegetation cover using remote sensing |
-
2022
- 2022-09-13 CN CN202211111362.9A patent/CN115479960B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN115479960A (en) | 2022-12-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115479960B (en) | Method for monitoring soil humidity in crop growth process by combining SAR (synthetic aperture radar) with optical remote sensing data | |
Jeong et al. | Predicting rice yield at pixel scale through synthetic use of crop and deep learning models with satellite data in South and North Korea | |
Betbeder et al. | Assimilation of LAI and dry biomass data from optical and SAR images into an agro-meteorological model to estimate soybean yield | |
Rahimzadegan et al. | Estimating evapotranspiration of pistachio crop based on SEBAL algorithm using Landsat 8 satellite imagery | |
CN108509836B (en) | Crop yield estimation method based on double-polarized synthetic aperture radar and crop model data assimilation | |
Kurum | Quantifying scattering albedo in microwave emission of vegetated terrain | |
CN108802728B (en) | Crop irrigation guiding method for assimilation of bipolar synthetic aperture radar and crop model | |
CN111751286B (en) | Soil moisture extraction method based on change detection algorithm | |
Yang et al. | In-season biomass estimation of oilseed rape (Brassica napus L.) using fully polarimetric SAR imagery | |
WO2018107245A1 (en) | Detection of environmental conditions | |
Togliatti et al. | Satellite L–band vegetation optical depth is directly proportional to crop water in the US Corn Belt | |
CN114065643B (en) | SAR and polarization decomposition-based plant soil water content estimation method and system | |
CN115526098B (en) | Remote sensing calculation method for leaf area index of surface vegetation in mining area and electronic equipment | |
Kong et al. | A coupling model for soil moisture retrieval in sparse vegetation covered areas based on microwave and optical remote sensing data | |
Luo et al. | Crop yield estimation based on assimilation of crop models and remote sensing data: A systematic evaluation | |
Nourani et al. | Hydrological model parameterization using NDVI values to account for the effects of land cover change on the rainfall–runoff response | |
Ouaadi et al. | C-band radar data and in situ measurements for the monitoring of wheat crops in a semi-arid area (center of Morocco) | |
Singh et al. | Incorporation of first-order backscattered power in Water Cloud Model for improving the Leaf Area Index and Soil Moisture retrieval using dual-polarized Sentinel-1 SAR data | |
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 | |
Das et al. | Ensemble surface soil moisture estimates at farm-scale combining satellite-based optical-thermal-microwave remote sensing observations | |
MAO et al. | The monitoring analysis for the drought in China by using an improved MPI method | |
Nasirzadehdizaji et al. | Application of sentinel-1 multi-temporal data for crop monitoring and mapping | |
CN114611699A (en) | Soil moisture downscaling method and device, electronic equipment and storage medium | |
Silvestro et al. | Development of an assimilation scheme for the estimation of drought-induced yield losses based on multi-source remote sensing and the AcquaCrop model | |
Humes et al. | High spatial resolution mapping of surface energy balance components with remotely sensed data |
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 |