A method, device and medium for determining the type of lodging of crops caused by wind

Abstract

The present invention relates to the field of agricultural production technology, and discloses a method, device and medium for determining the type of lodging of crops caused by wind disasters, including: obtaining meteorological data, soil data and crop data; determining the wind force, soil tension or crop bending stress of crops based on at least one of the meteorological data, soil data and crop data; determining whether crops have lodged and the type of lodging based on the comparison results of the wind force on crops and the soil tension and crop bending stress; wherein the meteorological data is data characterizing the influence of environmental factors on crop lodging, the soil data is data characterizing the influence of soil factors on crop lodging, and the crop data is data characterizing the influence of crops themselves on crop lodging; the lodging types include broken type lodging and non-broken type lodging. Based on the above scheme, it is possible to accurately determine whether crops have lodged and the type of lodging, which is conducive to taking targeted measures to reduce the impact of wind disasters on crops or formulating measures to save crop production.

Description

Technical Field

[0001] The present invention relates to the technical field of agricultural production, and in particular to a method, a device and a medium for determining the type of lodging of crops caused by wind. Background Art

[0002] During the growth process of crops, wind disasters, as a common natural disaster, have a huge impact on the normal growth of crops.

[0003] If we can accurately judge the damage and lodging of crops before or after a storm occurs, we can adopt certain rescue and emergency preparedness strategies in advance to avoid or reduce the impact of the storm on crop production. Summary of the invention

[0004] In order to solve the above technical problems, the present invention provides a method, device and medium for determining the type of lodging of crops caused by wind disasters.

[0005] In a first aspect, the present invention provides a method for determining the type of lodging of crops caused by wind disasters, comprising:

[0006] Obtain weather data, soil data, and crop data;

[0007] Determine, based on at least one of the meteorological data, the soil data and the crop data, whether the crop is subjected to wind force, soil tension or crop bending stress;

[0008] Based on the comparison result of the wind force on the crop, the soil tension and the crop bending stress, determining whether the crop has fallen and the type of lodging;

[0009] The meteorological data are data representing the influence of environmental factors on crop lodging, the soil data are data representing the influence of soil factors on crop lodging, and the crop data are data representing the influence of crops themselves on crop lodging;

[0010] The lodging types include breaking type lodging and non-breaking type lodging.

[0011] Optionally, the meteorological data includes air pressure, temperature and wind speed, and the crop data includes leaf thickness, leaf width, leaf length, leaf area, elastic modulus and friction coefficient;

[0012] Determining the wind force, soil tension or crop bending stress on the crop based on at least one of the meteorological data, the soil data and the crop data, comprising:

[0013] determining an air density based on the air pressure and the temperature, and determining a blade cross-sectional moment of inertia based on the blade thickness and the blade width;

[0014] determining a wind pressure based on the air density and the wind speed, and determining a deformed blade area based on the wind speed, the blade length, the elastic modulus, and the moment of inertia of the blade section;

[0015] determining wind friction based on the deformed blade area and the friction coefficient;

[0016] The wind force on the crop is determined based on the wind pressure, the blade area, the deformed blade area and the wind friction.

[0017] Optionally, the crop data include root dry matter weight, root dry matter density, taproot percentage, taproot radius, lateral root percentage, lateral root radius, fibrous root percentage, fibrous root radius, crop deadweight, effective root depth and root type aspect ratio; the soil data include soil type, soil type ratio, soil compaction, soil compaction ratio, soil moisture, soil water holding capacity, wilting point, soil bulk density and organic matter content;

[0018] Determining the wind force, soil tension or crop bending stress on the crop based on at least one of the meteorological data, the soil data and the crop data, comprising:

[0019] Determine the total root area based on the root dry matter weight, the root dry matter density, the taproot percentage, the taproot radius, the lateral root percentage, the lateral root radius, the fibrous root percentage, and the fibrous root radius;

[0020] determining an internal friction angle based on the soil type ratio, the wilting point, the soil moisture, and the soil water holding capacity;

[0021] Determine the total normal stress based on the crop deadweight, the soil bulk density, the effective root depth, the root type aspect ratio, the crop wind force and the horizontal angle of the terrain slope;

[0022] determining soil cohesion based on the soil type, the soil type ratio, the soil compaction, the soil compaction ratio, the wilting point, the soil moisture, the soil water holding capacity, and the organic matter content;

[0023] determining shear strength based on the soil cohesion, the total normal stress, and the internal friction angle;

[0024] The soil tension is determined based on the shear strength and the total root area.

[0025] Optionally, determining the total root area based on the root dry matter weight, the root dry matter density, the taproot percentage, the taproot radius, the lateral root percentage, the lateral root radius, the fibrous root percentage and the fibrous root radius comprises:

[0026] Determining the total root volume based on the root dry matter weight and the root dry matter density;

[0027] Determining the total area of ​​taproots based on the total root volume, the percentage of taproots, and the radius of taproots, determining the total area of ​​lateral roots based on the total root volume, the percentage of lateral roots, and the radius of lateral roots, and determining the total area of ​​fibrous roots based on the total root volume, the percentage of fibrous roots, and the radius of fibrous roots;

[0028] The total root area is determined based on the total main root area, the total lateral root area and the total fibrous root area.

[0029] Optionally, the total normal stress is determined based on the dead weight of the crop, the bulk density of the soil, the effective root depth, the aspect ratio of the root type, the wind force on the crop and the horizontal angle of the terrain slope, including:

[0030] The crop deadweight normal stress is determined based on the crop deadweight, the effective root depth and the root aspect ratio, the soil normal stress is determined based on the soil bulk density and the effective root depth, and the wind normal stress is determined based on the wind force on the crop and the horizontal angle of the terrain slope;

[0031] The total positive stress is determined based on the crop self-weight positive stress, the soil positive stress and the wind positive stress.

[0032] Optionally, determining soil cohesion based on the soil type, the soil type ratio, the soil compaction, the soil compaction ratio, the wilting point, the soil moisture, the soil water holding capacity and the organic matter content comprises:

[0033] determining foundation cohesion based on the soil type and the soil type ratio;

[0034] Determining a soil compactness correction coefficient based on the soil compactness and the soil compactness ratio;

[0035] determining a cohesion increment coefficient based on the wilting point, the soil moisture, the soil water holding capacity, the soil type ratio, and the soil compaction correction factor;

[0036] The soil cohesion is determined based on the basic cohesion, the cohesion increment coefficient and the organic matter content.

[0037] Optionally, the crop data at least includes the crop inner radius, the crop outer radius, the crop stem length, the elastic modulus and the distance from the crop center;

[0038] Determining the wind force, soil tension or crop bending stress on the crop based on at least one of the meteorological data, the soil data and the crop data, comprising:

[0039] Determining a stalk cross-sectional moment of inertia based on the crop inner radius and the crop outer radius;

[0040] Determining the crop bending moment based on the wind force on the crop and the length of the crop stem;

[0041] The crop bending stress is determined based on the stem cross-sectional moment of inertia, the crop bending moment and the distance from the crop center.

[0042] Optionally, determining whether the crop has fallen and the type of lodging based on the comparison result of the wind force on the crop, the soil tension, and the crop bending stress includes:

[0043] Based on the fact that the wind force on the crop is less than the crop bending stress and the soil tension, it is determined that the crop has not fallen over.

[0044] or,

[0045] Based on the fact that the wind force on the crop is greater than the crop bending stress and the soil tension, it is determined that the lodging type is the breakage type lodging;

[0046] or,

[0047] Based on the fact that the wind force on the crop is less than the soil tension and greater than the crop bending stress, it is determined that the lodging type is the breakage type lodging;

[0048] or,

[0049] Based on the fact that the wind force on the crop is greater than the soil tension and less than the crop bending stress, it is determined that the lodging type is non-breakage lodging.

[0050] In a second aspect, the present invention further provides a device for determining the type of lodging of crops caused by wind disasters, comprising:

[0051] A data acquisition module, used to acquire meteorological data, soil data and crop data;

[0052] A first determination module, configured to determine the wind force, soil tension or crop bending stress of the crop based on at least one of the meteorological data, the soil data and the crop data;

[0053] A second determination module is used to determine whether the crop has fallen and the type of lodging based on the comparison result of the wind force on the crop, the soil tension, and the crop bending stress;

[0054] The meteorological data are data representing the influence of environmental factors on crop lodging, the soil data are data representing the influence of soil factors on crop lodging, and the crop data are data representing the influence of crops themselves on crop lodging;

[0055] The lodging types include breaking type lodging and non-breaking type lodging.

[0056] In a third aspect, the present invention further provides a computer-readable storage medium, wherein the computer-readable storage medium comprises computer program instructions, wherein the computer program instructions enable a computer to execute the steps of the method as described in one of the first aspects.

[0057] The embodiments of the present invention have the following technical effects:

[0058] The method for determining the type of crop lodging caused by wind disasters provided by the present invention determines the wind force, soil tension or crop bending stress of crops by acquiring meteorological data, soil data and crop data before or after the occurrence of the wind disaster. Among them, the wind force on crops represents the force exerted by the wind disaster on crops, the soil tension represents the force required to blow down crops, and the crop bending stress represents the force required to break the stems of crops. During the occurrence of wind disasters, meteorological data, soil data and crop data will affect the wind force, soil tension and crop bending stress of crops. Therefore, the embodiment of the present invention can accurately determine the wind force, soil tension and crop bending stress of crops through the above scheme. After determining the wind force, soil tension and crop bending stress of crops, by comparing the wind force on crops with soil tension and crop bending stress, it is possible to accurately determine whether the crops have lodged and the type of lodging, such as non-lodging, broken type lodging and non-broken type lodging. After obtaining the above results, appropriate rescue measures can be made in a targeted manner to reduce the impact of wind disasters on crops. In summary, the embodiments of the present invention can accurately determine whether crops have fallen and the type of lodging based on the above scheme, which is conducive to taking targeted measures to reduce the impact of wind disasters on crops or formulating measures to save crop production. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0060] Figure 1 A schematic flow chart of a method for determining the type of lodging of crops caused by wind disasters provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0061] In order to make the purpose, technical solution and advantages of the present invention clearer, the technical solution of the present invention will be described clearly and completely below. Obviously, the described embodiments are only 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 belong to the scope of protection of the present invention.

[0062] Figure 1 A schematic flow chart of a method for determining the type of lodging of crops caused by wind disasters provided in an embodiment of the present invention, the method comprising:

[0063] S101. Acquire meteorological data, soil data, and crop data.

[0064] Specifically, the "crop" in the embodiment of the present invention refers to agricultural crops, such as corn, rice, and wheat, etc. Meteorological data is data about the environment in which the crop is located, such as wind speed, etc. Soil data is data about the soil in which the crop is located, such as soil type (sand, loam, clay), etc. Crop data is data about the crop itself, such as the length and diameter of the crop, etc.

[0065] During a wind disaster, meteorological data, soil data, and crop data will all have an impact on the type of lodging. For example, if the wind speed is high enough, the crops will fall over; if the soil is loose enough, the soil will not be able to fix the crops well, and the crops will fall over easily; if the diameter of the crops is small, the crops will have poor resistance to breaking, and they will break easily.

[0066] S102. Determine whether the crop is subjected to wind force, soil tension or crop bending stress based on at least one of meteorological data, soil data and crop data.

[0067] In the process of determining the lodging type, the embodiment of the present invention takes meteorological data, soil data and crop data into consideration. Through the above-mentioned multiple types of multiple data, the magnitude of the wind force, soil tension and crop bending stress on the crop can be accurately determined, and then the lodging type is determined based on the comparison results of the wind force on the crop with the soil tension and crop bending stress.

[0068] Among them, the wind force on crops refers to the force exerted by wind disasters on crops; soil tension refers to the force that the soil can use to fix crops to keep them upright, which can also be understood as the rooting force required to blow crops down; crop bending stress refers to the force required to break (blow off) the stems of crops.

[0069] S103: Determine whether the crop has fallen and the type of lodging based on the comparison results of the wind force on the crop, the soil tension, and the bending stress of the crop.

[0070] Among them, meteorological data are data representing the impact of environmental factors on crop lodging, soil data are data representing the impact of soil factors on crop lodging, and crop data are data representing the impact of crops themselves on crop lodging.

[0071] The lodging types include breaking lodging and non-breaking lodging.

[0072] In some embodiments, S103 includes:

[0073] The crop has not fallen over, based on the fact that the wind force on the crop is less than the crop bending stress and soil tension.

[0074] or,

[0075] The lodging type is determined to be the breaking type based on the fact that the wind force on the crop is greater than the crop bending stress and soil tension.

[0076] or,

[0077] The lodging type is determined to be the breaking type based on the fact that the wind force on the crop is less than the soil tension but greater than the crop bending stress.

[0078] or,

[0079] The lodging type is determined to be non-breakage type based on the fact that the wind force on the crop is greater than the soil tension but less than the crop bending stress.

[0080] Based on the above embodiments, it can be understood that if the wind force on the crops is greater than the soil tension and less than the crop bending stress, it means that the current wind disaster may blow down the crops, but cannot break the crops, so the lodging type is non-break lodging. If the wind force on the crops is greater than the crop bending stress and greater than the soil tension, it means that the wind disaster can both blow down the crops and break the crops, so the lodging type is break lodging. If the wind force on the crops is less than the crop bending stress but greater than the soil tension, it means that the wind disaster cannot blow down the crops, but can break the crops, so the lodging type is break lodging. If the wind force on the crops is less than the crop bending stress and less than the soil tension, it means that the wind disaster can neither blow down the crops nor break the crops, so the crops have not lodged.

[0081] In summary, through the above scheme, the embodiment of the present invention can accurately determine the magnitude of the wind force, soil tension and crop bending stress of the crop, and determine whether the crop has lodged or not based on the comparison results of the wind force, soil tension and crop bending stress of the crop. Furthermore, based on the lodging type, targeted rescue measures can be taken to reduce the impact of wind disasters on crops, or measures to rescue crop production can be formulated directly based on the prediction results after the wind disaster occurs.

[0082] In some embodiments, the meteorological data includes at least air pressure, temperature and wind speed, and the crop data includes at least leaf thickness, leaf width, leaf length, leaf area, elastic modulus and friction coefficient.

[0083] S102 includes:

[0084] The air density is determined based on the air pressure and temperature, and the blade section moment of inertia is determined based on the blade thickness and blade width.

[0085] The wind pressure is determined based on the air density and the wind speed, and the deformed blade area is determined based on the wind speed, the blade length, the elastic modulus and the moment of inertia of the blade section.

[0086] Wind friction is determined based on the deformed blade area and the friction coefficient.

[0087] The wind force on the crop is determined based on wind pressure, leaf area, deformed leaf area and wind friction.

[0088] Specifically, embodiments of the present invention can determine crop data such as leaf thickness, leaf width, leaf length, leaf area, elastic modulus and friction coefficient based on the crop model. The crop model represents an algorithm model used to determine the above-mentioned crop data, which can be obtained by training after a large number of sampling samples are input. Air pressure, temperature and wind speed can be measured, predicted or obtained through weather forecasts through corresponding tools. The elastic modulus represents the positive proportionality coefficient of the positive proportional relationship between the stress and strain of the material during the elastic deformation stage of the material. The friction coefficient represents the ratio of the friction between two surfaces to the vertical force acting on one of the surfaces. The moment of inertia of the blade section is a physical data that measures the bending resistance of the blade section. The deformed blade area represents the area of ​​the blade that is deformed under the action of wind disasters. Wind friction represents the force exerted on the crop by the wind disaster in the process of blowing over the side of the crop.

[0089] In the embodiment of the present invention, the air density can be determined by the following formula (1):

[0090] ; (1)

[0091] in, Indicates air density. The unit of measurement for air density is kg / m 3 , P represents air pressure, and the unit of air pressure is Pascal; T represents temperature, and the unit of temperature is Kelvin.

[0092] The blade section moment of inertia is determined by the following formula (2):

[0093] ; (2)

[0094] Among them, F1 represents the moment of inertia of the blade section, h1 represents the blade thickness, and w1 represents the blade width.

[0095] After calculating the moment of inertia of the blade section, the wind pressure is determined by the following formula (3):

[0096] ; (3)

[0097] Among them, F2 represents wind pressure and v1 represents wind speed.

[0098] The deformed blade area is determined by the following formula (4):

[0099] ; (4)

[0100] Among them, S1 represents the deformed blade area, l1 represents the blade length, and E1 represents the elastic modulus.

[0101] After determining the deformed blade area, the wind friction force is determined by the following formula (5):

[0102] ; (5)

[0103] Among them, F3 represents wind friction, Represents the friction coefficient.

[0104] After calculating the wind friction, the wind force on the crop is determined by the following formula (6):

[0105] ; (6)

[0106] Among them, F4 represents the wind force on the crop, and S2 represents the leaf area.

[0107] Based on the above scheme, the wind force on the crops can be determined. It should be noted that in the embodiment of the present invention, multiple formulas are involved, and letters are used to refer to noun definitions. For all letters in all formulas involved, they are case-sensitive. For example, the definitions of V and v are not the same.

[0108] In the above embodiment, the calculated wind force on the crop is the wind force on each leaf (when facing the wind direction). For a crop, it generally includes multiple leaves, and the orientation of each leaf is also different. Under different wind directions, the wind generates wind force on each leaf of the crop. For this, the embodiment of the present invention can obtain the number of leaves and the orientation of each leaf from the crop model, randomly generate the number of leaves and the orientation of each leaf, and determine the total wind force on all leaves of each column of crops based on the number of leaves and the orientation of the leaves. The embodiment of the present invention provides a simple calculation method. For example, when the leaves are facing the wind direction, the wind force on the crop is , if the acute angle between the plane where the blade is located (the blade is simply considered to be a plane) and the wind direction is , we can determine that the wind force on the crop corresponding to the leaf is ;in, Indicates that the crop corresponding to this leaf is subject to wind force.

[0109] The total wind force on all leaves is It is the sum of the wind forces acting on each leaf, and will not be shown here through a formula.

[0110] In some embodiments, crop data include root dry matter weight, root dry matter density, taproot percentage, taproot radius, lateral root percentage, lateral root radius, fibrous root percentage, fibrous root radius, crop deadweight, effective root depth and root type aspect ratio, and soil data include soil type, soil type ratio, soil compaction, soil compaction ratio, soil moisture, soil water holding capacity, wilting point, soil bulk density and organic matter content.

[0111] S102 includes:

[0112] The total root area was determined based on root dry matter weight, root dry matter density, taproot percentage, taproot radius, lateral root percentage, lateral root radius, fibrous root percentage, and fibrous root radius.

[0113] The internal friction angle is determined based on the proportion of soil types, wilting point, soil moisture, and soil water holding capacity.

[0114] The total normal stress is determined based on crop deadweight, soil bulk density, effective root depth, root aspect ratio, crop wind force and horizontal angle of terrain slope.

[0115] Soil cohesion is determined based on soil type, soil type ratio, soil compaction, soil compaction ratio, wilting point, soil moisture, soil water holding capacity, and organic matter content.

[0116] Determine the shear strength based on soil cohesion, total normal stress, and internal friction angle.

[0117] Determine soil tension based on shear strength and total root area.

[0118] Among them, the root dry matter weight indicates the weight of the material without water in the root of the crop. The root dry matter density indicates the density of the material without water in the root of the crop. The taproot percentage indicates the proportion of the taproot volume of the crop, and the taproot radius indicates the radius of the taproot of the crop. The lateral root percentage, lateral root radius, fibrous root percentage, fibrous root radius and taproot percentage and taproot radius are the same. The crop deadweight indicates the weight of the crop. The wilting point indicates the lowest soil water content that cannot be absorbed and utilized by the crop. The effective root depth indicates the root depth that the crop can effectively absorb water and nutrients, which can be obtained through the crop model. The root type aspect ratio indicates the ratio of the crop root radius to the root depth, which can be obtained through the crop model. The horizontal angle of the terrain slope indicates the angle between the terrain where the crop is located and the horizontal direction. Soil types can include sandy soil, loam and clay. The soil type ratio indicates the proportion of each soil type. The soil compaction indicates the looseness of the soil. The soil compaction ratio indicates the proportion of soils with different soil compactions. Soil moisture indicates the water content in the soil. Soil water holding capacity refers to the amount of water that the soil can stably maintain after sufficient rainfall or irrigation. Soil bulk density refers to the weight of one cubic meter of soil. Organic matter content refers to the proportion of organic matter in the soil. Shear strength refers to the shear strength when a substance yields or fails under shear force. Soil cohesion, also known as soil adhesion, refers to the binding force between soil particles and represents a measure of the strength of the association between soil particles. Normal stress refers to the internal force generated between the parts of an object when it is deformed due to external factors (force) to resist the effect of such external factors and try to restore the object from the deformed position to the position before deformation. The total normal stress represents the collection of normal stresses to which the crop is subjected. The internal friction angle is one of the shear strength indicators of soil or rock, reflecting the magnitude of the internal friction between the particles inside the soil or rock. The larger the internal friction angle, the higher the strength. Soil data can be obtained or calculated through public information, such as the "Chinese Land Surface Simulation Soil Property Dataset".

[0119] In some embodiments, determining the total root area based on root dry matter weight, root dry matter density, taproot percentage, taproot radius, lateral root percentage, lateral root radius, fibrous root percentage, and fibrous root radius comprises:

[0120] The total root volume was determined based on the root dry matter weight and root dry matter density.

[0121] The total area of ​​taproots was determined based on the total root volume, the percentage of taproots, and the radius of taproots; the total area of ​​lateral roots was determined based on the total root volume, the percentage of lateral roots, and the radius of lateral roots; and the total area of ​​fibrous roots was determined based on the total root volume, the percentage of fibrous roots, and the radius of fibrous roots.

[0122] The total root area was determined based on the total taproot area, total lateral root area, and total fibrous root area.

[0123] Specifically, the total root volume can be determined according to the following formula (7):

[0124] ; (7)

[0125] Where V1 represents the total volume of the root, W1 represents the dry matter weight of the root, Indicates the root dry matter density.

[0126] After determining the total root volume, the total taproot area is determined using the following formulas (8), (9), and (10):

[0127] ; (8)

[0128] ; (9)

[0129] ; (10)

[0130] Among them, V2 represents the total volume of the taproot, m1 represents the percentage of the taproot, l2 represents the total length of the taproot, r1 represents the radius of the taproot, and S2 represents the total area of ​​the taproot.

[0131] Similar to (8), (9) and (10), the total lateral root area S3 and the total fibrous root area S4 can be determined respectively, which will not be elaborated here.

[0132] After determining the total area of ​​the taproot, the total area of ​​the lateral roots and the total area of ​​the fibrous roots, the total root area is determined according to the following formula (11):

[0133] ; (11)

[0134] Where S5 represents the total root area.

[0135] In summary, the total root area can be determined based on the above scheme.

[0136] In some embodiments, determining the total normal stress based on the crop's own weight, soil bulk density, effective root depth, root aspect ratio, and the wind force on the crop comprises:

[0137] The crop deadweight positive stress is determined based on the crop deadweight, effective root depth and aspect ratio of the root type; the soil positive stress is determined based on the soil bulk density and effective root depth; and the wind positive stress is determined based on the wind force on the crop and the horizontal angle of the terrain slope.

[0138] The total normal stress is determined based on the crop self-weight normal stress, soil normal stress and wind normal stress.

[0139] Specifically, the crop self-weight positive stress can be determined by the following formula (12):

[0140] ; (12)

[0141] Among them, F5 represents the positive stress of crop deadweight, W2 represents the deadweight of crop, and S6 represents the root radius area.

[0142] The root radius area is the area calculated from the root radius, which can be obtained by the circle calculation formula and will not be described in detail. The root radius can be determined by the following formula (13):

[0143] ; (13)

[0144] Among them, r2 represents the root radius, h2 represents the effective root depth, and n represents the root aspect ratio.

[0145] The soil normal stress is determined by the following formula (14):

[0146] ; (14)

[0147] Among them, F6 represents soil normal stress and d represents soil bulk density.

[0148] The wind normal stress is determined by the following formula (15):

[0149] ; (15)

[0150] Among them, F7 represents the wind positive stress, Represents the horizontal angle of terrain slope.

[0151] After obtaining the crop self-weight normal stress, soil cohesion and wind normal stress, the total normal stress can be determined as:

[0152] ; (16)

[0153] Where F8 represents the total normal stress.

[0154] In some embodiments, determining soil cohesion based on soil type, soil type ratio, soil compaction, soil compaction ratio, wilting point, soil moisture, soil water holding capacity, and organic matter content comprises:

[0155] Determine the foundation cohesion based on soil type and soil type ratio.

[0156] The soil compaction correction factor is determined based on the soil compaction and the soil compaction ratio.

[0157] The incremental cohesion factor is determined based on the wilting point, soil moisture, soil water holding capacity, proportion of soil types and a correction factor for soil compaction.

[0158] Soil cohesion was determined based on basic cohesion, cohesion increment coefficient and organic matter content.

[0159] Specifically, the basic cohesion can be determined by the following formula (17):

[0160] ; (17)

[0161] Among them, F9 represents the basic cohesion, P1 represents the cohesion value of the loam silt in the first (highest) clod form, P2 represents the cohesion value of the loam silt in the second clod form, P3 represents the cohesion value of the clay in the first clod form, and P4 represents the cohesion value of the clay in the second clod form. Q1 represents the proportion of loam silt in the first clod form, Q2 represents the proportion of loam silt in the second clod form, Q3 represents the proportion of clay in the first clod form, and Q4 represents the proportion of clay in the second clod form.

[0162] SW1 represents the proportion of the first (highest) soil block shape, and SW2 represents the proportion of the second soil block shape. Soil block shapes can include many types, such as angular soil, angular block soil, block soil, large block soil, granular soil, columnar soil, etc., as shown in Table 1. Soil types are divided into sandy soil, loam and clay. When the soil block shape is the same, the assigned value of clay is the highest, the assigned value of loam is the second, and the assigned value of sand is 0, so it is not included in formula (16). On the basis of being loam or clay, the assigned values ​​of different soil block shapes may be the same or different. The cohesion correction data can be found based on the cohesion correction data table. Different soil block shapes correspond to different cohesion correction data.

[0163] Table 1 Cohesion value table of soil block shape

[0164]

[0165] As shown in Table 1, the column corresponding to sand is the cohesion assignment values ​​of different soil block shapes corresponding to sand, the column corresponding to loam is the cohesion assignment values ​​of different soil block shapes corresponding to loam, and the column corresponding to clay is the cohesion assignment values ​​of different soil block shapes corresponding to loam.

[0166] The soil compaction correction data is determined by the following formula (18):

[0167] ; (18)

[0168] Among them, df represents the soil compactness correction data, cv1 represents the soil compactness assignment according to the soil compactness ratio of the first, cv2 represents the soil compactness assignment according to the soil compactness ratio of the second, cw1 represents the proportion of the soil type according to the soil compactness ratio of the first, and cw2 represents the proportion of the soil type according to the soil compactness ratio of the second. The compactness assignment and soil type ratio can be found through the corresponding table. Different soil types have different compactness assignment values, different soil types have different soil type ratios, and different compactness.

[0169] The following is an exemplary soil compaction value table.

[0170] Table 2 Soil compactness assignment table

[0171]

[0172] Among them, c1 and c2 represent the soil compactness data of the top two soils.

[0173] The cohesion increment coefficient is determined by the following formula (19):

[0174] ; (19)

[0175] Among them, K represents the cohesion increment coefficient, silt represents the loam ratio, sand represents the sand ratio, clay represents the clay ratio, sm represents soil moisture, smw represents the wilting point, smfcf represents soil water holding capacity, and Abs represents the absolute value. Table 2 is a table of the cohesion increment coefficient range. If the calculated cohesion increment coefficient is within the limit value interval, the calculated value is taken; if it is greater than the limit value interval, the maximum value is taken; if it is less than the limit value interval, the minimum value is taken.

[0176] Table 3 Cohesion increment coefficient range

[0177]

[0178] After calculating the cohesion increment coefficient and the foundation cohesion, the soil cohesion is determined by the following formula (20):

[0179] ; (20)

[0180] Among them, F 10 represents soil cohesion, and q represents organic matter content.

[0181] In some embodiments, determining the internal friction angle based on soil type ratio, wilting point, soil moisture, and soil water holding capacity comprises:

[0182] Specifically, the internal friction angle is determined by the following equations (21)-(24):

[0183] (twenty one)

[0184] (twenty two)

[0185] (twenty three)

[0186] ; (twenty four)

[0187] in, represents the internal friction angle, is the internal friction angle of sand, is the internal friction angle of loam, Represents the internal friction angle of clay. After calculating the internal friction angle, the calculated result can be compared with the internal friction angle range obtained by multiplying the friction angle by the ratio of sand, loam and clay in the area. If the calculated internal friction angle is within the range, the calculated result is directly taken. If it is greater than the maximum value of the internal friction angle range, the maximum value of the internal friction angle range is taken. If it is less than the minimum value of the internal friction angle range, the minimum value of the internal friction angle range is taken. Table 4 is a table of the internal friction angle range.

[0188] Table 4 Internal friction angle range

[0189]

[0190] As shown in Table 4, the internal friction angle of loam is in the range of 15-25, the internal friction angle of sand is in the range of 25-35, and the internal friction angle of clay is in the range of 10-30. In some embodiments, the internal friction angle may also be calculated in other ways, which are not limited in the embodiments of the present invention.

[0191] After determining the soil cohesion, total normal stress, and internal friction angle, the shear strength can be determined as:

[0192] ; (25)

[0193] Where t represents the shear strength.

[0194] The soil tension is further determined based on the shear strength and the total root area:

[0195] ; (25)

[0196] Among them, F 11 Represents soil tension.

[0197] In some embodiments, the crop data includes at least the crop inner radius, the crop outer radius, the crop stem length, the elastic modulus, and the distance from the center of the crop.

[0198] S102 includes:

[0199] The stalk cross-sectional moment of inertia is determined based on the crop inner radius and the crop outer radius.

[0200] The crop bending moment is determined based on the wind force on the crop and the length of the crop stem.

[0201] The crop bending stress is determined based on the stalk cross-sectional moment of inertia, crop bending moment and distance from the crop center.

[0202] Specifically, the moment of inertia of the stem cross section can be determined by the following formula (26):

[0203] ; (26)

[0204] Where I represents the moment of inertia of the stalk cross section, r3 represents the inner radius of the crop, and r4 represents the outer radius of the crop.

[0205] The crop bending moment is determined by the following formula (27):

[0206] ; (27)

[0207] Where M represents the bending moment, H represents the stem length, and F 12 It represents the maximum bending stress of the crop and can be obtained by measurement.

[0208] The crop bending stress is determined by the following formula (28):

[0209] ; (28)

[0210] in, Indicates the crop bending stress, y is the distance from the center of the crop. Specifically, the distance y from the center of the crop can be simply understood as half of the length of the crop stem. For example, if the length of the crop stem is 1, then the distance from the center of the crop is 0.5; the crop inner radius and crop outer radius represent the inner radius and outer radius of the crop stem when the crop stem is hollow.

[0211] It has been explained in the above embodiments that determining the type of crop lodging only requires comparing the wind force on the crop with the crop bending stress and soil tension, and the above data have been obtained through calculation, so the type of lodging can be accurately determined.

[0212] Specifically, if and , the crop has not fallen. but , then the lodging type is broken lodging, if and , then the lodging type is broken lodging, if but , the lodging type is non-fracture lodging.

[0213] The embodiment of the present invention further provides a device for determining the type of lodging of crops caused by wind disasters, the device comprising:

[0214] The data acquisition module is used to acquire meteorological data, soil data and crop data.

[0215] The first determination module is used to determine whether the crop is subjected to wind force, soil tension or crop bending stress based on at least one of meteorological data, soil data and crop data.

[0216] The second determination module is used to determine whether the crop has fallen and the type of lodging based on the comparison results of the wind force on the crop, the soil tension, and the bending stress of the crop.

[0217] Among them, meteorological data are data representing the impact of environmental factors on crop lodging, soil data are data representing the impact of soil factors on crop lodging, and crop data are data representing the impact of crops themselves on crop lodging.

[0218] The lodging types include breaking lodging and non-breaking lodging.

[0219] A device for determining the type of lodging of crops caused by wind disasters provided in an embodiment of the present invention has the same technical features as the above method, and can therefore solve the same technical problems and achieve the same technical effects, which will not be elaborated herein.

[0220] An embodiment of the present invention further provides a computer-readable storage medium, which includes computer program instructions, and the computer program instructions enable a computer to execute the steps of the method provided by any embodiment of the present application.

[0221] The computer program instructions may be written in any combination of one or more programming languages ​​to form program codes for performing the operations of the embodiments of the present application, including object-oriented programming languages ​​such as Java, C++, etc., and conventional procedural programming languages ​​such as "C" or similar programming languages. The program code may be executed entirely on the user computing device, partially on the user computing device, as a separate software package, partially on the user computing device and partially on a remote computing device, or entirely on a remote computing device or server.

[0222] Computer readable storage media can adopt any combination of one or more readable media. The readable medium can be a readable signal medium or a readable storage medium. The readable storage medium can include, for example, but is not limited to, a system, device or device of electricity, magnetism, light, electromagnetic, infrared, or semiconductor, or any combination of the above. More specific examples (non-exhaustive list) of readable storage media include: an electrical connection with one or more wires, a portable disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the above.

[0223] 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 replace some or all of the technical features therein by equivalents. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the technical solutions of the embodiments of the present invention.

Claims

1. A method for determining the type of lodging of crops caused by wind disasters, characterized in that: include: S101, obtaining meteorological data, soil data and crop data; S102, determining whether the crop is subjected to wind force, soil tension or crop bending stress based on at least one of the meteorological data, the soil data and the crop data; S103, determining whether the crop has fallen and the type of lodging based on the comparison result of the wind force on the crop, the soil tension, and the crop bending stress; The meteorological data are data representing the influence of environmental factors on crop lodging, the soil data are data representing the influence of soil factors on crop lodging, and the crop data are data representing the influence of crops themselves on crop lodging; The lodging types include broken lodging and non-broken lodging; The crop data include root dry matter weight, root dry matter density, taproot percentage, taproot radius, lateral root percentage, lateral root radius, fibrous root percentage, fibrous root radius, crop deadweight, effective root depth and root type aspect ratio; the soil data include soil type, soil type ratio, soil compaction, soil compaction ratio, soil moisture, soil water holding capacity, wilting point, soil bulk density and organic matter content; Determining the wind force, soil tension or crop bending stress on the crop based on at least one of the meteorological data, the soil data and the crop data, comprising: Determine the total root area based on the root dry matter weight, the root dry matter density, the taproot percentage, the taproot radius, the lateral root percentage, the lateral root radius, the fibrous root percentage, and the fibrous root radius; determining an internal friction angle based on the soil type ratio, the wilting point, the soil moisture, and the soil water holding capacity; Determine the total normal stress based on the crop's deadweight, the soil bulk density, the effective root depth, the root type aspect ratio, the crop's wind force and the horizontal angle of the terrain slope; determining soil cohesion based on the soil type, the soil type ratio, the soil compaction, the soil compaction ratio, the wilting point, the soil moisture, the soil water holding capacity, and the organic matter content; determining shear strength based on the soil cohesion, the total normal stress, and the internal friction angle; The soil tension is determined based on the shear strength and the total root area.

2. The method according to claim 1, characterized in that The meteorological data include air pressure, temperature and wind speed, and the crop data include leaf thickness, leaf width, leaf length, leaf area, elastic modulus and friction coefficient; Determining the wind force, soil tension or crop bending stress on the crop based on at least one of the meteorological data, the soil data and the crop data, comprising: determining an air density based on the air pressure and the temperature, and determining a blade cross-sectional moment of inertia based on the blade thickness and the blade width; determining a wind pressure based on the air density and the wind speed, and determining a deformed blade area based on the wind speed, the blade length, the elastic modulus, and the moment of inertia of the blade section; determining wind friction based on the deformed blade area and the friction coefficient; The wind force on the crop is determined based on the wind pressure, the blade area, the deformed blade area and the wind friction.

3. The method according to claim 1, characterized in that Determining the total root area based on the root system dry matter weight, the root system dry matter density, the taproot percentage, the taproot radius, the lateral root percentage, the lateral root radius, the fibrous root percentage and the fibrous root radius comprises: Determining the total root volume based on the root dry matter weight and the root dry matter density; Determine the total area of ​​the taproot based on the total root volume, the taproot percentage, and the taproot radius; determine the total area of ​​the lateral roots based on the total root volume, the lateral root percentage, and the lateral root radius; and determine the total area of ​​the fibrous roots based on the total root volume, the fibrous root percentage, and the fibrous root radius; The total root area is determined based on the total main root area, the total lateral root area and the total fibrous root area.

4. The method according to claim 1, characterized in that: The total normal stress is determined based on the dead weight of the crop, the bulk density of the soil, the effective root depth, the aspect ratio of the root type, the wind force on the crop and the horizontal angle of the terrain slope, including: Determine the crop deadweight normal stress based on the crop deadweight, the effective root depth and the root type aspect ratio, determine the soil normal stress based on the soil bulk density and the effective root depth, and determine the wind normal stress based on the wind force on the crop and the horizontal angle of the terrain slope; The total positive stress is determined based on the crop self-weight positive stress, the soil positive stress and the wind positive stress.

5. The method according to claim 1, characterized in that Determining soil cohesion based on the soil type, the soil type ratio, the soil compaction, the soil compaction ratio, the wilting point, the soil moisture, the soil water holding capacity, and the organic matter content, comprises: determining foundation cohesion based on the soil type and the soil type ratio; Determining a soil compactness correction coefficient based on the soil compactness and the soil compactness ratio; determining a cohesion increment coefficient based on the wilting point, the soil moisture, the soil water holding capacity, the soil type ratio, and the soil compaction correction factor; The soil cohesion is determined based on the basic cohesion, the cohesion increment coefficient and the organic matter content.

6. The method according to claim 1, characterized in that The crop data at least includes the crop inner radius, the crop outer radius, the crop stem length, the elastic modulus and the distance from the crop center; Determining the wind force, soil tension or crop bending stress on the crop based on at least one of the meteorological data, the soil data and the crop data, comprising: Determining a stalk cross-sectional moment of inertia based on the crop inner radius and the crop outer radius; Determining the crop bending moment based on the wind force on the crop and the length of the crop stem; The crop bending stress is determined based on the stem cross-sectional moment of inertia, the crop bending moment and the distance from the crop center.

7. The method according to claim 1, characterized in that Based on the comparison result of the wind force on the crop, the soil tension and the crop bending stress, determining whether the crop has fallen and the type of the falling, including: Based on the fact that the wind force on the crop is less than the crop bending stress and the soil tension, it is determined that the crop has not fallen over. or, Based on the fact that the wind force on the crop is greater than the crop bending stress and the soil tension, it is determined that the lodging type is the breakage type lodging; or, Based on the fact that the wind force on the crop is less than the soil tension and greater than the crop bending stress, it is determined that the lodging type is the breakage type lodging; or, Based on the fact that the wind force on the crop is greater than the soil tension and less than the crop bending stress, it is determined that the lodging type is non-breakage lodging.

8. A device for determining the type of lodging of crops caused by windstorms, characterized in that: include: A data acquisition module, used to acquire meteorological data, soil data and crop data; A first determination module, configured to determine the wind force, soil tension or crop bending stress of the crop based on at least one of the meteorological data, the soil data and the crop data; A second determination module is used to determine whether the crop has fallen and the type of lodging based on the comparison result of the wind force on the crop, the soil tension, and the crop bending stress; The meteorological data are data representing the influence of environmental factors on crop lodging, the soil data are data representing the influence of soil factors on crop lodging, and the crop data are data representing the influence of crops themselves on crop lodging; The lodging types include broken lodging and non-broken lodging; The crop data include root dry matter weight, root dry matter density, taproot percentage, taproot radius, lateral root percentage, lateral root radius, fibrous root percentage, fibrous root radius, crop deadweight, effective root depth and root type aspect ratio; the soil data include soil type, soil type ratio, soil compaction, soil compaction ratio, soil moisture, soil water holding capacity, wilting point, soil bulk density and organic matter content; Determining the wind force, soil tension or crop bending stress on the crop based on at least one of the meteorological data, the soil data and the crop data, comprising: Determine the total root area based on the root dry matter weight, the root dry matter density, the taproot percentage, the taproot radius, the lateral root percentage, the lateral root radius, the fibrous root percentage, and the fibrous root radius; determining an internal friction angle based on the soil type ratio, the wilting point, the soil moisture, and the soil water holding capacity; Determine the total normal stress based on the crop's deadweight, the soil bulk density, the effective root depth, the root type aspect ratio, the crop's wind force and the horizontal angle of the terrain slope; determining soil cohesion based on the soil type, the soil type ratio, the soil compaction, the soil compaction ratio, the wilting point, the soil moisture, the soil water holding capacity, and the organic matter content; determining shear strength based on the soil cohesion, the total normal stress, and the internal friction angle; The soil tension is determined based on the shear strength and the total root area.

9. A computer-readable storage medium, characterized in that: The computer-readable storage medium comprises computer program instructions, which enable a computer to execute the steps of the method according to any one of claims 1 to 7.

Images