A method, apparatus and equipment for dynamic simulation of lettuce root zone

Abstract

The present invention relates to the field of root domain simulation technology, and specifically to a method, device and equipment for dynamic simulation of lettuce root domains, wherein the method first solves the root domain boundary gradient angle through specific formulas and parameters to determine the boundary coordinates, and then uses the lettuce root domain development calculation equation group to comprehensively consider the dynamic changes of accumulated temperature and growth coefficient based on the root domain background matrix, accumulated temperature correlation matrix and growth development coefficient matrix, and obtains the root domain development boundary and area. The application has a small amount of calculation and easy-to-obtain indicators, and can accurately present the root domain status of lettuce, which provides a scientific basis for lettuce planting management, is conducive to precise irrigation, fertilization and other operations, and has positive significance for improving the efficiency of lettuce planting and promoting the development of refined agriculture.

Description

Technical Field

[0001] The invention relates to the technical field of root zone simulation, and in particular to a method, a device and equipment for dynamic simulation of lettuce root zones. Background Art

[0002] Existing technologies often simulate lettuce roots based on fixed parameter models. First, the initial soil conditions of the lettuce roots are measured in detail, such as accurately analyzing the soil texture, measuring the ratio of sand, silt and clay, and determining the initial nutrient content, covering basic data such as the specific values ​​of major elements such as nitrogen, phosphorus, potassium, and some trace elements. Then, a lettuce root growth model is constructed based on empirical formulas. It is usually assumed that the root growth rate is constant, such as setting the root system to grow a fixed length value every day, and extending in the soil space according to a preset direction and fixed branching angle.

[0003] In addition, existing technologies are mainly aimed at simulating the depth and width of the root system, and the simulation of lateral roots and taproots, and few studies focus on the root area. For the simulation of root length, models are mostly constructed based on empirical formulas. The growth rate and depth extension law of the taproot are set according to the plant type, and the lateral root growth is generated from the taproot at fixed branch angles and intervals. Some studies use mathematical models, such as the L system, to describe the changes in the root topology through character rewriting rules and simulate branch morphology. In recent years, with the development of computer technology, numerical methods such as finite element analysis have emerged, which can divide soil space and simulate the growth direction of roots under different soil conditions and the impact of resource absorption on growth. However, they still face problems such as insufficient accuracy in complex environments and the need to improve the universality of the model. Summary of the invention

[0004] In view of this, the object of the present invention is to provide a method, device and equipment for dynamic simulation of lettuce root zone to overcome the problems existing in the current prior art.

[0005] To achieve the above objectives, the present invention adopts the following technical solutions:

[0006] In a first aspect, the present application provides a method for simulating the dynamics of a lettuce root zone, comprising:

[0007] Get the current root domain parameters of Lettuce;

[0008] According to the lettuce root zone parameters, the root zone boundary gradient angle solving formula is used to determine the root zone boundary coordinates of the current lettuce;

[0009] Obtain the root domain background matrix, the lettuce accumulated temperature relative index column matrix, the phenological period accumulated temperature threshold column matrix, the lettuce root domain lateral growth and development coefficient column matrix, and the lettuce root domain lateral growth and development coefficient column matrix;

[0010] According to the root domain boundary coordinates, the root domain background matrix, the lettuce accumulated temperature relative index column matrix, the phenological period accumulated temperature threshold column matrix, the lettuce root domain lateral growth and development coefficient column matrix and the lettuce root domain lateral growth and development coefficient column matrix, the current lettuce root domain is dynamically simulated through the lettuce root domain development calculation equation group to determine the current lettuce root domain development boundary and area.

[0011] Furthermore, in the method described above, the lettuce root zone parameters include: root zone lateral position coordinates, root zone longitudinal position coordinates and root zone boundary initial curvature.

[0012] Furthermore, the method described above, wherein the coordinates of the current root domain boundary of the lettuce are determined by solving a root domain boundary gradient angle formula according to the lettuce root domain parameters, comprises:

[0013] Solving the root zone parameters of the lettuce by using a root zone boundary gradient angle solution formula, and obtaining a solution result;

[0014] According to the solution results, iterative calculation is performed through a set of iterative equations to obtain a result matrix;

[0015] According to the result matrix, the boundary coordinates of the root domain of the current lettuce are determined.

[0016] Further, the method described above, according to the root domain boundary coordinates, the root domain background matrix, the lettuce accumulated temperature relative index column matrix, the phenological period accumulated temperature threshold column matrix, the lettuce root domain lateral growth and development coefficient column matrix and the lettuce root domain lateral growth and development coefficient column matrix, through the lettuce root domain development calculation equation group, dynamically simulates the current lettuce root domain to determine the root domain development boundary and area of ​​the current lettuce, including:

[0017] Determine the root system relative development rate coefficient through the lettuce root domain development calculation equation group according to the root domain boundary coordinates, the root domain background matrix, the lettuce accumulated temperature relative index column matrix, the phenological period accumulated temperature threshold column matrix, the lettuce root domain lateral growth development coefficient column matrix and the lettuce root domain lateral growth development coefficient column matrix;

[0018] Determining the root system growth position and daily growth position coordinates according to the root system relative development rate coefficient;

[0019] Determine a root zone development boundary coordinate set according to the root system growth position and daily growth position coordinates;

[0020] The root zone development area of ​​the current lettuce is determined according to the root zone development boundary coordinate set and the root zone area conversion constant.

[0021] In a second aspect, the present application provides a lettuce root zone dynamic simulation device, comprising:

[0022] A data acquisition module is used to obtain the current lettuce root domain parameters;

[0023] A root zone boundary gradient angle solving module is used to determine the root zone boundary coordinates of the current lettuce according to the root zone parameters of the lettuce and the root zone boundary gradient angle solving formula;

[0024] The data acquisition module is also used to obtain the root domain background matrix, the lettuce accumulated temperature relative index column matrix, the phenological period accumulated temperature threshold column matrix, the lettuce root domain lateral growth and development coefficient column matrix and the lettuce root domain lateral growth and development coefficient column matrix;

[0025] The lettuce root domain simulation module is used to dynamically simulate the current lettuce root domain according to the root domain boundary coordinates, the root domain background matrix, the lettuce accumulated temperature relative index column matrix, the phenological period accumulated temperature threshold column matrix, the lettuce root domain lateral growth and development coefficient column matrix and the lettuce root domain lateral growth and development coefficient column matrix, through the lettuce root domain development calculation equation group, to determine the root domain development boundary and area of ​​the current lettuce.

[0026] In a third aspect, the present application provides a lettuce root domain dynamic simulation device, including a processor and a memory, wherein the processor is connected to the memory:

[0027] Wherein, the processor is used to call and execute the program stored in the memory;

[0028] The memory is used to store the program, and the program is at least used to execute any one of the above methods for simulating the dynamics of the lettuce root zone.

[0029] The beneficial effects of the present invention are:

[0030] This application first solves the root domain boundary gradient angle through specific formulas and parameters to determine the boundary coordinates. Then, based on the root domain background matrix, accumulated temperature correlation matrix, and growth and development coefficient matrix, the lettuce root domain development calculation equation group is used to comprehensively consider the dynamic changes of accumulated temperature and growth coefficient to obtain the root domain development boundary and area. This application has a small amount of calculation and easy-to-obtain indicators, and can accurately present the root domain status of lettuce. This provides a scientific basis for lettuce planting management, is conducive to precise irrigation, fertilization and other operations, and has positive significance for improving the efficiency of lettuce planting and promoting the development of refined agriculture. BRIEF DESCRIPTION OF THE DRAWINGS

[0031] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only 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.

[0032] Figure 1 It is a flow chart provided by an embodiment of a method for dynamic simulation of lettuce root domain of the present invention;

[0033] Figure 2 It is a structural schematic diagram provided by an embodiment of a lettuce root zone dynamic simulation device of the present invention;

[0034] Figure 3 The present invention is a structural schematic diagram of a lettuce root zone dynamic simulation device provided in an embodiment. DETAILED DESCRIPTION

[0035] To make the purpose, technical solution and advantages of the present invention clearer, the technical solution of the present invention will be described in detail below. Obviously, the described embodiments are only part of the embodiments of the present invention, rather than all the embodiments. Based on the embodiments of the present invention, all other implementation methods obtained by ordinary technicians in this field without creative work belong to the scope of protection of the present invention.

[0036] Figure 1 This is a flow chart provided by an embodiment of a method for dynamically simulating the root zone of lettuce according to the present invention. Figure 1 , this embodiment may include the following steps:

[0037] S1. Get the current lettuce root domain parameters;

[0038] S2. According to the root domain parameters of the lettuce, the root domain boundary coordinates of the current lettuce are determined by solving the root domain boundary gradient angle formula;

[0039] S3, obtaining a root domain background matrix, a lettuce accumulated temperature relative index column matrix, a phenological period accumulated temperature threshold column matrix, a lettuce root domain lateral growth and development coefficient column matrix, and a lettuce root domain lateral growth and development coefficient column matrix;

[0040] S4. According to the root domain boundary coordinates, root domain background matrix, lettuce accumulated temperature relative index column matrix, phenological period accumulated temperature threshold column matrix, lettuce root domain lateral growth and development coefficient column matrix and lettuce root domain lateral growth and development coefficient column matrix, the current lettuce root domain is dynamically simulated through the lettuce root domain development calculation equation group to determine the current lettuce root domain development boundary and area.

[0041] It is understandable that in solving the gradient angle of the root domain boundary, the present application uses the horizontal and vertical position coordinates of the root domain and the initial curvature of the boundary as the basis, combined with the accuracy of the step matrix, and uses the cyclic iterative equation group for calculation. By setting the initial value and following a specific iteration method, the result of each cycle is recorded in the result storage matrix, and the root domain boundary coordinates are obtained after a series of complex calculations. This process lays a solid foundation for subsequent root domain simulation through precise mathematical models and rigorous calculation logic.

[0042] In the simulation of lettuce boundaries and root domains, several key matrices were introduced. The root domain background matrix constructs the basic framework of the simulation, while the lettuce accumulated temperature relative index column matrix reflects the impact of daily temperature and other data on lettuce growth, and the phenological accumulated temperature threshold column matrix defines the accumulated temperature requirements for lettuce at different growth stages. At the same time, the lettuce root domain lateral and longitudinal growth and development coefficient column matrices further reflect the growth characteristics of the root domain in different directions. With the help of the lettuce root domain development calculation equation group, this patent organically combines multiple factors such as the change in the relative index of accumulated temperature, the dynamic adjustment of the growth and development coefficient, and the change in the phenological accumulated temperature threshold. In the calculation process, the horizontal and vertical relative development rate coefficients of the root system, as well as parameters such as the daily longitudinal growth position coordinates are accurately considered to accurately determine the root domain development boundary and area.

[0043] From the perspective of computational characteristics, the amount of computation is relatively small, which has obvious advantages over traditional complex simulation technologies that consume large computing resources. In addition, the computational indicators are easy to obtain, which lowers the threshold for technology application. From the perspective of practical application in agricultural production, accurate root domain simulation results provide all-round support for lettuce planting. In terms of irrigation, the reasonable amount and frequency of irrigation water can be determined based on the root domain area, boundary and growth stage to avoid the adverse effects of water waste and insufficient irrigation on lettuce growth. In the fertilization stage, the type and amount of fertilizer can be accurately adjusted according to the development status of the root domain and the absorption characteristics of lettuce for nutrients, thereby improving fertilizer utilization and reducing environmental pollution. In addition, it also has important reference value for the selection and breeding of lettuce varieties and the optimization and regulation of the planting environment, which helps agricultural researchers to deeply study the growth mechanism of lettuce, promote the development of the lettuce planting industry in a refined and scientific direction, and enhance the economic benefits and sustainable development capabilities of the entire industry.

[0044] Preferably, the lettuce root zone parameters include: the root zone lateral position coordinates, the root zone longitudinal position coordinates and the initial curvature of the root zone boundary.

[0045] Preferably, step S2 comprises:

[0046] Solve the root zone parameters of lettuce through the root zone boundary gradient angle solution formula and obtain the solution result;

[0047] According to the solution results, the iterative calculation is performed through the iterative equation group to obtain the result matrix;

[0048] Based on the result matrix, determine the root domain boundary coordinates of the current lettuce.

[0049] Understandably,

[0050]

[0051]

[0052]

[0053] are the horizontal and vertical position coordinates of the root domain, is the initial curvature of the root domain boundary (unit: radians); To solve the step matrix accuracy, the default value is 500, that is, the accuracy is 1 / 500 (radians), To solve the step size matrix.

[0054]

[0055]

[0056]

[0057] is a set of iterative equations, is the current cycle number, is the intermediate iteration variable (unit: radians) and the iteration method is: To solve the angle value (unit: angle value), The initial value is 0 and the number of cycles is (Right now, ), For the last time The calculated value of the loop; during the loop iteration, each time ( )、 The results are recorded in , The results are stored in a matrix.

[0058]

[0059]

[0060]

[0061]

[0062]

[0063] The matrix is The result matrix is ​​subtracted from the constant The values ​​after that are taken as the matrix of absolute values, The meaning is to obtain Minimum value in a matrix , The function is to retrieve The line where , The current The boundary coordinates of the lower root domain.

[0064] Preferably, step S4 comprises:

[0065] According to the root zone boundary coordinates, root zone background matrix, lettuce accumulated temperature relative index column matrix, phenological period accumulated temperature threshold column matrix, lettuce root zone lateral growth and development coefficient column matrix and lettuce root zone lateral growth and development coefficient column matrix, the root system relative development rate coefficient is determined through the lettuce root zone development calculation equation group;

[0066] According to the root relative development rate coefficient, the root growth position and daily growth position coordinates are determined;

[0067] Determine the root zone development boundary coordinate set based on the root growth position and daily growth position coordinates;

[0068] The root domain development area of ​​the current lettuce is determined based on the root domain development boundary coordinate set and the root domain area conversion constant.

[0069] Understandably,

[0070]

[0071]

[0072]

[0073]

[0074]

[0075] Root domain background matrix, are the rows and columns of the matrix respectively. is the relative index matrix of accumulated temperature of lettuce (numeric value dimensionless), which is calculated from daily temperature and other data. There are many calculation methods in the public information, so this article does not give a calculation description. is the phenological period accumulated temperature threshold column matrix (value dimensionless), is the column matrix of the lateral growth coefficient of the lettuce root zone (numeric dimensionless), is the matrix of lateral growth coefficients of lettuce root zone (value dimensionless). The above matrix contains default parameters and needs to be adjusted according to the actual characteristics of lettuce when applied.

[0076]

[0077] is the lettuce root zone development calculation equation set, is the current accumulated temperature relative index (dimensionless), is the relative index of accumulated temperature of the previous day, is the step value of the accumulated temperature index (dimensionless); They are the current lettuce root zone lateral growth and development coefficient and the previous lettuce root zone lateral growth and development coefficient respectively; They are the current lettuce root zone longitudinal growth coefficient and the previous lettuce root zone longitudinal growth coefficient respectively; are the accumulated temperature threshold of the current phenological period and the accumulated temperature threshold of the previous phenological period respectively. Z( ) equation is the “root domain boundary gradient angle solution”, which is the input required for the S2 step.

[0078] and y are the relative root growth rate coefficients in the horizontal direction and in the vertical direction (both dimensionless), and are the root growth position and daily growth position coordinates, is the coordinate set of the root domain development boundary.

[0079]

[0080]

[0081] express Collection and Background Matrix The area formed by the boundary, i.e. the root domain . is the root area conversion constant, unit: mm2, is the converted root area of ​​lettuce (mm2). In actual application, this parameter needs to be adjusted according to different varieties, such as Salanova. , When 500 respectively, About 2.1mm2.

[0082] The present invention also provides a lettuce root zone dynamic simulation device for implementing the above method embodiment. Figure 21 is a schematic diagram of a structure of a dynamic simulation device for lettuce root zone according to an embodiment of the present invention. Figure 2 As shown, this embodiment includes:

[0083] Data acquisition module 1, used to obtain the current lettuce root zone parameters;

[0084] The root zone boundary gradient angle solving module 2 is used to determine the root zone boundary coordinates of the current lettuce according to the root zone parameters of the lettuce through the root zone boundary gradient angle solving formula;

[0085] The data acquisition module 1 is also used to obtain the root domain background matrix, the lettuce accumulated temperature relative index column matrix, the phenological period accumulated temperature threshold column matrix, the lettuce root domain lateral growth and development coefficient column matrix and the lettuce root domain lateral growth and development coefficient column matrix;

[0086] The lettuce root domain simulation module 3 is used to dynamically simulate the current lettuce root domain according to the root domain boundary coordinates, root domain background matrix, lettuce accumulated temperature relative index column matrix, phenological period accumulated temperature threshold column matrix, lettuce root domain lateral growth and development coefficient column matrix and lettuce root domain lateral growth and development coefficient column matrix, and determine the current root domain development boundary and area of ​​the lettuce through the lettuce root domain development calculation equation group.

[0087] Regarding the device in the above embodiment, the specific manner in which each module performs operations has been described in detail in the embodiment of the method, and will not be elaborated here.

[0088] The present invention also provides a lettuce root zone dynamic simulation device for implementing the above method embodiment. Figure 3 This is a schematic diagram of the structure of an embodiment of a lettuce root zone dynamic simulation device provided by the present invention. Figure 3 As shown, the lettuce root domain dynamic simulation device of this embodiment includes a processor 21 and a memory 22, and the processor 21 is connected to the memory 22. The processor 21 is used to call and execute the program stored in the memory 22; the memory 22 is used to store the program, and the program is at least used to execute the lettuce root domain dynamic simulation method in the above embodiment.

[0089] The specific implementation scheme of the lettuce root zone dynamic simulation device provided in the embodiment of the present application can refer to the implementation scheme of the lettuce root zone dynamic simulation method in any of the above embodiments, which will not be repeated here.

[0090] It can be understood that the same or similar parts of the above embodiments can be referenced to each other, and the contents not described in detail in some embodiments can refer to the same or similar contents in other embodiments.

[0091] It should be noted that, in the description of the present invention, the terms "first", "second", etc. are only used for descriptive purposes and cannot be understood as indicating or implying relative importance. In addition, in the description of the present invention, unless otherwise specified, the meaning of "plurality" refers to at least two.

[0092] Any process or method description in a flowchart or otherwise described herein may be understood to represent a module, segment or portion of code that includes one or more executable instructions for implementing the steps of a specific logical function or process, and the scope of the preferred embodiments of the present invention includes alternative implementations in which functions may not be performed in the order shown or discussed, including performing functions in a substantially simultaneous manner or in the reverse order depending on the functions involved, which should be understood by those skilled in the art to which the embodiments of the present invention belong.

[0093] It should be understood that the various parts of the present invention can be implemented by hardware, software, firmware or a combination thereof. In the above-mentioned embodiments, multiple steps or methods can be implemented by software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented by hardware, as in another embodiment, it can be implemented by any one of the following technologies known in the art or their combination: a discrete logic circuit having a logic gate circuit for implementing a logic function for a data signal, a dedicated integrated circuit having a suitable combination of logic gate circuits, a programmable gate array (PGA), a field programmable gate array (FPGA), etc.

[0094] A person skilled in the art may understand that all or part of the steps in the method for implementing the above-mentioned embodiment may be completed by instructing related hardware through a program, and the program may be stored in a computer-readable storage medium, which, when executed, includes one or a combination of the steps of the method embodiment.

[0095] In addition, each functional unit in each embodiment of the present invention may be integrated into a processing module, or each unit may exist physically separately, or two or more units may be integrated into one module. The above-mentioned integrated module may be implemented in the form of hardware or in the form of a software functional module. If the integrated module is implemented in the form of a software functional module and sold or used as an independent product, it may also be stored in a computer-readable storage medium.

[0096] The storage medium mentioned above can be a read-only memory, a magnetic disk or an optical disk, etc.

[0097] In the description of this specification, the description with reference to the terms "one embodiment", "some embodiments", "examples", "specific examples", or "some examples" means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present invention. In this specification, the schematic representation of the above terms does not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described may be combined in any one or more embodiments or examples in a suitable manner.

[0098] Although the embodiments of the present invention have been shown and described above, it is to be understood that the above embodiments are exemplary and are not to be construed as limitations of the present invention. A person skilled in the art may change, modify, replace and vary the above embodiments within the scope of the present invention.

Claims

1. A method for dynamic simulation of lettuce root zone, characterized in that: include: Obtaining current lettuce root zone parameters; wherein the lettuce root zone parameters include: root zone horizontal position coordinates, root zone vertical position coordinates and root zone boundary initial curvature; Solving the root zone parameters of the lettuce by using a root zone boundary gradient angle solution formula, and obtaining a solution result; According to the solution results, iterative calculation is performed through a set of iterative equations to obtain a result matrix; According to the result matrix, determining the boundary coordinates of the root domain of the current lettuce; Obtain the root domain background matrix, the lettuce accumulated temperature relative index column matrix, the phenological period accumulated temperature threshold column matrix, the lettuce root domain lateral growth and development coefficient column matrix, and the lettuce root domain longitudinal growth and development coefficient column matrix; Determine the root system relative development rate coefficient through the lettuce root domain development calculation equation group according to the root domain boundary coordinates, the root domain background matrix, the lettuce accumulated temperature relative index column matrix, the phenological period accumulated temperature threshold column matrix, the lettuce root domain lateral growth and development coefficient column matrix and the lettuce root domain longitudinal growth and development coefficient column matrix; Determining the root system growth position and daily growth position coordinates according to the root system relative development rate coefficient; Determine a root zone development boundary coordinate set according to the root system growth position and daily growth position coordinates; The root zone development area of ​​the current lettuce is determined according to the root zone development boundary coordinate set and the root zone area conversion constant.

2. A lettuce root zone dynamic simulation device, characterized in that: include: A data acquisition module is used to acquire the current lettuce root zone parameters; wherein the lettuce root zone parameters include: the root zone horizontal position coordinates, the root zone vertical position coordinates and the root zone boundary initial curvature; A root domain boundary gradient angle solving module is used to solve the root domain parameters of the lettuce through the root domain boundary gradient angle solving formula and obtain the solution result, and according to the solution result, perform iterative calculation through a cyclic iterative equation group to obtain a result matrix, and determine the root domain boundary coordinates of the current lettuce according to the result matrix; The data acquisition module is also used to obtain the root domain background matrix, the lettuce accumulated temperature relative index column matrix, the phenological period accumulated temperature threshold column matrix, the lettuce root domain lateral growth and development coefficient column matrix and the lettuce root domain longitudinal growth and development coefficient column matrix; The lettuce root domain simulation module is used to determine the root relative development rate coefficient according to the root domain boundary coordinates, the root domain background matrix, the lettuce accumulated temperature relative index column matrix, the phenological period accumulated temperature threshold column matrix, the lettuce root domain lateral growth and development coefficient column matrix and the lettuce root domain longitudinal growth and development coefficient column matrix, through the lettuce root domain development calculation equation group, determine the root system growth position and daily growth position coordinates according to the root system relative development rate coefficient, determine the root system growth position and daily growth position coordinates, determine the root domain development boundary coordinate set according to the root system growth position and daily growth position coordinates, and determine the current root domain development area of ​​the lettuce according to the root domain development boundary coordinate set and the root domain area conversion constant.

3. A lettuce root zone dynamic simulation device, characterized in that: The invention comprises a processor and a memory, wherein the processor is connected to the memory: Wherein, the processor is used to call and execute the program stored in the memory; The memory is used to store the program, and the program is at least used to execute the lettuce root zone dynamic simulation method described in claim 1.

Images