CN103927418A - Method for manufacturing meshed drainage channels of urban road drains on basis of DEM (digital elevation model) - Google Patents
Method for manufacturing meshed drainage channels of urban road drains on basis of DEM (digital elevation model) Download PDFInfo
- Publication number
- CN103927418A CN103927418A CN201410149299.7A CN201410149299A CN103927418A CN 103927418 A CN103927418 A CN 103927418A CN 201410149299 A CN201410149299 A CN 201410149299A CN 103927418 A CN103927418 A CN 103927418A
- Authority
- CN
- China
- Prior art keywords
- road
- node
- point
- urban
- dem
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Processing Or Creating Images (AREA)
- Road Paving Structures (AREA)
Abstract
The invention relates to a method for manufacturing meshed drainage channels of urban road drains on basis of a DEM (digital elevation model) and belongs to the crossing field of municipal engineering information technology and geographic information system technology. Aiming at the need of establishing a waterlogging simulation mode and on the basis of a high-precision DEM mode of urban roads, a vector surface layer area of the urban roads is divided into high-precision regular meshes, and node vector image layers are generated at joints; ground elevation values are extracted via the DEM of the roads, nodes are selected along middle lines of the roads, and an equation of road vertical section curves is fitted; slopes of the road vertical section curves are computed and fitted, and terrain feature points of road vertical sections are selected according to a set threshold value; road drain mesh sections are divided according to the terrain feature points; a terrain feature point of each road section is selected according to a terrain feather extracting algorithm, and finally construction of the meshed drainage channels of urban road drains is realized. The method is capable of quickly and accurately constructing one-dimensional drainage channels matched with road terrain changes, and implementation of urban waterlogging preventing adjustment work is greatly simplified while a requirement on hydraulic model simulation precision is met.
Description
Technical field
The present invention relates to a kind of urban road canal networking drainage channel method for making based on DEM, belong to municipal works infotech and geographic information system technology crossing domain.
Background technology
In recent years, the waterlogging disaster of urban area takes place frequently, and the waterlogging Numerical Simulation technology based on Geographic Information System and hydraulic model is being widely used aspect urban surface waterflooding Scene Simulation and drainage waterlogging prevention planning.The site engineering of City Terrain is the main part that city prevention waterlogging Engineering System is built, by landform, vertically adjust, make the effect of urban road performance drainage channel, the rainwash that guides the rainfall that exceeds standard to form enters downstream water body or storing facilities of water is one of conventional engineering measure of prevention waterlogging Engineering System.
At present, traditional urban road is to be all vertically described with landform grid digital elevation model (DEM), but the expression of planar rasterizing is difficult to be applied to its cross dimensions and vertical planning and adjustment, when the size of each fine setting road section or vertical elevation, if adopting grid DEM model expresses, all need to reappear extraction landform altitude point and carry out space interpolation calculating, this process is loaded down with trivial details and complicated, is unfavorable for planning efficiently carrying out of adjustment work.
The effect that urban road is brought into play in prevention waterlogging Engineering System is mainly the drainage channel that serves as canal networking, according to the terrain feature of real road, same section of road has identical cross-sectional form and different profile grade variations, and the two dimension flow motion on its earth's surface can be approximately one dimension Open Channels.Therefore, two-dimentional road digital elevation model is reduced to one dimension channel, is meeting under the requirement of hydraulic model simulation precision, will greatly simplify carrying out of city prevention waterlogging planning adjustment work.How according to high-precision road DEM model, construct quickly and accurately that to change with its landform the one dimension channel conforming to will be the content that this patent emphasis is set forth.
The demand that this patent builds from waterlogging realistic model, the high accuracy DEM model based on urban road, is divided into urban road vector surface layer region high-precision anyhow to regular grid, and in regular grid intersection point place, generates knot vector figure layer; By road DEM, extract node ground elevation value, along center line of road, select node matching profile of road curvilinear equation; Digital simulation profile of road slope of a curve, selects profile of road feature topographic(al) point according to setting threshold values; According to profile of road feature topographic(al) point, divide road canal net section; According to Extraction of Terrain Features, select the unique point of each section of road, finally realize the structure of one dimension road canal networking drainage channel.
Summary of the invention
The present invention is mainly and solves the current road DEM model based on traditional and carry out city prevention waterlogging planning to adjust work difficulty high, and the shortcoming of complex operation complexity, has proposed a kind of disposal route that goes out canal networking drainage channel based on urban road high accuracy DEM model construction.
Technical scheme of the present invention is as follows:
Urban road canal networking drainage channel method for making based on DEM, described method comprises following steps successively:
(1) based on urban road vector polygon surface layer, the regular grid that dividing precision is 0.1m～1m.Rule-based grid, builds knot vector figure layer in intersection point place, node diagram layer object need comprise attribute field: ground elevation, X coordinate and Y coordinate, and above-mentioned field is all set to double precision datum type;
(2) based on road high accuracy number elevation model, extract node diagram layer object ground elevation value, extracting method is the space X Y coordinate based on being sampled node, uses the elevation numerical value of bilinear interpolation interpolation node;
(3) matching of profile of road curve;
(3.1) based on road vectors polygon diagram layer, connecting joint point diagram layer object, the line layer VectorLayer of drafting road axis, to the situation that the forms of road cross section changes and multiple tracks road crosses, need be divided into a plurality of road axis object vectors;
(3.2) first the matching of profile of road curve need construct along the rectangular coordinate system of road direction and vertical elevation, the height value of the spot elevation that employing Local Polynomial method connects same road axis VectorLayer object is carried out regretional analysis and curve, and matched curve equation form is n algebraic polynomial H=h (X)=a
0+ a
1x+ ... + a
nx
n;
Being chosen as of n algebraic polynomial carried out from low order to high order preferably, and the Coefficient of determination of curve of take is greater than 0.9 as preferred standard.By building the normal equation group of least square method, Matrix Calculating solves polynomial parameters a
0, a
1.。。a
nthereby, obtain characterizing the polynomial of degree n equation that profile of road landform changes;
(4) interval △ X=0.1 is set, digital simulation profile of road slope of a curve h ' (X), set threshold values TR, the size of TR value derives from urban road Longitudinal Slope Design standard, its value is set to 0.01, according to threshold values, select profile of road feature topographic(al) point, vertical section feature topographic(al) point meets | h ' (X) | >=TR or h ' (X)=0;
(5) according to the profile of road feature topographic(al) point of selecting, make road axis vertical line numbering as the division of canal net section, based on canal net section, node diagram layer is carried out to Local Search, under node diagram layer is carried out, the grouping of section is divided;
(6), to same canal net section node, according to Douglas-Peucker algorithm, select the unique point of road Cross Section;
For arbitrary canal net section node, carry out following recursive calculation process:
(6.1) canal net section node is calculated to any two euclidean distance between node pair, two nodes of ultimate range, as P and Q are made as fringe node group, connect fringe node and make straight line PQ;
(6.2) calculate other nodes to the vertical range of PQ, the point of chosen distance maximum, if this point is P
i, its distance is D
i;
(6.3) according to the poor δ of determine precision limit, the setting of δ value derives from the horizontal Section Design standard of road, and it is 0.1 meter that its value is set, relatively δ and D
ivalue size, if D
ibe less than precision limit poor, feature point extraction unique point complete and that select is only fringe node P and Q, if D
ibe greater than precision limit poor, node P
iaugment as unique point, make straight-line segment PP
iwith straight-line segment P
iq, returns to respectively step 6.2, carries out recursive calculation process;
Beneficial effect:
1. adopt construction method of the present invention, can construct quickly and accurately with its landform and change the one dimension channel conforming to according to high-precision road DEM model, can meet the application demand of city waterlogging plan model.
2. conventional road grid DEM model is difficult to be applied to its cross dimensions and vertical planning and adjustment, this patent is reduced to one dimension channel by two-dimentional road digital elevation model, meeting under the requirement of hydraulic model simulation precision, will greatly simplify carrying out of city prevention waterlogging planning adjustment work.
3. this patent passes through Extraction of Terrain Features, farthest retaining on the basis of former terrain feature, reduce required data volume, and carried out spatial description and the analog computation to it of road drainage passage without expending great storage resources and computational resource.
Accompanying drawing explanation:
Fig. 1 is the schematic flow sheet of the present invention's " urban road canal networking drainage channel method for making based on DEM " work.
Embodiment:
Specific embodiment of the invention flow process as shown in Figure 1, comprises the following steps:
(1) under ArcGIS GIS platform, based on urban road vector polygon surface layer, to divide anyhow to high precision regular grid, grid precision setting is 0.1m～1m.Rule-based grid, builds knot vector figure layer in intersection point place, node diagram layer object need comprise attribute field: ground elevation, X coordinate and Y coordinate, and above-mentioned field is all set to double precision datum type;
(2) based on road high accuracy number elevation model, extract node diagram layer object ground elevation value, extracting method is the space X Y coordinate based on being sampled node, use the elevation numerical value of bilinear interpolation interpolation node, by " the Extract Values To Points " instrument under ArcGIS GIS platform, can realize fast;
(3) based on road vectors polygon diagram layer, connecting joint point diagram layer object, the line layer VectorLayer of drafting road axis, to the situation that the forms of road cross section changes and multiple tracks road crosses, is divided into a plurality of road axis object vectors;
(4) first the matching of profile of road curve constructs along the rectangular coordinate system of road direction and vertical elevation, the height value of the spot elevation that employing Local Polynomial method connects same road axis VectorLayer object is carried out regretional analysis and curve, and matched curve equation form is n algebraic polynomial H=h (X)=a
0+ a
1x+ ... + a
nxn;
Being chosen as of n algebraic polynomial carried out from low order to high order preferably, and the Coefficient of determination of curve of take is greater than 0.9 as preferred standard.By building the normal equation group of least square method, Matrix Calculating solves polynomial parameters a
0, a
1a
nthereby, obtain characterizing the polynomial of degree n equation that profile of road landform changes;
(5) interval △ X=0.1 is set, digital simulation profile of road slope of a curve h ' (X), set threshold values TR, the size of TR value is relevant to Road Longitudinal Slope Design standard, conventionally be set to 0.01, according to threshold values, select profile of road feature topographic(al) point, vertical section feature topographic(al) point meets | h ' (X) | >=TR or h ' (X)=0;
(6) according to the profile of road feature topographic(al) point of selecting, make road axis vertical line numbering as the division of canal net section, based on canal net section, node diagram layer is carried out to Local Search, under node diagram layer is carried out, the grouping of section is divided;
(7), to same canal net section node, according to Douglas-Peucker algorithm, select the unique point of road Cross Section;
For arbitrary canal net section node, carry out following recursive calculation process:
(7.1) canal net section node is calculated to any two euclidean distance between node pair, two nodes of ultimate range, as P and Q are made as fringe node group, connect fringe node and make straight line PQ;
(7.2) calculate other nodes to the vertical range of straight line PQ, the point of chosen distance maximum, if this point is P
i, its distance is D
i;
(7.3) according to the poor δ of determine precision limit, the selection of δ value is relevant to the horizontal Section Design standard of road, is conventionally set to 0.1 meter, relatively δ and D
ivalue size, if D
ibe less than precision limit poor, feature point extraction unique point complete and that select is only fringe node P and Q, if D
ibe greater than precision limit poor, node P
iaugment as unique point, make straight-line segment PP
iwith straight-line segment P
iq, returns to respectively step 6.2, carries out recursive calculation process.
Claims (1)
1. the urban road canal networking drainage channel method for making based on DEM, is characterized in that:
Described method comprises following steps successively:
(1) based on urban road vector polygon surface layer, the regular grid that dividing precision is 0.1m～1m; Rule-based grid, builds knot vector figure layer in intersection point place, node diagram layer object need comprise attribute field: ground elevation, X coordinate and Y coordinate, and above-mentioned field is all set to double precision datum type;
(2) based on road high accuracy number elevation model, extract node diagram layer object ground elevation value, extracting method is the space X Y coordinate based on being sampled node, uses the elevation numerical value of bilinear interpolation interpolation node;
(3) matching of profile of road curve;
(3.1) based on road vectors polygon diagram layer, connecting joint point diagram layer object, the line layer VectorLayer of drafting road axis, to the situation that the forms of road cross section changes and multiple tracks road crosses, need be divided into a plurality of road axis object vectors;
(3.2) first the matching of profile of road curve need construct along the rectangular coordinate system of road direction and vertical elevation, the height value of the spot elevation that employing Local Polynomial method connects same road axis VectorLayer object is carried out regretional analysis and curve, and matched curve equation form is n algebraic polynomial H=h (X)=a
0+ a
1x+ ... + a
nx
n;
Being chosen as of n algebraic polynomial carried out from low order to high order preferably, and the Coefficient of determination of curve of take is greater than 0.9 as preferred standard; By building the normal equation group of least square method, Matrix Calculating solves polynomial parameters a
0, a
1; ; ;a
nthereby, obtain characterizing the polynomial of degree n equation that profile of road landform changes;
(4) interval △ X=0.1 is set, digital simulation profile of road slope of a curve h ' (X), set threshold values TR, the size of TR value derives from urban road Longitudinal Slope Design standard, its value is set to 0.01, according to threshold values, select profile of road feature topographic(al) point, vertical section feature topographic(al) point meets | h ' (X) | >=TR or h ' (X)=0;
(5) according to the profile of road feature topographic(al) point of selecting, make road axis vertical line numbering as the division of canal net section, based on canal net section, node diagram layer is carried out to Local Search, under node diagram layer is carried out, the grouping of section is divided;
(6), to same canal net section node, according to Douglas-Peucker algorithm, select the unique point of road Cross Section;
For arbitrary canal net section node, carry out following recursive calculation process:
(6.1) canal net section node is calculated to any two euclidean distance between node pair, two nodes of ultimate range, as P and Q are made as fringe node group, connect fringe node and make straight line PQ;
(6.2) calculate other nodes to the vertical range of PQ, the point of chosen distance maximum, if this point is P
i, its distance is D
i;
(6.3) according to the poor δ of determine precision limit, the setting of δ value derives from the horizontal Section Design standard of road, and it is 0.1 meter that its value is set, relatively δ and D
ivalue size, if D
ibe less than precision limit poor, feature point extraction unique point complete and that select is only fringe node P and Q, if D
ibe greater than precision limit poor, node P
iaugment as unique point, make straight-line segment PP
iwith straight-line segment P
iq, returns to respectively step 6.2, carries out recursive calculation process.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410149299.7A CN103927418A (en) | 2014-04-12 | 2014-04-12 | Method for manufacturing meshed drainage channels of urban road drains on basis of DEM (digital elevation model) |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410149299.7A CN103927418A (en) | 2014-04-12 | 2014-04-12 | Method for manufacturing meshed drainage channels of urban road drains on basis of DEM (digital elevation model) |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN103927418A true CN103927418A (en) | 2014-07-16 |
Family
ID=51145636
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410149299.7A Pending CN103927418A (en) | 2014-04-12 | 2014-04-12 | Method for manufacturing meshed drainage channels of urban road drains on basis of DEM (digital elevation model) |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103927418A (en) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106503154A (en) * | 2016-10-21 | 2017-03-15 | 长江勘测规划设计研究有限责任公司 | The automation extracting method of batch river cross-section morphological data |
| CN107844636A (en) * | 2017-10-17 | 2018-03-27 | 河北工业大学 | A kind of method that double water system based on SWMM models determines depth of accumulated water |
| CN108319758A (en) * | 2017-12-29 | 2018-07-24 | 上海市政交通设计研究院有限公司 | Tunnel drainage Optimized System Design method is worn under a kind of city based on hydraulic model |
| CN110864682A (en) * | 2019-11-27 | 2020-03-06 | 成都纵横自动化技术股份有限公司 | Unmanned aerial vehicle safety return route planning method |
| CN111815117A (en) * | 2020-06-10 | 2020-10-23 | 河海大学 | Urban waterlogging tendency simulation evaluation method based on Grasshopper platform |
| CN112283595A (en) * | 2020-11-13 | 2021-01-29 | 中煤航测遥感集团有限公司 | Pipeline oil spilling path analysis method and device, computer equipment and storage medium |
| CN114140459A (en) * | 2021-12-09 | 2022-03-04 | 中铁二院工程集团有限责任公司 | Railway cross section measuring method based on original laser point cloud |
| CN114863385A (en) * | 2022-03-23 | 2022-08-05 | 禾多科技（北京）有限公司 | Road curved surface information generation method, device, equipment and computer readable medium |
| CN116341090A (en) * | 2023-05-31 | 2023-06-27 | 四川华恒升科技发展有限公司 | Automatic plotting method for farmland water conservancy channels |
-
2014
- 2014-04-12 CN CN201410149299.7A patent/CN103927418A/en active Pending
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106503154A (en) * | 2016-10-21 | 2017-03-15 | 长江勘测规划设计研究有限责任公司 | The automation extracting method of batch river cross-section morphological data |
| CN106503154B (en) * | 2016-10-21 | 2019-07-16 | 长江勘测规划设计研究有限责任公司 | The automation extracting method of batch river cross-section morphological data |
| CN107844636A (en) * | 2017-10-17 | 2018-03-27 | 河北工业大学 | A kind of method that double water system based on SWMM models determines depth of accumulated water |
| CN108319758A (en) * | 2017-12-29 | 2018-07-24 | 上海市政交通设计研究院有限公司 | Tunnel drainage Optimized System Design method is worn under a kind of city based on hydraulic model |
| CN110864682A (en) * | 2019-11-27 | 2020-03-06 | 成都纵横自动化技术股份有限公司 | Unmanned aerial vehicle safety return route planning method |
| CN110864682B (en) * | 2019-11-27 | 2023-04-07 | 成都纵横自动化技术股份有限公司 | Unmanned aerial vehicle safety return route planning method |
| CN111815117B (en) * | 2020-06-10 | 2022-08-26 | 河海大学 | Grasshopper platform-based urban waterlogging susceptibility simulation evaluation method |
| CN111815117A (en) * | 2020-06-10 | 2020-10-23 | 河海大学 | Urban waterlogging tendency simulation evaluation method based on Grasshopper platform |
| CN112283595A (en) * | 2020-11-13 | 2021-01-29 | 中煤航测遥感集团有限公司 | Pipeline oil spilling path analysis method and device, computer equipment and storage medium |
| CN114140459A (en) * | 2021-12-09 | 2022-03-04 | 中铁二院工程集团有限责任公司 | Railway cross section measuring method based on original laser point cloud |
| CN114140459B (en) * | 2021-12-09 | 2023-04-07 | 中铁二院工程集团有限责任公司 | Railway cross section measuring method based on original laser point cloud |
| CN114863385A (en) * | 2022-03-23 | 2022-08-05 | 禾多科技（北京）有限公司 | Road curved surface information generation method, device, equipment and computer readable medium |
| CN116341090A (en) * | 2023-05-31 | 2023-06-27 | 四川华恒升科技发展有限公司 | Automatic plotting method for farmland water conservancy channels |
| CN116341090B (en) * | 2023-05-31 | 2023-08-08 | 四川华恒升科技发展有限公司 | Automatic plotting method for farmland water conservancy channels |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103927418A (en) | Method for manufacturing meshed drainage channels of urban road drains on basis of DEM (digital elevation model) | |
| CN108399309B (en) | A kind of watershed partitioning method of large scale complex topographic area hydrological distribution model | |
| CN107832931B (en) | Modularized analysis method for waterlogging risk in plain water network region | |
| CN108010103B (en) | Rapid and fine generation method of complex river terrain | |
| CN103093114B (en) | A kind of distributed basin water deficit measuring method based on landform and soil characteristic | |
| CN102902844A (en) | Sub-water basin partitioning method based on DEM (Dynamic Effect Model) data with large data quantity | |
| CN101838958B (en) | Road gradient detection method | |
| CN105138722A (en) | Digital river-lake network based method for dividing water collection unit of river basin of plain river network region | |
| CN107704592A (en) | A kind of flood forecasting service construction method based on WebGIS | |
| CN102819568B (en) | Submarine topography data establishment method based on topographical sampling point positions | |
| CN103886152B (en) | Drainage basin water catchment range manufacturing method based on urban concave downward flyover area DEM | |
| CN107679498A (en) | A kind of airborne laser point cloud downtown roads recognition methods | |
| CN110232737B (en) | Urban catchment area division method | |
| CN104574512A (en) | Multi-scale DEM (digital elevation model) construction method considering topographical semantic information | |
| CN103927788A (en) | Building ground feature DEM manufacturing method based on city vertical planning | |
| CN110532641B (en) | Surface grid layered modeling method and system | |
| CN103970837A (en) | Discontinuous DEM classified manufacturing method based on urban land and vertical planning | |
| KR20230099317A (en) | Method for determining river zone using river topology information structure and boundary of inundation range, and recording medium thereof | |
| CN104750936A (en) | GIS (Geographic Information System)-based urban hydrological run-off calculating method and system | |
| CN114969944B (en) | High-precision road DEM construction method | |
| CN102930176B (en) | Nine-grid five-spline interpolation method for geoid model | |
| Naidu et al. | Geo-spatial database creation for Wazirabad canal command area | |
| Panagoulia et al. | Deciphering the floodplain inundation maps in Greece | |
| CN114840896A (en) | Waterlogging simulation analysis method based on urban road BIM | |
| KR102243453B1 (en) | Method for updating urban flood maps using river spatial topology, and recording medium thereof |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20140716 |
|
| WD01 | Invention patent application deemed withdrawn after publication |