CN110070619A - Take the river water surface three-dimensional visualization method of water surface gradient into account - Google Patents
Take the river water surface three-dimensional visualization method of water surface gradient into account Download PDFInfo
- Publication number
- CN110070619A CN110070619A CN201910278875.0A CN201910278875A CN110070619A CN 110070619 A CN110070619 A CN 110070619A CN 201910278875 A CN201910278875 A CN 201910278875A CN 110070619 A CN110070619 A CN 110070619A
- Authority
- CN
- China
- Prior art keywords
- water surface
- water level
- boundary point
- river
- website
- 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
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/20—Finite element generation, e.g. wire-frame surface description, tesselation
Abstract
The invention discloses a kind of river water surface three-dimensional visualization methods for taking water surface gradient into account.It includes the following steps, step 1: data preparation；Step 2: boundary point density is carried out to water surface boundary polygon；Step 3: polygon trigonometric ratio；Step 4: interpolation and adjustment are carried out based on elevation of the water level process data to boundary point；Step 5: water surface drafting is carried out based on the triangulation network adjusted.The present invention has the advantages of guaranteeing that water surface effect of visualization is truer, flexible, practical, strong applicability.
Description
Technical field
The present invention relates to three-dimensional artificials and visualization technique field, more specifically say it is the urban river water for taking water surface gradient into account
Face three-dimensional visualization method.
Background technique
Existing to be mainly divided to two classes about water surface three-dimensional visualization method: the first kind is texture water, mainly passes through dynamic line
The technologies such as reason, Z-Correct bump mapping Z-correct, projective textures, the visual effects such as reflection, refraction, the glistening light of waves of mimetic surface, this method efficiency compared with
Height, but the heaving of the sea of water body can not be embodied.Such as paper " realization of water body live effect in large-scale virtual scene ", " it is based on
The real-time water surface environmental simulation of Cg and OpenGL " etc..Second class is vertex water, mainly passes through wave spectrum simulation, Fourier transformation, letter
The mathematical methods such as array conjunction, fluid dynamics simulate the changing rule of water body grid vertex elevation, and then embody the wave of the water surface
Unrestrained fluctuation effect, this method fluctuation effect is true, but calculation amount is larger, inefficient, as paper " is based in 3d gaming
The dynamic water surface simulation algorithm of OGRE ", paper " a kind of sense of reality ripples face modeling method of Virtual Environment-Oriented ", patent " one
Real-time dynamic water surface analogy method of the kind based on GPU " etc..
These two kinds of methods all do not account for that natural river course upstream water level is high, the level of tail water is low, that is, there is the spy of water surface gradient
Property, thus it is difficult to the water level process three-dimensional visualization suitable for river basin level.
The fields such as the main Water Conservancy Information of water surface three-dimensional visualization, 3D game, digital city carry out three-dimensional visible to the water surface
Change the sense of reality that the water surface can be enhanced.At a wide range of (such as river basin level) if when river water surface three-dimensional visualization does not consider water surface ratio
Drop, will affect authenticity and visual effect that the river two sides water surface floods.
Therefore, now it is badly in need of researching and developing a kind of side of heaving of the sea water surface three-dimensional visualization that the effect is more real for embodying water body
Method.
Summary of the invention
The purpose of the invention is to provide a kind of river water surface three-dimensional visualization method for taking water surface gradient into account, guarantee water
Face the effect of visualization more true heaving of the sea of water body (can embody the effect is more real), flexible, practical, strong applicability.
To achieve the goals above, the technical solution of the present invention is as follows: taking the river water surface three-dimensional visualization of water surface gradient into account
Method, it is characterised in that: include the following steps,
Step 1: data preparation；
Step 2: boundary point density is carried out to water surface boundary polygon；
Step 3: polygon trigonometric ratio；
Step 4: interpolation and adjustment are carried out based on elevation of the water level process data to boundary point；
Step 5: water surface drafting is carried out based on the triangulation network adjusted.
In the above-mentioned technical solutions, in step 1, the data of preparation are river available data, the river available data packet
Include river water surface boundary polygon data, hydrographic station stage hydrograph data along river；Hydrographic station along river
Stage hydrograph data include that hydrographic station (the following abbreviation survey stations of hydrographic station) are numbered, x coordinate, y-coordinate, water level, when
Between；Wherein, survey station number is successively increased since 0 according to the sequence for being from upstream to downstream to n, and n is integer and is greater than or equal to
0。
In the above-mentioned technical solutions, in step 2, by the way that density point spacing is arranged, the boundary points after density are controlled
Amount.
In the above-mentioned technical solutions, in step 3, trigonometric ratio is carried out to the boundary polygon after density, it will be more after density
Side shape is converted into a series of triangle net collectives.
In the above-mentioned technical solutions, in step 4, the algorithm of boundary point water level elevation interpolation and adjustment includes the following steps:
S41: the Polygonal Boundary point P after traversal density；
S42: it calculates and current polygon boundary point P nearest survey station Si, number is i, if the water level at a certain moment is
li；
S43: the interpolation water level l of Polygonal Boundary point P is calculated by neighboring stations water levelp；Interpolation water level lpIn following
Hold:
1) as i=0, website SiTo originate survey station, no upstream survey station, downstream survey station Si+1；Website P is calculated separately to directly
Line SiSi+1Intersection point F0With website Si、Si+1Distance di、di+1；If Si+1Water level is li+1, the linear interpolation water level l of boundary point Pp
Meet (lp–li)/(lp–li+1)=di/di+1, then lp=(lidi+1-li+1di)/(di+1-di)；
2) as i=n, then website SiTo terminate survey station, no downstream survey station, upstream survey station Si-1；Website P is calculated separately to arrive
Straight line SiSi-1Intersection point F0With website Si、Si-1Distance di、di-1；If Si-1Water level is li-1, the linear interpolation water level of boundary point P
lpMeet (lp–li)/(lp–li-1)=di/di-1, then lp=(lidi-1-li-1di)/(di-1-di)；
3) as 0 < i < n, website SiFor intermediate stations；Upstream survey station Si-1, downstream survey station Si+1, water level is respectively li-1、
li+1；Website P is calculated separately to straight line SiSi-1, straight line SiSi+1Intersection point F1、F2；Including following content:
(1) work as F1、F2Fall in or do not fall in simultaneously simultaneously respectively line segment SiSi-1、SiSi+1When upper, l was takenp=li；
(2) work as F1Fall in line segment SiSi-1Upper and F2Line segment S is not fallen iniSi+1When upper, intersection point F1With website Si、Si-1Distance
Respectively d1i、di-1, meet (lp–li)/(li-1–lp)=d1i/di-1, then the linear interpolation water level l of boundary point Pp=(li-1d1i+
lidi-1)/(d1i+di-1)；
(3) work as F1Line segment S is not fallen iniSi-1Upper and F2In line segment SiSi+1When upper, intersection point F2With website Si、Si-1Distance point
It Wei not d2i、di+1, meet (li–lp)/(lp–li+1)=d2i/di+1, then the linear interpolation water level l of boundary point Pp=(li-1d2i+
lidi+1)/(d2i+di+1)；
S44: l is usedpValue update the elevation of boundary point P comprising the triangular apex elevation of boundary point P also it is corresponding more
Newly；
S45: step S41- step S44 is repeated, all Polygonal Boundary points are disposed after density.
In the above-mentioned technical solutions, in step 5, after the boundary point water level elevation interpolation and adjustment by step 4, originally
The coplanar triangulation network is changed into three-dimension curved surface；In each triangle of three-dimension curved surface, with including color translucent, texture
The conventional method of water or vertex water draws the water surface, can form the three-dimensional water surface of fitting natural river course water level.
The present invention has the advantage that
(1) present invention can guarantee that water surface effect of visualization is truer, flexible and practical；The prior art is not overcome generally not
The water surface gradient characteristic for taking natural river course into account can not embody lacking for the water level elevation progressive formation along river between discharge site
Point, while overcoming existing water surface method for visualizing and being difficult to be combined with the water level process data in true river, thus can not
The shortcomings that three-dimensional visualization applied to the river basin level river water surface；
(2) present invention is adjusted using polygon trigonometric ratio, water level interpolation and elevation, realizes three-dimension curved surface to natural river course water
The Continuous Approximation in face accounts for river upstream and downstream there are this objective fact of water surface gradient, and display effect is truer；
(3) present invention can encrypt the setting of spacing by Polygonal Boundary point, the number of Polygonal Boundary point after control encryption
Amount, and then the three-dimensional simulation water surface is adjusted flexibly to the approximation ratio of the true river water surface, facilitate user autonomous in efficiency and effect
Selection；
(4) the compatible traditional water surface method for visualizing of the present invention, pointedly solves the problems, such as river slope, to mention
The authenticity of the high river water surface, the present invention pass through the consecutive variations of interpolation calculation river water level, can be applied to flood analysis, dynamic
Emulation etc. has the occasion of required precision to water level expression, while being all suitable under small scene and river basin level large scene；
(5) strong applicability of the present invention can be applied to the industries such as Water Conservancy Information, water environment protection, geography information.
Detailed description of the invention
Fig. 1 is S of the present inventioniTo originate survey station schematic diagram.
Fig. 2 is S of the present inventioniTo terminate survey station schematic diagram.
Fig. 3 is that the situation that intersection point of the present invention is fallen in simultaneously respectively on the line segment of current survey station and upstream and downstream survey station line is illustrated
Figure.
Fig. 4 falls in current survey station for intersection point of the present invention respectively while not and shows with the situation on the line segment of upstream and downstream survey station line
It is intended to.
Fig. 5 is that intersection point of the present invention is only fallen in and the situation schematic diagram on the line segment of upstream survey station line.
Fig. 6 is that intersection point of the present invention is only fallen in and the situation schematic diagram on the line segment of downstream survey station line.
Fig. 7 is reservoir of embodiment of the present invention water surface boundary and survey station distribution map.
Fig. 8 is water level of embodiment of the present invention interpolation triangulation diagram adjusted.
Fig. 9 is the translucent water surface three-dimensional visualization effect picture of the embodiment of the present invention.
Figure 10 is process flow chart of the invention.
In Fig. 1, P is Polygonal Boundary point, lpFor the interpolation water level of Polygonal Boundary point, F0It is Polygonal Boundary point to directly
Line SiSi+1Intersection point, liFor the water level at the survey station a certain moment nearest with current border point, li+1For downstream survey station Si+1Water level,
SiFor the survey station nearest with current polygon boundary point (i=0, then SiTo originate survey station), Si+1For SiDownstream survey station, diIt is more
Side shape boundary point P to straight line SiSi+1Intersection point F0With website SiDistance, di+1For Polygonal Boundary point P to straight line SiSi+1Hang down
Sufficient F0With website Si+1Distance.
In Fig. 2, P is Polygonal Boundary point, lpFor the interpolation water level of Polygonal Boundary point, liIt is nearest with current border point
The survey station a certain moment water level, li-1For Si-1Water level, F0For Polygonal Boundary point to straight line SiSi-1Intersection point, SiFor with it is current
The nearest survey station of Polygonal Boundary point (i=n, then website SiTo terminate survey station), Si-1For SiUpstream survey station, diFor polygon edge
Boundary point P to straight line SiSi-1Intersection point F0With website SiDistance, di-1For Polygonal Boundary point P to straight line SiSi-1Intersection point F0With
Website Si-1Distance.
In Fig. 3, Fig. 4, Fig. 5, Fig. 6, P is Polygonal Boundary point, lpFor the interpolation water level of Polygonal Boundary point, liFor with work as
The water level at the nearest survey station of fore boundary point a certain moment, SiFor the survey station nearest with current polygon boundary point, (0 < i < n, then stand
Point SiFor intermediate stations), Si-1For SiUpstream survey station, Si+1For SiDownstream survey station, li-1For Si-1Water level, F1For polygon edge
Boundary point P to straight line SiSi-1Intersection point, F2For Polygonal Boundary point P to straight line SiSi+1Intersection point；
In Fig. 3 and Fig. 4, F1、F2Fall in or do not fall in simultaneously simultaneously respectively line segment SiSi-1、SiSi+1On, take lp=li；
In Fig. 5, F1Fall in line segment SiSi-1Upper and F2Line segment S is not fallen iniSi+1On, d1iFor intersection point F1With website SiDistance；
di-1For intersection point F1With website Si-1Distance；
In Fig. 6, F1Line segment S is not fallen iniSi-1Upper and F2In line segment SiSi+1On, d2iFor intersection point F2With website SiDistance, di+1
For intersection point F2With website Si-1Distance.
Fig. 7 is the distribution schematic diagram of reservoir of embodiment of the present invention water surface boundary and hydrographic survey website, and visual angle is positive
Top is overlooked.
Fig. 8 is the polygon triangle networking centre achievement that the embodiment of the present invention carries out water surface drafting using this method, visual angle
For above downstream towards upstream.
Fig. 9 is the effect picture that the embodiment of the present invention carries out the drafting of the color translucent water surface using this method, and visual angle is downstream
Top is towards upstream.
In Fig. 8 and Fig. 9, A is reservoir of embodiment of the present invention upstream survey station；B is mining under reservoir of the embodiment of the present invention
Certain survey station；In Fig. 8 and Fig. 9, survey station is set on the water surface and/or is set to water surface two sides.
Specific embodiment
The performance that the invention will now be described in detail with reference to the accompanying drawings, but they and do not constitute a limitation of the invention, only
It is for example.Keep advantages of the present invention more clear by explanation simultaneously and is readily appreciated that.
Refering to known to attached drawing: take the river water surface three-dimensional visualization method of water surface gradient into account, includes the following steps,
Step 1: data preparation；
Step 2: boundary point density is carried out to water surface boundary polygon；
Step 3: polygon trigonometric ratio；
Step 4: interpolation and adjustment are carried out based on elevation of the water level process data to boundary point；
Step 5: water surface drafting (as shown in Figure 10) is carried out based on the triangulation network adjusted.
In step 1, the data of preparation are river available data, and the river available data includes that river water surface boundary is more
Hydrographic station stage hydrograph data along side graphic data, river；Survey station stage hydrograph data include survey station along river
Number, x coordinate, y-coordinate, water level, time；Wherein, survey station number is successively increased since 0 according to the sequence for being from upstream to downstream
It is integer to n, n and is greater than or equal to 0.
In step 2, by the way that density point spacing is arranged, the boundary point quantity after density is controlled；Density back boundary point quantity
More, three-dimensional water surface flow display effect is better；Boundary point density method is the prior art；Spacing is that one of density method is defeated
Enter；Spacing setting is smaller, and the Polygonal Boundary point after density is more, and visualization expense is bigger, and water surface smooth effect is better；Density
Point spacing depends on the computer performance of user and the demand to effect of visualization；Computer performance is strong, effect of visualization requires
Lesser density point spacing is then arranged in height.
In step 3, trigonometric ratio is carried out to the boundary polygon after density, is converted the polygon after density to a series of
Triangle net collective；Polygon triangle turns to the prior art.
In step 4, the algorithm of boundary point water level elevation interpolation and adjustment includes the following steps:
S41: the Polygonal Boundary point P after traversal density；
S42: it calculates and current polygon boundary point P nearest survey station Si, number is i, if the water level at a certain moment is
li；
S43: the interpolation water level l of Polygonal Boundary point P is calculated by neighboring stations water levelp；Interpolation water level lpIn following
Hold:
1) as i=0, website SiTo originate survey station, no upstream survey station, downstream survey station Si+1, as shown in Figure 1；It calculates separately
Website P to straight line SiSi+1Intersection point F0With website Si、Si+1Distance di、di+1；If Si+1Water level is li+1, boundary point P's is linear
Interpolation water level lpMeet (lp–li)/(lp–li+1)=di/di+1, then lp=(lidi+1-li+1di)/(di+1-di)；
2) as i=n, then website SiTo terminate survey station, no downstream survey station, upstream survey station Si-1, as shown in Figure 2；It counts respectively
Calculate website P to straight line SiSi-1Intersection point F0With website Si、Si-1Distance di、di-1；If Si-1Water level is li-1, the line of boundary point P
Property interpolation water level lpMeet (lp–li)/(lp–li-1)=di/di-1, then lp=(lidi-1-li-1di)/(di-1-di)；
3) as 0 < i < n, website SiFor intermediate stations；Upstream survey station Si-1, downstream survey station Si+1, water level is respectively li-1、
li+1；Website P is calculated separately to straight line SiSi-1, straight line SiSi+1Intersection point F1、F2；Including following content:
(1) work as F1、F2Fall in or do not fall in simultaneously simultaneously respectively line segment SiSi-1、SiSi+1When upper, as shown in Figure 3, Figure 4, took
lp=li；
(2) work as F1Fall in line segment SiSi-1Upper and F2Line segment S is not fallen iniSi+1When upper, as shown in figure 5, intersection point F1With website Si、
Si-1Distance be respectively d1i、di-1, meet (lp–li)/(li-1–lp)=d1i/di-1, then the linear interpolation water level l of boundary point Pp=
(li-1d1i+lidi-1)/(d1i+di-1)；
(3) work as F1Line segment S is not fallen iniSi-1Upper and F2In line segment SiSi+1When upper, as shown in fig. 6, intersection point F2With website Si、
Si-1Distance be respectively d2i、di+1, meet (li–lp)/(lp–li+1)=d2i/di+1, then the linear interpolation water level l of boundary point Pp=
(li-1d2i+lidi+1)/(d2i+di+1)；
S44: l is usedpValue update the elevation of boundary point P comprising the triangular apex elevation of boundary point P also it is corresponding more
Newly；
S45: step S41- step S44 is repeated, all Polygonal Boundary points are disposed after density.
Above each step relies on computer and software to complete.
The present invention takes water surface gradient into account and is mainly reflected in: water level elevation after the water surface boundary point interpolation of step 4, is adjacent
The linear interpolation of upper and lower survey station water level elevation；It is higher using interpolation method precision of the present invention；It is inserted used by step 4 of the present invention
Value-based algorithm is not limited only to linear interpolation for the purpose of the smooth transition of guaranteed water level.
In step 5, after the boundary point water level elevation interpolation and adjustment by step 4, originally coplanar triangulation network transformation
For three-dimension curved surface；In each triangle of three-dimension curved surface, with conventional methods such as color translucent, texture water or vertex water
The water surface is drawn, the visualization of 3 d water surface of fitting natural river course water level can be formed.
Embodiment
Now apply the river water surface three-dimensional visualization method for taking water surface gradient into account of Mr. Yu's reservoir as embodiment using the present invention
It is described in detail, the river water surface three-dimensional visualization method for taking water surface gradient into account for being applied to other waters to the present invention is same
With directive function.
The present embodiment is using certain reservoir water surface boundary line and hydrographic station measured water level data as available data.Choose certain water
The water surface boundary line of two station of library the east of Sichuan to damsite Sandouping part, and 12 survey stations when having chosen the 8:00 on the 4th of August in 2012
Water level as input, swum over under upstream be respectively damsite Sandouping 158.68m, peaceful small stream 158.69m, osiery small stream 158.69m,
Nine curved small stream 158.71m, new beach 158.72m, fragrant small stream 158.73m, Zigui station 158.74m, Sha Zhenxi 158.75m, the town Guo Tan
158.75m, Xie Jiahe 158.77m, camel river 158.79m, two station 158.81m of the east of Sichuan, survey station distribution are as shown in Figure 7.
Step 1: data preparation,
Prepare the river water surface boundary polygon data measured, hydrographic station stage hydrograph data along river；
Survey station stage hydrograph data include along river: survey station number, xy coordinate, water level, time；Wherein, survey station number according to from
On swim over to the sequence in downstream and successively increased from 0 and be integer to n, n and be greater than or equal to 0；
Step 2: water surface boundary polygon boundary point density,
By the way that density point spacing is arranged, the boundary point quantity after density is controlled；Density back boundary point quantity is more, three-dimensional
Water surface flow display effect is better；Density point spacing depends on the computer performance of user and the demand to effect of visualization；Meter
Calculation machine performance is strong, effect of visualization requires height, then lesser density point spacing is arranged；
Step 3: polygon trigonometric ratio；A series of triangle net collectives are converted by the polygon after density；
Step 4: boundary point water level elevation interpolation and adjustment based on stage hydrograph,
Algorithm flow is as follows:
S41: the Polygonal Boundary point P after traversal density；
S42: it calculates and current border point P nearest survey station Si, number is i, the water level l at a certain momenti；
S43: the interpolation water level l of boundary point P is calculated by neighboring stations water levelp；Including following three categories situation:
1) as i=0, then website SiTo originate survey station, no upstream survey station, downstream survey station Si+1, as shown in Figure 1；It counts respectively
Calculate website P to straight line SiSi+1Intersection point F0With website Si、Si+1Distance di、di+1.If Si+1Water level is li+1, the line of boundary point P
Property interpolation water level lpMeet (lp–li)/(lp–li+1)=di/di+1, then lp=(lidi+1-li+1di)/(di+1-di)；
2) as i=n, then website SiTo terminate survey station, no downstream survey station, upstream survey station Si-1, as shown in Figure 2；It counts respectively
Calculate website P to straight line SiSi-1Intersection point F0With website Si、Si-1Distance di、di-1.If Si-1Water level is li-1, the line of boundary point P
Property interpolation water level lpMeet (lp–li)/(lp–li-1)=di/di-1, then lp=(lidi-1-li-1di)/(di-1-di)；
3) as 0 < i < n, then website SiFor intermediate stations；Upstream survey station Si-1, downstream survey station Si+1, water level is respectively li-1、
li+1.Website P is calculated separately to straight line SiSi-1, straight line SiSi+1Intersection point F1、F2, include the following three types situation:
If F1、F2Fall in or do not fall in simultaneously simultaneously respectively line segment SiSi-1、SiSi+1On, as shown in Figure 3, Figure 4, take lp=
li；
If F1Fall in line segment SiSi-1Upper and F2Line segment S is not fallen iniSi+1On, as shown in figure 5, intersection point F1With website Si、Si-1's
Distance d1i、di-1, meet (lp–li)/(li-1–lp)=d1i/di-1, then the linear interpolation water level l of boundary point Pp=(li-1d1i+
lidi-1)/(d1i+di-1)；
If F1Line segment S is not fallen iniSi-1Upper and F2In line segment SiSi+1On, as shown in fig. 6, intersection point F2With website Si、Si-1Away from
From d2i、di+1, meet (li–lp)/(lp–li+1)=d2i/di+1, then the linear interpolation water level l of boundary point Pp=(li-1d2i+
lidi+1)/(d2i+di+1)；
S44: l is usedpValue updates the elevation of boundary point P, and the triangular apex elevation comprising boundary point P also accordingly updates.
S45: step S41-S44 is repeated, all Polygonal Boundary points are disposed after density；
Step 5: the triangulation network water surface is drawn, the three-dimensional visualization surface chart of the continuous water surface of certain reservoir is obtained.
After boundary point water level elevation interpolation and adjustment, the originally coplanar triangulation network is changed into three-dimension curved surface, in three-dimensional
In each triangle of curved surface, the water surface is drawn with traditional method such as the methods of color translucent, texture water or vertex water,
The three-dimensional visualization surface chart of fitting natural river course water level can be formed.
Conclusion: encryption, trigonometric ratio and the height of water level interpolation on water surface boundary of the present embodiment Jing Guo the method for the present invention adjust
Later, original plane polygon becomes taking into account the three-dimension curved surface of river slope；The water level interpolation water surface adjusted of the present invention is shown
(the water level difference of each website has been done to be exaggerated processing (exaggerating processing is well-known technique, is by the water level of each hydrographic station effect
Data are greater than one amplification coefficient multiplied by one, so as to which the water level difference between adjacent sites is highlighted), to guarantee three
Dimension is shown intuitive), as shown in Figure 8, Figure 9, it is seen that the method for the present invention accounts for upstream to downstream, and there are the natural of water surface gradient
Characteristic can be realized the three-dimensional visualization of the continuous water surface.
Other unaccounted parts belong to the prior art.
Claims (6)
1. taking the river water surface three-dimensional visualization method of water surface gradient into account, it is characterised in that: include the following steps,
Step 1: data preparation；
Step 2: boundary point density is carried out to water surface boundary polygon；
Step 3: polygon trigonometric ratio；
Step 4: interpolation and adjustment are carried out based on elevation of the water level process data to boundary point；
Step 5: water surface drafting is carried out based on the triangulation network adjusted.
2. the river water surface three-dimensional visualization method according to claim 1 for taking water surface gradient into account, it is characterised in that: step
In one, the data of preparation are river available data, and the river available data includes river water surface boundary polygon data, river
Hydrographic station stage hydrograph data along the line；Hydrographic station stage hydrograph data include that hydrographic station is compiled along river
Number, x coordinate, y-coordinate, water level, the time；Wherein, hydrographic station number is according to being from upstream to the sequence in downstream since 0 successively
N is increased to, n is integer and is greater than or equal to 0.
3. the river water surface three-dimensional visualization method according to claim 1 or 2 for taking water surface gradient into account, it is characterised in that:
In step 2, by the way that density point spacing is arranged, the boundary point quantity after density is controlled.
4. the river water surface three-dimensional visualization method according to claim 3 for taking water surface gradient into account, it is characterised in that: step
In three, trigonometric ratio is carried out to the boundary polygon after density, converts a series of triangle net collectives for the polygon after density.
5. the river water surface three-dimensional visualization method according to claim 4 for taking water surface gradient into account, it is characterised in that: step
In four, the algorithm of boundary point water level elevation interpolation and adjustment includes the following steps:
S41: the Polygonal Boundary point P after traversal density；
S42: it calculates and current polygon boundary point P nearest hydrographic station Si, number is i, if the water level at a certain moment is
li；
S43: the interpolation water level l of Polygonal Boundary point P is calculated by adjacent hydrographic station water levelp；Interpolation water level lpIncluding as follows
Content:
1) as i=0, website SiTo originate hydrographic station, no upper water text measuring station, downstream hydrographic station Si+1；Respectively
Calculate website P to straight line SiSi+1Intersection point F0With website Si、Si+1Distance di、di+1；If Si+1Water level is li+1, boundary point P's
Linear interpolation water level lpMeet (lp–li)/(lp–li+1)=di/di+1, then lp=(lidi+1-li+1di)/(di+1-di)；
2) as i=n, then website SiTo terminate hydrographic station, no downstream hydrographic station, upstream hydrographic station Si-1；Point
It Ji Suan not website P to straight line SiSi-1Intersection point F0With website Si、Si-1Distance di、di-1；If Si-1Water level is li-1, boundary point P
Linear interpolation water level lpMeet (lp–li)/(lp–li-1)=di/di-1, then lp=(lidi-1-li-1di)/(di-1-di)；
3) as 0 < i < n, website SiFor intermediate stations；Upstream hydrographic station Si-1, downstream hydrographic station Si+1, water level difference
For li-1、li+1；Website P is calculated separately to straight line SiSi-1, straight line SiSi+1Intersection point F1、F2；Including following content:
(1) work as F1、F2Fall in or do not fall in simultaneously simultaneously respectively line segment SiSi-1、SiSi+1When upper, l was takenp=li；
(2) work as F1Fall in line segment SiSi-1Upper and F2Line segment S is not fallen iniSi+1When upper, intersection point F1With website Si、Si-1Distance difference
For d1i、di-1, meet (lp–li)/(li-1–lp)=d1i/di-1, then the linear interpolation water level l of boundary point Pp=(li-1d1i+
lidi-1)/(d1i+di-1)；
(3) work as F1Line segment S is not fallen iniSi-1Upper and F2In line segment SiSi+1When upper, intersection point F2With website Si、Si-1Distance be respectively
d2i、di+1, meet (li–lp)/(lp–li+1)=d2i/di+1, then the linear interpolation water level l of boundary point Pp=(li-1d2i+lidi+1)/
(d2i+di+1)；
S44: l is usedpValue updates the elevation of boundary point P comprising the triangular apex elevation of boundary point P also accordingly updates；
S45: step S41- step S44 is repeated, all Polygonal Boundary points are disposed after density.
6. the river water surface three-dimensional visualization method according to claim 5 for taking water surface gradient into account, it is characterised in that: step
In five, after the boundary point water level elevation interpolation and adjustment by step 4, the originally coplanar triangulation network is changed into three-dimension curved surface；?
In each triangle of three-dimension curved surface, the water surface is drawn with the conventional method for including color translucent, texture water or vertex water,
The three-dimensional water surface of fitting natural river course water level can be formed.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910278875.0A CN110070619B (en) | 2019-04-09 | 2019-04-09 | Three-dimensional visualization method for river water surface considering water surface gradient |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910278875.0A CN110070619B (en) | 2019-04-09 | 2019-04-09 | Three-dimensional visualization method for river water surface considering water surface gradient |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110070619A true CN110070619A (en) | 2019-07-30 |
| CN110070619B CN110070619B (en) | 2021-07-06 |
Family
ID=67367127
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910278875.0A Active CN110070619B (en) | 2019-04-09 | 2019-04-09 | Three-dimensional visualization method for river water surface considering water surface gradient |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110070619B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101359052B (en) * | 2008-09-02 | 2011-02-16 | 武汉大学 | Storage monitoring method |
| CN102982494A (en) * | 2012-10-31 | 2013-03-20 | 长江航道局 | Channel depth information generating method based on electronic navigation chart and water level data |
| CN103143465A (en) * | 2013-02-25 | 2013-06-12 | 中国水利水电科学研究院 | System and method for simulating regional rainfall process |
| CN103268633A (en) * | 2013-04-24 | 2013-08-28 | 中国科学院遥感与数字地球研究所 | Contour surface construction method for raster data |
| CN105741045A (en) * | 2016-02-05 | 2016-07-06 | 中国水利水电科学研究院 | Flood risk dynamic analysis and display system and method |
-
2019
- 2019-04-09 CN CN201910278875.0A patent/CN110070619B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101359052B (en) * | 2008-09-02 | 2011-02-16 | 武汉大学 | Storage monitoring method |
| CN102982494A (en) * | 2012-10-31 | 2013-03-20 | 长江航道局 | Channel depth information generating method based on electronic navigation chart and water level data |
| CN103143465A (en) * | 2013-02-25 | 2013-06-12 | 中国水利水电科学研究院 | System and method for simulating regional rainfall process |
| CN103268633A (en) * | 2013-04-24 | 2013-08-28 | 中国科学院遥感与数字地球研究所 | Contour surface construction method for raster data |
| CN105741045A (en) * | 2016-02-05 | 2016-07-06 | 中国水利水电科学研究院 | Flood risk dynamic analysis and display system and method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110070619B (en) | 2021-07-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102663827B (en) | Three-dimensional dynamic whole-process simulation method for storm surge and flood routing in complex flooding areas | |
| Hnaidi et al. | Feature based terrain generation using diffusion equation | |
| US7412362B2 (en) | River modeling | |
| CN106570929B (en) | Construction and drawing method of dynamic volume cloud | |
| Mei et al. | Fast hydraulic erosion simulation and visualization on GPU | |
| Peytavie et al. | Arches: a framework for modeling complex terrains | |
| Roudier et al. | Landscapes synthesis achieved through erosion and deposition process simulation | |
| CN107403459A (en) | Real terrain fast modeling method and landslide visualization technique | |
| Emilien et al. | Interactive procedural modelling of coherent waterfall scenes | |
| Anh et al. | Physically based hydraulic erosion simulation on graphics processing unit | |
| CN102663802A (en) | Game terrain road generating method and device | |
| CN104463937B (en) | Animation producing method and device based on irregularly shaped object | |
| CN102117501B (en) | Method and device for rendering three-dimensional lines | |
| CN104574489B (en) | Landform and motion vector integrated approach based on lamination quaternary tree atlas | |
| Semmo et al. | Real-time rendering of water surfaces with cartography-oriented design | |
| Rusnell et al. | Feature-rich distance-based terrain synthesis | |
| Trapp et al. | Interactive Rendering and Stylization of Transportation Networks using Distance Fields. | |
| CN110070619A (en) | Take the river water surface three-dimensional visualization method of water surface gradient into account | |
| Golubev et al. | Dijkstra-based terrain generation using advanced weight functions | |
| Wang et al. | Dynamic modeling and rendering of grass wagging in wind | |
| Thöny et al. | Large‐scale pixel‐precise deferred vector maps | |
| Teoh | Riverland: An efficient procedural modeling system for creating realistic-looking terrains | |
| Chang et al. | Hierarchical simplification of city models to maintain urban legibility. | |
| Toledo et al. | Visualizing 3D flow of black-oil reservoir models on arbitrary surfaces using projected 2D line integral convolution | |
| Grottel et al. | Real-Time Visualization of Urban Flood Simulation Data for Non-Professionals. |
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 |