A method for interpreting the topological relationship of water resource zoning codes with end-to-end connection

Abstract

The present invention provides a method for interpreting the topological relationship of water resource zoning codes in a head-to-tail manner. The method comprises: verifying the water resource zoning code structure, parsing the code expression, and counting the number of water resource zoning at each level to construct a multidimensional vector matrix for hierarchical storage; interpreting the topological relationship of water resource zoning codes at the same level, and based on the multidimensional vector matrix hierarchical storage combined with the water resource zoning code structure, using the lowest level water resource zoning to jump out step by step in a head-to-tail manner; connecting the water resource zoning at each level in series based on the head-to-tail manner to form a topological relationship interpretation code that increases sequentially from upstream to downstream, and using the confluence cumulative number method to count the size of the confluence cumulative number in the topological relationship interpretation code, and comparing the codes after the topological relationship interpretation. The present invention greatly simplifies the repetitive work in the hydrological model construction process and solves the problem of inconsistent boundaries between the statistical results of water resource zoning in regional water resource surveys and assessments and the simulation results of model processing.

Description

Technical Field

[0001] The present invention belongs to the technical field of hydrology and water resources, and in particular relates to a method for interpreting the topological relationship of end-to-end water resource zoning codes. Background Art

[0002] Distributed hydrological models are currently gaining popularity in hydrology and water resources research due to their consideration of basin-specific physical processes and their ability to apply large-scale spatial product data. They have significantly improved their feasibility in data-scarce regions. The upstream and downstream topological relationships between hydrological response units within a basin, a key component of distributed hydrological models, influence the flow patterns of rivers in these models, which in turn influences runoff and confluence characteristics, ultimately determining the accuracy of the model's simulated data. Currently, water resource assessments often utilize water resource zoning, which is uniquely coded. These assessments often focus on macroscopic, large-scale, and regularity-based summaries, complementing the detailed conclusions of distributed hydrological models. The zoning approach is similar to the division of hydrological response units in distributed hydrological models.

[0003] Therefore, when evaluating relevant watersheds, it is necessary not only to divide the hydrological response units within the distributed hydrological model but also to zonal code the water resources. When coding the topological relationships of the distributed hydrological model, it is necessary to again use raster files such as flow direction and flow rate to perform consecutive numbering to further explore the upstream and downstream confluence relationships. This is highly repetitive with the work done during water resource zoning coding, which is time-consuming and prone to problems such as inconsistent evaluation units. Code interpretation serves as a bridge to further connect the coding of watershed water resources regions with the coding of the topological relationships of the distributed hydrological model. To this end, a method for interpreting the topological relationships of water resource zoning codes has been developed. Summary of the Invention

[0004] In response to the above-mentioned deficiencies in the prior art, the present invention provides a method for interpreting the topological relationship of water resource zoning codes that is connected end to end. This method greatly reduces the time spent on repetitive work in large-scale basin or regional water resource evaluation and distributed hydrological model construction, and is conducive to improving the consistency of regional water resource zoning evaluation and distributed hydrological model evaluation scope.

[0005] In order to achieve the above purpose, the technical solution adopted by the present invention is:

[0006] This solution provides a method for interpreting the topological relationship of water resource zoning codes in a coherent manner, including the following steps:

[0007] S1. By verifying the water resource partition coding structure, parsing the coding expression, and counting the number of water resource partitions at each level, a multi-dimensional vector matrix is ​​constructed for hierarchical storage;

[0008] S2. Interpret the topological relationship of water resource zoning codes at the same level, and based on the hierarchical storage of multi-dimensional vector matrices combined with the water resource zoning code structure, use the lowest level water resource zoning end-to-end connection method to jump out step by step;

[0009] S3. Connect the water resource zones at all levels in series based on the head-to-tail connection method to form a topological relationship interpretation code that increases sequentially from upstream to downstream, and use the confluence cumulative number method to count the size of the confluence cumulative number in the topological relationship interpretation code, compare the codes after the interpretation of the topological relationship, determine the continuity and accuracy of the topological relationship, and complete the interpretation of the topological relationship of the water resource zone code.

[0010] The beneficial effects of the present invention are as follows: the present invention interprets the topological relationship code by connecting the water resource zoning code end to end, so that the water resource zoning evaluation unit is consistent with the scope of the hydrological response unit of the distributed hydrological model, which is helpful to combine the macro-regularity conclusion with the detailed hydrological simulation process; the end-to-end connection method reduces the process of performing conventional hydrological analysis again as much as possible, greatly reduces the time spent on repetitive work for large-scale basins or regions, and significantly improves the efficiency of hydrological work.

[0011] Furthermore, the step S1 includes the following steps:

[0012] S101. Verify the coding structure of water resource zoning and analyze the expressions for upstream and downstream, main and tributary rivers, and different levels of water resource zoning in the coding;

[0013] S102. If the water resource partition code does not have an upstream and downstream relationship or does not distinguish between different levels of water resource partitions, pre-process the data topology relationship of the resource partition code;

[0014] S103. Count the number of water resource partitions at each level and construct a multi-dimensional vector matrix for hierarchical storage.

[0015] The beneficial effect of the above further scheme is that the encoding method in the present invention is arranged according to a fixed logic. By parsing the encoding expression method and storing it in a hierarchical manner, the encoding can be effectively split. The hierarchical and block storage format is conducive to faster recognition and reading by the computer.

[0016] Furthermore, step S2 includes the following steps:

[0017] S201, using level n to represent the lowest level in the water resource zoning code, and constructing a topological relationship of the lowest level water resource zoning within all n-1 levels in the basin based on the parsing of the parsing code expression, where n represents the level of the water resource zoning;

[0018] S202, according to the hierarchical storage of the multidimensional vector matrix, locate the data positions of the most upstream and the most downstream in the multidimensional vector matrix step by step, and define the downstream of the water resource partition into 0;

[0019] S203, judging whether the n-1 level water resource partition where each n-level outlet is located is the outlet of the basin according to the difference in the parsing coding method; or

[0020] According to the difference in the parsing coding method, it is determined whether the coding area corresponds to the main stream;

[0021] S204. Based on the judgment result, combined with the river structure of the basin and the expression of the analytical code, the data position of the most upstream water resource partition where the n-1-level mainstream stream is located is imported. If there is a skipped level in the coding method, the data position of the most upstream water resource partition where the skipped level mainstream stream is located is imported;

[0022] S205. Based on the outlet of the n-th level water resource zone, connect the topological relationship of the water resource zones in the same level step by step and connect them in series. The downstream of the water resource zone is used as the water resource zone entering the sea as the end condition, and the topological relationship of the entire watershed is interpreted to complete the step-by-step exit.

[0023] The beneficial effect of the above-mentioned further scheme is that the present invention can directly interpret the topological relationship according to the digital changes corresponding to the construction characteristics of the coded topological relationship in a head-to-tail manner, without further relying on the flow direction, flow rate and other contents in conventional hydrological analysis to obtain the spatial position of each water resource partition through comparative judgment. This method is significantly faster than the latter in large-scale areas; and many conventional hydrological analysis tasks need to be performed when coding water resources. If this work still needs to be repeated again when coding the topological relationship of water resource partitions, it will undoubtedly increase the workload. This method significantly improves the efficiency of hydrological work and reduces repetitive work as much as possible.

[0024] Furthermore, the lowest level water resource zones are connected step by step, and the following conditions must be met:

[0025] The first condition: when multiple water resources flow into the sea, the location of the water resource divisions at the estuary should be clearly explored;

[0026] The second condition: whether there is a cross-level when the lowest-level water resource zone corresponding to the outlet jumps to the n-1 level, and combined with the river structure of the basin, determine the main stream water resource zone code connected when the lower-level water resource zone jumps in.

[0027] The beneficial effect of the above further scheme is: the present invention fully considers the advantages of block storage in step S1, and through the conditions of the water outlet in the block data, the topological relationship of the entire water resource partition can be connected in series; while jumping upward from the lowest level step by step to avoid omissions, it reduces repetitive work. BRIEF DESCRIPTION OF THE DRAWINGS

[0028] Figure 1Flow chart of the method of the present invention.

[0029] Figure 2 The European Danube region and its water resource zoning coding interpretation map used in the embodiments of the present invention.

[0030] Figure 3 This is a diagram showing the topological coding interpretation results in part of the study area of ​​the embodiment of the present invention.

[0031] Figure 4 This is a topological relationship diagram for topological interpretation in an embodiment of the present invention. DETAILED DESCRIPTION

[0032] The specific embodiments of the present invention are described below to facilitate understanding of the present invention by those skilled in the art. However, it should be clear that the present invention is not limited to the scope of the specific embodiments. For those skilled in the art, as long as various changes are within the spirit and scope of the present invention as defined and determined by the appended claims, these changes are obvious, and all inventions and creations utilizing the concepts of the present invention are protected.

[0033] Example

[0034] like Figure 1 As shown, the present invention provides a method for interpreting the topological relationship of water resource zoning codes, which is implemented as follows:

[0035] S1. By verifying the water resource partition coding structure, parsing the coding expression, and counting the number of water resource partitions at each level, a multi-dimensional vector matrix hierarchical storage is constructed. The implementation method is as follows:

[0036] S101. Verify the coding structure of water resource zoning and analyze the expressions for upstream and downstream, main and tributary rivers, and different levels of water resource zoning in the coding;

[0037] S102. If the water resource partition code does not have an upstream and downstream relationship or does not distinguish between different levels of water resource partitions, pre-process the data topology relationship of the resource partition code;

[0038] S103. Count the number of water resource partitions at each level and construct a multi-dimensional vector matrix for hierarchical storage.

[0039] In this embodiment, the water resource partition coding structure is verified, and the main method of distinguishing upstream and downstream (distinction methods such as partition codes gradually increasing or decreasing from upstream to downstream), main and tributary streams (distinction methods such as introducing new characters, odd and even numbers), and different levels of water resource partitions (distinction methods such as setting different lengths) in the coding expression is analyzed; if the water resource partition code does not have an upstream and downstream relationship or the water resource partitions of different levels have not yet been distinguished, the water resource partition code can be preprocessed as needed; the number of water resource partitions at each level is further counted, and a multi-dimensional vector matrix is ​​constructed for hierarchical storage.

[0040] S2. Interpret the topological relationship of water resource zoning codes at the same level, and based on the hierarchical storage of multi-dimensional vector matrices combined with the water resource zoning code structure, use the lowest level water resource zoning end-to-end connection method to jump out step by step. The implementation method is as follows:

[0041] S201, using level n to represent the lowest level in the water resource zoning code, and constructing a topological relationship of the lowest level water resource zoning within all n-1 levels in the basin based on the parsing of the parsing code expression, where n represents the level of the water resource zoning;

[0042] S202, according to the hierarchical storage of the multidimensional vector matrix, locate the data positions of the most upstream and the most downstream in the multidimensional vector matrix step by step, and define the downstream of the water resource partition into 0;

[0043] S203, judging whether the n-1 level water resource partition where each n-level outlet is located is the outlet of the basin according to the difference in the parsing coding method; or

[0044] According to the difference in the parsing coding method, it is determined whether the coding area corresponds to the main stream;

[0045] S204. Based on the judgment result, combined with the river structure of the basin and the expression of the analytical code, the data position of the most upstream water resource partition where the n-1-level mainstream stream is located is imported. If there is a skipped level in the coding method, the data position of the most upstream water resource partition where the skipped level mainstream stream is located is imported;

[0046] S205. Based on the outlet of the n-th level water resource zone, connect the topological relationship of the water resource zones in the same level step by step and connect them in series. The downstream of the water resource zone is used as the water resource zone entering the sea as the end condition, and the topological relationship of the entire basin is interpreted to complete the step-by-step exit.

[0047] The lowest level water resource zones are connected end to end in a step-by-step manner, and the following conditions must be met:

[0048] The first condition: when multiple water resources flow into the sea, the location of the water resource divisions at the estuary should be clearly explored;

[0049] The second condition: whether there is a cross-level when the lowest-level water resource zone corresponding to the outlet jumps to the n-1 level, and combined with the river structure of the basin, determine the main stream water resource zone code connected when the lower-level water resource zone jumps in.

[0050] In this embodiment, the topological relationship between water resource partitions of the same level is interpreted according to the preprocessing in step S1; then, based on the multidimensional vector matrix constructed by S1 and combined with the water resource partition coding structure, the lowest level water resource partition is used to jump out step by step in a head-to-tail manner; when jumping out in a head-to-tail manner, the following three aspects are mainly considered: (1) when multiple water resource areas in a region flow into the sea, the location of the water resource area at the estuary is explored (reflected in the water resource partition, mainly the breakpoint of the water resource partition coding); (2) whether there is a cross-level (jump to level n-2, n-3..., 1) when the lowest level water resource area corresponding to the outlet jumps out to level n-1; (3) combined with the river basin structure, the main stream water resource partition connected when the lower level water resource partition jumps in is identified.

[0051] In this embodiment, level n is used to represent the lowest level in the water resource partition coding, and level 1 represents the highest level; based on the analysis of the coding expression in S1, the topological relationship of the lowest level (level n) water resource partitions in all n-1 levels in the basin is constructed; based on the constructed multidimensional vector matrix, the data positions of the upstream and downstream in the step-by-step matrix are located, and the downstream of the water resource partition entering the sea is further defined as 0; based on the differences in the coding methods of upstream and downstream, main and tributary streams, etc. in step S1, it is determined whether the n-1 level water resource partition where each n-level outlet is located is the outlet of the basin or the coding partition in the basin. Whether the area corresponds to the main stream or tributary, it is further combined with the river structure of the basin and the coding expression method, and the data position of the upstream of the partition where the n-1 level main stream is located is merged. If there is a skipping situation in the coding method (the n-level partition is directly merged into the n-2,..., 1 level), the data position of the upstream of the partition where the skipped level main stream is located is merged; finally, based on the outlet of the n-th level water resource partition, the topological relationship of the water resource partitions within the level is connected step by step, and the downstream of the water resource partition is used as the sea water resource partition (the downstream matrix data is 0) as the end condition to realize the interpretation of the topological relationship of the entire basin.

[0052] S3. Connect the water resource zones at all levels in series based on the head-to-tail connection method to form a topological relationship interpretation code that increases sequentially from upstream to downstream, and use the confluence cumulative number method to count the size of the confluence cumulative number in the topological relationship interpretation code, compare the codes after the interpretation of the topological relationship, determine the continuity and accuracy of the topological relationship, and complete the interpretation of the topological relationship of the water resource zone code.

[0053] In this embodiment, the water resource zones at all levels are connected in series based on the end-to-end connection in step S2 to form a topological relationship interpretation code that increases sequentially from upstream to downstream; the confluence accumulation number in the topological code is counted by the confluence accumulation number method in conventional hydrological analysis, and the codes after the topological relationship interpretation are compared to ensure the continuity and accuracy of the topological relationship.

[0054] In this example, the Danube River in Europe is selected as the research area. Figure 2 As shown in the figure, the digital elevation data used is downloaded from the 90-meter SRTM (Shuttle Radar Topography Mission) data provided by NASA. The conventional hydrological analysis and corresponding preprocessing process are not described here.

[0055] A1. Based on conventional hydrological analysis, the Danube River basin water resources are zoned according to the relationship between the main and tributary streams, gradually increasing from upstream to downstream. (In this embodiment, the water resources are mainly divided into five levels, with odd-numbered levels representing catchment areas and even-numbered levels representing tributaries.) This preprocessing process will not be described in detail. According to the coding expression of the water resource zones, the number of water resource zones at each level is further counted (for example, the first level contains five second-level zones, the first second-level zone contains 13 third-level zones, the first third-level zone contains five fourth-level zones, and the first fourth-level zone contains one fifth-level zone). Statistics are then performed level by level in this manner to further construct a multidimensional vector matrix for hierarchical storage.

[0056] A2. At this point, the topological relationship between the five-level areas within each four-level area is clear. In this example, the topological relationship between the five-level areas within each four-level area is arranged in the same order of 1, 2, ...n. The main problem is that the start and end codes of each five-level area are unclear. Therefore, combining the characteristics of the beginning and end and the coding structure of the water resource zoning, the lowest level water resource zoning is connected to the end and jumped up step by step. The implementation method is as follows:

[0057] First, based on the multidimensional vector matrix constructed by A1, the downstream of the water resource partition (No. 8828) is defined as 0;

[0058] Secondly, based on the analysis of the coding methods of upstream and downstream, main and tributary rivers in A1, the topological coding of the 5-level water resource partitions except the outlet is determined; further, the situation of the water resource partition corresponding to the 5-level outlet merging into the previous level is determined one by one, that is, merging into the partition at the upstream of the main stream of the 4-level or even 3-level 1 water resource area.

[0059] Finally, starting from the uppermost reaches of the basin, the exit points are gradually linked to the uppermost partition points of the main stream. After linking, the codes corresponding to the upstream outlets can be superimposed according to the 1-n sequential codes of the five-level zones within the four-level zones, and then the water resource zones with upstream and downstream characteristics within the five-level zones are obtained in series. The downstream of the water resource zone is taken as the water resource zone entering the sea (the downstream matrix data is 0) as the end condition to realize the interpretation of the topological relationship of the entire basin.

[0060] A3: After connecting the water resource zones at all levels in series based on the end-to-end connection in A2, a topological relationship interpretation code is formed that increases in sequence from upstream to downstream (such as Figure 2 and Figure 3 The cumulative number of confluence in the topological code was calculated by the cumulative number of confluence in conventional hydrological analysis. The codes after the topological relationship interpretation all showed a trend of continuous rising of the main river channel, further ensuring the continuity and accuracy of the topological relationship (as shown in Figure 2). Figure 4 shown).

[0061] The present invention interprets the topological relationship code by connecting the water resource zoning codes end to end, so that the water resource zoning evaluation unit is consistent with the scope of the hydrological response unit of the distributed hydrological model, which is conducive to combining the macro-regularity conclusion with the detailed hydrological simulation process; the end-to-end connection method minimizes the process of performing conventional hydrological analysis again, greatly reduces the time spent on repetitive work for large-scale basins or regions, and significantly improves the efficiency of hydrological work.

Claims

1. A method for interpreting the topological relationship of water resource zoning codes that is connected end to end, characterized in that: The following steps are involved: S1. By verifying the water resource partition coding structure, parsing the coding expression, and counting the number of water resource partitions at each level, a multi-dimensional vector matrix is ​​constructed for hierarchical storage; S2. Interpret the topological relationship of water resource zoning codes at the same level, and based on the hierarchical storage of multi-dimensional vector matrices combined with the water resource zoning code structure, use the lowest level water resource zoning end-to-end connection method to jump out step by step; The step S2 comprises the following steps: S201, using level n to represent the lowest level in the water resource zoning code, and constructing a topological relationship of the lowest level water resource zoning within all n-1 levels in the basin based on the parsing of the parsing code expression, where n represents the level of the water resource zoning; S202, according to the hierarchical storage of the multidimensional vector matrix, locate the data positions of the most upstream and the most downstream in the multidimensional vector matrix step by step, and define the downstream of the water resource partition into 0; S203, judging whether the n-1 level water resource partition where each n-level outlet is located is the outlet of the basin according to the difference in the parsing coding method; or According to the difference in the parsing coding method, it is determined whether the coding area corresponds to the main stream; S204. Based on the judgment result, combined with the river structure of the basin and the expression of the analytical code, the data position of the most upstream water resource partition where the n-1-level mainstream stream is located is imported. If there is a skipped level in the coding method, the data position of the most upstream water resource partition where the skipped level mainstream stream is located is imported; S205, based on the outlet of the n-th level water resource zone, connect the topological relationships of the water resource zones in the same level step by step, and use the downstream of the water resource zone as the water resource zone entering the sea as the end condition to obtain the interpretation of the topological relationship of the entire watershed, and complete the step-by-step exit; The lowest level water resource zones are connected end to end in a step-by-step manner, and the following conditions must be met: The first condition: when multiple water resources flow into the sea, the location of the water resource divisions at the estuary should be clearly explored; The second condition is whether there is a cross-level when the lowest-level water resource zone corresponding to the outlet jumps out to the n-1 level, and combined with the river structure of the basin, the code of the main stream water resource zone connected when the lower-level water resource zone jumps in is determined; S3. Connect the water resource zones at all levels in series based on the head-to-tail connection method to form a topological relationship interpretation code that increases sequentially from upstream to downstream, and use the confluence cumulative number method to count the size of the confluence cumulative number in the topological relationship interpretation code, compare the codes after the interpretation of the topological relationship, determine the continuity and accuracy of the topological relationship, and complete the interpretation of the topological relationship of the water resource zone code.

2. The method for interpreting the topological relationship of water resource zoning codes according to claim 1 is characterized in that: The step S1 comprises the following steps: S101. Verify the coding structure of water resource zoning and analyze the expressions for upstream and downstream, main and tributary rivers, and different levels of water resource zoning in the coding; S102. If the water resource partition code does not have an upstream and downstream relationship or does not distinguish between different levels of water resource partitions, pre-process the data topology relationship of the resource partition code; S103. Count the number of water resource partitions at each level and construct a multi-dimensional vector matrix for hierarchical storage.

Images