A vertical elevation calculation method, device and calculation equipment for a megacity

Abstract

The application belongs to the technical field of city planning, and discloses a vertical elevation calculation method and device for super-large cities and a calculation equipment. The processed geographic spatial data in the research range is acquired to perform hydrological analysis and slope analysis, obtain catchment partition and slope analysis data, and determine the slope type of each catchment partition through the superposition of the catchment partition and the slope analysis. The rainwater discharge path of each catchment partition is drawn in combination with the slope type of each catchment partition, a calculation point is determined, and the current elevation value of each calculation point is extracted. The design elevation of each calculation point is determined in combination with the elevation of the rainwater discharge port and the current elevation of the calculation point, and finally the target terrain raster data of the design is generated. The city vertical planning control point is determined through the road traffic data, and the control elevation of the control point is determined. The example of the application can provide a calculation reference for the vertical planning of the whole city through the vertical planning control point, and solves the problem that the vertical planning lacks overall consideration.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of urban planning, and in particular to a vertical elevation calculation method, device and calculation equipment for a super-large city. BACKGROUND

[0002] Urban vertical planning plays an important and positive role in solving urban waterlogging problems. In recent years, urban waterlogging problems caused by frequent super-standard rainfall have caused serious waterlogging disasters in major cities, and urban vertical planning has gradually attracted widespread attention from society.

[0003] In the early stage of urbanization, vertical planning is mainly prepared separately by district. Due to the lack of overall planning and design at a higher level, the traditional vertical elevation calculation method has the shortcomings of non-uniform control element standards and insufficient overall consideration, and it is easy to cause problems of relatively low-lying local areas in later construction. Super-large cities generally have characteristics such as dense urban population, high urbanization level, and complex and diverse topography, and the vertical planning problem is more prominent, so an effective method is needed to comprehensively consider the topographic features, rainwater discharge demand, and present construction of the city, and to unify the vertical elevation calculation benchmark to solve the problem of lack of overall consideration in vertical planning. SUMMARY

[0004] Therefore, the embodiments of the present application provide a vertical elevation calculation method, device and calculation equipment for a super-large city, so that the calculated elevation value is more systematic and overall, and provides guidance and basis for urban vertical planning.

[0005] In a first aspect, the present application provides a vertical elevation calculation method for a super-large city.

[0006] The present application is realized by the following technical solutions:

[0007] A vertical elevation calculation method for a super-large city comprises:

[0008] Obtaining geographical space data of a research range, preprocessing the geographical space data to obtain preprocessed geographical space data; wherein the geographical space data comprises topographic data, river system data, road traffic data and rainwater pipe network data;

[0009] Analyzing the preprocessed geographical space data to obtain catchment partition results and slope analysis results, converting a grid map into slope vector surface data based on the slope analysis results, and analyzing the catchment partition results based on the slope vector surface data to obtain the slope type of each catchment partition;

[0010] According to the slope vector surface data, in combination with the slope type of each catchment area, the river system data and the road traffic data, the rainwater discharge path of each catchment area is drawn, the calculation point of the vertical elevation is determined based on the rainwater discharge path, and the present status elevation value of each calculation point is extracted;

[0011] According to the discharge port elevation value of the rainwater discharge port and the present status elevation value of the calculation point, the design elevation value of each calculation point along the rainwater discharge path is determined;

[0012] Based on the design elevation value, target terrain raster data is generated, the intersection points of the main roads and the secondary roads are determined through the road traffic data, the intersection points are taken as the control points of the urban vertical overall planning, and the control elevation value of the control points is determined through the target terrain raster data.

[0013] In a preferred example of the present application, the geographic spatial data of the research range can be further obtained, and the geographic spatial data is preprocessed to obtain preprocessed geographic spatial data, including:

[0014] Obtain terrain data, river system data, road traffic data and rainwater pipe network data in the same geographic coordinate system within the research range;

[0015] Erase the high bridge surface in the road traffic data and the river system surface data from the terrain data, and regenerate the terrain data in the research range by using the difference method;

[0016] Determine the starting point and the end point of the rainwater pipe network vector line from the rainwater pipe network data, wherein the starting point of the rainwater pipe network vector line is set as the discharge port of the rainwater pipe network, and the end point of the rainwater pipe network vector line is set as the starting point of the rainwater pipe network.

[0017] In a preferred example of the present application, the preprocessed geographic spatial data can be further analyzed to obtain catchment area results and slope analysis results, the raster map is converted into slope vector surface data based on the slope analysis results, and the catchment area results are analyzed based on the slope vector surface data to obtain the slope type of each catchment area, including:

[0018] The preprocessed geographic spatial data is subjected to depression filling processing, catchment analysis and flow direction analysis by using the hydrological analysis tool in Arcgis to determine the rainwater catchment area;

[0019] The preprocessed geographic spatial data is subjected to slope analysis by using the surface analysis tool in Arcgis to output the slope analysis results;

[0020] The slope analysis result is converted into integer slope data, the integer slope data is converted into slope vector surface data through a raster-to-surface tool in Arcgis, and the research range is divided into a plain area, a gentle slope area and a steep slope area based on the slope vector surface data;

[0021] The slope vector surface data of each rainwater catchment area is obtained, and the proportions of the plain area, the gentle slope area and the steep slope area of the rainwater catchment area are analyzed to obtain the slope type of the rainwater catchment area.

[0022] In a preferred example of the present application, the research range can be further divided into a plain area, a gentle slope area and a steep slope area based on the slope vector surface data, comprising:

[0023] When the slope vector surface data is located at 0-0.003, the corresponding research range is determined as the plain area;

[0024] When the slope vector surface data is located at 0.003-0.03, the corresponding research range is determined as the gentle slope area;

[0025] When the slope vector surface data is greater than 0.03, the corresponding research range is determined as the steep slope area.

[0026] In a preferred example of the present application, the calculation point of the vertical elevation can be further determined based on the rainwater discharge path, and the current elevation value of each calculation point is extracted, comprising:

[0027] A calculation point is generated every 50m along the rainwater discharge path by using the line generation point tool in Arcgis, and the starting point to the end point is numbered from small to large;

[0028] The current elevation value of each calculation point is extracted by using the sampling tool in Arcgis.

[0029] In a preferred example of the present application, the design elevation value of each calculation point along the rainwater discharge path can be further determined according to the discharge port elevation value of the rainwater discharge port and the current elevation value of the calculation point, comprising:

[0030] The calculation point of the vertical elevation is taken as the input data, and the urban space attribute is taken as the identification item to generate the vertical elevation calculation point with the spatial attribute by using the identification tool in Arcgis, the attribute field belonging to the built-up area is valued as 1, and the attribute field not belonging to the built-up area is valued as 0;

[0031] The vertical elevation calculation point with the spatial attribute is exported to an excel file, and the excel file is screened by Python;

[0032] Read the first value in the object column in the excel file, and subtract the first value from each value in the object column and assign it to the number column from 0;

[0033] Create a height difference column, multiply the value in the number column by 0.05 and assign it to the height difference column;

[0034] Create a discharge port elevation, assign the corresponding elevation value of the number column with number 0 to the discharge port elevation;

[0035] Create a calculation elevation column and a design elevation column, and determine the calculation elevation value and design elevation value of the calculation point according to the preset rules.

[0036] In a preferred example of the present application, it can be further provided that,

[0037] Determine the design elevation value of the calculation point according to the preset rules, including:

[0038] When the calculation elevation value is less than the existing elevation value, and the space attribute is 1, the design elevation value is equal to the existing elevation value;

[0039] When the calculation elevation value is less than the existing elevation value, and the space attribute is 0, the design elevation value is equal to the calculation elevation value;

[0040] When the calculation elevation value is greater than or equal to the existing elevation value, and the space attribute is 1, the design elevation value is equal to the calculation elevation value;

[0041] When the calculation elevation value is greater than or equal to the existing elevation value, and the space attribute is 0, the design elevation value is equal to the calculation elevation value.

[0042] In a preferred example of the present application, it can be further provided that,

[0043] The target terrain grid data is generated based on the design elevation value, the intersection of the main road and the secondary road is determined through the road traffic data, the intersection is taken as the control point of the urban vertical overall planning, the control elevation value of the control point is determined through the target terrain grid data, including:

[0044] The target terrain grid data is generated based on the design elevation value of the calculation point by using the triangular mesh method in Arcgis, and the filling and depression pre-processing is performed on the target terrain grid data;

[0045] The intersection of the main road and the secondary road in the road traffic data is extracted as the control point of the urban vertical planning in the research scope;

[0046] The design elevation value corresponding to the control point in the target terrain grid data after filling and depression pre-processing is extracted as the control elevation value of the control point, which is used to guide the urban vertical planning.

[0047] The second aspect of the present application provides a vertical elevation calculation device for a super large city.

[0048] The present application is realized by the following technical solutions:

[0049] A vertical elevation calculation device for a super large city is used to execute the method of the first aspect, and comprises:

[0050] A data acquisition module is configured to acquire geographic spatial data of a research scope, pre-process the geographic spatial data, and obtain pre-processed geographic spatial data; wherein the geographic spatial data comprises terrain data, river system data, road traffic data, and rainwater pipe network data.

[0051] A data analysis module is configured to analyze the pre-processed geographic spatial data, obtain catchment area results and slope analysis results, convert a raster map into slope vector surface data based on the slope analysis results, and analyze the catchment area results based on the slope vector surface data to obtain slope types of each catchment area.

[0052] A current elevation extraction module is configured to draw rainwater discharge paths of each catchment area according to the slope vector surface data, in combination with the slope types of each catchment area, the river system data, and the road traffic data, determine calculation points of the vertical elevation based on the rainwater discharge paths, and extract current elevation values of each calculation point.

[0053] A design elevation determination module is configured to determine design elevation values of each calculation point along the rainwater discharge path according to discharge port elevation values of rainwater discharge ports and the current elevation values of the calculation points.

[0054] A control elevation calculation module is configured to generate target terrain raster data based on the design elevation values, determine intersection points of main roads and secondary roads through the road traffic data, take the intersection points as control points of urban vertical overall planning, and determine control elevation values of the control points through the target terrain raster data.

[0055] The third aspect of the present application is realized by the following technical solutions:

[0056] A vertical elevation calculation device for a super large city comprises a memory, a processor, and a computer program stored in the memory and executable on the processor, and the processor executes the computer program to realize the steps of any one of the methods of the first aspect.

[0057] Compared with the prior art, the technical solutions provided in the embodiments of the present application have at least the following beneficial effects:

[0058] This application obtains geospatial data within the research scope, preprocesses the geospatial data to obtain preprocessed geospatial data, analyzes the preprocessed geospatial data to obtain water catchment zoning results and slope analysis results, converts the raster map into slope vector surface data based on the slope analysis results, and then analyzes the water catchment zoning results based on the slope vector surface data to obtain the slope type of each water catchment zone. Based on the slope vector surface data, combined with the slope type of each water catchment zone, river system data, and road traffic data, the rainwater discharge path of each water catchment zone is drawn, the calculation points for vertical elevation are determined based on the rainwater discharge path, and the current elevation value of each calculation point is extracted. Based on the discharge outlet elevation value of the rainwater discharge outlet and the current elevation value of the calculation point, the design elevation value of each calculation point along the rainwater discharge path is determined. Based on the design elevation value, target terrain raster data is generated, the intersection points of main roads and secondary roads are determined through road traffic data, the intersection points are used as control points for urban vertical planning, and the control elevation value of the control points is determined through the target terrain raster data. After obtaining the geospatial data for the research area, this application first preprocesses the data to remove areas that do not reflect the true terrain, thereby improving calculation accuracy. Through detailed analysis of terrain slope zoning and rainwater drainage paths, it can more accurately determine the direction of rainwater flow and accumulation areas, and more accurately determine the vertical elevation requirements of different areas. Compared to traditional vertical planning methods that are limited to point data, this application generates two-dimensional surface data based on the calculated design elevation values, avoiding the problem of relatively low-lying areas in local areas, and providing an accurate basis for lower-level vertical planning. Attached Figure Description

[0059] Figure 1 A flowchart illustrating a method for calculating vertical elevation in megacities, provided as an embodiment of this application;

[0060] Figure 2 A schematic diagram of the structure of a vertical elevation calculation device for megacities provided in an embodiment of this application;

[0061] Figure 3 A schematic diagram of the structure of a vertical elevation calculation device for megacities provided in another embodiment of this application;

[0062] Explanation of reference numerals in the attached figures:

[0063] Data acquisition module-01, data analysis module-02, current elevation extraction module-03, design elevation determination module-04, control elevation calculation module-05. Detailed Implementation

[0064] The specific embodiments are merely explanatory of the present application, and are not intended to limit the present application. Those skilled in the art can make modifications to the embodiments according to the needs without creative contribution, and as long as the modifications are within the scope of the claims of the present application, they are protected by the patent law.

[0065] To make the objectives, technical solutions and advantages of the embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be described clearly and completely below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are only some, but not all of the embodiments of the present application. Based on the embodiments in the present application, all other embodiments obtained by those of ordinary skill in the art without creative contribution fall within the scope of the present application.

[0066] In addition, the term "and / or" in the present application merely describes the association relationship of the associated objects, and indicates that there can be three relationships, for example, A and / or B can represent the three cases of A alone, A and B together, and B alone. In addition, the character " / " in the present application generally represents an "or" relationship between the front and rear associated objects unless otherwise specified.

[0067] The terms "first", "second", and the like in the present application are used to distinguish the same items or similar items with basically the same function, and it should be understood that there is no logical or time sequence dependence between "first", "second", and "nth", and the quantity and execution order are not limited.

[0068] The embodiments of the present application will be described in further detail below with reference to the drawings of the specification. The vertical elevation calculation method provided by the present application is applicable to the vertical planning of a super large city, where the super large city refers to a city including multiple municipal districts, counties or larger areas. As shown in the drawings, the vertical elevation calculation method for a super large city provided by the present application includes the following steps: Figure 1

[0069] S10: Obtain geographic spatial data of the research scope, pre-process the geographic spatial data, and obtain pre-processed geographic spatial data; wherein the geographic spatial data includes terrain data, river system data, road traffic data, and rainwater pipe network data.

[0070] The terrain data is terrain raster data of the research scope, the resolution of the terrain raster data is 2m x 2m, and the terrain raster data can be obtained through aerial photogrammetry, satellite remote sensing, ground laser scanning, etc.; the river system data is river system surface data composed of river system embankments under the current status; the road traffic data includes elevated road distribution surface data, bridge distribution surface data, and urban road distribution data, as well as road width and grade information; the rainwater pipe network data includes drainage pipe discharge direction data, discharge port distribution data, and discharge port water level data.​

[0071] S20: Analyzing the pre-processed geospatial data to obtain a catchment partition result and a slope analysis result, converting the raster map into a slope vector surface data based on the slope analysis result, and analyzing the catchment partition result based on the slope vector surface data to obtain a slope type of each catchment partition.

[0072] Specifically, a series of spatial analysis functions are provided in the spatial analysis toolbox of the Arcgis geographic information system, and the geographic spatial data can be effectively processed by the Arcgis to perform hydrological analysis and slope analysis. The slope analysis result obtained by the slope analysis is recorded as "SLOPE", and the catchment partition result obtained by the hydrological analysis is recorded as "JSFQ".

[0073] S30: Drawing a rainwater discharge path of each catchment partition according to the slope vector surface data, combining the slope type of each catchment partition, river system data and road traffic data, determining a calculation point of the vertical elevation based on the rainwater discharge path, and extracting a current elevation value of each calculation point.

[0074] Specifically, the slope type of each catchment partition is determined according to the slope vector surface data, different slope types reflect different characteristics of the terrain, and the calculation point of the vertical elevation is determined differently for different types of catchment partitions, and the current elevation value of each calculation point is extracted, which can ensure that the vertical design not only conforms to the natural terrain, but also can efficiently discharge water.

[0075] S40: Determining a design elevation value of each calculation point along the rainwater discharge path according to the discharge elevation value of the rainwater discharge port and the current elevation value of the calculation point.

[0076] The application considers both rainwater surface runoff and the comparison with the current elevation value when calculating the design elevation value of the calculation point of the vertical elevation. When the current elevation value can meet the minimum requirement of rainwater surface runoff, the calculation point of the vertical elevation is considered as a current limiting factor, and the current elevation value is maintained unchanged. Moreover, the application standardizes the calculation process of the calculation point of the vertical elevation, providing a basis for dynamic update of planning data.

[0077] S50: Generating target terrain raster data based on the design elevation value, determining the intersection points of the main road and the secondary road through the road traffic data, taking the intersection points as control points of the urban vertical overall planning, and determining the control elevation value of the control points through the target terrain raster data.

[0078] Compared with the traditional vertical planning method which is limited to the calculation of the vertical elevation of the road intersection, the application generates two-dimensional DEM data from the calculated vertical elevation calculation points, converts the traditional one-dimensional point data of vertical planning into two-dimensional surface data, and then performs depression filling processing to solve the problem of relative low-lying in some areas. The intersection points of the main road and the secondary road are identified as the vertical planning control points, and the elevation value is assigned to the control points, which can provide accurate calculation basis for the lower level vertical detailed vertical planning and solve the problem of inconsistent calculation basis for the lower level vertical planning.

[0079] In some embodiments, the geographic spatial data of the study area is obtained in step S10, and the geographic spatial data is preprocessed to obtain preprocessed geographic spatial data, specifically including: obtaining the shapefile data of the terrain data, river system data, road traffic data and rainwater pipe network data in the same geographic coordinate system within the study area; using the erase tool in the Arcgis analysis toolbox to erase the elevated bridge surface in the road traffic data and the river system surface data from the terrain data, and using the difference method to regenerate the geographic spatial data of the study area; combining the rainwater discharge direction, determining the starting point and the end point of the rainwater pipe network vector line from the rainwater pipe network data, wherein the starting point of the rainwater pipe network vector line is set as the discharge port of the rainwater pipe network, and the end point of the rainwater pipe network vector line is set as the starting point of the rainwater pipe network. It needs to be noted that the discharge port is the outlet of rainwater in the rainwater pipe network system. Usually, these discharge ports are located at the downstream end of the pipe network, connected to natural water bodies (such as rivers, lakes, oceans) or artificial drainage facilities (such as drainage ditches, water storage tanks); the starting point of the rainwater pipe network is the entrance of rainwater into the pipe network, which can be a rainwater inlet, a rainwater well or other collection facilities. In actual implementation process, it is found that the water system and the elevated bridge road surface cannot reflect the true ground elevation, which will affect the calculation result, and the water system and the elevated bridge road surface are erased to improve the accuracy of subsequent elevation calculation.

[0080] In some embodiments, the preprocessed geographic spatial data is analyzed in step S20 to obtain catchment partition results and slope analysis results, the raster map is converted into slope vector surface data based on the slope analysis results, and the slope type of each catchment partition is obtained by analyzing the catchment partition results based on the slope vector surface data, specifically including:

[0081] The terrain data of the preprocessed geographic spatial data is analyzed by using the hydrological analysis tool in Arcgis to fill the depression, catchment area and flow direction, determine the rainwater catchment partition, and name it as "JSFQ". The terrain data of the preprocessed geographic spatial data is analyzed by using the slope analysis in the surface analysis tool in Arcgis, and the slope analysis result with the measurement unit of PERCENT_RISE is output, named as "SLOPE". The slope analysis result is converted into integer slope data, and the integer slope data is converted into slope vector surface data by using the raster to surface tool in Arcgis. Specifically, the "SLOPE" data is multiplied by 10 times to generate "SLOPE_time10" data, the "SLOPE_time10" data is converted into integer, and "INT_SLOPE_time10" data is output. The "INT_SLOPE_time10" data is converted into surface data by using the raster to surface tool, and the field is the raster value "value". The terrain slope value after being enlarged by 10 times is obtained, and the terrain slope surface data is named as "XZDXPD", which is the present situation terrain slope. The field "SLOPE" is added in the vertical table of "XZDXPD", and the field calculator tool is used, and the formula is "SLOPE=[value]", so that the slope vector surface data with the present situation terrain slope attribute is obtained, and the unit is ‰. Based on the slope vector surface data, the research range is divided into plain area, gentle slope area and steep slope area. The slope vector surface data of each rainwater catchment partition is obtained, and the proportion of the plain area, the gentle slope area and the steep slope area of the rainwater catchment partition is analyzed to obtain the slope type of the rainwater catchment partition. Specifically, the rainwater partition result "JSFQ" (rainwater catchment partition) is taken as an input factor, the slope partition result is taken as an identification factor, the area proportion of each slope partition in each rainwater catchment partition is counted, and when the area proportion of a certain slope partition is greater than 80%, the rainwater partition can be regarded as the terrain area. For example, when 85% of the area of the rainwater catchment partition JSFQ-1 belongs to the gentle slope area, 5% belongs to the plain area, and 10% belongs to the steep slope area, the slope type of the rainwater catchment partition JSFQ-1 is the gentle slope area.

[0082] The step of dividing the research range into the plain area, the gentle slope area and the steep slope area includes: when the slope vector surface data is located at 0-0.003, the corresponding area is determined as the plain area; when the slope vector surface data is located at 0.003-0.03, the corresponding area is determined as the gentle slope area; and when the slope vector surface data is greater than 0.03, the corresponding area is determined as the steep slope area. In the actual implementation process, the area of the terrain slope located at 0-0.003, 0.003-0.03 and greater than 0.03 in the slope vector surface data is divided by using the attribute selection tool in Arcgis.

[0083] For different slope types of rainwater catchment partition, differential classification calculation can be performed, and the calculation target of vertical elevation is as follows:

[0084] For plain areas, the vertical planning should be based on the principle of meeting rainwater discharge to ensure that the overall terrain is sloped to the rainwater discharge area with a slope of not less than 0.001; for gentle slope areas, the existing terrain basically meets the rainwater discharge requirements, and the vertical planning should strengthen the connection with the existing area to avoid local low-lying areas; for steep slope areas, due to the steep terrain slope, rainwater discharge is not the main limiting factor, and the vertical planning should be based on maintaining the status quo while meeting the reasonable planning under the road design requirements.

[0085] After obtaining the slope type of the rainwater catchment partition, further combined with the river system distribution, road traffic distribution, rainwater pipe network, etc. in each rainwater catchment partition, the rainwater discharge path is drawn, named "YSJL", i.e. rainwater runoff, and numbered. At the same time, the drawing principles of rainwater discharge path of each slope type are as follows: for the rainwater catchment partition of the plain area, the terrain slope is less than 0.003, and the hydrological analysis result is difficult to reflect the true rainwater catchment, the vertical planning should be based on the purpose of meeting the reasonable rainwater discharge of the area, and the rainwater pipe network distribution should be taken as the basis to draw the rainwater discharge path vector line data in each partition from the rainwater discharge port as the starting point in the reverse drainage direction; for the rainwater catchment partition of the gentle slope area, the terrain slope is between 0.003 and 0.03, and the regional rainwater can basically meet the natural discharge condition, the vertical planning should be based on the hydrological analysis result, combined with the traffic road distribution and river system distribution, to draw the rainwater discharge path vector line data in each partition; for the rainwater catchment partition of the steep slope area, the terrain slope is greater than 0.03, and the regional rainwater discharge condition is good, basically maintaining the status quo.

[0086] In some embodiments, the step of determining the calculation points of the vertical elevation based on the rainwater discharge path in step S30 and extracting the existing elevation value of each calculation point includes:

[0087] The along-line point generation tool in Arcgis is used to generate a calculation point every 50 m along the rainwater discharge path, and the starting point to the end point is numbered from small to large. The current elevation value of each calculation point is extracted by using the sampling tool in the extraction analysis of Arcgis. In the actual implementation process, the along-line point generation tool in the sampling tool of Arcgis is used to generate vertical elevation calculation points along the rainwater discharge path at a distance of 50 m, and the end points are included. In this operation process, the points generated from the same runoff path have the same number as in the previous step, and the "OBJECTID" field in the point data attribute table is automatically sorted from small to large from the starting point to the end point. The current elevation value is extracted from the current terrain raster data according to the generated vertical elevation calculation point by using the sampling tool in the extraction analysis of Arcgis. The current elevation value is named "RASTERVALU". The "XZBG" field is created. The "RASTERVALU" value is assigned to the "XZBG" field. The X and Y fields with horizontal and vertical coordinate attributes are added, and the coordinate values are assigned to the X and Y fields respectively.

[0088] In some embodiments, the step of determining the design elevation value of each calculation point along the rainwater discharge path in step S40 according to the discharge port elevation value of the rainwater discharge port and the current elevation value of the calculation point includes: using the identification tool in Arcgis, taking the vertical elevation calculation point as the input data, and taking the urban space attribute as the identification item to generate the vertical elevation calculation point with the space attribute. The attribute field belonging to the built-up area is assigned a value of 1, and the attribute field not belonging to the built-up area is assigned a value of 0. The vertical elevation calculation point with the space attribute is exported to an excel file, and the excel file is data filtered by Python. The first value in the object column of the excel file is read, and each value in the object column is subtracted from the first value and assigned to the number column for numbering from 0. The difference column is created, and the value of the number column is multiplied by 0.05 to assign to the difference column. The discharge port elevation is created, and the elevation value corresponding to the number 0 in the number column is assigned to the discharge port elevation. The calculation elevation column and the design elevation column are created, and the design elevation value of the calculation point is determined row by row according to the preset rule.

[0089] Specifically, the identification tool in Arcgis is used, the calculation point of vertical elevation is taken as the input element, the built-up area and the urban space attribute not belonging to the built-up area are taken as the identification elements, the calculation point of vertical elevation with the urban space attribute is generated, the field "JCQ" (built-up area) and "FJCQ" (non-built-up area) are added, the "JCQ" field corresponding to the calculation point located in the built-up area is assigned a value of "1", otherwise the "JCQ" field corresponding to the calculation point located in the non-built-up area is assigned a value of "0"; the attribute table of the calculation point of vertical elevation with complete attribute information is exported to an excel file, the "pandas" table reading and processing module in Python is used, data filtering is performed according to the number of each rainwater discharge path, and a separate sheet table is created; the first number in the "OBJECTID" column of each sheet table is read, and each value in "OBJECTID" is subtracted from the first number and assigned to a new column "BH", that is, the value "0" indicates that the point is a rainwater pipe discharge port, and the number value of the calculation point is increased in the reverse direction of the rainwater discharge; a new column "GC" (i.e. height difference) is created, the value "BH" x 0.05 is assigned to the column "GC", and the column "GC" represents the minimum height difference requirement from the calculation point to the rainwater discharge point; a new column "PFKBG" (i.e. discharge port elevation) is created, the value of "XZBG" in the first row "BH" = 0 is assigned to the column; new columns "JSBG" and "SJBG" (i.e. calculation elevation and design elevation) are created, and the calculation elevation value and the design elevation value are calculated row by row according to the number sorting.

[0090] The step of determining the design elevation value of the calculation point according to the preset rule row by row includes:

[0091] When the calculation elevation value is less than the current elevation value, and the space attribute is 1, the design elevation value is equal to the current elevation value;

[0092] When the calculation elevation value is less than the current elevation value, and the space attribute is 0, the design elevation value is equal to the calculation elevation value;

[0093] When the calculation elevation value is greater than or equal to the current elevation value, and the space attribute is 1, the design elevation value is equal to the calculation elevation value;

[0094] When the calculation elevation value is greater than or equal to the current elevation value, and the space attribute is 0, the design elevation value is equal to the calculation elevation value.

[0095] 1) For the first row "BH" = 0, that is, the rainwater discharge port, generally the river road elevation or the river embankment elevation is taken, so as to maintain the current elevation, and "JSBG" = "XZBG" = "SJBG", that is, the calculation elevation value is equal to the current elevation value and the design elevation value;

[0096] 2) For "BH" = 1, the "JSBG" value is equal to the "SJBG" at "BH" = 0 in the previous row plus a height difference of 0.05 m, and the design elevation value for this row needs to be determined according to the following judgment conditions:

[0097] 1. Judgment condition one:

[0098] If "JSBG" < "XZBG" and "JCQ" = 1;

[0099] then "SJBG" = "XZBG".

[0100] 2. Judgment condition two:

[0101] If "JSBG" < "XZBG" and "JCQ" = 0;

[0102] then "SJBG" = "JSBG".

[0103] 3. Judgment condition three:

[0104] If "JSBG" >= "XZBG" and "JCQ" = 1;

[0105] then "SJBG" = "JSBG".

[0106] 4. Judgment condition four:

[0107] If "JSBG" >= "XZBG" and "JCQ" = 0;

[0108] then "SJBG" = "JSBG".

[0109] 3) For "BH" = 2, the above processing flow is looped.

[0110] In some embodiments, step S50 specifically comprises: generating target terrain raster data based on the design elevation value of the calculation point using the triangular grid method in Arcgis, and performing depression filling preprocessing on the target terrain raster data;

[0111] Extracting the intersection points of the main road and the secondary road in the road traffic data as the control points of the urban vertical planning of the research scope;

[0112] Extracting the design elevation value corresponding to the control point in the target terrain raster data after depression filling preprocessing as the control elevation value of the control point, which is used to guide the urban vertical planning.

[0113] The technical scheme provided in the application solves the problem that the original technical method is not universally applicable to complex and large cities, and at the same time, due to more limiting factors in the original technical method, it is difficult to be applicable to the calculation of large-scale vertical planning elevation of a city. The application mainly provides a vertical planning elevation calculation method capable of guaranteeing the safety bottom line of regional flood in a large-scale city. Compared with the original vertical planning method which is limited to the calculation of vertical elevation of road intersections, the application is based on hydrological analysis, considers rainwater surface runoff, current construction and other limiting factors, and finely calculates the elevation along the rainwater runoff line, converts the traditional vertical planning point data into two-dimensional surface data, and then performs depression processing, solves the problem of relative low-lying in some areas, and realizes the change from one-dimensional control to two-dimensional plane control of vertical planning. Finally, the intersection points of the main road and the secondary road are identified and extracted as the vertical planning control points, the elevation values are assigned to the control points, and calculation basis is provided for the lower-level vertical detailed planning, solving the problem of non-uniform calculation basis of lower-level vertical planning elevation.

[0114] Another embodiment of the application also provides a vertical elevation calculation device for a super large city, which is used to execute the above-mentioned calculation method, as shown in Figure 2 The calculation device comprises:

[0115] A data acquisition module 01 is configured to acquire geographic spatial data of a research range, pre-process the geographic spatial data, and obtain pre-processed geographic spatial data; wherein the geographic spatial data comprises terrain data, river system data, road traffic data and rainwater pipe network data.

[0116] A data analysis module 02 is configured to analyze the pre-processed geographic spatial data, obtain catchment partition results and slope analysis results, convert a grid map into a slope vector surface data based on the slope analysis results, and analyze the catchment partition results based on the slope vector surface data to obtain the slope type of each catchment partition.

[0117] A current elevation extraction module 03 is configured to draw a rainwater discharge path of each catchment partition according to the slope vector surface data, in combination with the slope type of each catchment partition, the river system data and the road traffic data, determine a calculation point of the vertical elevation based on the rainwater discharge path, and extract a current elevation value of each calculation point.

[0118] A design elevation determination module 04 is configured to determine a design elevation value of each calculation point along the rainwater discharge path according to a discharge port elevation value of a rainwater discharge port and the current elevation value of the calculation point.

[0119] The control elevation calculation module 05 is used to generate target terrain raster data based on the design elevation value, determine the intersection of main roads and secondary roads through road traffic data, use the intersection as the control point of the city's vertical master plan, and determine the control elevation value of the control point through the target terrain raster data.

[0120] The specific limitations of the vertical elevation calculation device for megacities provided in this embodiment can be found in the embodiment of the vertical elevation calculation method for megacities described above, and will not be repeated here. Each module in the above-described vertical elevation calculation device for megacities can be implemented entirely or partially through software, hardware, or a combination thereof. These modules can be embedded in or independent of the processor in a computer device in hardware form, or stored in the memory of a computer device in software form, so that the processor can call and execute the operations corresponding to each module.

[0121] In some embodiments, such as Figure 3 As shown in the illustration, this application also provides a vertical elevation calculation device for megacities. This device may include a processor, memory, network interface, and database connected via a system bus. The processor provides computational and control capabilities. The memory includes a non-volatile storage medium and internal memory. The non-volatile storage medium stores an operating system, computer programs, and a database. The internal memory provides an environment for the operation of the operating system and computer programs in the non-volatile storage medium. The network interface is used to communicate with external terminals via a network connection. When the computer program is executed by the processor, it causes the processor to perform the steps of the vertical elevation calculation method for megacities as described in any of the above embodiments.

[0122] The working process, working details, and technical effects of the computer equipment provided in this embodiment can be found in the embodiment of the vertical elevation calculation method for megacities described above, and will not be repeated here.

[0123] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program. The computer program can be stored in a non-volatile computer readable storage medium, and when the computer program is executed, the processes of the above-mentioned embodiment methods can be included. Any reference to memory, storage, database or other medium used in each embodiment provided in the present application can include non-volatile and / or volatile memory. Non-volatile memory can include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM) or flash memory. Volatile memory can include random access memory (RAM) or external cache memory. As an illustration but not limitation, RAM is available in various forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDR SDRAM), enhanced SDRAM (ESDRAM), synchronous link (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), etc.

[0124] The technical features of the above embodiments can be combined in any manner. In order to make the description concise, not all possible combinations of the technical features in the above embodiments are described, but as long as the combinations of the technical features do not exist contradictions, they should be considered as the scope of the present disclosure.

[0125] Those skilled in the art can clearly understand that, for the convenience and brevity of description, only the division of the above functional units and modules is exemplified, and in actual application, the above functions can be completed by different functional units and modules according to needs, that is, the internal structure of the system is divided into different functional units or modules to complete all or part of the functions described above.

Claims

1. A method for calculating vertical elevation for mega-cities, characterized by, The method comprises the following steps: acquiring geographic spatial data in a research range, preprocessing the geographic spatial data, and obtaining preprocessed geographic spatial data; wherein the geographic spatial data comprises terrain data, river system data, road traffic data, and rainwater pipe network data; analyzing the preprocessed geographic spatial data to obtain catchment area results and slope analysis results, converting a raster map into slope vector surface data based on the slope analysis results, and analyzing the catchment area results based on the slope vector surface data to obtain slope types of each catchment area; drawing rainwater discharge paths of each catchment area according to the slope vector surface data, combining the slope types of each catchment area, the river system data, and the road traffic data, determining calculation points of vertical elevations based on the rainwater discharge paths, and extracting current elevation values of each calculation point; determining design elevation values of each calculation point along the rainwater discharge paths according to discharge elevation values of rainwater discharge outlets and the current elevation values of the calculation points; generating target terrain raster data based on the design elevation values, determining intersection points of main roads and secondary roads through the road traffic data, taking the intersection points as control points of urban vertical overall planning, and determining control elevation values of the control points through the target terrain raster data.

2. The method for calculating the vertical elevation for mega cities as claimed in claim 1, wherein, The method of acquiring geographic spatial data in a research range, preprocessing the geographic spatial data, and obtaining preprocessed geographic spatial data comprises the following steps: acquiring terrain data, river system data, road traffic data, and rainwater pipe network data in the same geographic coordinate system within the research range; erasing high bridge surfaces in the road traffic data and river system surface data from the terrain data, and regenerating terrain data within the research range by using a difference method; determining starting points and ending points of rainwater pipe network vector lines from the rainwater pipe network data, wherein the starting points of the rainwater pipe network vector lines are set as rainwater discharge outlets of the rainwater pipe network, and the ending points of the rainwater pipe network vector lines are set as starting points of the rainwater pipe network.

3. The method for calculating the vertical elevation of mega cities as claimed in claim 2, wherein, The method of analyzing the preprocessed geographic spatial data to obtain catchment area results and slope analysis results, and converting a raster map into slope vector surface data based on the slope analysis results, and analyzing the catchment area results based on the slope vector surface data to obtain slope types of each catchment area comprises the following steps: using a hydrological analysis tool in Arcgis to perform depression filling processing, catchment analysis, and flow direction analysis on the preprocessed geographic spatial data to determine rainwater catchment areas; using a surface analysis tool in Arcgis to perform slope analysis on the preprocessed geographic spatial data to output slope analysis results; converting the slope analysis results into integer slope data, converting the integer slope data into slope vector surface data by using a raster-to-surface tool in Arcgis, and dividing the research range into plain areas, gentle slope areas, and steep slope areas based on the slope vector surface data; acquiring slope vector surface data of each rainwater catchment area, and analyzing proportions of the plain areas, the gentle slope areas, and the steep slope areas of the rainwater catchment areas to obtain slope types of the rainwater catchment areas.

4. The method for calculating the vertical elevation of mega-cities according to claim 3, characterized in that, The method of dividing the research range into plain areas, gentle slope areas, and steep slope areas based on the slope vector surface data comprises the following steps: When the slope vector surface data is located at 0-0.003, the corresponding research range is determined as a plain area; When the slope vector surface data is located at 0.003-0.03, the corresponding research range is determined as a gentle slope area; When the slope vector surface data is greater than 0.03, the corresponding research range is determined as a steep slope area.

5. The method for calculating the vertical elevation of mega cities as claimed in claim 1 wherein, The calculation points of the vertical elevation are determined based on the rainwater discharge path, and the current elevation values of the calculation points are extracted, including: A calculation point is generated every 50 m along the rainwater discharge path by using the line generation point tool in Arcgis, and the starting point to the end point is numbered from small to large; The current elevation value of each calculation point is extracted by using the sampling tool in the extraction analysis of Arcgis.

6. The method for calculating the vertical elevation of mega cities as claimed in claim 1 wherein, The design elevation value of each calculation point along the rainwater discharge path is determined according to the discharge port elevation value of the rainwater discharge port and the current elevation value of the calculation point, including: The calculation point of the vertical elevation is used as the input data, and the urban space attribute is used as the identification item to generate the vertical elevation calculation point with the space attribute, and the attribute field of the built-up area is assigned a value of 1, and the attribute field of the non-built-up area is assigned a value of 0; The vertical elevation calculation point with the space attribute is exported to an excel file, and the excel file is screened by Python; The first value in the object column of the excel file is read, and each value in the object column is subtracted from the first value and assigned to the number column from 0; A difference column is created, and the value of the number column is multiplied by 0.05 to assign to the difference column; The corresponding elevation value of the number column with a number of 0 is assigned to the discharge port elevation; The calculation elevation column and the design elevation column are created, and the calculation elevation value and the design elevation value of the calculation point are determined row by row according to the preset rule.

7. The method for calculating the vertical elevation of mega-cities according to claim 6, characterized in that, The design elevation value of the calculation point is determined row by row according to the preset rule, including: When the calculation elevation value is less than the current elevation value, and the space attribute is 1, the design elevation value is equal to the current elevation value; When the calculation elevation value is less than the current elevation value, and the space attribute is 0, the design elevation value is equal to the calculation elevation value; When the calculation elevation value is greater than or equal to the current elevation value, and the space attribute is 1, the design elevation value is equal to the calculation elevation value; When the calculation elevation value is greater than or equal to the current elevation value, and the space attribute is 0, the design elevation value is equal to the calculation elevation value.

8. The method for calculating the vertical elevation of mega cities as claimed in claim 1 wherein, The target terrain raster data is generated based on the design elevation value, the intersection points of the main road and the secondary road are determined through the road traffic data, the intersection points are used as the control points of the urban vertical overall planning, the control elevation value of the control point is determined through the target terrain raster data, including: The target terrain raster data is generated based on the design elevation value of the calculation point by using the triangular mesh method in Arcgis, and the target terrain raster data is preprocessed by filling the pits; The intersection points of the main road and the secondary road in the road traffic data are extracted as the control points of the urban vertical planning in the research range; The design elevation value corresponding to the control point in the target terrain raster data after the pit filling preprocessing is extracted as the control elevation value of the control point, which is used to guide the urban vertical planning.

9. A vertical elevation calculating device for mega cities, characterized in that, A computer program product for performing the method of any one of claims 1 to 8, comprising: a data acquisition module configured to acquire geospatial data of a study area, and to preprocess the geospatial data to obtain preprocessed geospatial data, wherein the geospatial data comprises terrain data, river system data, road traffic data, and rainwater pipe network data; a data analysis module configured to analyze the preprocessed geospatial data to obtain catchment area results and slope analysis results, to convert a raster map into slope vector surface data based on the slope analysis results, and to analyze the catchment area results based on the slope vector surface data to obtain slope types of the catchment areas; a current elevation extraction module configured to draw rainwater discharge paths of the catchment areas based on the slope vector surface data, the slope types of the catchment areas, the river system data, and the road traffic data, to determine calculation points of the vertical elevations based on the rainwater discharge paths, and to extract current elevation values of the calculation points; a design elevation determination module configured to determine design elevation values of the calculation points along the rainwater discharge paths based on discharge elevation values of rainwater discharge outlets and the current elevation values of the calculation points; a control elevation calculation module configured to generate target terrain raster data based on the design elevation values, to determine intersection points of main roads and secondary roads based on the road traffic data, to take the intersection points as control points of a city vertical overall plan, and to determine control elevation values of the control points based on the target terrain raster data.

10. A vertical elevation calculation device for mega cities, characterized by, A computer program product for performing the method of any one of claims 1 to 8, comprising: a data acquisition module configured to acquire geospatial data of a study area, and to preprocess the geospatial data to obtain preprocessed geospatial data, wherein the geospatial data comprises terrain data, river system data, road traffic data, and rainwater pipe network data; a data analysis module configured to analyze the preprocessed geospatial data to obtain catchment area results and slope analysis results, to convert a raster map into slope vector surface data based on the slope analysis results, and to analyze the catchment area results based on the slope vector surface data to obtain slope types of the catchment areas; a current elevation extraction module configured to draw rainwater discharge paths of the catchment areas based on the slope vector surface data, the slope types of the catchment areas, the river system data, and the road traffic data, to determine calculation points of the vertical elevations based on the rainwater discharge paths, and to extract current elevation values of the calculation points; a design elevation determination module configured to determine design elevation values of the calculation points along the rainwater discharge paths based on discharge elevation values of rainwater discharge outlets and the current elevation values of the calculation points; a control elevation calculation module configured to generate target terrain raster data based on the design elevation values, to determine intersection points of main roads and secondary roads based on the road traffic data, to take the intersection points as control points of a city vertical overall plan, and to determine control elevation values of the control points based on the target terrain raster data.

Images