Water supply network subarea metering area division method and device and electronic equipment

Abstract

The invention relates to the technical field of municipal pipe networks, in particular to a water supply pipe network zoning metering area dividing method and device and electronic equipment. The method includes: acquiring digital elevation model data of a target area; performing depression filling on the digital elevation model data by adopting a mathematical morphology method to generate a depression-free elevation model; performing watershed division on the depression-free elevation model based on a hydrological analysis algorithm to obtain river network distribution characteristics and sub-basin division results; combining the river network distribution characteristics and the sub-basin boundaries to generate a preliminary partition metering region; performing spatial overlay analysis on the sub-basin division result and geographic information data to obtain water supply network management parameters; and adjusting the preliminary zoning metering area according to the water supply network management parameters, and outputting a final zoning result. The reliability of the division result is high, and support is provided for construction of urban water supply networks.

Description

Technical Field

[0001] The present invention relates to the technical field of municipal pipe networks, and particularly to a method, device and electronic equipment for dividing metering areas of a water supply pipe network. Background Art

[0002] DMA (District Metered Area) zonal metering plays an important role in the management of water supply systems. DMA zonal metering can accurately monitor the water supply and consumption in each zone. On the one hand, it can help detect leakage phenomena by comparing water supply and consumption, reducing water resource waste. On the other hand, it helps formulate reasonable water price policies to promote the sustainable utilization of water resources. Through DMA zoning, managers can conduct refined management of the water supply system, analyze the data of each zone, which can help identify abnormal situations, optimize water supply strategies, and improve the operation efficiency of the water supply system. In addition, the detailed data of DAM zoning can provide a scientific basis for the planning and expansion of urban water supply systems and rationally plan the water supply network.

[0003] There are various ways to divide metering zones. According to different objectives and methods, they can be divided into methods based on geographical regions, methods based on water use types, methods based on pipe network structures, methods based on pressure zones, methods based on water consumption, and methods that combine multiple factors. Among them, the division based on geographical regions, as the name implies, is to divide by combining geographical information such as administrative divisions. This method has a good guiding role in the project pre - plan.

[0004] However, the existing division method based on geographical information systems only considers factors such as administrative divisions and ignores factors such as terrain distribution, resulting in unreliable results of directly dividing DMA zones relying on geographical administrative boundaries, and a large amount of workload for subsequent changes. Summary of the Invention

[0005] In view of this, the embodiments of the present application provide a method, device and electronic equipment for dividing metering areas of a water supply pipe network to solve the problems of low reliability and large subsequent change workload of using traditional methods for dividing metering areas.

[0006] The first aspect of the embodiments of the present application provides a method for dividing metering areas of a water supply pipe network, including:

[0007] Obtain digital elevation model data of a target area;

[0008] Use the mathematical morphology method to fill depressions in the digital elevation model data to generate a depression - free elevation model;

[0009] Based on a hydrological analysis algorithm, perform watershed division on the depression - free elevation model to obtain river network distribution characteristics and sub - basin division results;

[0010] Generate a preliminary zoned metering area by combining the described river network distribution characteristics and sub - basin boundaries;

[0011] Perform a spatial overlay analysis on the sub - basin division result and the geographic information data to obtain water supply network management parameters;

[0012] Adjust the preliminary zoned metering area according to the water supply network management parameters and output the final zoning result.

[0013] The second aspect of the embodiments of the present application provides a device for dividing a zoned metering area of a water supply network, including:

[0014] A data acquisition module for acquiring digital elevation model data of a target area;

[0015] A reconstruction module for filling depressions in the digital elevation model data using mathematical morphology methods to generate a depression - free elevation model;

[0016] A watershed division module for performing watershed division on the depression - free elevation model based on a hydrological analysis algorithm to obtain river network distribution characteristics and sub - basin division results;

[0017] A preliminary zoning module for generating a preliminary zoned metering area by combining the river network distribution characteristics and sub - basin boundaries;

[0018] An overlay analysis module for performing a spatial overlay analysis on the sub - basin division result and the geographic information data to obtain water supply network management parameters;

[0019] An adjustment and output module for adjusting the preliminary zoned metering area according to the water supply network management parameters and outputting the final zoning result.

[0020] The third aspect of the embodiments of the present application provides an electronic device, including a processor, a memory, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the electronic device implements the method for dividing a zoned metering area of a water supply network provided in the first aspect of the embodiments of the present application.

[0021] The fourth aspect of the embodiments of the present application provides a computer program product, including a computer program. When the computer program is run, the method described in the first aspect of the embodiments of the present application is executed.

[0022] The method for dividing the water supply pipe network into partition metering areas provided by the first aspect of the embodiments of the present application includes obtaining the digital elevation model data of the target area; using the mathematical morphology method to fill the depressions in the digital elevation model data to generate a depression-free elevation model; based on the hydrological analysis algorithm, performing watershed division on the depression-free elevation model to obtain the river network distribution characteristics and the sub-watershed division results; combining the river network distribution characteristics and the sub-watershed boundaries to generate preliminary partition metering areas; performing spatial overlay analysis on the sub-watershed division results and the geographic information data to obtain the water supply pipe network management parameters; and adjusting the preliminary partition metering areas according to the water supply pipe network management parameters to output the final partition results. This method accurately eliminates the false depressions in the DEM data through mathematical morphology, avoiding subsequent hydrological analysis errors; uses the hydrological analysis method to divide the sub-watersheds of the DEM after geodesic reconstruction, generates the initial partitions based on the natural terrain features (river network, sub-watersheds), improving the scientificity and reliability of the partition scheme; fuses the geographic information data to realize the quantification of management parameters, ensuring that the partition results meet the actual management requirements; the full-process automatic processing significantly improves the efficiency of the water supply pipe network planning, is more efficient than the traditional manual division method, and no subsequent changes are required, providing support for the construction of the urban water supply network.

[0023] It can be understood that the beneficial effects of the second to fourth aspects above can be referred to the relevant descriptions in the first aspect above, and will not be repeated here. BRIEF DESCRIPTION OF THE DRAWINGS

[0024] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the following will briefly introduce the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings in the following description are only some embodiments of the present application. For those of ordinary skill in the art, without creative efforts, other drawings can be obtained based on these drawings.

[0025] Figure 1 is a flowchart of a method for dividing the water supply pipe network into partition metering areas provided by an embodiment of the present application;

[0026] Figure 2 is a schematic diagram of morphological erosion reconstruction depression filling of a one-dimensional signal;

[0027] Figure 3 is a flowchart of a method for dividing the water supply pipe network into partition metering areas provided by another embodiment of the present application;

[0028] Figure 4 is a schematic diagram of the DEM depression filling process based on geodesic erosion reconstruction;

[0029] Figure 5 is a schematic diagram of the DEM data after geodesic reconstruction using mathematical morphology;

[0030] Figure 6 It is a schematic diagram of the sub - watershed division result obtained by the hydrological analysis algorithm;

[0031] Figure 7 It is a schematic diagram of the preliminary DMA zoning divided according to the sub - watershed division result;

[0032] Figure 8 It is a distribution map of administrative divisions and settlements within the design scope;

[0033] Figure 9 It is a schematic diagram of the final DMA zoning;

[0034] Figure 10 It is a schematic diagram of the structure of the device for dividing the metering area of the water supply network provided by the embodiment of the present application;

[0035] Figure 11 It is a schematic diagram of the structure of the electronic device provided by the embodiment of the present application. Detailed implementation manners

[0036] In the following description, for the purpose of illustration rather than limitation, specific details such as specific system architectures, technologies, etc. are presented in order to thoroughly understand the embodiments of the present application. However, those skilled in the art should clearly understand that the present application can also be implemented in other embodiments without these specific details. In other cases, detailed descriptions of well - known systems, devices, circuits, and methods are omitted to avoid unnecessary details from interfering with the description of the present application.

[0037] It should be understood that when used in the specification of the present application and the appended claims, the term "comprising" indicates the presence of the described features, wholes, steps, operations, elements, and / or components, but does not exclude the presence or addition of one or more other features, wholes, steps, operations, elements, components, and / or their combinations.

[0038] It should also be understood that the term "and / or" used in the specification of the present application and the appended claims refers to any combination and all possible combinations of one or more of the associated listed items, and includes these combinations.

[0039] References to "one embodiment" or "some embodiments" etc. described in the specification of this application mean that specific features, structures, or characteristics described in connection with that embodiment are included in one or more embodiments of this application. Thus, statements such as "in one embodiment", "in some embodiments", "in other some embodiments", "in still other embodiments", etc. that appear in different places in this specification do not necessarily all refer to the same embodiment, but mean "one or more but not all embodiments", unless otherwise specifically emphasized in other ways. The terms "comprising", "including", "having" and their variants all mean "including but not limited to", unless otherwise specifically emphasized in other ways.

[0040] As Figure 1 shown, the method for dividing the water supply pipe network zoning metering area provided by the embodiment of this application includes the following steps S101 to S106:

[0041] Step S101, obtain the digital elevation model data (Digital Elevation Model, DEM) of the target area.

[0042] Step S102, use the mathematical morphology method to fill the depressions in the digital elevation model data to generate a depression-free elevation model.

[0043] Step S103, based on the hydrological analysis algorithm, perform watershed division on the depression-free elevation model to obtain the river network distribution characteristics and the sub-watershed division results.

[0044] Step S104, combine the river network distribution characteristics and the sub-watershed boundaries to generate a preliminary zoning metering area.

[0045] Step S105, perform a spatial overlay analysis on the sub-watershed division results and the geographic information data to obtain the water supply pipe network management parameters.

[0046] Step S106, adjust the preliminary zoning metering area according to the water supply pipe network management parameters, and output the final zoning result.

[0047] In application, this method can be executed by a Geographic Information System (GIS) platform, and can be specifically implemented through software such as ArcGIS and QGIS.

[0048] In applications, digital elevation model data can be surface elevation raster data obtained through satellite remote sensing or aerial photogrammetry techniques, and its spatial resolution is selected according to the scale of the target area (such as 30 meters × 30 meters). Depression filling uses the geodesic erosion reconstruction algorithm in mathematical morphology to eliminate the closed depression areas in the DEM through iterative calculations. Hydrological analysis algorithms include flow direction calculation, flow accumulation analysis, and watershed division, among which the D8 algorithm is used to determine the water flow direction of each raster for flow direction calculation. Spatial overlay analysis involves vector layer overlay and raster weighted operations. For example, a spatial intersection operation is performed between the sub-watershed boundary and the administrative region boundary, and the zoning range is corrected by combining the settlement density data.

[0049] In this embodiment, the false depressions in the DEM data are accurately eliminated through mathematical morphology, ensuring the accuracy of surface runoff simulation in subsequent hydrological analysis and avoiding zoning boundary deviations caused by data errors. The initial zoning is generated by combining natural terrain features (river networks, sub-watersheds), effectively utilizing the law of terrain continuity, making the zoning boundary highly consistent with natural hydrological units, and reducing the need for manual intervention. Geographical information data is quantified into management parameters through spatial overlay analysis, realizing the scientific transformation from natural geography to management strategies. Compared with traditional empirical zoning methods, the implementability of the zoning scheme is improved, with high reliability and no need for subsequent changes. The final zoning result has both terrain adaptability and management operability, providing support for the construction of urban water supply networks.

[0050] In one embodiment, obtaining the digital elevation model data of the target area includes:

[0051] Downloading or real-time collecting the original elevation data through a geospatial data cloud platform;

[0052] Performing cropping and stitching operations on the original elevation data according to the design scope to generate the digital elevation model data of the target area.

[0053] In applications, DEM data can be downloaded on the Internet such as through a geospatial data cloud platform, or obtained through remote sensing technical means such as oblique photography, and then operations such as cropping and stitching are performed on the DEM data according to the design scope to obtain the DEM data within the design scope as the input for the next step.

[0054] In applications, geospatial data cloud platforms include but are not limited to public data sources such as NASA Earthdata and USGS EarthExplorer, and real-time collection can be achieved by using a LiDAR device carried by an unmanned aerial vehicle. The cropping process uses the mask extraction function of GIS software with the administrative boundary or planning red line as the cropping range. The stitching process is for cross-sheet data, and the seams are eliminated through coordinate system unification and edge matching algorithms.

[0055] The cropping and splicing process in the embodiments of this application adopts an adaptive boundary matching algorithm to eliminate the seam misalignment problem caused by traditional manual splicing, improving the data splicing efficiency. The generated standardized DEM can be directly connected to the subsequent processing module, reducing the data preprocessing time, and is particularly suitable for large-scale cross-regional pipe network planning projects.

[0056] In one embodiment, step S102 includes:

[0057] Step S201: Generate a marked image with the same data range as the digital elevation model data. The boundary pixel values of the marked image are the same as the corresponding pixel values of the digital elevation model data, and the non-boundary pixel values are set to the maximum elevation value in the digital elevation model data.

[0058] Step S202: Take the marked image as the current processing image, and compare the minimum pixel value in the eight-neighborhood of each non-boundary pixel in the current processing image with the corresponding pixel elevation value of the digital elevation model data through iterative calculation, and take the larger value to update the current pixel value until all pixel values converge.

[0059] In application, the marked image is raster data with the same size as the DEM, and the boundary is defined as the outermost circle of pixels in the DEM data matrix. The iterative calculation adopts a row-first traversal strategy, and the eight-neighborhood of each pixel refers to its 8 adjacent pixels (up, down, left, right, and 4 diagonal directions). The convergence determination criterion is that the change amount of all pixel values in two consecutive iterations is less than 0.01 meters. This algorithm gradually raises the elevation of the depression area while keeping the boundary elevation unchanged, and finally generates a continuous terrain surface without closed depressions.

[0060] In application, step S102 is implemented by the geodesic reconstruction method of mathematical morphology to obtain a DEM without false depressions. Depressions in the DEM refer to sets of pixels in the DEM that do not have non-ascending and boundary-connected pixels, usually the regional minimum values in the DEM. Mathematical morphology provides a complete framework for the depression filling algorithm, and can effectively handle even nested depressions and real natural depressions. Geodesic reconstruction means that a bounded image will always converge after a finite number of geodesic transformation iterative calculations until the dilation or contraction of the marked image is completely blocked by the mask image. If the loop iteration continues, the pixel values of the marked image will no longer change. Geodesic reconstruction includes dilation reconstruction and erosion reconstruction.

[0061] In application, taking a one-dimensional signal as an example, the display of geodesic erosion is as Figure 2 shown. Figure 2 is a comparison chart of the geodesic erosion reconstruction results of a one-dimensional signal. Figure 2 (a) The gray rectangular histogram in it is the one-dimensional original signal f, and the solid line in the figure is the marked image signal f m , after morphological erosion reconstruction, the erosion result is shown in Figure 2(as shown in (b). Among them Figure 2 the gray rectangular histogram in (b) is the one-dimensional original signal, and the dashed line is the result after erosion reconstruction. It can be seen from Figure 2 (b) that the depressions located within the original signal are all filled.

[0062] The operation symbol FILL is used for DEM depression filling based on morphology and is defined by erosion reconstruction. The mask image is the original DEM, denoted as f; the marker image is f m , whose boundary pixels are equal to the original DEM data f, and the values of non-boundary pixels are uniformly set to the maximum value h of the pixels in the original DEM data f max .

[0063]

[0064] Among them, when x is on the DEM boundary: f m (x) = f(x); otherwise: f m (x) = h max .

[0065] Using the above formula, read the DEM data, set a marker DEM that is the same as the DEM range, set all pixels of the marker DEM except the boundary pixels to the maximum value of the DEM elevation value, and then start from the non-boundary pixels of the marker DEM, traverse each grid in the order from left to right and from top to bottom, and perform geodesic erosion calculations in turn. Specifically, compare the minimum value of the 8-neighborhood of the center grid of the marker DEM with the elevation of the corresponding grid of the original DEM, and assign the larger value to the center grid. When all grids are traversed and the pixel values no longer change, the geodesic erosion reconstruction ends.

[0066] The embodiment of the present application uses the mathematical morphology method for depression filling. By generating a marker image with the same range as the DEM, the boundary pixel values of which are the same as the DEM, and the non-boundary pixels are set to the maximum value of the DEM; by iteratively calculating and comparing the minimum value of the eight neighborhoods of the non-boundary pixels of the marker image with the original DEM value, and taking the larger value for update until convergence. It effectively eliminates the false depressions caused by the DEM data acquisition error and retains the true terrain features (such as lakes, basins); the filling result meets the requirements of surface runoff simulation in hydrological analysis.

[0067] In one embodiment, as Figure 3 shown, step S202 includes:

[0068] Step S2021: Take the marker image as the current processing image;

[0069] Step S2022: Traverse all non-boundary pixels in the current processing image in row-major order;

[0070] Step S2023: Calculate the minimum pixel value of the eight-neighborhood for each current pixel;

[0071] Step S2024: Compare the minimum pixel value with the elevation value of the corresponding pixel in the digital elevation model data;

[0072] Step S2025: Update the pixel value of the current pixel to the larger value between the minimum pixel value and the elevation value of the corresponding pixel in the digital elevation model data;

[0073] Step S2026: After completing the traversal of the entire image, detect the change status of the pixel values of the currently processed image. If there are pixel value changes, use the updated image as the new currently processed image and return to execute the step of traversing all non-boundary pixels in the currently processed image in row-major order;

[0074] Step S2027: If all pixel values remain unchanged, terminate the iteration and output the processing result as a depression-free elevation model.

[0075] In the application, row-major traversal means accessing each pixel from left to right and top to bottom in turn. The calculation of the eight-neighborhood minimum value uses the sliding window method, taking the minimum value for each 3×3 neighborhood (excluding the central pixel). The pixel update rule is: if the neighborhood minimum value is greater than the original DEM value, update it to the neighborhood minimum value; otherwise, keep the original DEM value. Change detection is achieved by comparing the Euclidean distance of the image difference matrix before and after iteration. When the distance is less than the threshold ε (usually taken as 0.01), it is determined to converge.

[0076] Figure 4 Shows the process of eliminating depressions by DEM erosion reconstruction. Figure 4 (a) is the original DEM, where the numbers represent elevation values, Figure 4 (b) is the labeled DEM, obtained by calculating from the original DEM ( Figure 4 (a)), and the elevation of other grids except the boundary is set to Figure 4 the maximum value of the elevation values in (a), for example, Figure 4 the gray area in (b) is assigned a value of 15. The specific steps of the process of eliminating depressions by DEM erosion reconstruction are as follows:

[0077] ① Starting from the non-boundary pixels of the labeled DEM, traverse each grid in the order from left to right and top to bottom, and perform geodesic erosion calculations in turn.

[0078] ② Read the first pixel in the labeled DEM in the above order, use it as the central grid, compare the minimum value of the 8-neighborhood of the central grid of the labeled DEM with the elevation of the corresponding grid in the original DEM, and assign the larger value to the central grid.

[0079] ③ After traversing all grids, the geodesic erosion reconstruction ends.

[0080] Figure 4 (c) to Figure 4 (e) shows the entire calculation process. As can be seen from the figure, the grid elevation value of C3 in the original DEM is raised from 8 to 9, as Figure 4 shown in (e). It can be seen that the depression is correctly filled to form a depression-free DEM.

[0081] Using this geodesic reconstruction method, after geodesic reconstruction of the original DEM obtained in step S101, the reconstructed DEM obtained is shown in Figure 5 the figure.

[0082] In the embodiment of the present application, the marked image is set as the current processed image; the non-boundary pixels are traversed row by row and the minimum value of the eight-neighborhood is calculated; the minimum value is compared with the corresponding value of the original DEM and the pixel value is updated; the change of the whole image is detected to determine whether to continue the iteration. The traversal order ensures that the calculation process is reproducible, and the pixel-level comparison retains the topographic detail features.

[0083] In one embodiment, step S103 includes:

[0084] Step S301, calculating the water flow direction of each pixel in the depression-free elevation model;

[0085] Step S302, generating a flow accumulation distribution map according to the water flow direction;

[0086] Step S303, extracting the river network distribution characteristics from the flow accumulation distribution map based on a preset flow threshold;

[0087] Step S304, dividing the adjacent sub-watershed boundaries according to the river network distribution characteristics to obtain the sub-watershed division result.

[0088] In application, hydrological analysis is to study the hydrological characteristics of the watershed and simulate the surface hydrological process based on the elevation model. Hydrological analysis mainly studies the movement law of water bodies in the watershed. Common research contents include the source, distribution and movement law of water. Using the hydrological analysis algorithm to divide the watershed of the DEM data without false depressions, the specific operations include flow direction calculation, flow accumulation setting, watershed division, etc. Finally, the sub-watershed division result within the design range is obtained.

[0089] In application, the hydrological analysis operation can be realized by using GIS software. Taking ArcGIS software as an example, using ArcGIS to load the reconstructed DEM, the hydrological analysis toolbox in ArcGIS realizes the work of DEM flow direction calculation, flow accumulation calculation and watershed division. After this work, the sub-watershed division result can be obtained, as shown in Figure 6 the figure.

[0090] In the embodiments of the present application, the water flow direction of each pixel in the DEM without depressions is calculated; a confluence accumulation distribution map is generated; river network features are extracted based on a confluence threshold; and sub-basin boundaries are divided according to the river network. The natural river network distribution is accurately extracted, and the sub-basin division reflects the true hydrological units.

[0091] In one embodiment, step S104 may further include:

[0092] Extracting river network confluence nodes from the sub-basin division results as partition control points;

[0093] Merging the regions connected to the same river network tributary in adjacent sub-basins;

[0094] Determining the natural boundary based on the river network distribution characteristics, and generating a preliminary partition measurement area in combination with the merged sub-basin boundaries.

[0095] In applications, river network confluence nodes are automatically identified through topological analysis and defined as the intersection positions of two or more river channels. The regional merging rule is based on water flow connectivity analysis, and only adjacent sub-basins sharing the same river network tributary are merged. The generation of the natural boundary can be divided based on the river network distribution characteristics or controlled by the confluence threshold. Among them, the river network distribution characteristics (such as the main river channel and tributary positions) can be extracted through the watershed algorithm. For example, the main river channel between two sub-basins is the natural boundary, naturally separating the catchment areas on both sides. Using the confluence threshold for control, a natural boundary is formed by setting a threshold for the confluence accumulation on the tributary to divide the sub-basins, such as Figure 7 the sub-basin corresponding to the middle tributary. Taking the main river network of the trunk stream as the natural boundary framework and the sub-basin boundaries as filling units, a preliminary partition measurement area is generated in combination with the merged sub-basin boundaries as shown in Figure 7 shown.

[0096] In the embodiments of the present application, the control point positioning of the river network confluence nodes ensures the coverage of the hub for partition management; the regional merging rule maintains the balance of the partition scale and avoids the imbalance of resource allocation; the natural boundary highly fits the terrain features and reduces the number of cross-regional connection points of the pipe network.

[0097] In applications, the administrative division and residential area distribution map (i.e., geographical information data) within the design scope are as shown in Figure 8 shown. By performing a spatial overlay analysis on the layer of Figure 8 and the sub-basin division result of Figure 6 with the geographical information data, water supply pipe network management parameters can be obtained, including management authority division data, management area boundary data, end-user water consumption scale data, etc. Finally, the preliminary partition measurement area is checked and adjusted according to the water supply pipe network management parameters, and the final partition result is output as shown in Figure 9 shown. Comparing the positions corresponding to the 1# circle and 2# circle in Figure 9 with Figure 7By comparing the corresponding positions, it is found that there are obvious differences at these two positions. Combining the analysis of the water supply network management parameters, Figure 9 the final sub-region division results obtained are more in line with the requirements of water supply management. The results show that the use of water supply network management parameters can optimize and adjust the water supply sub-region data.

[0098] In the application, the permission division data uses spatial intersection operations to determine the overlapping relationship between sub-watersheds and administrative regions, and generates a responsibility attribution mapping table. The regional boundary data defines different management priority regions through the weighted overlay of density grids and sub-watersheds. The water consumption scale data extracts the distribution information of pipe network nodes, and combines the number of users and service radius to establish a demand prediction model.

[0099] It should be understood that the magnitudes of the sequence numbers of the steps in the above embodiments do not imply the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation to the implementation process of the embodiments of the present application.

[0100] The embodiments of the present application further provide a device for dividing water supply network partition metering areas, which is used to execute the steps in the method embodiments for dividing water supply network partition metering areas. The device for dividing water supply network partition metering areas can be a virtual appliance in an electronic device, which is run by the processor of the electronic device, or the electronic device itself.

[0101] As Figure 10 shown, the device 100 for dividing water supply network partition metering areas provided by the embodiments of the present application includes:

[0102] A data acquisition module 101, configured to acquire digital elevation model data of a target area;

[0103] A reconstruction module 102, configured to fill depressions in the digital elevation model data by using mathematical morphology methods to generate a depression-free elevation model;

[0104] A watershed division module 103, configured to perform watershed division on the depression-free elevation model based on a hydrological analysis algorithm to obtain river network distribution characteristics and sub-watershed division results;

[0105] A preliminary partition module 104, configured to generate preliminary partition metering areas by combining river network distribution characteristics and sub-watershed boundaries;

[0106] An overlay analysis module 105, configured to perform spatial overlay analysis on the sub-watershed division results and geographic information data to obtain water supply network management parameters;

[0107] An adjustment and output module 106, configured to adjust the preliminary partition metering areas according to the water supply network management parameters and output the final partition results.

[0108] As shown Figure 11 in the figure, an embodiment of the present application further provides an electronic device 200, including: at least one processor 201( Figure 11 only one processor is shown in the figure), a memory 202, and a computer program 203 stored in the memory 202 and executable on at least one processor 201. When the processor 201 executes the computer program 203, the steps in the above method embodiments are implemented.

[0109] In applications, the electronic device may include, but is not limited to, a processor and a memory. Those skilled in the art can understand that Figure 11 this is only an example of the electronic device and does not constitute a limitation on the electronic device. It may include more or fewer components than shown in the figure, or combine certain components, or different components.

[0110] In applications, the processor may be a central processing unit (CPU), and the processor may also be other general-purpose processors, digital signal processors (DSPs), application specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. The general-purpose processor may be a microprocessor or the processor may also be any conventional processor, etc.

[0111] In applications, the memory may be an internal storage unit of the electronic device in some embodiments, such as the hard disk or memory of the electronic device. The memory may also be an external storage device of the electronic device in other embodiments, such as a plug-in hard disk, a smart media card (SMC), a secure digital (SD) card, a flash card, etc. equipped on the electronic device. Further, the memory may also include both the internal storage unit and the external storage device of the electronic device. The memory is used to store an operating system, application programs, a boot loader, data, and other programs, such as the program code of the computer program. The memory may also be used to temporarily store data that has been output or will be output.

[0112] It should be noted that for the content such as information interaction and execution process between the above devices / units, since it is based on the same concept as the method embodiments of the present application, for its specific functions and the technical effects brought, please refer to the method embodiment part for details and will not be elaborated here.

[0113] Those skilled in the art can clearly understand that, for the convenience and conciseness of description, only the above division of each functional unit and module is used as an example. In actual applications, the above functions can be allocated to different functional units and modules as needed, that is, the internal structure of the device can be divided into different functional units or modules to complete all or part of the functions described above. Each functional unit and module in the embodiment can be integrated into a processing unit, or each unit can exist physically alone, or two or more units can be integrated into one unit. The above integrated unit can be implemented in the form of hardware or in the form of a software functional unit. In addition, the specific names of each functional unit and module are only for the convenience of mutual distinction and do not limit the protection scope of this application. The specific working process of the units and modules in the above system can refer to the corresponding process in the foregoing method embodiment and will not be elaborated herein.

[0114] The embodiment of this application also provides a computer-readable storage medium. The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps in the foregoing method embodiments can be implemented.

[0115] The embodiment of this application provides a computer program product, including a computer program. When the computer program product runs on an electronic device, the electronic device can execute the steps in the foregoing method embodiments.

[0116] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on such an understanding, to implement all or part of the processes in the above method embodiments of this application, a computer program can be used to instruct relevant hardware to complete. The computer program can be stored in a computer-readable storage medium. When the computer program is executed by a processor, the steps in the foregoing method embodiments can be implemented. Among them, the computer program includes computer program code, and the computer program code can be in the form of source code, object code, executable file or some intermediate form, etc. The computer-readable medium can at least include: any entity or device that can carry the computer program code to the device / electronic device, recording medium, computer memory, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), electrical carrier signal, telecommunication signal, and software distribution medium. For example, a USB flash drive, a mobile hard disk, a magnetic disk or an optical disc, etc. In some jurisdictions, according to legislation and patent practice, the computer-readable medium cannot be an electrical carrier signal and a telecommunication signal.

[0117] In the above embodiments, the descriptions of the various embodiments each have their own emphasis. For parts not described in detail or recorded in a certain embodiment, reference may be made to the relevant descriptions of other embodiments.

[0118] Those of ordinary skill in the art can realize that the units and algorithm steps of the examples described in combination with the embodiments disclosed herein can be implemented by electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are executed in a hardware or software manner depends on the specific application and design constraints of the technical solution. A professional technician can use different methods to implement the described functions for each specific application, but such implementation should not be considered to exceed the scope of this application.

[0119] In the embodiments provided in this application, it should be understood that the disclosed devices and methods can be implemented in other ways. For example, the device embodiments described above are merely illustrative. For example, the division of modules or units is only a logical function division. In actual implementation, there may be other division methods. For example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the displayed or discussed couplings or direct couplings or communication connections to each other can be through some interfaces. The indirect couplings or communication connections of devices or units can be in electrical, mechanical or other forms.

[0120] The units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they can be located in one place, or can be distributed to multiple network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0121] The above embodiments are only used to illustrate the technical solutions of this application, rather than to limit them; although this application has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that they can still modify the technical solutions recorded in the foregoing embodiments, or perform equivalent replacements for some of the technical features; and these modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of this application, and should all be included within the protection scope of this application.

Claims

1. A method for dividing water supply network into zones and metering areas, characterized in that: include: Obtain digital elevation model data of the target area; Using mathematical morphology method to fill depressions in the digital elevation model data to generate a depression-free elevation model; Based on the hydrological analysis algorithm, the watershed division is performed on the depression-free elevation model to obtain the river network distribution characteristics and sub-basin division results; Generating preliminary zoning measurement areas in combination with the river network distribution characteristics and sub-basin boundaries; Perform spatial overlay analysis on the sub-basin division results and geographic information data to obtain water supply network management parameters; The preliminary zoning metering area is adjusted according to the water supply network management parameters, and the final zoning result is output.

2. The method for dividing water supply network into zone metering areas according to claim 1, characterized in that: The step of obtaining digital elevation model data of the target area includes: Download or collect raw elevation data in real time through the geospatial data cloud platform; The original elevation data is clipped and spliced ​​according to the design range to generate digital elevation model data of the target area.

3. The method for dividing water supply network into zone metering areas according to claim 1, characterized in that: The method of using a mathematical morphology method to fill depressions in the digital elevation model data to generate a depression-free elevation model includes: Generate a marker image with the same range as the digital elevation model data, wherein the boundary pixel value of the marker image is the same as the corresponding pixel value of the digital elevation model data, and the non-boundary pixel value is set to the maximum elevation value in the digital elevation model data; The marked image is used as the current processing image, and the minimum pixel value of the eight-neighborhood of each non-boundary pixel in the current processing image is compared with the corresponding pixel elevation value of the digital elevation model data through iterative calculation, and the larger value is taken to update the current pixel value until all pixel values ​​converge.

4. The method for dividing water supply network into zone metering areas according to claim 3, characterized in that: The method of using the marked image as the current processing image, comparing the minimum pixel value of the eight neighborhoods of each non-boundary pixel in the current processing image with the corresponding pixel elevation value of the digital elevation model data through iterative calculation, and taking the larger value to update the current pixel value until all pixel values ​​converge, includes: Using the marked image as the current processing image; Traverse all non-boundary pixels in the current processing image in row-major order; For each current pixel, calculate the minimum pixel value of its eight neighborhoods; Comparing the minimum pixel value with the elevation value of the corresponding pixel in the digital elevation model data; Updating the pixel value of the current pixel to the larger value between the minimum pixel value and the elevation value of the corresponding pixel in the digital elevation model data; After the whole image traversal is completed, the pixel value change state of the current processing image is detected. If there is a pixel value change, the updated image is used as the new current processing image, and the step of traversing all non-boundary pixels in the current processing image in row priority order is returned to be executed; If all pixel values ​​remain unchanged, the iteration is terminated and the processing result is output as a depression-free elevation model.

5. The method for dividing water supply network into zone metering areas according to claim 1, characterized in that: The watershed division of the non-depression elevation model based on the hydrological analysis algorithm to obtain the river network distribution characteristics and sub-basin division results include: Calculating the water flow direction of each pixel in the depression-free elevation model; Generate a flow accumulation distribution diagram according to the water flow direction; Extracting river network distribution characteristics from the runoff accumulation distribution map based on a preset runoff threshold; The adjacent sub-basin boundaries are divided according to the river network distribution characteristics to obtain the sub-basin division results.

6. The method for dividing water supply network into zone metering areas according to claim 1, characterized in that: The generation of preliminary zoning measurement areas by combining the river network distribution characteristics with the sub-basin boundaries includes: Extracting river network intersection nodes from the sub-basin division results as partition control points; Merge the areas in adjacent sub-basins that are connected to the same river network tributaries; The river network distribution characteristics are used as natural boundaries, combined with the merged sub-basin boundaries to generate preliminary zoning measurement areas.

7. The method for dividing water supply network into zone metering areas according to claim 1, characterized in that: The water supply network management parameters include the management authority division data, the management area boundary data and the end-user water consumption scale data.

8. The water supply network zone metering area division device according to claim 7, characterized in that: include: A data acquisition module is used to acquire digital elevation model data of a target area; A reconstruction module, used for filling depressions in the digital elevation model data using a mathematical morphology method to generate an elevation model without depressions; A watershed division module is used to perform watershed division on the depression-free elevation model based on a hydrological analysis algorithm to obtain river network distribution characteristics and sub-watershed division results; A preliminary zoning module, for generating preliminary zoning metering areas by combining the river network distribution characteristics and sub-basin boundaries; An overlay analysis module is used to perform spatial overlay analysis on the sub-basin division results and geographic information data to obtain water supply network management parameters; The adjustment output module is used to adjust the preliminary partition metering area according to the water supply network management parameters and output the final partition result.

9. An electronic device, characterized in that: The electronic device comprises a processor, a memory, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, the electronic device implements the method as claimed in any one of claims 1 to 7.

10. A computer program product, characterized in that The invention comprises a computer program, which, when being executed, enables the method according to any one of claims 1 to 7 to be performed.

Images