A state-owned asset supervision method based on geographic information layering

Abstract

This invention discloses a method for supervising state-owned assets based on geographic information layering, belonging to the field of supervision and management technology. The method includes: generating a digital elevation model based on the centerline of a flood control dike; generating an elevation function based on the digital elevation model; calculating intermediate elevations based on the elevation function; generating candidate shoulder stages based on the intermediate elevations; generating planar coordinates based on the candidate shoulder stages; generating a shoulder centerline based on the planar coordinates; calculating the average shoulder width based on the shoulder centerline; and calculating the total land area occupied by the hidden inner shoulders of the flood control dike based on the average shoulder width. This invention improves the completeness and precision of state-owned asset supervision.

Description

Technical Field

[0001] This invention relates to the field of regulatory management technology, and in particular to a method for supervising state-owned assets based on geographic information layering. Background Technology

[0002] In the field of flood control engineering and water conservancy infrastructure management, flood control dikes, as important flood control and disaster reduction engineering facilities, are usually linearly distributed along rivers, lakes, or coastal areas. Their spatial morphology is complex, their structural scale is large, and their coverage is wide. During long-term operation, they not only undertake public safety functions such as flood control and drainage, but also constitute an important component of state-owned assets. With the development of refined land space management and digital governance, spatial representation, area measurement, and asset supervision of flood control dikes and their ancillary structures based on geographic information data has become an important application scenario in water conservancy engineering management and state-owned asset management. Especially on the backwater side of the dike, due to historical construction conditions, topographical evolution, and multiple reinforcement and renovations, there are often a large number of irregularly shaped, slowly changing elevation areas on the inner dike shoulder that have long occupied land resources. These areas exist in the engineering entity and are spatially continuous, but are difficult to fully represent in conventional drawings and ledgers. Their land area and location information lack unified, precise, and verifiable quantitative methods, thus posing a significant challenge to the accurate registration, dynamic supervision, and value assessment of flood control engineering assets.

[0003] In existing technologies, when calculating and supervising the area of ​​flood control dikes and related assets, the methods usually rely on manual experience, simplified geometric assumptions, or estimation based on a single cross-sectional parameter. This makes it difficult to fully reflect the topographical changes along the dike and the continuous spatial form of the inner shoulder. In particular, for structures such as the hidden inner shoulder that are not easy to identify directly, their area is often ignored or roughly merged, resulting in a discrepancy between the asset measurement results and the actual situation. Summary of the Invention

[0004] The purpose of this invention is to address the shortcomings of existing technologies that fail to adequately reflect the topographical changes along the embankment and the continuous spatial morphology of the inner embankment shoulders, and to propose a method for supervising state-owned assets based on geographic information layering.

[0005] To address the problems existing in the prior art, the present invention adopts the following technical solution:

[0006] A geographic information-based hierarchical method for supervising state-owned assets includes:

[0007] S1. Generate a digital elevation model based on the centerline of the flood control embankment, and generate an elevation function based on the digital elevation model;

[0008] S2. Calculate intermediate elevation based on elevation function, and generate candidate abutment platforms based on intermediate elevation.

[0009] S3. Generate planar coordinates based on the candidate platform stage of the embankment shoulder, and generate the centerline of the embankment shoulder based on the planar coordinates;

[0010] S4. Calculate the average shoulder width based on the centerline of the shoulder, and calculate the total land area occupied by the hidden inner shoulder in the flood control dike based on the average shoulder width.

[0011] S5. Update the monitoring and measurement of flood control dike assets based on the total land area of ​​the hidden inner shoulder.

[0012] Preferably, generating a digital elevation model based on the centerline of the flood control dike includes:

[0013] Obtain the centerline of the flood control dike;

[0014] Locate the adjacent surface area on the backwater side of the dike's centerline;

[0015] Obtain surveying data of the dike centerline and adjacent surface areas;

[0016] Elevation data was collected from the spatial area corresponding to the centerline of the dike based on the survey data to obtain elevation sampling data.

[0017] Spatial interpolation and rasterization are performed on the elevation sampling data to obtain a digital elevation model.

[0018] Preferably, generating an elevation function based on a digital elevation model includes:

[0019] The centerline of the dike was divided into segments, resulting in multiple mileage points;

[0020] Determine the tangent direction of the dike's centerline at the kilometer point;

[0021] Construct profile lines at each mileage point in the normal direction of the tangent and facing the backwater side;

[0022] Multiple positional parameters are obtained by sampling the profile lines at equal intervals.

[0023] Based on the digital elevation model, the elevation of the location parameters is queried to obtain the elevation value of each location parameter;

[0024] A functional relationship model is performed on the location parameters and their elevation values ​​to obtain the elevation function of the profile line.

[0025] Preferably, calculating the intermediate elevation based on the elevation function includes:

[0026] Based on the elevation function, the elevation values ​​corresponding to the position parameters on the side closest to the center line of the dike are extracted to obtain the first set of elevation values;

[0027] The maximum value of the first set of elevation values ​​is calculated to obtain the top elevation of the embankment of the profile line;

[0028] Based on the elevation function, the elevation values ​​corresponding to the position parameters on the side away from the center line of the dike are extracted to obtain the second set of elevation values;

[0029] The average value of the second set of elevation values ​​is calculated to obtain the backwater average elevation of the profile line.

[0030] The average elevation of the dike crest and the average elevation of the backwater are calculated to obtain the intermediate elevation.

[0031] Preferably, the stage of generating candidate abutment platforms based on intermediate elevations includes:

[0032] The elevation function values ​​are divided into elevation intervals based on the intermediate elevation to obtain the intermediate elevation level.

[0033] Based on the intermediate elevation, all positional parameters of the profile line are filtered to obtain a set of candidate positional parameters;

[0034] Based on the elevation function, the continuity analysis of the function values ​​corresponding to adjacent candidate position parameters in the candidate position parameter set is performed to obtain the candidate platform stage of the embankment shoulder.

[0035] Preferably, generating planar coordinates based on the candidate platform stage of the embankment includes:

[0036] The starting position is extracted for the candidate platform stage of the embankment.

[0037] The termination position is extracted from the candidate platform stage of the embankment.

[0038] The midpoint position is obtained by calculating the midpoint between the starting and ending positions.

[0039] Perform a planar coordinate transformation on the midpoint position to obtain the planar coordinates.

[0040] Preferably, generating the centerline of the embankment shoulder based on planar coordinates includes:

[0041] The average elevation of the embankment shoulder is obtained by averaging the elevation values ​​corresponding to all position parameters during the candidate platform stage of the embankment shoulder.

[0042] Use plane coordinates and average elevation of the embankment shoulder as attribute information;

[0043] Based on attribute information, spatial points corresponding to planar coordinates are marked to obtain candidate points for the embankment shoulder;

[0044] Sort all the profile lines to obtain an ordered sequence of profiles;

[0045] Adjacent association processing is performed on candidate points of the embankment shoulder based on the ordered sequence of the profile to obtain a sequence of candidate points of adjacent embankment shoulders;

[0046] Connect the candidate point sequences of adjacent embankment shoulders to obtain the centerline of the embankment shoulders.

[0047] Preferably, the average shoulder width is calculated based on the centerline of the shoulder, including:

[0048] The point sequence along the centerline of the embankment shoulder is paired up to obtain multiple candidate point pairs for adjacent embankment shoulders.

[0049] The adjacent candidate points of the embankment shoulder are divided into line segments to obtain multiple smaller segments;

[0050] Obtain the first shoulder width of the candidate shoulder point corresponding to one of the endpoints of a small segment of adjacent shoulder candidate points;

[0051] Obtain the second shoulder width of the shoulder candidate point corresponding to the other endpoint in the adjacent shoulder candidate point pair;

[0052] The average shoulder width is obtained by averaging the widths of the first and second shoulders.

[0053] Preferably, the total land area occupied by the concealed inner shoulder of the flood control dike is calculated based on the average shoulder width, including:

[0054] The distance between each pair of adjacent candidate points on the embankment shoulder is calculated to obtain the length of the segment.

[0055] Multiply the length of the short segment by the average width of the embankment shoulder to obtain the area occupied by the embankment shoulder.

[0056] The total area occupied by the hidden inner shoulder is obtained by summing up the areas occupied by all the shoulders.

[0057] Compared with the prior art, the beneficial effects of the present invention are:

[0058] 1. This invention constructs a digital elevation model based on the centerline of the dike and generates a profile elevation function, enabling the terrain on the back side of the flood control dike to be expressed in a continuous and computable manner. This transforms the terrain information that was originally scattered in the survey data into an analytical object with spatial continuity, allowing structures such as the hidden inner dike shoulder, which are difficult to identify intuitively in conventional management, to be stably extracted at the profile scale, providing a reliable data foundation for subsequent area calculations.

[0059] 2. This invention introduces an intermediate elevation calculation method, using the intermediate height layer between the top of the dike and the surface height on the back side as the basis for analysis. Within this height layer, the profile elevation changes are screened and the continuity is analyzed, thereby automatically obtaining the candidate stage of the dike shoulder. This makes the identification process of the dike shoulder structure consistent and repeatable, effectively avoiding the uncertainty caused by experience judgment in the prior art.

[0060] 3. This invention further constructs a centerline of the embankment shoulder by associating multiple candidate points of the embankment shoulder on multiple cross sections along the embankment line. Based on this, the area of ​​each segment is calculated by combining the width of the embankment shoulder with the length along the line. Finally, the total area occupied by the hidden inner embankment shoulder is obtained. This allows the long-neglected area of ​​the inner embankment shoulder in flood control embankments to be quantified and incorporated into the state-owned asset measurement system. This achieves consistent updating of flood control embankment asset data with actual spatial form, and improves the completeness and precision of state-owned asset supervision. Attached Figure Description

[0061] The accompanying drawings, which are included to provide a further understanding of the invention and form part of this application, illustrate exemplary embodiments of the invention and, together with their description, serve to explain the invention and do not constitute an undue limitation thereof. In the drawings:

[0062] Figure 1 This is a flowchart illustrating a method for supervising state-owned assets based on geographic information layering, as provided in an embodiment of the present invention. Detailed Implementation

[0063] The technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings of the embodiments of the present invention. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments.

[0064] Example: This example provides a method for supervising state-owned assets based on geographic information layering. See [link to example]. Figure 1 Specifically, including:

[0065] S1. Generate a digital elevation model based on the centerline of the flood control embankment, and generate an elevation function based on the digital elevation model;

[0066] In an embodiment of the present invention, generating a digital elevation model based on the centerline of a flood control embankment includes:

[0067] Obtain the centerline of the flood control dike;

[0068] Specifically, obtaining the centerline of a flood control dike includes the process of determining a linear element representing the overall orientation of the dike based on its engineering structure and spatial location characteristics. This process first identifies the continuous extension direction of the dike in planar space based on its actual distribution on the ground surface, and selects a series of spatial location points along this direction to reflect the axial position of the dike. These spatial location points can originate from dike engineering design data, as-built data, or spatial measurement results of the existing dike entity. Subsequently, these spatial location points are connected sequentially along the dike's extension direction to form a continuous linear element located in the middle of the dike entity, reflecting the overall orientation and location distribution of the flood control dike in space. After forming the linear element, a spatial consistency check is performed to ensure its continuity and smoothness throughout the entire flood control dike area, thereby obtaining a centerline that stably represents the spatial location and extension direction of the flood control dike. This centerline serves as the basic spatial data for subsequent location of the backwater area, construction of profile lines, and topographic analysis.

[0069] Locate the adjacent surface area on the backwater side of the dike's centerline;

[0070] Specifically, locating the adjacent surface area on the backwater side of the dike's centerline involves determining the dike's spatial orientation based on the acquired dike centerline, using this orientation as a reference direction to distinguish the upstream and downstream sides of the dike. After clarifying the direction of the downstream side, a surface area directly connected to the dike entity and spatially continuous is determined along the downstream side direction from the dike centerline. This surface area covers the toe of the dike's downstream slope and its outer extension area. The extension area is bounded by a boundary that can fully reflect the topographic changes on the downstream side of the dike. By advancing segment by segment along the dike centerline, a continuous downstream adjacent surface area corresponding to the dike centerline is formed for subsequent topographic elevation analysis.

[0071] Obtain surveying data of the dike centerline and adjacent surface areas;

[0072] Specifically, acquiring surveying data for the dike centerline and adjacent surface areas includes collecting measurement results describing the spatial location and elevation changes of the area for the determined dike centerline and its adjacent surface areas on the backwater side. The surveying data includes plane coordinate data to express the location of the dike centerline and elevation data to express the topographic relief of the adjacent surface areas. The surveying data is obtained by spatially measuring the dike entity and its backwater surface, so that each measurement location has corresponding plane position and elevation information. The surveying data is then uniformly organized into a data set that can simultaneously reflect the spatial location of the dike centerline and the topographic features of the adjacent surface areas, serving as the basis for subsequent elevation data acquisition and digital elevation model construction.

[0073] Specifically, the centerline of a flood control dike refers to a linear spatial element set along the overall direction of the dike. This line reflects the location and extension direction of the dike in planar space, usually corresponding to the location of the dike crest or axis, and is used to describe the overall geometric shape of the dike project. The backwater side of the dike centerline refers to the spatial direction located inside the dike or on the land side relative to the water-facing side of the dike. The area in this direction is usually not in direct contact with the water body. The adjacent surface area refers to the surface range that is adjacent to the backwater side of the dike centerline and is spatially continuous with the physical structure of the flood control dike. The topography of this area is closely related to the evolution of the dike structure. Surveying data refers to the set of information obtained through measurement methods to describe the spatial location and topographic elevation of the dike centerline and adjacent surface areas. This data can reflect the real surface undulation and spatial distribution characteristics of the area, and is used for subsequent construction of digital elevation models and analysis of the dike physical structure.

[0074] Elevation data was collected from the spatial area corresponding to the centerline of the dike based on the survey data to obtain elevation sampling data.

[0075] Specifically, using the spatial range of the dike centerline and adjacent surface areas obtained in the previous step, the elevation information recorded in the surveying data is extracted point by point within this spatial range. The elevation information is the surface height value obtained through field measurement. Taking the dike centerline and its adjacent area on the backwater side as the boundary, the planar position and corresponding elevation value of all measurement points covering this area are read together. After reading, all measurement points are sorted according to their planar position, so that each measurement point forms a discrete set of elevation values ​​in space that can reflect the topographic relief of the area. This discrete set is used as elevation sampling data. The elevation sampling data contains the coordinates and corresponding elevation values ​​of a series of specific measurement points, which are used to construct an elevation model that continuously expresses the topography.

[0076] Spatial interpolation and rasterization are performed on the elevation sampling data to obtain a digital elevation model.

[0077] Specifically, by utilizing the planar positions and corresponding elevation values ​​of each measurement point in the elevation sampling data, the elevations of spatial locations without measurement points are extrapolated, transforming the originally discrete elevation distribution into a continuous spatial representation. Spatial interpolation calculates the elevations of surrounding locations based on the actual height of the measurement points. After spatial interpolation, the entire spatial area is divided into multiple grid cells with fixed planar ranges according to a fixed grid size. For each grid cell, the interpolated continuous elevation data is used to determine its corresponding elevation value. The planar positions and elevation values ​​of all grid cells are combined to form a digital elevation model that represents the topographic relief of the entire dike centerline and adjacent surface areas, serving as the foundational data for subsequent profile construction and hidden dike shoulder identification.

[0078] Specifically, elevation data acquisition refers to the process of obtaining elevation information at various locations on the ground surface within the spatial area corresponding to the centerline of the dike based on surveying data. This elevation information reflects the undulation of the ground surface relative to a unified reference surface. Elevation sampling data refers to the set of discrete elevation information obtained through the elevation data acquisition. This set consists of multiple sampling points with spatial locations and corresponding elevation values, used to describe the topography of the dike and its adjacent areas. Spatial interpolation refers to the process of extrapolating the elevation of locations not directly sampled based on existing elevation sampling data, so that the discrete elevation information forms a continuous distribution in space. Rasterization refers to dividing the continuous spatial elevation distribution into regular grid units, so that each grid unit corresponds to a specific elevation value. The digital elevation model refers to the data model formed after spatial interpolation and rasterization of the elevation sampling data. This model is used to express the overall topographic elevation changes of the dike centerline and adjacent surface areas.

[0079] In an embodiment of the present invention, generating an elevation function based on a digital elevation model includes:

[0080] The centerline of the dike was divided into segments, resulting in multiple mileage points;

[0081] Specifically, segmenting the dike centerline involves selecting a series of discrete points representing the dike's location along the acquired centerline based on spatial distance. The selection method involves advancing along the dike centerline from a starting position at fixed spatial intervals, recording the spatial coordinates of the dike centerline at each interval, thus forming multiple mileage points arranged along the dike's extension direction. These mileage points serve as reference positions for subsequent profile construction. This method allows the continuous spatial path of the dike centerline to be represented by a sequence of discrete points usable for calculation and analysis, ensuring that each mileage point is arranged sequentially along the dike's direction, thereby completely covering the entire range of the dike centerline.

[0082] Determine the tangent direction of the dike's centerline at the kilometer point;

[0083] Specifically, using the obtained multiple mileage points, according to the arrangement order of the mileage points on the center line of the dike, for each mileage point, select its two adjacent mileage points before and after it, and calculate the connection direction between the mileage point and its two adjacent mileage points to obtain the extension direction of the dike center line at that position. This extension direction is taken as the tangent direction at that mileage point. At the starting point and ending point of the dike center line, the direction is calculated by using its adjacent mileage points to obtain the corresponding tangent direction.

[0084] Construct profile lines at each mileage point in the normal direction of the tangent and facing the backwater side;

[0085] Specifically, at each mileage point, using the previously determined tangent direction, the normal direction is obtained by calculating the direction perpendicular to the tangent direction. Based on the backwater direction determined in the previous step, the normal direction is adjusted to point towards the backwater side of the dike. After obtaining the corresponding normal direction, a linear path with a length sufficient to cover the range of topographic changes on the backwater side of the dike is generated from the spatial position of the mileage point along the normal direction towards the backwater side. This path is used as the profile line at the mileage point. The profile line is used to represent the spatial sampling direction that gradually extends from the dike towards the backwater side, laying the foundation for subsequent elevation sampling along the profile line and the construction of the profile elevation function.

[0086] Specifically, segmenting the dike centerline involves selecting multiple points representing the dike centerline's location at certain spatial intervals along the established centerline. These selected points are called mileage points, which are linearly distributed in space to represent the discrete spatial locations of the dike centerline at different positions. The tangent direction of the dike centerline at each mileage point refers to the direction vector formed at that mileage point along the dike centerline. This direction vector reflects the dike's extension direction at that location and is used to distinguish the left and right spatial ranges of the dike. The normal direction of the tangent direction is the direction vector perpendicular to the tangent direction, used to construct an analysis direction consistent with the dike centerline direction and capable of unfolding a profile in space. Constructing a profile line along the normal direction of the tangent direction and towards the backwater side involves extending a linear spatial path along the normal direction and pointing towards the backwater side of the dike at each mileage point. This path represents the elevation sampling direction extending from the dike body towards the backwater side, so as to generate continuous profile data in this direction, thereby describing the topographic changes in the dike body's transverse direction.

[0087] Multiple positional parameters are obtained by sampling the profile lines at equal intervals.

[0088] Specifically, on each profile line constructed along the backwater side of the dike, the starting point of the profile line is first determined, which is located at the spatial coordinates of the corresponding mileage point. Then, the profile line is extended segment by segment towards the backwater side at fixed spatial intervals, and the distance value of that position on the profile line is recorded at each interval. This distance value is used as a position parameter, so that the continuous spatial path on the profile line is represented as a sequence of multiple discrete sampling points with clear distance information. Through this sampling method, the extension process of the profile from the dike to the backwater side is discretized and transformed into multiple spatial locations with directly queryable elevations, laying the foundation for extracting the height values ​​of each location based on the digital elevation model.

[0089] Based on the digital elevation model, the elevation of the location parameters is queried to obtain the elevation value of each location parameter;

[0090] Specifically, the spatial location of each sampling point on the profile line, represented by the location parameters, is used to find the corresponding ground elevation in the digital elevation model (DEM). The DEM consists of regular grid cells, each with a unique range of planar coordinates and a corresponding elevation value. Therefore, given the planar coordinates corresponding to the location parameters, these coordinates are located within the grid of the DEM, and the ground elevation corresponding to that location is directly obtained from the grid cell. This elevation is then used as the elevation value of the location parameter. By performing the above operation on all location parameters one by one, a complete discrete elevation dataset consisting of location parameters and their corresponding elevation values ​​is obtained, thus fully reflecting the elevation changes of the profile line in the backwater direction.

[0091] A functional relationship model is performed on the location parameters and their elevation values ​​to obtain the elevation function of the profile line.

[0092] Specifically, based on the arrangement order of the various location parameters along the profile line, each location parameter and its corresponding elevation value are sequentially associated in order from the profile starting point to the backwater side. The location parameter is used as the independent variable, and the corresponding elevation value as the dependent variable. By combining all the associations, a continuous expression describing the trend of profile line elevation changes with location is formed. This continuous expression is the profile line elevation function, used to reflect the height changes at various points along the profile line, mathematically transforming the originally discrete elevation sampling points into a continuous height change curve. Through this functional expression, further analysis of profile elevation changes can be conducted to identify potential abutment platform structures.

[0093] Specifically, isometric sampling of the profile line refers to selecting multiple points representing different positions on each profile line constructed along the normal direction of the dike centerline at fixed spatial intervals. The positions of these points are represented by their distance from the starting position of the profile line, used to describe the discrete sampling positions of the profile in the backwater direction. The position parameter refers to the distance value used to identify each selected sampling point on the profile line, reflecting the spatial position of the point on the profile line. Elevation query based on the position parameter using the digital elevation model refers to using the spatial position determined by the position parameter to read the corresponding ground height from the digital elevation model, thereby obtaining the elevation value of each position parameter. The elevation value of the position parameter refers to the ground height information corresponding to each sampling point on the profile line, used to describe the topographic changes of the profile line in the backwater direction. Functional relationship modeling of the position parameter and its elevation value refers to sequentially associating the corresponding elevation values ​​according to the arrangement order of each position parameter on the profile line, forming an elevation function that can describe the height change trend of the profile line along the sampling direction.

[0094] Specifically, the elevation function is a mathematical expression used to describe the height variation of a profile line along its sampling direction. This function uses the positional parameters on the profile line as independent variables and the height value obtained by the digital elevation model at that location as the dependent variable. By associating the height information of each sampling location in spatial order, the topographic undulation of the profile line in the direction of the backwater side can be expressed as a continuous curve. The elevation function can reflect the height variation characteristics of the profile as it extends from the dike body to the backwater side, and is an important data basis for identifying hidden dike shoulder structures and analyzing the lateral morphology of dikes.

[0095] S2. Calculate intermediate elevation based on elevation function, and generate candidate abutment platforms based on intermediate elevation.

[0096] In an embodiment of the present invention, calculating the intermediate elevation based on the elevation function includes:

[0097] Based on the elevation function, the elevation values ​​corresponding to the position parameters on the side closest to the center line of the dike are extracted to obtain the first set of elevation values;

[0098] Specifically, the range of location parameters for the starting point of the profile line is determined in the constructed profile line elevation function. The starting point is located at the mileage point and is spatially aligned with the center line of the dike. Then, each location parameter is traversed sequentially from the starting point of the profile line towards the backwater side. A continuous range of location parameters near the starting point is selected as the location parameter set on the side closer to the center line of the dike. The corresponding elevation values ​​of each location parameter set in the elevation function are read one by one. The read elevation values ​​are then collected to form a first elevation value set, which represents the height distribution of the profile line in the vicinity of the dike.

[0099] The maximum value of the first set of elevation values ​​is calculated to obtain the top elevation of the embankment of the profile line;

[0100] Specifically, all elevation values ​​in the first set of elevation values ​​are compared one by one, and the elevation value with the largest value is determined as the maximum elevation value. This maximum elevation value is then used as the top elevation of the embankment corresponding to the profile line. At the same time as determining the top elevation, the corresponding positional parameters that generate the maximum elevation value are recorded, so that the top elevation corresponds to the specific spatial position on the profile line, thereby enabling the top elevation to represent the highest height of the profile when it passes through the embankment area.

[0101] Based on the elevation function, the elevation values ​​corresponding to the position parameters on the side away from the center line of the dike are extracted to obtain the second set of elevation values;

[0102] Specifically, a continuous range of position parameters is determined at the end of the profile line. The end is located in the direction away from the mileage point of the profile line and corresponds to the surface area outside the backwater side. Then, the corresponding elevation values ​​in the elevation function are read one by one from the continuous range of position parameters near the end, and the elevation values ​​are collected to form a second set of elevation values. The second set of elevation values ​​represents the height distribution of the profile line in the surface area outside the backwater side, providing input data for the subsequent calculation of the average elevation of the backwater side.

[0103] Specifically, the location parameters on the side closer to the dike centerline refer to a set of location parameters near the starting point of the profile line. The spatial location corresponding to these location parameters is near the dike body and adjacent to the dike crest area. The location parameters on the side farther from the dike centerline refer to a set of location parameters near the end of the profile line. The spatial location corresponding to these location parameters is further out on the backwater side and closer to the natural ground area. The first set of elevation values ​​refers to the set of elevation values ​​corresponding to the location parameters on the side closer to the dike centerline extracted from the elevation function. This set reflects the height distribution of the profile line in the area near the dike body. The dike crest elevation of the profile line refers to the elevation value obtained by taking the maximum value in the first set of elevation values. This elevation value corresponds to the highest position of the profile line when it passes through the dike body area and represents the height of the dike crest. The second set of elevation values ​​refers to the set of elevation values ​​corresponding to the location parameters on the side farther from the dike centerline extracted from the elevation function. This set reflects the height distribution of the profile line in the surface area outside the backwater side and is used to characterize the height level of the ground on the backwater side.

[0104] The average value of the second set of elevation values ​​is calculated to obtain the backwater average elevation of the profile line.

[0105] Specifically, each elevation value in the second set of elevation values ​​is summed, the number of elevation values ​​in the second set of elevation values ​​is recorded, and the summation result is divided by the number to obtain the average elevation value. This average elevation value is used as the backwater average elevation of the profile line. In the calculation process, the elevation values ​​in the second set of elevation values ​​are all derived from the elevation values ​​corresponding to the continuous position parameters on the side of the profile line away from the center line of the dike, so that the backwater average elevation can characterize the overall height level of the profile line in the surface area outside the backwater side.

[0106] The average elevation of the dike crest and the average elevation of the backwater are calculated to obtain the intermediate elevation.

[0107] Specifically, the top elevation of the dike and the average elevation of the backwater of the profile line are read, and the two are added together to obtain the elevation sum. The elevation sum is then divided by two to obtain the intermediate elevation. The intermediate elevation is located between the top elevation of the dike and the average elevation of the backwater, and is used as a reference height value for screening candidate position parameters in the intermediate height layer in the profile line elevation function.

[0108] Specifically, the average elevation of the backwater refers to the height value calculated by averaging the various elevation values ​​in the second set of elevation values, used to characterize the overall height level of the profile line on the outer side of the ground. The top elevation of the dike refers to the height value representing the top position of the dike body determined from the elevation values ​​on the side of the profile line closest to the dike centerline, and this height value corresponds to the highest point of the dike body on this profile. The intermediate elevation refers to the height value calculated by averaging the top elevation of the dike body and the average elevation of the backwater, and this height value is between the height of the top of the dike body and the height of the ground on the backwater side, used to characterize the height levels of the dike shoulder stage that may appear on the profile line.

[0109] In an embodiment of the present invention, the stage of generating candidate abutment platforms based on intermediate elevations includes:

[0110] The elevation function values ​​are divided into elevation intervals based on the intermediate elevation to obtain the intermediate elevation level.

[0111] Specifically, the intermediate elevation of the profile line and the elevation values ​​of the profile line elevation function at each position parameter are read. All position parameters on the profile line are traversed sequentially. The elevation value corresponding to each position parameter is compared with the intermediate elevation. The difference between the elevation value and the intermediate elevation is calculated. Position parameters with smaller differences and at the same height level as the intermediate elevation are grouped into the same elevation interval. This forms an intermediate elevation level corresponding to the intermediate elevation, so that the intermediate elevation level can represent the set of spatial positions on the profile line within the intermediate height range between the top elevation of the dike and the average elevation of the backwater.

[0112] Specifically, the intermediate elevation of the profile line and the elevation values ​​of the profile line elevation function at each position parameter are read. All position parameters on the profile line are traversed sequentially. The elevation value corresponding to each position parameter is compared with the intermediate elevation. The difference between the elevation value and the intermediate elevation is calculated. Position parameters with smaller differences and at the same height level as the intermediate elevation are grouped into the same elevation interval. This forms an intermediate elevation level corresponding to the intermediate elevation, so that the intermediate elevation level can represent the set of spatial positions on the profile line within the intermediate height range between the top elevation of the dike and the average elevation of the backwater.

[0113] Based on the intermediate elevation, all positional parameters of the profile line are filtered to obtain a set of candidate positional parameters;

[0114] Specifically, using the location parameters corresponding to the intermediate elevation level as the screening criteria, all location parameters of the profile line are judged one by one. The location parameters belonging to the intermediate elevation level are selected and collected to form a candidate location parameter set. At the same time, the order of the location parameters in the candidate location parameter set is kept consistent with their spatial order on the profile line. This allows the candidate location parameter set to correspond to one or more continuous or discontinuous location sequences in the intermediate height range along the profile line direction, providing basic data for subsequent continuity analysis of adjacent locations.

[0115] Specifically, using the location parameters corresponding to the intermediate elevation level as the screening criteria, all location parameters of the profile line are judged one by one. The location parameters belonging to the intermediate elevation level are selected and collected to form a candidate location parameter set. At the same time, the order of the location parameters in the candidate location parameter set is kept consistent with their spatial order on the profile line. This allows the candidate location parameter set to correspond to one or more continuous or discontinuous location sequences in the intermediate height range along the profile line direction, providing basic data for subsequent continuity analysis of adjacent locations.

[0116] Based on the elevation function, the continuity analysis of the function values ​​corresponding to adjacent candidate position parameters in the candidate position parameter set is performed to obtain the candidate platform stage of the embankment shoulder.

[0117] Specifically, according to the order of the position parameters in the candidate position parameter set, the corresponding elevation function values ​​are extracted for any two adjacent candidate position parameters. The change in the two adjacent elevation function values ​​is calculated, and the parts with small changes and adjacent position parameters that are closely distributed on the profile line are grouped into the same continuous segment. The continuous segment is extended segment by segment along the candidate position parameter set towards the backwater side until the elevation change is obvious or the position parameters are no longer continuous, thereby obtaining one or more candidate segments. The starting position parameters and ending position parameters of each candidate segment are recorded as the abutment candidate platform stage of the profile line.

[0118] Specifically, the process of obtaining the abutment candidate stage by performing continuity analysis on the function values ​​corresponding to adjacent candidate position parameters in the candidate position parameter set based on the elevation function includes: extracting the corresponding elevation function values ​​for any two adjacent candidate position parameters according to the order of the position parameters in the candidate position parameter set; calculating the change in the elevation function values ​​of the two adjacent positions; grouping the parts with small changes and adjacent position parameters that are closely distributed on the profile line into the same continuous segment; extending the continuous segment segment by segment along the candidate position parameter set towards the backwater side until there is a significant change in elevation or the position parameters are no longer continuous, thereby obtaining one or more candidate segments; and recording the starting and ending position parameters of each candidate segment as the abutment candidate stage of the profile line.

[0119] Specifically, the intermediate elevation refers to the height value between the top elevation of the dike and the average elevation of the backwater, used to represent the intermediate height level between the dike body and the backwater ground on the profile line; the function value of the elevation function refers to the surface height value corresponding to each position parameter on the profile line, used to describe the topographic undulation of the profile along the backwater side; elevation interval division refers to classifying the elevation values ​​of the elevation function at each position parameter according to their relationship with the intermediate elevation, using the intermediate elevation as a reference, thereby determining the intermediate elevation level corresponding to the intermediate elevation; the intermediate elevation level refers to a group of positions on the profile line whose elevation values ​​are close to the intermediate elevation. The height levels corresponding to the parameters are used to reflect the spatial distribution of the profile within the intermediate height range; the location parameters refer to the distance markers extending from the starting point to the backwater side along the profile line, used to determine the spatial location of each sampling point on the profile line; the candidate location parameter set refers to the set of multiple location parameters selected within the intermediate elevation level, and the spatial locations corresponding to these location parameters are located within the intermediate height range of the profile; the continuity analysis refers to comparing the elevation function values ​​corresponding to adjacent location parameters in the candidate location parameter set to determine the segments on the profile line where the height change between adjacent positions is small and the segments are spatially continuous.

[0120] Specifically, the candidate platform stage of the embankment shoulder refers to a segment within the intermediate height range between the top elevation of the embankment and the average elevation of the backwater on a single profile line. It consists of multiple adjacent position parameters and is spatially continuous. The elevation function value corresponding to this segment changes little within the segment and exhibits a relatively gentle height continuity characteristic. It is used to characterize the spatial range where there may be inner embankment shoulder steps or platforms on the profile line. The candidate platform stage of the embankment shoulder is defined by its starting position parameter and ending position parameter, and may include statistical information on elevation values ​​within the segment as the basis for subsequent association and construction of the embankment shoulder centerline along the embankment line.

[0121] S3. Generate planar coordinates based on the candidate platform stage of the embankment shoulder, and generate the centerline of the embankment shoulder based on the planar coordinates;

[0122] In an embodiment of the present invention, generating planar coordinates based on the candidate platform stage of the embankment includes:

[0123] The starting position is extracted for the candidate platform stage of the embankment.

[0124] Specifically, the starting position extraction for the candidate abutment stage includes reading all position parameters contained in the candidate abutment stage, maintaining the arrangement order of the position parameters in the profile line direction, determining the minimum position parameter by comparing the magnitude of each position parameter, using the minimum position parameter as the starting position parameter, and using the spatial position of the profile line corresponding to the starting position parameter as the starting position, so that the starting position can characterize the starting point of the candidate abutment stage on the profile line.

[0125] The termination position is extracted from the candidate platform stage of the embankment.

[0126] Specifically, the process of extracting the termination position of the candidate abutment stage includes reading all the position parameters contained in the candidate abutment stage, determining the maximum position parameter by comparing the magnitude of each position parameter, using the maximum position parameter as the termination position parameter, and using the spatial position of the profile line corresponding to the termination position parameter as the termination position, so that the termination position can characterize the end point of the candidate abutment stage on the profile line.

[0127] The midpoint position is obtained by calculating the midpoint between the starting and ending positions.

[0128] Specifically, the starting position parameter corresponding to the starting position and the ending position parameter corresponding to the ending position are read. The starting position parameter and the ending position parameter are summed and divided by two to obtain the midpoint position parameter. The spatial position corresponding to the midpoint position parameter on the profile line is determined as the midpoint position, so that the midpoint position can serve as the representative position of the candidate platform stage of the embankment on the profile line.

[0129] Perform a planar coordinate transformation on the midpoint position to obtain the planar coordinates.

[0130] Specifically, the starting plane coordinates of the profile line and the direction vector information of the profile line are obtained. The starting point of the profile line is located at the mileage point and corresponds to the spatial position of the dike centerline. Based on the distance along the profile line represented by the midpoint position parameter, the distance is displaced from the starting point of the profile line along the direction of the profile line to obtain the two-dimensional position of the midpoint position in geographic space, thereby obtaining the plane coordinates corresponding to the midpoint position. These coordinates are used to express the spatial position of the dike shoulder candidate platform stage in plane space and to provide a coordinate basis for the subsequent construction of dike shoulder candidate points.

[0131] Specifically, the starting position refers to the initial position of the candidate abutment stage along the position parameter direction on the profile line, corresponding to the spatial point at the beginning of the segment; the ending position refers to the ending position of the candidate abutment stage along the position parameter direction on the profile line, corresponding to the spatial point at the end of the segment; the midpoint position refers to the center position of the segment determined by the starting and ending positions, located in the middle of the segment's length direction and used to represent the typical position of the segment on the profile line; the planar coordinate transformation refers to mapping the position of the midpoint position in the profile line coordinate expression to a planar position in geographic space, so that the midpoint position can be represented by a unified planar coordinate; the planar coordinate refers to the two-dimensional position representation of the midpoint position in geographic space, used to express the spatial position corresponding to the candidate abutment stage on a planar map and to provide a coordinate basis for subsequent association and connection along the embankment line direction.

[0132] In an embodiment of the present invention, generating the centerline of the embankment shoulder based on planar coordinates includes:

[0133] The average elevation of the embankment shoulder is obtained by averaging the elevation values ​​corresponding to all position parameters during the candidate platform stage of the embankment shoulder.

[0134] Specifically, the starting and ending position parameters of the candidate abutment stage are read, and the range of position parameters covered by the candidate abutment stage is determined accordingly. All position parameters within this range are traversed sequentially, and the corresponding elevation value is obtained from the elevation function of the profile line for each position parameter. The obtained elevation values ​​are summed one by one, and the number of elevation values ​​involved in the summation is counted. After the traversal is completed, the summation result is divided by the number to obtain the average elevation value. This average elevation value is determined as the average elevation of the abutment in the candidate abutment stage, so that the average elevation of the abutment can characterize the overall height level of the candidate abutment stage in the profile direction.

[0135] Use plane coordinates and average elevation of the embankment shoulder as attribute information;

[0136] Specifically, the use of planar coordinates and average elevation of the embankment shoulder as attribute information includes obtaining the planar coordinates corresponding to the midpoint of the candidate embankment stage, and associating the planar coordinates with the aforementioned average elevation of the embankment shoulder to form attribute information describing the spatial location and height characteristics of the candidate embankment stage. The planar coordinates are used to determine the location of the candidate embankment stage in the geographic plane, and the average elevation of the embankment shoulder is used to determine the representative value of the candidate embankment stage in the height dimension, so that the attribute information can simultaneously express the planar location and height level of the candidate embankment stage.

[0137] Based on attribute information, spatial points corresponding to planar coordinates are marked to obtain candidate points for the embankment shoulder;

[0138] Specifically, the spatial location indicated by the plane coordinates is taken as a spatial point, and the attribute information is written into the record content of the spatial point, so that the spatial point carries the average elevation of the embankment shoulder and retains its plane position identifier, thereby marking the spatial point as the spatial point corresponding to the embankment shoulder candidate stage. The marked spatial point is output as the embankment shoulder candidate point, which is used to associate and connect the embankment shoulder candidate points on adjacent profiles along the embankment direction.

[0139] Specifically, the elevation values ​​corresponding to all position parameters in the candidate platform stage refer to the set of surface height values ​​corresponding to each sampling position along the profile line within the candidate platform stage. This set reflects the spatial height distribution of the candidate platform stage. The average elevation of the embankment refers to the height value obtained by averaging the set of elevation values, used to characterize the overall height level of the candidate platform stage in the profile direction. Attribute information refers to the set of descriptive information associated with spatial points corresponding to planar coordinates, including at least planar coordinates and the average elevation of the embankment, used to express the spatial location and height characteristics of the spatial point. Marking refers to the process of attaching the attribute information to the spatial point corresponding to the planar coordinates, so that the spatial point is identified as a point with embankment characteristics in subsequent calculations. The candidate embankment point refers to the spatial point after the marking process is completed. This spatial point has a clear planar location and carries the average elevation attribute of the embankment, used as the basic element for subsequently associating adjacent candidate points along the embankment centerline and constructing the embankment centerline.

[0140] Sort all the profile lines to obtain an ordered sequence of profiles;

[0141] Specifically, each profile line is assigned a corresponding mileage point, and the mileage points are sorted according to their position on the dike centerline. All mileage points are traversed sequentially from the start to the end of the dike centerline. The profile lines corresponding to each mileage point are arranged in the traversal order to form a sequential list of profile lines. This sequential list is then defined as the ordered sequence of profiles, which reflects the spatial sequence of profile lines in the direction of dike extension and provides an ordered reference for establishing the correspondence between candidate dike shoulder points between adjacent profiles.

[0142] Adjacent association processing is performed on candidate points of the embankment shoulder based on the ordered sequence of the profile to obtain a sequence of candidate points of adjacent embankment shoulders;

[0143] Specifically, two adjacent profile lines are selected sequentially along the ordered profile sequence. The planar coordinates and average elevation of the marked candidate shoulder points on the two profile lines are read respectively. The relative proximity of the candidate points between the two profile lines in planar position is calculated, and the difference in the average elevation of the candidate shoulders is compared. In this way, the adjacent candidate points corresponding to the candidate points of the previous profile line are determined on the next profile line. The above process is repeated for each pair of adjacent profiles in the ordered profile sequence, so that the candidate shoulder points on each profile form a continuous pairing relationship with the corresponding candidate points on the next profile, thereby obtaining a sequence of adjacent candidate shoulder points arranged in the profile order. The sequence of adjacent candidate shoulder points is used to characterize the set of discrete corresponding points of the same shoulder structure along the direction of the embankment.

[0144] Connect the candidate point sequences of adjacent embankment shoulders to obtain the centerline of the embankment shoulders.

[0145] Specifically, according to the order of the adjacent shoulder candidate point sequence, two adjacent shoulder candidate points are selected in sequence, and their planar coordinates are connected by line segments to form a linear element. The connection operation is repeated for all adjacent point pairs in the sequence, so that multiple line segments are connected end to end to form a continuous broken line. The continuous broken line is determined as the shoulder centerline. After the shoulder centerline is formed, the average elevation of the shoulder carried by each shoulder candidate point in the sequence is stored as the height attribute information corresponding to the shoulder centerline, so that the shoulder centerline has both planar position expression and height feature expression, which is used for subsequent construction of the shoulder area and updating of state-owned assets measurement.

[0146] Specifically, sorting all profile lines refers to arranging the profile lines according to the sequential positions of their corresponding mileage points on the dike centerline, so that the profile lines form a continuous sequence along the dike's extension direction; the ordered profile sequence refers to the arrangement of profile lines obtained after the sorting, which reflects the spatial sequence of profile lines along the dike's direction; the adjacent association processing refers to pairing candidate points of the dike shoulder on adjacent profiles according to the adjacent profile relationship expressed by the ordered profile sequence, so that candidate points representing the same dike shoulder structure in adjacent profiles form a continuous spatial correspondence; the adjacent dike shoulder candidate point sequence refers to the set of dike shoulder candidate points arranged in profile order obtained through adjacent association processing, which is continuously distributed in space along the dike direction; and the connection refers to connecting the candidate points in the adjacent dike shoulder candidate point sequence in their order as line segments to form linear elements.

[0147] Specifically, the shoulder centerline refers to a linear spatial element formed by connecting candidate shoulder points on adjacent cross-sections along the dike's direction in spatial order. This linear element is continuously distributed in planar space and reflects the overall orientation and position of the shoulder structure on the backwater side of the dike. Each candidate shoulder point represents a typical position of the candidate shoulder platform stage on the corresponding cross-section. By connecting these typical positions in the direction of dike extension, the shoulder centerline transforms the discrete expression of the shoulder structure at the cross-sectional scale into a continuous expression along the dike line, which is used for subsequent calculation of the shoulder's occupied area and for carrying out the measurement and registration of state-owned assets.

[0148] S4. Calculate the average shoulder width based on the centerline of the shoulder, and calculate the total land area occupied by the hidden inner shoulder in the flood control dike based on the average shoulder width.

[0149] In an embodiment of the present invention, calculating the average shoulder width based on the shoulder centerline includes:

[0150] The point sequence along the centerline of the embankment shoulder is paired up to obtain multiple candidate point pairs for adjacent embankment shoulders.

[0151] The adjacent candidate points of the embankment shoulder are divided into line segments to obtain multiple smaller segments;

[0152] Specifically, the point sequence of the embankment centerline refers to an ordered set of points formed by arranging multiple embankment candidate points in sequence according to the embankment extension direction, where each point corresponds to a representative spatial position of the embankment structure at a certain cross-sectional location; the adjacent pairing process refers to pairing two adjacent points in the point sequence as a group, so that each pair of paired points corresponds to a continuous connection unit on the embankment centerline; the adjacent embankment candidate point pair refers to the combination of two adjacent embankment candidate points obtained by the adjacent pairing process, which includes the previous candidate point and the next candidate point on the embankment centerline, used to characterize the spatial connection relationship of the embankment centerline within a local range; the line segment division process refers to combining the endpoints of each adjacent embankment candidate point pair as a straight line connecting line segment, thereby forming a set of line segments composed of multiple endpoint connecting line segments; the small segment refers to a single line segment determined by the adjacent embankment candidate point pair, which corresponds to a local length of the embankment centerline in planar space.

[0153] Obtain the first shoulder width of the candidate shoulder point corresponding to one of the endpoints of a small segment of adjacent shoulder candidate points;

[0154] Specifically, obtaining the first shoulder width of a candidate point corresponding to one endpoint of a short segment of an adjacent shoulder candidate point pair includes reading the adjacent shoulder candidate point pair of the short segment, determining the first endpoint of the adjacent shoulder candidate point pair, extracting the shoulder width attribute from the attribute record of the candidate point corresponding to the first endpoint, and using the extracted shoulder width attribute as the first shoulder width. The shoulder width attribute is derived from the distance calculation result between the start position parameter and the end position parameter of the candidate shoulder platform stage corresponding to the candidate shoulder candidate point on the profile line.

[0155] Obtain the second shoulder width of the shoulder candidate point corresponding to the other endpoint in the adjacent shoulder candidate point pair;

[0156] Specifically, obtaining the second shoulder width of the shoulder candidate point corresponding to the other endpoint in the adjacent shoulder candidate point pair includes determining a second endpoint that is different from the first endpoint in the adjacent shoulder candidate point pair, extracting the shoulder width attribute from the attribute record of the shoulder candidate point corresponding to the second endpoint, and using the extracted shoulder width attribute as the second shoulder width, so that the second shoulder width characterizes the lateral scale of the shoulder candidate platform stage at the other endpoint of the segment in the profile direction.

[0157] The average shoulder width is obtained by averaging the widths of the first and second shoulders.

[0158] Specifically, the widths of the first and second embankments are read, and the two are added together to obtain the sum of the widths. The sum of the widths is then divided by two to obtain the average embankment width. This average embankment width is used as the representative value of the embankment width corresponding to the small segment, and is used together with the length of the small segment to participate in the segmented calculation of the embankment shoulder area.

[0159] Specifically, an adjacent shoulder candidate point segment refers to a single line segment formed by connecting two adjacent shoulder candidate points on the centerline of the shoulder, which corresponds to a local area of ​​the centerline of the shoulder in planar space; an endpoint refers to any one of the two shoulder candidate points constituting the segment, which corresponds to a representative position of the shoulder structure at a certain cross-sectional location in space; the shoulder width of a shoulder candidate point refers to the lateral scale of the shoulder candidate stage corresponding to that shoulder candidate point in the cross-sectional direction, which reflects the spatial span of the candidate stage from the starting position to the ending position in the cross-sectional line; the first shoulder width refers to the shoulder width attribute carried by the shoulder candidate point corresponding to one of the endpoints of the segment, which is used to characterize the lateral scale of the shoulder structure at that endpoint; the second shoulder width refers to the shoulder width attribute carried by the shoulder candidate point corresponding to the other endpoint of the segment, which is used to characterize the lateral scale of the shoulder structure at the other endpoint; the average shoulder width refers to the width value obtained by averaging the first shoulder width and the second shoulder width, which is used to characterize the representative value of the lateral scale of the shoulder structure within the segment.

[0160] In embodiments of the present invention, the total land area occupied by the concealed inner shoulder of the flood control dike is calculated based on the average shoulder width, including:

[0161] The distance between each pair of adjacent candidate points on the embankment shoulder is calculated to obtain the length of the segment.

[0162] Specifically, the planar coordinates of each adjacent candidate point of the embankment shoulder to the candidate points at both ends are read sequentially, and the planar coordinates at both ends are used as the input of the two point positions. The straight-line distance between the two points in the planar space is calculated, and the straight-line distance is recorded as the length of the corresponding small segment of the adjacent candidate point of the embankment shoulder. The length of each small segment can characterize the spatial extension scale of the embankment shoulder centerline in this local area, and the lengths of all small segments are formed into a length sequence in point-to-point order for subsequent area calculation.

[0163] Multiply the length of the short segment by the average width of the embankment shoulder to obtain the area occupied by the embankment shoulder.

[0164] Specifically, the average shoulder width and length of each segment are read, and the length and average shoulder width are multiplied to obtain the shoulder area of ​​the segment. The shoulder area is then recorded in relation to the segment number, so that each segment has a traceable shoulder area result, thus forming an area set composed of the shoulder areas of multiple segments.

[0165] The total area occupied by the hidden inner shoulder is obtained by summing up the areas occupied by all the shoulders.

[0166] Specifically, the area occupied by each shoulder in the area set is summed up to obtain the summed area of ​​the shoulder, and the summed result is determined as the total area occupied by the hidden inner shoulder, so that the total area can represent the overall planar area covered by the hidden inner shoulder structure on the back side of the dike, and serve as the area input data for subsequent measurement and updating of state-owned assets of flood control dikes.

[0167] Specifically, an adjacent shoulder candidate point pair refers to a pair of two sequentially adjacent shoulder candidate points in the shoulder centerline point sequence, which corresponds to a local connecting unit on the shoulder centerline; distance calculation refers to obtaining the straight-line distance between the two points based on the planar coordinates of the two endpoints of the adjacent shoulder candidate point pair, thereby obtaining the length of the local connecting unit in planar space; the length of a segment refers to the straight-line distance between the two endpoints of the adjacent shoulder candidate point pair, which is used to characterize the spatial extension scale of the shoulder centerline within this local area; the shoulder land area refers to the area value obtained by multiplying the length of the segment by the average shoulder width, which is used to characterize the land area occupied by the shoulder structure corresponding to the segment in planar space; summation refers to the calculation process of adding up the shoulder land areas corresponding to multiple segments to obtain the total area; the total land area of ​​the hidden inner shoulder refers to the area value obtained by summing the shoulder land areas of all segments, which is used to characterize the overall planar land area covered by the hidden inner shoulder structure on the backwater side of the dike.

[0168] S5. Update the monitoring and measurement of flood control dike assets based on the total land area of ​​the hidden inner shoulder.

[0169] Specifically, the process begins by retrieving data from the existing flood control dike asset layer, specifically the dike segment number information, terrain and digital elevation model layer corresponding to the hidden inner dike shoulder area's geometric description, and the state-owned asset ledger layer. Simultaneously, the total area and total volume of the hidden inner dike shoulder, previously obtained through geographic information layering data processing, are extracted. Then, using cross-layer data association technology within the geographic information layering framework, the spatial information of the hidden inner dike shoulder's area is precisely matched with the spatial location information of the existing dike segment, ensuring accurate spatial correspondence within the geographic information layering framework. Subsequently, in accordance with the compliance management requirements for state-owned asset supervision, the source legality of the measurement data, including the total area, total volume, incremental area, and incremental volume, is verified, the calculation process is traced, and the spatial range is compared for consistency. After successful verification, an independent asset supervision object for the hidden inner dike shoulder is established in the state-owned asset management system based on the spatial indexing mechanism of geographic information layering. Its core attribute information, such as spatial coordinates, geometric morphology, total area, and total volume, is fully entered and associated with the corresponding data in the geographic information layering framework. The layer identification and data source information clearly define the state-owned asset attributes, ownership, and management responsibility boundaries of the asset under supervision. Then, the key measurement indicators such as the total area and total volume of the corresponding dike section in the geographic information layer are updated. This ensures real-time synchronization between the ledger layer data and the spatial data in the geographic information layer representing the actual asset status of the dike. A dedicated asset registration document is generated and archived in conjunction with the geographic information layer data. Simultaneously, all updated asset data is synchronized to the state-owned asset supervision platform, achieving two-way data linkage between the supervision platform and the geographic information layer system. Subsequently, in accordance with the standardized management of state-owned assets and the dynamic update cycle of the geographic information layer data, the spatial location of the hidden inner dike shoulder assets is periodically reviewed using geographic information layering technology to verify the measurement data. Based on the updated asset data, targeted maintenance plans are developed, and the asset usage and maintenance status are tracked and recorded. Through the technical support of geographic information layering combined with the standardized management of state-owned assets, it is ensured that this portion of state-owned assets, which might otherwise be in a regulatory blind spot, remains in a clear, traceable, manageable, and controllable regulatory state.

[0170] The above description is only a preferred embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any equivalent substitutions or modifications made by those skilled in the art within the scope of the technology disclosed in the present invention, based on the technical solution and inventive concept of the present invention, should be covered within the scope of protection of the present invention.

Claims

1. A state-owned asset supervision method based on geographic information layering, characterized in that, Includes the following steps: S1. Generate a digital elevation model based on the centerline of the flood control embankment, and generate an elevation function based on the digital elevation model; The specific steps for generating an elevation function based on a digital elevation model are as follows: The centerline of the dike was divided into segments, resulting in multiple mileage points; Determine the tangent direction of the dike's centerline at the kilometer point; Construct profile lines at each mileage point in the normal direction of the tangent and facing the backwater side; Multiple positional parameters are obtained by sampling the profile lines at equal intervals. Based on the digital elevation model, the elevation of the location parameters is queried to obtain the elevation value of each location parameter; Model the functional relationship between the location parameters and their elevation values ​​to obtain the elevation function of the profile line; S2. Calculate intermediate elevation based on elevation function, and generate candidate abutment platforms based on intermediate elevation. The specific steps for calculating intermediate elevations based on elevation functions are as follows: Based on the elevation function, the elevation values ​​corresponding to the position parameters on the side closest to the center line of the dike are extracted to obtain the first set of elevation values; The maximum value of the first set of elevation values ​​is calculated to obtain the top elevation of the embankment of the profile line; Based on the elevation function, the elevation values ​​corresponding to the position parameters on the side away from the center line of the dike are extracted to obtain the second set of elevation values; The average value of the second set of elevation values ​​is calculated to obtain the backwater average elevation of the profile line. The average elevation of the dike crest and the average elevation of the backwater are averaged to obtain the intermediate elevation. S3. Generate planar coordinates based on the candidate platform stage of the embankment shoulder, and generate the centerline of the embankment shoulder based on the planar coordinates; S4. Calculate the average shoulder width based on the centerline of the shoulder, and calculate the total land area occupied by the hidden inner shoulder in the flood control dike based on the average shoulder width. S5. Update the monitoring and measurement of flood control dike assets based on the total land area of ​​the hidden inner shoulder.

2. The state-owned asset supervision method based on geographic information layering according to claim 1, characterized in that, A digital elevation model is generated based on the centerline of the flood control dike, including: Obtain the centerline of the flood control dike; Locate the adjacent surface area on the backwater side of the dike's centerline; Obtain surveying data of the dike centerline and adjacent surface areas; Elevation data was collected from the spatial area corresponding to the centerline of the dike based on the survey data to obtain elevation sampling data. Spatial interpolation and rasterization are performed on the elevation sampling data to obtain a digital elevation model.

3. The state-owned asset supervision method based on geographic information layering according to claim 1, characterized in that, The stage of generating candidate abutment platforms based on intermediate elevations includes: The elevation function values ​​are divided into elevation intervals based on the intermediate elevation to obtain the intermediate elevation level. Based on the intermediate elevation, all positional parameters of the profile line are filtered to obtain a set of candidate positional parameters; Based on the elevation function, the continuity analysis of the function values ​​corresponding to adjacent candidate position parameters in the candidate position parameter set is performed to obtain the candidate platform stage of the embankment shoulder.

4. The state-owned asset supervision method based on geographic information layering according to claim 1, characterized in that, Based on the candidate platform stage of the embankment, the planar coordinates are generated, including: The starting position is extracted for the candidate platform stage of the embankment. The termination position is extracted from the candidate platform stage of the embankment. The midpoint position is obtained by calculating the midpoint between the starting and ending positions. Perform a planar coordinate transformation on the midpoint position to obtain the planar coordinates.

5. A method for supervising state-owned assets based on geographic information layering according to claim 1, characterized in that, Generating the centerline of the embankment shoulder based on planar coordinates includes: The average elevation of the embankment shoulder is obtained by averaging the elevation values ​​corresponding to all position parameters during the candidate platform stage of the embankment shoulder. Use plane coordinates and average elevation of the embankment shoulder as attribute information; Based on attribute information, spatial points corresponding to planar coordinates are marked to obtain candidate points for the embankment shoulder; Sort all the profile lines to obtain an ordered sequence of profiles; Adjacent association processing is performed on candidate points of the embankment shoulder based on the ordered sequence of the profile to obtain a sequence of candidate points of adjacent embankment shoulders; Connect the candidate point sequences of adjacent embankment shoulders to obtain the centerline of the embankment shoulders.

6. A method for supervising state-owned assets based on geographic information layering as described in claim 5, characterized in that, The average shoulder width is calculated based on the centerline of the shoulder, including: The point sequence along the centerline of the embankment shoulder is paired up to obtain multiple candidate point pairs for adjacent embankment shoulders. The adjacent candidate points of the embankment shoulder are divided into line segments to obtain multiple smaller segments; Obtain the first shoulder width of the candidate shoulder point corresponding to one of the endpoints of a small segment of adjacent shoulder candidate points; Obtain the second shoulder width of the shoulder candidate point corresponding to the other endpoint in the adjacent shoulder candidate point pair; The average shoulder width is obtained by averaging the widths of the first and second shoulders.

7. A method for supervising state-owned assets based on geographic information layering according to claim 6, characterized in that, The total land area occupied by the concealed inner shoulder of the flood control dike is calculated based on the average shoulder width, including: The distance between each pair of adjacent candidate points on the embankment shoulder is calculated to obtain the length of the segment. Multiply the length of the short segment by the average width of the embankment shoulder to obtain the area occupied by the embankment shoulder. The total area occupied by the hidden inner shoulder is obtained by summing up the areas occupied by all the shoulders.

Images