CN112014323B - Remote sensing inversion method and system for surface evapotranspiration - Google Patents
Remote sensing inversion method and system for surface evapotranspiration Download PDFInfo
- Publication number
- CN112014323B CN112014323B CN202010895830.0A CN202010895830A CN112014323B CN 112014323 B CN112014323 B CN 112014323B CN 202010895830 A CN202010895830 A CN 202010895830A CN 112014323 B CN112014323 B CN 112014323B
- Authority
- CN
- China
- Prior art keywords
- soil
- temperature
- vegetation
- relative
- calculating
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
- G01W1/02—Instruments for indicating weather conditions by measuring two or more variables, e.g. humidity, pressure, temperature, cloud cover or wind speed
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N2021/1793—Remote sensing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N2021/1793—Remote sensing
- G01N2021/1795—Atmospheric mapping of gases
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Atmospheric Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Ecology (AREA)
- Environmental Sciences (AREA)
- Radiation Pyrometers (AREA)
Abstract
The invention discloses a remote sensing inversion method of surface evapotranspiration, which comprises the following steps: acquiring meteorological data and remote sensing data of the earth surface; respectively calculating the soil temperature and the vegetation temperature relative to the soil moisture of 0 and the soil temperature and the vegetation temperature relative to the soil moisture of 1 according to the meteorological data and the remote sensing data; further calculating the soil dry edge temperature, the soil humidity limiting factor, the vegetation humidity limiting factor and the vegetation temperature limiting factor; and calculating the surface evaporation capacity by using a Priestley-Taylor formula according to the soil humidity limiting factor, the vegetation humidity limiting factor and the temperature limiting factor. According to the invention, the calculation limiting factors of the soil temperature and the vegetation temperature relative to the soil moisture of 0 and the calculation limiting factors of the soil temperature and the vegetation temperature relative to the soil moisture of 1 are introduced, so that the water change of the soil surface layer and the vegetation root area on the daily scale can be reflected, and the accuracy of remote sensing inversion of the earth surface evapotranspiration is improved by adopting the inversion method of daily scale earth surface evapotranspiration.
Description
Technical Field
The invention relates to the technical field of surface evapotranspiration calculation, in particular to a remote sensing inversion method and system for surface evapotranspiration.
Background
The surface evapotranspiration is an important component of surface energy balance and water circulation, and is also an important index of vegetation growth condition and crop yield. The method has the advantages that regional or global-scale surface evapotranspiration can be accurately estimated, and the method has important significance for global climate evolution research, water resource planning and management, agricultural water conservation and crop yield estimation. However, the evapotranspiration has strong spatial variability, and the remote sensing technology has the characteristics of high speed, high space-time resolution and suitability for large-area long-term observation, and is considered as the most effective means for acquiring the evapotranspiration distribution of the earth surface area or the global scale with high time efficiency and high precision.
Due to the characteristics of instantaneous transit of the remote sensing satellite, the earth surface evapotranspiration obtained based on the thermal infrared remote sensing inversion model is generally an instantaneous value of the transit time of the satellite, and the change rule of the earth surface evapotranspiration in a daily scale and a longer time scale cannot be reflected. In practical application, the daily scale and longer time scale evapotranspiration are more valuable. The ratio of variables closely related to evapotranspiration under clear sky conditions to the evapotranspiration is always assumed to be basically unchanged in the daytime, and the time scale expansion is carried out on the instantaneous evapotranspiration obtained by the inversion of the satellite transit time to obtain the daily scale evapotranspiration. However, the daily scale expansion method is based on the instantaneous evapotranspiration of remote sensing inversion, and the accuracy and the application condition of different remote sensing inversion methods are different. Therefore, the estimation result of the earth surface evapotranspiration of the current day scale mostly has two uncertainties: (1) uncertainty of the remote sensing inversion model; (2) uncertainty of the daily scale extension method. In order to avoid the error accumulation of uncertainty of the two methods, a daily scale earth surface evapotranspiration remote sensing inversion model is directly constructed, and the method has important research significance for improving the accuracy of the daily scale earth surface evapotranspiration of remote sensing inversion.
Disclosure of Invention
The invention aims to provide a remote sensing inversion method and a remote sensing inversion system for ground evapotranspiration, so as to avoid error accumulation of uncertainty of a remote sensing inversion model and uncertainty of a daily scale expansion method.
In order to achieve the purpose, the invention provides the following scheme:
a remote sensing inversion method of ground surface evapotranspiration comprises the following steps:
acquiring meteorological data and remote sensing data of the earth surface; the meteorological data comprise air temperature, atmospheric pressure, air relative humidity, wind speed, uplink short-wave radiation, downlink short-wave radiation and the like; the remote sensing data comprises earth surface temperature, earth surface emissivity, reflectivity, leaf area index, enhanced vegetation index, photosynthetically active radiation and the like;
respectively calculating the soil temperature and the vegetation temperature relative to the soil moisture of 0 and the soil temperature and the vegetation temperature relative to the soil moisture of 1 according to the meteorological data and the remote sensing data of the earth surface;
calculating and distinguishing the soil dry edge temperature of the sufficient water state of the root zone and the water-free state of the soil surface layer according to the soil temperature of 0 relative to the soil moisture, the vegetation temperature of 1 relative to the soil moisture and remote sensing data;
calculating a soil humidity limiting factor, a vegetation humidity limiting factor and a vegetation temperature limiting factor according to the soil temperature and the vegetation temperature relative to the soil moisture of 0, the soil temperature and the vegetation temperature relative to the soil moisture of 1 and the soil dry edge temperature;
and calculating the surface evaporation capacity by using a Priestley-Taylor (Priesterli-Taylor) formula according to the soil humidity limiting factor, the vegetation humidity limiting factor and the vegetation temperature limiting factor.
Optionally, according to meteorological data and the remote sensing data of earth's surface respectively calculate soil temperature and vegetation temperature that relative soil moisture is 0, soil temperature and vegetation temperature that relative soil moisture is 1, specifically include:
using a formula based on the meteorological data and the remote sensing dataCalculating the soil temperature relative to the soil moisture of 0;
using a formula based on the meteorological data and the remote sensing dataCalculating the vegetation temperature relative to the soil moisture of 0;
using a formula based on the meteorological data and the remote sensing dataCalculating the soil temperature relative to the soil moisture of 1;
using a formula based on the meteorological data and the remote sensing dataCalculating the vegetation temperature relative to the soil moisture of 1;
wherein, Ts0And Tv0Respectively representing the soil temperature and vegetation temperature, T, relative to the soil moisture of 0s1And Tv1Respectively representing the soil temperature and the vegetation temperature relative to the soil moisture of 1; ρ is the air density; cpIs the specific heat at constant pressure; e.g. of the typesSaturated water vapor pressure; e.g. of the typeaThe actual water vapor pressure of the atmosphere; t isaThe temperature near the surface is obtained; gvwVegetation canopy conductivity sufficient for water supply; gavAnd gasThe aerodynamic conductivity of the vegetation and the upper layer of the soil are respectively; rn,sAnd Rn,vRespectively the soil component net radiation and the vegetation component net radiation; gsIs the soil heat flux.
Optionally, the soil dry edge temperature for distinguishing the sufficient water state of the root zone and the moisture-free state of the soil surface layer is calculated according to the soil temperature relative to the soil moisture of 0 and the vegetation temperature relative to the soil moisture of 1 and the remote sensing data, and specifically includes:
Wherein EVI is Enhanced Vegetation Index (EVI)minAnd EVImaxRespectively the minimum enhancement type vegetation index corresponding to bare soil and the maximum enhancement type vegetation index corresponding to full vegetation coverage;
according to the soil temperature with the relative soil moisture of 0 and the vegetation temperature with the relative soil moisture of 1 and the vegetation proportion, utilizing a formulaCalculating and distinguishing soil dry edge temperature T of sufficient water state and no water state on soil surface layer of root zone*;
Wherein, Ts0And Tv1The soil temperature is 0 relative to the soil moisture and the vegetation temperature is 1 relative to the soil moisture.
Optionally, the soil temperature and vegetation temperature that is 0 according to relative soil moisture, relative soil moisture and the soil temperature and vegetation temperature that is 1 and soil dry edge temperature calculate soil humidity limiting factor, vegetation humidity limiting factor and vegetation temperature limiting factor, specifically include:
judging whether the surface temperature of the remote sensing image pixel is greater than the dry edge temperature of the soil or not to obtain a judgment result;
if the judgment result shows no, the formula f is usedVMCalculating a vegetation moisture limiting factor f 1VM(ii) a Using formulasCalculating the soil humidity limiting factor fSM;
If the judgment result shows yes, the formula f is utilizedSMCalculating the soil humidity limiting factor f as 0SM(ii) a Using formulasCalculating vegetation humidity limiting factor fVM;
Wherein, Ts0And Ts1Respectively, a soil temperature of 0 relative to the soil moisture and a soil temperature of 1 relative to the soil moisture, TsoilWhich is indicative of the temperature of the soil components,TRfor remote-sensing the surface temperature of image elements, FvIs the vegetation proportion; t isv0And Tv1Respectively vegetation temperature of 0 relative to the soil moisture and vegetation temperature of 1 relative to the soil moisture, TvegIs the temperature of the components of the vegetation,
Wherein, TmaxThe maximum atmospheric temperature in one day; t isoptThe temperature is the optimum temperature for vegetation growth.
Optionally, calculating the surface evaporation capacity by using a Priestley-Taylor formula according to the soil humidity limiting factor, the vegetation humidity limiting factor and the vegetation temperature limiting factor specifically includes:
according to the vegetation temperature limiting factor and the vegetation humidity limiting factor, using a formulaCalculating the transpiration amount of the vegetation;
according to the soil humidity limiting factor, using the formulaCalculating the evaporation capacity of the soil;
according to the vegetation transpiration amount, the soil evaporation amount and the cut-off evaporation amount, utilizing a formula ET ═ ETv+ETs+ETiCalculating the evaporation capacity of the earth surface;
wherein, ETv、ETsAnd ETiThe vegetation transpiration amount, the soil evaporation amount and the closure evaporation amount are respectively, and ET is the earth surface evaporation amount; f. ofwetIs the proportion of the wet leaves of the vegetation,RH is the relative humidity of the atmosphere; f. ofgThe ratio of vegetation to green leaves;fAPARphotosynthetically active radiation absorbed for green plants; f. ofIPARPhotosynthetically active radiation intercepted for green vegetation; f. ofTIs a vegetation temperature limiting factor; f. ofSMIs a soil humidity limiting factor; f. ofVMIs vegetation humidity limiting factor; alpha is the coefficient of a Priestley-Taylor formula, and alpha is 1.26; e.g. of the typesSaturated water vapor pressure; t isaThe temperature near the surface is obtained; rn,vFor canopy net radiation, Rn,sNet radiation for soil; gsIs the soil heat flux.
A remote sensing inversion system of surface evapotranspiration, the inversion system comprising, for example:
the data acquisition module is used for acquiring meteorological data and remote sensing data of the earth surface; the meteorological data comprise air temperature, atmospheric pressure, air relative humidity, wind speed, uplink short-wave radiation, downlink short-wave radiation and the like; the remote sensing data comprises earth surface temperature, earth surface emissivity, reflectivity, leaf area index, enhanced vegetation index, photosynthetically active radiation and the like;
the soil temperature and vegetation temperature calculation module is used for calculating the soil temperature and vegetation temperature relative to the soil moisture of 0 and the soil temperature and vegetation temperature relative to the soil moisture of 1 respectively according to the meteorological data and the remote sensing data of the earth surface;
the soil dry edge temperature calculation module is used for calculating and distinguishing the soil dry edge temperature of the sufficient water state of the root area and the no water state of the soil surface layer according to the soil temperature relative to the soil water content of 0, the vegetation temperature relative to the soil water content of 1 and remote sensing data;
the limiting factor calculation module is used for calculating a soil humidity limiting factor, a vegetation humidity limiting factor and a vegetation temperature limiting factor according to the soil temperature and the vegetation temperature relative to the soil moisture of 0, the soil temperature and the vegetation temperature relative to the soil moisture of 1 and the soil dry edge temperature;
and the earth surface evaporation amount calculation module is used for calculating the earth surface evaporation amount by utilizing a Priestley-Taylor formula according to the soil humidity limiting factor, the vegetation humidity limiting factor and the temperature limiting factor.
Optionally, the soil temperature and vegetation temperature calculation module specifically includes:
a soil temperature calculation submodule with a relative soil moisture of 0 for utilizing a formula according to the meteorological data and the remote sensing dataCalculating the soil temperature relative to the soil moisture of 0;
a vegetation temperature calculation submodule with a relative soil moisture of 0 for utilizing public service according to the meteorological data and the remote sensing dataFormula (II)Calculating the vegetation temperature relative to the soil moisture of 0;
a soil temperature calculation submodule with a relative soil moisture of 1 for utilizing a formula according to the meteorological data and the remote sensing dataCalculating the soil temperature relative to the soil moisture of 1;
a vegetation temperature calculation submodule with relative soil moisture of 1 for utilizing a formula according to the meteorological data and the remote sensing dataCalculating the vegetation temperature relative to the soil moisture of 1;
wherein, Ts0And Tv0Respectively representing the soil temperature and vegetation temperature, T, relative to the soil moisture of 0s1And Tv1Respectively representing the soil temperature and the vegetation temperature relative to the soil moisture of 1; ρ is the air density; cpIs the specific heat at constant pressure; e.g. of the typesSaturated water vapor pressure; e.g. of the typeaThe actual water vapor pressure of the atmosphere; t isaThe temperature near the surface is obtained; gvwVegetation canopy conductivity sufficient for water supply; gavAnd gasThe aerodynamic conductivity of the vegetation and the upper layer of the soil are respectively; rn,sAnd Rn,vRespectively the soil component net radiation and the vegetation component net radiation; gsIs the soil heat flux.
Optionally, the soil dry edge temperature calculation module specifically includes:
a vegetation coverage calculation submodule for utilizing a formula according to the earth surface remote sensing dataCalculating the vegetation ratio Fv;
Wherein EVI is Enhanced Vegetation Index (EVI)minAnd EVImaxRespectively corresponding to the most bare soilA small enhanced vegetation index and a maximum enhanced vegetation index corresponding to a full vegetation coverage;
a soil dry edge temperature operator module for utilizing a formula according to the soil temperature of 0 relative to the soil moisture and the vegetation temperature and vegetation proportion of 1 relative to the soil moistureCalculating and distinguishing soil dry edge temperature T of sufficient water state and no water state on soil surface layer of root zone*;
Wherein, Ts0And Tv1The soil temperature is 0 relative to the soil moisture and the vegetation temperature is 1 relative to the soil moisture.
Optionally, the limiting factor calculating module specifically includes:
the judgment submodule is used for judging whether the surface temperature of the remote sensing image pixel is greater than the dry edge temperature of the soil or not to obtain a judgment result;
a first limiting factor calculating submodule for utilizing the formula f if the judgment result shows noVMCalculating vegetation humidity limiting factor f 1VM(ii) a Using formulasCalculating the soil humidity limiting factor fSM;
A second limiting factor calculating submodule for using the formula f if the judgment result shows yesSMCalculating a vegetation humidity limiting factor fVM(ii) a Using formulasCalculating vegetation humidity limiting factor fVM;
Wherein, Ts0And Ts1The soil temperature is 0 relative to the soil moisture and the soil temperature is 1 relative to the soil moisture; t issoilWhich is indicative of the temperature of the soil components,TRas remote sensing imagesSurface temperature of the picture element, FvIs the vegetation proportion; t isv0And Tv1Respectively vegetation temperature of 0 relative to the soil moisture and vegetation temperature of 1 relative to the soil moisture, TvegIs the temperature of the components of the vegetation,
a third limiting factor calculation submodule for utilizing a formulaCalculating a vegetation temperature limiting factor fT;
Wherein, TmaxThe maximum atmospheric temperature in one day; t isoptThe temperature is the optimum temperature for vegetation growth.
Optionally, the surface evaporation amount calculation module specifically includes:
a vegetation transpiration meter operator module for utilizing a formula according to the vegetation temperature limiting factor and the vegetation humidity limiting factorCalculating the transpiration amount of the vegetation;
a soil evaporation meter operator module for utilizing a formula according to the soil humidity limiting factorCalculating the evaporation capacity of the soil;
a cut-off evaporation amount calculation operator module for utilizingCalculating the interception evaporation capacity;
the earth surface evapotranspiration meter operator module is used for utilizing a formula ET-ET according to the vegetation evapotranspiration, the soil evaporation and the interception evaporationv+ETs+ETiCalculating the evapotranspiration amount of the earth surface;
wherein, ETv、ETsAnd ETiThe vegetation transpiration amount, the soil evaporation amount and the cut-off evaporation amount are respectivelyQuantity, ET is the surface evaporation capacity; f. ofwetIs the proportion of the wet leaves of the vegetation,RH is the relative humidity of the atmosphere; f. ofgThe ratio of vegetation to green leaves;fAPARphotosynthetically active radiation absorbed for green plants; f. ofIPARPhotosynthetically active radiation intercepted for green vegetation; f. ofTIs a vegetation temperature limiting factor; f. ofSMIs a soil humidity limiting factor; f. ofVMIs vegetation humidity limiting factor; alpha is the coefficient of a Priestley-Taylor formula, and alpha is 1.26; e.g. of the typesSaturated water vapor pressure; t isaThe temperature near the surface is obtained; rn,vFor canopy net radiation, Rn,sNet radiation for soil; gsIs the soil heat flux.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
the invention discloses a remote sensing inversion method of surface evapotranspiration, which comprises the following steps: acquiring meteorological data and remote sensing data of the earth surface; respectively calculating the soil temperature and the vegetation temperature relative to the soil moisture of 0 and the soil temperature and the vegetation temperature relative to the soil moisture of 1 according to the meteorological data and the remote sensing data of the earth surface; calculating and distinguishing the soil dry edge temperature of the sufficient water state of the root zone and the water-free state of the soil surface layer according to the soil temperature of 0 relative to the soil moisture, the vegetation temperature of 1 relative to the soil moisture and the remote sensing data of the earth surface; calculating a soil humidity limiting factor, a vegetation humidity limiting factor and a vegetation temperature limiting factor according to the soil temperature and the vegetation temperature relative to the soil moisture of 0, the soil temperature and the vegetation temperature relative to the soil moisture of 1 and the soil dry edge temperature; and calculating the surface evaporation capacity by using a Priestley-Taylor formula according to the soil humidity limiting factor, the vegetation humidity limiting factor and the temperature limiting factor. According to the invention, the calculation limiting factors of the soil temperature and the vegetation temperature relative to the soil moisture of 0 and the calculation limiting factors of the soil temperature and the vegetation temperature relative to the soil moisture of 1 are introduced, so that the water change of the soil surface layer and the vegetation root area on the daily scale can be reflected, and the accuracy of remote sensing inversion of the earth surface evapotranspiration is improved by adopting the inversion method of daily scale earth surface evapotranspiration.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without inventive exercise.
Fig. 1 is a flowchart of a remote sensing inversion method of surface evapotranspiration provided by the present invention.
Detailed Description
The invention aims to provide a remote sensing inversion method and a remote sensing inversion system for ground evapotranspiration, so as to avoid error accumulation of uncertainty of a remote sensing inversion model and uncertainty of a daily scale expansion method.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
As shown in fig. 1, the present invention provides a remote sensing inversion method of surface evapotranspiration, which comprises the following steps:
And 102, respectively calculating the soil temperature and the vegetation temperature relative to the soil moisture of 0 and the soil temperature and the vegetation temperature relative to the soil moisture of 1 according to the meteorological data and the remote sensing data of the earth surface.
Step 102 said method is based on said earth's surfaceThe meteorological data and the remote sensing data respectively calculate the soil temperature and vegetation temperature that relative soil moisture is 0, the soil temperature and vegetation temperature that relative soil moisture is 1, specifically include: using a formula based on the meteorological data and the remote sensing dataCalculating the soil temperature relative to the soil moisture of 0; using a formula based on the meteorological data and the remote sensing dataCalculating the vegetation temperature relative to the soil moisture of 0; using a formula based on the meteorological data and the remote sensing dataCalculating the soil temperature relative to the soil moisture of 1; using a formula based on the meteorological data and the remote sensing dataCalculating the vegetation temperature relative to the soil moisture of 1; wherein, Ts0And Tv0Respectively representing the soil temperature and vegetation temperature, T, relative to the soil moisture of 0s1And Tv1Respectively representing the soil temperature and the vegetation temperature relative to the soil moisture of 1; ρ is the air density; cpIs the specific heat at constant pressure; e.g. of the typesSaturated water vapor pressure; e.g. of the typeaThe actual water vapor pressure of the atmosphere; t isaThe temperature near the surface is obtained; gvwVegetation canopy conductivity sufficient for water supply; gavAnd gasThe aerodynamic conductivity of the vegetation and the upper layer of the soil are respectively; rn,sAnd Rn,vRespectively the soil component net radiation and the vegetation component net radiation; gsIs the soil heat flux.
103, calculating and distinguishing the soil dry edge temperature of the sufficient water state of the root area and the water-free state of the soil surface layer according to the soil temperature relative to the soil water content of 0, the vegetation temperature relative to the soil water content of 1 and the remote sensing data of the earth surface;
step 103 is as followsThe soil dry edge temperature for distinguishing the sufficient water state of the root zone and the moisture-free state of the soil surface layer is calculated according to the soil temperature relative to the soil moisture of 0, the vegetation temperature relative to the soil moisture of 1 and the remote sensing data of the earth surface, and the method specifically comprises the following steps: according to the surface remote sensing data, using a formulaCalculating the vegetation ratio Fv(ii) a Wherein EVI is Enhanced Vegetation Index (EVI)minAnd EVImaxRespectively the minimum enhancement type vegetation index corresponding to bare soil and the maximum enhancement type vegetation index corresponding to full vegetation coverage; according to the soil temperature with the relative soil moisture of 0 and the vegetation temperature with the relative soil moisture of 1 and the vegetation proportion, utilizing a formulaCalculating and distinguishing soil dry edge temperature T of sufficient water state and no water state on soil surface layer of root zone*(ii) a Wherein, Ts0And Tv1The soil temperature is 0 relative to the soil moisture and the vegetation temperature is 1 relative to the soil moisture.
And 104, calculating a soil humidity limiting factor, a vegetation humidity limiting factor and a vegetation temperature limiting factor according to the soil temperature and the vegetation temperature relative to the soil moisture of 0, the soil temperature and the vegetation temperature relative to the soil moisture of 1 and the soil dry edge temperature.
Step 104 the soil humidity limiting factor, vegetation humidity limiting factor and vegetation temperature limiting factor are calculated according to the soil temperature and vegetation temperature relative to the soil moisture of 0, the soil temperature and vegetation temperature relative to the soil moisture of 1 and the soil dry edge temperature, and specifically include: judging whether the surface temperature of the remote sensing image pixel is greater than the dry edge temperature of the soil or not to obtain a judgment result; if the judgment result shows no, the formula f is usedVMCalculating a vegetation moisture limiting factor f 1VM(ii) a Using formulasCalculating the soil humidity limiting factor fSM(ii) a If the judgment result shows yes, the formula f is utilizedSMCalculating the soil humidity limiting factor f as 0SM(ii) a Using formulasCalculating vegetation humidity limiting factor fVM(ii) a Wherein, Ts0And Ts1Respectively, a soil temperature of 0 relative to the soil moisture and a soil temperature of 1 relative to the soil moisture, TsoilWhich is indicative of the temperature of the soil components,TRfor remote-sensing the surface temperature of image elements, FvIs the vegetation proportion; t isv0And Tv1Respectively vegetation temperature of 0 relative to the soil moisture and vegetation temperature of 1 relative to the soil moisture, TvegIs the temperature of the components of the vegetation,using formulasCalculating a vegetation temperature limiting factor fT(ii) a Wherein, TmaxThe maximum atmospheric temperature in one day; t isoptThe temperature is the optimum temperature for vegetation growth.
And 105, calculating the surface evaporation capacity by using a Priestley-Taylor formula according to the soil humidity limiting factor, the vegetation humidity limiting factor and the temperature limiting factor.
105, calculating the surface evaporation capacity by using a Priestley-Taylor formula according to the soil humidity limiting factor, the vegetation humidity limiting factor and the temperature limiting factor, wherein the method specifically comprises the following steps: according to the vegetation temperature limiting factor and the vegetation humidity limiting factor, using a formulaCalculating the transpiration amount of the vegetation; according to the soil humidity limiting factor, using the formulaCalculating the evaporation capacity of the soil; by usingCalculating the interception evaporation capacity; according to the vegetation transpiration amount, the soil evaporation amount and the cut-off evaporation amount, utilizing a formula ET ═ ETv+ETs+ETiCalculating the evaporation capacity of the earth surface; wherein, ETv、ETsAnd ETiThe vegetation transpiration amount, the soil evaporation amount and the closure evaporation amount are respectively, and ET is the earth surface evaporation amount; f. ofwetIs the proportion of the wet leaves of the vegetation,RH is the relative humidity of the atmosphere; f. ofgThe ratio of vegetation to green leaves;fAPARphotosynthetically active radiation absorbed for green plants; f. ofIPARPhotosynthetically active radiation intercepted for green vegetation; f. ofTIs a vegetation temperature limiting factor; f. ofSMIs a soil humidity limiting factor; f. ofVMIs vegetation humidity limiting factor; alpha is the coefficient of a Priestley-Taylor formula, and alpha is 1.26; e.g. of the typesSaturated water vapor pressure; t isaThe temperature near the surface is obtained; rn,vFor canopy net radiation, Rn,sNet radiation for soil; gsIs the soil heat flux.
A remote sensing inversion system of surface evapotranspiration, the inversion system comprising:
the data acquisition module is used for acquiring meteorological data and remote sensing data of the earth surface; the meteorological data comprise air temperature, atmospheric pressure, air relative humidity, wind speed, uplink short-wave radiation, downlink short-wave radiation and the like; the remote sensing data comprises earth surface temperature, earth surface emissivity, reflectivity, leaf area index, enhanced vegetation index, photosynthetically active radiation and the like;
the soil temperature and vegetation temperature calculation module is used for calculating the soil temperature and vegetation temperature relative to the soil moisture of 0 and the soil temperature and vegetation temperature relative to the soil moisture of 1 respectively according to the meteorological data and the remote sensing data of the earth surface;
soil temperature and vegetation temperature calculation module specifically includes: a soil temperature calculation submodule with a relative soil moisture of 0 for utilizing a formula according to the meteorological data and the remote sensing dataCalculating the soil temperature relative to the soil moisture of 0; a vegetation temperature calculation submodule with a relative soil moisture of 0 for utilizing a formula according to the meteorological data and the remote sensing dataCalculating the vegetation temperature relative to the soil moisture of 0; a soil temperature calculation submodule with a relative soil moisture of 1 for utilizing a formula according to the meteorological data and the remote sensing dataCalculating the soil temperature relative to the soil moisture of 1; a vegetation temperature calculation submodule with relative soil moisture of 1 for utilizing a formula according to the meteorological data and the remote sensing dataCalculating the vegetation temperature relative to the soil moisture of 1; wherein, Ts0And Tv0Respectively representing the soil temperature and vegetation temperature, T, relative to the soil moisture of 0s1And Tv1Respectively representing the soil temperature and the vegetation temperature relative to the soil moisture of 1; ρ is the air density; cpIs the specific heat at constant pressure; e.g. of the typesSaturated water vapor pressure; e.g. of the typeaThe actual water vapor pressure of the atmosphere; t isaThe temperature near the surface is obtained; gvwVegetation canopy conductivity sufficient for water supply; gavAnd gasThe aerodynamic conductivity of the vegetation and the upper layer of the soil are respectively; rn,sAnd Rn,vRespectively the soil component net radiation and the vegetation component net radiation; gsIs the soil heat flux.
And the soil dry edge temperature calculation module is used for calculating and distinguishing the soil dry edge temperature of the sufficient water state of the root area and the soil surface layer water-free state according to the soil temperature relative to the soil water content of 0, the vegetation temperature relative to the soil water content of 1 and remote sensing data.
The soil dry edge temperature calculation module specifically comprises: a vegetation coverage calculation submodule for utilizing a formula according to the earth surface remote sensing dataCalculating the vegetation ratio Fv(ii) a Wherein EVI is Enhanced Vegetation Index (EVI)minAnd EVImaxRespectively the minimum enhancement type vegetation index corresponding to bare soil and the maximum enhancement type vegetation index corresponding to full vegetation coverage; a soil dry edge temperature operator module for utilizing a formula according to the soil temperature of 0 relative to the soil moisture and the vegetation temperature and vegetation proportion of 1 relative to the soil moistureCalculating and distinguishing soil dry edge temperature T of sufficient water state and no water state on soil surface layer of root zone*(ii) a Wherein, Ts0And Tv1The soil temperature is 0 relative to the soil moisture and the vegetation temperature is 1 relative to the soil moisture.
And the limiting factor calculation module is used for calculating a soil humidity limiting factor, a vegetation humidity limiting factor and a vegetation temperature limiting factor according to the soil temperature and the vegetation temperature relative to the soil moisture of 0 and the soil temperature and the vegetation temperature relative to the soil moisture of 1.
The limiting factor calculating module specifically comprises: the judgment submodule is used for judging whether the surface temperature of the remote sensing image pixel is greater than the dry edge temperature of the soil or not to obtain a judgment result; a first limiting factor calculating submodule for utilizing the formula f if the judgment result shows noVMCalculating vegetation humidity limiting factor f 1VM(ii) a Using formulasCalculating soil moisture limitsFactor fSM(ii) a A second limiting factor calculating submodule for using the formula f if the judgment result shows yesSMCalculating a vegetation humidity limiting factor fVM(ii) a Using formulasCalculating vegetation humidity limiting factor fVM(ii) a Wherein, Ts0And Ts1The soil temperature is 0 relative to the soil moisture and the soil temperature is 1 relative to the soil moisture; t issoilWhich is indicative of the temperature of the soil components,TRfor remote-sensing the surface temperature of image elements, FvIs the vegetation proportion; t isv0And Tv1Respectively vegetation temperature of 0 relative to the soil moisture and vegetation temperature of 1 relative to the soil moisture, TvegIs the temperature of the components of the vegetation,a third limiting factor calculation submodule for utilizing a formulaCalculating a vegetation temperature limiting factor fT(ii) a Wherein, TmaxThe maximum atmospheric temperature in one day; t isoptThe temperature is the optimum temperature for vegetation growth.
And the earth surface evaporation amount calculation module is used for calculating the earth surface evaporation amount by utilizing a Priestley-Taylor formula according to the soil humidity limiting factor, the vegetation humidity limiting factor and the temperature limiting factor.
The surface evaporation amount calculation module specifically comprises: a vegetation transpiration meter operator module for utilizing a formula according to the vegetation temperature limiting factor and the vegetation humidity limiting factorCalculating the transpiration amount of the vegetation; a soil evaporation meter operator module for utilizing a formula according to the soil humidity limiting factorCalculating the evaporation capacity of the soil; a cut-off evaporation amount calculation operator module for utilizingCalculating the interception evaporation capacity; the earth surface evapotranspiration meter operator module is used for utilizing a formula ET-ET according to the vegetation evapotranspiration, the soil evaporation and the interception evaporationv+ETs+ETiCalculating the evapotranspiration amount of the earth surface; wherein, ETv、ETsAnd ETiThe vegetation transpiration amount, the soil evaporation amount and the closure evaporation amount are respectively, and ET is the earth surface evaporation amount; f. ofwetIs the proportion of the wet leaves of the vegetation,RH is the relative humidity of the atmosphere; f. ofgThe ratio of vegetation to green leaves;fAPARphotosynthetically active radiation absorbed for green plants; f. ofIPARPhotosynthetically active radiation intercepted for green vegetation; f. ofTIs a vegetation temperature limiting factor; f. ofSMIs a soil humidity limiting factor; f. ofVMIs vegetation humidity limiting factor; alpha is the coefficient of a Priestley-Taylor formula, and alpha is 1.26; e.g. of the typesSaturated water vapor pressure; t isaThe temperature near the surface is obtained; rn,vFor canopy net radiation, Rn,sNet radiation for soil; gsIs the soil heat flux.
The invention provides a remote sensing inversion method of daily scale earth surface evapotranspiration, and belongs to the technical field of earth surface evapotranspiration remote sensing inversion. The method comprises the following steps: (A) collecting meteorological data and remote sensing data, and constructing an input data set of a daily scale earth surface evapotranspiration remote sensing inversion method; (B) calculating a soil humidity limiting factor and a vegetation humidity limiting factor by calculating the soil temperature and the vegetation temperature of 0 dry relative to the soil moisture, the soil temperature and the vegetation temperature and the soil dry edge temperature of 1 relative to the soil moisture; (C) calculating vegetation green leaf proportion, vegetation wet leaf proportion and vegetation temperature limiting factor according to the vegetation index, relative humidity, atmospheric temperature and vegetation growth optimum temperature; (D) and calculating vegetation transpiration, soil evaporation and cut-off evaporation according to a Priestley-Taylor formula, a vegetation green leaf proportion, a vegetation wet leaf proportion, a vegetation temperature limiting factor, a vegetation humidity limiting factor and a soil humidity limiting factor. The invention can reflect the water change of the soil surface layer and the vegetation root area on the daily scale by calculating the limiting factors of the soil temperature and the vegetation temperature relative to the soil moisture of 0 and the soil temperature and the vegetation temperature relative to the soil moisture of 1, directly provides the inversion method of daily scale surface evapotranspiration, and has important effect on improving the accuracy of remote sensing inversion of the surface evapotranspiration.
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other.
The principle and the implementation manner of the present invention are explained by applying specific examples, the above description of the embodiments is only used to help understanding the method of the present invention and the core idea thereof, the described embodiments are only a part of the embodiments of the present invention, not all embodiments, and all other embodiments obtained by a person of ordinary skill in the art based on the embodiments of the present invention without creative efforts belong to the protection scope of the present invention.
Claims (2)
1. A remote sensing inversion method of ground surface evapotranspiration is characterized by comprising the following steps:
acquiring meteorological data and remote sensing data of the earth surface;
respectively calculating the soil temperature and the vegetation temperature relative to the soil moisture of 0 and the soil temperature and the vegetation temperature relative to the soil moisture of 1 according to the meteorological data and the remote sensing data of the earth surface; the method specifically comprises the following steps: using a formula based on the meteorological data and the remote sensing dataCalculating the soil temperature relative to the soil moisture of 0; using a formula based on the meteorological data and the remote sensing dataCalculating the vegetation temperature relative to the soil moisture of 0; using a formula based on the meteorological data and the remote sensing dataCalculating the soil temperature relative to the soil moisture of 1; using a formula based on the meteorological data and the remote sensing dataCalculating the vegetation temperature relative to the soil moisture of 1; wherein, Ts0And Tv0Respectively representing the soil temperature and vegetation temperature, T, relative to the soil moisture of 0s1And Tv1Respectively representing the soil temperature and the vegetation temperature relative to the soil moisture of 1; ρ is the air density; cpIs the specific heat at constant pressure; e.g. of the typesSaturated water vapor pressure; e.g. of the typeaThe actual water vapor pressure of the atmosphere; t isaThe temperature near the surface is obtained; gvwVegetation canopy conductivity sufficient for water supply; gavAnd gasThe aerodynamic conductivity of the vegetation and the upper layer of the soil are respectively; rn,sAnd Rn,vRespectively the soil component net radiation and the vegetation component net radiation; gsIs the soil heat flux;
calculating and distinguishing the soil dry edge temperature of the sufficient water state of the root zone and the water-free state of the soil surface layer according to the soil temperature of 0 relative to the soil moisture, the vegetation temperature of 1 relative to the soil moisture and remote sensing data; the method specifically comprises the following steps: according to the surface remote sensing data, using a formulaCalculating the vegetation ratio Fv(ii) a Wherein EVI is Enhanced Vegetation Index (EVI)minAnd EVImaxMinimum enhanced vegetation index for the corresponding bare soil, respectivelyCorresponding to the maximum enhanced vegetation index of the fully planted and covered vegetation; according to the soil temperature with the relative soil moisture of 0 and the vegetation temperature with the relative soil moisture of 1 and the vegetation proportion, utilizing a formulaCalculating and distinguishing soil dry edge temperature T of sufficient water state and no water state on soil surface layer of root zone*(ii) a Wherein, Ts0And Tv1Respectively the soil temperature relative to the soil moisture of 0 and the vegetation temperature relative to the soil moisture of 1;
calculating a soil humidity limiting factor, a vegetation humidity limiting factor and a vegetation temperature limiting factor according to the soil temperature and the vegetation temperature relative to the soil moisture of 0, the soil temperature and the vegetation temperature relative to the soil moisture of 1 and the soil dry edge temperature; the method specifically comprises the following steps: judging whether the surface temperature of the remote sensing image pixel is greater than the dry edge temperature of the soil or not to obtain a judgment result; if the judgment result shows no, the formula f is usedVMCalculating a vegetation moisture limiting factor f 1VM(ii) a Using formulasCalculating the soil humidity limiting factor fSM(ii) a If the judgment result shows yes, the formula f is utilizedSMCalculating the soil humidity limiting factor f as 0SM(ii) a Using formulasCalculating vegetation humidity limiting factor fVM(ii) a Wherein, Ts0And Ts1Respectively, a soil temperature of 0 relative to the soil moisture and a soil temperature of 1 relative to the soil moisture, TsoilWhich is indicative of the temperature of the soil components,TRfor remote-sensing the surface temperature of image elements, FvIs the vegetation proportion; t isv0And Tv1Respectively vegetation temperature of 0 relative to soil moisture and vegetation temperature of 1 relative to soil moistureVegetation temperature, TvegIs the temperature of the components of the vegetation,using formulasCalculating a vegetation temperature limiting factor fT(ii) a Wherein, TmaxThe maximum atmospheric temperature in one day; t isoptThe optimum temperature for vegetation growth;
calculating the surface evaporation capacity by using a Priestley-Taylor formula according to the soil humidity limiting factor, the vegetation humidity limiting factor and the vegetation temperature limiting factor; the method specifically comprises the following steps: according to the vegetation temperature limiting factor and the vegetation humidity limiting factor, using a formulaCalculating the transpiration amount of the vegetation; according to the soil humidity limiting factor, using the formulaCalculating the evaporation capacity of the soil; by usingCalculating the interception evaporation capacity; according to the vegetation transpiration amount, the soil evaporation amount and the cut-off evaporation amount, utilizing a formula ET ═ ETv+ETs+ETiCalculating the evaporation capacity of the earth surface; wherein, ETv、ETsAnd ETiThe vegetation transpiration amount, the soil evaporation amount and the closure evaporation amount are respectively, and ET is the earth surface evaporation amount; f. ofwetIs the proportion of the wet leaves of the vegetation,RH is the relative humidity of the atmosphere; f. ofgThe ratio of vegetation to green leaves;fAPARis a green vegetationAbsorbed photosynthetically active radiation; f. ofIPARPhotosynthetically active radiation intercepted for green vegetation; f. ofTIs a vegetation temperature limiting factor; f. ofSMIs a soil humidity limiting factor; f. ofVMIs vegetation humidity limiting factor; alpha is the coefficient of a Priestley-Taylor formula, and alpha is 1.26; e.g. of the typesSaturated water vapor pressure; t isaThe temperature near the surface is obtained; rn,vFor canopy net radiation, Rn,sNet radiation for soil; gsIs the soil heat flux.
2. A remote sensing inversion system of surface evapotranspiration, the inversion system comprising:
the data acquisition module is used for acquiring meteorological data and remote sensing data of the earth surface;
the soil temperature and vegetation temperature calculation module is used for calculating the soil temperature and vegetation temperature relative to the soil moisture of 0 and the soil temperature and vegetation temperature relative to the soil moisture of 1 respectively according to the meteorological data and the remote sensing data of the earth surface;
soil temperature and vegetation temperature calculation module specifically includes: a soil temperature calculation submodule with a relative soil moisture of 0 for utilizing a formula according to the meteorological data and the remote sensing dataCalculating the soil temperature relative to the soil moisture of 0; a vegetation temperature calculation submodule with a relative soil moisture of 0 for utilizing a formula according to the meteorological data and the remote sensing dataCalculating the vegetation temperature relative to the soil moisture of 0; a soil temperature calculation submodule with a relative soil moisture of 1 for utilizing a formula according to the meteorological data and the remote sensing dataCalculating the soil temperature relative to the soil moisture of 1; vegetation with relative soil moisture of 1A temperature calculation submodule for utilizing a formula according to the meteorological data and the remote sensing dataCalculating the vegetation temperature relative to the soil moisture of 1; wherein, Ts0And Tv0Respectively representing the soil temperature and vegetation temperature, T, relative to the soil moisture of 0s1And Tv1Respectively representing the soil temperature and the vegetation temperature relative to the soil moisture of 1; ρ is the air density; cpIs the specific heat at constant pressure; e.g. of the typesSaturated water vapor pressure; e.g. of the typeaThe actual water vapor pressure of the atmosphere; t isaThe temperature near the surface is obtained; gvwVegetation canopy conductivity sufficient for water supply; gavAnd gasThe aerodynamic conductivity of the vegetation and the upper layer of the soil are respectively; rn,sAnd Rn,vRespectively the soil component net radiation and the vegetation component net radiation; gsIs the soil heat flux;
the soil dry edge temperature calculation module is used for calculating and distinguishing the soil dry edge temperature of the sufficient water state of the root area and the no water state of the soil surface layer according to the soil temperature relative to the soil water content of 0, the vegetation temperature relative to the soil water content of 1 and remote sensing data;
the soil dry edge temperature calculation module specifically comprises: a vegetation coverage calculation submodule for utilizing a formula according to the earth surface remote sensing dataCalculating the vegetation ratio Fv(ii) a Wherein EVI is Enhanced Vegetation Index (EVI)minAnd EVImaxRespectively the minimum enhancement type vegetation index corresponding to bare soil and the maximum enhancement type vegetation index corresponding to full vegetation coverage; a soil dry edge temperature operator module for utilizing a formula according to the soil temperature of 0 relative to the soil moisture and the vegetation temperature and vegetation proportion of 1 relative to the soil moistureCalculating and distinguishing root zone water sufficient stateAnd the dry edge temperature T of the soil in a state of no moisture on the surface layer of the soil*(ii) a Wherein, Ts0And Tv1Respectively the soil temperature relative to the soil moisture of 0 and the vegetation temperature relative to the soil moisture of 1;
the limiting factor calculation module is used for calculating a soil humidity limiting factor, a vegetation humidity limiting factor and a vegetation temperature limiting factor according to the soil temperature and the vegetation temperature relative to the soil moisture of 0, the soil temperature and the vegetation temperature relative to the soil moisture of 1 and the soil dry edge temperature;
the limiting factor calculating module specifically comprises: the judgment submodule is used for judging whether the surface temperature of the remote sensing image pixel is greater than the dry edge temperature of the soil or not to obtain a judgment result; a first limiting factor calculating submodule for utilizing the formula f if the judgment result shows noVMCalculating vegetation humidity limiting factor f 1VM(ii) a Using formulasCalculating the soil humidity limiting factor fSM(ii) a A second limiting factor calculating submodule for using the formula f if the judgment result shows yesSMCalculating a vegetation humidity limiting factor fVM(ii) a Using formulasCalculating vegetation humidity limiting factor fVM(ii) a Wherein, Ts0And Ts1The soil temperature is 0 relative to the soil moisture and the soil temperature is 1 relative to the soil moisture; t issoilWhich is indicative of the temperature of the soil components,TRfor remote-sensing the surface temperature of image elements, FvIs the vegetation proportion; t isv0And Tv1Respectively vegetation temperature of 0 relative to the soil moisture and vegetation temperature of 1 relative to the soil moisture, TvegIs the temperature of the components of the vegetation,a third limiting factor calculation submodule for utilizing a formulaCalculating a vegetation temperature limiting factor fT(ii) a Wherein, TmaxThe maximum atmospheric temperature in one day; t isoptThe optimum temperature for vegetation growth;
the earth surface evaporation amount calculation module is used for calculating earth surface evaporation amount by using a Priestley-Taylor formula according to the soil humidity limiting factor, the vegetation humidity limiting factor and the temperature limiting factor;
the surface evaporation amount calculation module specifically comprises: a vegetation transpiration meter operator module for utilizing a formula according to the vegetation temperature limiting factor and the vegetation humidity limiting factorCalculating the transpiration amount of the vegetation; a soil evaporation meter operator module for utilizing a formula according to the soil humidity limiting factorCalculating the evaporation capacity of the soil; a cut-off evaporation amount calculation operator module for utilizingCalculating the interception evaporation capacity; the earth surface evapotranspiration meter operator module is used for utilizing a formula ET-ET according to the vegetation evapotranspiration, the soil evaporation and the interception evaporationv+ETs+ETiCalculating the evapotranspiration amount of the earth surface; wherein, ETv、ETsAnd ETiThe vegetation transpiration amount, the soil evaporation amount and the closure evaporation amount are respectively, and ET is the earth surface evaporation amount; f. ofwetIs the proportion of the wet leaves of the vegetation,RH is the relative humidity of the atmosphere; f. ofgThe ratio of vegetation to green leaves;fAPARphotosynthetically active radiation absorbed for green plants; f. ofIPARPhotosynthetically active radiation intercepted for green vegetation; f. ofTIs a vegetation temperature limiting factor; f. ofSMIs a soil humidity limiting factor; f. ofVMIs vegetation humidity limiting factor; alpha is the coefficient of a Priestley-Taylor formula, and alpha is 1.26; e.g. of the typesSaturated water vapor pressure; t isaThe temperature near the surface is obtained; rn,vFor canopy net radiation, Rn,sNet radiation for soil; gsIs the soil heat flux.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010895830.0A CN112014323B (en) | 2020-08-31 | 2020-08-31 | Remote sensing inversion method and system for surface evapotranspiration |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010895830.0A CN112014323B (en) | 2020-08-31 | 2020-08-31 | Remote sensing inversion method and system for surface evapotranspiration |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112014323A CN112014323A (en) | 2020-12-01 |
CN112014323B true CN112014323B (en) | 2021-07-20 |
Family
ID=73503957
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010895830.0A Active CN112014323B (en) | 2020-08-31 | 2020-08-31 | Remote sensing inversion method and system for surface evapotranspiration |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112014323B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112733445A (en) * | 2021-01-05 | 2021-04-30 | 四川大学 | Large-area-scale soil moisture inversion method based on evapotranspiration vegetation index spatial characteristics |
CN113075012B (en) * | 2021-03-29 | 2021-11-30 | 南京晓庄学院 | Moisture evaporation device for remotely sensing and monitoring near-surface soil humidity |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106771073A (en) * | 2016-12-28 | 2017-05-31 | 中国科学院地理科学与资源研究所 | A kind of method that soil and vegetation evapotranspiration are estimated based on end member information model |
CN109142674A (en) * | 2018-08-02 | 2019-01-04 | 中国科学院地理科学与资源研究所 | The remote sensing inversion method that root zone soil moisture opposite with surface layer is estimated simultaneously |
CN110059362A (en) * | 2019-03-22 | 2019-07-26 | 兰州大学 | Consider the construction method of the regional scale double source evapotranspiration model of vegetation dynamic changes |
WO2019145895A1 (en) * | 2018-01-24 | 2019-08-01 | The State Of Israel, Ministry Of Agriculture & Rural Development Agricultural Research Organization | Method and system for estimating crop coefficient and evapotranspiration of crops based on remote sensing |
CN110321784A (en) * | 2019-05-08 | 2019-10-11 | 中国科学院地理科学与资源研究所 | Method, apparatus, electronic equipment and the computer media of soil moisture estimation |
CN110334381A (en) * | 2019-05-08 | 2019-10-15 | 中国科学院地理科学与资源研究所 | The method, apparatus and computer media of evapotranspiration estimation |
-
2020
- 2020-08-31 CN CN202010895830.0A patent/CN112014323B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106771073A (en) * | 2016-12-28 | 2017-05-31 | 中国科学院地理科学与资源研究所 | A kind of method that soil and vegetation evapotranspiration are estimated based on end member information model |
WO2019145895A1 (en) * | 2018-01-24 | 2019-08-01 | The State Of Israel, Ministry Of Agriculture & Rural Development Agricultural Research Organization | Method and system for estimating crop coefficient and evapotranspiration of crops based on remote sensing |
CN109142674A (en) * | 2018-08-02 | 2019-01-04 | 中国科学院地理科学与资源研究所 | The remote sensing inversion method that root zone soil moisture opposite with surface layer is estimated simultaneously |
CN110059362A (en) * | 2019-03-22 | 2019-07-26 | 兰州大学 | Consider the construction method of the regional scale double source evapotranspiration model of vegetation dynamic changes |
CN110321784A (en) * | 2019-05-08 | 2019-10-11 | 中国科学院地理科学与资源研究所 | Method, apparatus, electronic equipment and the computer media of soil moisture estimation |
CN110334381A (en) * | 2019-05-08 | 2019-10-15 | 中国科学院地理科学与资源研究所 | The method, apparatus and computer media of evapotranspiration estimation |
Non-Patent Citations (1)
Title |
---|
"遥感反演蒸散发的日尺度扩展方法研究进展";王桐等;《遥感学报》;20191231;第813-830页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112014323A (en) | 2020-12-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ca et al. | Reductions in air conditioning energy caused by a nearby park | |
Yang et al. | Study on the local climatic effects of large photovoltaic solar farms in desert areas | |
Molion | Amazonia rainfall and its variability | |
Ramos et al. | Evaluation of satellite evapotranspiration estimates using ground-meteorological data available for the Flumen District into the Ebro Valley of NE Spain | |
CN112014323B (en) | Remote sensing inversion method and system for surface evapotranspiration | |
Nicholson et al. | Micrometeorological conditions and surface mass and energy fluxes on Lewis Glacier, Mt Kenya, in relation to other tropical glaciers | |
Rorison et al. | Local climate, topography and plant growth in Lathkill Dale NNR. I. A twelve‐year summary of solar radiation and temperature | |
CN111814317B (en) | Surface energy balance component estimation method and system based on remote sensing | |
Ali Khan Niazi et al. | Comparative analysis of photovoltaic configurations for agrivoltaic systems in Europe | |
Linacre | Net radiation to various surfaces | |
CN115796659A (en) | Method for estimating urban water storage capacity based on regional scale evaporation degradation | |
Stanhill et al. | The climate of the cotton crop: Physical characteristics and microclimate relationships | |
Iziomon et al. | Empirical models for estimating net radiative flux: A case study for three mid-latitude sites with orographic variability | |
Graser et al. | Within-canopy temperature patterns of sorghum at two row spacings | |
Dickinson | Evapotranspiration in global climate models | |
LeDrew et al. | A comparison of the radiation and energy balance during the growing season for an arctic and alpine tundra | |
Kim et al. | Combined effects of wind and solar irradiance on the spatial variation of midday air temperature over a mountainous terrain | |
Hojdová et al. | Microclimate of a peat bog and of the forest in different states of damage in the Šumava National Park | |
Nullet | Water balance of pacific atolls 1 | |
Udo et al. | New empirical relationships for determining global PAR from measurements of global solar radiation, infrared radiation or sunshine duration | |
Decker | Precision of Estimates of Evapotranspiration in Missouri Climate 1 | |
Lockwood | The influence of temperature variations on interception loss and water storage in vegetation canopies | |
Chaowiwat et al. | Prediction of future agriculture water demand in Thailand using multi bias corrected climate models | |
Gajere et al. | Estimation of Global Solar Radiation of Plateau state, Nigeria | |
Boyce | An historical comparison between the Met Office Surface Exchange Scheme-Probability Distributed Model (MOSES-PDM) and the Met Office Rainfall and Evaporation Calculation System (MORECS) |
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 |