Method for monitoring soil erosion rate in wind- water complex erosion region

Patent metadata

Field	Value
Publication number	US-12541971-B2
Application number	US-202519215867-A
Country	US
Kind code	B2
Filing date	May 22, 2025
Priority date	Jun 27, 2024
Publication date	Feb 3, 2026
Grant date	Feb 3, 2026

Abstract

Disclosed is a method for monitoring soil erosion rate in a wind-water complex erosion region, which includes dividing a soil erosion monitoring area into a bare area and a vegetated area, performing grid division, and obtaining soil erodibility of each grid in the bare area and the vegetated area; obtaining a bare area reference grid and a vegetated area reference grid; obtaining image data of both the bare area and the vegetated area at monitoring starting and ending moments; calculating soil erosion volumes of both the bare area reference grid and the vegetated area reference grid; calculating a soil erosion volume corrected by soil erodibility; and calculating the soil erosion rate according to the soil erosion volume corrected by soil erodibility, and completing the monitoring of the soil erosion rate. The problem that the workload and monitoring accuracy of the existing monitoring methods cannot be balanced is solved.

First claim

What is claimed is:

1 . A method for monitoring soil erosion rate in a wind-water complex erosion region, comprising the following steps:
S 1 : determining a soil erosion monitoring area; S 2 : at a monitoring starting moment, dividing the soil erosion monitoring area into a bare area and a vegetated area based on a surface cover condition, performing grid division on both the bare area and the vegetated area, and randomly sampling soil in each grid to obtain soil erodibility of each grid in the bare area and the vegetated area; S 3 : according to the soil erodibility of each grid in the bare area and the vegetated area, taking a grid with maximum soil erodibility of the bare area and a grid with maximum soil erodibility of the vegetated area as reference grids to obtain a bare area reference grid and a vegetated area reference grid, and inserting erosion pins into the bare area reference grid and the vegetated area reference grid; S 4 : acquiring image data of the soil erosion monitoring area at the monitoring starting moment and a monitoring ending moment by using an unmanned aerial vehicle, and segmenting the image data of the soil erosion monitoring area at the monitoring starting moment and the monitoring ending moment based on the bare area and the vegetated area to respectively obtain image data of the bare area at the monitoring starting moment and the monitoring ending moment and image data of the vegetated area at the monitoring starting moment and the monitoring ending moment; S 5 : calculating a soil erosion volume of the bare area reference grid and a soil erosion volume of the vegetated area reference grid according to the erosion pins in the bare area reference grid and the vegetated area reference grid; S 6 : calculating a soil erosion volume corrected by soil erodibility according to the image data of the bare area at the monitoring starting moment and the monitoring ending moment, the image data of the vegetated area at the monitoring starting moment and the monitoring ending moment, the soil erosion volume of the bare area reference grid and the soil erosion volume of the vegetated area reference grid; wherein the S 6 specifically is as follows: S 601 : performing gray processing on the image data of the bare area at the monitoring starting moment and the monitoring ending moment and the image data of the vegetated area at the monitoring starting moment and the monitoring ending moment to obtain gray images of the bare area at the monitoring starting moment and the monitoring ending moment and gray images of the vegetated area at the monitoring starting moment and the monitoring ending moment; S 602 : according to a gray image of the bare area at the monitoring starting moment and a gray image of the bare area at the monitoring ending moment, subtracting an average gray value at a grid ending moment from an average gray value at a grid starting moment to obtain a gray difference of each grid in the bare area; S 603 : according to a gray image of the vegetated area at the monitoring starting moment and a gray image of the vegetated area at the monitoring ending moment, subtracting an average gray value at a grid ending moment from an average gray value at a grid starting moment to obtain a gray difference of each grid in the vegetated area; S 604 : respectively obtaining a gray difference of a reference grid in the bare area and a gray difference of a reference grid in the vegetated area according to the gray difference of each grid in the bare area and the gray difference of each grid in the vegetated area, and calculating a soil erosion volume of a single gray value change according to the soil erosion volume of the reference grid in the bare area, the soil erosion volume of the reference grid in the vegetated area, the gray difference of the reference grid in the bare area and the gray difference of the reference grid in the vegetated area; wherein an expression for calculating the soil erosion volume of the single gray value change in the S 604 is as follows:

E
s

=

A

s
,
O

H

s
,
O

E
c

=

A

c
,
O

H

c
,
O

wherein E s is a soil erosion volume of a single gray value change in the bare area; A s,O is the soil erosion volume of the reference grid in the bare area; H s,O is the gray difference of the reference grid in the bare area; E c is a soil erosion volume of a single gray value change in the vegetated area; A c,O is the soil erosion volume of the reference grid in the vegetated area; and H c,O is the gray difference of the reference grid in the vegetated area;
S 605 : calculating the soil erosion volume corrected by soil erodibility according to the soil erosion volume of the single gray value change, the gray difference of each grid in the bare area and the gray difference of each grid in the vegetated area; wherein an expression for calculating the soil erosion volume corrected by soil erodibility in the S 605 is as follows:

A
=

∑

H

s
,
i

⁢

E
s

⁢

S
i
s

S
O
s

+

∑

H

c
,
j

⁢

E
c

⁢

S
j
c

S
O
c

wherein A is the soil erosion volume corrected by soil erodibility; H s,i is a gray difference of an i th grid in the bare area; E s is the soil erosion volume of the single gray value change in the bare area; S i s is soil erodibi

Assignees

• China Inst Water Resources & Hydropower Res

Inventors

• Liu Guangquan
• Tu Xiaoning
• Ai Ning
• Wang Xueyan
• Li Puhang
• Du Pengfei
• Qi Chunyu
• Zhang Bin
• Qi Yingyin
• Hou Jian

Classifications

• G06T2207/30188

  Vegetation; Agriculture · CPC title

• G06T2207/10032

  Satellite or aerial image; Remote sensing · CPC title

• G06V20/188

  Vegetation · CPC title

• G06V10/26

  Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion · CPC title

• G06T7/62

  of area, perimeter, diameter or volume · CPC title