Soil erosion intensity measuring and calculating method and device based on high-resolution data driving

Abstract

The invention relates to the technical field of soil erosion monitoring and evaluation, and discloses a soil erosion intensity measuring and calculating method and device based on high-resolution data driving, and the method comprises the steps: dividing water and soil loss data according to an engineering disturbance type and a catchment type, and obtaining a soil loss type; unmanned aerial vehicle operation data and disturbance survey data corresponding to the target production construction project are acquired, and land utilization type data, water and soil conservation measure data, rainfall data and surface change data are determined based on the unmanned aerial vehicle operation data and the disturbance survey data; calculating water and soil loss influence factors based on the land utilization type data, the water and soil conservation measure data, the rainfall data and the earth surface change data; and respectively calculating the soil erosion intensity based on the water and soil loss influence factors according to the soil loss type. According to the method, the precision and spatial resolution of soil erosion intensity measurement and calculation are remarkably improved, and a scientific basis is provided for soil erosion control and ecological environment construction in a small-scale range.

Description

Technical Field

[0001] This invention relates to the field of soil erosion monitoring and assessment technology, specifically to a method and apparatus for calculating soil erosion intensity based on high-resolution data. Background Technology

[0002] Soil erosion has a profound impact on the stability of ecosystems, the sustainable use of land resources, and the healthy development of agricultural production. Accurate and efficient measurement of soil erosion intensity is a key prerequisite for scientifically preventing and controlling soil erosion, rationally planning land resources, and protecting the ecological environment.

[0003] However, existing soil erosion monitoring methods cannot reflect changes in soil erosion and are insufficient to meet the needs for real-time monitoring and accurate assessment. Summary of the Invention

[0004] This invention provides a method and apparatus for calculating soil erosion intensity based on high-resolution data, in order to solve the problem that related soil erosion monitoring methods cannot reflect changes in soil erosion and are difficult to meet the needs of real-time monitoring and accurate assessment.

[0005] In a first aspect, the present invention provides a method for calculating soil erosion intensity based on high-resolution data, the method comprising: Obtain soil and water loss data corresponding to the target production and construction project, and classify the soil and water loss data according to the engineering disturbance type and water catchment type to obtain the soil loss type; Obtain drone operation data and disturbance survey data corresponding to the target production and construction project, and determine land use type data, soil and water conservation measures data, rainfall data and surface change data based on the drone operation data and disturbance survey data; The soil erosion influencing factors are calculated based on land use type data, soil and water conservation measures data, rainfall data, and surface change data. Soil erosion intensity was calculated based on soil loss type and water and soil loss influencing factors.

[0006] This embodiment provides a method for calculating soil erosion intensity based on high-resolution data. It acquires UAV operation data and disturbance survey data corresponding to the target production and construction project, utilizes UAV remote sensing technology to predict and monitor soil erosion at the project scale, and then determines land use type data, soil and water conservation measures data, rainfall data, and surface change data based on the UAV operation data and disturbance survey data. Based on these data, it calculates soil erosion influencing factors. By fully utilizing the rich information in high-resolution data, it can more precisely depict the spatial variations of various soil erosion influencing factors, significantly improve the accuracy and spatial resolution of soil erosion intensity calculation, and more accurately reflect the actual situation of soil erosion, providing a scientific basis for serving small-scale soil erosion prevention and control and ecological environment construction.

[0007] In one optional implementation, land use type data, soil and water conservation measures data, rainfall data, and surface change data are determined based on UAV operation data and disturbance survey data, including: Based on UAV operational data, determine UAV imagery data, multispectral data, and laser point cloud data; The UAV image data, multispectral data, and laser point cloud data were preprocessed to obtain digital surface model data, multispectral image data, and digital elevation model data, respectively. Land use type data, soil and water conservation measures data, rainfall data, and surface change data were determined based on digital surface model data, multispectral image data, digital elevation model data, and disturbance survey data, respectively.

[0008] This embodiment provides a method for calculating soil erosion intensity based on high-resolution data. It uses UAV remote sensing technology to predict and monitor soil and water loss at the project scale, overcoming the limitations of insufficient accuracy, high difficulty, and high cost of monitoring methods for large-scale production and construction projects.

[0009] In one optional implementation, soil erosion data are classified according to engineering disturbance type and water catchment type to obtain soil erosion types, including: Based on soil erosion data, determine the current geomorphological data and the water catchment data above the target area; based on the current geomorphological data, determine the type of engineering disturbance; and based on the water catchment data above the target area, determine the type of water catchment. Soil loss types are classified based on engineering disturbance type and water catchment type. Soil loss types include disturbance without water catchment encroachment, disturbance without water catchment excavation, disturbance without water catchment deposition, disturbance with water catchment encroachment, disturbance with water catchment excavation, and disturbance with water catchment deposition.

[0010] This embodiment provides a method for calculating soil erosion intensity based on high-resolution data. It divides soil erosion data according to engineering disturbance type and water catchment type, realizing accurate classification of soil erosion data at the project scale, and laying the foundation for subsequent calculation of soil erosion intensity corresponding to different soil loss types.

[0011] In one optional implementation, the soil erosion influencing factors include rainfall erosion factors, runoff scour factors, soil and rock factors, sedimentary mass factors, slope factors, slope length factors, vegetation cover factors, and soil and water conservation measures factors; the soil erosion influencing factors are calculated based on land use type data, soil and water conservation measures data, rainfall data, and surface change data, including: Calculate rainfall erosion factor based on rainfall data; Obtain the runoff coefficient, determine the catchment area and measured rainfall based on rainfall data, and calculate the runoff scour factor based on the runoff coefficient, catchment area and measured rainfall; Based on surface change data, the content of sand, silt, clay and organic carbon were determined, and soil and rock factors were calculated based on the content of sand, silt, clay and organic carbon. Based on surface change data, the depositional material and depositional morphology are determined, and the depositional body factors are determined based on the depositional material and depositional morphology. Slope gradient is determined based on surface change data, and slope factor is calculated based on slope gradient. The slope length is determined based on surface change data, and the slope length factor is calculated based on the slope length and slope gradient. The disturbed surface type is determined based on land use type data, and the vegetation cover factor is determined based on the disturbed surface type. The types of soil and water conservation measures are determined based on the data of soil and water conservation measures, and the factors of soil and water conservation measures are determined based on the types of soil and water conservation measures.

[0012] This embodiment provides a high-resolution data-driven method for calculating soil erosion intensity. Based on different types of soil erosion, it calculates different soil erosion influencing factors, namely rainfall erosion factor, runoff scour factor, soil and rock factor, sedimentary mass factor, slope factor, slope length factor, vegetation cover factor, and soil and water conservation measures factor. This enables multi-scale and multi-dimensional monitoring of surface changes caused by soil erosion, more accurately reflects the actual situation of soil erosion, and improves the accuracy of soil erosion intensity calculation.

[0013] In one optional implementation, soil erosion intensity is calculated based on soil erosion influencing factors according to soil loss type, including: The soil erosion intensity corresponding to the non-catchment occupancy type disturbance is calculated based on rainfall erosion factor, soil and rock factor, slope factor, slope length factor, vegetation cover factor and soil and water conservation measure factor. The soil erosion intensity corresponding to the excavation disturbance without water catchment is calculated based on the rainfall erosion factor, soil and rock factor, slope factor, slope length factor and soil and water conservation measure factor. The soil erosion intensity corresponding to the non-catchment depositional disturbance is calculated based on rainfall erosion factor, soil and rock factor, deposit factor, slope factor, slope length factor and soil and water conservation measure factor. The soil erosion intensity corresponding to the water catchment-occupying disturbance is calculated based on rainfall erosion factor, runoff scour factor, soil and rock factor, slope factor, slope length factor, vegetation cover factor and soil and water conservation measure factor. The soil erosion intensity corresponding to the water catchment excavation type disturbance is calculated based on rainfall erosion factor, runoff scour factor, soil and rock factor, slope factor, slope length factor and soil and water conservation measure factor. Soil erosion intensity corresponding to water-catchment-deposition type disturbances is calculated based on rainfall erosion factor, runoff scour factor, soil and rock factor, deposit factor, slope factor, slope length factor, and soil and water conservation measure factor.

[0014] This embodiment provides a high-resolution data-driven method for calculating soil erosion intensity. It calculates soil erosion intensity for different types of soil loss, improving the accuracy of soil erosion intensity calculation and more accurately reflecting the actual situation of soil erosion. This provides a scientific basis for serving small-scale soil erosion prevention and control and ecological environment construction.

[0015] In one alternative implementation, it further includes: Soil erosion intensity is used to determine the amount of soil loss, and soil erosion control strategies are optimized based on the amount of soil loss.

[0016] This embodiment provides a method for calculating soil erosion intensity based on high-resolution data. The method has good scalability and versatility, and can be applied to soil erosion intensity monitoring and assessment in different regions, providing reliable technical support for soil erosion prevention and control decisions.

[0017] Secondly, the present invention provides a soil erosion intensity measurement device based on high-resolution data, the device comprising: The segmentation module is used to obtain soil and water loss data corresponding to the target production and construction project, and to segment the soil and water loss data according to the engineering disturbance type and the water catchment type to obtain the soil loss type; The determination module is used to acquire UAV operation data and disturbance survey data corresponding to the target production and construction project, and to determine land use type data, soil and water conservation measures data, rainfall data and surface change data based on the UAV operation data and disturbance survey data; The first calculation module is used to calculate the soil erosion impact factor based on land use type data, soil and water conservation measures data, rainfall data, and surface change data. The second calculation module is used to calculate soil erosion intensity according to soil loss type and based on soil erosion influencing factors.

[0018] Thirdly, the present invention provides an electronic device, comprising: a memory and a processor, wherein the memory and the processor are communicatively connected to each other, the memory stores computer instructions, and the processor executes the computer instructions to perform a method for calculating soil erosion intensity based on high-resolution data driven by the first aspect or any corresponding embodiment described above.

[0019] Fourthly, the present invention provides a computer-readable storage medium storing computer instructions for causing a computer to execute a method for calculating soil erosion intensity based on high-resolution data in the first aspect or any corresponding embodiment described above.

[0020] Fifthly, the present invention provides a computer program product, including computer instructions for causing a computer to execute a high-resolution data-driven method for calculating soil erosion intensity based on the first aspect or any corresponding embodiment described above. Attached Figure Description

[0021] To more clearly illustrate the specific embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the specific embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of the present invention. For those skilled in the art, other drawings can be obtained from these drawings without creative effort.

[0022] Figure 1 This is a schematic diagram of an application scenario according to an embodiment of the present invention; Figure 2 This is a schematic diagram of the first process of a method for calculating soil erosion intensity based on high-resolution data according to an embodiment of the present invention. Figure 3 This is a schematic diagram of a second process for a soil erosion intensity calculation method based on high-resolution data according to an embodiment of the present invention; Figure 4 This is a schematic diagram of the third process of a soil erosion intensity measurement method based on high-resolution data according to an embodiment of the present invention. Figure 5 This is a schematic diagram of the fourth process of a soil erosion intensity measurement method based on high-resolution data according to an embodiment of the present invention; Figure 6 This is a schematic diagram of the process for calculating human-induced soil erosion according to an embodiment of the present invention; Figure 7 This is a structural block diagram of a soil erosion intensity measurement device based on high-resolution data according to an embodiment of the present invention; Figure 8 This is a schematic diagram of the hardware structure of an electronic device according to an embodiment of the present invention. Detailed Implementation

[0023] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0024] It is understood that before using the technical solutions disclosed in the various embodiments of the present invention, users should be informed of the types, scope of use, and usage scenarios of the personal information involved in the present invention and their authorization should be obtained in accordance with relevant laws and regulations through appropriate means.

[0025] The terms "first" and "second" are used for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include one or more of that feature. In the description of this invention, "a plurality of" means two or more, unless otherwise explicitly specified.

[0026] As an optional application scenario of this invention, such as Figure 1 As shown, the soil erosion intensity measurement device based on high-resolution data may include at least one terminal device and at least one server. Figure 1 The system is illustrated in the example, which includes a computer 101, a mobile terminal 102, and a server 103, and the terminal devices such as the computer 101 and the mobile terminal 102 are connected to the server 103 through a network 110.

[0027] Specifically, the terminal device can be a smartphone, tablet, laptop, PDA, desktop computer, game console, smart TV, smart wearable device, in-vehicle terminal, VR (Virtual Reality) device, AR (Augmented Reality) device, etc. Server 103 can be a standalone physical server, a server cluster, a distributed system, or a cloud server providing cloud services. Network 110 can be a wired or wireless network, examples of which include, but are not limited to, the Internet, corporate intranet, local area network, wide area network, mobile communication network, and combinations thereof.

[0028] Some dynamic monitoring and assessment systems for soil erosion only focus on the main watershed and county scale (usually larger than the scale of the main watershed and county). ), neglecting the project scale (usually smaller than) Development and application of soil and water loss models.

[0029] In project-scale soil erosion monitoring, relevant monitoring methods mainly rely on a limited number of ground monitoring stations, acquiring data through on-site sampling, manual observation, and simple instrument measurements. However, this approach has serious limitations. On the one hand, due to the limited number of monitoring stations, it is difficult to comprehensively cover the monitoring area, resulting in spatially unrepresentative data that cannot accurately reflect the spatial distribution characteristics of regional soil erosion. On the other hand, the data collection frequency of these monitoring methods is low, failing to capture the dynamic changes in soil erosion intensity over time, and thus failing to meet the needs of real-time monitoring and accurate assessment. Furthermore, these monitoring methods are also limited by natural conditions such as terrain and transportation. In areas with complex terrain and inconvenient transportation, monitoring work is difficult to carry out, and data acquisition costs are high.

[0030] Regarding the process of soil erosion measurement, the method based on soil erosion prediction guidelines also has many problems. This method usually requires a series of tedious and complicated steps, involving a large amount of data collection, processing and analysis. This not only wastes a lot of manpower, material resources and time, but also easily introduces human error in the data processing process, which greatly reduces the accuracy of the measurement results. At the same time, since soil erosion prediction guidelines are often based on certain assumptions and empirical models, they are not adaptable enough to factors such as soil type, topography, climate conditions and land use in different regions, and it is difficult to accurately reflect the actual soil erosion situation in the local area.

[0031] With the continuous advancement of technology, the integration of drones with multiple sensors and remote sensing monitoring technologies has opened up new avenues for refined monitoring of soil erosion. Drones acquire high-resolution imagery data at low altitudes, clearly reflecting subtle changes in the land surface. Simultaneously, they can carry multispectral sensors to accurately identify the spectral characteristics of different land features, providing data for analyzing key indicators such as vegetation cover and soil moisture content. LiDAR sensors can quickly acquire high-precision topographic data, effectively supplementing the deficiencies of Digital Elevation Model (DEM) data and assisting in the identification of soil erosion risk areas. The data acquired by these sensors, combined with the rich surface information obtained from macroscopic monitoring of large areas by remote sensing technology, enables multi-scale, multi-dimensional monitoring of soil erosion, providing comprehensive and abundant data support for the accurate calculation of soil erosion intensity. However, currently, there is no mature method for calculating soil erosion intensity that fully utilizes the advantages of drone multi-sensor and high-resolution remote sensing data. How to effectively integrate and utilize this high-resolution data to establish an accurate and efficient soil erosion intensity calculation model has become a current research hotspot and challenge in this field.

[0032] This invention provides a method for calculating soil erosion intensity based on high-resolution data. It utilizes a drone equipped with different payloads to acquire high-resolution drone imagery, multispectral data, and point cloud data. Through remote sensing interpretation, multispectral data calculation, and spatial triangulation, basic geographic information data, vegetation cover index, and DEM are obtained. Multiple factors, such as rainfall, topography, vegetation, and soil and water conservation measures, are determined and generated. Soil erosion intensity is calculated based on CSLE (the Coupled Soil Landscape Evolution Equation, a model framework for studying the development and changes of land surface landscapes). This method improves the accuracy and spatial resolution of soil erosion intensity calculation, more accurately reflects the actual situation of soil erosion, and provides a scientific basis for serving small-scale soil erosion prevention and control and ecological environment construction.

[0033] According to an embodiment of the present invention, a method for measuring soil erosion intensity based on high-resolution data is provided. It should be noted that the steps shown in the flowchart in the accompanying drawings can be executed in a computer system such as a set of computer-executable instructions. Furthermore, although a logical order is shown in the flowchart, in some cases, the steps shown or described may be executed in a different order than that shown here.

[0034] This embodiment provides a method for calculating soil erosion intensity based on high-resolution data, which can be used in the aforementioned terminal equipment. Figure 2 This is a flowchart of a soil erosion intensity measurement method based on high-resolution data according to an embodiment of the present invention, as shown below. Figure 2 As shown, the process includes the following steps: Step S201: Obtain soil and water loss data corresponding to the target production and construction project, and classify the soil and water loss data according to the engineering disturbance type and the water catchment type to obtain the soil loss type.

[0035] Specifically, retrieve data on the scope of responsibility for soil and water conservation, prevention and control zones, and soil and water conservation measures (SHP data, a vector data format of geographic information systems) for production and construction projects, and divide the disturbance units according to the soil loss type of the production and construction projects.

[0036] Step S202: Obtain UAV operation data and disturbance survey data corresponding to the target production and construction project, and determine land use type data, soil and water conservation measures data, rainfall data and surface change data based on the UAV operation data and disturbance survey data.

[0037] Specifically, UAV operational data mainly includes high-resolution UAV images obtained by using UAVs equipped with oblique photography payloads, multispectral data obtained by using UAVs equipped with multispectral remote sensing payloads, and laser point cloud data obtained by using lidar payloads.

[0038] Furthermore, the disturbance investigation work includes investigating the soil texture, vegetation cover, forest canopy closure, type of soil and water conservation measures, and morphology of engineering accumulation (determining whether it is an inclined slope or a conical accumulation) of each disturbance unit; taking soil samples from each disturbance unit on the excavation face to determine the material composition and properties such as soil density, mechanical composition, and gravel content; and when the land use type of the disturbance unit is forest land, using photographic or visual estimation methods to measure the canopy closure on-site.

[0039] Furthermore, land use type data, soil and water conservation measures data, rainfall data, and surface change data are the basic data for calculating the influencing factors of soil and water loss.

[0040] Step S203: Calculate the soil erosion influencing factors based on land use type data, soil and water conservation measures data, rainfall data, and surface change data.

[0041] Step S204: Calculate the soil erosion intensity based on the soil loss influencing factors according to the soil loss type.

[0042] Specifically, the amount of soil loss is determined based on the soil erosion intensity, and the soil loss control strategy is optimized based on the amount of soil loss.

[0043] Furthermore, soil erosion intensity is used to estimate soil and water loss, or to determine whether an area needs to be protected against water loss by judging whether the soil erosion intensity exceeds a certain threshold.

[0044] This embodiment provides a method for calculating soil erosion intensity based on high-resolution data. It acquires UAV operation data and disturbance survey data corresponding to the target production and construction project, utilizes UAV remote sensing technology to predict and monitor soil erosion at the project scale, and then determines land use type data, soil and water conservation measures data, rainfall data, and surface change data based on the UAV operation data and disturbance survey data. Based on these data, it calculates soil erosion influencing factors. By fully utilizing the rich information in high-resolution data, it can more precisely depict the spatial variations of various soil erosion influencing factors, significantly improve the accuracy and spatial resolution of soil erosion intensity calculation, and more accurately reflect the actual situation of soil erosion, providing a scientific basis for serving small-scale soil erosion prevention and control and ecological environment construction.

[0045] This embodiment provides a method for calculating soil erosion intensity based on high-resolution data, which can be used in the aforementioned terminal equipment. Figure 3 This is a flowchart of a soil erosion intensity measurement method based on high-resolution data according to an embodiment of the present invention, as shown below. Figure 3 As shown, the process includes the following steps: Step S301: Obtain soil erosion data corresponding to the target production and construction project. Classify the soil erosion data according to the type of engineering disturbance and the type of water catchment to obtain the soil erosion type. For details, please refer to [link to relevant documentation]. Figure 2 Step S201 of the illustrated embodiment will not be described again here.

[0046] Step S302: Obtain UAV operation data and disturbance survey data corresponding to the target production and construction project, and determine land use type data, soil and water conservation measures data, rainfall data and surface change data based on the UAV operation data and disturbance survey data.

[0047] Specifically, step S302 includes: Step S3021: Determine UAV image data, multispectral data, and laser point cloud data based on UAV operation data.

[0048] Step S3022: Preprocess the UAV image data, multispectral data and laser point cloud data respectively to obtain digital surface model data, multispectral image data and digital elevation model data.

[0049] Specifically, the preprocessing of UAV imagery data, multispectral data, and laser point cloud data includes: obtaining true radiometric DSM (Digital Surface Model) data from UAV imagery data through differential calculation and aerial triangulation; obtaining multispectral images through geometric correction, atmospheric correction, image registration, and band fusion; and obtaining DEM (Digital Elevation Model) data from laser point cloud data.

[0050] Furthermore, the acquired UAV image data undergoes preprocessing operations such as distortion correction, coordinate transformation, and image stitching to improve image quality; the point cloud data acquired by the UAV-mounted LiDAR is processed by filtering and denoising, ground point extraction, and digital surface model construction, and then combined with high-precision digital elevation model data to extract terrain change information through differential calculation and other methods.

[0051] Step S3023: Based on digital surface model data, multispectral image data, digital elevation model data, and disturbance survey data, determine land use type data, soil and water conservation measures data, rainfall data, and surface change data, respectively.

[0052] Specifically, the calculation steps for land use type shp (i.e., land use type data) include: using true orthophoto maps and visual remote sensing interpretation methods to sequentially interpret land use types such as farmland, trees, shrubs, grassland, hard ground, water bodies, and bare land; at the same time, referring to the national crop rotation area to assign values ​​to tillage factors; and filling in the vegetation cover and canopy closure field information according to the survey results of disturbance units. The attribute table of land use type shp fields is shown in Table 1 below.

[0053] Table 1:

[0054] Furthermore, the calculation steps for soil and water conservation measures shp (i.e., soil and water conservation measures data) include: interpreting the types of soil and water conservation measures using true orthophoto maps with visual remote sensing interpretation combined with the results of disturbance unit surveys. The types of soil and water conservation measures can be divided into sediment trapping projects, slope protection projects, land reclamation projects, flood control and drainage projects, rainfall storage and infiltration projects, temporary protection projects, vegetation construction projects, and windbreak and sand fixation projects. The attribute table of soil and water conservation measures shp fields is shown in Table 2 below.

[0055] Table 2:

[0056] Furthermore, the rainfall Excel (i.e., rainfall data) is compiled according to the daily rainfall.

[0057] Furthermore, the calculation unit shapefile (i.e., surface change data) adds information fields such as soil texture and sediment morphology to the disturbance unit. According to the disturbance survey data, the relevant field information is filled in. The field attribute table of the calculation unit shapefile is shown in Report 3 below.

[0058] Table 3:

[0059] Step S303: Calculate the soil erosion influencing factors based on land use type data, soil and water conservation measures data, rainfall data, and surface change data. For details, please refer to [link to relevant documentation]. Figure 2 Step S203 of the illustrated embodiment will not be described again here.

[0060] Step S304: Calculate soil erosion intensity based on soil erosion influencing factors according to soil loss type. For details, please refer to [link to relevant documentation]. Figure 2 Step S204 of the illustrated embodiment will not be described again here.

[0061] This embodiment provides a method for calculating soil erosion intensity based on high-resolution data. It uses UAV remote sensing technology to predict and monitor soil and water loss at the project scale. For large-scale production and construction projects, it overcomes the limitations of insufficient accuracy, high difficulty, and high cost of monitoring methods.

[0062] This embodiment provides a method for calculating soil erosion intensity based on high-resolution data, which can be used in the aforementioned terminal equipment. Figure 4 This is a flowchart of a soil erosion intensity measurement method based on high-resolution data according to an embodiment of the present invention, as shown below. Figure 4 As shown, the process includes the following steps: Step S401: Obtain soil erosion data corresponding to the target production and construction project. Classify the soil erosion data according to the type of engineering disturbance and the type of water catchment to obtain the soil erosion type. For details, please refer to [link to relevant documentation]. Figure 3 Step S301 of the illustrated embodiment will not be described again here.

[0063] Step S402: Obtain UAV operation data and disturbance survey data corresponding to the target production and construction project. Based on the UAV operation data and disturbance survey data, determine land use type data, soil and water conservation measures data, rainfall data, and surface change data. For details, please refer to [link to relevant documentation]. Figure 3 Step S302 of the illustrated embodiment will not be described again here.

[0064] Step S403: Calculate the soil erosion influencing factors based on land use type data, soil and water conservation measures data, rainfall data, and surface change data.

[0065] Specifically, the factors influencing soil and water loss include rainfall erosion factors, runoff scour factors, soil and rock factors, sedimentary mass factors, slope factors, slope length factors, vegetation cover factors, and soil and water conservation measures factors; the above step S403 includes: Step S4031: Calculate the rainfall erosion factor based on rainfall data.

[0066] Specifically, the rainfall erosion factor is calculated using measured daily rainfall data. By accumulating the daily rainfall erosion factors, rainfall erosion factors for different time scales, such as multi-day, monthly, multi-month, and yearly, are obtained. The calculation formula is as follows: (1) (2) (3) in, Indicates the first Daily rainfall erosion factor, in units of ; Indicates the first Daily rainfall, in units The daily rainfall must be ≥12mm; otherwise, it will be calculated as 0. and Represents statistical coefficients, dimensionless; This refers to the average daily rainfall of ≥12mm, expressed in units of... ; This represents the annual average rainfall with a daily rainfall of ≥12mm, expressed in units of... .

[0067] Step S4032: Obtain the runoff coefficient, determine the catchment area and measured rainfall based on the rainfall data, and calculate the runoff scour factor based on the runoff coefficient, catchment area and measured rainfall.

[0068] Specifically, the runoff scour factor is a function of three parameters: catchment area, runoff coefficient, and rainfall. The calculation formula is as follows: (4) in, Indicates the catchment area. This represents the runoff coefficient, which is determined by the underlying surface type. This represents the measured rainfall amount, obtained through observation.

[0069] Furthermore, by replacing the concept of catchment area with the cumulative step size in the raster data, a two-dimensional calculation of the runoff scour factor is achieved.

[0070] Step S4033: Determine the sand content, silt content, clay content, and organic carbon content based on surface change data, and calculate the soil-rock factors based on the sand content, silt content, clay content, and organic carbon content.

[0071] Specifically, the soil-rock factor is calculated using the EPIC (Erosion and Soil Productivity Model) formula. The calculation formula is as follows: (5) in, The percentage of sand particles. The content of powder particles (%). The content of clay particles (%). Organic carbon content (%) The content of non-sand particles is SN1 = 1 - SAN / 100.

[0072] Step S4034: Determine the deposited material and deposit morphology based on surface change data, and determine the deposited body factor based on the deposited material and deposit morphology.

[0073] Specifically, the packing factor is determined by comprehensively considering both the packing morphology and the packing material. The reference values ​​are shown in Table 4 below.

[0074] Table 4:

[0075] Step S4035: Determine the slope gradient based on surface change data, and calculate the slope factor based on the slope gradient.

[0076] Specifically, slope factor The calculation formula is as follows: (6) in, This is for calculating the slope of the slope in the unit shapep.

[0077] Step S4036: Determine the slope length based on surface change data, and calculate the slope length factor based on the slope length and slope gradient.

[0078] Specifically, the formula for calculating the slope length factor L is as follows: L (7) (8) in, The horizontal projection slope length, The length of the slope. This is the slope length index.

[0079] Step S4037: Determine the disturbed surface type based on land use type data, and determine the vegetation cover factor based on the disturbed surface type.

[0080] Specifically, vegetation cover factor Interpolation methods were used to obtain the vegetation factors. When the general disturbed surface type was grassland or shrubland, vegetation factors were determined by linear interpolation using vegetation cover. When the general disturbed surface was arbor forest, vegetation factors were determined by linear interpolation using both vegetation cover and forest canopy closure.

[0081] Step S4038: Determine the type of soil and water conservation measures based on the soil and water conservation measures data, and determine the soil and water conservation measures factors based on the type of soil and water conservation measures.

[0082] Specifically, soil and water conservation measures factors The value is obtained by taking the values. According to Table 5 below, the total number of values ​​corresponding to the measure types in the target area is calculated. The measure factor is equal to 1 minus the total number of values. If there are no soil and water conservation measures, the factor value is 1.

[0083] Table 5:

[0084] Step S404: Calculate soil erosion intensity based on soil erosion influencing factors according to soil loss type. For details, please refer to [link to relevant documentation]. Figure 3 Step S304 of the illustrated embodiment will not be described again here.

[0085] This embodiment provides a high-resolution data-driven method for calculating soil erosion intensity. Based on different types of soil erosion, it calculates different soil erosion influencing factors, namely rainfall erosion factor, runoff scour factor, soil and rock factor, sedimentary mass factor, slope factor, slope length factor, vegetation cover factor, and soil and water conservation measures factor. This enables multi-scale and multi-dimensional monitoring of surface changes caused by soil erosion, more accurately reflects the actual situation of soil erosion, and improves the accuracy of soil erosion intensity calculation.

[0086] This embodiment provides a method for calculating soil erosion intensity based on high-resolution data, which can be used in the aforementioned terminal equipment. Figure 5 This is a flowchart of a soil erosion intensity measurement method based on high-resolution data according to an embodiment of the present invention, as shown below. Figure 5 As shown, the process includes the following steps: Step S501: Obtain soil and water loss data corresponding to the target production and construction project, and classify the soil and water loss data according to the engineering disturbance type and the water catchment type to obtain the soil loss type.

[0087] Specifically, step S501 includes: Step S5011: Determine the current geomorphological data and the water catchment data above the target area based on the soil erosion data; determine the engineering disturbance type based on the current geomorphological data; and determine the water catchment type based on the water catchment data above the target area.

[0088] Specifically, based on the current micro-topography (i.e., current geomorphological data), the type of engineering disturbance is determined to be either encroachment disturbance, excavation disturbance, or accumulation disturbance. The specific determination process is shown in Table 6 below.

[0089] Table 6:

[0090] Furthermore, the type of water catchment is determined by whether there is water flowing from above or not, based on whether there is a catchment area above the target area; where the target area is the area involved in the target production and construction project.

[0091] Step S5012: Soil loss types are classified based on engineering disturbance type and water catchment type. Soil loss types include disturbance without water catchment occupancy, disturbance without water catchment excavation, disturbance without water catchment deposition, disturbance with water catchment occupancy, disturbance with water catchment excavation, and disturbance with water catchment deposition.

[0092] Specifically, soil loss is classified into six types according to the type of engineering disturbance and whether there is a water catchment area, as shown in Table 7 below.

[0093] Table 7:

[0094] Step S502: Obtain UAV operation data and disturbance survey data corresponding to the target production and construction project. Based on the UAV operation data and disturbance survey data, determine land use type data, soil and water conservation measures data, rainfall data, and surface change data. For details, please refer to [link to relevant documentation]. Figure 4 Step S402 of the illustrated embodiment will not be described again here.

[0095] Step S503: Calculate the soil erosion influencing factors based on land use type data, soil and water conservation measures data, rainfall data, and surface change data. For details, please refer to [link to relevant documentation]. Figure 4 Step S403 of the illustrated embodiment will not be described again here.

[0096] Step S504: Calculate the soil erosion intensity according to the soil loss type and based on the soil erosion influencing factors.

[0097] Specifically, step S504 includes: Step S5041: Calculate the soil erosion intensity corresponding to the disturbance without water catchment area based on rainfall erosion factor, soil and rock factor, slope factor, slope length factor, vegetation cover factor and soil and water conservation measures factor.

[0098] Specifically, soil erosion intensity was calculated for three types of engineering disturbance without a catchment area. The soil erosion intensity corresponding to the non-catchment occupancy type disturbance was calculated. The calculation formula is as follows: (9) Step S5042: Calculate the soil erosion intensity corresponding to the excavation disturbance without water catchment based on rainfall erosion factor, soil and rock factor, slope factor, slope length factor and soil and water conservation measure factor.

[0099] Specifically, the soil erosion intensity corresponding to excavation-type disturbance without water catchment. The calculation formula is as follows: (10) Step S5043: Calculate the soil erosion intensity corresponding to the non-catchment depositional disturbance based on rainfall erosion factor, soil and rock factor, depositional mass factor, slope factor, slope length factor, and soil and water conservation measures factor.

[0100] Specifically, soil erosion intensity corresponding to disturbances without water catchment and deposition. The calculation formula is as follows: (11) Step S5044: Calculate the soil erosion intensity corresponding to the water catchment-occupying disturbance based on rainfall erosion factor, runoff scour factor, soil and rock factor, slope factor, slope length factor, vegetation cover factor, and soil and water conservation measures factor.

[0101] Specifically, considering the water catchment area, the soil erosion of the above three types of engineering disturbances needs to be increased by the runoff scouring factor.

[0102] Furthermore, the soil erosion intensity corresponding to water catchment-occupancy type disturbance. The calculation formula is as follows: (12) Step S5045: Calculate the soil erosion intensity corresponding to the water catchment excavation type disturbance based on rainfall erosion factor, runoff scour factor, soil and rock factor, slope factor, slope length factor and soil and water conservation measure factor.

[0103] Specifically, there is soil erosion intensity corresponding to water catchment excavation disturbance. The calculation formula is as follows: (13) Step S5046: Calculate the soil erosion intensity corresponding to the water-catchment-deposition type disturbance based on rainfall erosion factor, runoff scour factor, soil and rock factor, deposit factor, slope factor, slope length factor, and soil and water conservation measures factor.

[0104] Specifically, there is soil erosion intensity corresponding to water-collecting and depositional disturbances. The calculation formula is as follows: (14) This embodiment provides a high-resolution data-driven method for calculating soil erosion intensity. It calculates soil erosion intensity for different types of soil loss, improving the accuracy of soil erosion intensity calculation and more accurately reflecting the actual situation of soil erosion. This provides a scientific basis for serving small-scale soil erosion prevention and control and ecological environment construction.

[0105] The following specific embodiment illustrates the steps of a method for calculating soil erosion intensity based on high-resolution data.

[0106] Example 1: like Figure 6 As shown, the specific steps of a soil erosion intensity measurement method based on high-resolution data include: S1. Preparation work: retrieve the shp of the responsibility area for soil and water conservation of the production and construction project, divide the disturbance units according to the soil loss type of the production and construction project, and determine the scope of drone operation and field operation.

[0107] S2. Fieldwork: including drone operations and disturbance unit surveys.

[0108] S3. UAV Data Production and Preprocessing: This includes UAV remote sensing data production and preprocessing. Preprocessing operations such as distortion correction, coordinate transformation, and image stitching are performed on the acquired UAV image data to improve image quality. Point cloud data acquired by the UAV-mounted LiDAR undergoes filtering and noise reduction, ground point extraction, and digital surface model (DSM) construction. This data is then combined with high-precision digital elevation model (DEM) data to extract terrain change information through differential calculations and other methods.

[0109] S4. Erosion Factor Calculation: Calculate the various erosion factors for different soil loss types. Common erosion factors include rainfall erosion factor, soil and rock factor, slope length factor, slope factor, vegetation cover factor, and soil and water conservation measure factor. For engineering accumulations, the accumulation factor needs to be calculated. In the case of water catchment, the runoff scour factor needs to be calculated.

[0110] S5. Modulus Calculation: Calculate the soil erosion intensity for different types of soil loss, including disturbance without catchment area encroachment, disturbance without catchment area excavation, disturbance without catchment area deposition, disturbance with catchment area encroachment, disturbance with catchment area excavation, and disturbance with catchment area deposition.

[0111] This embodiment also provides a soil erosion intensity calculation device based on high-resolution data, which is used to implement the above embodiments and preferred embodiments; details already described will not be repeated. As used below, the term "module" can refer to a combination of software and / or hardware that performs a predetermined function. Although the device described in the following embodiments is preferably implemented in software, hardware implementation, or a combination of software and hardware, is also possible and contemplated.

[0112] This embodiment provides a soil erosion intensity measurement device based on high-resolution data, such as... Figure 7 As shown, it includes: The segmentation module 701 is used to obtain soil and water loss data corresponding to the target production and construction project, and to segment the soil and water loss data according to the engineering disturbance type and the water catchment type to obtain the soil loss type; The determination module 702 is used to acquire UAV operation data and disturbance survey data corresponding to the target production and construction project, and to determine land use type data, soil and water conservation measures data, rainfall data and surface change data based on the UAV operation data and disturbance survey data; The first calculation module 703 is used to calculate the soil erosion impact factor based on land use type data, soil and water conservation measures data, rainfall data and surface change data; The second calculation module 704 is used to calculate the soil erosion intensity according to the soil loss type and based on the soil erosion influencing factors.

[0113] In some alternative implementations, the determining module 702 includes: The first determining unit is used to determine UAV image data, multispectral data, and laser point cloud data based on UAV operation data; The preprocessing unit is used to preprocess UAV image data, multispectral data and laser point cloud data respectively to obtain digital surface model data, multispectral image data and digital elevation model data. The second determining unit is used to determine land use type data, soil and water conservation measures data, rainfall data, and surface change data based on digital surface model data, multispectral image data, digital elevation model data, and disturbance survey data, respectively.

[0114] In some alternative implementations, the partitioning module 701 includes: The third determining unit is used to determine the current geomorphological data and the water catchment data above the target area based on the soil erosion data, determine the engineering disturbance type based on the current geomorphological data, and determine the water catchment type based on the water catchment data above the target area. The division unit is used to classify soil loss types based on engineering disturbance type and water catchment type. Soil loss types include disturbance without water catchment occupancy, disturbance without water catchment excavation, disturbance without water catchment deposition, disturbance with water catchment occupancy, disturbance with water catchment excavation, and disturbance with water catchment deposition.

[0115] In some alternative implementations, the first computing module 703 includes: The first calculation unit is used to calculate the rainfall erosion factor based on rainfall data; The second calculation unit is used to obtain the runoff coefficient, determine the catchment area and measured rainfall based on the rainfall data, and calculate the runoff scour factor based on the runoff coefficient, catchment area and measured rainfall. The third calculation unit is used to determine the sand content, silt content, clay content and organic carbon content based on surface change data, and to calculate soil and rock factors based on the sand content, silt content, clay content and organic carbon content; The fourth determining unit is used to determine the depositional material and depositional morphology based on surface change data, and to determine the depositional body factor based on the depositional material and depositional morphology. The fourth calculation unit is used to determine the slope gradient based on surface change data and to calculate the slope factor based on the slope gradient. The fifth calculation unit is used to determine the slope length based on surface change data, and to calculate the slope length factor based on the slope length and slope gradient; The fifth determining unit is used to determine the disturbed surface type based on land use type data, and to determine the vegetation cover factor based on the disturbed surface type; The sixth determining unit is used to determine the type of soil and water conservation measures based on the data of soil and water conservation measures, and to determine the soil and water conservation measure factors based on the type of soil and water conservation measures.

[0116] In some alternative implementations, the second computing module 704 includes: The sixth calculation unit is used to calculate the soil erosion intensity corresponding to the disturbance without water catchment area based on rainfall erosion factor, soil and rock factor, slope factor, slope length factor, vegetation cover factor and soil and water conservation measure factor. The seventh calculation unit is used to calculate the soil erosion intensity corresponding to the excavation disturbance without water catchment based on the rainfall erosion factor, soil and rock factor, slope factor, slope length factor and soil and water conservation measure factor. The eighth calculation unit is used to calculate the soil erosion intensity corresponding to non-catchment depositional disturbance based on rainfall erosion factor, soil and rock factor, deposit factor, slope factor, slope length factor and soil and water conservation measure factor. The ninth calculation unit is used to calculate the soil erosion intensity corresponding to the water catchment-occupying type disturbance based on rainfall erosion factor, runoff scour factor, soil and rock factor, slope factor, slope length factor, vegetation cover factor and soil and water conservation measure factor. The tenth calculation unit is used to calculate the soil erosion intensity corresponding to the water catchment excavation type disturbance based on rainfall erosion factor, runoff scour factor, soil and rock factor, slope factor, slope length factor and soil and water conservation measure factor. The eleventh calculation unit is used to calculate the soil erosion intensity corresponding to water-catchment-deposition type disturbances based on rainfall erosion factor, runoff scour factor, soil and rock factor, deposit factor, slope factor, slope length factor, and soil and water conservation measure factor.

[0117] In some alternative implementations, it also includes: The optimization module is used to determine the amount of soil loss based on soil erosion intensity and optimize soil loss control strategies based on the amount of soil loss.

[0118] The soil erosion intensity calculation device based on high-resolution data provided in this embodiment of the invention can execute the soil erosion intensity calculation method based on high-resolution data provided in any embodiment of the invention, and has the corresponding functional modules and beneficial effects of the method. Further functional descriptions of the above modules and units are the same as those in the corresponding embodiments above, and will not be repeated here.

[0119] Figure 8 This is a schematic diagram of the structure of an electronic device provided in an embodiment of the present invention.

[0120] The following is a detailed reference. Figure 8 This diagram illustrates a suitable structural schematic for implementing an electronic device according to embodiments of the present invention. The electronic device may include a processor (e.g., a central processing unit, graphics processor, etc.) 801, which can perform various appropriate actions and processes based on a program stored in read-only memory (ROM) 802 or a program loaded from memory 808 into random access memory (RAM) 803. The RAM 803 also stores various programs and data required for the operation of the electronic device. The processor 801, ROM 802, and RAM 803 are interconnected via a bus 804. An input / output (I / O) interface 805 is also connected to the bus 804.

[0121] Typically, the following devices can be connected to I / O interface 805: input devices 806 including, for example, touchscreens, touchpads, keyboards, mice, cameras, microphones, accelerometers, gyroscopes, etc.; output devices 807 including, for example, liquid crystal displays (LCDs), speakers, vibrators, etc.; memory devices 808 including, for example, magnetic tapes, hard disks, etc.; and communication devices 809. Communication device 809 allows electronic devices to communicate wirelessly or wiredly with other devices to exchange data. Although Figure 8 Electronic devices with various devices are shown, but it should be understood that it is not required to implement or have all of the devices shown, and more or fewer devices may be implemented or have instead.

[0122] In particular, according to embodiments of the present invention, the processes described above with reference to the flowcharts can be implemented as computer software programs. For example, embodiments of the present invention include a computer program product comprising a computer program carried on a non-transitory computer-readable medium, the computer program containing program code for performing the methods shown in the flowcharts. In such embodiments, the computer program can be downloaded and installed from a network via a communication device 809, or installed from a memory 808, or installed from a ROM 802. When the computer program is executed by a processor 801, it performs the functions defined in the high-resolution data-driven soil erosion intensity measurement method of the embodiments of the present invention.

[0123] Figure 8 The electronic device shown is merely an example and should not be construed as limiting the functionality and scope of use of the embodiments of the present invention.

[0124] This invention also provides a computer-readable storage medium. The methods described above according to embodiments of the invention can be implemented in hardware or firmware, or implemented as computer code that can be recorded on a storage medium, or implemented as computer code downloaded via a network and originally stored on a remote storage medium or a non-transitory machine-readable storage medium and then stored on a local storage medium. Thus, the methods described herein can be processed by software stored on a storage medium using a general-purpose computer, a dedicated processor, or programmable or dedicated hardware. The storage medium can be a magnetic disk, optical disk, read-only memory, random access memory, flash memory, hard disk, or solid-state drive, etc.; further, the storage medium can also include combinations of the above types of memory. It is understood that computers, processors, microprocessor controllers, or programmable hardware include storage components capable of storing or receiving software or computer code. When the software or computer code is accessed and executed by the computer, processor, or hardware, it implements the high-resolution data-driven soil erosion intensity calculation method shown in the above embodiments.

[0125] A portion of this invention can be applied as a computer program product, such as computer program instructions, which, when executed by a computer, can invoke or provide the methods and / or technical solutions according to the invention through the operation of the computer. Those skilled in the art will understand that the forms in which computer program instructions exist in a computer-readable medium include, but are not limited to, source files, executable files, installation package files, etc. Correspondingly, the ways in which computer program instructions are executed by a computer include, but are not limited to: the computer directly executing the instructions, or the computer compiling the instructions and then executing the corresponding compiled program, or the computer reading and executing the instructions, or the computer reading and installing the instructions and then executing the corresponding installed program. Here, the computer-readable medium can be any available computer-readable storage medium or communication medium accessible to a computer.

[0126] Although embodiments of the invention have been described in conjunction with the accompanying drawings, those skilled in the art can make various modifications and variations without departing from the spirit and scope of the invention, and such modifications and variations all fall within the scope defined by the appended claims.

Claims

1. A high-resolution data-driven based soil erosion intensity estimation method, characterized in that, The method comprises: obtaining water and soil loss data corresponding to the target production and construction project, dividing the water and soil loss data according to engineering disturbance types and catchment types to obtain soil loss types; obtaining unmanned aerial vehicle operation data and disturbance investigation data corresponding to the target production and construction project, and determining land use type data, water and soil conservation measure data, rainfall data and surface change data based on the unmanned aerial vehicle operation data and the disturbance investigation data; calculating water and soil loss influence factors based on the land use type data, the water and soil conservation measure data, the rainfall data and the surface change data; calculating soil erosion intensity based on the water and soil loss influence factors according to the soil loss types.

2. The method of claim 1, wherein, The method comprises: determining unmanned aerial vehicle image data, multispectral data and laser point cloud data based on the unmanned aerial vehicle operation data; respectively pre-processing the unmanned aerial vehicle image data, the multispectral data and the laser point cloud data to obtain digital surface model data, multispectral image data and digital elevation model data; respectively determining the land use type data, the water and soil conservation measure data, the rainfall data and the surface change data based on the digital surface model data, the multispectral image data, the digital elevation model data and the disturbance investigation data.

3. The method of claim 1, wherein, The method comprises: determining current landform data and catchment data above a target region based on the water and soil loss data, determining the engineering disturbance types based on the current landform data, and determining the catchment types based on the catchment data above the target region; dividing the soil loss types based on the engineering disturbance types and the catchment types, wherein the soil loss types include no-catchment pressure-type disturbance, no-catchment excavation-type disturbance, no-catchment accumulation-type disturbance, catchment pressure-type disturbance, catchment excavation-type disturbance and catchment accumulation-type disturbance.

4. The method of claim 3, wherein, The water and soil loss influence factors include rainfall erosion factor, runoff scouring factor, soil and rock quality factor, accumulation body factor, slope factor, slope length factor, vegetation coverage factor and water and soil conservation measure factor; the method comprises: calculating the rainfall erosion factor based on the rainfall data; obtaining a runoff coefficient, determining catchment area and measured rainfall based on the rainfall data, and calculating the runoff scouring factor based on the runoff coefficient, the catchment area and the measured rainfall; determining sand content, silt content, clay content and organic carbon content based on the surface change data, and calculating the soil and rock quality factor based on the sand content, the silt content, the clay content and the organic carbon content; determine accumulation material and accumulation form based on the ground surface change data, and determine the accumulation body factor based on the accumulation material and the accumulation form; determine slope gradient based on the ground surface change data, and calculate the slope factor based on the slope gradient; determine slope length based on the ground surface change data, and calculate the slope length factor based on the slope length and the slope gradient; determine disturbed ground surface type based on the land use type data, and determine the vegetation coverage factor based on the disturbed ground surface type; determine water and soil conservation measure type based on the water and soil conservation measure data, and determine the water and soil conservation measure factor based on the water and soil conservation measure type.

5. The method of claim 4, wherein, The soil erosion intensity is calculated based on the water and soil erosion influencing factor according to the soil erosion type, including: calculating the soil erosion intensity corresponding to the no-confluence pressure type disturbance based on the rainfall erosion factor, the soil and rock quality factor, the slope factor, the slope length factor, the vegetation coverage factor and the water and soil conservation measure factor; calculating the soil erosion intensity corresponding to the no-confluence excavation type disturbance based on the rainfall erosion factor, the soil and rock quality factor, the slope factor, the slope length factor and the water and soil conservation measure factor; calculating the soil erosion intensity corresponding to the no-confluence accumulation type disturbance based on the rainfall erosion factor, the soil and rock quality factor, the accumulation body factor, the slope factor, the slope length factor and the water and soil conservation measure factor; calculating the soil erosion intensity corresponding to the confluence pressure type disturbance based on the rainfall erosion factor, the runoff scouring factor, the soil and rock quality factor, the slope factor, the slope length factor, the vegetation coverage factor and the water and soil conservation measure factor; calculating the soil erosion intensity corresponding to the confluence excavation type disturbance based on the rainfall erosion factor, the runoff scouring factor, the soil and rock quality factor, the slope factor, the slope length factor and the water and soil conservation measure factor; calculating the soil erosion intensity corresponding to the confluence accumulation type disturbance based on the rainfall erosion factor, the runoff scouring factor, the soil and rock quality factor, the accumulation body factor, the slope factor, the slope length factor and the water and soil conservation measure factor.

6. The method of claim 1, wherein, Further comprising: determining the amount of water and soil loss based on the soil erosion intensity, and optimizing the water and soil loss treatment strategy based on the amount of water and soil loss.

7. A device for measuring soil erosion intensity based on high-resolution data driving, characterized in that, The device comprises: a division module configured to obtain water and soil loss data corresponding to a target production and construction project, divide the water and soil loss data according to engineering disturbance types and confluence types, and obtain soil erosion types; a determination module configured to obtain unmanned aerial vehicle operation data and disturbance investigation data corresponding to the target production and construction project, determine land use type data, water and soil conservation measure data, rainfall data and ground surface change data based on the unmanned aerial vehicle operation data and the disturbance investigation data; a first calculation module configured to calculate water and soil erosion influencing factors based on the land use type data, the water and soil conservation measure data, the rainfall data and the ground surface change data; The second calculation module is configured to calculate soil erosion intensity based on the water and soil loss influencing factors according to the soil loss types.

8. An electronic device, comprising: The application further provides a computer readable storage medium having computer instructions stored thereon, wherein the computer instructions are used to make a computer execute the soil erosion intensity calculation method based on high-resolution data driving. The memory and the processor are in communication connection with each other, and the memory stores computer instructions.

9. A computer-readable storage medium, characterized in that, The computer readable storage medium has computer instructions stored thereon, and the computer instructions are used to make a computer execute the soil erosion intensity calculation method based on high-resolution data driving.

10. A computer program product, characterised in that, The computer readable storage medium has computer instructions stored thereon, and the computer instructions are used to make a computer execute the soil erosion intensity calculation method based on high-resolution data driving.

Images