This application is based upon and claims priority to Chinese Patent Application No. 202111088757.7, filed on Sep. 16, 2021, the entire contents of which are incorporated herein by reference
The present invention relates to the technical field of plant growth analysis, and in particular to a dynamic growth model construction method and system for a crop canopy structure, and a storable medium.
The crop model quantitatively describes the complex interaction relationship of elements such as meteorological factors, soil moisture, nutrients and temperature in the crop ecosystem in the processes of crop growth, development and the like by the description of the characteristics of all components in the farmland ecosystem and the influence of the components on the system, and provides an important means for optimizing cultivation, reasonable moisture supply, reasonable resource utilization, planting benefit optimization and the like. In recent years, with the development of computer technology and the continuous improvement of computer performance, scholars at home and abroad have made great progress in the research and application of crop models.
Currently, at least 250 crop simulation models have been developed based on incomplete statistics. According to the development characteristics of the crop models, the crop models can be divided into the following categories: category one is a mechanism model, such as a computer model of maize photosynthetic production constructed by de Wit, a Dutch scholar, a model of influence of maize leaf area and leaf blade angle on group photosynthesis developed by Duncan et al., an ELCROS (Elementary Crop Simulator) model established by the school of de Wit based on system dynamics as a theoretical basis, and a C3 plant photosynthetic biochemical model (referred to as FvCB model for short) provided by Farquhar et al., and a C4 plant photosynthetic model established by colatz et al., which can simulate photosynthetic biochemical reactions inside leaf blades; category two is an empirical model, in order to improve the applicability of the crop model, the crop mechanism model is connected with multidisciplinary, the complex crop model is subjected to parameter simplification, such as a CERES-Maize model and a GROPGRO-Soybean model, the connection is established between a crop growth model and a geographic information system ArcView, and the spatial-temporal change of soybeans and maize on the field scale can be simulated; the United States WGEN and the United Kingdom LARS-WG establish the combination of a stochastic weather generator and a crop model so as to form an agricultural production decision system; the Lobell team in the United States proposed an SCYM (Scalable Crop Yield Mapper) framework that couples multispectral satellite remote sensing data to crop growth estimates; and category three is a comprehensive model, which can not only simulate the growth, the development, the dry matter, the yield and the like of crops, but also simulate the habitat conditions of the crops, and there are early models, for example, the BACROS model (basic crop simulation) and the SUCROS (simple and universal crop simulator) model developed by the team from Wageningen University & Research in the Netherlands on the basis of the ELCROS model; the SUCROS model is based on the simulation of the photosynthesis of the crop canopies, and can be applied to different environmental conditions and different crop species only by changing the crop parameters. Currently, the comprehensive models widely used include a WOFOST (World Food Studies) model developed by Hijmans, an APSIM (Agricultural Production Systems slMulator) model developed by CSIRO of Australia, a DSSAT (Decision Support System for Agrotechnology Transer) model developed by University of Hawaii, an AquaCrop model developed by the Food and Agriculture Organization of the United Nations, and the like. However, most of the current crop models are one-dimensional vertical models of a soil-crop-atmosphere system, and characterize the conditions of crop canopies by mainly adopting parameters such as leaf area indexes and canopy coverage. Due to insufficient attention on the characteristics of crop canopy structure, the simulation error of the hydrothermal carbonization processes inside the system is large and the universality of the model is insufficient, and particularly the existing crop models generate remarkable deviation due to remarkably change of the canopy structure caused by the adjusted planting mode.
Although some scholars implement the expansion of the crop model in the regional scale by means of the GIS technology, the model is actually “pseudo-three-dimensional”, because this model cannot solve the problem of uneven distribution of components in the spatial ecosystem, for example, the problem of competition of different crop growth (including illumination, moisture, nutrients, land resources and the like) and horizontal migration of substances (including moisture, nutrients, pesticides and the like) on the region. Therefore, the establishment of a three-dimensional spatial crop growth simulation system is an important direction for breaking through the one-dimensional spatial limitation of a crop model and expanding the crop model to a region, and is a basis for the application of the model expanding to the regional scale.
With respect to the limitations of one-dimensional models, numerous scholars have attempted to characterize the canopy structures. Greyson et al. studied the relationship between leaf blade width and leaf position; Hesketh et al. established the relationship between leaf blades, leaf sheaths and internodes; Girardin et al. studied the azimuth angle distribution law of leaf blade of maize and the like. Although many scholars have established models of spatial crop structures, these models are mostly static simulations, and thus have not been coupled with dynamic crop growth models. Researches find that the crop canopy structure finally causes different photosynthetic yields by influencing the light interception quantity inside the canopy, and the photosynthetic yields influence the development and morphogenesis of each organ of a plant and finally cause the change of the canopy structure. Therefore, a crop model is urgently needed to fully reflect the complex physiological and ecological change process.
In view of this, the present invention provides a dynamic growth model construction method and system for a crop canopy structure, and a storable medium. Through full research of the characteristics of a crop canopy structure (canopy type, leaf type, spike type, internode, leaf sheath and the like) and the variation laws in the crop development processes, the present invention establishes a three-dimensional crop canopy structure model, establishes a dynamic maize growth simulation by fusing the three-dimensional crop canopy structure model, and finally establishes a characteristic platform of a phenotype structure of a three-dimensional spatial crop canopy structure, thus realizing dynamic simulation of a water-carbon process of a soil-crop-atmospheric system, and effectively improving the simulation capability and universality of the crop model.
In order to achieve the above objective, the present invention adopts the following technical solutions:
Provided is a dynamic crop growth model construction method for a three-dimensional canopy structure, comprising the following steps:
Optionally, the method further comprises establishing a dynamic model of a water-carbon process of a soil-crop-atmospheric system according to the three-dimensional dynamic crop growth model, specifically as follows:
On the basis of quantifying a three-dimensional dynamic single-plant growth simulation model, constructing a deduction method for crop leaf blades to a single plant canopy and then to a spatial canopy structure of a crop group, constructing a database based on data, realizing topological expansion of the spatial structure, and meanwhile, integrating a time dimension, and constructing a crop water-carbon flux multi-spatial-temporal cooperative coupling quantitative characterization method.
Optionally, the crop parameters comprise canopy structure characteristics and variation laws in crop development processes; the canopy structure characteristics comprise: a canopy type, a leaf type, a spike type, an internode, and a leaf sheath.
Optionally, the two-dimensional canopy structure model describes crop growth and development stages by adopting growing degree-days, and carries out dynamic simulation on a two-dimensional crop canopy structure by adopting a Logistic equation.
Optionally, the step of constructing a three-dimensional dynamic crop growth simulation model specifically comprises the following steps: on the basis of the two-dimensional canopy structure simulation model, characterizing veins and leaf margins by adopting a spatial topological structure, and constructing a three-dimensional dynamic growth simulation model.
Provided is a dynamic crop growth model construction system for a three-dimensional canopy structure, comprising:
Optionally, the system further comprises a dynamic model establishment module of a water-carbon process of a soil-crop-atmospheric system used for establishing a dynamic model of a water-carbon process of a soil-crop-atmospheric system according to the three-dimensional dynamic crop growth model.
Provided is a computer storage medium having a computer program stored thereon, wherein the computer program, when being executed by a processor, implements the steps of the dynamic crop growth model construction method for a three-dimensional canopy structure.
According to the above-mentioned technical solutions, compared with the prior art, the dynamic crop growth model construction method and system and the storable medium disclosed herein solve the problem that the deviation is generated due to the fact that the existing crop models characterize the crop canopy conditions by only adopting parameters such as leaf area indexes and canopy coverage, and particularly the precision of the crop growth simulation is greatly improved after the canopy structure is remarkably changed due to the adjusted planting mode, and the universality of the crop model is greatly improved. In addition, the present invention provides a novel means for quantitative precision agriculture in the future, and has a wide prospect for understanding the interrelation among all components of a farmland ecosystem, improving the utilization efficiency of agricultural irrigation water and realizing the expansion and application of a model on a regional scale.
In order to more clearly illustrate the technical solutions in the embodiments of the present invention or in the prior art, the drawings required to be used in the description of the embodiments or the prior art are briefly introduced below. It is obvious that the drawings in the description below are some embodiments of the present invention, and those of ordinary skilled in the art can obtain other drawings according to the drawings provided herein without creative efforts.
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention but not all of them. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skilled in the art without creative efforts fall within the protection scope of the present invention.
An embodiment of the present invention discloses a dynamic crop growth model construction method for a three-dimensional canopy structure, which comprises the following steps:
The method further comprises establishing a dynamic model of a water-carbon process of a soil-crop-atmospheric system according to the three-dimensional dynamic crop growth model, specifically as follows:
On the basis of quantifying a three-dimensional dynamic single-plant growth simulation model, constructing a deduction method for crop leaf blades to a single plant canopy and then to a spatial canopy structure of a crop group, constructing a database based on data, realizing topological expansion of the spatial structure, and meanwhile, integrating a time dimension, and constructing a crop water-carbon flux multi-spatial-temporal cooperative coupling quantitative characterization method.
Specifically, the method is as follows:
The crop model usually describes crop growth and development stages by adopting growing degree-days (GDDs), and if TMAX and TMIN meet the following conditions, then a calculation formula of GDD is as follows:
wherein To is an optimum temperature for maize growth, Tb is a base temperature for maize growth, and TMAX and TMIN are a daily maximum temperature and a daily minimum temperature, respectively.
If the daily maximum temperature and the daily minimum temperature do not meet the above-mentioned conditions, then the influence of the day-night temperature difference on the development of crops shall be considered. This means that the influence of the day-night temperature difference on the crops needs to be considered. The following equations can be used:
wherein Tfac(i) (i = 1, 2, 3,......24) is a temperature change factor of each period, T(i) is a temperature of each period, Tm is a maximum temperature of the crop growth period, and SUMGDD is a sum of GDDs of each development period.
The maize development period is divided into 7 stages according to the genetic characteristics and environmental factors of maize. Quantitative development of thermal effect, photoperiod effect and genetic effect is considered, and a unified physiological development time scale is established. The maize development period is as follows:
The simulation of the entire maize development stage ends when the cumulative GDDs reach P5 at stages 5, 6 and 7. S7 is the maize mature stage. The parameters of the maize development period are shown in Table 1.
Research shows that the elongation of the internodes of maize is a process from slow to fast to slow along with the change of GDD, the internodes grow in an upward elongation way from the morphological lower end, and the growth has certain sequentiality and overlapping property; the thickness of each internode of maize gradually increases along with the growth of maize. Thus, the process of internode elongation and thickening in maize can be quantitatively described by equations (13, 14):
wherein INLen (GDD) and INWen (GDD) are the length (cm) and width (cm) of the nth internode at the time of GDD, respectively; INGDDn is GDD of the nth internode which begins to elongate and thicken, and INPa and INPb are crop genetic parameters and are used for calibrating the model; s1, s2, m1 and m2 are model parameters and are used for calibrating the model and are related to the characteristics of the variety; BUMDTT is GDD of maize at the jointing stage; a1, b1 and c1 are model parameters obtained when a potential maximum value is taken as the internode lengthand are related to the characteristics of the variety; a2, b2 and c2 are model parameters obtained when a potential maximum value is taken as the internode width and are related to the characteristics of the variety.
The elongation process of the leaf sheath of maize is a process from slow to fast to slow, and the growth curve of the leaf sheath is an S-shaped growth curve. Therefore, the present study adopts Logistic equations (18, 19) to describe the maize leaf sheath growth process:
wherein LSLen (GDD) and LSWen (GDD) are the length (cm) and width (mm) of the nth leaf sheath at the time of GDD, respectively; LSLen and LSWen are potential maximum values (cm) of the length and width of the nth leaf sheath, respectively; LSGDDn is the initial GDD when the nth leaf sheath begins to elongate or thicken; GDD1 is GDD required by maize from sowing to seedling emergence, and LSPa and LSPb are crop genetic parameters and are used for calibrating the model; s3, s4, m3 and m4 are model parameters and are used for calibrating the model and are related to the characteristics of the variety; a3, b3 and c3 are model parameters obtained when a potential maximum value is taken as the length of the leaf sheath and are related to the characteristics of the variety; a4, b4 and c4 are model parameters obtained when a potential maximum value is taken as the width of the leaf sheath and are related to the characteristics of the variety.
The elongation process of the leaf blade of maize is a process from slow to fast along with the change of GDD, the growth curve of the leaf blade is an S-shaped increase, namely, a Logistic equation, so that the elongation process of the leaf blade of maize can be quantitatively described by equations (23 and 24):
wherein LELen (GDD) and LEWen (GDD) are the length (cm) and width (mm) of the nth leaf sheath at the time of GDD, respectively; LEGDDn is the initial GDD when the nth leaf sheath begins to elongate or thicken; GDD1 is GDD required by maize from sowing to seedling emergence, and LEPa and LEPb are crop genetic parameters and are used for calibrating the model; s5, s6, m5 and m6 are model parameters and are used for calibrating the model and are related to the characteristics of the variety; a5, b5 and c5 are model parameters obtained when a potential maximum value is taken as the length of the leaf blade and are related to the characteristics of the variety; a6, b6 and c6 are model parameters obtained when a potential maximum value is taken as the width of the leaf blade and are related to the characteristics of the variety.
The maize tassel begins to grow after the maize enters into the jointing stage, the length of the tassel is mainly the elongation of the internode below the tassel, and the elongation and the thickening of the maize tassel can be quantitatively described by equations (28, 29):
wherein TALen (GDD) and TAWen (GDD) are the length (cm) and width (mm) of tassel at the time of GDD, respectively; TALemax and TAWemax are the potential length and potential width of the maize tassel, respectively; TAGDDn is the initial GDD when the tassel begins to elongate or thicken; TAB is the GDD required by the growth of the maize tassel, and TAPa and TAPb are crop genetic parameters and are used for calibrating the model; m7 and m8 are model parameters obtained when potential maximum values are taken as the length and width of the tassel, respectively and are related to the characteristics of the variety.
The growth of the maize silking is slightly slower than that of the tassel, and when the heading stage ends, the silking growth is rapid, and the elongation and thickening of the maize silking can be quantitatively described by equations (30, 31):
wherein SLen (GDD) and SWen (GDD) are the length (cm) and width (mm) of the silking at the time of GDD, respectively; SLemax and SWemax are the potential length and the potential width of the maize silking, respectively; SGDDn is the initial GDD when the silking begins to elongate or thicken; SB is the GDD required by the growth of the maize silking, and Spaa and SPb are crop genetic parameters and are used for calibrating the model; m9 and m10 are model parameters and are related to the characteristics of the variety.
In the process of describing the three-dimensional spatial structure of maize, the spatial variation characteristics of each organ of the maize need to be described, such as the spatial orientation of veins, the characteristics of leaf margins, the characteristics of leaf shape growth, the morphology of leaf sheath, leaf azimuth angle and the like.
The veins run through the vascular bundles and other tissue leaf blades in the mesophyll to play the roles of conduction and support, and extend from the petiole to the main veins in the leaf blades. The veins directly determine the spatial orientation of leaf blades, and the spatial orientation of the veins is often described by means of “oblique ejection track” in research. Assuming that the vein curve is in a two-dimensional plane, then the motion track equation of the growth point of the vein curve can be deduced from the oblique ejection motion law as follows:
wherein i is any point of the growth time sequence of the veins (0 ≤ i ≤ N), and N is the total time of the “oblique ejection track” motion of the veins and is determined by the length of the veins. g is the gravitational acceleration, α is the initial inclination of the leaf blades, and h is the vertical maximum height of the veins. The vein curve is determined by the parameters a and h, and the values of a and h are determined by the growth time of the veins and change along with the change of the growth time of the veins.
The leaf margin of the leaf blade of maize consists of an unfolding part and a folding part of a leaf blade. As the leaf blades unfold, the leaf margins gradually appear wavy folds, while the leaf margins of the folding part are closed, which can be approximately seen as a curled conical surface.
Assuming that the coordinates of point P are (Xp, Yp, Zp), d=
wherein θ is the relative positions of point S, point P and each plane, and I is the leaf width of point S at point P on the vein.
The folding part of the leaf margin may be approximately a curled conical surface, the cross-section of which is a spiral. α is the polar angle of the spiral, r is the minimum radius vector, R is the maximum radius vector, and P (Xp, Yp, Zp) is the center of the spiral. Assuming that the radius vector of the spiral is in linear relation with the polar angle, the spiral is divided into n segments according to angle, so as to form a sequence {Pi} of n +1 points, and if Pi (xi, yi, zi) exists, then the following relation exists:
The point corresponding to the polar angle α/2 on the spiral is superposed with the vein, then the spiral is attached to the vein; the point on the vein is taken as the point corresponding to the polar angle α/2 on the spiral, and the point on each spiral is respectively solved by combining the relationship between the radius vector and the center of the spiral and by using the above-mentioned formula, so that the folding leaf blades can be obtained.
Research shows that the leaf blade length of maize has a significant correlation with leaf arrangement, and the maize changes in a “V” shape from the plant base to the top of the canopy.
Usually, the number of leaf blades of maize plants is determined by genetic characteristics and is related to crop varieties, and the length of leaf blades of most maize varieties is longer in the middle and shorter in both sides. The simulation of leaf blade length can be expressed by the following formula:
wherein LDLM is the maximum length of a single leaf, and NLDLM is the leaf arrangement of the longest leaf blade of maize. LDLa and LDLb are model parameters and are determined by genetic characteristics.
The variation laws in the width and length of leaf blades of maize are similar, i.e., the leaf blade is wider in the middle and slightly narrower in the upper and lower parts. The simulation of leaf blade width and leaf arrangement can be expressed by the following formula:
wherein LDWM is the maximum width of a single leaf, and NLDWM is the leaf arrangement of the widest leaf blade of maize. LDWa and LDWb are model parameters and are determined by genetic characteristics.
The leaf shape of maize is a dynamic change process from the beginning of leaf blade growth to the unfolding and shaping of the leaves, and is determined by the change of leaf width in the elongation direction. Therefore, a functional relationship between leaf width and leaf length is constructed to simulate the change in leaf shape:
wherein LDL is the length of a leaf blade, LDW is the maximum width of a leaf blade, lw is the width of a leaf at the length ll, and n1, n2 and n3 are model parameters and are related to crop varieties.
The leaf sheaths of maize are wrapped on the outer sides of the internodes, playing the roles of supporting the leaf blades, enhancing the plant stalks, and meanwhile, achieving the lodging resistance effect. Research shows that the length of the leaf sheath gradually increases with the elevation of the leaf position, while it begins to shorten at ear-leaf. The relationship between the leaf sheath length and leaf arrangement of unfolding leaves of maize can be expressed by the following formula:
wherein LN is the total number of leaves, LSLM is the maximum length of a leaf sheath between 1 to Ni of the leaf arrangement, NLSLM is the leaf arrangement of the longest leaf sheath between 1 to Ni of the leaf arrangement, and n4, n5 and n6 are model parameters and are related to crop varieties.
Internode elongation in maize is usually in the form of a slow-fast-slow law, with a regular sequence of internode elongation, the lower part of the plant elongating first and then growing gradually. The relationship between the length of the elongated internode and the leaf arrangement of maize can be expressed by the following formula:
wherein LN is the total number of leaves, INLM is the maximum length of a leaf sheath between 1 to Ni of the leaf arrangement, NINLM is the leaf arrangement of the longest leaf sheath between 1 to Ni of the leaf arrangement, and n7, n8 and n9 are model parameters and are related to crop varieties.
On the basis of quantifying a three-dimensional spatial structure of a single plant, a deduction method for crop leaf blades to a single plant canopy and then to a spatial canopy structure of a crop group is systematically constructed, the existing data are combined to construct a database, topological expansion of the spatial structure is realized, and meanwhile, a time dimension is integrated, and a crop water-carbon flux multi-spatial-temporal cooperative coupling quantitative characterization method is constructed. The specific implementation process is shown in
In the present embodiment, also included is a dynamic crop growth model construction system for a three-dimensional canopy structure, comprising:
The system further comprises a dynamic model establishment module of a water-carbon process of a soil-crop-atmospheric system used for establishing a dynamic model of a water-carbon process of a soil-crop-atmospheric system according to the three-dimensional dynamic crop growth model.
Provided is a computer storage medium having a computer program stored thereon, wherein the computer program, when being executed by a processor, implements the steps of the dynamic crop growth model construction method for a three-dimensional canopy structure.
The above description of the disclosed embodiments enables those skilled in the art to implement or use the present invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the present invention. Thus, the present invention is not intended to be limited to these embodiments shown herein but is to accord with the broadest scope consistent with the principles and novel features disclosed herein.
Number | Date | Country | Kind |
---|---|---|---|
202111088757.7 | Sep 2021 | CN | national |