CN109472868B - A kind of watershed partitioning method in continental river basin - Google Patents
A kind of watershed partitioning method in continental river basin Download PDFInfo
- Publication number
- CN109472868B CN109472868B CN201811313462.3A CN201811313462A CN109472868B CN 109472868 B CN109472868 B CN 109472868B CN 201811313462 A CN201811313462 A CN 201811313462A CN 109472868 B CN109472868 B CN 109472868B
- Authority
- CN
- China
- Prior art keywords
- waterways
- network
- grid
- basin
- virtual
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2210/00—Indexing scheme for image generation or computer graphics
- G06T2210/32—Image data format
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Remote Sensing (AREA)
- Computer Graphics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of watershed partitioning methods in continental river basin, comprising: sets in continental river basin low-lying land or lake center as its virtual exit point;Height value at exit point grid virtual in DEM is set as Null, obtains DEM1;The virtual network of waterways of vector is drawn, so that the network of waterways is finally connected to virtual exit point, and change into reference to network of waterways grid;Continental river basin inner and outer boundary is changed into grid;It is modified to elevation at the network of waterways and basin inner and outer boundary grid positions is referred in DEM1, obtains DEM2;The simulation network of waterways is extracted from DEM2;Lake range is encoded and is converted to grid;Traversal of tracing to the source is carried out to simulation network of waterways grid, and considers that lake range carries out network of waterways coding;To simulate network of waterways grid as starting point, traversal of tracing to the source is carried out to non-network of waterways grid, the network of waterways coding of the network of waterways grid is set as all non-network of waterways grid sub-basins for importing current network of waterways grid and is encoded.The present invention can continental river basin to no river outlet carry out watershed partitioning, and ensure that lake is located in the same sub-basin.
Description
Technical field
The present invention relates to hydrological model fields, and in particular to a kind of watershed partitioning method in continental river basin.
Background technique
In hydrological model, conventional extract basin and the method for sub-basin include using GIS (Geographic
Information System, GIS-Geographic Information System) software hydrological analysis module from DEM (Digital Elevation
Model, digital elevation model) start to undergo and fills out low-lying area, flow to calculating, the calculating of confluence cumulative number, simulate extraction of drainage, sub-basin stroke
Divide and etc.;The premise that this method can correctly carry out watershed partitioning, which is the need for one or more, to flow out through a series of grids
The basin on the boundary DEM exports.
For continental river basin, basin outlet is normally at basin inside center, and periphery grid is higher by one big section,
And a big continental river basin often by multiple relatively independent small interior stream district's groups at.If mentioned using above-mentioned conventional method
Basin and sub-basin are taken, is had the following problems:
(1) when carrying out filling out low-lying area, the outlet inside basin can be filled and led up, and from the lower region exit boundary in boundary,
When watershed partitioning be easy make mistake upstream and downstream confluence relationship and by the small interior stream region merging technique of multiple independences be an area
Domain;
(2) for continental river basin, mountain area water abundance, usually perennial stream, middle and lower reaches are generally risen in upstream
Typically seasonal river, least significant end are usually wetland or lower reaches lake, and moisture, which is basically used for evaporation, will not flow out outside basin;
In order to facilitate Water Balance of Lakes analysis, need for lake range to be divided into the same sub-basin;Conventional method is without in satisfaction
Requirement is stated, the problem of drawing in lake to multiple sub-basins often occurs;
(3) for the middle and lower reaches of continental river basin, since height difference is relatively small in region, the son divided according to a conventional method
With water resources etc., other have been fruitful and often differ larger basin perimeter, are not easy in model analog analysing unified to phase
With unit analyzed, and will lead to borderline region divide it is abnormal so that the non-basin region in periphery draw to current basin or
The current basin of person lacks a part of area.
Summary of the invention
The technical problem to be solved by the present invention is to the existing conventional methods for extracting basin and sub-basin, are internally flowing
When river valley carries out sub-basin stroke, merge there are sub-basin upstream and downstream confluence relationship exception, by the independent area little Nei Liu, lake divides
The problems such as to multiple sub-basins and sub-basin range partition error.
In order to solve the above-mentioned technical problem, the technical scheme adopted by the invention is that providing a kind of subflow in continental river basin
Domain classification method, comprising the following steps:
Step S10, the lake center position of the local elevation minimum point or No way out in continental river basin to be divided is set
For the virtual exit point in continental river basin, using the virtual exit point as continental river basin river direction terminal;
Step S20, virtual exit point position grid height value in corresponding basin DEM is emptied using GIS software, is obtained
To an amendment DEM figure layer DEM1;
Step S30, according to the DEM elevation in continental river basin and actual measurement network of waterways information, the virtual network of waterways is drawn, and make the river
Net can finally be connected to some virtual exit point, and change into grid according to DEM grid positions and size, obtain with reference to network of waterways grid
Lattice;
Step S40, the continental river Basin Boundary of collection is converted to grid according to DEM grid positions and size;
Step S50, it will be burnt in DEM1 with reference to network of waterways grid and continental river Basin Boundary grid information, and respectively to ginseng
It examines network of waterways grid and continental river Basin Boundary grid position carries out elevation amendment, obtain continental river Basin Boundary grid positions
It is higher by the certain elevation of original DEM, the amendment DEM figure layer DEM2 of the certain elevation of original DEM is lower than with reference to network of waterways grid positions elevation;
Step S60, the simulation network of waterways is extracted using the hydrological analysis module of GIS software to DEM2;
Step S70, unique and unduplicated lake is carried out to each lake range in the continental river basin of collection to encode,
And grid is converted it into according to DEM grid positions and size;
Step S80, successively since each virtual exit point, simulation network of waterways grid trace to the source time according to grid flow direction
It goes through, unique network of waterways that do not repeat is carried out to each section and is encoded, and makes the network of waterways with the simulation network of waterways grid of identical lake coding
It encodes identical;
Step S90, using each simulation network of waterways grid as starting point, non-analog network of waterways grid is carried out according to grid flow direction
It traces to the source traversal, searches all non-network of waterways grids for importing present day analog network of waterways grids, and set the non-network of waterways grid found
Sub-basin is encoded to the network of waterways coding of the simulation network of waterways grid, encodes internally stream river valley progress sub-basin according still further to sub-basin and draws
Point.
In the above-mentioned methods, continental river basin to be divided area information is first collected before setting virtual exit point;
Continental river basin to be divided area information includes the continental river basin DEM elevation;Lake number in continental river basin
Amount, range and its location information;Continental river Basin Boundary information;
Wherein, continental river Basin Boundary information includes inner boundary information and outer boundary information;Outer boundary is comprising entire interior
The boundary of river valley range is flowed, inner boundary is the boundary of water resources regionalization in continental river basin.
In the above-mentioned methods, in step slo, in entire continental river basin perimeter, when presence is independent, it is more not contact
When a water resources regionalization, virtual outlet is respectively set in each independence, the water resources regionalization not contacted.
In the above-mentioned methods, in step s 30, the virtual network of waterways of drafting includes two kinds of forms:
The first is to extend to have the virtual network of waterways that the actual measurement network of waterways is virtually exported to some;
Second for large area without actual measurement network of waterways region, is newly drawn according to DEM elevation and with survey the network of waterways it is relatively independent,
It is set finally to be connected to the virtual network of waterways that some is virtually exported.
In the above-mentioned methods, in step s 30, the virtual network of waterways and the actual measurement network of waterways are first merged into one by the virtual network of waterways of drafting
A linear vector figure layer, then the linear vector figure layer after merging is changed into grid according to DEM grid positions and size.
In the above-mentioned methods, when drawing the virtual network of waterways, when the river lake Shang You, internal drawing one connection in lake
The virtual network of waterways of upstream and downstream section.
In the above-mentioned methods, in step s 50, it according to reference network of waterways grid and continental river Basin Boundary grid information, presses
Elevation amendment is carried out to reference network of waterways grid and continental river Basin Boundary grid position respectively according to priority, specifically:
The DEM1 elevation numerical value of continental river Basin Boundary grid position is increased into a fixed value;Correspondingly, will ginseng
The DEM1 elevation numerical value for examining network of waterways grid position reduces a fixed value;
If a grid is both outer boundary grid and refers to network of waterways grid, handled as boundary raster, by place
The DEM1 elevation numerical value of position increases a fixed value;
If a grid is both inner boundary grid and refers to network of waterways grid, handled as network of waterways grid, by place
The DEM1 elevation numerical value of position reduces a fixed value.
In the above-mentioned methods, lake coding, the network of waterways are encoded to the coding of the consecution natural number since 1.
In the above-mentioned methods, step S80 specifically:
Traversal of tracing to the source is carried out to simulation network of waterways grid according to grid flow direction, assigns a unique volume to each network of waterways grid
Code;
If not encountering the simulation network of waterways grid or lake grid of fork, next network of waterways Lattice encoding is equal to current
The network of waterways of grid encodes;
If encountering the simulation network of waterways grid or lake grid of fork, the network of waterways coding of next network of waterways grid is equal to
The network of waterways coding of current grid plus 1;
If encountering lake grid, and also lake having the same encodes present day analog network of waterways grid, then whether divides
The network of waterways coding of trouble, next simulation network of waterways grid is equal to present day analog network of waterways raster river network coding;
Network of waterways coding identical simulation network of waterways grid belongs to the same section;
Traversal of tracing to the source is carried out after the completion of some virtual outlet upstream river traversal, then to other virtual outlets, this is stylish
The starting network of waterways for traversal of tracing to the source is encoded to previous stage maximum coding plus 1.
Compared with prior art, the present invention is based on DEM, in conjunction with letters such as the actual measurement network of waterways, Basin Boundary, lake ranges
Breath, different sub-basin ranges is extracted using multiple virtual exit points, and on the basis of considering lake range respectively, carries out river
Net coding can guarantee that dividing obtained sub-basin substantially conforms to upstream-downstream relationship in continental river basin, and make a lake only
In a sub-basin, guarantee that sub-basin range divides correctness.
Detailed description of the invention
Fig. 1 is a kind of flow chart of the watershed partitioning method in continental river basin provided by the invention;
Fig. 2 is that there are two the specific embodiment schematic diagrames in lake and two subregions for tool in the present invention;
Fig. 3 is to carry out the schematic diagram after virtual outlet setting to Fig. 2 in the present invention;
Fig. 4 is that the schematic diagram after virtual network of waterways drafting is carried out to Fig. 3 in the present invention;
Fig. 5 is that the schematic diagram after simulation extraction of drainage is carried out to Fig. 4 in the present invention;
Fig. 6 is to carry out the schematic diagram after the coding of the network of waterways in the present invention to Fig. 5;
Fig. 7 is to realize showing for watershed partitioning after carrying out sub-basin coding to the non-analog network of waterways grid of Fig. 6 in the present invention
It is intended to.
Specific embodiment
The present invention provides a kind of watershed partitioning methods in continental river basin, based on DEM, in conjunction with actual measurement the network of waterways,
The information such as Basin Boundary, lake range extract different sub-basin ranges, Neng Goubao using multiple virtual exit points respectively
Card divides obtained sub-basin and substantially conforms to upstream-downstream relationship in continental river basin, and a lake is made to be only positioned at a sub-basin
It is interior, guarantee that sub-basin range divides correctness.The present invention is made in detail with specific embodiment with reference to the accompanying drawings of the specification
Explanation.
As shown in Figure 1, a kind of watershed partitioning method in continental river basin provided by the invention, comprising the following steps:
Step S10, setting low-lying land in continental river basin to be divided, (the local elevation in entire continental river basin perimeter is most
Low spot) or No way out lake center position be the virtual exit point in continental river basin, using the virtual exit point as continental river
Basin river direction terminal;
Continental river basin to be divided area information, including continental river basin DEM are first collected before setting virtual exit point
Elevation, lake quantity, range and its location information in continental river basin, the information such as continental river Basin Boundary;Wherein, continental river
Basin Boundary information includes inner boundary information and outer boundary information, and outer boundary is the boundary comprising entire continental river basin perimeter
(peripheral solid line in such as Fig. 3), inner boundary is the boundary (inner dotted line in such as Fig. 3) of water resources regionalization in continental river basin;Then,
According to DEM elevation in each water resources regionalization, it is arranged in the minimum point of elevation (the lake center position of low-lying land or No way out)
Virtual exit point;Due to there is the independent water resources regionalization not contacted in a continental river basin, so needing respectively above-mentioned
The virtual exit point of setting in the different water resources regionalizations of compartment, therefore according to drainage characteristics concrete condition, a continental river
Multiple virtual exit points, at least 1 virtual exit point can be set in basin.As shown in Fig. 2, in this embodiment, collection
The continental river basin area information includes 2 lakes (black region), the lake a being located on perennial stream and river lower reaches
Lake b, there are two water resources regionalization in basin, water resources regionalization boundary is dotted line in Fig. 2, and entire continental river basin is divided into point
Two water resources regionalizations of area I and subregion II;According to DEM elevation, virtual exit point is set up in subregion I and subregion II respectively,
As shown in figure 3, being respectively the circle A and circle B at partition boundaries.
Step S20, virtual exit point position grid height value in corresponding basin DEM is emptied (i.e. using GIS software
It is set as Null value), obtain the amendment DEM figure layer emptied in virtual outlet point position, referred to as DEM1.
Step S30, it according to the DEM elevation in continental river basin and the actual measurement network of waterways (linear vector) information, draws (manual)
The virtual network of waterways (linear vector), and so that the network of waterways is finally connected to some virtual exit point, and according to DEM grid positions and greatly
It is small to change into grid, it obtains with reference to network of waterways grid.
In the present invention, the virtual network of waterways of drafting includes two kinds of forms, a kind of extension to survey the network of waterways extends actual measurement
The network of waterways is to some virtual exit point;Another kind is the independent virtual network of waterways of individualism, does not have the region in river, root in large area
The network of waterways relatively independent with the actual measurement network of waterways is newly drawn according to DEM elevation, it is made finally to be connected to some virtual exit point.It has drawn virtual
Behind the network of waterways, after needing its being merged into a linear vector figure layer with the actual measurement network of waterways, then by the linear vector figure layer root after merging
Grid is changed into according to DEM grid positions and size.When the river lake Shang You, then needs to draw a virtual network of waterways in lake and connect
Connect lower reache.If Fig. 4 is that the virtual network of waterways carried out to embodiment illustrated in fig. 3 is drawn, drawn in the lake a on perennial stream
One virtual network of waterways ensures that lower reache is connected thereon, extends to the actual measurement network of waterways of subregion I and subregion II further according to DEM elevation
Corresponding virtual exit point;And in order to there is no river region to carry out watershed partitioning the upper left corner subregion II large area, at this
It is added to an independent virtual network of waterways, and is connected to correspond to virtual exit point B, thus the region can be independently divided into one
(i.e. No. 2 sections in Fig. 6, if not adding the virtual section, No. 2 and 4 work song basins will merge into 1 son to sub-basin in 7
Basin).
Step S40, the continental river Basin Boundary (linear vector) of collection is converted to grid according to DEM grid positions and size
Lattice;
Step S50, it will be burnt with reference to network of waterways grid and continental river Basin Boundary (including inner boundary and outer boundary) grid information
It records in DEM1, and elevation amendment is carried out to reference network of waterways grid and continental river Basin Boundary grid position respectively, obtain
One is emptied in virtual outlet point position, and continental river Basin Boundary grid positions are higher by the certain elevation of original DEM, with reference to the network of waterways
Grid positions elevation is lower than the amendment DEM figure layer of the certain elevation of original DEM, referred to as DEM2.
In step S50 of the present invention, according to reference network of waterways grid and continental river Basin Boundary grid information, according to priority
The elevation of DEM1 internal reference network of waterways grid, outer boundary grid, inner boundary grid position is modified, specifically:
The DEM1 elevation numerical value of continental river Basin Boundary (inner boundary and outer boundary) grid position is increased by one admittedly
Definite value (such as 1000m);Correspondingly, the DEM1 elevation numerical value of reference network of waterways grid position is reduced a fixed value (example
Such as 1000m);
If a grid is that outer boundary grid is to handle (elevation as outer boundary grid with reference to network of waterways grid again
Numerical value increases a fixed value);If a grid is that inner boundary grid is to be used as network of waterways grid with reference to network of waterways grid again
Processing (elevation numerical value reduces a fixed value).
Step S60, the simulation network of waterways is extracted using the hydrological analysis module of GIS software to DEM2, i.e., by filling out hollow, flow direction meter
Calculate, confluence accumulative total calculate, simulation the network of waterways threshold value setting and extract and etc. obtain simulation the network of waterways, as Fig. 5 chain lines be pair
The simulation network of waterways that Fig. 4 is extracted;The simulation river density generally extracted is greater than the actual measurement network of waterways, and river length is not yet
Unanimously.
Step S70, it carries out unique to each lake range (planar vector) in the continental river basin of collection and does not repeat
Lake coding, and convert it into grid according to DEM grid positions and size;Wherein, lake is encoded to continuous since 1
Natural number coding.
Step S80, successively started with each virtual exit point, simulation network of waterways grid trace to the source time according to grid flow direction
It goes through, network of waterways coding is carried out to each section, and make have the network of waterways coding of the simulation network of waterways grid of identical lake coding identical, i.e.,
The influence for considering lake region range in an encoding process, so that the lake region of one and same coding is located in identical sub-basin.In this hair
In bright, step S80 specifically:
Traversal of tracing to the source is carried out to simulation network of waterways grid according to grid flow direction, assigns a unique volume to each network of waterways grid
Code;
If not encountering the simulation network of waterways grid or lake grid of fork, next network of waterways Lattice encoding is equal to current
The network of waterways of grid encodes;
If encountering the simulation network of waterways grid or lake grid of fork, the network of waterways coding of next network of waterways grid is equal to
The network of waterways coding of current grid plus 1;
If encountering lake grid, and also lake having the same encodes present day analog network of waterways grid, then whether divides
The network of waterways coding of trouble, next simulation network of waterways grid is equal to present day analog network of waterways raster river network coding, so that being located at same
Network of waterways grid within the scope of lake is encoded with the identical network of waterways;
Network of waterways coding identical simulation network of waterways grid belongs to the same section, and the network of waterways coding of section is the company since 1
Continuous natural number coding;
Traversal of tracing to the source is carried out after the completion of some virtual outlet upstream river traversal, then to other virtual outlets, this is stylish
The starting network of waterways for traversal of tracing to the source is encoded to previous stage maximum coding plus 1.
By taking Fig. 5 as an example, the simulation network of waterways shown in fig. 5 grid is successively traced to the source since 2 virtual exit points traversal, is first traversed
Subregion II traverses subregion I again, encodes to each section, totally 18 sections, as shown in Figure 6;Wherein, the network of waterways coding of section
Object for 1 is that contraposition has in 4 sections within the scope of lower reaches water body in lake because these sections are located in the same lake
Identical lake coding, so network of waterways coding is identical, remaining network of waterways coded object is then the independent section of no fork.
Step S90, using each simulation network of waterways grid as starting point, non-analog network of waterways grid is carried out according to grid flow direction
It traces to the source traversal, finds all non-network of waterways grids for importing present day analog network of waterways grids, and set the non-network of waterways grid found
Sub-basin be encoded to the simulation network of waterways grid the network of waterways coding, finally according to sub-basin encode internally stream river valley carry out subflow
Domain divides, i.e., identical sub-basin is encoded to same sub-basin, and different sub-basins are encoded to different sub-basins, to realize continental river
Basin watershed partitioning obtains 18 sub-basins as shown in fig. 7, dividing altogether.
The invention is not limited to above-mentioned preferred forms, and anyone should learn that is made under the inspiration of the present invention
Structure change, the technical schemes that are same or similar to the present invention are fallen within the scope of protection of the present invention.
Claims (8)
1. a kind of watershed partitioning method in continental river basin, which comprises the following steps:
Step S10, it is interior for setting the lake center position of local elevation minimum point or No way out in continental river basin to be divided
The virtual exit point of river valley is flowed, using the virtual exit point as continental river basin river direction terminal;
Step S20, virtual exit point position grid height value in corresponding continental river basin DEM is emptied using GIS software,
Obtain an amendment DEM figure layer DEM1;
Step S30, according to the DEM elevation in continental river basin and actual measurement network of waterways information, the virtual network of waterways is drawn, and enable the network of waterways
It is enough to be finally connected to some virtual exit point, and grid is changed into according to DEM grid positions and size, it obtains with reference to network of waterways grid;
Step S40, the continental river Basin Boundary of collection is converted to grid according to DEM grid positions and size;
Step S50, it will be burnt in DEM1 with reference to network of waterways grid and continental river Basin Boundary grid information, and respectively to reference river
Grid lattice and continental river Basin Boundary grid position carry out elevation amendment, obtain continental river Basin Boundary grid positions and are higher by
The original certain elevation of DEM is lower than the amendment DEM figure layer DEM2 of the certain elevation of original DEM with reference to network of waterways grid positions elevation;
Step S60, the simulation network of waterways is extracted using the hydrological analysis module of GIS software to DEM2;
Step S70, it carries out unique and unduplicated lake to each lake range in the continental river basin of collection to encode, and root
Grid is converted it into according to DEM grid positions and size;
Step S80, successively since each virtual exit point, traversal of tracing to the source is carried out to simulation network of waterways grid according to grid flow direction,
Unique network of waterways that do not repeat is carried out to each section to encode, and makes the network of waterways coding with the simulation network of waterways grid of identical lake coding
It is identical;
Step S90, it using each simulation network of waterways grid as starting point, is traced to the source according to grid flow direction non-analog network of waterways grid
Traversal searches all non-network of waterways grids for importing present day analog network of waterways grid, and sets the subflow of the non-network of waterways grid found
Domain is encoded to the network of waterways coding of the simulation network of waterways grid, encodes internally stream river valley according still further to sub-basin and carries out watershed partitioning;
Step S80 specifically:
Traversal of tracing to the source is carried out to simulation network of waterways grid according to grid flow direction, assigns a unique coding to each network of waterways grid;
If not encountering the simulation network of waterways grid or lake grid of fork, next network of waterways Lattice encoding is equal to current grid
The network of waterways coding;
If encountering the simulation network of waterways grid or lake grid of fork, the network of waterways coding of next network of waterways grid is equal to current
The network of waterways coding of grid plus 1;
If encountering lake grid, and also lake having the same encodes present day analog network of waterways grid, then whether diverges, under
The network of waterways coding of one simulation network of waterways grid is equal to present day analog network of waterways raster river network coding;
Network of waterways coding identical simulation network of waterways grid belongs to the same section;
Traversal of tracing to the source is carried out after the completion of some virtual outlet upstream river traversal, then to other virtual outlets, this is stylish to trace to the source
The starting network of waterways of traversal is encoded to previous stage maximum coding plus 1.
2. the method according to claim 1, wherein first collecting continental river to be divided before setting virtual exit point
Basin area information;
Continental river basin to be divided area information includes the continental river basin DEM elevation;Lake quantity in continental river basin, model
It encloses and its location information;Continental river Basin Boundary information;
Wherein, continental river Basin Boundary information includes inner boundary information and outer boundary information;Outer boundary is to include entire continental river
The boundary of basin perimeter, inner boundary are the boundary of water resources regionalization in continental river basin.
3. according to the method described in claim 2, it is characterized in that, in step slo, in entire continental river basin perimeter,
When there is multiple water resources regionalizations that are independent, not contacting, to be respectively set in each independence, the water resources regionalization not contacted
Virtual outlet.
4. according to the method described in claim 2, it is characterized in that, in step s 30, the virtual network of waterways of drafting includes two kinds of shapes
Formula:
The first is to extend to have the virtual network of waterways that the actual measurement network of waterways is virtually exported to some;
Second for large area without actual measurement network of waterways region, is newly drawn according to DEM elevation and with survey the network of waterways it is relatively independent, make it
Finally it is connected to some virtual network of waterways virtually exported.
5. according to the method described in claim 4, it is characterized in that, in step s 30, the virtual network of waterways of drafting is first by virtual river
Net and the actual measurement network of waterways are merged into a linear vector figure layer, then by the linear vector figure layer after merging according to DEM grid positions and
Size changes into grid.
6. according to the method described in claim 4, it is characterized in that, when drawing the virtual network of waterways, when the river lake Shang You,
The virtual network of waterways of lake internal drawing one connection upstream and downstream section.
7. according to the method described in claim 2, it is characterized in that, in step s 50, according to reference network of waterways grid and continental river
Basin Boundary grid information respectively carries out reference network of waterways grid and continental river Basin Boundary grid position according to priority
Elevation amendment, specifically:
The DEM1 elevation numerical value of continental river Basin Boundary grid position is increased into a fixed value;Correspondingly, river will be referred to
The DEM1 elevation numerical value of grid lattice position reduces a fixed value;
If a grid is both outer boundary grid and refers to network of waterways grid, handled as boundary raster, by position
DEM1 elevation numerical value increase a fixed value;
If a grid is both inner boundary grid and refers to network of waterways grid, handled as network of waterways grid, by position
DEM1 elevation numerical value reduce a fixed value.
8. according to the method described in claim 2, it is characterized in that, lake coding, the network of waterways are encoded to the continuous nature since 1
Number encoder.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811313462.3A CN109472868B (en) | 2018-11-06 | 2018-11-06 | A kind of watershed partitioning method in continental river basin |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811313462.3A CN109472868B (en) | 2018-11-06 | 2018-11-06 | A kind of watershed partitioning method in continental river basin |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109472868A CN109472868A (en) | 2019-03-15 |
CN109472868B true CN109472868B (en) | 2019-10-18 |
Family
ID=65672294
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811313462.3A Active CN109472868B (en) | 2018-11-06 | 2018-11-06 | A kind of watershed partitioning method in continental river basin |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109472868B (en) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110147423B (en) * | 2019-05-21 | 2021-06-25 | 中国科学院南京地理与湖泊研究所 | Automatic division method for basin area of internal flow lake based on lake confluence relation |
CN110717251B (en) * | 2019-09-19 | 2020-11-24 | 中国水利水电科学研究院 | Distributed hydrological model sub-basin partitioning method considering multiple elements |
CN112650827B (en) * | 2020-12-31 | 2023-11-21 | 北京师范大学 | River basin and river coding method and device based on Strahler classification |
CN113128009B (en) * | 2021-04-27 | 2021-12-07 | 中国水利水电科学研究院 | Sub-basin unit division method considering mountain area plain landform difference |
CN112906252B (en) * | 2021-04-27 | 2021-12-07 | 中国水利水电科学研究院 | Sub-basin dividing method considering lake reservoir range |
CN116126983A (en) * | 2022-12-10 | 2023-05-16 | 黄河水利委员会黄河水利科学研究院 | Method for determining multi-stage drainage basin outlet position based on DEM |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102915547A (en) * | 2012-10-12 | 2013-02-06 | 吴立新 | Catchment area dividing method based on triangle patch |
CN102915227A (en) * | 2012-09-03 | 2013-02-06 | 南京师范大学 | Parallel method for large-area drainage basin extraction |
CN105160658A (en) * | 2015-08-10 | 2015-12-16 | 南京师范大学 | Ridge line extraction method based on subbasin boundaries and flow path characteristics |
CN108090164A (en) * | 2017-12-13 | 2018-05-29 | 武汉大学 | A kind of basin large scale calculates the coding method of magnanimity grid |
CN108399309A (en) * | 2018-03-16 | 2018-08-14 | 中国水利水电科学研究院 | A kind of watershed partitioning method of large scale complex topographic area hydrological distribution model |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7249257B2 (en) * | 2001-03-05 | 2007-07-24 | Digimarc Corporation | Digitally watermarked maps and signs and related navigational tools |
US9117305B2 (en) * | 2012-06-29 | 2015-08-25 | Weyerhaeuser Nr Company | Method and apparatus for deriving stream network vectors from digital elevation map data |
-
2018
- 2018-11-06 CN CN201811313462.3A patent/CN109472868B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102915227A (en) * | 2012-09-03 | 2013-02-06 | 南京师范大学 | Parallel method for large-area drainage basin extraction |
CN102915547A (en) * | 2012-10-12 | 2013-02-06 | 吴立新 | Catchment area dividing method based on triangle patch |
CN105160658A (en) * | 2015-08-10 | 2015-12-16 | 南京师范大学 | Ridge line extraction method based on subbasin boundaries and flow path characteristics |
CN108090164A (en) * | 2017-12-13 | 2018-05-29 | 武汉大学 | A kind of basin large scale calculates the coding method of magnanimity grid |
CN108399309A (en) * | 2018-03-16 | 2018-08-14 | 中国水利水电科学研究院 | A kind of watershed partitioning method of large scale complex topographic area hydrological distribution model |
Non-Patent Citations (3)
Title |
---|
《Coding of watershed and river hierarchy to support GIS-based hydrological analyses at different scales》;HORHAN T;《Computers & Geosciences》;20091231;第35卷(第3期);第688-696页 * |
《基于DEM 河网干支拓扑关系的子流域编码规则》;刘佳嘉 等;《河海大学学报(自然科学版)》;20130731;第41卷(第4期);第288-293页 * |
《基于DEM与实测河网的流域编码方法》;罗翔宇 等;《水科学进展》;20060331;第17卷(第2期);第259-264页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109472868A (en) | 2019-03-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109472868B (en) | A kind of watershed partitioning method in continental river basin | |
CN110717251B (en) | Distributed hydrological model sub-basin partitioning method considering multiple elements | |
CN103092572B (en) | The parallel method of distributing numerical control under a kind of cluster environment | |
De Jager et al. | Development and demonstration of a structured hydrological feature coding system for Europe | |
CN101864750A (en) | Multi-model meta-synthesis flood forecasting system and forecasting method thereof | |
CN112231993B (en) | Channel cross section generalization method and device based on Classification-Puck | |
CN109285219B (en) | A kind of grid type hydrological model grid calculation order encoding method based on DEM | |
CN108510008B (en) | Road network extraction method based on floating car track point spatial relationship and distribution | |
CN108875815A (en) | Feature Engineering variable determines method and device | |
CN104820724A (en) | Method for obtaining prediction model of knowledge points of text-type education resources and model application method | |
Qin et al. | OpenGF: An ultra-large-scale ground filtering dataset built upon open ALS point clouds around the world | |
CN109033322A (en) | A kind of test method and device of multidimensional data | |
CN109902133B (en) | Multi-source data error correction processing method and system based on arbitrary partition area of power grid | |
CN114661744A (en) | Terrain database updating method and system based on deep learning | |
CN113963003A (en) | Regional surface water body rapid extraction method based on Google Earth Engine | |
CN104331870A (en) | Structured multi-scale reconstruction method suitable for regular grid DEM | |
CN110737931B (en) | ArcGIS-based extraction method for key parameters of railway bridge crossing hydrologic | |
Schreiber et al. | Regionalization of low flows in southwest Germany | |
CN113159044A (en) | Deep learning-based road material identification method for convolutional neural network | |
CN108345871A (en) | A kind of cervical carcinoma slice recognition methods | |
CN115953687A (en) | Small and micro water body damage grade division method and device based on remote sensing technology | |
CN106777724B (en) | Basic flow segmentation method for semi-humid and semi-arid regions and humid regions | |
CN112508315A (en) | Flood prevention consultation system | |
CN110532969A (en) | Ramp unit division methods based on multi-scale image segmentation | |
CN115688247A (en) | Loess plateau channel area extraction and calculation method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |