Crop evapotranspiration prediction method, device, equipment, storage medium and program product

Abstract

The invention provides a crop evapotranspiration prediction method and device, equipment, a storage medium and a program product. The method comprises the steps of determining reference evapotranspiration data based on evapotranspiration input data; determining the relationship between the soil water potential and the soil moisture stress coefficient based on the evapotranspiration measured data and the reference evapotranspiration data; and based on the reference evapotranspiration data and the relationship, predicting the crop evapotranspiration of the to-be-predicted region to obtain a crop evapotranspiration prediction result. According to the method, reference evapotranspiration data is determined through evapotranspiration input data obtained from a to-be-predicted area, evapotranspiration measured data obtained from the to-be-predicted area are used for jointly determining the relation between the soil water potential and the soil moisture stress coefficient, and finally the crop evapotranspiration amount of the to-be-predicted area is predicted based on the reference evapotranspiration data and the determined relation. And obtaining a crop evapotranspiration prediction result. According to the method, the abiotic stress condition of the soil is effectively judged through the soil water potential, and the accuracy of crop evapotranspiration prediction is improved.

Description

Technical Field

[0001] The present invention relates to the technical field of crop monitoring, and in particular to a crop evapotranspiration prediction method, device, equipment, storage medium and program product. Background Art

[0002] Evapotranspiration is the total flux of water vapor transported from vegetation and the ground to the atmosphere. Evapotranspiration is an important part of energy balance and water cycle. Evapotranspiration data not only affects the growth and development of plants, but also affects atmospheric circulation and thus the climate. Therefore, the prediction of crop evapotranspiration helps to understand the various aspects of the ecosystem.

[0003] The existing crop reference evapotranspiration prediction methods are usually based on physical process models and statistical models. The estimation of crop actual evapotranspiration is based on the reference evapotranspiration, and the crop coefficient and soil moisture stress coefficient must also be considered. Existing methods usually try to establish the relationship between soil moisture and soil moisture stress coefficient, and then further estimate the actual crop evapotranspiration. However, soil moisture is not sensitive enough to soil changes. Therefore, when facing abiotic stresses of soil moisture (such as drought, freeze-thaw and salinity, etc.), the actual crop evapotranspiration cannot be accurately estimated. Summary of the invention

[0004] The present invention provides a crop evapotranspiration prediction method, device, equipment, storage medium and program product, which are used to solve the defect of inaccurate prediction of existing crop actual evapotranspiration prediction schemes and achieve more accurate prediction of crop actual evapotranspiration.

[0005] The invention provides a crop evapotranspiration prediction method, comprising the following steps.

[0006] Determining reference evapotranspiration data based on evapotranspiration input data; the evapotranspiration input data is obtained from the area to be predicted; Determining the relationship between soil water potential and soil water stress coefficient based on the evapotranspiration measured data and the reference evapotranspiration data; the evapotranspiration measured data is obtained from the area to be predicted; Based on the reference evapotranspiration data and the relationship, the crop evapotranspiration of the area to be predicted is predicted to obtain a crop evapotranspiration prediction result.

[0007] According to a crop evapotranspiration prediction method provided by the present invention, determining reference evapotranspiration data based on evapotranspiration input data comprises: Obtain evapotranspiration input data; the evapotranspiration input data includes crop surface net radiation , soil heat flux , average temperature , the slope of the relationship between saturated water vapor pressure and temperature , wind speed , saturated vapor pressure and the actual vapor pressure ; Based on the evapotranspiration input data and formula Determine the reference evapotranspiration data; where, is a reference evapotranspiration rate included in the reference evapotranspiration data; is a constant; a, b, c and d are parameters.

[0008] According to a crop evapotranspiration prediction method provided by the present invention, the relationship between soil water potential and soil water stress coefficient determined based on the evapotranspiration measured data and the reference evapotranspiration data includes: Determine the crop coefficient; Based on the crop coefficient and formula Determine the soil water stress coefficient, where is the measured evapotranspiration rate included in the evapotranspiration measured data; is the soil moisture stress coefficient; is the crop coefficient; The relationship between soil water potential and the soil water stress coefficient is determined by symbolic regression based on genetic programming.

[0009] According to a crop evapotranspiration prediction method provided by the present invention, the relationship between soil water potential and the soil water stress coefficient determined by the symbolic regression based on genetic programming includes: Determine the initial population and fitness function; Determining the fitness of each individual in the initialized population based on the fitness function; Updating the initialized population based on the fitness of each individual in the initialized population to obtain an updated population; In the case where the renewal population meets the termination condition, the relationship between the soil water potential and the soil water stress coefficient is determined based on the fitness of each individual of the renewal population.

[0010] According to a crop evapotranspiration prediction method provided by the present invention, the updating of the initialization population based on the fitness of each individual in the initialization population to obtain the updated population includes: Based on the fitness of each individual in the initialized population, an update operation is performed on the initialized population to obtain an updated population; the update operation includes a selection operation, a crossover operation and a mutation operation; wherein, a selection operation is performed on the initialized population to obtain parent individuals; Performing a crossover operation on the parent individuals to obtain new offspring individuals; A mutation operation is performed on the individual part of the initialized population to obtain new individuals.

[0011] According to a crop evapotranspiration prediction method provided by the present invention, the crop evapotranspiration of the to-be-predicted area is predicted based on the reference evapotranspiration data and the relationship, and obtaining the crop evapotranspiration prediction result includes: Based on the formula The crop evapotranspiration of the area to be predicted is predicted to obtain a crop evapotranspiration prediction result; wherein, is the actual evapotranspiration; is the soil water stress coefficient determined based on the relationship; The actual evapotranspiration is evaluated based on the regression prediction model evaluation index to obtain a crop evapotranspiration prediction result including the evaluation result.

[0012] The present invention also provides a crop evapotranspiration prediction device, comprising the following modules: A reference evapotranspiration data determination module, used to determine reference evapotranspiration data based on evapotranspiration input data; the evapotranspiration input data is obtained from a region to be predicted; A relationship determination module, used to determine the relationship between soil water potential and soil water stress coefficient based on evapotranspiration measured data and the reference evapotranspiration data; the evapotranspiration measured data is obtained from the area to be predicted; The crop evapotranspiration prediction module is used to predict the crop evapotranspiration of the area to be predicted based on the reference evapotranspiration data and the relationship to obtain a crop evapotranspiration prediction result.

[0013] The present invention also provides an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, any of the crop evapotranspiration prediction methods described above is implemented.

[0014] The present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, and when the computer program is executed by a processor, the crop evapotranspiration prediction method as described in any one of the above is implemented.

[0015] The present invention also provides a computer program product, comprising a computer program, wherein when the computer program is executed by a processor, the crop evapotranspiration prediction method as described above is implemented.

[0016] The crop evapotranspiration prediction method, device, equipment, storage medium and program product provided by the present invention determine the reference evapotranspiration data by using the evapotranspiration input data obtained from the area to be predicted, and the evapotranspiration measured data obtained from the area to be predicted, and jointly determine the relationship between the soil water potential and the soil water stress coefficient, and finally predict the crop evapotranspiration of the area to be predicted based on the reference evapotranspiration data and the determined relationship, and obtain the crop evapotranspiration prediction result. The present application effectively judges the abiotic stress of the soil through the soil water potential, thereby improving the accuracy of the crop evapotranspiration prediction. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] In order to more clearly illustrate the technical solutions in the present invention or the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0018] Figure 1 This is one of the flow charts of the crop evapotranspiration prediction method provided by the present invention.

[0019] Figure 2 This is the second flow chart of the crop evapotranspiration prediction method provided by the present invention.

[0020] Figure 3 It is a structural schematic diagram of the crop evapotranspiration prediction device provided by the present invention.

[0021] Figure 4 It is a structural schematic diagram of the electronic device provided by the present invention. DETAILED DESCRIPTION

[0022] In order to make the purpose, technical solution and advantages of the present invention clearer, the technical solution of the present invention will be clearly and completely described below in conjunction with the drawings of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0023] Combine the following Figure 1-Figure 4 The present invention describes the crop evapotranspiration prediction method, device, equipment, storage medium and program product.

[0024] Figure 1 FIG. 1 is one of the flow charts of the crop evapotranspiration prediction method provided by the present invention, such as Figure 1 As shown, the method includes the following: Step 100: Determine reference evapotranspiration data based on evapotranspiration input data; the evapotranspiration input data is obtained from a region to be predicted; The crop evapotranspiration prediction method provided in the embodiment of the present application may further include: Obtain evapotranspiration input data; the evapotranspiration input data includes crop surface net radiation , soil heat flux , average temperature , the slope of the relationship between saturated water vapor pressure and temperature , wind speed , saturated vapor pressure and the actual vapor pressure ; Based on the evapotranspiration input data and formula Determine the reference evapotranspiration data; where, is a reference evapotranspiration rate included in the reference evapotranspiration data; is a constant; a, b, c and d are parameters.

[0025] Specifically, the input data of the relevant evapotranspiration prediction model of the area to be predicted and the measured evapotranspiration value of the area to be predicted are obtained and preprocessed. The evapotranspiration input data to be obtained mainly include geographic information, meteorological data, crop factors, soil factors, etc.

[0026] The evapotranspiration input data further include the altitude of the study area, the longitude and latitude of the study area, solar radiation, wind speed at a height of 2 meters, daily maximum temperature, daily minimum temperature, daily maximum humidity, daily minimum humidity, leaf area index, and soil water potential, etc. In addition, the corresponding evapotranspiration measured values ​​need to be obtained.

[0027] The reference evapotranspiration is calculated by the Penman-Monteith formula (evapotranspiration calculation formula), and the specific calculation is shown in formula 1, where: represents the reference evaporation rate, which can be in millimeters per day; is the net radiation on the crop surface, in megajoules per square meter per day; is the soil heat flux in megajoules per square meter per day; It can be the daily average temperature at a height of 2 meters above the ground, in degrees Celsius; It can be the wind speed at a height of 2 meters above the ground, in meters per second; is the saturated water vapor pressure, the unit can be kilopascal; is the actual water vapor pressure, the unit can be kilopascal; It is the slope of the relationship between saturated water vapor pressure and air temperature, and the unit can be kilopascals per degree Celsius.

[0028] ; (1) Among them, a, b, c and d are parameters. Specifically, a can be 0.408, which is an empirical coefficient determined through experiments and data analysis; b can be 900, 900 represents the relationship between saturated water vapor pressure and temperature at standard atmospheric pressure (1013.25 hPa) and 25°C; c can be 273, T+273 is the conversion formula for absolute temperature; d can be 0.34, which is an empirical coefficient determined through experiments and data analysis.

[0029] Step 200: determining the relationship between soil water potential and soil water stress coefficient based on the evapotranspiration measured data and the reference evapotranspiration data; the evapotranspiration measured data is obtained from the area to be predicted; The crop evapotranspiration prediction method provided in the embodiment of the present application may further include: Step 210, determining the crop coefficient; Step 220: Based on the crop coefficient and formula Determine the soil water stress coefficient, where is the measured evapotranspiration rate included in the evapotranspiration measured data; is the soil moisture stress coefficient; is the crop coefficient; Step 230: Determine the relationship between soil water potential and the soil water stress coefficient based on symbolic regression of genetic programming.

[0030] Specifically, the relationship between soil water potential and soil water stress coefficient was established by using symbolic regression based on genetic programming, and the soil water stress coefficient was estimated by soil water potential.

[0031] It specifically includes the following steps: Step 1, calculate the soil moisture stress coefficient according to the evapotranspiration formula; Step 2: Input data set division and preprocessing; Step 3: Implement symbolic regression based on genetic programming to fit the relationship between soil water potential and soil water stress coefficient.

[0032] Among them, the specific contents of step 1 include: the corresponding crop coefficient , the measured evapotranspiration and the reference evapotranspiration calculated above Substitute into the evapotranspiration formula The soil moisture stress coefficient is calculated ; Step 2 specifically includes: defining the data set, given the data of soil water potential and the corresponding soil water stress coefficient, and preprocessing the original data. A part of the data (e.g. 80%) is used as the training set, and the other part (e.g. 20%) is used as the validation set. The training goal is to fit the relationship between soil water potential and soil water stress coefficient and output the corresponding relationship function.

[0033] Step 300: Based on the reference evapotranspiration data and the relationship, the crop evapotranspiration of the area to be predicted is predicted to obtain a crop evapotranspiration prediction result.

[0034] The crop evapotranspiration prediction method provided in the embodiment of the present application may further include: Step 310: Based on the formula The crop evapotranspiration of the area to be predicted is predicted to obtain a crop evapotranspiration prediction result; wherein, is the actual evapotranspiration; is the soil water stress coefficient determined based on the relationship; Step 320: Evaluate the actual evapotranspiration based on the regression prediction model evaluation index to obtain a crop evapotranspiration prediction result including the evaluation result.

[0035] Specifically, according to the formula obtained from the above actual evapotranspiration calculation , predict the actual evapotranspiration of crops in the region and analyze the model prediction performance.

[0036] in is the reference evaporation rate obtained above, is the soil moisture stress coefficient obtained by the above symbolic regression (determined based on the relationship), is the crop coefficient. The mean square error (MSE), root mean square error (RMSE), and mean absolute error (MAE) can be used to evaluate the crop evapotranspiration prediction results. As shown in formula 2-4, is the true value, is the predicted value.

[0037] ; (2) ; (3) ; (4) This embodiment determines reference evapotranspiration data by using evapotranspiration input data obtained from the area to be predicted, and evapotranspiration measured data obtained from the area to be predicted, jointly determines the relationship between soil water potential and soil water stress coefficient, and finally predicts the crop evapotranspiration of the area to be predicted based on the reference evapotranspiration data and the determined relationship, and obtains the crop evapotranspiration prediction result. This application effectively judges the abiotic stress of the soil through soil water potential, thereby improving the accuracy of crop evapotranspiration prediction.

[0038] Figure 2 FIG. 2 is a flow chart of the crop evapotranspiration prediction method provided by the present invention. Figure 2 As shown, the method may also include: Step 231, determine the initialization population and fitness function; Step 232: determining the fitness of each individual of the initialized population based on the fitness function; Step 233: updating the initialized population based on the fitness of each individual in the initialized population to obtain an updated population; Step 234: When the renewal population meets the termination condition, determine the relationship between the soil water potential and the soil water stress coefficient based on the fitness of each individual in the renewal population.

[0039] The crop evapotranspiration prediction method provided in the embodiment of the present application may further include: Step 233-1: Based on the fitness of each individual in the initialized population, an update operation is performed on the initialized population to obtain an updated population; the update operation includes a selection operation, a crossover operation, and a mutation operation; wherein, a selection operation is performed on the initialized population to obtain parent individuals; Performing a crossover operation on the parent individuals to obtain new offspring individuals; A mutation operation is performed on the individual part of the initialized population to obtain new individuals.

[0040] Specifically, the content of the above step 3 mainly includes: Step 3.1, generate the initial population, and then set the parameters of the initial population, including population size, maximum genetic generation, crossover probability and mutation probability, etc.

[0041] Step 3.2, define the fitness function, each individual in the population is a mathematical expression; the fitness function can be selected from the mean square error formula , as shown in Formula 2.

[0042] Step 3.3: Perform selection, crossover and mutation operations on individuals in the population to update the population.

[0043] Step 3.4, continuously perform the above selection operation, crossover operation and mutation operation on the population to update the population. After the update, determine whether the updated population meets the preset termination condition. The preset termination condition can be set to 500 training times and the updated fitting error is less than 0.01. When the preset termination condition is met after the update, the individual with the best fitness is selected as the best mathematical expression obtained by training.

[0044] This embodiment establishes the relationship between soil water potential and soil water stress coefficient through symbolic regression based on genetic programming, and accurately estimates the soil water stress coefficient through soil water potential.

[0045] The crop evapotranspiration prediction device provided by the present invention is described below. The crop evapotranspiration prediction device described below and the crop evapotranspiration prediction method described above can be referred to each other.

[0046] Please refer to Figure 3 The present invention also provides a crop evapotranspiration prediction device, comprising: A reference evapotranspiration data determination module 301 is used to determine reference evapotranspiration data based on evapotranspiration input data; the evapotranspiration input data is obtained from a region to be predicted; A relationship determination module 302 is used to determine the relationship between soil water potential and soil water stress coefficient based on the evapotranspiration measured data and the reference evapotranspiration data; the evapotranspiration measured data is obtained from the area to be predicted; The crop evapotranspiration prediction module 303 is used to predict the crop evapotranspiration of the to-be-predicted area based on the reference evapotranspiration data and the relationship to obtain a crop evapotranspiration prediction result.

[0047] Optionally, the reference evapotranspiration data determination module includes: The evapotranspiration input data acquisition unit is used to acquire the evapotranspiration input data; the evapotranspiration input data includes the crop surface net radiation , soil heat flux , average temperature , the slope of the relationship between saturated water vapor pressure and temperature , wind speed , saturated vapor pressure and the actual vapor pressure ; Reference evapotranspiration data determination unit for inputting data and formula based on the evapotranspiration Determine the reference evapotranspiration data; where, is a reference evapotranspiration rate included in the reference evapotranspiration data; is a constant; a, b, c and d are parameters.

[0048] Optionally, the relationship determination module includes: A crop coefficient determination unit, used for determining a crop coefficient; A soil moisture stress coefficient determination unit is used to determine the soil moisture stress coefficient based on the crop coefficient and formula Determine the soil water stress coefficient, where is the measured evapotranspiration rate included in the evapotranspiration measured data; is the soil moisture stress coefficient; is the crop coefficient; The relationship determination unit between soil water potential and the soil water stress coefficient is used to determine the relationship between soil water potential and the soil water stress coefficient based on symbolic regression of genetic programming.

[0049] Optionally, the unit for determining the relationship between the soil water potential and the soil water stress coefficient includes: An initialization population and fitness function determination unit, used to determine the initialization population and fitness function; An individual fitness determination unit, used to determine the fitness of each individual of the initialized population based on the fitness function; An initialization population updating unit, used for updating the initialization population based on the fitness of each individual in the initialization population to obtain an updated population; A relationship determination unit is used to determine the relationship between the soil water potential and the soil moisture stress coefficient based on the fitness of each individual of the renewal population when the renewal population meets the termination condition.

[0050] Optionally, the initialization population update unit includes: An initialization population update operation unit, used to perform an update operation on the initialization population based on the fitness of each individual in the initialization population to obtain an updated population; the update operation includes a selection operation, a crossover operation and a mutation operation; wherein, a selection operation is performed on the initialized population to obtain parent individuals; Performing a crossover operation on the parent individuals to obtain new offspring individuals; A mutation operation is performed on the individual part of the initialized population to obtain new individuals.

[0051] Optionally, the crop evapotranspiration prediction module includes: The crop evapotranspiration prediction result determination unit is used to determine the crop evapotranspiration prediction result based on the formula The crop evapotranspiration of the area to be predicted is predicted to obtain a crop evapotranspiration prediction result; wherein, is the actual evapotranspiration; is the soil water stress coefficient determined based on the relationship; The actual evapotranspiration evaluation unit is used to evaluate the actual evapotranspiration based on the regression prediction model evaluation index to obtain a crop evapotranspiration prediction result including the evaluation result.

[0052] Figure 4 An example of a physical structure diagram of an electronic device is shown in FIG. Figure 4As shown, the electronic device may include: a processor 410, a communication interface 420, a memory 430 and a communication bus 440, wherein the processor 410, the communication interface 420 and the memory 430 communicate with each other through the communication bus 440. The processor 410 may call the logic instructions in the memory 430 to execute the crop evapotranspiration prediction method, the method comprising: determining reference evapotranspiration data based on evapotranspiration input data; the evapotranspiration input data is obtained from the area to be predicted; determining the relationship between soil water potential and soil water stress coefficient based on evapotranspiration measured data and the reference evapotranspiration data; the evapotranspiration measured data is obtained from the area to be predicted; based on the reference evapotranspiration data and the relationship, predicting the crop evapotranspiration of the area to be predicted, and obtaining the crop evapotranspiration prediction result.

[0053] In addition, the logic instructions in the above-mentioned memory 430 can be implemented in the form of a software functional unit and can be stored in a computer-readable storage medium when it is sold or used as an independent product. Based on this understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art or the part of the technical solution, can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including a number of instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk, etc. Various media that can store program codes.

[0054] On the other hand, the present invention also provides a computer program product, which includes a computer program. The computer program can be stored on a non-transitory computer-readable storage medium. When the computer program is executed by a processor, the computer can execute the crop evapotranspiration prediction method provided by the above methods, the method including: determining reference evapotranspiration data based on evapotranspiration input data; the evapotranspiration input data is obtained from the area to be predicted; determining the relationship between soil water potential and soil moisture stress coefficient based on evapotranspiration measured data and the reference evapotranspiration data; the evapotranspiration measured data is obtained from the area to be predicted; based on the reference evapotranspiration data and the relationship, predicting the crop evapotranspiration of the area to be predicted to obtain a crop evapotranspiration prediction result.

[0055] On the other hand, the present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, is implemented to execute the crop evapotranspiration prediction method provided by the above-mentioned methods, the method comprising: determining reference evapotranspiration data based on evapotranspiration input data; the evapotranspiration input data is obtained from the area to be predicted; determining the relationship between soil water potential and soil moisture stress coefficient based on the evapotranspiration measured data and the reference evapotranspiration data; the evapotranspiration measured data is obtained from the area to be predicted; based on the reference evapotranspiration data and the relationship, predicting the crop evapotranspiration of the area to be predicted to obtain a crop evapotranspiration prediction result.

[0056] The device embodiments described above are merely illustrative, wherein the units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the scheme of this embodiment. Ordinary technicians in this field can understand and implement it without paying creative labor.

[0057] Through the description of the above implementation methods, those skilled in the art can clearly understand that each implementation method can be implemented by means of software plus a necessary general hardware platform, and of course, can also be implemented by hardware. Based on this understanding, the above technical solution is essentially or the part that contributes to the prior art can be embodied in the form of a software product, and the computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, a disk, an optical disk, etc., including a number of instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to execute the methods described in each embodiment or some parts of the embodiments.

[0058] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A crop evapotranspiration prediction method, characterized in that: include: Determining reference evapotranspiration data based on evapotranspiration input data; The evapotranspiration input data is obtained from the area to be predicted; Determining the relationship between soil water potential and soil water stress coefficient based on the evapotranspiration measured data and the reference evapotranspiration data; The evapotranspiration measured data is obtained from the area to be predicted; Based on the reference evapotranspiration data and the relationship, the crop evapotranspiration of the area to be predicted is predicted to obtain a crop evapotranspiration prediction result.

2. The crop evapotranspiration prediction method according to claim 1, characterized in that: Determining reference evapotranspiration data based on evapotranspiration input data comprises: Obtain evapotranspiration input data; the evapotranspiration input data includes crop surface net radiation , soil heat flux , average temperature , the slope of the relationship between saturated water vapor pressure and temperature , wind speed , saturated vapor pressure and the actual vapor pressure ; Based on the evapotranspiration input data and formula Determine the reference evapotranspiration data; where, is a reference evapotranspiration rate included in the reference evapotranspiration data; is a constant; a, b, c and d are parameters.

3. The crop evapotranspiration prediction method according to claim 2, characterized in that: Determining the relationship between soil water potential and soil water stress coefficient based on the evapotranspiration measured data and the reference evapotranspiration data includes: Determine the crop coefficient; Based on the crop coefficient and formula Determine the soil water stress coefficient, where is the measured evapotranspiration rate included in the evapotranspiration measured data; is the soil moisture stress coefficient; is the crop coefficient; The relationship between soil water potential and the soil water stress coefficient is determined by symbolic regression based on genetic programming.

4. The crop evapotranspiration prediction method according to claim 3, characterized in that: The relationship between soil water potential and soil water stress coefficient determined by the symbolic regression based on genetic programming includes: Determine the initial population and fitness function; Determining the fitness of each individual in the initialized population based on the fitness function; Updating the initialized population based on the fitness of each individual in the initialized population to obtain an updated population; When the renewal population meets the termination condition, the relationship between the soil water potential and the soil water stress coefficient is determined based on the fitness of each individual of the renewal population.

5. The crop evapotranspiration prediction method according to claim 4, characterized in that: The updating of the initialized population based on the fitness of each individual in the initialized population to obtain an updated population comprises: Based on the fitness of each individual in the initialized population, an update operation is performed on the initialized population to obtain an updated population; the update operation includes a selection operation, a crossover operation and a mutation operation; wherein, a selection operation is performed on the initialized population to obtain parent individuals; Performing a crossover operation on the parent individuals to obtain new offspring individuals; A mutation operation is performed on the individual part of the initialized population to obtain new individuals.

6. The crop evapotranspiration prediction method according to claim 3, characterized in that: The step of predicting the crop evapotranspiration of the area to be predicted based on the reference evapotranspiration data and the relationship to obtain a crop evapotranspiration prediction result includes: Based on the formula The crop evapotranspiration of the area to be predicted is predicted to obtain a crop evapotranspiration prediction result; wherein, is the actual evapotranspiration; is the soil water stress coefficient determined based on the relationship; The actual evapotranspiration is evaluated based on the regression prediction model evaluation index to obtain a crop evapotranspiration prediction result including the evaluation result.

7. A crop evapotranspiration prediction device, characterized in that: include: A reference evapotranspiration data determination module, used for determining reference evapotranspiration data based on evapotranspiration input data; The evapotranspiration input data is obtained from the area to be predicted; A relationship determination module, used to determine the relationship between soil water potential and soil water stress coefficient based on the evapotranspiration measured data and the reference evapotranspiration data; The evapotranspiration measured data is obtained from the area to be predicted; The crop evapotranspiration prediction module is used to predict the crop evapotranspiration of the area to be predicted based on the reference evapotranspiration data and the relationship to obtain a crop evapotranspiration prediction result.

8. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the computer program, the crop evapotranspiration prediction method according to any one of claims 1 to 6 is implemented.

9. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the crop evapotranspiration prediction method according to any one of claims 1 to 6 is implemented.

10. A computer program product, comprising a computer program, characterized in that When the computer program is executed by a processor, the crop evapotranspiration prediction method according to any one of claims 1 to 6 is implemented.

Images