Method for automatically judging incidence relation between water intaking user and river section

Abstract

The invention provides a method for automatically judging the incidence relation between a water intaking user and a river section, and relates to the technical field of water conservancy and water resources, and the method comprises the steps: obtaining digital elevation data covering a research area, and carrying out the preprocessing of the digital elevation data; the flow direction of the grid units is calculated, the grid units are traversed, the convergence range of the section is automatically generated, a unique identifier is distributed, and grid data of the convergence range of the section are generated; a coordinate system of the water intake data is converted to be consistent with the digital elevation data of the research area, a unique identifier is allocated to each water intake, the water intake point data is converted into grid unit data with the same resolution as the grid unit, and the grid unit is assigned according to the ID of the water intake; and according to the generated raster data of the section convergence range and the assignment result, judging the space association relationship between the water taking user position and the section convergence range. According to the invention, the water intaking users distributed in the section convergence range can be judged in batches.

Description

Technical Field

[0001] The invention belongs to the technical field of water conservancy and water resources, and in particular relates to a method for automatically distinguishing the correlation between a water user and a river section. Background Art

[0002] With the development of economy and society, the number of water users (such as enterprises, agricultural irrigation areas, water plants, etc.) who directly take water from natural water bodies has increased significantly. Under the background of refined management of water resources, the management of water users has been gradually strengthened. The locations of water users in the basin / region are relatively scattered, but their water intake, water use and drainage processes have a great impact on the runoff of the basin / regional control section. Therefore, it is very important to determine the number and amount of water users within the confluence range of the control section when constructing the basin / regional water resources allocation model, which is conducive to improving the accuracy of the model calculation of the section runoff. At the same time, when there is water shortage or insufficient ecological flow in the section, it is necessary to regulate the water users with related relationships in the river section in a targeted manner.

[0003] At present, there is no targeted method to automatically identify the spatial relationship between water users and key river sections. It is mainly based on GIS software to generate watershed partitions (but the section location cannot be considered) or to determine the confluence range of river sections based on topographic maps, and to determine the confluence range of the section where the water users are located based on manual naked eye.

[0004] The manual judgment method will cause certain errors. For example, it is difficult to distinguish the relationship between water users and watersheds around the watershed or in the plains, so the water users in watershed A are judged to be in watershed B. When there are many self-provided water users in a certain area or the relationship between river systems is complex, manual judgment will increase a lot of work time. In addition, the existing watershed generation method is basically automatically generated based on DEM data, and it is impossible to generate the confluence range of the section based on the section location alone. Summary of the invention

[0005] In view of the above-mentioned deficiencies in the prior art, the present invention provides an automatic determination method for the relationship between water users and river sections, which solves the problem that manual determination of the relationship between water users and sections is time-consuming, has large errors and is highly subjective.

[0006] In order to achieve the above purpose, the technical solution adopted by the present invention is: a method for automatically determining the association between water users and river sections, comprising the following steps:

[0007] S1, obtain digital elevation data covering the study area and preprocess it;

[0008] S2. Calculate the flow direction of the grid cell according to the preprocessed digital elevation data, traverse the grid cell, automatically generate the confluence range of the section, assign a unique identifier, and generate the confluence range grid data of the section;

[0009] S3, convert the coordinate system of the water intake data to be consistent with the digital elevation data of the study area, assign a unique identifier to each water intake, convert the water intake data to grid cell data with the same resolution as the grid cell, and assign values ​​to the grid cells according to the water intake ID;

[0010] S4. According to the generated cross-section runoff range grid data and the assignment result of step S3, the spatial correlation between the water user location and the cross-section runoff range is determined.

[0011] The beneficial effects of the present invention are as follows: the present application first generates the confluence range of the cross section based on the digital elevation data and cross section location of the study area; secondly, collects the data of water users and converts the coordinate information of the water users into raster data; finally, the raster data of the cross section confluence range and the raster data of the water user location are superimposed to identify their spatial correlation. The present invention realizes the rapid identification of the spatial correlation between water users and the confluence range of the watershed cross section, which can provide a technical basis for evaluating the impact of dispersed water users on cross section runoff, provide technical support for the construction of watershed / regional water resources allocation models and parameter calibration, and can provide assistance for the refinement and intelligent management of watershed / regional water resources.

[0012] Furthermore, the digital elevation data covering the study area is obtained and preprocessed, which is specifically as follows:

[0013] According to the longitude and latitude coordinates of the study area boundary, the digital elevation model is used to extract the digital elevation data covering the study area;

[0014] The extracted digital elevation data is preprocessed by removing outliers and filling missing values.

[0015] The beneficial effect of the above further scheme is that the present invention corrects the possible abnormal value and missing value problems in the original digital elevation data, which provides the basic working conditions for extracting the confluence range of the following section and avoids abnormal results.

[0016] Furthermore, the step S2 is specifically as follows:

[0017] According to the pre-processed digital elevation data, for any central grid, the slopes of the eight adjacent grids are calculated respectively;

[0018] Based on the slope, the convergence direction of the central grid is determined and the flow direction is coded;

[0019] Construct a two-dimensional relationship table of grid flow direction, and traverse all grid cells in turn, and determine the inflow and outflow relationship of adjacent grids according to the flow direction code. In the two-dimensional relationship table of grid flow direction, the vertical direction represents the outflow grid, and the horizontal direction represents the inflow grid;

[0020] The calculation starts from the downstream section, and searches upward for grids with hydraulic connections based on the two-dimensional flow direction relationship table. All grids with the same hydraulic connection are the confluence range of the section, and a unique identifier ID is assigned. BASIN , generate section confluence range grid data, where the sections with upstream and downstream relationships on the main stream are traversed to the sections, and the tributaries are traversed to the grids without inflow relationships. The value of the generated section confluence range grid is ID BASIN , representing the cross-section catchment range to which the grid belongs.

[0021] Furthermore, the slope is expressed as follows:

[0022]

[0023] Among them, S i Indicates the slope, Z c Indicates the elevation value of the center grid, Z i Represents the elevation value of the adjacent grid, d i Represents the distance between two grids.

[0024] Furthermore, the expression of the flow direction coding is as follows:

[0025]

[0026] Among them, DIR represents the flow direction coding, argmax(S i ) indicates the number of the surrounding adjacent grid with the largest slope.

[0027] The beneficial effect of the above further scheme is that the present invention determines the flow direction based on the elevation through specific flow direction calculation and coding rules, determines the hydraulic connection between grids, and provides a basis for determining the confluence range of the section.

[0028] Furthermore, the step S3 is specifically as follows:

[0029] The coordinate system of the water intake data is converted to be consistent with the digital elevation data of the study area, and a unique identifier ID is assigned to each water intake. INT ;

[0030] Initialize the water intake grid data so that its number of rows and columns and grid resolution are consistent with the cross-section confluence range grid data. The latitude range of the generated grid data is [lat j |j=1,2,…,n], the longitude value range is [lon i |i=1,2,…,m],lat j Indicates the latitude range of the raster data, lon iIndicates the range of longitude values, j indicates the grid row index in the latitude direction, from 1 to n, indicating the current row; n indicates the total number of grid rows in the latitude direction; i indicates the grid column index in the longitude direction, from 1 to m, indicating the current column; m indicates the total number of grid columns in the longitude direction;

[0031] Traverse all water users, and determine the grid position of the water user and the corresponding grid data row and column number according to the longitude and latitude of the water user; if the longitude and latitude of a water user is (lat0, lon0), the longitude and latitude distance of the water user from the grid point is:

[0032] DIS lat =Abs(lat0-[lat j |i=1,2,…,n])

[0033] DIS lon =Abs(lon0-[lon i |i=1,2,…,m]

[0034] The water user is closest to the j0th column of the grid and the i0th row of the grid, and the grid (i0, j0) in the grid data is determined as the grid where the water user is located. lat and DIS lon They represent the latitude and longitude distances of the water user from the grid point, Abs() represents the absolute value function, and j0 represents Min(DIS lat ) Index value, i0 means Min(DIS lon )'s Index value;

[0035] The unique identifier of the water intake ID INT Mapped to the grid cell (i0, j0) of the water intake, the grid cell is assigned ID INT Water user ID.

[0036] The beneficial effect of the above further scheme is that the present invention generates water intake grid data that is consistent with the cross-section confluence range grid data range, so as to perform matching calculations on the corresponding grids of the two files in step S4, and assign a unique water user ID to each grid to support the calculation in step S4.

[0037] Furthermore, the step S4 is specifically as follows:

[0038] The following formula is used to overlay and combine the cross-section confluence range grid data and the water intake user grid data to generate a new grid unit, where each pixel value in the new grid unit is a combination of the cross-section confluence range ID and the water intake ID:

[0039] R combined (i,j)=(R INT(i,j)<<4)+R BASIN (i,j)

[0040] Among them, R combined (i,j) represents the combined value of grid c(i,j), R BASIN (i,j) represents the grid value corresponding to the i-th row and j-th column in the cross-section confluence range grid data, and its value is the unique identifier ID of the corresponding cross-section confluence range BASIN , R INT (i,j) represents the grid value corresponding to the i-th row and j-th column in the water intake grid data, and its value is the unique identifier ID of the corresponding water intake. INT , R INT (i, j) << 4 means shifting the water intake ID 4 bits to the left;

[0041] The new grid cells are traversed, and the spatial correlation between the section runoff range and the water users is determined based on the newly generated pixel values. Based on the determination results, a data set of the spatial correlation between the water users and the section runoff range is generated.

[0042] The beneficial effect of the above further scheme is: the present invention generates new grid units based on the superposition and combination of section confluence range raster data and water intake raster data, and judges the spatial correlation between the section confluence range and the water user according to the newly generated pixel value, and can automatically generate the spatial correlation between the water user and the section confluence range. BRIEF DESCRIPTION OF THE DRAWINGS

[0043] Figure 1 The figure is a flow chart of the method of the present invention. DETAILED DESCRIPTION

[0044] The specific implementation modes of the present invention are described below so that those skilled in the art can understand the present invention. However, it should be clear that the present invention is not limited to the scope of the specific implementation modes. For those of ordinary skill in the art, as long as various changes are within the spirit and scope of the present invention as defined and determined by the attached claims, these changes are obvious, and all inventions and creations utilizing the concept of the present invention are protected.

[0045] Example

[0046] like Figure 1 As shown, the present invention provides a method for automatically determining the association between water users and river sections, and the implementation method is as follows:

[0047] S1. Obtain digital elevation data covering the study area and preprocess it, specifically:

[0048] According to the longitude and latitude coordinates of the study area boundary, the digital elevation model is used to extract the digital elevation data covering the study area;

[0049] The extracted digital elevation data is preprocessed by removing outliers and filling missing values.

[0050] In this embodiment, elevation data covering the study area is extracted from the digital elevation model (DEM) database based on the longitude and latitude coordinates of the study area boundary, and preprocessed, including removing outliers and filling missing data, to avoid errors caused by local "depressions" during flow direction analysis.

[0051] S2. Calculate the flow direction of the grid cells according to the preprocessed digital elevation data, traverse the grid cells, automatically generate the confluence range of the section, assign a unique identifier, and generate the confluence range grid data of the section, which is specifically:

[0052] According to the pre-processed digital elevation data, for any central grid, the slopes of the eight adjacent grids are calculated respectively;

[0053] Based on the slope, the convergence direction of the central grid is determined and the flow direction is coded;

[0054] Construct a two-dimensional relationship table of grid flow direction, and traverse all grid cells in turn, and determine the inflow and outflow relationship of adjacent grids according to the flow direction code. In the two-dimensional relationship table of grid flow direction, the vertical direction represents the outflow grid, and the horizontal direction represents the inflow grid;

[0055] The calculation starts from the downstream section, and searches upward for grids with hydraulic connections based on the two-dimensional flow direction relationship table. All grids with the same hydraulic connection are the confluence range of the section, and a unique identifier ID is assigned. BASIN , generate section confluence range grid data, where the sections with upstream and downstream relationships on the main stream are traversed to the sections, and the tributaries are traversed to the grids without inflow relationships. The value of the generated section confluence range grid is ID BASIN , representing the cross-section catchment range to which the grid belongs.

[0056] In this embodiment, the slope calculation is as follows: for the central grid c, the slopes of the eight surrounding adjacent grids are calculated respectively:

[0057]

[0058] In this embodiment, the flow direction is determined. Based on the slope, the confluence direction of the central grid c is determined, and the flow direction is coded:

[0059]

[0060] Among them, DIR represents the flow direction code, which corresponds to the surrounding adjacent grid numbers as 1, 2, 4, 8, 16, 32, 64, and 128, respectively. i) represents the number of the adjacent grid with the largest slope, S i Indicates the slope, Z c Indicates the elevation value of the center grid, Z i Represents the elevation value of the adjacent grid, d i Represents the distance between two grids.

[0061] In this embodiment, a grid flow direction relationship table is established: a two-dimensional grid flow direction relationship table is established, in which the vertical direction represents the outflow grid and the horizontal direction represents the inflow grid; all grid cells are traversed in sequence, and the inflow and outflow relationship of adjacent grids is determined according to the flow direction coding. If there is a confluence relationship, it is marked as 1, and if there is no confluence relationship, it is marked as 0.

[0062] In this embodiment, the confluence range of a section is determined by traversing and calculating from the downstream section, searching upward for grids with hydraulic connections based on the flow direction relationship table, wherein the sections with upstream and downstream relationships on the main stream are traversed to the sections; the tributaries are traversed to grids without inflow relationships. All grids with hydraulic connections on the same section are the confluence range of the section, and are assigned a unique ID, for example, the ID format is "BASIN_" plus four digits (such as BASIN_0001).

[0063] S3, convert the coordinate system of the water intake data to be consistent with the digital elevation, assign a unique identifier to each water intake, convert the water user data to grid cell data with the same resolution as the grid cell, and assign values ​​to the grid cells according to the water intake ID, which is specifically:

[0064] The coordinate system of the water intake data is converted to be consistent with the digital elevation data of the study area, and a unique identifier ID is assigned to each water intake. INT ;

[0065] Initialize the water intake grid data so that its number of rows and columns and grid resolution are consistent with the cross-section confluence range grid data. The latitude range of the generated grid data is [lat j |j=1,2,…,n], the longitude value range is [lon i |i=1,2,…,m],lat j Indicates the latitude range of the raster data, lon i Indicates the range of longitude values, j indicates the grid row index in the latitude direction, from 1 to n, indicating the current row; n indicates the total number of grid rows in the latitude direction; i indicates the grid column index in the longitude direction, from 1 to m, indicating the current column; m indicates the total number of grid columns in the longitude direction;

[0066] Traverse all water users, and determine the grid position of the water user and the corresponding grid data row and column number according to the longitude and latitude of the water user; if the longitude and latitude of a water user is (lat0, lon0), the longitude and latitude distance of the water user from the grid point is:

[0067] DIS lat =Abs(lat0-[lat j |i=1,2,…,n])

[0068] DIS lon =Abs(lon0-[lon i |i=1,2,…,m]

[0069] The water user is closest to the j0th column of the grid and the i0th row of the grid, and the grid (i0, j0) in the grid data is determined as the grid where the water user is located. lat and DIS lon They represent the latitude and longitude distances of the water user from the grid point, Abs() represents the absolute value function, and j0 represents Min(DIS lat ) Index value, i0 means Min(DIS lon )'s Index value;

[0070] The unique identifier of the water intake ID INT Mapped to the grid cell (i0, j0) of the water intake, the grid cell is assigned ID INT Water user ID.

[0071] In this embodiment, determine the water intake identification number: determine the basic information of the water intake such as the name, longitude, and latitude, and assign a unique ID to it. The ID format is "INT_" plus four digits (such as INT_0001)

[0072] In this embodiment, the conversion into raster data is as follows: the water intake raster data is initialized so that the number of rows and columns and the grid resolution are consistent with the cross-sectional confluence range unit raster data. All water users are traversed and mapped to the cross-sectional confluence range grid cells according to the number of rows and columns, and the grid cells are assigned to the water user ID. That is, the data element encoding of the water user is performed; an empty raster data (with longitude and latitude information) is generated, called the water intake raster; based on the location information, the water intake is traversed and calculated, and based on the attribute information, the raster is matched to add data attributes. Among them, the matching process includes: 1. Calculating the spatial distance between the longitude and latitude of the water user and the longitude and latitude of all grid cells; 2. Assigning the water intake information to the grid cell with the minimum distance.

[0073] S4. According to the generated cross-section confluence range grid data and the assignment result of step S3, the spatial correlation between the water user location and the cross-section confluence range is determined. The implementation method is as follows:

[0074] The following formula is used to overlay and combine the cross-section confluence range grid data and the water intake user grid data to generate a new grid unit, where each pixel value in the new grid unit is a combination of the cross-section confluence range ID and the water intake ID:

[0075] R combined (i,j)=(R INT (i,j)<<4)+R BASIN (i,j)

[0076] Among them, R combined (i,j) represents the combined value of grid c(i,j), R BASIN (i,j) represents the grid value corresponding to the i-th row and j-th column in the cross-section confluence range grid data, and its value is the unique identifier ID of the corresponding cross-section confluence range BASIN , R INT (i,j) represents the grid value corresponding to the i-th row and j-th column in the water intake grid data, and its value is the unique identifier ID of the corresponding water intake. INT , R INT (i, j) << 4 means shifting the water intake ID 4 bits to the left;

[0077] The new grid cells are traversed, and the correlation between the water users and the cross-section confluence range is judged according to the newly generated pixel values. Based on the judgment results, a data set of the spatial correlation between the water users and the cross-section confluence range is generated.

[0078] In this embodiment, the water user is superimposed with the section confluence range: the bitwise operation algorithm is used to superimpose and combine the section confluence range unit grid data and the water user grid data to generate a new grid data, in which each pixel value is a combination of the section confluence range ID and the water intake ID.

[0079] In this embodiment, the newly generated grid cells are traversed, and the spatial association between the cross-section flow range and the water users is determined based on the newly generated pixel values. For example, the pixel value 00010010 means that the 0010 water intake is located within the 0001 cross-section flow range.

[0080] In this embodiment, a spatial correlation relationship dataset between water users and cross-section confluence ranges is generated, including the correspondence between water intake IDs and cross-section IDs, water intake information, cross-section information, and other contents.

[0081] In summary, the present invention realizes the rapid identification of the spatial correlation between water users and the confluence range of river basin sections, which can provide a technical basis for evaluating the impact of dispersed water users on cross-sectional runoff, provide technical support for the construction of river basin / regional water resources allocation models and parameter calibration, and provide assistance for the refined and intelligent management of river basin / regional water resources.

Claims

1. A method for automatically determining the relationship between water users and river sections, characterized in that: The following steps are involved: S1, obtain digital elevation data covering the study area and preprocess it; S2. Calculate the flow direction of the grid cell according to the preprocessed digital elevation data, traverse the grid cell, automatically generate the confluence range of the section, assign a unique identifier, and generate the confluence range grid data of the section; S3, convert the coordinate system of the water intake data to be consistent with the digital elevation data of the study area, assign a unique identifier to each water intake, convert the water intake data to grid cell data with the same resolution as the grid cell, and assign values ​​to the grid cells according to the water intake ID; S4. According to the generated cross-section runoff range grid data and the assignment result of step S3, the spatial correlation between the water user location and the cross-section runoff range is determined.

2. The method for automatically determining the relationship between water users and river sections according to claim 1, characterized in that: The digital elevation data covering the study area is obtained and preprocessed, which is specifically as follows: According to the longitude and latitude coordinates of the study area boundary, the digital elevation model is used to extract the digital elevation data covering the study area; The extracted digital elevation data is preprocessed by removing outliers and filling missing values.

3. The method for automatically determining the relationship between water users and river sections according to claim 1, characterized in that: The step S2 is specifically as follows: According to the pre-processed digital elevation data, for any central grid, the slopes of the eight adjacent grids are calculated respectively; Based on the slope, the convergence direction of the central grid is determined and the flow direction is coded; Construct a two-dimensional relationship table of grid flow direction, and traverse all grid cells in turn, and determine the inflow and outflow relationship of adjacent grids according to the flow direction code. In the two-dimensional relationship table of grid flow direction, the vertical direction represents the outflow grid, and the horizontal direction represents the inflow grid; The calculation starts from the downstream section, and searches upward for grids with hydraulic connections based on the two-dimensional relationship table of flow direction. All grids with the same hydraulic connection are the confluence range of the section, and a unique identifier ID is assigned. BASIN , generate section confluence range grid data, where the sections with upstream and downstream relationships on the main stream are traversed to the sections, and the tributaries are traversed to the grids without inflow relationships. The value of the generated section confluence range grid is ID BASIN , representing the cross-section catchment range to which the grid belongs.

4. The method for automatically determining the relationship between water users and river sections according to claim 3, characterized in that: The slope expression is as follows: Among them, S i Indicates the slope, Z c Indicates the elevation value of the center grid, Z i Represents the elevation value of the adjacent grid, d i Represents the distance between two grids.

5. The method for automatically determining the relationship between water users and river sections according to claim 4, characterized in that: The expression of the flow direction coding is as follows: Among them, DIR represents the flow direction coding, argmax(S i ) indicates the number of the surrounding adjacent grid with the largest slope.

6. The method for automatically determining the relationship between water users and river sections according to claim 1, characterized in that: The step S3 is specifically as follows: The coordinate system of the water intake data is converted to be consistent with the digital elevation data of the study area, and a unique identifier ID is assigned to each water intake. INT ; Initialize the water intake grid data so that its number of rows and columns and grid resolution are consistent with the cross-section confluence range grid data. The latitude range of the generated grid data is [lat j |j=1,2,…,n], the longitude value range is [lon i |i=1,2,…,m],lat j Indicates the latitude range of the raster data, lon i Indicates the range of longitude values, j indicates the grid row index in the latitude direction, from 1 to n, indicating the current row; n indicates the total number of grid rows in the latitude direction; i indicates the grid column index in the longitude direction, from 1 to m, indicating the current column; m indicates the total number of grid columns in the longitude direction; Traverse all water users, and determine the grid position of the water user and the corresponding grid data row and column number according to the longitude and latitude of the water user; if the longitude and latitude of a water user is (lat0, lon0), the longitude and latitude distance of the water user from the grid point is: DIS lat =Abs(lat0-[lat j |i=1,2,…,n]) DIS lon =Abs(lon0-[lon i |i=1,2,…,m] The water user is closest to the j0th column of the grid and the i0th row of the grid, and the grid cell (i0, j0) in the grid data is determined to be the grid where the water user is located. lat and DIS lon They represent the latitude and longitude distances of the water user from the grid point, Abs() represents the absolute value function, and j0 represents Min(DIS lat ) Index value, i0 means Min(DIS lon )'s Index value; The unique identifier of the water intake ID INT Mapped to the grid cell (i0, j0) of the water intake, where the grid cell (i0, j0) is assigned ID INT Water user ID.

7. The method for automatically determining the relationship between water users and river sections according to claim 1, characterized in that: The step S4 is specifically as follows: The following formula is used to overlay and combine the cross-section confluence range grid data and the water intake user grid data to generate a new grid unit, where each pixel value in the new grid unit is a combination of the cross-section confluence range ID and the water intake ID: R combined (i,j)=(R INT (i,j)<<4)+R BASIN (i,j) Among them, R combined (i,j) represents the combined value of grid c(i,j), R BASIN (i,j) represents the grid value corresponding to the i-th row and j-th column in the cross-section confluence range grid data, and its value is the unique identifier ID of the corresponding cross-section confluence range BASIN , R INT (i,j) represents the grid value corresponding to the i-th row and j-th column in the water intake grid data, and its value is the unique identifier ID of the corresponding water intake INT , R INT (i, j)<<4 means shifting the water intake ID left by 4 bits; The new grid cells are traversed, and the spatial correlation between the section runoff range and the water users is determined based on the newly generated pixel values. Based on the determination results, a data set of the spatial correlation between the water users and the section runoff range is generated.

Images