CN102915547A - Catchment area dividing method based on triangle patch - Google Patents

Catchment area dividing method based on triangle patch Download PDF

Info

Publication number
CN102915547A
CN102915547A CN2012103848051A CN201210384805A CN102915547A CN 102915547 A CN102915547 A CN 102915547A CN 2012103848051 A CN2012103848051 A CN 2012103848051A CN 201210384805 A CN201210384805 A CN 201210384805A CN 102915547 A CN102915547 A CN 102915547A
Authority
CN
China
Prior art keywords
point
water
catchment area
triangle
minimum
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.)
Granted
Application number
CN2012103848051A
Other languages
Chinese (zh)
Other versions
CN102915547B (en
Inventor
吴立新
李志锋
张振鑫
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Individual
Original Assignee
Individual
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Priority to CN201210384805.1A priority Critical patent/CN102915547B/en
Publication of CN102915547A publication Critical patent/CN102915547A/en
Application granted granted Critical
Publication of CN102915547B publication Critical patent/CN102915547B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Sewage (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a catchment area dividing method based on a triangle patch, which comprises: 1) a confluence path extracting method based on a plane-point-side-point confluence mode that water of the triangle patch converges at the lowest point and then converges at the maximum side along slope; 2) a dividing method for marking a first-level catchment area of an upstream triangle based on the plane-point-side-point confluence mode; and 3) a multi-level catchment area dividing method for filling low-lying vertex in the catchment area based on the elevation of a catchment area water outlet. The catchment area dividing method based on the triangle patch can be applied to basin characteristic analysis of confluence path extraction, catchment area (multi-level) division and water outlet judgment of digital earth surface (comprising complex earth surface based on CD-triangulated irregular network (CD-TIN)) represented by the TIN, provides a new method for TIN-based digital elevation model (DEM) basin characteristic analysis, and can develop the corresponding software for the digital hydrological field (applied to urban waterlogging flood process simulation and simulated analysis).

Description

A kind of water catchment area division methods based on triangle surface
Technical field
The present invention relates to based on the digital elevation model water catchment area division methods of (Digital Elevation Model, DEM divide TIN (Triangulated Irregular Network, TIN) and two kinds of expression waies of grid).More specifically, relate to a kind of water catchment area division methods based on triangle surface.
Background technology
The water catchment area, or claim catchment area, retaining basin, drainage basin, referring to that rainwash or other materials converge to the ground table section of flowing through in the process of a common water delivering orifice, it is the zone of a sealing.Water delivering orifice refers to that current leave the point of watershed, is the borderline minimum point of watershed.Exist the landform more low-lying than peripheral ground to be called the depression in the water catchment area.
The water catchment area is the basis of hydrological analysis, also is the elementary cell of the area researches such as basin and flood.Current, based on the water catchment area division methods of DEM two large types are arranged: based on grid with based on TIN.Wherein, because raster data structure is simple, algorithm process is convenient, and is more based on the water catchment area partitioning algorithm of grid.Determine the depression method such as propositions such as Band (1986) and Qian (1990) with rectangular window scanning DEM matrix; O ' Callaghan and Mark (1984) propose based on rainwash cross flow model, simulated flow and then definite watershed line are divided watershed.And wherein the water (flow) direction algorithm has: based on D8 method, DEMON method and the D ∞ method of unimodal current, based on the Quinn method of polymodal current and Freeman method etc.But, based on the digital earth's surface of triangle surface approaching to reality landform more, especially better to complicated earth surface TIN expression effect.Therefore the hydrological analysis of being correlated with based on TIN can access more accurately result.Liu Xuejun has proposed a kind of river network extraction algorithm based on TIN in " based on the extraction of drainage algorithm of TIN " on " Journal of Hydraulic Engineering " the 1st phase of the 39th volume in 2008, and Zhou Qiming has simulated surface water flow path at " Water Resource Research " " Estimating surface flowpaths on a digital elevation model using a triangular facet network " on the 47th phase based on leg-of-mutton DEM in 2011.But below all do not provide the water catchment area generating algorithm based on triangle surface.
To sum up, the method that is based on grid DEM is divided in the water catchment area more at present, and divides comparatively rare based on the water catchment area of triangle surface.Patent of the present invention is taken the current rule into account, has proposed the face-point-Bian-point sink stream mode on a kind of triangle surface, has provided based on this pattern that flow concentration path extracting method, water catchment area (multistage) are divided and the technical method of water delivering orifice judgement.
Summary of the invention
(1) technical matters that solves
Current hydrological analysis is many based on grid DEM, and accuracy and dirigibility are not good enough, and based on the DEM of TIN approaching to reality earth's surface preferably, but divide for its water catchment area (multistage), do not provide corresponding solution.Patent of the present invention has provided the technical method that flow concentration path extraction, water catchment area (multistage) division and water catchment area water delivering orifice are judged based on triangle surface, and the method can be applicable to the digital earth's surface (comprising the complicated earth surface that Constraint-based sealene triangle net (CD-TIN) is expressed) that TIN represents.
(2) technical scheme
For ease of hereinafter setting forth, now relational language is explained as follows:
The upstream triangle: refer to according to the pattern of confluxing of technical scheme hereinafter, all triangles that are flowed to certain summit along flow concentration path are called the upstream triangle on this summit; Wherein, the upstream triangle on each minimum summit of triangle elevation initially is denoted as the triangle at its place.
One-level water catchment area: refer to based on the TIN earth's surface to divide first the sealing ground table section that forms according to the pattern of confluxing of technical scheme hereinafter.
Water catchment area water delivering orifice: refer to the minimum point of height value in the marginal point of water catchment area.
Low-lying summit: refer to that height value in the water catchment area is lower than the summit of water delivering orifice; Wherein, the low-lying summit that local elevation is minimum in the water catchment area is called local depression point.
Parent water catchment area: refer to the division with respect to the upper level water catchment area, fill and lead up low-lying summit in the upper level water catchment area according to the water delivering orifice height value of upper level water catchment area, then according to the pattern of confluxing of technical scheme hereinafter, again divide the sealing ground table section that forms.As, the parent water catchment area after the one-level water catchment area is upwards divided is the secondary water catchment area, and the sub-water catchment area of secondary water catchment area is the one-level water catchment area, and the rest may be inferred.
Multistage water catchment area: refer to for the one-level water catchment area, according to the division result of one-level water catchment area, according to the secondary water catchment area of the hereinafter multistage water catchment area division methods formation of technical scheme, three grades of water catchment areas etc., be referred to as multistage water catchment area.
1, the flow concentration path extraction is the prerequisite that the water catchment area is divided, and is flowed to each minimum summit of triangle elevation according to the water on the triangle surface, is flowed to Close edge gradient steepest summit along triangle edges, until the terminal point that confluxes; The water yield of confluxing on each summit of triangle that indicates according to this mode of confluxing is followed the trail of the extraction flow concentration path from low to high.The water based on triangle surface that the technical program relates to is flowed to the flow concentration path extracting method of face-point of confluxing along gradient steepest limit after the minimum point-Bian-point sink stream mode, mainly comprises following step:
(1) all triangular apex in the initialization area, the unified initial water yield that indicates each summit is a;
(2) each triangle in the traversal zone is pooled to the minimum summit of this triangle elevation with each leg-of-mutton water, and namely the sign water yield of each minimum point is 3a, and the sign water yield on all the other two summits does not become a; Store each leg-of-mutton minimum point in chained list L, and sort from high to low by its elevation;
(3) get certain leg-of-mutton minimum point P among the chained list L, the limit abutment points of search TIN net mid point P.Be lower than the limit abutment points such as fruit dot P height value, then the water yield does not flow out; If the height value of P is higher than its limit abutment points, then get the limit abutment points S of gradient steepest, with the water yield inflow S of minimum point P, the water yield of namely putting S increases the water yield of P, and the water yield of some P indicates constant.Gauge point P and the searched mistake of some S; Get a S and repeat this step.
(4) to all unlabelled some repeating steps (3) among the chained list L, obtain the final each point water yield of confluxing and indicate;
(5) water yield of confluxing according to each point indicates, and water withdrawal indicates the point greater than a, by the limit neighbouring relations, indicates ascending tracking according to the water yield and connects each point formation flow concentration path.Wherein, the water yield indicates larger expression the closer to the terminal point that confluxes, and the water yield sign value on the main flow concentration path is greater than the water yield sign value on the flow concentration path, and it is local depression point that the local water yield indicates maximum point.
2, according to the flow concentration path of face-point-Bian-point sink stream mode, indicate the upstream triangle, according to the one-level water catchment area in the final sign zoning.The division methods based on face-point-Bian-point sink stream mode sign leg-of-mutton one-level water catchment area, upstream that the technical program relates to mainly comprises following step:
(1) all leg-of-mutton minimum points in the zoning are deposited among the chained list Q, and initially indicate the upstream triangle of each minimum point;
(2) if chained list Q non-NULL is got the minimum point P among the chained list Q successively, do not search for if P is labeled, then entered (3), otherwise got the next one; If chained list Q is sky then withdraws from;
(3) the limit abutment points of search point P deposits among the chained list N;
(4) compare with the height value of chained list N mid point such as the height value of fruit dot P and be not minimum, enter (5); Otherwise enter (7);
(5) value of slope of each point among calculating P and the chained list N; Get some M corresponding to gradient steepest limit and be deposited among the flow concentration path chained list R, increase the upstream triangle that some P indicates in the upstream triangle that some M indicates, indicate the searched mistake of P;
(6) will put M and be assigned to P, forward (3) to;
(7) record this point and be the terminal point that confluxes (local depression point), indicate the upstream triangle that flows into this point;
(8) sign of basis (7), the one-level water catchment area in the formation zone.
3, the water catchment area water delivering orifice is the minimum point of height value in the marginal point of water catchment area.Have number of times according to the quilt of triangle edges and judge marginal edge, get that the elevation minimum point is water delivering orifice in the summit of marginal edge, the summit that height value is lower than water delivering orifice in the water catchment area is called low-lying summit, and water delivering orifice is judged and mainly comprised following step:
(1) all limits in the search water catchment area deposit among the chained list C;
(2) get limit among the chained list C, judge whether this limit is had by certain triangle in the water catchment area, record it and had number of times;
(3) if its number of times that is had is 2, then this limit is the internal edges of water catchment area; If its number of times that is had is 1, then this limit of mark is the marginal edge of water catchment area;
(4) point of getting on the marginal edge is marginal point, and getting the point that elevation is minimum in the marginal point is water delivering orifice.
4, extract the water catchment area water delivering orifice according to summary of the invention 3, fill and lead up low-lying summit in the water catchment area according to the water delivering orifice elevation of water catchment area, then indicate the multistage water catchment area of triangle recursive generation, upstream based on face-point-Bian-point sink stream mode.Fill and lead up the multistage water catchment area division methods on low-lying summit in the water catchment area based on water catchment area water delivering orifice elevation, mainly comprise following step:
(1) water delivering orifice of extraction water catchment area;
(2) fill and lead up low-lying summit in the water catchment area according to the water delivering orifice elevation of water catchment area, make it a little more than the water delivering orifice elevation;
(3) according to the method for technical scheme 2, generate the parent water catchment area of this water catchment area;
(4) recursion step (1)~(3) can obtain the multistage water catchment area in the zone.
(3) beneficial effect
1. patent of the present invention adopts more accurately approaching to reality earth's surface of TIN data structure, proposed innovatively based on face-point-Bian-the pattern of confluxing, invented flow concentration path extracting method, water catchment area division methods based on this.
2. patent of the present invention has number of times according to the quilt of triangle edges and judges marginal edge, gets that the elevation minimum point is water delivering orifice in the summit of marginal edge; Invented according to water catchment area water delivering orifice elevation and filled and led up low-lying summit in the water catchment area, based on face-point-Bian-point sink stream mode mark upstream triangle, the method for the multistage water catchment area of recursive generation.
3. the flow concentration path of patent of the present invention extracts, water delivering orifice is judged, water catchment area (multistage) division methods, for the DEM drainage characteristics analysis based on TIN provides new method, Facing Digital hydrology field (as be applied to the city waterlogging flood process simulation and simulation analysis) exploitation corresponding software accordingly.
Description of drawings
Fig. 1 is certain regional TIN net schematic diagram, letter representation Ground Point, the figure denote height value in the bracket.
Fig. 2 is flow concentration path and sign water yield schematic diagram,
Figure BSA00000787945600041
Indicate the leg-of-mutton water yield and flow into minimum point,
Figure BSA00000787945600042
Indicate flow concentration path,
Figure BSA00000787945600043
Indicate local depression point,
Figure BSA00000787945600044
Indicate the water yield of this point.
Fig. 3 is that schematic diagram is divided in the one-level water catchment area, and distinct symbols has indicated different water catchment areas:
Figure BSA00000787945600045
Indicate water catchment area J, Indicate water catchment area F,
Figure BSA00000787945600051
Indicate water catchment area D.
Fig. 4 is the water delivering orifice schematic diagram of one-level water catchment area,
Figure BSA00000787945600052
Indicate the water delivering orifice of water catchment area.
Fig. 5 is that schematic diagram is divided in the secondary water catchment area, and distinct symbols has indicated different water catchment areas:
Figure BSA00000787945600053
Indicate water catchment area H, Indicate local depression point,
Figure BSA00000787945600055
Indicate the water delivering orifice of water catchment area.
Embodiment
In order to be illustrated more clearly in the present invention, below in conjunction with case study on implementation and accompanying drawing the present invention is described further.Below specifically described content be illustrative and nonrestrictive, should not limit protection scope of the present invention with this.
What need to further specify is: relate to the comparison of height value and value of slope among the present invention, when the situation that the multiple spot height value equates or polygon value of slope equates occurring, can sort according to certain rule (as pressing the ID size of period).Other special circumstances processing mode does not affect whole content of the present invention, should not replace, limit or improve the present invention with the processing mode of special circumstances.
Adopt the foregoing invention content can obtain the flow concentration path of target area and the division result of water catchment area (multistage).This result is the earth's surface elevation feature of target area and the true reflection of surface configuration feature, can be the digital hydrology field earth's surface hydrology characteristic information on basis is provided, can be the city waterlogging that the heavy rain charge for remittance causes such as the water catchment area, city and flood the analytic unit that provides basic.Hereinafter provide specific embodiment, the summary of the invention of patent of the present invention is described.Take certain regional TIN net shown in Figure 1 as example, the implementation step is as follows:
1. the flow concentration path that is flowed to face-point of confluxing along gradient steepest limit after the minimum point-Bian-point sink stream mode based on the water of triangle surface extracts:
(1) A is to the P point in the initialization area, and the unified sign each point water yield is a, therefore each leg-of-mutton initial water yield is denoted as 3a;
(2) each triangle in the traversal zone, each leg-of-mutton sign water yield is pooled to the minimum summit of this triangle elevation, as shown in Figure 2, point D is the minimum summit of △ BDC, △ CDK, △ KDE, △ EDB, therefore it indicates the water yield is 12a, the sign water yield that in like manner can get other minimum summit F, H, I, J, M is respectively 15a, 3a, 3a, 18a, 6a, and the sign water yield on all the other summits is constant (to be a); To put D, F, H, I, J, M are deposited among the chained list L, and just sort according to elevation and to be M, I, H, J, F, D; Annotate: the elevation of this routine mid point F, D is identical, by further specifying above, gets the more a little bigger F ordering of period ID front; Vice versa;
(3) get first chained list L intermediate cam shape minimum point M, its limit abutment points I, J, L, P, N in the search TIN net, some I, J height value all are lower than a M, but the gradient steepest of limit MJ flows into J therefore will put the sign water yield (for 6a) of M.The sign water yield of point J increases 6a (becoming 24a), and the sign water yield of some M is constant.The searched mistake of gauge point M, J;
(4) get successively chained list L intermediate cam shape minimum point I, search for its limit abutment points G, H, J, M, N, some G, H, J elevation all are lower than an I, but the gradient steepest of limit IJ, therefore the sign water yield (being 3a) of I is flowed into J.The sign water yield of J increases 3a (becoming 27a), the sign water yield constant (being 3a) of some I.The searched mistake of gauge point I;
(5) get successively chained list L intermediate cam shape minimum point H, search for its limit abutment points E, F, G, I, J, some F, J elevation all are lower than a H, but the gradient steepest of limit HF, therefore the sign water yield (being 3a) of H is flowed into some F.The sign water yield of point F increases 3a (becoming 18a), the sign water yield constant (being 3a) of some H.The searched mistake of gauge point F, H;
(6) get successively the minimum point J of chained list L intermediate cam shape, mark is skipped;
(7) get successively the minimum point F of chained list L intermediate cam shape, mark is skipped;
(8) get successively the minimum point D of chained list L intermediate cam shape, search for its limit abutment points B, C, E, K, its elevation all is higher than the D point, and D water does not flow out, the searched mistake of gauge point D;
(9) finally conflux each point the sign water yield as shown in Figure 2, water withdrawal indicates the point greater than a, by the limit neighbouring relations, indicate ascending connection each point according to the water yield and form 3 flow concentration path: H → F, I → J, M → J, 3 local depression points: some D (the sign water yield is 12a), some F (the sign water yield is 18a) and some J (the sign water yield is 27a).
2. indicate the division of leg-of-mutton one-level water catchment area, upstream based on face-point-Bian-point sink stream mode:
(1) all leg-of-mutton minimum point D, F, H, I, J, M in the zoning are deposited among the chained list Q; And initially indicate the upstream triangle of each minimum point: the upstream triangle of some D is △ BDC, △ CDK, △ KDE, △ EDB, the upstream triangle of point F is △ AFB, △ BFE, △ EFH, △ HFG, △ GFA, the upstream triangle of point H is △ GHI, the upstream triangle of point I is △ MIN, the upstream triangle of point J is △ EJK, △ KJL, △ LJM, △ MJI, △ IJH, △ HJE, and the upstream triangle of some M is △ LMP, △ PMN;
(2) get first chained list Q mid point D, be not labeled and search for, searching for its neighborhood point B, C, E, K deposits among the chained list N, the height value of this point is minimum to be local depression point, this is searched for mark, mark upstream triangle is attributed to water catchment area D (being convenient expression, with local depression point name water catchment area);
(3) then get chained list Q mid point F, be not labeled and search for, search for its neighborhood point A, B, E, H, G and deposit among the chained list N, the height value of some F is minimum to be local depression point, and this is searched for mark, and mark upstream triangle is attributed to water catchment area F;
(4) then get chained list Q mid point H, be not labeled and search for, search for its neighborhood point E, F, G, I, J, point F, J elevation all are lower than a H, but limit HF gradient steepest, store in the flowing water path, increase the upstream triangle that some H indicates in the upstream triangle that some F indicates, the searched mistake of gauge point H; Get a F and repeat step, F is local depression point, and then mark upstream triangle is attributed to water catchment area F;
(5) then get chained list Q mid point I, search neighborhood point G, H, J, M, N, some G, H, J elevation all are lower than an I, but the gradient steepest of limit IJ, therefore the searched mistake of gauge point I stores in the flowing water path.Get a J and repeat step, J is local depression point, and then mark upstream triangle is attributed to water catchment area J;
(6) then get a J, searched mistake, mark upstream triangle is attributed to water catchment area J;
(7) then get chained list Q mid point M, search for its neighborhood point I, J, L, P, N, point I, J elevation all are lower than a M, but the gradient steepest of limit MJ, the searched mistake of gauge point M, store in the flowing water path, increase the upstream triangle that some M indicates in the upstream triangle that some J indicates, the searched mistake of gauge point M; Get a J and repeat step, J is local depression point, and then mark upstream triangle is attributed to water catchment area J;
(8) indicate difference according to leg-of-mutton ownership water catchment area, triangle is belonged in the different water catchment areas, the one-level water catchment area effect of generation as shown in Figure 3.
3. have number of times according to the quilt of triangle edges and judge marginal edge, extract that the elevation minimum point is water delivering orifice in the summit of marginal edge.Water catchment area D in Fig. 3 is as example, and these water catchment area water delivering orifice computing method are as follows:
(1) all limit BC, CD, BD, EK, ED, DK, CK, the BE in the search water catchment area D deposits among the chained list C;
(2) get limit BC among the chained list C, this limit is only had by △ BCD, and being had number of times is 1;
(3) limit CK, BE, EK situation are with limit BC, and being had number of times is 1;
(4) get limit CD among the chained list C, this limit is had by △ BCD and △ CDK, and being had number of times is 2;
(5) limit CD, BD, ED, DK situation are with limit CD, and being had number of times is 2;
(6) to have number of times be 1 to the quilt of limit BC, CK, BE, EK, then is the marginal edge of water catchment area D;
Get (7) that the summit on marginal edge is marginal point in (6), get that elevation minimum point E is water delivering orifice in the marginal point.
All the other water catchment area water delivering orifice computing method are similar, generate the result as shown in Figure 4.Point E is the water delivering orifice of water catchment area D, and some H is the water delivering orifice of water catchment area F and water catchment area J.
4. fill and lead up the multistage water catchment area division on low-lying summit in the water catchment area based on water catchment area water delivering orifice elevation:
(1) generates the one-level water catchment area according to embodiment 2, such as Fig. 3;
(2) the water catchment area water delivering orifice that calculates according to embodiment 3 is such as Fig. 4;
(3) (value of this small quantity only otherwise the one-piece construction that changes this regional DEM gets final product the elevation of the water outlet stomion E of water catchment area D to be added small quantity 0.01, should not limit the invention accordingly equally, hereinafter similar) be assigned to local depression point D, the elevation of the water outlet stomion H of water catchment area F adds small quantity 0.01 and is assigned to local depression point F, and the elevation of the water outlet stomion H of water catchment area J adds small quantity 0.01 and is assigned to local depression point J;
(4) copy the step of embodiment 2, generate the parent water catchment area, i.e. secondary water catchment area, as shown in Figure 5;
Above-described specific embodiment; purpose of the present invention, technical scheme and beneficial effect are further described; institute is understood that; the above only is specific embodiments of the invention; be not limited to the present invention; within the spirit and principles in the present invention all, any modification of making, be equal to replacement, improvement etc., all should be included within protection scope of the present invention.

Claims (5)

1. the water catchment area division methods based on triangle surface is characterized in that, the method comprises:
(1) is flowed to the flow concentration path extracting method of face-point of confluxing along gradient steepest limit after the minimum point-Bian-point sink stream mode based on the water of triangle surface;
(2) indicate the division methods of leg-of-mutton one-level water catchment area, upstream based on face-point-Bian-point sink stream mode;
(3) fill and lead up the multistage water catchment area division methods on low-lying summit in the water catchment area based on water catchment area water delivering orifice elevation.
2. the water based on triangle surface described in according to claim 1 is flowed to the flow concentration path extracting method of face-point of confluxing along gradient steepest limit after the minimum point-Bian-point sink stream mode, it is characterized in that, water on the triangle surface is flowed to each minimum summit of triangle elevation, be flowed to Close edge gradient steepest summit along triangle edges, until the terminal point that confluxes; The water yield of confluxing on each summit of triangle that indicates according to this mode of confluxing is followed the trail of the extraction flow concentration path from low to high, mainly comprises following step:
(1) all triangular apex in the initialization district, the unified initial water yield that indicates each summit is a;
(2) each triangle in the traversal district is pooled to the minimum summit of this triangle elevation with each leg-of-mutton water, and namely the sign water yield on each minimum summit is 3a, and the sign water yield on all the other two summits does not become a; Store each leg-of-mutton minimum point in chained list L, and sort from high to low by its elevation;
(3) get certain leg-of-mutton minimum point P among the chained list L, the limit abutment points of search TIN net mid point P.Be lower than its limit abutment points such as fruit dot P height value, then the water yield does not flow out; If the height value of P is higher than its limit abutment points, then get the limit abutment points S of gradient steepest, with the water yield inflow S of minimum point P, the water yield of namely putting S increases the water yield of P, and the water yield of some P indicates constant.Gauge point P and the searched mistake of some S; Get a S and repeat this step;
(4) to all unlabelled some repeating steps (3) among the chained list L, obtain the final each point water yield of confluxing and indicate;
(5) water yield of confluxing according to each point indicates, and water withdrawal indicates the point greater than a, by the limit neighbouring relations, indicates ascending tracking according to the water yield and connects each point formation flow concentration path.Wherein, the larger expression of the water yield is the closer to the terminal point that confluxes, and the water yield sign value on the main flow concentration path is greater than the water yield sign value on the flow concentration path, and it is local depression point that the local water yield indicates maximum point.
3. indicate the division methods of leg-of-mutton one-level water catchment area, upstream based on face-point-Bian-point sink stream mode described in according to claim 1, it is characterized in that, according to the flow concentration path of face-point-Bian-point sink stream mode, indicate the upstream triangle, mainly comprise following step:
(1) all leg-of-mutton minimum points in the zoning are deposited among the chained list Q, and initially indicate the upstream triangle of each minimum point;
(2) if chained list Q non-NULL is got the minimum point P among the chained list Q successively, do not search for if P is labeled, then entered (3), otherwise got the next one; If chained list Q is sky then withdraws from;
(3) the limit abutment points of search point P deposits among the chained list N;
(4) compare with the height value of chained list N mid point such as the height value of fruit dot P and be not minimum, enter (5); Otherwise enter (7);
(5) value of slope of each point among calculating P and the chained list N; Get some M corresponding to gradient steepest limit and be deposited among the flow concentration path chained list R, increase the upstream triangle that some P indicates in the upstream triangle that some M indicates, indicate the searched mistake of some P;
(6) will put M and be assigned to P, forward (3) to;
(7) record this point and be the terminal point that confluxes, indicate the upstream triangle that flows into this point;
(8) sign of basis (7), the one-level water catchment area in the formation zone.
4. fill and lead up the multistage water catchment area division methods on low-lying summit in the water catchment area based on water catchment area water delivering orifice elevation described in according to claim 1, it is characterized in that, extract the water catchment area water delivering orifice, fill and lead up low-lying summit in the water catchment area according to the water delivering orifice elevation of water catchment area, make it a little more than the water delivering orifice elevation, indicate the multistage water catchment area of leg-of-mutton method recursive generation, upstream based on face-point-Bian-point sink stream mode.
5. the extraction water catchment area water delivering orifice described in according to claim 4 is characterized in that, having number of times according to the quilt of triangle edges in the water catchment area is 1 to be judged to be marginal edge, gets that the elevation minimum point is water delivering orifice in the summit of marginal edge.
CN201210384805.1A 2012-10-12 2012-10-12 A kind of water catchment area division methods based on triangle surface Expired - Fee Related CN102915547B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210384805.1A CN102915547B (en) 2012-10-12 2012-10-12 A kind of water catchment area division methods based on triangle surface

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210384805.1A CN102915547B (en) 2012-10-12 2012-10-12 A kind of water catchment area division methods based on triangle surface

Publications (2)

Publication Number Publication Date
CN102915547A true CN102915547A (en) 2013-02-06
CN102915547B CN102915547B (en) 2016-03-02

Family

ID=47613899

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210384805.1A Expired - Fee Related CN102915547B (en) 2012-10-12 2012-10-12 A kind of water catchment area division methods based on triangle surface

Country Status (1)

Country Link
CN (1) CN102915547B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103345589A (en) * 2013-07-19 2013-10-09 吴立新 Urban area catchment unit partition method with consideration of constraint characteristics
CN105279317A (en) * 2015-09-30 2016-01-27 西北农林科技大学 DEM based flat ground river network water flow direction estimation method
CN105740812A (en) * 2016-01-29 2016-07-06 武汉大学 Digital surface model based city catchment area extraction method
CN109472868A (en) * 2018-11-06 2019-03-15 中国水利水电科学研究院 A kind of watershed partitioning method in continental river basin
CN113128009A (en) * 2021-04-27 2021-07-16 中国水利水电科学研究院 Sub-basin unit division method considering mountain area plain landform difference

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040215428A1 (en) * 2003-04-28 2004-10-28 Massachusetts Institute Of Technology Method for producing property-preserving variable resolution models of surfaces
CN102722909A (en) * 2012-05-28 2012-10-10 武汉大学 Assembly line topology network dynamic simulation method based on adaptive-dimensional DEM (dynamic effect model)

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040215428A1 (en) * 2003-04-28 2004-10-28 Massachusetts Institute Of Technology Method for producing property-preserving variable resolution models of surfaces
CN102722909A (en) * 2012-05-28 2012-10-10 武汉大学 Assembly line topology network dynamic simulation method based on adaptive-dimensional DEM (dynamic effect model)

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
刘学军等: "基于不规则三角网的河网提取算法", 《水利学报》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103345589A (en) * 2013-07-19 2013-10-09 吴立新 Urban area catchment unit partition method with consideration of constraint characteristics
CN105279317A (en) * 2015-09-30 2016-01-27 西北农林科技大学 DEM based flat ground river network water flow direction estimation method
CN105279317B (en) * 2015-09-30 2018-04-03 西北农林科技大学 A kind of level land Flow of River direction estimating method based on DEM
CN105740812A (en) * 2016-01-29 2016-07-06 武汉大学 Digital surface model based city catchment area extraction method
CN105740812B (en) * 2016-01-29 2018-10-26 武汉大学 A kind of city water catchment area extracting method based on digital surface model
CN109472868A (en) * 2018-11-06 2019-03-15 中国水利水电科学研究院 A kind of watershed partitioning method in continental river basin
CN109472868B (en) * 2018-11-06 2019-10-18 中国水利水电科学研究院 A kind of watershed partitioning method in continental river basin
CN113128009A (en) * 2021-04-27 2021-07-16 中国水利水电科学研究院 Sub-basin unit division method considering mountain area plain landform difference

Also Published As

Publication number Publication date
CN102915547B (en) 2016-03-02

Similar Documents

Publication Publication Date Title
CN103236086B (en) One takes the contextual multiple dimensioned DEM modeling method of the earth's surface hydrology into account
Gholami et al. Modeling of groundwater level fluctuations using dendrochronology in alluvial aquifers
Hadded et al. A decision support system to manage the groundwater of the zeuss koutine aquifer using the WEAP-MODFLOW framework
Zhang et al. A remote sensing method for estimating regional reservoir area and evaporative loss
Chen et al. Exploring the spatio-temporal interrelation between groundwater and surface water by using the self-organizing maps
CN102915547B (en) A kind of water catchment area division methods based on triangle surface
Nalbantis et al. Holistic versus monomeric strategies for hydrological modelling of human-modified hydrosystems
Kenny et al. Routing overland flow through sinks and flats in interpolated raster terrain surfaces
Van Zijl et al. Functional digital soil mapping: a case study from Namarroi, Mozambique
CN104200387A (en) Land use zoning method giving consideration to geographic entity semantic similarity
Liang et al. Connecting urban area with rural hinterland: A stepwise ecological security network construction approach in the urban–rural fringe
CN112199456B (en) Basin water resource management method based on supply type hydrological ecosystem service
Villarreal et al. The influence of wastewater subsidy, flood disturbance and neighbouring land use on current and historical patterns of riparian vegetation in a semi‐arid watershed
Zhu et al. Loess terrain segmentation from digital elevation models based on the region growth method
Haag et al. Development of a data model to facilitate rapid watershed delineation
Wang et al. Multi-objective optimization of the coastal groundwater abstraction for striking the balance among conflicts of resource-environment-economy in Longkou City, China
García et al. Deriving spatially explicit water uses from land use change modelling results in four river basins across Europe
Zhang et al. An enhanced approach for surface flow routing over drainage‐constrained triangulated irregular networks
Wu et al. Effect of wetland hydrological connectivity on runoff based on a basin comparison
Zhang et al. Construction of multi-level ecological security network in fragmented nature landscape using the three-dimensional framework of ecological adaptability
CN114004141A (en) River basin groundwater super-mining area prediction decision platform
CN114648617A (en) Water system extraction method based on digital elevation model DEM
Sharma et al. Long-term sustainability of groundwater resources in the central Ganga Alluvial Plain, India: Study from Gomti River Basin
CN105808947A (en) Intuitive analysis method of ecological water demand monthly scale on the basis of water demand for sediment transport
O'Neil et al. Quantifying local flow direction uncertainty

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160302

Termination date: 20161012

CF01 Termination of patent right due to non-payment of annual fee